模型1.m

来自「时变电磁场的数值分析与仿真」· M 代码 · 共 20 行

clear;clc;
e0=1;
A=0:10;
Z=-50:50;
[a,z]=meshgrid(A,Z);
for d=1:32
    if z>0
        E=d.*(1-z./(z.^2+a.^2).^0.5)/(2*e0);
    else E=-d.*(1+z./(z.^2+a.^2).^0.5)/(2*e0);
    end
    surf(a,z,E);
    axis([0 10 -50 50 -30 30]);
    view([230 40]);
    xlabel('a'), 
    ylabel('z'),
    zlabel('E(z)'),
    m(d)=getframe;
end
    movie(m,2);