给出绘制垂直极化波由空气(z<0)斜入射理想导体(z>0)平面(z=0，媒质分界面位于xoy面)的电场矢量图和磁场矢量图的Matlab代码

以下是绘制垂直极化波由空气(z<0)斜入射理想导体(z>0)平面(z=0，媒质分界面位于xoy面)的电场矢量图和磁场矢量图的Matlab代码：

% 定义入射波的频率、波长、振幅、极化方向和入射角度
f = 1e9; % 1GHz
lambda = 3e8 / f;
E0 = 1;
theta_i = 45 * pi / 180;

% 定义导体的位置
z0 = 0;

% 定义网格范围和步长
x = -lambda:lambda/50:lambda;
z = 0:lambda/50:lambda;
[X,Z] = meshgrid(x,z);

% 定义介电常数
epsilon1 = 8.854e-12; % 真空的介电常数
epsilon2 = Inf; % 无限大的介电常数，即理想导体

% 计算电场和磁场
E = zeros(size