4_8.rar_matlab 电动力学_场分布_电力线_电动力学_电动力学_matlab

thi=[1:0.1:180]*pi/180; dthi=thi(2)-thi(1); Le=length(thi); N=10; kd=2*pi*9/40; beta=11/20*pi; F=zeros(1,Le); U0=0; Nor=1/N*sin(pi/2)/sin(pi/(2*N)); for i=1:Le %numerator Nu=sin(N/2*(kd*cos(thi(i))+beta)); %demomenator De=sin(1/2*(kd*cos(thi(i))+beta)); F(i)=(1/N)*Nu/De/Nor; U0=U0+0.5*F(i)^2*sin(thi(i))*dthi; end 1/U0