電磁數(shù)值計(jì)算作業(yè)4玉海燕

1、題四：用矩量法（拋物線(xiàn)法）分析導(dǎo)電圓環(huán)的電容變化 玉海燕（Y002090833）一：已知條件和計(jì)算目標(biāo) 已知導(dǎo)電圓環(huán)，半徑是3m ,導(dǎo)線(xiàn)半徑是0.05m 用矩量法分析計(jì)算圓環(huán)的電容的變化和線(xiàn)電荷密度（C=1.7*e-10F）二：圓環(huán)結(jié)構(gòu)示意圖圖一：圓環(huán)的結(jié)構(gòu)示意圖圖二：電容隨分段數(shù)的變化曲線(xiàn)圖三，電容隨圓環(huán)半徑的變化曲線(xiàn)圖四：電容隨導(dǎo)線(xiàn)半徑的變化曲線(xiàn)圖五： 電荷密度分布四：結(jié)果分析分析圖二，隨著分段數(shù)的增加，圓環(huán)的電容持續(xù)增加。分析圖三，圓環(huán)電容隨半徑的增加而增大。分析圖四，圓環(huán)電容隨導(dǎo)線(xiàn)半徑的增加呈上升。分析圖五，不管分段數(shù)如何變化，圓環(huán)電荷密度為一條直線(xiàn)五：源程序計(jì)算電容主程序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）電容隨分段數(shù)的變化y=;for k=1:100 y(k,1)=dianrong(k,3,0.05);end;k=1:1:100;plot(k,y,r)grid onxlabel(分段數(shù)N)ylabel(電容值/ F)title(電容隨分段數(shù)變化的曲線(xiàn))（3）電容隨圓環(huán)半徑的變化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( 圓環(huán)半徑/m )ylabel(電容值/F)title(電容隨圓環(huán)半徑的變化曲線(xiàn))（4）電容隨導(dǎo)

4、線(xiàn)半徑的變化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(導(dǎo)線(xiàn)半徑/m)ylabel(電容值/ F )title(電容隨導(dǎo)線(xiàn)半徑變化的曲線(xiàn))（5）圓環(huán)線(xiàn)電荷密度分布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