OFDM 系统中常见星座图畸变机理分析与数学建模(含 CFO 载波频偏,SFO 采样频偏,CTO 粗定时偏移,STO 细定时偏移,IQ 不平衡情况讨论)

前言

在数字通信系统中,星座图以直观方式呈现调制符号在复平面上的统计分布,其形态不仅反映了调制方式与噪声水平,也集中体现了同步误差、射频前端非理想与信道畸变等因素对接收基带信号的综合影响。对于正交频分复用(Orthogonal Frequency Division Multiplexing, OFDM)系统而言,星座图畸变具有更强的“结构性”:由于多载波正交性的存在,频率同步、采样同步与定时同步的细微偏差均可能诱发子载波间干扰(Inter-Carrier Interference, ICI)或相位误差累积,从而显著劣化误码率与误差矢量幅度(Error Vector Magnitude, EVM)。因此,围绕星座图畸变的机理进行系统化建模与分析,不仅有助于解释仿真或实测现象,也为同步算法、补偿结构与射频校准策略的设计提供理论依据。

本文以 AWGN 信道下的 OFDM 基带模型为主线,依次讨论理想加性高斯白噪声(Additive White Gaussian Noise, AWGN)条件下的星座分布特性,并在此基础上引入并分析若干最常见且对 OFDM 影响显著的非理想因素,包括载波频偏(Carrier Frequency Offset, CFO)、采样频偏(Sampling Frequency Offset, SFO)、粗定时偏移(Coarse Timing Offset, CTO)、细定时偏移(Symbol/Fine Timing Offset, STO)以及 IQ 不平衡(I/Q Imbalance)。行文结构上,本文首先给出 OFDM 系统的时域与频域数学描述以及 AWGN 条件下的接收模型,随后分别针对上述各类非理想因素建立等效基带模型,推导其对频域子载波符号的影响形式,并阐明由此导致的星座图畸变特征及其系统危害,最后对全文进行总结归纳。

1 仅 AWGN(理想)条件下的 OFDM 星座模型

