4. 电路仿真基础与射频前端行为建模

射频前端同时包含线性无源网络与强非线性有源器件：滤波器、双工器、匹配网络、传输线属于线性网络，功率放大器、混频器、开关、LNA 的压缩与交调属于非线性行为。若只用一种仿真引擎，要么无法处理非线性失真，要么在长调制包络上耗时过长。因此本节的核心不是“学会点哪个按钮”，而是理解四类方法各自回答什么问题、在什么条件下成立，以及如何把它们串成一条可信的仿真链。

一、四类仿真方法的分工

S 参数级联处理的是线性、时不变、小信号条件下的网络组合。每个二端口用散射矩阵描述入射波与反射波的关系，级联时把各段 S 矩阵转换为传输矩阵（T 矩阵或 ABCD 矩阵）相乘，再转回 S 参数。它回答的是：匹配网络加滤波器加封装寄生后，整体插损、回波损耗和群时延是多少。前提是各端口阻抗一致（通常 50 Ω）、网络线性、激励足够小。

谐波平衡（HB）处理周期稳态下的强非线性。它把时域波形展开为有限次谐波的频域未知量，在时域计算非线性电流、在频域施加线性网络约束，通过牛顿迭代使两者自洽。它回答的是：给定输入功率和频率，PA 的 1 dB 压缩点、二次/三次谐波、双音交调 IM3 落在哪里。HB 的代价是谐波数越多、非线性阶数越高，雅可比矩阵越大。

包络仿真处理“高频载波 + 慢速调制”的双时间尺度问题。载波仍用 HB 或谐波域描述，包络（幅度与相位随时间的慢变化）用时间步进积分。它回答的是：WCDMA/LTE 信号激励下，PA 的邻道功率比 ACLR、记忆效应引起的 AM-PM 随包络变化如何。相比纯时域 SPICE 逐周期积分，包络仿真把步长从载波周期放宽到包络带宽尺度，速度提升常达数十到数百倍。

非线性器件模型是上述三种方法的共同底座。模型不准，HB 和包络仿真给出的压缩点与交调就是“精确地算错”。

二、S 参数级联的实际做法

以一条接收链为例：开关（插损 0.5 dB）→ 滤波器（插损 1.8 dB）→ LNA（增益 15 dB、噪声系数 1.2 dB）→ 匹配网络（插损 0.4 dB）。级联计算的关键是先算前级、再算后级，因为后级噪声会被前级增益“稀释”。用 Friis 公式：

F_total = F1 + (F2-1)/G1 + (F3-1)/(G1·G2) + ...
线性值：F1=1.122, G1=0.891
        F2=1.514, G2=0.661
        F3=15.85, G3=31.62
        F4=1.096
F_total = 1.122 + (1.514-1)/0.891 + (15.85-1)/(0.891×0.661)
        + (1.096-1)/(0.891×0.661×31.62)
        ≈ 1.122 + 0.577 + 25.22 + 0.005 ≈ 26.92
NF_total ≈ 10·log10(26.92) ≈ 14.3 dB

这个结果说明：滤波器插损 1.8 dB 直接恶化噪声系数，而 LNA 之后的匹配网络插损几乎不影响总 NF。若把 LNA 前移，NF 可降到约 2.5 dB。级联时还要注意：每段 S 参数必须在同一参考阻抗下测得，否则需先做阻抗重归一化；有源器件在小信号下才可用 S 参数，一旦进入压缩区，级联结果不再成立。

三、谐波平衡的输入与判读

HB 的典型设置包括：基波频率 f0、谐波阶数（如 5 阶）、本振或双音间隔 Δf、输入功率扫描范围。以双音测试为例，f1=2.40 GHz、f2=2.41 GHz，Δf=10 MHz，输入双音总功率从 -10 dBm 扫到 +10 dBm。仿真输出的是各阶交调分量：IM3 落在 2f1-f2=2.39 GHz 和 2f2-f1=2.42 GHz。判读时看两条曲线：基波输出随输入线性上升（斜率 1），IM3 随输入以斜率 3 上升，两者交点外推即为三阶截点 OIP3。若仿真中 IM3 曲线斜率明显偏离 3，通常意味着谐波阶数不足或模型外推越界。

四、包络仿真与记忆效应

调制信号下，PA 的增益和相位不仅取决于当前包络幅度，还受历史包络影响，这就是记忆效应。包络仿真把载波周期与包络周期分开：载波用谐波域描述，包络用时间步进。设置时需给出包络采样率（通常为信号带宽的 4～8 倍）、载波频率、以及调制信号的 I/Q 数据。输出 ACLR 时，观察主信道两侧相邻信道的积分功率比；输出 AM-AM/AM-PM 曲线时，观察包络幅度变化引起的增益压缩与相位旋转。若 ACLR 在包络峰值处急剧恶化，说明 PA 已接近饱和，需要回退功率或调整匹配。

五、非线性器件模型的选用

模型选择取决于器件类型与工作区。GaN HEMT 常用 Angelov 或 EEHEMT 模型，前者对跨导和栅电容的拟合较平滑，适合 HB 与包络仿真；GaAs FET 可用 Curtice 立方模型；LDMOS 常用 MET 模型。模型参数通常来自 IV 曲线、S 参数和负载牵引数据的分步提取。需要特别注意的是：模型的有效区。若器件实际工作在膝点电压附近或深饱和区，而模型参数只拟合了线性区，HB 给出的 P1dB 会偏乐观。工程上常用“模型验证三件套”：直流 IV 对比、小信号 S 参数对比、单音功率扫描对比，三者都吻合才用于系统级仿真。

六、一条可复现的仿真链

综合以上，推荐按以下顺序搭建射频前端仿真：先用 S 参数级联完成线性链路预算，确定各级增益、插损与噪声系数；再把 PA 单独放入 HB，用双音扫描提取 P1dB 与 OIP3；随后用包络仿真加载实际调制信号，评估 ACLR 与 EVM；最后把 PA 的行为模型（如记忆多项式）回嵌到系统级链路，验证整机指标。每一步都要保留输入条件、模型版本和输出曲线，便于回溯。这样得到的仿真结果才既有物理依据，又能支撑工程决策。