6_12_directivity.m

来自「层结构的电力线分布图」· M 代码 · 共 20 行

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