用matlab的语言编程，要求表示出带电粒子在匀强磁场中的运动规律，分别对带电粒子垂直入射磁场方向，平行磁场方向，任意角度进入磁场进行讨论以及作图（不要动画）

在MATLAB中，你可以使用数值积分和向量数学来模拟带电粒子在均匀磁场中的运动。首先，我们假设粒子的质量、电荷和速度都是已知的，磁场B也是常数。以下是针对三种情况的基本步骤：

1. **垂直入射**：
当粒子速度v垂直于磁场B时，粒子将做圆周运动（洛伦兹力提供向心力）。可以建立二维坐标系，x轴代表粒子轨迹，y轴代表垂直速度分量。通过牛顿第二定律和洛伦兹力公式，你可以求解角速度，并计算每个时间步的位置。绘制出粒子在y轴上的轨迹。

% 假设初速度 vy = v0 * sin(theta) (theta为入射角)
vy = v0 * sind(angle);
dx = zeros(t_max, 1); % x轴位置
dy = zeros(t_max, 1); % y轴位置

for t = 0:dt:t_max
    r = vy * dt; % 半径
    dx(end+1) = dx(end) + r * cosd(theta); % 更新x轴位置
    dy(end+1) = dy(end) + r; % 更新y轴位置
end
plot(dx, dy);
xlabel('x');
ylabel('y');

2. **平行入射**：
粒子将不受磁场影响，沿直线匀速移动。只需记录初始速度并随时间增加即可。

vx = v0 * cosd(angle); % 水平速度
dx = vx * t;
plot(dx, ones(size(dx)) * initial_y_position, 'r'); % 直线轨迹

3. **任意角度入射**：
这需要分解粒子速度到垂直和水平分量，然后分别应用上述两个情况的算法。可以使用三角函数来计算分量。

v_parallel = v0 * cosd(angle);
v_perpendicular = v0 * sind(angle);

dx = v_parallel * t;
dy = zeros(1, t_max); % 初始化y轴
for t = 1:length(t)
    dy(t) = dy(t-1) + v_perpendicular * t; % 圆周运动累积
end
plot(dx, dy);
hold on;