%導體圓環的電容計算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('導體圓環的電容隨圓環半徑的變化的曲線圖')