jfk4.asv

来自「电磁辐射的数值分析与仿真（基于MATLAB语言）。包括两个点电荷的电场」· ASV 代码 · 共 14 行

clear;clc;
I=10;l=1;omega=10*pi;epsilon=0.7;lambda=10;z=100;
theta=0:pi/10:pi*2;r=100:50:700;
c2=I*l/(2*pi*lambda);
[Theta,R]=meshgrid(theta,r);
Hphi=c2.*sin(Theta)./(R.^2*2);
Er=c2.*cos(Theta)./(R.^3*omega*epsilon);
Etheta=c2.*z.*sin(Theta)./(R.^3*2*omega*);
subplot(1,3,1),surf(Theta,R,Hphi);axis([0 6 0 700 -5e-3 5e-3])
xlabel('\theta'),ylabel('r'),zlabel('H_{\phi}'),title('远区场H_{\phi}')
subplot(1,3,2),surf(Theta,R,Er);axis([0 6 0 700 -6e-5 6e-5])
xlabel('\theta'),ylabel('r'),zlabel('E_r'),title('远区场E_r')
subplot(1,3,3),surf(Theta,R,Etheta);axis([0 6 0 700 -0.5 0.5])
xlabel('\theta'),ylabel('r'),zlabel('E_{\theta}'),title('远区场E_{\theta}')