电磁场的Matlab仿真.pdf

.\ Matlab 与电磁场模拟 一 单电荷的场分布： 单电荷的外部电位计算公式： q 4 0 r 等位线就是连接距离电荷等距离的点， 在图上表示就是一圈一圈的圆， 而电力线就是由点向 外辐射的线。 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); 同号电荷的静电场图像为： 1 1.5 2 50 2.5