電磁數(shù)值計算作業(yè)_第1頁
電磁數(shù)值計算作業(yè)_第2頁
電磁數(shù)值計算作業(yè)_第3頁
電磁數(shù)值計算作業(yè)_第4頁
已閱讀5頁，還剩8頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、作業(yè) 15 地面上方的單導體傳輸線：一已知條件和計算目標已知條件：導線的直徑 d=0.2m ，距離地面的高度 H=2m計算目標：用矩量法計算與地面平行的圓柱形長導線的長度電容、電感、特征阻抗隨高度和半徑的變化，并將計算結果與理論結果相比較。三源程序(1)clear;長度電容 C 與高度 H 的關系mu=4*pi*10(-7); ep=(10(-9)/(36*pi); d=0.2;H=0.2:0.1:2; D=2*H;L=;N=40;l=L/N; t=2*pi/N;for i=1:length(H)the/2:t:2*pi;x=cos(theta)*d/2;y=sheta)*d/2+H(i);t

2、heta1=0:t:2*pi-t;x1=cos(theta1)*d/2;y1=sheta1)*d/2+H(i);for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(y(m)+y(n)*yy(n); c1=l2;a2=(x(m)-x(n)2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*

3、xx(n)-2*(y(m)-y(n)*yy(n); c2=l2;s(m,n)=l*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-ln(a 2,b2,c2,-1/2);endend LL=ones(1,N)*l;B0=4*pi*ep*ones(N,1); C(i)=LL*inv(s)*B0;C1(i)=pi*ep/log(D(i)+sqrt(D(i)2-d2)/d);L(i)=ep*mu/C(i);L1=ep*mu/C1(i); L2(i)=L1;Zc(i)=sqrt(ep*mu)/C(i);Zc1=sqrt(L1/C1(i); Zc

4、2(i)=Zc1;endplot(H,C,-*r,LineWidth,1.5); hold on;plot(H,C1,-,LineWidth,1.5); grid onxlabel(導線離地面的距離H(m);ylabel(電容C(F/m);title(長度電容C與高度H的關系);legend(矩量法,理論值)（2）clear;長度電感 L 與高度 H 的關系mu=4*pi*10(-7); ep=(10(-9)/(36*pi); d=0.4;H=0.2:0.1:2; D=2*H;L=;N=40;l=L/N; t=2*pi/N;for i=1:length(H)the/2:t:2*pi;x=cos

5、(theta)*d/2;y=sheta)*d/2+H(i);theta1=0:t:2*pi-t;x1=cos(theta1)*d/2;y1=sheta1)*d/2+H(i);for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(y(m)+y(n)*yy(n); c1=l2;a2=(x(m)-x(n)

6、2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*xx(n)-2*(y(m)-y(n)*yy(n); c2=l2;s(m,n)=l*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-ln(a 2,b2,c2,-1/2);endend LL=ones(1,N)*l;B0=4*pi*ep*ones(N,1); C(i)=LL*inv(s)*B0;C1(i)=pi*ep/log(D(i)+sqrt(D(i)2-d2)/d);L(i)=ep*mu/C(i); L1=ep*mu/C1(i);L2(i)=L1;Zc(i)=sqrt(ep

7、*mu)/C(i);Zc1=sqrt(L1/C1(i); Zc2(i)=Zc1;endplot(H,L,-*r,LineWidth,1.5); hold on;plot(H,L2,-,LineWidth,1.5); grid onxlabel(導體離地面的距離H(m);ylabel(電感L(H/m);title(長度電感L與高度H的關系); legend(矩量法,理論值)（3）直導線特征阻抗 Zc 與高度 H 的關系clear; mu=4*pi*10(-7);ep=(10(-9)/(36*pi); d=0.4;H=0.2:0.1:2; D=2*H;L=;N=40;l=L/N; t=2*pi/N

8、;for i=1:length(H)the/2:t:2*pi;x=cos(theta)*d/2;y=sheta)*d/2+H(i);theta1=0:t:2*pi-t;x1=cos(theta1)*d/2;y1=sheta1)*d/2+H(i);for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(

9、y(m)+y(n)*yy(n); c1=l2;a2=(x(m)-x(n)2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*xx(n)-2*(y(m)-y(n)*yy(n); c2=l2;s(m,n)=l*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-ln(a 2,b2,c2,-1/2);endend LL=ones(1,N)*l;B0=4*pi*ep*ones(N,1); C(i)=LL*inv(s)*B0;C1(i)=pi*ep/log(D(i)+sqrt(D(i)2-d2)/d);L(i)=ep*mu/C(i);L1

10、=ep*mu/C1(i); L2(i)=L1;Zc(i)=sqrt(ep*mu)/C(i); Zc1=sqrt(L1/C1(i);Zc2(i)=Zc1;endplot(H,Zc,-*r,LineWidth,1.5); hold on;plot(H,Zc2,-,LineWidth,1.5);grid onxlabel(導線離地面的距離H(m);ylabel(特征阻抗Zc();title(直導線特征阻抗Zc與高度H的關系); legend(矩量法,理論值)（4）clear;長度電容 C 與導線直徑 d 的關系mu=4*pi*10(-7); ep=(10(-9)/(36*pi); d=0.1:0.1

11、:1;H=2; D=2*H; L=;N=40;l=L/N; t=2*pi/N;for i=1:length(d)the/2:t:2*pi;x=cos(theta)*d(i)/2;y=sheta)*d(i)/2+H;theta1=0:t:2*pi-t;x1=cos(theta1)*d(i)/2;y1=sheta1)*d(i)/2+H;for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-

12、x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(y(m)+y(n)*yy(n); c1=l(i)2;a2=(x(m)-x(n)2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*xx(n)-2*(y(m)-y(n)*yy(n); c2=l(i)2;s(m,n)=l(i)*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-l n(a2,b2,c2,-1/2);endend LL=ones(1,N)*l(i);B0=4*pi*ep*ones(N,1); C(i)=LL*inv(s)*B0;C

13、1(i)=pi*ep/log(D+sqrt(D2-d(i)2)/d(i);L(i)=ep*mu/C(i);L1=ep*mu/C1(i); L2(i)=L1;Zc(i)=sqrt(ep*mu)/C(i); Zc1=sqrt(L1/C1(i);Zc2(i)=Zc1;endplot(d,C,-*r,LineWidth,1.5); hold on;plot(d,C1,-,LineWidth,1.5); grid onxlabel(導線直徑d(m);ylabel(電容C(F/m);title(長度電容C與導線直徑d關系);legend(矩量法,理論值)（5）clear;長度電感 L 與導線直徑 d 的關

14、系mu=4*pi*10(-7); ep=(10(-9)/(36*pi); d=0.1:0.1:1;H=2; D=2*H; L=;N=40;l=L/N; t=2*pi/N;for i=1:length(d)the/2:t:2*pi;x=cos(theta)*d(i)/2;y=sheta)*d(i)/2+H;theta1=0:t:2*pi-t;x1=cos(theta1)*d(i)/2;y1=sheta1)*d(i)/2+H;for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(

15、N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(y(m)+y(n)*yy(n); c1=l(i)2;a2=(x(m)-x(n)2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*xx(n)-2*(y(m)-y(n)*yy(n); c2=l(i)2;s(m,n)=l(i)*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-l n(a2,b2,c2,-1/2);endend LL=ones(1,N)*

16、l(i); B0=4*pi*ep*ones(N,1); C(i)=LL*inv(s)*B0;C1(i)=pi*ep/log(D+sqrt(D2-d(i)2)/d(i);L(i)=ep*mu/C(i);L1=ep*mu/C1(i); L2(i)=L1;Zc(i)=sqrt(ep*mu)/C(i); Zc1=sqrt(L1/C1(i);Zc2(i)=Zc1;endplot(d,L,-*r,LineWidth,1.5); hold on;plot(d,L2,-,LineWidth,1.5); grid onxlabel(導線直徑d(m);ylabel(電感L(H/m);title(長度電感L與導線直

17、徑d的關系); legend(矩量法,理論值)（6）直導線特征阻抗 Zc 與導線直徑 d 的關系clear; mu=4*pi*10(-7);ep=(10(-9)/(36*pi); d=0.1:0.1:1;H=2;D=2*H; L=;N=40;l=L/N; t=2*pi/N;for i=1:length(d)the/2:t:2*pi;x=cos(theta)*d(i)/2;y=sheta)*d(i)/2+H;theta1=0:t:2*pi-t;x1=cos(theta1)*d(i)/2;y1=sheta1)*d(i)/2+H;for j=1:N-1xx(j)=abs(x1(j+1)-x1(j);

18、yy(j)=abs(y1(j+1)-y1(j);endxx(N)=abs(x1(N)-x1(1);yy(N)=abs(y1(N)-y1(1);for m=1:Nfor n=1:Na1=(x(m)-x(n)2+(y(m)+y(n)2;b1=-2*(x(m)-x(n)*xx(n)+2*(y(m)+y(n)*yy(n); c1=l(i)2;a2=(x(m)-x(n)2+(y(m)-y(n)2;b2=-2*(x(m)-x(n)*xx(n)-2*(y(m)-y(n)*yy(n); c2=l(i)2;s(m,n)=l(i)*(ln(a1,b1,c1,1/2)-ln(a1,b1,c1,-1/2)-(ln(a2,b2,c2,1/2)-l n(a2,b2,c2,-1/2);endend LL=ones(1,N)*l(i); B0=4*pi*ep*ones(N,1); C(i)=LL*inv

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論