相位-频率传递模型

相位-频率传递 (Phase-Frequency Transfer) 模型

相位-频率传递 (Phase-Frequency Transfer, PFT) 模型 是 Martens & Gielen《High-Level Modeling and Synthesis of Analog Integrated Systems》(2008) Ch.5 提出的全新频域建模框架。核心创新是用 Polyphase Harmonic Transfer Matrix (PHTM) 统一表示多相（多分支）RF 前端的线性/非线性行为，支持下变频/上变频操作，并自然处理频率偏移和非理性效应。

---

动机与需求

传统的 RF 前端模型使用等效低通(equivalent low-pass)方法，对窄带系统有效但对宽带/多频系统有误差。PFT 模型的设计需求（§5.3.1）：

1. 允许多相(多分支)信号表示 — RF 中常用 I/Q、多相滤波
2. 支持下变频/上变频操作 — 混频器核心功能
3. 自然处理频率偏移 — 本振产生的频谱搬移
4. 支持级联、并联、反馈互联 — 常见架构组合方式
5. 可扩展至弱非线性 — 覆盖交叉调制和互调产物

---

Polyphase Harmonic Signal (PHS)

定义

$N$ 相多相信号表示为向量：

$$s(t)=[s_0(t),s_1(t),...,s_{N-1}(t){]}^T$$

多相初等信号 (Polyphase Elementary Signals)

模式	信号	用途
差分模式	$s=[1,-1,1,-1,...{]}^T$	相邻反相
共模	$s=[1,1,...,1{]}^T$	同相、偏置
正序	$s_k=e^{j2\pi k/N}$	顺时针旋转相量
负序	$s_k=e^{-j2\pi k/N}$	逆时针旋转相量

PHS 谱表示

Polyphase Harmonic Signal 的频谱表示为（Fig. 5.2）：

$$S(f)=\sum _{a=0}^{A-1}\sum _{n=0}^{N-1}{s}_{a,n}\cdot \delta (f-f_a-n\cdot f_s/N)$$

其中 $N$ 相位数，$f_s$ 采样频率，$A$ 频率数，$s_{a,n}$ 为各频率-相位分量。

基函数：PHS 由 $(N\times A)$ 个基函数张成，每个基对应一个特定相位模式和特定频率。

---

Polyphase Harmonic Transfer Matrix (PHTM)

定义

PHTM 为线性变换算子，将输入多相信号的相位-频率平面映射到输出多相信号的相位-频率平面：

$${\tilde{H}}_L(f)=\left[\begin{array}{cccc}{H}_L^{0,0}(f)& H_L^{0,1}(f)& \cdots & H_L^{0,M-1}(f)\\ H_L^{1,0}(f)& H_L^{1,1}(f)& \cdots & H_L^{1,M-1}(f)\\ ⋮& ⋮& \ddots & ⋮\\ H_L^{N-1,0}(f)& H_L^{N-1,1}(f)& \cdots & H_L^{N-1,M-1}(f)\end{array}\right]$$

其中：

• $M$ 为输入相位数，$N$ 为输出相位数
• $H_L^{i,j}(f)$ 表示从输入第 $j$ 相位—频率平面到输出第 $i$ 相位—频率平面的传递函数
• 频率维度隐含在 $f$ 参数中

Box Plot 可视化（Fig. 5.8-5.9）

PHTM 可以可视化为 $N\times M$ 的网格，每个网格元素 $H_L^{i,j}(f)$ 为一复数传递函数：

输出相位 i →
   H^0,0(f)  H^0,1(f)  ...  H^0,M-1(f)
   H^1,0(f)  H^1,1(f)  ...  H^1,M-1(f)
   ...       ...       ...  ...
   H^N-1,0(f) H^N-1,1(f) ... H^N-1,M-1(f)

Box Plot 用灰度表示传递幅度（暗=大传递），简明显示系统中不同相位模式间的耦合强度。

相位变换器 (Phase-converter)（Fig. 5.3）

特殊类型 PHTM，完成相位模式间的转换：

• 差分→正序：将 $[1,-1{]}^T$ 转换为 $[1,j{]}^T$
• 正序→差分：将 $[1,j{]}^T$ 转换为 $[1,-1{]}^T$
• N→M 转换：转换相位模式维度（完成匹配/不匹配转换）

三种基本连接类型（Fig. 5.7）

PHTM 支持代数组合，类似线性系统块图：

1. 级联 (Cascade)：${\tilde{H}}_{out}(f)={\tilde{H}}_2(f)\cdot {\tilde{H}}_1(f)\cdot {\tilde{H}}_{in}(f)$
   • 矩阵乘法，注意维度兼容
