dantiaodai.m

导体的电感 电容 特性阻抗曲线 的ＭＡＴＬＡＢ模拟

clear all;close all;
format long;
e0=8.854187818e-12;
u=pi*4e-7;
W=1;
N=30;
deltl=W/N;
c1=deltl^2;
c2=deltl^2;
for tt=1:25
    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)=e0*W/h;
    l(tt)=e0*u/c(tt);
    l1(tt)=u*h/W;
    z(tt)=sqrt(l(tt)/c(tt));
    z1(tt)=sqrt(u/e0)*(h/W);
end
tt=(1:25)';
plot(tt,c,'r',tt,cc1,'bo')
xlabel('距离h(m)');
ylabel('电容 C(F/m)');
legend('矩量法','解析方法');
title('电容随距离h变化曲线图');
grid on;
figure;
hold on;
tt=(1:25)';
plot(tt,l,'r',tt,l1,'bo')
xlabel('距离h(m)');
ylabel('电感L(H/m)');
legend('矩量法','解析方法');
title('电感随距离h变化曲线图')
grid on;
figure;
hold on;
tt=(1:25)';
plot(tt,z,'r',tt,z1,'bo')
xlabel('距离h(m)');
ylabel('特性阻抗Zc(Ω)');
legend('矩量法','解析方法');
title('特性阻抗随距离h变化曲线图')
grid on;