雷达 · 通信 · 电子战 / 基础精学 001
一条波形,一个相量,一组相位递推
在通信和雷达处理中,同一个正弦信号常有三种写法:时域里的余弦曲线、复平面上的旋转相量,以及数字电路中的相位累加器。把这三种表示对应起来,I/Q、相位补偿和数字本振就有了共同的数学起点。
核心关系:正弦波是旋转相量在实轴上的投影;频率决定相位增长的速度,初始相位决定起点。
一个连续时间正弦信号可以写成:
A 为正的峰值振幅,f 的单位为 Hz,φ 的单位为 rad;本文取 A > 0、f > 0。
振幅决定波形有多高,频率决定每秒重复多少次,相位决定观察开始时信号处于周期的哪个位置。角频率 ω 表示每秒转过的弧度数,与频率 f 相差一个 2π。
例如,A = 1、f = 1 kHz、φ = π/4 时,周期为 1 ms,初始相位为 45°。在 t = 0 时,信号已经处于 cos 45° ≈ 0.7071 的位置,并不从零开始。
相位比较需要共同的时间参考。若两个信号的频率不同,它们的相位差会随时间变化;只有同频时,相位差才可以保持为常数。
欧拉公式把三角函数与复指数连接起来。工程中通常用 j 表示虚数单位,满足 j² = −1:
Re{·} 表示取实部。
z(t) 可以看成复平面上一根长度为 A 的向量。时间增加时,它以角速度 ω 转动;向量在水平轴上的投影,正是余弦波。频率越高,转动越快;振幅越大,向量越长。
同样,由正转和反转两个复指数相加,可以得到实余弦:
这也解释了实余弦为什么同时具有正、负频率分量。负频率在复数表示中对应反向旋转,是描述信号的数学组成部分。
把复信号的实部记作 I,虚部记作 Q,就得到同相与正交分量(In-phase / Quadrature,I/Q):
图 1 相量的水平、垂直投影分别形成 I 和 Q。算例初始相位为 45°,两个初始分量均约为 0.7071。
两个坐标共同确定复平面上的位置。对本文的非零单频复信号,令 θ(t) = ωt + φ,可由 I、Q 求得幅度与主值相位 θp:
atan2 区分四个象限,返回 [−π, π] 范围内的主值相位;θp 与 θ 按模 2π 等价。
主值相位不记录已转过的整周数,不能直接当作随时间连续增长的相位。连续相位需要在采样相位变化可判别等条件下进行相位展开;当 I = Q = 0 时,幅度为零,物理相位没有定义。
若只知道某一时刻的 I 值,通常无法区分相位 θ 与 −θ,因为 cos θ = cos(−θ)。保留 Q 后,两个位置就能区分。对于随时间变化的信号,复数表示也便于统一完成混频、滤波和相位旋转。
符号约定:本文采用 z = I + jQ。有些接收机采用相反的 Q 极性;处理链中的约定必须一致,否则相位和频率方向的解释会发生变化。
将信号延迟 τ,意味着把时间变量 t 换成 t − τ。代入后可直接看到相位变化:
在本文的正频率和复指数约定下,正时延引起负相位变化,也就是相位滞后。以 1 kHz 信号为例,延迟 125 μs,相当于延迟了八分之一个周期:
图 2 橙色虚线相对蓝色实线向右移动 125 μs,对应同频信号的 45° 相位滞后。
这一关系有两个重要前提。其一,单频相位仅在模 2π 的意义下确定:对单频信号,τ 与 τ + T 产生相同的相位结果,单频相位不能唯一确定绝对时延。
其二,同一个时延对不同频率产生不同相移。多频或宽带信号的真实延迟,对应随频率线性变化的相位;一个对所有频率施加相同角度的移相器,不能普遍替代真实时延。
采样率为 Fs 时,第 n 个样本对应 t = n/Fs。连续时间模型变为:
Ω 是离散时间角频率,单位为 rad/sample。
用 16 kHz 的采样率采集 1 kHz 正弦信号,每个周期得到 16 个样本,相邻样本的相位增加 π/8,也就是 22.5°。
图 3 连续曲线与离散样本的对应关系。横轴保留实际时间,蓝色离散点每隔 62.5 μs 取一个样本。
频率 f 描述每秒转几圈,Ω 描述每个样本转多少弧度。比较不同采样率的系统时,需要同时注明实际频率和归一化频率;相同的离散序列,在不同采样时钟下可以对应不同的实际频率。
数控振荡器(Numerically Controlled Oscillator,NCO)可以用整数相位累加器实现等步进旋转。设累加器有 W 位,相位步进字为 K:
K 为整数,Fs 为相位更新率;这里取 0 < K < 2W−1,输出频率低于奈奎斯特频率。
每次加 K,就让相量向前转一小步;累加器溢出回绕,相当于相位绕完一周。将相位映射到余弦查找表,就得到离散波形。
当 W = 16、K = 4096 时,每个采样周期前进一周的 1/16。若时钟为 16 kHz,输出就是 1 kHz;若时钟为 160 MHz,同样的整数步进对应 10 MHz。
上述步进字决定频率,初相位则由累加器初值决定。为复现前文 φ = π/4 的波形,还需设置 p[0] = 216φ/(2π) = 8192,并以更新前的 p[n] 生成第 n 个样本。相位映射这里按理想余弦值计算;实际查表和幅度量化会引入误差。
三个位宽的作用需要分开:累加器位宽决定频率步进分辨率,查表地址位宽决定相位网格,输出位宽决定幅度量化精度。增加其中一个位宽,不能自动消除另外两类误差。
以下结果沿用第 1 课固定参数:A = 1、f = 1 kHz、Fs = 16 kHz、φ = π/4、τ = 125 μs。原课 MATLAB 数值结果与解析值一致;本次图文版按相同参数重新绘制了示意图。
| 检查项 | 解析关系 | 数值 |
|---|---|---|
| 周期 | 1/f | 1 ms |
| 每周期样本数 | Fs/f | 16 |
| 每样本相位增量 | 2πf/Fs | π/8 |
| 延迟相移 | −2πfτ | −45° |
| 初始 I / Q | cos φ / sin φ | 均约 0.7071 |
三种表示都在描述同一件事:相位如何随时间或采样序号增长。连续时间用 ωt,离散时间用 Ωn,数字电路用相位累加值 p[n]。把单位、正负号和参考时刻对齐,三者就可以互相核对。
从正弦波走向 I/Q、数字混频和相干处理,首先需要确认的始终是:相位相对于谁定义,每次更新又前进了多少。
内容依据:每日精学第 01 课《正弦信号、复指数与相位》及配套 MATLAB 结果。图为固定参数的数值示意,不代表设备实测。