mom1.m

矩量法计算长直线导体电容的计算matlab程序

%长直线导体电容的计算
function  C=NO1                
L=3;  
% N=30;                   
a=0.01;                         
e=1/(36*pi)*(1e-9);          
k=1;
%for N=3:1:50
for L=1:10
    m=N;
    n=N;
    S=rand(m,n);                      
    A=rand(m,1);            
    D=ones(m,1); 
syms l                      
for j=1:m
   for i=1:n
    if i==j
S(i,i)=double(int(1/sqrt(a^2+l^2),-L/(2*N),L/(2*N)));
    else   
S(j,i)=1/abs(j-i);
    end
   end
end
A=double(4*pi*e*inv(S)*D);    
Q=sum(A);
C(k)=double(Q*(L/N))
k=k+1;
end
%N=3:1:50;
L=1:10
%plot(N,C(),'r')
grid on
%xlabel('x轴（分段的段数N）’);
%ylabel('y轴（电容值）')
%legend('电容随分段段数N的变化曲线 ')
%title('直线导体的电容分布图')

plot(N,C(),'r')
grid on
xlabel('x轴（直线导体的长度L(m)）’);
ylabel('y轴（电容值）')
legend('电容随直线导体的长度L的变化曲线 ')
title('直线导体的电容随直线导体的长度L的变化曲线图')