111.txt

用矩量法计算单条带传输线的电容

interg.m
function  y=interg(a,b,c,x)
z=b^2-4*a*c;
if z>=0
y=(x+b/(2*c))*log(a+b*x+c*x^2)-2*x+sqrt(z+eps)/c*atanh((2*c*x+b)/sqrt(z+eps));
else
y=(x+b/(2*c))*log(a+b*x+c*x^2)-2*x+sqrt(-z)/c*atan((2*c*x+b)/sqrt(-z));
end
程序一：
clear;
clf;
format long;
e0=8.854187818e-12;
W=1;
N=10;
deltl=W/N;
c1=deltl^2;
c2=deltl^2;
for tt=1:20
    h=0.01*tt;
for m=1:N
    for n=1:N
        xm=(m-1)*deltl+deltl/2;  %每段中点的坐标
        ym=h;
        xn=(n-1)*deltl+deltl/2;
        yn=h;
        dxn=deltl;
        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;
s(m,n)=deltl*(interg(a1,b1,c1,0.5)-interg(a1,b1,c1,-0.5)-interg(a2,b2,c2,0.5)+interg(a2,b2,c2,-0.5))
    end
end
A=s\ones(N,1);
C=deltl*sum(A)*(4*pi*e0);
c(tt)=C;
cc1(tt)=2*e0*W/(2*h);
end
tt=(1:20)';
plot(tt,c,'r',tt,cc1)
title('单条带传输线的电容随其与地面的距离的变化曲线')
xlabel('distance(cm)')
ylabel('Capacity(F)')
grid on
legend('数值结果','解析结果')