原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Matlab与电磁场模拟
一 单电荷的场分布：
单电荷的外部电位计算公式：
等位线就是连接距离电荷等距离的点，在图上表示就是一圈一圈的圆，而电力线就是由点向外辐射的线。
MATLAB程序：
theta=[0:.01:2*pi];
r=0:10;
x=sin(theta)*r;
y=cos(theta)*r;
plot(x,y,b)
x=linspace(-5,5,100);
for theta=[-pi/4 0 pi/4]
y=x*tan(theta);
hold on;
plot(x,y);
end
grid on
单电荷的等位线和电力线分布图：
二 多个点电荷的电场情况：
模拟一对同号点电荷的静电场
设有两个同号点电荷,其带电量分别为 +Q1和+Q2(Q1、Q20 )距离为 2a则两电荷在点P(x, y)处产生的电势为:
由电场强度可得E = -?U,在xOy平面上,电场强度的公式为:
为了简单起见,对电势U做如下变换:
。
Matlab程序：
q=1;
xm=2.5;
ym=2;
x=linspace(-xm,xm);
y=linspace(-ym,ym);
[X,Y]=meshgrid(x,y);
R1=sqrt((X+1).^2+Y.^2);
R2=sqrt((X-1).^2+Y.^2);
U=1./R1+q./R2;
u=1:0.5:4;
figure
contour(X,Y,U,u)
grid on
legend(num2str(u))
hold on
plot([-xm;xm],[0;0])
plot([0;0],[-ym;ym])
plot(-1,0,o,MarkerSize,12)
plot(1,0,o,MarkerSize,12)
[DX,DY] = gradient(U);
quiver(X,Y,-DX,-DY);
surf(X,Y,U);
同号电荷的静电场图像为：
同理，将程序稍作修改，便可以得到异号电荷的静电场图像：
三、线电荷产生的电位：
设电荷均匀分布在从z=-L到z=L,通过原点的线段上，其密度为q(单位C/m),求在xy平面上的电位分布。
点电荷产生的电位可表示为 是一个标量。其中r为电荷到测量点的距离。线电荷所产生的电位可用积分或叠加的方法来求。为此把线电荷分为N段，每段长为dL。每段上电荷为q*dL,看作集中在中点的点电荷，它产生的电位为
然后对全部电荷求和即可。
把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;
L0=linspace(-L,L,N+1);
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),grid
线电荷产生的静电位分布图：
四 计算平面上N个电荷之间的库伦引力
1 建模：
由库仑定律：
先输入电荷的数目，各电荷的坐标及电荷量，再选一个电荷，求其它电荷对它的作用力，叠加求合力。再选下一个电荷，依次类推。
Matlab程序：
clear all;
N = input(ê?è?μ?oéêy??N=:);
for ic = 1:N
fprintf(----/n??μ?oé#%g\n,ic);
rc = input(ê?è?μ?oé????[x,y]￡¨?×￡?:);
x(ic) = rc(1);
y(ic) = rc(2);
q(ic) = input(ê?è?μ?oéá?￡¨?a??￡?￡o);
end
E0 = 8.85e-12;
C0 = 1/(4*pi*E0);
for ic = 1:N
Fx = 0.0;Fy = 0.0;
for jc
文档评论（0）