利用matlab实现光波（电磁波）的仿真

根据电磁波余弦曲线公式

用matlab进行三维做图，做出电磁波仿真。

%Electromagnetic wave simulation
%假定电磁波沿x轴方向传播，电磁波在y轴方向和z轴方向有分量
Ey_m=6;  %电磁波y分量幅值
Ez_m=8;  %电磁波z分量幅值
E=(Ey_m^2+Ez_m^2)^0.5;   %求出电磁波的总幅值
w=10;  %电磁波的角速度为10
a_1=pi/4;  %y分量电磁波的初相位
a_2=pi/2;  %z分量电磁波的初相位
%a_1和a_2的大小差值，决定了相位差
%位相差与分量的幅值进一步决定了电磁波为何种程度的极化电磁波（线极化，圆极化，椭圆极化）
x=0:0.01:3;  %设置沿x轴方向传播的点位
m_0=zeros(size(x));  %m_0为与x矢量同尺寸的0坐标参考点位

figure
for t=0:200
    Ey=Ey_m*cos(w*x+w*t*1e-2+a_1)  %沿x轴方向传播的波，在y-x平面上的波函数
    Ez=Ez_m*cos(w*x+w*t*1e-2+a_2)  %沿x轴方向传播的波，在z-x平面上的波函数
    
    plot3(x,m_0,m_0,'r','linewidth',2.5)  %绘制参考轴线
    hold on
    plot3(x,Ey,m_0,'b','linewidth',1.0)  %绘制y-x平面的波形图
    hold on
    plot3(x,m_0,Ez,'k','linewidth',1.0)  %绘制z-x平面的波形图
    hold on
    plot3(x,Ey,Ez,'g','linewidth',0.8)  %绘制出空间三维的波形图
    hold off
    
    title('光波（电磁波）空间传播图')
    xlabel('光波传播方向x')
    ylabel('光波电矢量分量Ey')
    zlabel('光波磁矢量分量Bz')
    
    drawnow
end