matlab带电粒子在正交电磁场的螺旋运动，入射角为锐角，带注释

下面是一个简单的 Matlab 代码，模拟了一个带电粒子在正交电磁场中的螺旋运动。注释已经添加在代码中。

% 定义常数
q = 1.6e-19;    % 电荷量
m = 9.11e-31;   % 粒子质量
B = 1;          % 磁场强度
E = 1;          % 电场强度

% 定义初值
x0 = 0;
y0 = 0;
z0 = 0;

vx0 = 1;
vy0 = 0;
vz0 = 1;

% 定义时间步长和总时间
dt = 0.01;      % 时间步长
t = 0:dt:10;    % 总时间

% 定义数组
x = zeros(size(t));
y = zeros(size(t));
z = zeros(size(t));

vx = zeros(size(t));
vy = zeros(size(t));
vz = zeros(size(t));

% 初值赋值
x(1) = x0;
y(1) = y0;
z(1) = z0;

vx(1) = vx0;
vy(1) = vy0;
vz(1) = vz0;

% 计算粒子位置和速度
for i = 2:length(t)
    % 计算电磁力和加速度
    F = q * (E*[1 0 0]' + cross([vx(i-1) vy(i-1) vz(i-1)]', [0 0 B]'));
    ax = F(1)/m;
    ay = F(2)/m;
    az = F(3)/m;
    
    % 计算速度和位置
    vx(i) = vx(i-1) + ax * dt;
    vy(i) = vy(i-1) + ay * dt;
    vz(i) = vz(i-1) + az * dt;
    
    x(i) = x(i-1) + vx(i) * dt;
    y(i) = y(i-1) + vy(i) * dt;
    z(i) = z(i-1) + vz(i) * dt;
end

% 绘图
plot3(x, y, z);
xlabel('x');
ylabel('y');
zlabel('z');
title('带电粒子在正交电磁场中的螺旋运动');

这个程序中，我们首先定义了粒子的常数、初值、时间步长和总时间，并且定义了用来存储位置和速度的数组。然后，在每个时间步长内，我们计算了电磁力和加速度，并且使用欧拉法计算了速度和位置。最后，我们绘制了三维图形来展示粒子的螺旋运动。