-94- /2012.11/ 科研发展 》》 MATLAB仿真带电粒子在磁场中磁镜现象 上海海事大学文理学院　方瑞银国家知识产权局专利局专利审查协作江苏中心　戚海洋 【摘要】借助MATLAB数学工具软件，使用数值计算的方法计算带电粒子磁场中的运动方程，对非均匀磁场的磁镜现象进行仿真，使人们对带电粒子在磁场磁镜现象有个更直观的了解。 【关键词】MATLAB；带电粒子；磁场；磁镜现象 1.引言 对于带电粒子在磁场中的磁镜现象，人们都有了一定的认识，但是却不能对其 有比较直观形象的想象。为了比较直观形象的认识带电粒子在磁场中的磁镜现象，本文利用MATLAB对其常见运动进行数值计算，并进行仿真，方便人们得到直观印象，为进一步工作学习提供方便。 2.带电粒子在均匀磁场中的运动方程 在均匀磁场中，忽略带电粒子所受重力的影响，带电粒子在磁场中只受到洛伦兹力的作用而使得运动轨迹发生相应的变化。根据牛顿第二定律，质量为m，电荷量为q，初速度为v 的带电粒子在磁感应强 度为B 的磁场中满足的运动方程为： 将 kzjyixr   和 kBjBiBB zyx   代入 上式，可将上式分解在直角坐标系展成标量式： (1) 为便于求解，可以将其化为一阶一阶方程组。令得到一组一阶微分方程组： (2) 3.磁场中的磁镜现象 在非均匀磁场中，粒子发生旋转的半径与磁感应强度成反比，图1示带电粒子在递增磁场的运动轨迹[3]。若带电粒子以一定角度入射到磁感应强度递增地磁场区，带电粒子受到的洛伦兹力可以分解为 1F 、 2 F ，如图2。当带电粒子从弱磁场区向强磁场区运动时，如图2恒有分力 1 F 阻止带电粒子向强磁场区的运动，于是当带电粒子运动一段时间后，粒子速度减小至零，粒子像被反弹一样开始反向运动，这种现象被称为磁镜现象，如图3。 4.MATLAB数值求解及仿真演示 设定磁场为沿z轴的轴对称场，如图 4， 0,0,0  zr BBB  。由 0 B  ，可得： (3) 在一个拉莫尔半径内，近似均匀磁场。由 |||| zr BB ，(3)可得： z Br B z r    2 - [5]。利用均匀磁场中直角坐 标系的方程形式(1)、(2)，假定带电粒子的电量q=1.6e-19，质量m=1.6725e-27，初速度V0=1e6，B0=1T， T/m；方程组 (2)的初始条件为： 0)0(,1)0(,0)0(  zyx ， 其中 4    。 MATLAB仿真源程序：clear; q=1.6e-19; m=1.6725e-27; V0=1e6;sita=pi/4; Vy=0,Vz=V0*cos(sita),Vx=V0*sin(s ita); c=[0,1.1,0,Vx,Vy,Vz];%给定初始量[t,y]=ode45(‘cjxxfun’,[0:1e- 7:1.5e-6],c,[],q,m);%解微分方程组plot3(y(:,1),y(:,2),y(:,3),’lin ewidth’,2);%绘制三维粒子运动轨迹 xlabel('z');ylabel('y');zlabel('x'); grid on;title(‘磁镜现象’) figure subplot(2,1,1); plot(y(:,1),y(:,2),’linewid th’,2);%从yoz面看到的轨迹 xlabel('z');ylabel('y'); grid on;title(‘yoz平面’) subplot(2,1,2); plot(y(:,1),y(:,