综述：无线电路仿真中多时间偏微分-代数方程的高级数值方法

《Axioms》：Advanced Numerical Methods for Multitime Partial Differential–Algebraic Equations in Wireless Circuit Simulation

　　

现代无线通信电路的仿真仍然具有挑战性，因为存在非线性行为、异构子系统以及广泛分离的时间尺度。本综述对求解射频（RF）和微波系统电路级建模中出现的多时间偏微分-代数方程（MPDAEs）的高级数值方法进行了结构化概述。与之前的综述论文相比，本工作的主要贡献在于根据底层数值策略对文献进行组织，区分了纯时域、混合时频、多维频域和电路块划分方法。综述的方法表明，多时间公式可以在传统仿真技术上实现显著的计算收益，特别是对于多速率和多尺度电路。时域技术在强非线性条件下通常更稳健，而频域和混合方法在波形可以用有限数量的谐波表示时往往更高效。电路块划分通过利用活跃变量和潜在变量进一步提高了效率，但MPDAE方法的计算复杂度随时间尺度数量迅速增加，并且它们对非周期或高度通用的多速率激励的适用性变得更加有限。总体而言，本综述强调了当前基于MPDAE的数值方法的优势和实践局限性，并指出了未来研究的开放挑战。

论文主体部分总结如下：

2. Theoretical Background Material

2.1. Mathematical Model of an Electronic Circuit—General DAE Formulation
非线性微分-代数方程（DAE）系统被广泛用于建模电子电路的动态行为，形式为 q′(x(t)) + f(x(t)) = s(t)，其中 s(t) 为独立激励向量，x(t) 为状态变量（电压、电流）向量。通过改进节点分析（MNA）结合基尔霍夫电流定律（KCL）导出，函数 q 描述动态元件（电容、电感），f 描述无记忆元件（电阻、非线性电流源）。以包含电流源、线性电感、电导、非线性电容和非线性电流源的示例电路说明，节点电压和电感电流构成的 DAE 系统可划分为记忆函数和无记忆函数。

2.2. Initial Value Problems (Transient Simulation)
瞬态仿真求解给定初始条件 x(0) = x₀ 下的初始值问题，常用 Runge–Kutta（RK）方法或线性多步法（LMM）进行时间步积分。SPICE 类工具普遍采用二阶 LMM（如 Gear-2、梯形法）。此技术直接追踪电路时域响应，构成模拟的基本方法。

2.3. Boundary Value Problems (Steady-State Simulation)
稳态仿真关注周期性稳态行为，需满足周期边界条件 x(T) = x(0)，转化为两点边值问题。数值求解主要有时域法和频域法。时域法采用 shooting-Newton 迭代，频域法采用谐波平衡（HB）方法，两者分别适用于不同场景。

2.3.1. Time-Domain Technique
Shooting 法通过状态转移函数 Φ 将边值问题转化为求解 Φ(x(0)) = x(0)，应用牛顿法更新初值。每步迭代需计算灵敏度矩阵，它由各时间步的偏导矩阵链式乘积得到，这些矩阵依赖于具体时间步积分格式（如 Gear-2）。迭代至满足误差容限为止。

2.3.2. Frequency-Domain Technique
HB 法将周期解展开为傅里叶级数，将 DAE 转换为频域代数方程组。通过牛顿–拉夫逊迭代求解，Jacobian 矩阵由 Toeplitz 转换矩阵组成，涉及离散傅里叶变换（DFT）。对于大规模系统，常使用 GMRES 等 Krylov 子空间迭代法求解稀疏线性系统。

3. Problems Modeled by Partial Differential Systems

3.1. Multitime Formulation
针对包络调制信号（如 A(t₁)cos(ωt₂)），引入两个人工时间尺度：慢时间 t₁ 和快时间 t₂，将一维信号映射为二维函数 ?x(t₁, t₂)。DAE 系统转化为多时间偏微分-代数方程（MPDAE）：?q/?t₁ + ?q/?t₂ + f(?x) = ?s(t₁, t₂)。该表示法可更紧凑地描述多尺度信号，减少采样点数，且一维解可通过取对角线 t₁ = t, t₂ = t mod T 恢复。对于 m 个不同变化率的时间尺度，可推广至 m 维 MPDAE。

3.2. Initial and Boundary Conditions for the Multitime System
对于包络调制激励，需要慢时间方向的初始条件和快时间方向的周期边界条件，构成初始-边值问题。对于准周期激励（如多音信号），两个方向均需周期边界条件，形成双周期边值问题。通过求解这些条件，可在矩形域上计算二维解，并通过锯齿路径恢复一维解。

