### 如何用Python实现一个简单的水平层状模型 在处理水平层状模型时，通常会涉及光在不同介质中的传播问题。以下是一个基于传输矩阵法（TMM）的简单水平层状模型示例代码，用于计算光通过多层介质时的反射率和透射率[^4]。 ```python import numpy as np # 定义传输矩阵法函数 def transfer_matrix_method(n_list, d_list, wavelength, theta_i): """ 使用传输矩阵法计算多层水平结构的反射率和透射率。 参数: n_list: 各层折射率列表 (包括入射介质和透射介质) d_list: 各层厚度列表 (单位为纳米)，首尾为无穷大 wavelength: 光波长 (单位为纳米) theta_i: 入射角 (单位为弧度) 返回: r: 反射率 t: 透射率 """ # 初始化变量 k_list = [2 * np.pi * n / wavelength for n in n_list] cos_theta_list = [np.sqrt(1 - (n_list[0] / n) ** 2 * np.sin(theta_i) ** 2) for n in n_list] # 构建传输矩阵 M = np.eye(2, dtype=complex) for i in range(1, len(n_list) - 1): delta = k_list[i] * d_list[i] * cos_theta_list[i] phi = np.array([[np.exp(-1j * delta), 0], [0, np.exp(1j * delta)]]) T_layer = 1 / (2 * cos_theta_list[i]) * np.array([ [cos_theta_list[i-1] + cos_theta_list[i], cos_theta_list[i-1] - cos_theta_list[i]], [cos_theta_list[i-1] - cos_theta_list[i], cos_theta_list[i-1] + cos_theta_list[i]]]) M = np.dot(np.dot(T_layer, phi), M) # 计算反射率和透射率 r = M[1, 0] / M[0, 0] t = (2 * cos_theta_list[0]) / M[0, 0] R = abs(r) ** 2 T = abs(t) ** 2 * (cos_theta_list[-1] / cos_theta_list[0]) return R, T # 示例参数 n_list = [1.0, 1.5, 2.5, 1.0] # 折射率列表 (空气、第一层、第二层、空气) d_list = [np.inf, 100, 200, np.inf] # 厚度列表 (单位为纳米) wavelength = 500 # 波长 (单位为纳米) theta_i = 0 # 入射角 (单位为弧度) # 调用函数 R, T = transfer_matrix_method(n_list, d_list, wavelength, theta_i) print(f"反射率: {R}, 透射率: {T}") ``` 上述代码使用传输矩阵法（TMM）来模拟光通过多层水平结构的情况。用户可以调整 `n_list` 和 `d_list` 来定义不同的层状结构，并通过改变波长或入射角进一步分析其光学特性。 ---