Matlab 与电磁场模拟一单电荷的场分布：单电荷的外部电位计算公式：4 兀 sr0等位线就是连接距离电荷等距离的点，在图上表示就是一圈一圈的圆，而电力线就是由点向外辐射的线。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);fortheta=[-pi/40pi/4]y=x*tan(theta);holdon;plot(x,y);endgridon单电荷的等位线和电力线分布图：二多个点电荷的电场情况：模拟一对同号点电荷的静电场设有两个同号点电荷，其带电量分别为+Q1 和+Q2(Q1、Q2>0)距离为 2a 则两电荷在点P(x,y)处产生的电势为：q_也甩也,©片 5血估卜 bJa-a)Sb由电场强度可得 E=-VU,在 xOy 平面上，电场强度的公式为:为了简单起见,对电势 U 做如下变换:"现呱+：『+”中氓Matlab 程序：q=1;xm=2.5;ym=2;.\10-2.\同理，将程序稍作修改，便可以得到异号电荷的静电场图像：三、线电荷产生的电位：设电荷均匀分布在从 z=-L 到 z=L,通过原点的线段上，其密度为 q（单位 C/m），求在xy 平面上的电位分布。点电荷产生的电位可表示为 VFQfE^n 是一个标量。其中 r 为电离。线电荷所产生的电位可用积分或叠加的方法来求。为此把线电荷分为 N 段，每段长为dL。每段上电荷为 q*dL，看作集中在中点的点电荷，它产生的电位为然后对全部电荷求和即可。把 xy 平面分成网格，因为 xy 平面上的电位仅取决于离原点的垂直距离 R,所以可以省略一维，只取 R 为自变量。把 R 从 0 到 10 米分成 Nr+1 点，对每一点计算其电位。Matlab 程序：clearall;L=input（'线电荷长度 L='）;N=input（'分段数 N='）;Nr=input（'分段数 Nr='）;q=input（'电荷密度 q='）;29.\0652123466a91fiedHs：11810^—J 口凶igureC:'■■.Pt-ogjraTiFilczSNAllJlBT]kvc-iCormaiidWindowLOe-KJlO♦:9.31990„B654线电荷怅度 L=・E芬段数 N=:50付段炯 r=i50电荷 i?l^q=；:FileEdit\旦 3Zn.£4r1T^ols实=lc+aFHdp□QB；Q'K□S□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);fork=1:Nr+1Rk=sqrt(Lm.A2+R(k)A2);Vk=C0*dL*q./Rk;V(k)=sum(Vk);end[max(V),min(V)]plot(R,V),grid线电荷产生的静电位分布图：.\.\C0=1/(4*pi*E0);foric=1:NFx=0.0;Fy=0.0;forjc=1:Nif(ic~=jc)xij=x(ic)-x(jc);yij=y(ic)-y(jc);Rij=sqrt(xijA2+yijA2);Fx=Fx+CO*q(ic)*q(jc)*xi...
