电磁数值计算作业4玉海燕

1、题四：用矩量法（抛物线法）分析导电圆环的电容变化 玉海燕（Y002090833）一：已知条件和计算目标 已知导电圆环，半径是3m ,导线半径是0.05m 用矩量法分析计算圆环的电容的变化和线电荷密度（C=1.7*e-10F）二：圆环结构示意图图一：圆环的结构示意图图二：电容随分段数的变化曲线图三，电容随圆环半径的变化曲线图四：电容随导线半径的变化曲线图五： 电荷密度分布四：结果分析分析图二，随着分段数的增加，圆环的电容持续增加。分析图三，圆环电容随半径的增加而增大。分析图四，圆环电容随导线半径的增加呈上升。分析图五，不管分段数如何变化，圆环电荷密度为一条直线五：源程序计算电容主程序functi

2、onC=dianrong(N,b,r)s=zeros();d=ones(N,1);E0=(1/36/pi)*1e-9;k=0;t=2*pi/N;dz=t*b;for i= 1:N for j=1:N if i=j s(i,j)=1/4/pi/E0*log(dz/2+sqrt(dz/2)2+r2)/(-dz/2+sqrt(dz/2)2+r2); else x=cos(t*(i-1)-cos(t*(j-1); y=sin(t*(i-1)-sin(t*(j-1); R=b*sqrt(x2+y2); s(i,j)=11*dz/(12*R); end endend; a=4*pi*E0*inv(s)*d

3、; k=sum(a); C=k*dz;（2）电容随分段数的变化y=;for k=1:100 y(k,1)=dianrong(k,3,0.05);end;k=1:1:100;plot(k,y,r)grid onxlabel(分段数N)ylabel(电容值/ F)title(电容随分段数变化的曲线)（3）电容随圆环半径的变化y=;for k=1:300 b=k/100; y(k,1)=dianrong(100,b,0.05);end;k=1:1:300;plot(k,y,r)grid onxlabel( 圆环半径/m )ylabel(电容值/F)title(电容随圆环半径的变化曲线)（4）电容随导

4、线半径的变化y=;for k=1:300 n=k/100; y(k,1)=dianrong(100,3,n);end;n=0.01:0.01:3;plot(n,y,r)grid onxlabel(导线半径/m)ylabel(电容值/ F )title(电容随导线半径变化的曲线)（5）圆环线电荷密度分布functiona=dianrong2(N,b,r)s=zeros();d=ones(N,1);E0=(1/36/pi)*1e-9;k=0;t=2*pi/N;dz=t*b;for i= 1:N for j=1:N if i=j s(i,j)=1/4/pi/E0*log(dz/2+sqrt(dz/2)2+r2)/(-dz/2+sqrt(dz/2)2+r2); else p=cos(t*(i-1)-cos(t*(j-1); q=sin(t*(i-1)-sin(t*(j-1); R=b*sqrt(p2+q2); s(i,j)=11*dz/(12*R); end endend; a=4*pi*E0*inv(s)*d;y=zeros();y=dianr