Matlab在电磁场中的应用(2)要点.pdf
《Matlab在电磁场中的应用(2)要点.pdf》由会员分享,可在线阅读,更多相关《Matlab在电磁场中的应用(2)要点.pdf(17页珍藏版)》请在三一文库上搜索。
1、Matlab 在电磁场中的 应用 专业 : 电气信息与自动化 班级:2012 级自动化 3 班 学号:12012242065 学院:物电学院 指导老师:李虹 完成日期:2013 年 12 月 15 日 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬 1 Matlab 在电磁场中的应用 摘要 Matlab 是美国 Mathworks 公司于 80 年代推出的大型数学软件, 通过多年的 升级换代, 现在已发展成为集数值计算、符号计算、 可视化功能以及诸多的工具 箱为一体的大型科学计算软件, 它已广泛应用于科研院所、 工程技术等各个部门, 并成为大学生、研究生必备的工具软件。 电磁学是物理
2、学的一个分支,是研究电场和电磁的相互作用现象。电磁 学从原来互相独立的两门科学( 电学、磁学) 发展成为物理学中一个完整的 分支学科,主要是基于电流的磁效应和变化的磁场的电效应的发现。这两 个实验现象,加上麦克斯韦关于变化电场产生磁场的假设,奠定了电磁学 的整个理论体系,发展了对现代文明起重大影响的电工和电子技术。 针对电磁场学习理论性强、概念抽象等特点,利用Matlab 强大的数值计算 和图形技术,通过具体实例进行仿真,绘制相应的图形,使其形象化,便于对其 的理解和掌握。将Matlab 引入电磁学中,利用其可视化功能对电磁学实验现象 进行计算机模拟,可以提高学习效率于学习积极性,使学习效果明
3、显。 本文通过 Matlab 软件工具,对点电荷电场、线电荷产生的电位、平面上N 个电荷之间的库仑引力、 仿真电荷在变化磁场中的运动等问题分别给出了直观形 象的的仿真图, 形实现了可视化学习, 丰富了学习内容, 提高了对电磁场理论知 识的兴趣。 关键词 :Matlab 电磁学仿真计算机模拟 一、点电荷电场 问题描述: 真空中,两个带正电的点电荷,在电量相同和电量不同情况下的电场分布。 根据电学知识, 若电荷在空间激发的电势分布为V,则电场强度等于电势梯度的 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬 2 负值,即: 根据题意,真空中若以无穷远为电势零点,则在两个点电荷的电场中,
4、空间 的电势分布为: 程序实现: clear all ep0=8.85*1e-12; c0=1/(4*pi*ep0); e=1.60e-10; h=0.018; x=-0.5:h:0.5; y=-0.5:h:0.5; X,Y=meshgrid(x,y); q=e;1.9*e; for i=1:2 V=c0*e./sqrt(X+0.2).2+Y.2)+c0.*q(i)./sqrt(X-0.2).2+Y .2); Ex,Ey=gradient(-V,h); figure(i) contour(X(:,:,1),Y(:,:,1),V,. 20,-20,19,-19,18,-18,17,-17,. 1
5、6,-16,15,-15,14,-14,13,-13,. 12,-12,11,-11,10,-10); axis(-0.38,0.38,-0.28,0.28) hold on phi=0:pi/17:2*pi; sx1=0.2+0.01*cos(phi); sy1=0.01*sin(phi); streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx1,sy1); hold on sx2=-0.2+0.01*cos(phi); sy2=0.01*sin(phi); streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx2,sy2); title(stri)
6、 text(-0.212,0,+,fontsize,20); text(0.187,0,+,fontsize,20); end EV 12 12 0 10 2 44 qq VVV rr 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬 3 图 1-1 两个同号等量电荷的电场分布图 1-2 两个同号不等量电荷的电场分布 二、线电荷产生的电位 设电荷均匀分布在从z=-L 到 z=L,通过原点的线段上,其密度为 q(单位 C/m), 求在 xy 平面上的电位分布。 点电荷产生的电位可表示为 0 /4VQr是一个标量。其中 r 为电荷到 测量点的距离。 线电荷所产生的电位可用积分或叠加的方法来
7、求。为此把线电荷 分为 N 段,每段长为 dL。每段上电荷为 q*dL, 看作集中在中点的点电荷,它产 生的电位为 0 4 qdL dV r 然后对全部电荷求和即可。 把 xy 平面分成网格,因为xy 平面上的电位仅取决于离原点的垂直距离R, 所以可以省略一维,只取R 为自变量。把 R 从 0 到 10 米分成 Nr+1 点,对每一 点计算其电位。 matlab 程序 clear all; L=input(线电荷长度 L:); N=input( 分段数 N:); Nr=input(分段数 Nr:); q=input(电荷密度 q=: ); E0=8.85e-12; C0=1/4/pi/E0;
8、L0=linspace(-L,L,N+1); 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬 4 L1=L0(1:N);L2=L0(2:N+1); Lm=(L1+L2)/2;dL=2*L/N; R=linspace(0,10,Nr+1); for k=1:Nr+1 Rk=sqrt(Lm.2+R(k)2); Vk=C0*dL*q./Rk; V(k)=sum(Vk); end max(V),min(V) plot(R,V),grad 输入: 线电荷长度 L:5 分段数 N:50 分段数 Nr:50 电荷密度 q=:1 可得最大值和最小值为 : ans = 1.0e+010 *9.3199
9、 0.8654 图(2-1)线电荷产生的静电位分布图 三、平面上 N 个电荷之间的库仑引力 建模: 由库仑定律: 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬 5 3 120 / 4Fq qr 其分量的公式可以写成: 3 12210 3 12210 22 2121 ()/ 4 () / 4 ()() x y Fq qxxr Fq qyyr rxxyy 编写程序时,先输入电荷的数目,各电荷的坐标及电荷量,再选一个电荷,求其 它电荷对它的作用力,叠加求合力。再选下一个电荷,依次类推。 Matlab 程序: clear all; N = input(输入电荷数目 N=:); for ic
10、 = 1:N %输入给定条件 fprintf(-/n 对电荷 #%gn,ic); rc = input(输入电荷位置 x,y(米) :); x(ic) = rc(1); %电荷 ic 的 x 坐标 y(ic) = rc(2); %电荷 ic 的 y 坐标 q(ic) = input(输入电荷量(库仑):); end E0 = 8.85e-12; %真空中的常数 C0 = 1/(4*pi*E0); %合并常数 for ic = 1:N %循环计每个电荷所受的力 Fx = 0.0;Fy = 0.0; for jc = 1:N if(ic = jc) xij = x(ic)-x(jc);yij =
11、y(ic)-y(jc); Rij = sqrt(xij2+yij2); Fx = Fx+C0*q(ic)*q(jc)*xij/Rij3; Fy = Fy+C0*q(ic)*q(jc)*yij/Rij3; end end fprintf( 其它电荷作用在电荷 #%g 上的合力为: n,ic); fprintf(x- 分量:%gNn,Fx); fprintf(y- 分量:%gNn,Fy); end 本程序注意学会循环提示并输入参数的方法,以及用双循环解决较复杂的 计算过程的编程问题。 输入已知条件 : 输入电荷数目 N=3 -对电荷 #1 工程物理作业 Matlab 在电磁场中的应用-电力赵洪彬
12、6 输入电荷位置 x,y(m):1 2 输入电荷量(库仑) :2 -对电荷 #2 输入电荷位置 x,y(m):1 1 输入电荷量(库仑) :1 -对电荷 #3 输入电荷位置 x,y(m):3 3 输入电荷量(库仑) :3 计算结果 : 其它电荷作用在# 1 上的合力为 : X-分量为 :-9.65102e+009N Y-分量为 1.31581e+010 其它电荷作用在# 2 上的合力为 : X-分量为 :-2.38431e+009N Y-分量为 -2.03679e+010 其它电荷作用在# 3 上的合力为 : X-分量为 :1.20353e+010N Y-分量为 7.20982e+009 利用
13、 matlab 软件仿真电荷在变化磁场中的运动 程序一 % 电荷在非均匀磁场中的运动 v=10;sita=pi/6; %设定带电粒子的初速度及入射角 v=v*cos(sita); u=v*sin(sita); %计算 x,y 方向的初速度 w=0; t,y = ode23(yy,0:0.002:2,0,v,0,u,0,w); %求解名为“ yy”的微 分方程组 figure %描绘运动轨迹 plot(t,y(:,1); %绘制一般二维曲线 %comet(t,y(:,1); %绘制二维动态曲线 xlabel(t);ylabel(x); figure plot(t,y(:,3); %comet(t
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 电磁场 中的 应用 要点
链接地址:https://www.31doc.com/p-5197503.html