带电粒子在电磁场中运动的MATLAB仿真Word下载.docx
《带电粒子在电磁场中运动的MATLAB仿真Word下载.docx》由会员分享，可在线阅读，更多相关《带电粒子在电磁场中运动的MATLAB仿真Word下载.docx（20页珍藏版）》请在冰点文库上搜索。
help
帮助
linspace
区间等分
lookfor
查找
length
获取数组长度
type
列程序清单
size
矩阵大小
which
查找文件目录
max
最大值
double
双精度
min
最小值
str2num
字符串转化为数值
sum
求和
num2str
数值转化为字符串
2.2.2数据和变量
1.表达式
在命令窗口作一些简单的计算，就如同使用一个功能强大的计算器，使用变量无须预先定义类型。
几个表达式可以写在一行，用分号（；
）或逗号（，）分割，用分号（；
）使该表达式运算结果不显示，而逗号（，）则显示结果。
也可以将一个长表达式分在几行上写，用三点（…）续行。
若需要修改已执行过的命令行，可以在命令历史中找到该命令行复制，再粘贴至命令窗口修改。
也可以直接使用键盘↑↓调出已执行过的命令行修改。
2.数据显示格式
MATLAB默认的数据显示格式为短格式（short）：
当结果为整数，就作为整数显示；
当结果是实数，以小数点后四位的长度显示。
若结果的有效数字超出一定范围，以科学计数法显示（如3.2000e-006表示
）。
数据显示格式可使用命令Format改变。
3.预定义变量
MATLAB有一些预定义变量（表2-2），启动时就已赋值，可以直接使用，如前我们使用的圆周率pi和虚数单位i。
表2-2常用预定义变量
变量名
说明
i或j
虚数单位
pi
圆周率3.14159…
eps
浮点数识别精度2^（-52）=
realmin
最小正实数
realmax
最大正实数
inf
无穷大
NaN
没有意义的数
4.用户变量
MATLAB变量名总以字母开头，以字母、数字或下划线组成，区分大小写，有效字符长度为63个。
如A，a，a1，a_b都是合法的，且a与A表示不同变量。
在CommandWindow中使用的变量一旦被赋值，就会携带这个值存在于工作空间，直到被清除或被赋予新的值。
2.2.3数组和矩阵运算
MATLAB基本数据单元是无需指定维数的数组。
数组运算是MATLAB最鲜明的特点，一方面可以使得计算程序简明易读，另一方面可以提高计算速度。
1.数组的输入
最常用的数组是双精度数值数组（doublearray）。
一维数组相当于向量，二维数组相当于矩阵，一维数组可以视为二维数组的特例。
二维数组的第一维称为“行”，第二维称为“列”。
MATLAB数组无需预先定义维数。
直接输入数组的元素，用中括号（[]）表示一个数组，同行元素间用空格或逗号分隔，不同行间用分号或回车分隔。
2.数组运算
数组运算是指数组对应元素之间的运算，也称点运算。
矩阵的乘法、乘方和除法有特殊的数学含义，并不是数组对应元素的运算，所数组乘法、乘方和除法的运算符前特别加了一个点。
特别要区分数组运算在乘法、乘方和除法上的意义和表示上与矩阵运算的不同[3]。
表2-3数组运算符
运算
符号
说明
数组加与减
A+B与A-B
对应元素之间加减
数乘数组
k*A或A*k
k乘A的每个元素
数与数组加减
k+A或k-A
k加（减）A的每个元素
数组乘数组
A.*B
点运算只有点乘、点乘方、点除三个，表示对应元素之间的运算；
（.*）是一个整体，点（.）不能漏掉，（.）和（*）之间也不能有空格
数组乘方
A.^k,k.^A
数除以数组
k./A
数组除法
左除A.\B,右除B./A
3.数学函数
数组的数学函数也是按每个元素的运算，使用通常的函数符号，常用数学函数见表2-4。
表2-4数学函数
函数
意义
sin
正弦
fix
向0取整
cos
余弦
mod
模余
tan
正切
rem
除法余数
cot
余切
abs
绝对值（模）
asin
反正弦
exp
指数函数
acos
反余弦
log
自然对数
sqrt
开方
log10
以10为底的对数
2.3利用MATLAB作图
1.常用作图命令和函数（表2-5）。
表2-5常用作图命令和函数
含义
plot
基本二维图形
clabel
等高线高度标志
fplot
一元函数图像
grid
格栅
ezplot
画二维曲线的符号命令
hold
图形保持
plot3
空间曲线
axis
定制坐标轴
meshgrid
网格数据生成
view
改变视点
mesh
网面图
subplot
子图
surf
曲面图
figure
新图形窗口
contour
等高线图
clf
清除图形
contour3
三维等高线图
close
关闭图形窗口
title
标题
ylabel
y轴说明
xlabel
x轴说明
zlabel
z轴说明
2.曲线图
plot（x,y）作出以数据（x（i）,y（i））为节点的折线图，其中x,y为同长度的向量
fplot（'
fun'
[a,b]）作出函数fun在区间[a,b]上的函数图，fun可以是M函数主名，也可是字符串
ezplot（fun,lims）绘制字符串fun（可是显函数、隐函数或参数方程）指定的函数
plot3（x,y,z）空间曲线图，其中x，y，z为同长度向量
图形显示在图形窗口。
在图形窗口可以使用File菜单保存（Save）为M文件，导出（Export）为图形文件。
也可利用图形窗口Edit菜单Copyfigure作为图片复制到剪贴板，从而进一步粘贴到Word或其他应用程序中。
图形的线型、标记、颜色均可根据要求设定。
常用的见表2-6。
表2-6图形元素设定
颜色
线型
标记
b
蓝（默认）
-
实线（默认）
无标记（默认）
g
绿
：
虚线
.
点
r
红
-.
点划线
o
圈
m
洋红
--
划线
x
叉
y
青
+
十字
k
黑
*
星
3.曲面图
[x,y]=meshgrid（xa,ya）当xa,ya分别为m维和n维行向量，得到x和y均为n行m列矩阵。
meshgrid常用于生成X-Y平面上的网格数据
mesh（x,y,z）绘制网面图，是最基本的曲面图形命令，其中x,y,z是同阶矩阵，表示曲面三维数据
surf（x,y,z）绘制曲面图，与mesh用法类似
contour（x,y,z）绘制等高线图，与mesh用法类似
contour3（x,y,z）绘制三维等高线图，与mesh用法类似
4.图形说明和定制
title（'
字符串'
）图形标题说明
xlabel,ylabel,zlabel用法类似于title，分别说明坐标轴x,y,z
gridon/off显示/不显示格栅
boxon/off使用/不使用坐标框
axis（[a,b,c,d]）定制二维坐标轴范围a<
x<
b,c<
y<
d
axis（[a,b,c,d,e,f]）定制三维坐标轴范围a<
d,e<
z<
f
figure开一个新图形窗口
close关闭现有图形窗口
legend（str1,str2,……）图例，字符串str1,str2,……依次为各图形对象说明
第3章实验原理及仿真结果分析
3.1带电粒子在电磁场中运动的原理
带电粒子在磁场中运动时会受到洛伦兹力的作用，且随着初始运动方向和磁
场分布的不同，其运动轨迹会发生不同的变化。
由洛伦兹力的推导公式可知，它垂直于粒子的运动速度，不对运动粒子作功，只改变其运动方向，其大小为：
（3—1）
因此，综合牛顿运动定律就可以精确确定带电粒子在磁场中的运动轨迹。
3.2质量较大的带电微粒在复合场中的运动
这里我们只研究垂直射入磁场的带电微粒在垂直磁场的平面内的运动，并分几种情况进行讨论。
（1）只受重力和洛仑兹力：
此种情况下，要使微粒在垂直磁场的平面内运动，磁场方向必须是水平的。
微粒所受的合外力就是重力与洛仑兹力的合力。
在此合力作用下，微粒不可能再做匀速圆周运动，也不可能做与抛体运动类似的运动。
在合外力不等于零的情况下微粒将做一般曲线运动，其运动加速度遵从牛顿第二定律；
在合外力等于零的情况下，微粒将做匀速直线运动。
无论微粒在垂直匀强磁场的平面内如何运动，由于洛仑兹力不做功，只有重力做功，因此微粒的机械能守恒，即：
（3—2）
（2）微粒受有重力、电场力和洛仑兹力：
此种情况下。
要使微粒在垂直磁场的平面内运动，匀强磁场若沿水平方向，则所加的匀强电场必须与磁场方向垂直。
在上述复合场中，带电微粒受重力、电场力和洛仑兹力。
这三种力的矢量和即是微粒所受的合外力，其运动加速度遵从牛顿第二定律。
如果微粒所受的重力与电场力相抵消，微粒相当于只受洛仑兹力，微粒将以洛仑兹力为向心力，以射入时的速率做匀速圆周运动。
若重力与电场力不相抵，微粒不可能再做匀速圆周运动，也不可能做与抛体运动类似的运动，而只能做一般曲线运动。
如果微粒所受的合外力为零，即所受的三种力平衡，微粒将做匀速直线运动。
无论微粒在复合场中如何运动，洛仑兹力对微粒不做功。
若只有重力对微粒做功，则微粒的机械能守恒；
若只有电场力对微粒做功，则微粒的电势能和动能的总和守恒；
若重力和电场力都对微粒做功，则微粒的电势能与机械能的总和守恒，用公式表示为：
（3—3）
3.3带电粒子垂直射入E和B正交的叠加场的运动分析
带电粒子垂直射入E和B正交的叠加场常见的模型就是速度旋转器。
如图3-1（a）所示，粒子受力特点——电场力F与洛仑兹力f方向相反。
粒子匀速通过速度选择器的条件——带电粒子从小孔S1水平射入，匀速通过叠加场，并从小孔S2水平射出，从不同角度看有三种等效条件：
从力的角度—电场力与洛仑兹力平衡，即：
；
从速度角度——v0的大小等于E与B的比值，即：
从功的角度——电场力对粒子不做功，即：
两个重要的功能关系——当粒子进入速度选择器时速度
，粒子将因侧移而不能通过选择器。
如图3-1（b），设在电场方向侧移
后粒子速度为v，当
时：
粒子向f方向侧移，F做负功—粒子动能减少电势能增加，有：
当
时，粒子向F方向侧移，F做正功——粒子动能增加，电势能减少，有：
（a）（b）
图3-1速度选择器原理图
3.4实验内容
1、1．用Matlab数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动，即
带电粒子进入磁场的方向与磁场方向的角度
（0<
<
90︒）。
（1）带电粒子质量为m,电量为q（q>
0）的带电粒子在均匀稳定的电磁场中的
运动微分方程为：
（3-4）
以电磁场中某点为原点，以E为Oy方向，B为Oz方向建立坐标系O-xyz。
由于
，则上式的投影方程为：
（3-5）
将其转换为一阶微分线性方程组，以便用差分替代微分作数值计算，令
，
，则上式化简为：
（3-6）
下面使用MATLABR2009b编写程序，分三种情况考虑:
①电场强度和磁场强度都不为零；
②电场强度为零,磁场强度不为零；
③电场强度不为零,磁场强度为零。
程序见附录1。
运行程序得到该带电粒子在均匀稳定的电磁场中的运动轨迹如图3-2所示：
图3-2带电粒子在均匀稳定的电磁场中的运动轨迹
（2）当v与B有一定夹角
时，可将v分解为：
（即
若只有分量
，带电粒子将在垂直的平面内作匀速圆周运动，若只有分量
，粒子将沿磁场B方向作匀速直线运动，当两分量同时存在时，带电粒子的运动轨迹将是一条螺旋线。
用MATLABR2009b编写程序，源程序见附录2.
运行程序，得到带电粒子的运动轨迹如图3-3所示：
图3-3带电粒子的运动轨迹（螺旋线）
第4章MATLAB仿真的应用
4.1用MATLAB仿真带电粒子在磁场中运动的优点
由于电磁场看不见，摸不着，它不像普通的“三态”物质那样由原子、分子
构成，也没有可见的形态，但其具有可以被检测的运动速度、能量和动量，占有
空间，是一种真实的客观存在。
实验中通过MATLAB数值模拟模拟带电粒子在恒定磁场中的螺旋运动和磁聚焦现象，让我们对带电粒子在磁场中的运动有了更深的感性认识，用MATLAB画出的立体图也更有利于对此的理解，对于对应知识的理解和吸收有很大的帮助。
运用数值模拟方法作出的螺旋线运动和磁聚焦现象的轨迹不仅生动形象,而且结果是建立在定量计算的基础上的,不是简单的定性的模拟,因此更具有科学性与说服力。
4.2用MATLAB仿真在物理实验中的应用
MATLAB现已成为集数值分析、符号运算、图形处理、系统仿真等功能于一体的科学与工程计算软件平台。
MATLAB的最大特点在于功能强大而界面友好,其丰富的库函数和各种专用工具箱,将使用者从繁琐的底层编程中解放出来,使他们有更多的时间和精力去探究科学问题本身；
它对科学计算结果迅捷而准确的可视化能力,有助于使用者化抽象思维为形象思维,从而更好地洞察含义、理解概念、发现规律。
目前,MATLAB已成为工程界最流行的软件工具,在大学理工实验课程教学中的应用亦渐成热点[4]。
在构建物理模型的前提下,利用MATLAB的计算、声音、图像、动画等功能来模拟物理现象,诠释物理规律,演示运动过程,并由此构建学生自主探索问题和解决问题的平台。
利用它开展大学物理计算机模拟实验,一方面使抽象的问题形象化,加深学生对物理现象和规律的理解,提高教学效率,另一方面则有利于激发学的学习兴趣,发展学生的想象力,开发学生的创造力[5]。
在物理实验教学过程中,教师可利用MATLAB绘图和动画功能,选择其中合适的部分用在自己的讲解和演示文稿及多媒体课件中,以便说明讲解的结构,形象地演示物理教学中某些难以理解的内容;
或用模拟的图表、动画等展示动态的变化过程和理论模型等,帮助学生了解和理解所学的知识。
形象而生动,非常便于学生接受。
在信息技术、计算机辅助教学迅速发展的时代,上述工作显示了比较大的发展空间和潜力,已经广泛地引起了学生的各方面兴趣,正在成为学生自我学习的平台,已显示出在提高学生素质、提高学生学习质量以及在教学方法和教学内容改革等方面所具有的优势。
第5章总结
基于MATLAB计算机物理模拟方法实现的物理图像和物理过程,教学意图明确,操作简单,人机交互性强,交互参数调节方便,响应速度快,具有很高的容错能力。
在物理实验教学中,以MATLAB为实验演示平台,对物理实验现象进行模拟,作为演示配合物理理论的讲授,能使学生加深对物理现象和规律的理解,提高教学效率。
MATLAB绘图和动画功能,能将计算结果或外部文件数据以二维、三维图形呈现而且提供了方便的图形导出功能,其图形可以jpg、tif等多种格式导出,在物理实验教学中开发利用,可作为物理实验CAI课件的素材图[6]。
通过这次试验，进一步学习和掌握了MATLAB的使用，为今后的学习和应用MATLAB打下了更好的基础。
另外通过这次使用MATLAB仿真带电粒子在不同条件下在磁场中的运动，深入且形象的认识了带电粒子在磁场中的运动方式，而不在限于之前的理论认识。
同时，在实验时也遇到了不少的问题。
1、还是像上次实验一样不能熟练使用MATLAB软件，导致花费时间较多。
2、模拟效果感觉不是很好。
希望通过这次的学习及今后的进一步学习使用MATLAB，能更好的使用MATLAB解决更多的问题。
参考文献
[1]StevenE.Koonin，秦克诚译.计算物理学.北京：
高等教育出版社，1993.
[2]马文淦等。
计算物理学.合肥：
中国科学技术大学出版社，1992.
[3]张志涌。
精通MATLAB6.5.北京：
北京航空航天大学出版社，2003.
[4]罗志恒,陆建隆.MATLAB在普通物理教学中的应用[J].物理通报,2003.
[5]梁齐,许晓琳、王明虎.光干涉实验的计算机模拟[J].大学物理验,2003.
[6]张星辉,钟晓春,张世昌.大学物理教学中引用Matlab图像和动画的实践与思考[J].贵州大报,2004.
附录
附录1
%用Matlab数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动
functionlx
globalqmbe
q=1.6e-2;
m=0.02;
b=[2;
1;
0];
e=[1;
0;
1];
strd{1}='
e\neq0,b\neq0'
;
strd{2}='
e=0,b\neq0'
strd{3}='
e\neq0,b=0'
fori=1:
3
[t,w]=ode23（@ddlzfun,[0:
0.01:
20],[0,0.01,0,6,0,0.01],[],q,m,b（i）,e（i））;
axes（'
unit'
'
normalized'
position'
[0.045+（i-1）*0.350.0620.27860.6583]）;
plot3（w（:
1）,w（:
3）,w（:
5）,'
linewidth'
2）;
gridon
title（strd{i},'
fontsize'
12,'
fontweight'
demi'
）;
xlabel（'
x'
ylabel（'
y'
zlabel（'
z'
view（[-51,18]）;
end
functionwdot=ddlzfun（t,w,q,m,b,e）%该函数实现（3）式的差分迭代
symsqmbe
wdot=[w
（2）;
q.*b.*w（4）./m;
w（4）;
q.*e./m-q.*b.*w
（2）./m;
w（6）;
0]
附录2：
globalqmB
[t,y]=ode45（@ddlzfun,[0:
20],[0,0.01,0,6,0,0.01],[],q,m,B,0）;
%用ode23解微分方程组，时间设为20s
%指定初始条件，传递相关参数
comet3（y（:
1）,y（:
3）,y（:
5））;
plot3（y（:
gridon%开启坐标网格线
模拟带电粒子在恒定磁场中的螺旋运动'
functionydot=ddlzfun（t,y,q,m,B,E）%q,m,B,E为参量
symsqmB
B=2;
ydot=[y
（2）;
q*B*y（4）/m;
y（4）;
-q*B*y
（2）/m;
y（6）;
0]
升级会员 