4. Numerical Methods for Solving the Multitime System

4.1. Time-Domain Envelope Following
该方法使用后向微分公式（如 Gear-2）离散慢时间导数，将 MPDAE 转化为一系列关于快时间周期的边界值问题。每个边界值问题通过 shooting 法求解，快时间积分采用 Gear-2 等 LMM。在每步 shooting 迭代中，需通过链式求导计算灵敏度矩阵。该方法适用于强非线性电路。

4.2. Hybrid Time–Frequency Envelope Following
慢时间方向同样用 BDF 离散，但每个快时间周期问题由 HB 法求解。在每个慢时间点，构建 HB 系统并利用牛顿法迭代，Jacobian 矩阵包含随时间变化的 Toeplitz 子块。该方法在波形谐波成分较少时高效，适用于中度非线性场景。

4.3. Multidimensional Frequency-Domain Computation
对于双周期问题，对两个时间方向同时应用 HB，形成二维 HB 系统。状态变量和激励表示为二维傅里叶级数，MPDAE 转化为频域代数方程组。Jacobian 矩阵由 Toeplitz 矩阵构成，通过牛顿–拉夫逊迭代求解。该方法在弱非线性条件下具有谱精度。

4.4. Advanced Techniques Based on Circuit-Block Partitioning

4.4.1. Time-Domain Latency Envelope Following
将电路变量分为活跃（快变）和潜在（慢变）子集，在快时间方向上采用多速率 Runge–Kutta 方案（如基于 Bogacki–Shampine 的公式），活跃变量用小步长积分，潜在变量用大步长积分，从而显著降低计算量。该技术比标准时域包络跟随更高效，且当潜在块大于活跃块时优势更明显。

4.4.2. Mixed Time–Frequency Latency Envelope Following
在慢时间方向上采用 HB，同时利用变量活跃/潜在划分：活跃变量以 K 个谐波表示，潜在变量仅需直流分量（一个系数）。这大幅缩减了 HB 系统的维度和 Jacobian 矩阵尺寸，减少了时频域转换次数。该技术比混合时频包络跟随更具效率，尤其适用于潜在变量数量远大于活跃变量的场景。

4.5. Order and Stability

4.5.1. The BDF2 (Gear-2) Scheme: Order and Stability
BDF2 是二阶隐式多步法，具有 A-稳定性（覆盖整个左半平面）和 L-稳定性，能快速阻尼高频振荡，适合刚性系统。

4.5.2. Harmonic Balance: Spectral Accuracy and Numerical Robustness
HB 具有谱精度，截断误差随谐波数 K 指数衰减。数值稳定性取决于多频 Jacobian 矩阵的条件数和混叠误差控制；在强非线性下，收敛速率可能降至代数级。

4.5.3. Multirate Runge–Kutta Schemes: Order and Stability Regions
多速率 RK 方案基于 Bogacki–Shampine 方法，需满足耦合条件以保持三阶精度。稳定性分析通过耦合刚度测试方程，绘制稳定区域图，方法在耦合参数和刚度比下保持稳定。

5. Numerical Results

5.1. Highly Nonlinear Illustrative Example
以 RF 极性功率放大器为例，比较时域包络跟随（TDEF）和时域延迟包络跟随（TDLEF）与单变量 Gear-2 方法。混合时频方法因谐波过多而失效。结果显示 TDEF 和 TDLEF 的 CPU 时间远低于单变量方法，且 TDLEF 因块划分进一步加速，精度损失极小（归一化均方误差约 10^?8）。

5.2. Moderately Nonlinear Illustrative Example
以 RF 混频器为例，比较 TDEF、混合时频包络跟随（HTFEF）、TDLEF 和混合时频延迟包络跟随（MTFLEF）与单变量参考方法。HTFEF 和 MTFLEF 优于纯时域方法，且块划分技术（TDLEF vs TDEF, MTFLEF vs HTFEF）均提升效率，最大误差保持在 10^?7 以下。

6. Discussion
综述总结了 MPDAE 方法的优势：可大幅加速多速率、多尺度电路仿真；时域方法适用于强非线性，频域/混合方法适用于中度非线性；电路块划分进一步利用变量时间尺度差异提高效率。主要挑战包括：MPDAE 计算复杂度随时间尺度数量指数增长，限制了超过两三个尺度的应用；以及处理非周期多速率激励的数值策略尚待开发。