电磁场与电磁波matlab大作业.docx
- 文档编号:117580
- 上传时间:2022-10-03
- 格式:DOCX
- 页数:5
- 大小:615.26KB
电磁场与电磁波matlab大作业.docx
《电磁场与电磁波matlab大作业.docx》由会员分享,可在线阅读,更多相关《电磁场与电磁波matlab大作业.docx(5页珍藏版)》请在冰豆网上搜索。
电磁场与电磁波Matlab大作业
姓名:
张佳辰学号:
13020188003
班级:
1302018提交日期:
2015年4月27日
一、程序概述程序开始时,提示用户输入Q1和Q2两电荷的电荷量,程序读取后,算出Q1和Q2两电荷之比
lamda。
然后根据公式,画出电场线和电势线分布。
二、Matlab程序
clear;
q1=input('请输入Q1的电荷量:
');
q2=input('请输入Q2的电荷量');
lambda=q1./q2;%计算电荷比
xm=4;ym=4;
x=linspace(-xm,xm);
y=linspace(0,ym);
[X,Y]=meshgrid(x,y);
R1=sqrt((X+1).^2+Y.^2);%
第一个电荷到场点的距离
R2=sqrt((X-1).^2+Y.^2);%
第二个电荷到场点的距离
U=2*(log(1./R1)+lambda*log(1./R2));%计算电势
u=-5:
1:
5;
contour(X,Y,U,u,'r','LineWidth',1);%画上半平面等势线
gridon,axisequal,holdon
contour(X,-Y,U,u,'r','LineWidth',1)%画下半平面等势线
plot([-xm;xm],[0;0],[0;0],[-ym;ym],'LineWidth',1)
titleTxt=['µçºÉ±ÈΪ',num2str(lambda),'µÄµç³¡Ïߺ͵çÊÆÏß·Ö²¼£º'];%标题文本
title(titleTxt,'FontSize',12)%显示标题
C=atan((X+1)./Y)+lambda*atan((X-1)./Y);
r0=0.1;
dth=20;
th=dth:
dth:
180-dth;
th=th*pi/180;
x0=r0*cos(th)-1;
y0=r0*sin(th);
c=atan((x0+1)./y0)+lambda*atan((x0-1)./y0);
contour(X,Y,C,c,'g','LineWidth',1)%画上半平面电场线
contour(X,-Y,C,c,'g','LineWidth',1)%画下半平面电场线
iflambda>0%如果电性相同
dth=dth*abs(1/lambda);
th=dth:
dth:
180-dth;
th=th*pi/180;
x0=r0*cos(th)+1;
y0=r0*sin(th);
c=atan((x0+1)./y0)+lambda*atan((x0-1)./y0);
contour(X,Y,C,c,'g','LineWidth',1)%画上半平面电场线
contour(X,-Y,C,c,'g','LineWidth',1)%画下半平面电场线
end
三、程序测试
1)Q1=Q2=1时
2)Q1=1,Q2=-1时
3)Q1=1,Q2=2时
3)Q1=-3,Q2=2时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电磁场 电磁波 matlab 作业