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

以下是一个基于有限差分法（FDTD）的动态Matlab代码，用于模拟一维行波电磁波的传播和反射。其中，假设空间中只存在电场，而磁场可以通过安培定理计算得出。

% 模拟参数
c = 3e8;         % 光速
f = 1e9;         % 频率
lambda = c/f;    % 波长
dx = lambda/10;  % 空间步长
dt = dx/c/2;     % 时间步长
t_max = 2/f;     % 模拟总时间
x_max = 10*lambda;% 模拟空间范围

% 网格初始化
x = (-x_max:dx:x_max);   % 空间网格
t = (0:dt:t_max);        % 时间网格
E = zeros(size(x));      % 电场

% 初始条件
E(x <= 0) = 1;           % 左侧为驻波源

% 模拟循环
for n = 1:length(t)
    % 计算电场在时刻n+1的值
    E(2:end-1) = E(2:end-1) + (c*dt/dx)^2*(E(3:end) - 2*E(2:end-1) + E(1:end-2));
    
    % 添加边界条件：电场为0
    E(1) = 0;
    E(end) = 0;
    
    % 绘制电场随时间和空间的变化
    plot(x, E);
    axis([-x_max x_max -2 2]);
    xlabel('位置 (m)');
    ylabel('电场强度');
    title(sprintf('时刻 %.2f s', t(n)));
    drawnow;
end

在这个代码中，我们使用了一个二阶有限差分法来计算电场在下一个时间步长的值，并添加了边界条件来模拟电磁波在两端的反射。我们还用`plot`函数在每个时间步长绘制了电场在空间中的分布图，以便观察电磁波的传播和反射。