作业陈陈兵辉.pdfVIP

第二题 用矩量法计算圆环 陈兵辉 Y060439 1 题目: 一个金属导线圆环，其圆环半径为 3m，其导线半径为 0.05m ， 求这圆环的电容？ 2 已知条件: 圆环半径b=3m ，导线半径a=0.05m 3 计算目标 圆环的电容c 4 分析对象的结构示意图 5 计算曲线 5.1Pulse expansion plus point matching 圆环半径为3m 导线半径0.05m 圆环半径为3m 导线半径0.05m 1--300 分11 段时 电容c=1.7287 10−10 F 分21 段时 电容c=1.727910−10 F −10 分30 段时 电容c= 1.726910 F 圆环半径20m 导线半径0.01—3m 分100 段 5.2Parabola expansion plus point matching 圆环半径为3m 导线半径0.05m 1--300 分11 段时 电容c=1.826910−10 F 分21 段时 电容c=1.845210−10 F 分30 段时 电容c= 1.852910−10 F 圆环半径20m 导线半径0.01—3m 分100 段 六源程序 源程序 1 function nj=eee(n,r,a) k=zeros(n,n); sc=ones(n,1); qiuhe=0; c=pi/n; gh=log(pi/n*3+sqrt(a^2+(pi/n*3)^2))-log(-pi/n*3+sqrt(a^2+(pi/n*3)^2)); for i=1:n for j=1:n if i==j k(i,j)=gh; else k(i,j)=c/sin(abs(i-j)*c);%1/abs(i -j); end end end niju=inv(k); xc=niju*sc; for h=1:n qiuhe=qiuhe+xc(h,1); end nj=qiuhe*4*2*pi*r/n/36/10^9; 源程序2 function nj=fff(n,r,a) k=zeros(n,n); sc=ones(n,1); qiuhe=0; c=pi/n; lnln=2*c*r; %gh=log(pi/n*3+sqrt(a^2+(pi/n*3)^2))-log(-pi/n*3+sqrt(a^2+(pi/n*3)^2)); gh=(1+a^2/2/lnln^2)*(log(lnln/2+sqrt(lnln^2/4+a^2))-log(-lnln/2+sqrt(lnln^2/4+a^2)))-sqrt(lnln^2 /4+a^2)/2/lnln; for i=1:n for j=1:n if i==j k(i,j)=gh; else k(i,j)=11*c/sin(abs(i-j)*c)/12;%1/abs(i -j); end end end niju=inv(k); xc=niju*sc; for h=1:n qiuhe=qiuhe+xc(h,1); end nj=qiuhe*4*2*pi*r/n/36/10^9; 画图程序 aa=zeros(1,300); bb=zeros(1,300); for ff=1:300 aa(1,ff)=ff; end for gg=1:300 bb(1,gg)=fff(gg,3,0.05); end plot(aa,bb,r) xlabel(parabola