matlab磁镜实验代码

以下是一个简单的 Matlab 磁镜实验代码：

% 定义常数
m = 9.10938356e-31; % 电子质量
q = -1.60217662e-19; % 电子电荷
B = 0.5; % 磁场强度

% 定义初始条件
x0 = 0;
y0 = 0;
z0 = 0;
vx0 = 1e5;
vy0 = 1e5;
vz0 = 0;

% 定义时间步长和总时间
dt = 1e-10;
t = 1e-5;

% 初始化位置和速度数组
x = zeros(1, t/dt);
y = zeros(1, t/dt);
z = zeros(1, t/dt);
vx = zeros(1, t/dt);
vy = zeros(1, t/dt);
vz = zeros(1, t/dt);

% 初始条件
x(1) = x0;
y(1) = y0;
z(1) = z0;
vx(1) = vx0;
vy(1) = vy0;
vz(1) = vz0;

% 计算电子运动轨迹
for i = 1:t/dt-1
    ax = q*B/m*vy(i);
    ay = -q*B/m*vx(i);
    az = 0;
    vx(i+1) = vx(i) + ax*dt;
    vy(i+1) = vy(i) + ay*dt;
    vz(i+1) = vz(i) + az*dt;
    x(i+1) = x(i) + vx(i+1)*dt;
    y(i+1) = y(i) + vy(i+1)*dt;
    z(i+1) = z(i) + vz(i+1)*dt;
end

% 绘制电子运动轨迹
figure;
plot(x, y);
axis equal;
title('Electron Trajectory in Magnetic Mirror');
xlabel('x');
ylabel('y');

这个代码使用了欧拉法（Euler Method）来计算电子在磁场中的运动轨迹。在这个简单的磁镜实验中，我们假设电子的初始速度沿着 x 轴方向，且电子在 y-z 平面内运动。运行代码后，程序将绘制出电子在磁镜中的运动轨迹。