MATLAB语言课程论文 基于MATLAB的电磁场数值图像分析文档格式.docx
- 文档编号:22232109
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:16
- 大小:245.71KB
MATLAB语言课程论文 基于MATLAB的电磁场数值图像分析文档格式.docx
《MATLAB语言课程论文 基于MATLAB的电磁场数值图像分析文档格式.docx》由会员分享,可在线阅读,更多相关《MATLAB语言课程论文 基于MATLAB的电磁场数值图像分析文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
根据麦克斯韦方程组和结构方程组!
在静电场中,以电位为求解对象,在各向同性介质中,电位满足泊松方程:
(4)
在恒定磁场中,取矢量磁位为求解对象,令有
(5)
考虑电磁波动方程:
(6)
在正弦稳态条件下,上式可分别导出亥姆霍兹方程:
(7)
如上述几个例子,对于不同的电磁场实际应用问题求解可以得到对应的电磁场偏微分方程,直接用解析法求解这些方程组往往会遇到很多困难甚至无法求解,电磁场数值分析方法已成为求解电磁场问题的重要方法。
数值分析方法将原来连续的场域离散化求解,再用离散点的结果近似逼近连续场域的解&
常用的电磁场数值分析方法包括有限差分法、边界元法和有限元法。
有限差分法是以差分原理为基础的一种数值计算方法,即用差分方程代替偏微分方程,把要求解的边值问题转化为一组相应的差分问题,将求解区域划分,求解差分方程组从而得出各网格单元的场值。
这种方法的特点是方法简单,网格划分容易,但对不规则边界处理不便,网格划分缺乏灵活性。
边界元法以麦克斯韦积分方程为基础,它采用分步积分如格林定理等在一定条件下把该积分方程转化为关于边界的积分方程,并据此进行离散获得对应的代数方程!
求出场域中变量的数值。
它的特点是将数值法与解析法相结合,在数学上起到降维的作用,减少了计算量,但对非线性情况失去了高精度特点,有局限性。
有限元法由于单元定义灵活,处理边界条件容易,具有正定对称系数矩阵而占据主导地位&
有限元法是根据变分原理和离散化而取得近似解的方法。
是先从偏微分方程边值问题出发,找出一个能量泛函的积分式,并令其在满足第一类边界条件的前提下取极值,即构成条件变分问题&
例如与上式对应的二维泊松场第三类边界条件下的泛函极值问题为:
(8)
与上式对应的亥姆霍兹方程泛函为:
(9)
这个条件变分问题是和偏微分方程边值问题等价的。
有限元法便是以条件变分问题为对象来求解电磁场问题。
在求解过程中,将场的求解区域剖分成有限个单元,因此在单元中构造出插值函数,将插值函数代入能量泛函的积分式,再把泛函离散化成多元函数。
通过多元函数极值求极值方法得到一个代数方程组。
最后由此方程组求解得到数值解。
下面举几例论证:
例证1
设圆线圈的中心为O,半径为R,放置于y-z平面,线圈通过的电流为I0,如右图所示。
用毕奥-萨伐尔定律计算载流圆线圈在z=0处x-y平面上的磁场分布。
解题分析
根据毕奥-萨伐尔定律,
(10)
线圈上任一点处的电流元在x-y平面上一点P产生的元磁场为dB。
在编制程序时,将电流环分为N段,每一小段视为一电流元,然后求出每一电流元在观察点处的磁场分量,求出总磁场,最后叠加。
MATLAB程序如下:
clearall%清除
R=input('
请输入圆环半径,R='
);
%定义输入变量
I0=input('
请输入电流,I0='
%定义输出变量
mu0=4*pi*1e-7;
C0=mu0/(4*pi);
%归并常数
N=20;
%电流环分段数
x=linspace(-3,3,N);
y=x;
%确定观测点范围
theta0=linspace(0,2*pi,N+1);
%环的圆周角分段
theta1=theta0(1:
N);
y1=R*cos(theta1);
z1=R*sin(theta1);
%环各段矢量的起始坐标y1,z1
theta2=theta0(2:
N+1);
y2=R*cos(theta2);
z2=R*sin(theta2);
%环各段矢量的终点坐标y2,z2
xc=0;
yc=(y2+y1)./2;
zc=(z2+z1)./2;
%计算环各段矢量中点的三个坐标分量xc,yc,zc
dlx=0;
dly=y2-y1;
dlz=z2-z1;
%计算环各段矢量dl的三个长度分量,其中x1=x2=0。
NGx=N;
NGy=NGx;
%网格线数
fori=1:
NGy%循环计算各网点上的B(x,y)值
forj=1:
NGx
rx=x(j)-xc;
ry=y(i)-yc;
rz=0-zc;
%计算径矢r的3个长度分量,r在z=0的面上。
r3=sqrt(rx.^2+ry.^2+rz.^2).^3;
%计算r3
dlXr_x=dly.*rz-dlz.*ry;
%计算叉乘dl×
r的x和y分量,z分量为0
dlXr_y=dlz.*rx-dlx.*rz;
Bx(i,j)=sum(C0*I0.*dlXr_x./r3);
%把环各段产生的磁场分量累加
By(i,j)=sum(C0*I0.*dlXr_y./r3);
B=(Bx.^2+By.^2).^0.5;
%计算B的大小
end
subplot(1,2,1),quiver(x,y,Bx,By),%画矢量场图
holdon
plot(0,1,'
ro'
0,-1,'
bo'
),
xlabel('
x'
),ylabel('
y'
),%修饰图形,标注坐标轴
axis([-3,3,-3,3]),
subplot(1,2,2)
mesh(x,y,B);
axis([-3,3,-3,3,0,1e-4])%画磁场大小分布图
),zlabel('
B'
)
运行该程序,例如,在命令窗中的R和I0的提示后分别键入1和100,运行结果如图1所示
图1载流圆线圈的磁场分布度及位置曲线
例证2
一对相同的圆形线圈,彼此平行而共轴。
设两线圈内的电流都是I,且回绕方向一致,线圈的半径为R,二者的间距为a(当a=R时,称为亥姆霍兹线圈),求轴线附近的磁场分布。
本题是把观测区域取在两线圈之间的小范围内。
线圈B生成的左边的磁场等于线圈A的左边磁场。
因此,A、B两线圈在中间部分的合成磁场等于A线圈的右磁场与其左磁场平移R后的和。
clearall%清除
C0=mu0/(4*pi);
%归并常数
I0=5.0;
R=1;
NGx=21;
NGy=21;
%设定网格线数
x=linspace(-1,1,NGx);
y=linspace(-1,1,NGy);
%电流环分段数
%环的圆周角分段
%计算环各段矢量dl的三个长度分量,其中x1=x2=0。
yc=(y2+y1)/2;
zc=(z2+z1)/2;
NGy
rx=x(j)-xc;
%计算径矢r的3个长度分量,r在z=0的面上。
r3=sqrt(rx.^2+ry.^2+rz.^2).^3;
dlXr_y=dly.*rx-dlx.*rz;
Bx(i,j)=sum(C0*I0*dlXr_x./r3);
%把环各段产生的磁场分量累加
By(i,j)=sum(C0*I0*dlXr_y./r3);
Bax=Bx(:
11:
21)+Bx(:
1:
11);
%把x<
0区域内的磁场平移,叠加到x>
0区域
Bay=By(:
21)+By(:
subplot(1,2,1),mesh(x(11:
21),y,Bax);
ylabel('
zlabel('
subplot(1,2,2),quiver(x,y,Bx,By,1.5),axis('
square'
),axis([-1,1,-1,1]),
运行结果如图2所示。
可以看出,在轴线附近磁场大小均匀且沿x方向。
图2亥姆霍兹线圈轴线附近Bx在x-y平面上的分布及矢量场
例证3
设带电粒子质量为m,带电量为,电场强度E沿方向,磁感应强度B沿qyz方向.则带电粒子在均匀电磁场中的运动微分方程为
(11)
令,则上面微分方程可化作:
(12)
选择和为参量,就可以分别研究0≠E,0=B和,等情况.编写MATLAB程序如下:
clear%清除
symswxyztBEmq;
%定义变量
E=input('
E='
B=input('
B='
%输入E和B值
[x,y,z]=dsolve('
D2x=q*B/m*Dy'
'
D2y=q*E/m-q*B/m*Dx'
D2z=0'
x(0)=0'
y(0)=0'
z(0)=0'
Dx(0)=0.01'
Dy(0)=6'
Dz(0)=0.01'
);
%初始条件取x(0)=y(0)=z(0)=0,Dx(0)=0.01,Dy(0)=6,Dz(0)=0.01
q=1.6e-2;
m=0.02;
%赋值
X=subs([xyz]);
x=X
(1),y=X
(2),z=X(3),
ezplot3(X
(1),X
(2),X(3))%绘图函数调用
运行上述程序,例如,取E=4,B=8可得下列特解并给出图运行结果如图3所示。
研究时可以采用不同的初始条件和不同的参量观察不同的现象。
运算特解如下:
x=
(t*exp
(1))/8-(100*i*exp
(1)-8*i+4800)/(10240*exp((32*i*t)/5))-(exp((32*i*t)/5)*(8*i-100*i*exp
(1)+4800))/10240+15/16
y=
(5*exp
(1))/256+(i*(100*i*exp
(1)-8*i+4800))/(10240*exp((32*i*t)/5))-(i*exp((32*i*t)/5)*(8*i-100*i*exp
(1)+4800))/10240-1/640
z=
t/100
图3现有E=4,B=8参数运行结果
例证4
如右图所示,求垂直于无限长载流直导线的平面内磁感应强度的分布。
设场点P的位置为
,电流元位置为
,电流元矢量为
。
由此,场点P相对于电流元的位置矢量为
利用行列式计算Idl×
r,可写为
(13)
也可利用MATLAB中的det命令函数来求该行列式,MATLAB程序如下:
symsdxdydzx0xy0yz0z;
dl=[dx,dy,dz];
%定义行列式
r=[x0-x,y0-y,z0-z];
d1cr=cross(dl,r)%求dl×
r的积
运行结果为
d1cr=[dy*(z0-z)-dz*(y0-y),dz*(x0-x)-dx*(z0-z),dx*(y0-y)-dy*(x0-x)]
即
又,r的大小为
设载流导体通过坐标原点垂直于x-y平面放置,电流元Idl沿z轴正向,场点P位于x-y平面上。
对本题目而言,dx=dy=0,x=y=0,z0=0,矢量叉乘积为
,r的大小为
由毕奥-萨伐尔定律
(14)
有
;
(15)
MATLAB程序
1、用符号运算求B的表达式
symsC0Izxyrr0;
Bx=C0.*I.*int(-y./(x.^2+y.^2+z.^2).^(3/2),'
z'
-inf,inf)
By=C0.*I.*int(x./(x.^2+y.^2+z.^2).^(3/2),'
B=(Bx.^2+By.^2).^0.5
运行结果:
Bx=-2*C0*I*y/(x^2+y^2)^(3/2)/(1/(x^2+y^2))^(1/2)
By=2*C0*I*x/(x^2+y^2)^(3/2)/(1/(x^2+y^2))^(1/2)
B=(4*C0^2*I^2*y^2/(x^2+y^2)^2+4*C0^2*I^2*x^2/(x^2+y^2)^2)^(1/2)
即
(16)
2、绘制磁场大小分布图和矢量场图
x=-0.5:
0.05:
0.5;
y=x;
I=input('
请输入电流I='
%设置输入
%设置函数
[X,Y]=meshgrid(x,y);
Bx=-2.*C0.*I.*Y./(X.^2+Y.^2).^(3./2.)./(1./(X.^2+Y.^2)).^(1./2);
By=2.*C0.*I.*X./(X.^2+Y.^2).^(3./2)./(1./(X.^2+Y.^2)).^(1./2);
B=(4.*C0.^2.*I.^2.*Y.^2./(X.^2+Y.^2).^2+4.*C0.^2.*I.^2.*X.^2./(X.^2+Y.^2).^2).^(1./2);
subplot(1,2,1)%选择1×
2区域中的1号区
quiver(X,Y,Bx,By,2),axis([-0.5,0.5,-0.5,0.5]),axis('
subplot(1,2,2)%选择1×
2区域中的2号区
mesh(X,Y,B)%绘图
运行该程序,在命令窗中的提示后键入I0值(例如,取I0=100A),便得到图4所示图形。
图4长载流导线的磁场在x-y平面上的分布
三、结论
从以上利用MATLAB语言对几种电磁场模型的分析我们不难的出以下结论:
电磁场与电磁波理论作为电子信息类专业的一门重要基础课程,该课程的特点
是:
一是理论性强,课程中涉及许多高等数学方面的计算,特别是需要应用矢量分析和微分方程,致使电磁场问题的求解十分复杂;
二是理论模型抽象,各种电磁场与电磁波现象都具有复杂的空间分布,而电磁场与电磁波又都具有不可见和不可触摸的特性,只能进行抽象的想象或通过仪器进行数据测量;
三是电磁波是动态的,电磁波是电磁场相互激发的结果,可以通过麦克斯韦方程和边界条件来定量描述,它在空间的传播每时每刻的位置和状态都在发生变化。
四、课程体会
通过在本次写论文的过程中,我受益匪浅。
在我选择MATLAB论文的题目的过程中,我通过查找大量的资料,发现MATLAB涉及的范围非常的大,例如在物理学中、电磁学中、线性代数等中的应用,认识到MATLAB对我们是如此的实用。
其次,在初学习的MATLAB的时候,觉得它和我们大一第一学期学的C语言有很多相似之处,在老师讲的过程中,我争取好好的听每一节课,每次在实验课上我也认真的做好每一个实验的题目,每次在程序结果运行出来的时候,我觉得特有成就感,同时我感觉到了MATLAB强大功能。
经过一学期紧张而有序的课程学习,在忙碌之余也得到了颇多的收获。
我深深体会到MATLAB语言相对于同类程序语言更方便更简洁易懂。
利用Matlab可以设计大量仿真实验,把理论学习和仿真实验教学有效结合起来,加深了对理论学习的理解,自己动手利用MATLAB解决电磁场与电磁波一些实际问题,这也是培养了我们综合素质。
[参考文献]
[1]刘卫国.MATLAB程序设计与应用(第二版)[M].北京:
高等教育出版社,2006.
[2]马文蔚.物理学(中册)(第四版)[M],北京:
高等教育出版社,1999.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB语言课程论文 基于MATLAB的电磁场数值图像分析 MATLAB 语言 课程 论文 基于 电磁场 数值 图像 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)