仿真平面电磁波的传播、反射及透射代码

下面是一个简单的MATLAB代码，用于仿真平面电磁波在介质分界面上的传播、反射和透射。

% 定义常数
c = 3e8;         % 光速
mu0 = pi*4e-7;   % 真空磁导率
eps0 = 8.85e-12; % 真空介电常数

% 定义介质的参数
n1 = 1;          % 第一种介质的折射率
n2 = 1.5;        % 第二种介质的折射率

% 定义平面波参数
f = 1e9;         % 频率
lambda = c/f;    % 波长
k = 2*pi/lambda; % 波数
omega = 2*pi*f;  % 角频率

% 定义模拟区域
xMin = -0.05;    % 区域最小值
xMax = 0.05;     % 区域最大值
yMin = -0.05;    
yMax = 0.05;
z = 0;

% 定义网格
dx = lambda/10;  % 网格间距
dy = lambda/10;
Nx = round((xMax-xMin)/dx);
Ny = round((yMax-yMin)/dy);

% 创建网格
[x, y] = meshgrid(xMin:dx:xMax, yMin:dy:yMax);

% 计算电场和磁场
E1 = exp(1i*(k*n1*z - omega*t));
E2 = exp(1i*(k*n2*z - omega*t));
H1 = E1/(n1*mu0*c);
H2 = E2/(n2*mu0*c);

% 计算反射和透射系数
r = (n1-n2)/(n1+n2);
t = 2*n1/(n1+n2);

% 计算反射和透射电场和磁场
Er = r*E1;
Et = t*E1;
Hr = Er/(n1*mu0*c);
Ht = Et/(n2*mu0*c);

% 计算总电场和磁场
E = E1 + Er;
H = H1 + Hr;

% 绘制电场图像
figure;
quiver(x, y, real(E), imag(E));
title('Electric Field');

% 绘制磁场图像
figure;
quiver(x, y, real(H), imag(H));
title('Magnetic Field');

上述代码中，我们首先定义了常数和介质的参数，然后定义了平面波的参数和模拟区域。接着，我们创建了网格，并计算了电场和磁场。然后，我们计算了反射和透射系数，并计算了反射和透射电场和磁场。最后，我们计算了总电场和磁场，并绘制了电场和磁场的图像。

需要注意的是，这只是一个简单的示例代码，实际的仿真过程可能需要更复杂的计算和处理。同时，仿真结果也需要经过验证才能得出可靠的结论。