利用Matlab模拟带电粒子在磁场中的运动

1、利用Matlab模拟带电粒子在磁场中的运动一、实验目的1、理解数值模拟研究物理问题的思路，能独立地运用此方法研究物理问题,掌握数值模拟的编程。2、运用Matlab数值模拟的方法研究三维空间中带电粒子在复杂磁场环境下的运动行为。二、实验原理带电粒子在磁场中运动时会受到洛伦兹力的作用，且随着初始运动方向和磁场分布的不同，其运动轨迹会发生不同的变化。由洛伦兹力的推导公式可知，它垂直于粒子的运动速度，不对运动粒子作功,只改变其运动方向，其大小为：F=qvBsin&；因此，综合牛顿运动定律就可以精确确定带电粒子在磁场中的运动轨迹。三、实验内容用Matlab数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动

2、，即带电粒子进入磁场的方向与磁场方向的角度&(Ov&v9O。)。带电粒子质量为m,电量为q(q0)的带电粒子在均匀稳定的电磁场中的运动微分方程为：md2r=qE+qvxB(1)以电磁场中某点为原点，以E为Oy方向,B为Oz方向建立坐标系0-粧。由于a)=qB/mf则(1)式的投影方程为:d2z=0drindt将其转换为一阶微分线性方程组，以便用差分替代微分作数值计算，令w=x,w2=dy/dx,w3=y,w4=dy/dt,巴=Z，%=ck/dt,则(2)式化简为:dwdtdwy苻=7dw.二=w.dt4dw.qE=-cow.dtmd叫页=%叫=0.dt下面使用MATLABR2009b编写程序，

3、分三种情况考虑：电场强度和磁场强度都不为零；电场强度为零，磁场强度不为零；电场强度不为零，磁场强度为零。源程序如下：治分3种情况模拟带电粒子在恒定磁场中的螺旋运动globalqmbeq=l6e-2;m=002;b=2;l;0;e=l；0;l;figurestrd1=1eneq0,bneq01;strd2=1e=0,bneq0;strd3=,eneq0,b=01;fori=l:3t/w=ode23(ddlzfun/0:0.01:20,0,0.01,0,6,0,0.01,qzmzb(i),e(iaxes(1unit,1normalized1,1position,0.045+(i-1)*0350.0

4、62027860.6583);plot3(w(:,1)zw(:z3)zw(:,5)z1linewidth1,2);gridontitle(strdiffontsize,12,1fontweight1fdemi);xlabel(1x1);ylabel(1y1);zlabel(1z);view(-51,18);endfunctionwdot=ddlzfun%该函数实现(3)式的差分迭代symsqmbeq=l6e-2;m=002;b=2;l;0J;e=l；0;l;wdot=w(2);q.*b.*w(4)./m;w(4);q.*e./m-q.*b.*w(2)./m;w(6);0运行以上程序得到该带电粒

5、子在均匀稳定的电磁场中的运动轨迹如图1所示。图1带电粒子在均匀稳定的电磁场中的运动轨迹（2）当v与B有一定夹角&时，可将V分解为：v/z=vcos,v=vsin（即匕工0,比HO,冬h0）。若只有分量u，带电粒子将在垂直的平面内作匀速圆周运动，若只有分量，粒子将沿磁场B方向作匀速直线运动，当两分量同时存在时，带电粒子的运动轨迹将是一条螺旋线。使用MATLABR2009b数值模拟的方法编写程序，源代码如下：治用Matlab数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动functionlxglobalqmBt,y=ode45(ddlzfunz0:0.01:20z0,0.010,6,0,0.01z

6、zqzmzBz0);治用。de23解微分方程组，时间设为20s治指定初始条件，传递相关参数comet3(y(:zl)zy(:,3),y(:,5);plot3(y(:,1)zy(:z3)zy(:,5);gridon殆开启坐标网格线xlabel(1x1);ylabel(y);zlabel(1z);titled模拟带电粒子在恒定磁场中的螺旋运动，)；functionydot=ddlzfun(t,yrqzm,BrE)E为参量symsqmBq=l6e-2;B=2;m=002;ydot=y(2);q*B*y(4)/m;y(4);-q*B*y(2)/m;y(6);0;运行以上程序，得到带电粒子的运动轨迹(螺

7、旋线)如图2所示。模拟带电粒子在恒定磁场中的蝮旋运动图2带电粒子的运动轨迹(螺旋线)用Matlab数值模拟的方法模拟磁聚焦现象，即在均匀磁场中某点引入一发散角&不大的带电粒子束，并使束中粒子的速度v大致相同。对图2所示的螺旋线进行简单的分析：其螺距为:h=qT=汕,带电qB粒子运动一周所前进的距离与v无关，所以若从磁场中某点A发射出一束很窄的电子流，使他们的速度很接近，并且与的夹角都很小，则V/y=VCOS-Vo它们具有近似相同的螺距力，尽管它们的v=vsin-v不同，各粒子会沿着不同的半径作螺旋线运动，但各粒子在经过距离后乂会重新聚在一起，这就是磁聚焦现象。用Matlab数值模拟的方法编写如

8、下的源程序模拟磁聚焦现象：治用Matlab数值模拟的方法模拟磁聚焦现象t=0:0.01:2*pi;al=O.5.*(t-pi);form=-16:2:10*pi/180;axis(07-11-0.40.4);gridon;view(12,18);holdon;comet3(cos(m)*tt2*sin(m)*cos(al)A292*sin(m)*cos(al)*sin(al);plot3(cos(m)*tz2*sin(m)*cos(al)A2/2*sin(m)*cos(al)*sin(al);endxlabel(1xf);ylabel(1y1);zlabel(fz1);title(模拟磁聚焦现

9、象)以上程序中默认粒子入射速度相同，运行结果是一束带电粒子做螺旋运动的三维动画，我们可以从不同的视角进行观察，并可发现当发散角不大时粒子确实会聚到同一点，如图3；在平行磁场方向可以看到粒子做圆周运动，如图4所示。!1!图4磁聚焦现象（Y-Z方向）由于电磁场看不见，摸不着，它不像普通的“三态”物质那样由原子、分子构成，也没有可见的形态，但其具有可以被检测的运动速度、能量和动量，占有空间，是一种真实的客观存在。实验中通过MATLAB数值模拟模拟带电粒子在恒定磁场中的螺旋运动和磁聚焦现象，让我们对带电粒子在磁场中的运动有了更深的感性认识，用MATLAB画出的立体图也更有利于对此的理解，对于对应知识的理解和吸收有很大的帮助。运用数值模拟方法作出的螺旋线运动和磁聚焦现象的轨迹不仅生动形象,而且结果是建立在定量计算的基础上的，不是简单的定性的模拟，因此更具有科学性与说服力。五、参考文献宿刚，