徐紅梅師兄作業(yè)

1、線問題直導(dǎo)線電容值的計算（pulse展開）一 計算目的：其中 L=常數(shù)，a=常數(shù),（把直導(dǎo)線等分為n求直導(dǎo)線的電容值LL2aPulse展開脈沖函數(shù)將任意形狀的導(dǎo)線L分為N段：。Pulse展開的方程為1、L=3m，a=0.01m，改變分段數(shù)n，從3到50，觀察電容值的變化如圖一。圖一源程序1a=0.01;l=3;o=8.854e-12;i=50;C=ones(1,i-2)syms xfor K=3:1:it=1:1:K;v(t)=(l/K)*t;V=zeros(K,K);for m=1:K for n=1:K if m=n Y=(int(1/sqrt(a2+x2),x,-(l/K)/2,(l/K

2、)/2); V(m,n)=vpa(Y); else V(m,n)=(l/K)/abs(v(m)-v(n); end end U=Vones(K,1); u=U*(l/K); s1=sum(u); c=4*s1*pi*o;endC(1,K-2)=c;endK=3:1:iplot(K,C);hold on text(30,3.17e-011,leftarrowitfontnamearialC=3.17e-011)title(L=3m,a=0.01m,change n)xlabel(n/)ylabel(C/F)grid on2、a=0.01，n=30，改變導(dǎo)線L長度從3m到30m；觀察電容值的變化如

3、圖二。圖二程序2a=0.01;i=30;o=8.854e-12;C=ones(1,i-2);syms xfor l=3:1:30K=30;t=1:1:K;v(t)=(l/K)*t;V=zeros(K,K);for m=1:K for n=1:K if m=n Y=(int(1/sqrt(a2+x2),x,-(l/K)/2,(l/K)/2); V(m,n)=vpa(Y); else V(m,n)=(l/K)/abs(v(m)-v(n); end end U=Vones(K,1); u=U*(l/K); s1=sum(u); c=4*s1*pi*o;endC(1,l-2)=c;endl=3:1:i

4、plot(l,C);hold onstairs(l,C,:r);hold offtext(3,3.17e-011,leftarrowitfontnamearialC=3.17e-011)title(n=3,a=0.01m,change L)xlabel(L/m)ylabel(C/F)grid on3、L=3m，n=30，改變半徑a從0.01m到0.06m圖三程序l=3;i=11;o=8.854e-12;C=ones(1,i);syms x,r=0for a=0.01:0.005:0.06K=30;t=1:1:K;v(t)=(l/K)*t;V=zeros(K,K);for m=1:K for n

5、=1:K if m=n Y=(int(1/sqrt(a2+x2),x,-(l/K)/2,(l/K)/2); V(m,n)=vpa(Y); else V(m,n)=(l/K)/abs(v(m)-v(n); end end U=Vones(K,1); u=U*(l/K); s1=sum(u); c=4*s1*pi*o;endr=r+1;C(1,r)=c;endfor i=1:1:11 p(1,i)=0.01+0.005*(i-1);endplot(p,C);hold onstem(p,C,*r);hold offtitle(n=3,L=3m,change a)xlabel(a/m)ylabel(C/F)grid on4、導(dǎo)線L上電荷分布情況如圖四所示。圖四程序a=0.01;i=30;o=8.854e-12;C=ones(1,i);syms xl=3;K=30;t=1:1:K;v(t)=(l/K)*t;V=zeros(K,K);U=zeros(1.30)for m=1:K for n=1:K if m=n Y=(int(1/sqrt(a2+x2),x,-(l/K)/2,(l/K)/2); V(m,n)=vpa(Y); else V(m,n)=(l/K)/abs(v(m)-v(n); end end U=Vones(K,1)