考虑一个基带等效的 OFDM 系统,其每个 OFDM 符号由 N N N 点 IFFT 生成,子载波间隔为 Δ f = F s N \Delta f = \frac{F_s}{N} Δf=NFs,其中 F s F_s Fs 为采样率。设第 m m m 个 OFDM 符号在频域上第 k k k 个子载波 ( k ∈ K ) (k\in\mathcal{K}) (kK)承载的调制符号为 X m [ k ] X_m[k] Xm[k],其中 K \mathcal{K} K 表示有效子载波集合, X m [ k ] X_m[k] Xm[k] 为复数星座点(例如 M − Q A M M-QAM MQAM),并满足 E ∣ X m [ k ] ∣ 2 = P X \mathbb{E}{|X_m[k]|^2}=P_X EXm[k]2=PX。将频域符号映射到长度为 N N N 的频域向量 X m \mathbf{X}_m Xm,其第 k k k 个分量为 X m [ k ] X_m[k] Xm[k],未使用子载波置零。忽略窗函数与脉冲成形影响时,第 m m m 个 OFDM 符号的时域采样序列可写为
x m [ n ] = 1 N ∑ k = 0 N − 1 X m [ k ] e j 2 π N k n , n = 0 , 1 , … , N − 1 , x_m[n]=\frac{1}{\sqrt{N}}\sum_{k=0}^{N-1} X_m[k] e^{j\frac{2\pi}{N}kn}, \quad n=0,1,\dots,N-1, xm[n]=N 1k=0N1Xm[k]ejN2πkn,n=0,1,,N1,
其中 n n n 表示离散时间采样索引, 1 N \frac{1}{\sqrt{N}} N 1 用以保证能量归一化。为抵抗多径引起的符号间干扰,OFDM 引入循环前缀(Cyclic Prefix, CP),设 CP 长度为 N cp N_{\text{cp}} Ncp,则带 CP 的发送序列为
x ~ m [ n ] = x m [ ( n   m o d   N ) ] , n = − N cp , … , N − 1. \tilde{x}_m[n]=x_m[(n \bmod N)],\quad n=-N_{\text{cp}},\dots,N-1. x~m[n]=xm[(nmodN)],n=Ncp,,N1.
在纯 AWGN 信道下,假定无载波与采样同步误差且定时对齐准确,接收端采样得到的序列可表示为
r ~ m [ n ] = x ~ m [ n ] + w m [ n ] , \tilde{r}_m[n]=\tilde{x}_m[n]+w_m[n], r~m[n]=x~m[n]+wm[n],
其中 w m [ n ] w_m[n] wm[n]为复高斯白噪声序列,满足 w m [ n ] ∼ C N ( 0 , σ w 2 ) w_m[n]\sim\mathcal{CN}(0,\sigma_w^2) wm[n]CN(0,σw2),且在不同 n n n 上独立同分布。接收端去除 CP 后得到
r m [ n ] = x m [ n ] + w m [ n ] , n = 0 , 1 , … , N − 1. r_m[n]=x_m[n]+w_m[n],\quad n=0,1,\dots,N-1. rm[n]=xm[n]+wm[n],n=0,1,,N1.
r m [ n ] r_m[n] rm[n] 进行 N N N 点 FFT 得到频域观测
Y m [ k ] = 1 N ∑ n = 0 N − 1 r m [ n ] e − j 2 π N k n = X m [ k ] + W m [ k ] , Y_m[k]=\frac{1}{\sqrt{N}}\sum_{n=0}^{N-1} r_m[n] e^{-j\frac{2\pi}{N}kn} =X_m[k]+W_m[k], Ym[k]=N 1n=0N1rm[n]ejN2πkn=Xm[k]+Wm[k],
其中
W m [ k ] = 1 N ∑ n = 0 N − 1 w m [ n ] e − j 2 π N k n W_m[k]=\frac{1}{\sqrt{N}}\sum_{n=0}^{N-1} w_m[n] e^{-j\frac{2\pi}{N}kn} Wm[k]=N 1n=0N1wm[n]ejN2πkn
仍为零均值复高斯随机变量,且 W m [ k ] ∼ C N ( 0 , σ w 2 ) W_m[k]\sim\mathcal{CN}(0,\sigma_w^2) Wm[k]CN(0,σw2)(在单位ary 变换下方差保持不变)。因此,在理想 AWGN 条件下,星座图可视为以 X m [ k ] X_m[k] Xm[k] 为中心、方差由 σ w 2 \sigma_w^2 σw2 决定的二维圆对称高斯散布,其畸变不具有结构性,主要由信噪比 SNR = P X / σ w 2 \text{SNR}=P_X/\sigma_w^2 SNR=PX/σw2 控制。
在这里插入图片描述

2 CFO 载波频偏的模型与危害

