带电粒子在磁场中的运动matlab代码

当一个带电粒子在磁场中运动时，它会受到洛伦兹力的作用。洛伦兹力的方向垂直于粒子的速度和磁场方向，大小与速度、电荷量以及磁场强度有关。下面是一个用MATLAB实现的简单示例代码，演示了带电粒子在磁场中的运动：

% 输入参数
q = 1.6e-19; % 粒子电荷量（库仑）
m = 1e-27; % 粒子质量（千克）
v0 = [0, 1e4, 0]; % 初始速度（米/秒）
B = [0, 0, 1]; % 磁场强度（特斯拉）

% 设置模拟参数
dt = 1e-8; % 时间步长（秒）
num_steps = 10000; % 模拟步数

% 初始化位置和速度
r = [0, 0, 0]; % 初始位置（米）
v = v0; % 初始速度（米/秒）

% 模拟带电粒子的运动
for i = 1:num_steps
    % 计算洛伦兹力
    F = q * cross(v, B);
    
    % 计算加速度
    a = F / m;
    
    % 更新速度和位置
    v = v + a * dt;
    r = r + v * dt;
    
    % 绘制粒子轨迹
    plot3(r(1), r(2), r(3), 'bo'); % 绘制蓝色点
    
    % 设置坐标轴范围
    xlim([-1e-5, 1e-5]);
    ylim([-1e-5, 1e-5]);
    zlim([-1e-5, 1e-5]);
    
    % 添加标签和标题
    xlabel('x');
    ylabel('y');
    zlabel('z');
    title('带电粒子在磁场中的运动');    
    % 更新绘图
    drawnow;
end

你可以根据需要修改输入参数、模拟参数以及绘图设置。这个代码将产生一个动画，显示带电粒子在三维空间中的运动轨迹。