2. 并联 (Parallel)：${\tilde{H}}_{out}(f)={\tilde{H}}_1(f)+{\tilde{H}}_2(f)$
   • 矩阵加法，要求维度相同
3. 反馈 (Feedback)：${\tilde{H}}_{cl}(f)=(I-{\tilde{H}}_1(f)\cdot {\tilde{H}}_2(f){)}^{-1}\cdot {\tilde{H}}_1(f)$
   • 类似反馈系统传递函数，$I$ 为单位矩阵

---

噪声模型（§5.3.4, Fig. 5.10）

PHTM 框架中噪声建模：

• 多相噪声向量 $\mathbf{n}(t)$
• 输入噪声 ${\mathbf{n}}_{in}$ 通过 PHTM 传递：

- 支持周期平稳噪声：噪声谱密度在相位-频率平面中随相位模式变化 --- ## 非线性扩展：Distortion Tensor（§5.3.5） ### 定义 $K$ 阶失真张量 (Distortion Tensor of Order $K$)：

\tilde{H}_{NL}^{(K)}(f_1, …, f_K): \mathbb{C}^{N \times 1} \times \cdots \times \mathbb{C}^{N \times 1} \to \mathbb{C}^{M \times 1}

将 $K$ 个输入多相信号映射到输出，完整非线性行为为各阶贡献之和（Fig. 5.11）：

\tilde{Y}(f) = \sum_{K=1}^\infty \int_{-\infty}^{\infty} \cdots \int_{-\infty}^{\infty} \tilde{H}_{NL}^{(K)}(f_1, …, f_K) \cdot \tilde{X}(f_1) \cdots \tilde{X}(f_K) \cdot \delta(f - \sum f_i) , df_1 \cdots df_K

其中： - $K=1$：线性项 — 对应 $H_L(f)$ - $K=2$：二阶失真（HD2、IM2） - $K=3$：三阶失真（HD3、IM3） ### 三阶非线性实例（Fig. 5.12） 当设置 $K_3=1$，三阶失真张量可展开为三维矩阵 Box Plot，显示： - 频谱分量间的耦合：$f_1 + f_2 + f_3$ 产生新的频率分量 - 镜像频率折叠：干扰信号通过三阶非线性折叠到目标频带 - 相位模式混叠：不同相位模式的输入产生交叉相位输出 ### 与 Volterra 级数的关系 Distortion Tensor 是 **Volterra 级数在 Polyphase 域的多维扩展**： - Volterra 核 $h_K(t_1, ..., t_K)$ → $H_K(f_1, ..., f_K)$ → $\tilde{H}_{NL}^{(K)}(f_1, ..., f_K)$ - 额外维度来源于相位模式间耦合 - 复杂度：$K$ 阶张量的维度为 $M \times N^K$（线性时 $K=1$ 为 $M \times N$） --- ## RF 前端通用行为模型（§5.4） PFT 模型作为 RF 前端通用行为模型的数学基础（Fig. 5.14）： **RF 前端标准分解**： ``` LNA (PHTM_LNA) → Mixer (PHTM_mix) → Filter (PHTM_filt) → Sampler (PHTM_samp) ``` **性能评估实例**： - **IRR (Image Rejection Ratio)**：通过混频器 PHTM 计算，增益失配 $\sigma_{gain}$ 影响 IRR（Fig. 5.18） - **NF (Noise Figure)**：结合各模块 PHTM 和噪声模型计算 - **非线性产物**：Distortion Tensor 计算 IM3/HD3 和频谱再生（Fig. 5.20-5.21） --- ## 核心优势总结 | 方面 | 传统方法 | PFT 模型 | |------|---------|----------| | 多相信号 | 独立处理各分支 | 统一向量表示 | | 频率搬移 | 需等效低通近似 | 自然包含于 PHTM | | 镜像折叠 | 需单独分析 | PHTM 自动计算 | | 非线性 | Volterra 级数（无多相） | Distortion Tensor（含多相） | | 互联组合 | 手动推导 | 矩阵代数自动组合 | | 计算效率 | 时域仿真慢 | 频域直接计算快 | --- ## 与其他概念的关系 - [[通用行为建模方法]] — PFT 是频域 GBM 的核心 - [[描述与抽象层设计框架]] — PFT 位于行为描述层 - [[自顶向下异构综合优化]] — PFT 模型可嵌入性能计算器 - [[CT-ΔΣ调制]] — 对比：时域 GBM 用于 ΔΣ，频域 GBM/PFT 用于 RF