couple_line.m

用数值算法

function y1=couple_line( )
clear;
clf;
format long
e0=8.854187818e-12;
u0=pi*(4e-7);
Wt=1;
W=0.1;
%S=0.1;
for S=0.01:0.01:0.5
b=1;
N1=10;
N2=10;
Nt=100;
deltl=W/N1;%%%deltl=0.01%%%
c1=deltl^2;
c2=c1;
V11=ones(N1,1);
V10=zeros(N2,1);
V1=[V11;V10];
V2=[zeros(N1,1);ones(N2,1)];
for m=1:N1
    for n=1:N1
        x(m)=-(W+S/2)+(m-1)*deltl+deltl/2;%条带1每段中点的坐标
        y(m)=b/2;
        x(n)=-(W+S/2)+(n-1)*deltl+deltl/2;
        y(n)=b/2;
        dxn=deltl;
        dyn=0;
    end
end
for m=N1+1:N1+N2
    for n=N1+1:N1+N2
        x(m)=S/2+(m-1)*deltl+deltl/2;%条带2每段中点的坐标
        y(m)=b/2;
        x(n)=S/2+(n-1)*deltl+deltl/2;
        y(n)=b/2;
        dxn=deltl;
        dyn=0;
    end
end
for m=N1+N2+1:N1+N2+Nt
    for n=N1+N2+1:N1+N2+Nt
        x(m)=-Wt/2+(m-1)*deltl+deltl/2;
        y(m)=b;
        x(n)=-Wt/2+(n-1)*deltl+deltl/2;
        y(n)=b;
        dxn=deltl;
        dyn=0;
    end
end
for m=1:N1+N2
    for n=1:N1+N2
        a1=(x(m)-x(n))^2+(y(m)+y(n))^2;
        b1=-2*(x(m)-x(n))*dxn+2*(y(m)+y(n))*dyn;
        a2=(x(m)-x(n))^2+(y(m)-y(n))^2;
        b2=-2*(x(m)-x(n))*dxn-2*(y(m)-y(n))*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
A1=s\V1;
A2=s\V2;
%c11=(4*pi*e0)*sum(A1(1:N1))*deltl;
%c12=(4*pi*e0)*abs(sum(A2(1:N1)))*deltl;
for i=1:5
    if(i/100==S)            c11(i)=(4*pi*e0)*sum(A1(1:N1))*deltl;%%%%%%%%%%%%%c11=c22
    %c22=(4*pi*e0)*sum(A2(N1+1:N1+N2))*deltl;
    c12(i)=(4*pi*e0)*abs(sum(A2(1:N1)))*deltl;
    %c21=(4*pi*e0)*abs(sum(A1(N1+1:N1+N2)))*deltl;
    end
end
end
plot(c11,'r','LineWidth',2)
hold on
plot(c12,'b','LineWidth',2)
grid on
title('耦合带状线')
xlabel('distance(cm)')
ylabel('Capacity(F)')
legend('自电容值','互电容值')