带电粒子在均匀电磁场运动.docVIP

实验五 带电粒子在均匀电磁场中运动 实验目的 掌握m文件，进行编程绘制带电粒子在均匀电磁场中的运动轨迹图； 掌握利用函数ode23( ) 求解微分方程组； 掌握三维曲面绘图和彗星轨迹图的语句。 实验原理 （一）三维曲线的作图命令 三维曲线的作图命令 plot3(x,y,z,s) 其中x，y，z是同维数的矢量或矩阵，每组x，y，z构成一个点的坐标，各点依次相连，形成一条曲线。如果是矩阵，则它们相应的列构成一条三维曲线的数据点坐标，所以用矩阵可以同时画多条空间曲线。S是线型、颜色和标志的参数。 （二）彗星轨迹图命令 彗星轨迹图可以绘出运动的轨迹，其使用格式为： comet3(z); comet3(x,y,z); 三维空间彗星轨迹，x，y，z的含义与plot3命令中的x，y，z的含义，绘出的轨迹慧长默认为0.1 。 （三）解常微分方程组 MATLAB只能解一阶的常微分方程组，高价的常微分方程需要转化成一阶方程组才能求解。对于二阶常微分方程，首先需要化成显式形式，然后令，则二阶常微分方程化为两个一阶常微分方程组成的方程组，从而使问题得到解决。 （四）建立带电粒子在均匀电磁场中的运动方程 1 带电粒子在均匀稳定电磁场中受力分析: 2 带电粒子在均匀稳定电磁场中的运动微分方程为: ① 可将上式分解在直角坐标系展成标量式： ② 令 则化简为: ③ 令 则得出程序认可方程： ④ 3． 设带电粒子质量为m ，带电量为q ，电场强度E沿y方向，磁感应强度B沿z方向，则带电粒子在均匀电磁场中的运动微分方程为 令，则上面微分方程可化作： 实验仪器 序 号 名 称 主要用途 1 一台安装Windows2000的pc机。计算机的具体要求：⒈Pentium3以上的CPU；⒉建议至少256MB的内存；⒊建议硬盘至少20GB4.安装MATLAB仿真软件。 运行 MATLAB仿真软件。 实验内容 利用MATLAB的函数，编写m文件，用函数comet3(x,y,z)观察轨迹；用函数plot3(x,y,z)绘出图形。 附代码 BX=input(‘请输入BX=’); BY=input(‘请输入BY=’); BZ=input(‘请输入BZ=’); EX=input(‘请输入EX=’); EY=input(‘请输入EY=’); EZ=input(‘请输入EZ=’); q=1.6e-2; m=0.02; [t,y]=ode23(‘fun1’,[0:0.01:20],[0,0.6,0,6,0,6],[],q); figure(1); plot3(y(:,1),y(:,3),y(:,5),’linewidth’,2); grid on xlabel(‘x’); ylabel(‘y’); zlabel(‘z’); figure(2); comet3(y(:,1),y(:,3),y(:,5)); q=1.6e-2; m=0.02; B=3; E=0; [t,y]=ode(‘fun’,[0:0.2:50],[0,0.01,0,6,0,0.01],[],q,m,B,E); figure(1); plot3(y(:,1),y(:,3),y(:,5),’linewidth’,2); grid on xlabel(‘x’); ylabel(‘y’); zlabel(‘z’); figure(2); comet3(y(:,1),y(:,3),y(:,5)); function YMA=fun(t,y,flag,q,m,B,E) YMA=[y(2); q*B*y(4)/m; y(4); q*E/m-q*B*y(2)/m; y(6); 0];