3. 光电路仿真引擎与数值方法

光电路级仿真的核心任务，是在给定器件模型和激励条件下，求解光信号在由波导、耦合器、谐振腔、调制器、探测器等构成的网络中的传输与响应。与电子电路仿真类似，光电路仿真也需要一套“引擎”来组织模型、建立方程并求解；但光信号具有复数幅度与相位、宽频带、多模、双向传输等特性，因此其数值方法比纯电阻网络复杂得多。理解频域、时域和频时混合三类引擎的适用边界，以及S参数、传输矩阵、耦合模理论三种模型语言如何嵌入引擎，是正确选择仿真策略的前提。

一、三类仿真引擎的数学本质

频域仿真把光场写成复振幅随频率变化的函数，对线性时不变光网络直接求解稳态响应。其基本形式是：给定输入频谱 E_in(ω)，通过系统传递函数 H(ω) 得到输出 E_out(ω)=H(ω)E_in(ω)。频域引擎的优点是精度高、速度快，尤其适合分析滤波器、微环谐振器、波分复用器等窄带或周期响应器件；缺点是难以直接处理非线性、时变调制和宽带瞬态。

时域仿真把光场写成随时间演化的复包络 E(t)，通过求解微分方程或差分方程推进。对于包含载流子色散、热光调谐、饱和吸收等动态效应的器件，时域引擎是自然选择。例如，硅基微环调制器的载流子寿命约纳秒量级，必须用时域方法才能捕捉调制边沿和瞬态啁啾。时域仿真的代价是步长受限于最小时间常数，长脉冲或高Q腔需要大量迭代。

频时混合仿真则把网络拆成“线性频域部分”和“非线性/时变时域部分”，在两者之间做域变换。典型做法是：无源波导、耦合器、延迟线用频域S参数描述，有源调制器或探测器用时域状态方程描述，通过傅里叶变换和逆变换在接口处交换复包络。混合引擎兼顾了线性部分的效率和有源部分的动态精度，是硅光链路级验证中最常用的折中方案。

二、S参数：多端口网络的“黑盒”描述

S参数是频域引擎最自然的模型语言。对于一个 N 端口光网络，定义入射波列向量 a 和出射波列向量 b，则 b=Sa。在硅光中，常见端口包括输入/输出波导、上/下路微环端口、调制器射频端口（若做光电协同）。S参数矩阵的每个元素 S_ij(ω) 表示从端口 j 入射到端口 i 出射的复传输系数，包含幅度和相位。

实际实现时，S参数通常由全波仿真（FDTD、FDFD、EME）或测试提取，然后以Touchstone或CSV格式导入电路级仿真器。仿真器在频域对每个频点做矩阵乘法，再通过逆傅里叶变换得到时域波形。需要注意的是，S参数必须满足无源性和因果性；若提取频带过窄或端口去嵌入不干净，会在时域产生非物理的前沿响应。一个实用的检查方法是：对 S₁₁(ω) 做逆变换，观察其冲激响应是否在 t<0 时为零。

三、传输矩阵：级联结构的“乘法链”

传输矩阵（TMM）把一段波导或一个器件的输入/输出场用2×2矩阵联系起来。对于布拉格光栅、多层抗反射膜、级联微环等周期或分层结构，TMM的优势在于级联只需矩阵相乘，且能直接给出反射谱和透射谱。以均匀布拉格光栅为例，其传输矩阵可写为：

M = [ cosh(γL) - j(Δβ/γ)sinh(γL)    -j(κ/γ)sinh(γL) ]
    [ j(κ/γ)sinh(γL)                 cosh(γL) + j(Δβ/γ)sinh(γL) ]

其中 κ 为耦合系数，Δβ 为相位失配，γ²=κ²−Δβ²。在电路级仿真中，TMM通常用于频域引擎，因为矩阵元素是频率的函数。若要在时域使用，需要先拟合为有理函数或状态空间模型，再转为时域卷积或微分方程。

四、耦合模理论：谐振与动态的“状态方程”

耦合模理论（CMT）把谐振腔或耦合区用少数几个模式幅度描述，得到一阶微分方程组。以单微环为例，忽略反向耦合时：

da/dt = (jω₀ - γ_total)a + jμ·s_in
s_out = s_in + jμ·a

其中 a 为腔模复幅度，ω₀ 为谐振频率，γ_total=γ_i+γ_e 为总衰减率，μ 为输入耦合系数。CMT的突出优点是物理图像清晰、参数少、可直接嵌入时域引擎，并且能自然扩展到载流子色散、热光调谐等动态效应。例如，在微环调制器中，把 γ_total 写成载流子浓度 N 的函数，即可模拟调制响应。CMT的局限是假设模式正交、弱耦合、慢变包络，对强耦合或宽带场景需要修正。

五、如何选择：一个可操作的决策流程

面对具体仿真任务，可按以下顺序判断：

1. 看信号带宽与器件Q值。若信号带宽远小于器件自由光谱范围，且器件为线性无源，优先频域+S参数或TMM。
2. 看是否存在时变或非线性。若包含调制器、探测器、热光开关，必须引入时域或混合引擎，有源部分用CMT或速率方程。
3. 看是否需要链路级指标。眼图、误码率、星座图需要时域波形，建议混合引擎：无源链路用S参数频域加速，有源器件用CMT时域推进。
4. 看模型来源。若已有全波提取的宽带S参数，优先频域；若只有CMT参数，优先时域或混合。
5. 做收敛性检查。频域检查频带边缘是否衰减足够；时域检查步长是否小于最小时间常数；混合检查FFT窗口是否覆盖脉冲响应长度。

总结来说，频域、时域、频时混合三类引擎并非互斥，而是对应不同的物理假设和计算代价。S参数、传输矩阵、耦合模理论则是三种模型语言，分别擅长多端口黑盒、级联周期结构和谐振动态。实际仿真中，先明确信号特征和器件物理，再选择引擎与模型语言的组合，并通过收敛性检查验证结果，才能既保证精度又控制计算量。