电磁数值计算作业自己要交的徐红梅

1、作业七1.题目：用矩量法计算圆板电容徐红梅Y0020908142.已知条件：一个很薄的圆盘导体板，其半径为 R。计算目标：用矩量法计算其电容值3.分析对象的结构示意图：图一：圆盘结构示意图4，计算曲线：图二：圆盘电容与等面积方板电容随 N 变化比较曲线分析：从图二的曲线变化趋势可以看出，由此可见两者误差很小。圆盘电容可以由方板近似求得。5.源程序：function Capacitance=disk(R,N)epsino=1/(36*pi*1e9);deltaR=R/N;SN=N2;% the disk is segmentedto SN equal partsdeltaS=eltaR2;x(1

2、)=0;y(1)=0;for m=2:SNRingNum(m)=ceil(sqrt(m);%the ring numberNumOfS=RingNum(m)2-(RingNum(m)-1)2;%total number of sectionshis ringdelheta(m)=2*pi/NumOfS;%the angle of eachsection inthis ringr(m)=(RingNum(m)-0.5)*deltaR;theta(m)=(m-(RingNum(m)-1)2-0.5)*delheta(m);x(m)=r(m)*cos(theta(m);y(m)=r(m)*sheta

3、(m);end%if m=nform=1:SNS(m,m)=0.282*sqrt(deltaS)/epsino;end%ifmnform=1:SNfor n=1:(m-1)S(m,n)=deltaS/(4*pi*epsino*sqrt(x(m)-x(n)2+(y(m)-y(n)2);S(n,m)=S(m,n);endendD=one,1);invS=inv(S);C=invS*D;P=one,1);Q=P*C*deltaS;Capacitance=Q;function Capacitanquareplate(a,N)epsino=1/(36*pi*1e9);sqrtN=sqrt(N);b=a/

4、sqrtN;deltaS=4*b2;S=ones(N);form=1:Nif mod(m,sqrtN)0 x1(m)=mod(m,sqrtN);elsex1(m)=sqrtN;endy1(m)=ceiqrtN);x(m)=(2*x1(m)-1)*b;y(m)=(2*y1(m)-1)*b;end% m=nfor n=1:NS(n,n)=2*b*log(1+sqrt(2)/(pi*epsino);end% m=/nfor m=1:Nfor n=1:(m-1)S(m,n)=b2/(pi*epsino*sqrt(x(m)-x(n)2+(y(m)-y(n)2);S(n,m)=S(m,n);endendS;D=ones(N,1);invS=inv(S);C=invS*D;P=ones(N,1);Q=P*C*deltaS;Capacitance=Q;主函数clear;clf;for n=3:30diskC(n-2)=disk(2/sqrt(pi),n);squareplateC(n-2)=squareplate(1,n2);endn1=3:30;plot(n1,diskC,b-d)hold