给出绘制xoz面合成电场幅度图和合成磁场矢量图的matlab代码

以下是绘制xoz面合成电场幅度图和合成磁场矢量图的Matlab代码：

% 定义空气和导体的坐标范围
xmin = -1; xmax = 1;
zmin = -1; zmax = 1;

% 定义空气和导体的网格数量
Nx = 101;
Nz = 101;

% 定义空气和导体的网格步长
dx = (xmax - xmin) / (Nx - 1);
dz = (zmax - zmin) / (Nz - 1);

% 定义空气和导体的网格点坐标
x_air = linspace(xmin, xmax, Nx);
z_air = linspace(zmin, zmax, Nz);

% 定义入射波的波长和角频率
lambda = 3e8 / f;
omega = 2 * pi * f;

% 计算入射波的波矢和相位差
kz = 2 * pi / lambda * sin(theta_i);
delta = kz * z_air;

% 计算反射波的相位差和振幅
delta_r = -delta;
E_r = -E0;

% 计算空气中的合成电场和磁场
E = E0 * cos(theta_i) + E_r * cos(theta_r) * exp(1i * 2 * delta_r);
H = E / sqrt(epsilon0);

% 计算空气中的电场和磁场分量
Ex = E * sin(theta_i);
Ez = E * cos(theta_i);
Hx = H * sin(theta_i);
Hz = H * cos(theta_i);

% 绘制xoz面合成电场幅度图
figure;
imagesc(x_air, z_air, abs(E));
set(gca, 'YDir', 'normal');
colorbar;
xlabel('x (m)');
ylabel('z (m)');
title('合成电场幅度图');

% 绘制xoz面合成磁场矢量图
figure;
quiver(x_air, z_air, Hx, Hz);
xlim([xmin, xmax]);
ylim([zmin, zmax]);
xlabel('x (m)');
ylabel('z (m)');
title('合成磁场矢量图');

需要注意的是，由于是垂直极化波，电场和磁场的方向都垂直于入射面。因此，在绘制磁场矢量图时，只需要考虑磁场在 xoz 平面上的分量，即 $H_x$ 和 $H_z$。