calc.m

来自「计算9格均匀电流面元阵方向图的matlab程序」· M 代码 · 共 11 行

[theta,phi]=meshgrid(0:0.02*pi:2*pi,0:pi/12:1.1*pi);
the=pi/2-theta;
for kk=1:100
    aaa=-0.8+kk/500;
F=@(the,phi)cell3(the,phi,aaa).^2.*sin(the);
A=dblquad(F,0,pi,0,2*pi);
B=max(max(cell3(the,phi,aaa)));
q(kk)=4*pi/A*B*B
log10(q)*10;
end