xizhidaoxiandianrong.m

用matlab编写的关于天线问题的程序

clear all;close all;

e=8.854e-12;

L=3;

a=0.05;

N=30;

for mm=1:N;
    
    nn=mm;
    
    delte=L/nn;
    
    O=ones(nn,1);
    
    for ii=1:nn;
        
        S(ii,ii)=log(delte/2+sqrt(a^2+(delte/2)^2))-log(-delte/2+sqrt(a^2+(-delte/2)^2));
        
    end
    
    for ii=1:nn;
        
        for jj=1:nn;
            
            if ii~=jj;
                
                S(ii,jj)=1/(abs(ii-jj));
                
            end
            
        end
        
    end
    
    P=inv(S)*O;
    
    A=sum(P,1);
    
    C(mm,1)=4*pi*e*A*delte;
    
end

C_temp=zeros(N+1,1);

for nn=1:N;
    
    C_temp(nn+1,1)=C(nn,1);
    
end

mm=1:N;

plot(mm,C,'b:');

xlabel('分段段数') ;

ylabel('电容值C') ;