mom2.m

矩量法计算导体圆环的电容计算matlab程序

%导体圆环的电容计算matlab程序
function  C=NO3
a=0.05; 
b=3;
%N=30;
e=1/(36*pi)*(1e-9);
k=1;
for N=3:1:50
%for b=1:0.1:5
f=2*pi*b/N;  
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)=1/(4*pi*e)*[(1+a^2/(2*f^2))*log((f+sqrt(f^2+a^2))/(-f+sqrt(f^2+a^2)))-sqrt(f^2+a^2)/f];
    else   
S(j,i)=f/(3*pi*e*(2*abs(sin((j-i)*pi/N))));
    end
   end
end
A=double(inv(S)*D);
Q=sum(A);
C(k)=double(2*pi*b*Q/N)
k=k+1;
end
N=3:1:50
plot(N,C(),'r')
%b=1:0.1:5
%plot(b,C(),'r')
grid on
xlabel('分段的段数N’);
ylabel('电容值C(F)')
legend('电容随分段段数N的变化曲线 ')
title('导体圆环的电容分布图
%xlabel('导体圆环的半径b(m)’);
%ylabel('电容值C(F)')
%legend('电容随圆环半径的变化曲线 ')
%title('导体圆环的电容随圆环半径的变化的曲线图')