根据电磁波余弦曲线公式
用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