用matlab写出完整代码实现上述功能

在MATLAB中实现电磁波的反射和折射，需要用到一些专门的工具箱，如`pdepe`用于求解偏微分方程，`wavefun`用于创建和操作波函数。由于这里涉及到的内容较多，我会提供一个简单的示例，实际应用可能需要更复杂的设置。以下是一个基于波动方程的一维简化版本：

% 定义介质参数
medium1 = [1, 0]; % 第一种介质，假设为理想导体
medium2 = [1, 0.001]; % 第二种介质，例如空气

% 设定边界条件
boundaries = [-L, L]; % 分界面位置
boundaryconds = @(x,D) [D(1); 0]; % 边界上E和H的边界条件（无反射）

% 创建PDE对象并设置方程
pde = pdepe('laplace', 'u', boundaryconds, [], medium1, boundaries);

% 定义初始条件和边界
ic = u; % 初始条件，可以根据实际情况设定
bc = {[-L, 0], [L, 0]}; % 左右边界条件

% 求解波函数
[t, x, u] = pdepe(pde, ic, bc);

% 计算反射和折射
reflected_phase = -atan(1/medium2(2)); % 反射相位
incident_phase = atan(1/medium1(2)); % 入射相位
reflected_field = exp(-1i * (t + incident_phase)) * u(L);
transmitted_field = exp(-1i * (t + reflected_phase)) * u(-L);

% 输出或可视化结果
disp(['Reflection coefficient: ', num2str(abs(reflected_field)^2)]);
disp(['Transmission coefficient: ', num2str(abs(transmitted_field)^2)]);

% 简单的波动模式可视化，这取决于实际的波函数形状
plot(x, real(u));
xlabel('Distance');
ylabel('Electric Field');