载波频偏(CFO)来源于发射端与接收端本振频率不一致或多普勒频移等因素。设发射端载波频率为 f c f_c fc,接收端下变频本振为 f c + Δ f c f_c+\Delta f_c fc+Δfc,其中 Δ f c \Delta f_c Δfc 为频偏,则理想基带等效下,接收采样序列将叠加一个随时间线性增长的相位项。若以采样周期 T s = 1 / F s T_s=1/F_s Ts=1/Fs 表示离散时间,则 CFO 引起的复包络旋转可写为
r ~ m [ n ] = x ~ m [ n ] e j 2 π Δ f c ( n T s ) + w m [ n ] = x ~ m [ n ] e j 2 π N ε n + w m [ n ] , \tilde{r}_m[n]=\tilde{x}_m[n] e^{j2\pi \Delta f_c (nT_s)} + w_m[n] =\tilde{x}_m[n] e^{j\frac{2\pi}{N}\varepsilon n}+w_m[n], r~m[n]=x~m[n]ej2πΔfc(nTs)+wm[n]=x~m[n]ejN2πεn+wm[n],
其中引入归一化频偏 ε = Δ f c / Δ f \varepsilon=\Delta f_c/\Delta f ε=Δfcf,而 Δ f = F s / N \Delta f=F_s/N Δf=Fs/N 为子载波间隔。去 CP 并取 FFT 后,第 m m m 个 OFDM 符号频域输出可表示为
Y m [ k ] = ∑ ℓ = 0 N − 1 X m [ ℓ ] C ε ( k − ℓ ) e j ϕ m + W m [ k ] , Y_m[k]=\sum_{\ell=0}^{N-1} X_m[\ell] C_{\varepsilon}(k-\ell) e^{j\phi_m} + W_m[k], Ym[k]==0N1Xm[]Cε(k)ejϕm+Wm[k],
其中 C ε ( ⋅ ) C_{\varepsilon}(\cdot) Cε() 为由 CFO 导致的 ICI 系数,通常可写为
C ε ( p ) = 1 N ∑ n = 0 N − 1 e j 2 π N ( ε − p ) n = 1 N e j π ( ε − p ) N − 1 N sin ⁡ ( π ( ε − p ) ) sin ⁡ ( π ( ε − p ) / N ) , C_{\varepsilon}(p)=\frac{1}{N}\sum_{n=0}^{N-1} e^{j\frac{2\pi}{N}(\varepsilon-p)n} =\frac{1}{N}e^{j\pi(\varepsilon-p)\frac{N-1}{N}}\frac{\sin\big(\pi(\varepsilon-p)\big)}{\sin\big(\pi(\varepsilon-p)/N\big)}, Cε(p)=N1n=0N1ejN2π(εp)n=N1e(εp)NN1sin(π(εp)/N)sin(π(εp)),
而公共相位项
ϕ m = 2 π N ε ( m ( N + N cp ) + N cp ) \phi_m = \frac{2\pi}{N}\varepsilon (m(N+N_{\text{cp}})+N_{\text{cp}}) ϕm=N2πε(m(N+Ncp)+Ncp)
描述了 CFO 在不同 OFDM 符号之间累积的相位旋转。由上式可见,当 ε ≠ 0 \varepsilon\neq 0 ε=0 时, C ε ( 0 ) C_{\varepsilon}(0) Cε(0) 不再等于 1 且 C ε ( p ≠ 0 ) C_{\varepsilon}(p\neq 0) Cε(p=0) 不再为 0,导致期望子载波分量 X m [ k ] X_m[k] Xm[k] 的幅度衰减并伴随来自其他子载波的泄漏,即 ICI。该效应在星座图上体现为两类叠加畸变:其一为符号级公共相位误差(Common Phase Error, CPE)使星座整体随 OFDM 符号索引旋转;其二为 ICI 使星座点呈现扩散与模糊化,且扩散程度与 ∣ ε ∣ |\varepsilon| ε 增大显著增强。系统层面上,CFO 会破坏 OFDM 子载波正交性,造成有效信噪比降低与误码率劣化,并且在高阶 QAM 下更易表现为 EVM 超标与误码地板。
在这里插入图片描述

3 SFO 采样频偏的模型与危害

