calc3.m

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

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