单条带传输线的矩量法分析.docVIP

作业八 单条带传输线的矩量法分析一、条件和计算目标已知：条带的宽度为W，宽面与地面平行，且与地面距离为h，周围空间介质参数为和，则该导线与它的镜像构成了一对间距为的平行板传输线。求：（1）矩量法求单条带传输线的单位长度电容C、单位长度电感L、特性阻抗 （2）解析法求平行板传输线单位长度电容C、单位长度电感L、特性阻抗（3）比较计算结果与理论公式计算结果二、示意图图1 单导体传输线三、矩量法计算和理论公式(a)矩量法解此问题的关键是根据矩阵方程求出矩阵元素，然后利用公式，求出电容值，然后再求出电感和特性阻抗值。(b)理论公式(F/M)(H/m)四、结果与分析图2 矩量法与解析法单位长度电容结果比较图3 矩量法与解析法单位长度电感结果比较图4 矩量法与解析法单位长度电容结果比较附程序:clear;clf;format long;e0=8.854187818e-12;u0=4*pi*10^(-7);w=0.1;%h=0.1;%d=2*h;N=10;len=w/N;c1=len^2;c2=len^2;syms t;for kk=1:20 h=0.01*kk;for m=1:N for n=1:N xm=(m-1)*len+len/2; %每段中点的坐标 ym=h; xn=(n-1)*len+len/2; yn=h; dxn=len; dyn=0; a1=(xm-xn)^2+(ym+yn)^2; b1=-2*(xm-xn)*dxn+2*(ym+yn)*dyn; a2=(xm-xn)^2+(ym-yn)^2; b2=-2*(xm-xn)*dxn-2*(ym-yn)*dyn; smn(m,n)=int(log((a1+b1*t+c1*t^2)/(a2+b2*t+c2*t^2)),-1/2,1/2); smn(m,n)=len*double(smn(m,n)); endendA=inv(smn)*4*pi*e0*ones(N,1);c(kk)=len*sum(A);l(kk)=u0*e0/c(kk);z(kk)=sqrt(l(kk)/c(kk));CC(kk)=e0*w/(2*h);LL(kk)=u0*e0/CC(kk);ZZ(kk)=sqrt(u0/e0)*2*h/w;endkk=1:20;figure(1)plot(kk,double(c),r-,linewidth,1.5);hold ontitle(单条带传输线的单位长度电容C)xlabel(高度h(m)),ylabel(电容(F))plot(kk,3*CC,b--,linewidth,1.5);legend(矩量法计算结果,理论公式计算结果);grid onfigure(2)plot(kk,double(l),r--,linewidth,1.5);hold onplot(kk,1/3*LL,b--,linewidth,1.5);title(单条带传输线的单位长度电感L)xlabel(高度h(m)),ylabel(电感(H))legend(矩量法计算结果,理论公式计算结果);figure(3)%plot(kk,eval(vpa(z)),k--,linewidth,2);plot(kk,double(z),r-,linewidth,1.5);hold onplot(kk,1/3*ZZ,b--,linewidth,1.5);title(单条带传输线的特性阻抗Zc(\Omega))xlabel(高度h(m)),ylabel(特性阻抗(\Omega))legend(矩量法计算结果,理论公式计算结果);分析:由上述几个图可知，当取时，矩阵法计算得出的单位长度电容等于理论公式电容值的两倍，而计算得出的单位长度电感及特性阻抗等于理论公式电感和特性阻抗计算值的一半。附程序:%--------------单导体传输线----------------clear all%clfe0=8.854*10^(-12);u0=4*pi*10^(-7);%h=0.2;%离地面高度d=0.01;%直径kk=1;N=10;for h=1.8:0.1:3circle=d*pi;%周长len=circle/N;%每段的长度syms t;D=2*h;for m=1:N for n=1:N x1=d/2*cos(pi/N+(m-1)*2*pi/N); y1=d/2*sin(pi/N+(m-1)*2*pi/N)+h; x2=d/2*cos(pi/N+(n-1)*2*pi/N);