采样频偏(SFO)是指接收端采样时钟与发射端采样时钟存在微小比例误差,典型来源为晶振精度不足、温漂或时钟分配不一致。设发射端采样率为 F s F_s Fs,接收端采样率为 F s ( 1 + ρ ) F_s(1+\rho) Fs(1+ρ),其中 ρ \rho ρ 为无量纲的相对采样频偏(通常以 ppm 表示)。在连续时间层面,接收机在时刻 t n = n / F s ( 1 + ρ ) t_n=n/F_s(1+\rho) tn=n/Fs(1+ρ) 对信号采样,相当于对理想离散序列在非整数索引处重采样。若忽略插值误差,将 SFO 视为时间轴伸缩,则可表示为
r ~ m [ n ] ≈ x ~ ( n 1 + ρ ) + w [ n ] , \tilde{r}_m[n] \approx \tilde{x}\left(\frac{n}{1+\rho}\right) + w[n], r~m[n]x~(1+ρn)+w[n],
其中 x ~ ( ⋅ ) \tilde{x}(\cdot) x~()表示连续时间或带限重构意义下的发送信号。对 OFDM 而言,SFO 的主要危害在于其引入的相位误差同时依赖符号索引 m m m 与子载波索引 k k k。在常用的小偏差近似 ∣ ρ ∣ ≪ 1 |\rho|\ll 1 ρ1 下,SFO 在频域可近似导致第 m m m 个 OFDM 符号的第 k k k 个子载波出现附加相位项
Y m [ k ] ≈ X m [ k ] e j 2 π k ρ m ( N + N cp ) N + ICI ρ ( m , k ) + W m [ k ] , Y_m[k]\approx X_m[k] e^{j 2\pi k \rho \frac{m(N+N_{\text{cp}})}{N}} + \text{ICI}_{\rho}(m,k) + W_m[k], Ym[k]Xm[k]ej2πkρNm(N+Ncp)+ICIρ(m,k)+Wm[k],
其中指数项体现了相位随符号索引 m m m 线性累积、并随子载波索引 k k k 成比例放大的特性, ICI ρ ( m , k ) \text{ICI}_{\rho}(m,k) ICIρ(m,k) 表示由时间伸缩破坏正交性引起的 ICI 残差。与 CFO 相比,SFO 的突出特征在于其相位误差不仅在时间上累积,而且在频率维度上呈线性斜坡,使不同子载波的星座旋转角不同,进而在汇聚绘制时表现为更明显的环状或带状结构。系统危害方面,SFO 会导致跨符号的相位漂移与子载波间干扰逐步增强,使得导频辅助的信道估计与相位跟踪难以保持一致性,并在长帧传输、窄子载波间隔或高阶调制条件下显著劣化性能。
在这里插入图片描述

4 CTO 粗定时偏移的模型与危害

粗定时偏移(CTO)通常指接收端在帧同步或符号定界时产生的整数采样级偏移,即 FFT 窗口的起始位置相对于理想位置偏移 Δ n ∈ Z \Delta n\in\mathbb{Z} ΔnZ。在 OFDM 系统中,若 Δ n \Delta n Δn仍落在循环前缀覆盖范围内,即 ∣ Δ n ∣ ≤ N cp |\Delta n|\le N_{\text{cp}} ∣ΔnNcp,则不会引入来自相邻 OFDM 符号的符号间干扰,但会在频域造成确定性的线性相位斜坡。设理想情况下去 CP 后的序列为 x m [ n ] x_m[n] xm[n],若 FFT 窗起点相对理想位置提前 Δ n > 0 \Delta n>0 Δn>0个采样,则实际输入 FFT 的序列可记为
r m [ n ] = x m [ ( n − Δ n )   m o d   N ] + w m [ n ] , n = 0 , 1 , … , N − 1. r_m[n]=x_m[(n-\Delta n)\bmod N] + w_m[n],\quad n=0,1,\dots,N-1. rm[n]=xm[(nΔn)modN]+wm[n],n=0,1,,N1.
其 FFT 输出为
Y m [ k ] = X m [ k ] e − j 2 π N k Δ n + W m [ k ] , Y_m[k]=X_m[k] e^{-j\frac{2\pi}{N}k\Delta n}+W_m[k], Ym[k]=Xm[k]ejN2πkΔn+Wm[k],
其中 Δ n \Delta n Δn 的整数性使相位项在每个子载波上为固定旋转,且旋转角随 k k k 线性变化。由此可见,在 CP 内 CTO 不会产生 ICI,但会造成子载波间的相位一致性被破坏,尤其在需要多子载波联合处理的场景(如信道估计插值、MIMO 预编码与解码)中,该线性相位斜坡若不被补偿将导致等效信道相位响应产生系统性偏移。若 ∣ Δ n ∣ > N cp |\Delta n|>N_{\text{cp}} ∣Δn>Ncp,则 FFT 窗将截入相邻符号成分,引发符号间干扰与 ICI,使星座显著扩散并出现误码地板。因此,CTO 的危害呈现“门限性”:CP 内主要体现为相位斜坡,CP 外则导致严重的 ISI/ICI 与性能塌陷。
在这里插入图片描述

