Matlab 磁镜现象数值模拟：带电粒子束运动轨迹

以下是一个简单的 Matlab 程序，模拟磁镜现象：

% 磁镜模拟程序

% 定义常数
q = 1.6e-19;    % 电荷量
m = 9.11e-31;   % 粒子质量
B0 = 1;         % 磁场中心强度
L = 1e-6;       % 磁镜长度
v0 = 1e7;       % 初始速度

% 定义时间步长和模拟时间
dt = 1e-11;
tmax = 1e-8;
t = 0:dt:tmax;

% 定义初始位置和速度
x0 = -L/2;
y0 = 0;
vx0 = v0;
vy0 = 0;

% 定义粒子的运动方程
fx = @(t,x,y,vx,vy) q*B0*vy;
fy = @(t,x,y,vx,vy) -q*B0*vx;
f = @(t,y) [y(4); y(5); fx(t,y(1),y(2),y(4),y(5))/m; fy(t,y(1),y(2),y(4),y(5))/m; 0];

% 解微分方程
y0 = [x0; y0; vx0; vy0; 0];
options = odeset('RelTol',1e-6,'AbsTol',1e-9);
[t,y] = ode45(f,t,y0,options);

% 绘制粒子轨迹图
figure;
plot(y(:,1),y(:,2),'LineWidth',2);
xlabel('x (m)');
ylabel('y (m)');
title('磁镜模拟');

上述程序中，我们使用了 ODE45 函数来求解微分方程，其中 fx 和 fy 分别表示粒子在 x 和 y 方向上的受力，f 则是整个粒子运动方程。程序中的其它参数可以根据需要自行修改。运行程序后，可以得到类似下图所示的粒子轨迹图：

从轨迹图中可以看出，粒子在磁场中心周围呈圆形轨迹，而在两端则被强磁场弯曲，形成了‘磁镜’效应。