带电粒子在互相垂直的恒定电磁场下的漂移.docxVIP

带电粒子在互相垂直的恒定电磁场下的漂移 1,实验目的： 通过粒子运动轨迹，分析带点粒子运动规律。 2,实验内容： 利用数值模拟计算粒子运动轨迹，并实现可视化，并分析带电粒子运动规律。 3,实验过程： （1），设定参数：荷质比Q/m=1 E(V/m) B(T) V(m/s) x 0 0 0 y 0.1 0 0 z 0 1 0 (2),运算公式： 引入磁标势A以及电势φ： B=?× -?φ= 设广义势能为： U=q 则拉格朗日函数为： L= 使电场E=Ej,磁场 A 得拉格朗日函数为： L= 代入拉格朗日方程ddt x 解此微分方程组即得结果。 四，结果： （1），轨迹： 带电粒子运动轨迹如图： 可见，带电粒子大致朝着x轴方向螺旋前进，可认为粒子轨迹与y(x)为周期函数。 如右图所示，为t-x\y图： 可见粒子沿着y轴方向做简谐振动，并且以恒定速率沿着x轴方向运动时做简谐振动；对x和y进行傅里叶分析，如图所示： 可见，x和y有相同的固定单一的频率，则以上推断应为正确的。 如右图为t-Vx\Vy图， 可见，推测vx(t)和vy(t)为有固定相位差的正弦函数。 如右图为Vx-Vy相图： 粒子相图为椭圆，则粒子做速度周期正弦振动，以上估计正确。 附： Matlab程序： functionimp globalEB E=0.1; B=1; %q/m=1; [t,r]=ode45(@f,[0:0.1:50],[010000]); x=r(:,1); y=r(:,3); z=r(:,5); vx=r(:,2); vy=r(:,4); figure(1) plot(x,y) axis([-1,6,-2,0.5]) gridon xlabel(x) ylabel(y) title(x-y图) figure(2) plot(t,x,t,y) gridon xlabel(t) ylabel(x\y) title(t-x\y图) figure(3) plot(vx,vy) gridon xlabel(vx) ylabel(vy) title(vx-vy图) figure(4) plot(t,vx,t,vy) gridon xlabel(t) ylabel(vx\vy) title(t-vx\vy图) figure(5) X=fft(x,512); Pxx=X.*conj(X)/512; Y=fft(y,512); Pyy=Y.*conj(Y)/512; f=1000/256*(0:255); plot(f,Pxx(1:256),f,Pyy(1:256)) title(FFT-x\y) xlabel(Frequency(Hz(x\y))) axis([0,50,0,1000]) end functionrdot=f(t,r) globalEB rdot=[ r(2); B.*r(4); r(4); E-B.*r(2); r(6); 0 ]; end