5 STO 细定时偏移的模型与危害

细定时偏移(STO)通常指小数采样级的定时误差,即 FFT 窗口起点相对理想位置存在非整数偏移 δ ∈ R \delta\in\mathbb{R} δR,其来源包括定时估计器分辨率不足、采样相位误差以及残余时钟抖动等。与 CTO 的整数移位不同,STO 在带限重构意义下对应一个分数延时算子。若假设 ∣ δ ∣ ≤ N cp |\delta|\le N_{\text{cp}} δNcp且信号在 CP 区间内循环延拓成立,则去 CP 后输入 FFT 的等效序列可表示为
r m [ n ] = x m [ n − δ ] + w m [ n ] , n = 0 , 1 , … , N − 1 , r_m[n]=x_m[n-\delta]+w_m[n],\quad n=0,1,\dots,N-1, rm[n]=xm[nδ]+wm[n],n=0,1,,N1,
其中 x m [ n − δ ] x_m[n-\delta] xm[nδ] 表示对 x m [ n ] x_m[n] xm[n] 的分数延时重采样。理想带限情况下,分数延时在频域对应乘以线性相位
Y m [ k ] ≈ X m [ k ] e − j 2 π N k δ + W m [ k ] . Y_m[k]\approx X_m[k] e^{-j\frac{2\pi}{N}k\delta} + W_m[k]. Ym[k]Xm[k]ejN2πkδ+Wm[k].
由该式可见,在理想条件下 STO 与 CTO 在形式上相似,均引入随子载波索引 k k k 线性变化的相位斜坡,不同之处在于 δ \delta δ 为实数,导致相位更细粒度地变化。然而在实际实现中, x m [ n − δ ] x_m[n-\delta] xm[nδ] 需由插值或分数延时滤波器逼近,若插值器频响不平坦或相位非线性,则会额外引入幅相畸变,进而产生 ICI 残差,使星座图呈现结构性扩散。若 ∣ δ ∣ > N cp |\delta|>N_{\text{cp}} δ>Ncp,则同样会引入 ISI/ICI,且其劣化程度通常比纯整数偏移更难补偿,因为分数误差会与均衡与相位跟踪耦合。系统危害上,STO 直接影响子载波相位一致性与正交性维持,导致信道估计偏差、导频相位不一致以及高阶调制下 EVM 增大,并在存在频偏或相噪时进一步加重性能退化。
在这里插入图片描述

6 IQ 不平衡的模型与危害

