写出行驻波电磁波的动态matlab代码

以下是一份简单的行波电磁波的动态MATLAB代码：

%% 行波电磁波的动态MATLAB代码
clear all;
close all;
clc;

%% 参数设置
c = 3e8; % 光速
f = 2e9; % 频率
lambda = c/f; % 波长
k = 2*pi/lambda; % 波数
w = 2*pi*f; % 圆频率
Vp = 1; % 正向波振幅
Vm = 1; % 反向波振幅
d = lambda/2; % 距离
t = linspace(0, 4/f, 400); % 时间

%% 计算电场强度和磁场强度
E = Vp*cos(w*t-k*d) + Vm*cos(w*t+k*d); % 电场强度
H = (Vp/Z0)*cos(w*t-k*d) - (Vm/Z0)*cos(w*t+k*d); % 磁场强度

%% 绘图
figure(1);
subplot(2, 1, 1);
plot(t, E, 'r', 'linewidth', 2);
grid on;
xlabel('时间 (s)', 'fontsize', 14);
ylabel('电场强度 (V/m)', 'fontsize', 14);
title('行波电磁波的动态演示', 'fontsize', 16);
subplot(2, 1, 2);
plot(t, H, 'b', 'linewidth', 2);
grid on;
xlabel('时间 (s)', 'fontsize', 14);
ylabel('磁场强度 (A/m)', 'fontsize', 14);

这段代码可绘制出时间轴上电场强度和磁场强度的动态变化，从而形象地展示出行波电磁波的传播特性。