IQ 不平衡是直接变频(Direct-Conversion)或低中频接收机中常见的射频前端非理想,其根源在于 I/Q 两路的增益不匹配与正交相位误差,使得复基带表示不再严格满足理想的正交解调条件。设理想复基带信号为 x ( t ) x(t) x(t),IQ 不平衡可用广泛采用的广义线性模型表示为
y ( t ) = α x ( t ) + β x ∗ ( t ) + v ( t ) , y(t)=\alpha x(t)+\beta x^{*}(t)+v(t), y(t)=αx(t)+βx(t)+v(t),
其中 α \alpha α, β ∈ C \beta\in\mathbb{C} βC 为由增益失配与相位失配决定的复系数, ( ⋅ ) ∗ (\cdot)^{*} () 表示共轭, v ( t ) v(t) v(t) 为噪声项。若用增益不平衡系数 g g g 与相位误差 ϕ \phi ϕ 参数化,则一类常用形式为
α = 1 2 ( ( 1 + g ) e − j ϕ / 2 + ( 1 − g ) e j ϕ / 2 ) , β = 1 2 ( ( 1 + g ) e − j ϕ / 2 − ( 1 − g ) e j ϕ / 2 ) . \alpha=\frac{1}{2}\left((1+g)e^{-j\phi/2}+(1-g)e^{j\phi/2}\right),\\\quad \beta=\frac{1}{2}\left((1+g)e^{-j\phi/2}-(1-g)e^{j\phi/2}\right). α=21((1+g)ejϕ/2+(1g)ejϕ/2),β=21((1+g)ejϕ/2(1g)ejϕ/2).
g = 0 g=0 g=0 ϕ = 0 \phi=0 ϕ=0 时有 β = 0 , α = 1 \beta=0,\alpha=1 β=0,α=1,退化为理想接收。对 OFDM 而言,将 y ( t ) y(t) y(t) 采样并做 FFT 后, β x ∗ ( t ) \beta x^{*}(t) βx(t) 项在频域中体现为镜像子载波耦合,即第 k k k 个子载波不仅包含 X m [ k ] X_m[k] Xm[k]的线性变换,还混入了镜像子载波 X m ∗ [ − k ] X_m^{*}[-k] Xm[k] 的泄漏分量,可表示为
Y m [ k ] = α X m [ k ] + β X m ∗ [ − k ] + W m [ k ] , Y_m[k]=\alpha X_m[k]+\beta X_m^{*}[-k]+W_m[k], Ym[k]=αXm[k]+βXm[k]+Wm[k],
其中 − k -k k 表示相对于直流的镜像索引(具体取值依赖子载波编号约定)。该耦合破坏了子载波间的独立性,使得即便在 AWGN 信道下仍出现结构性失真,星座图上常表现为点云形状拉伸、旋转不一致以及由镜像干扰引起的“重影”扩散。系统危害方面,IQ 不平衡会等效降低信干噪比,尤其在高阶 QAM、频谱利用率高且子载波占用对称时更为显著,并会显著增加接收机补偿复杂度,通常需要数字域 IQ 校准或自适应镜像抑制算法予以修正。
在这里插入图片描述

代码开源:hjHe-ee/OFDM-constellation-chart-distortion-CFO-SFO-CTO-STO-and-IQ-imbalance-: Analysis and Mathematical Modeling of Common Constellation Diagram Distortion Mechanisms in OFDM Systems (including CFO carrier frequency offset, SFO sampling frequency offset, CTO coarse timing offset, STO fine timing offset, and discussion of IQ imbalance)

总结

本文围绕 OFDM 系统在 AWGN 信道下的星座图畸变问题,建立了从理想噪声模型到典型同步误差与射频非理想的统一数学描述。在仅 AWGN 条件下,FFT 输出满足 Y m [ k ] = X m [ k ] + W m [ k ] Y_m[k]=X_m[k]+W_m[k] Ym[k]=Xm[k]+Wm[k],星座散布呈圆对称高斯特征且由 SNR 决定;CFO 通过引入随时间增长的相位旋转与 ICI 系数破坏子载波正交性,使星座表现为公共相位旋转叠加扩散;SFO 由于采样率失配导致相位在时间与子载波维度双重累积,并诱发 ICI 残差,进一步强化结构性环状畸变;CTO 在 CP 覆盖范围内主要表现为随子载波线性变化的确定性相位斜坡,而当偏移超出 CP 则引发严重 ISI/ICI;STO 作为分数采样级误差在理想带限条件下亦对应线性相位斜坡,但其实际实现常因插值逼近而引入幅相畸变与额外 ICI;IQ 不平衡则通过镜像子载波耦合 Y m [ k ] = α X m [ k ] + β X m ∗ [ − k ] + W m [ k ] Y_m[k]=\alpha X_m[k]+\beta X_m^{*}[-k]+W_m[k] Ym[k]=αXm[k]+βXm[k]+Wm[k] 破坏频域独立性,导致星座出现结构性拉伸与重影扩散。上述分析表明,OFDM 星座畸变的本质往往不是简单噪声扩散,而是多维同步与前端非理想共同作用下的正交性破坏与相位一致性失效。基于本文的模型推导,可进一步面向实际系统开展 CFO/SFO 联合估计、定时同步精细化以及 IQ 失配自适应补偿等研究,从而在高阶调制与高频谱效率场景下实现更稳健的接收性能。

Logo

DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。

更多推荐