做通信物理层仿真的人大概都有这个体会教科书里的调制解调链路跑起来很容易MATLAB里几行代码就能画出漂亮的BER曲线。但一旦往信道模型里加入“多普勒频移”四个字整条链路就开始表演花式翻车。我第一次在8-PSK仿真里加入1kHz频偏时星座图直接转成了陀螺误码率从10^-6量级一路崩到接近0.15当时第一反应以为是代码写错了。后来才意识到问题不在代码而在同步算法缺位。这篇文章就围绕这个场景完整走一遍MATLAB下的8-PSK调制解调、频移影响分析、载波同步和符号定时同步的设计思路适合正在做无线通信物理层仿真、或者准备把高阶PSK搬进移动场景的同学参考。1. 为什么8-PSK仿真首先绕不开多普勒频移1.1 多普勒频移在复基带里到底是什么先明确一个基本概念。多普勒频移不是信道里额外加的一段噪声而是收发两端存在相对径向运动时接收信号载波频率发生的偏移。它的物理公式是fd v * cosθ * fc / c其中v是相对运动速度fc是载波频率c是光速。举个具体例子载波2.4GHz终端以120km/h约33m/s接近基站fd大约是267Hz。如果符号速率只有100kSps这个频偏看似只有符号率的0.27%但问题在于它会让星座图持续旋转。在MATLAB里做复基带仿真时实际GHz量级的载波不需要显式建模。我们用复包络表示信号那么多普勒频移在接收端就等效为乘上一个时变的复指数rx(t) tx(t) * exp(j * 2π * fd * t)这个表达式直观上就说明频偏不是给星座图加噪声而是给每个采样点乘一个持续旋转的单位向量。相位会随时间累积这也是为什么频偏的杀伤力远比同等信噪比损失更隐蔽——你很难靠提高发射功率来对抗它。1.2 8-PSK的相位间隔对频偏有多敏感8-PSK每符号携带3比特8个星座点均匀分布在单位圆上相邻星座点的相位间隔是45度。解调时最常用的硬判决准则是把接收符号判到最近的那个星座点上所以判决边界正好落在相邻星座点的中间也就是22.5度。这里有一个非常关键的量化关系每个符号周期内频偏造成的相位旋转量是Δφ 2π * fd * Ts 2π * fd / Rs。当Δφ接近判决边界22.5度时判决几乎必然出错。换算成归一化频偏fd / Rs ≈ 22.5 / 360 6.25%这是理论上的绝对崩溃点。但实际工程中要保证BER退化不明显会让fd / Rs远小于这个值一般要求在一个数量级以上——比如1%以内对应每符号3.6度旋转就已经开始明显劣化0.1%以内才相对安全。QPSK为什么对频偏容忍度高一些因为它的相位间隔是90度判决边界45度同样绝对频偏下比8-PSK从容得多。8-PSK处于一个尴尬位置它的频谱效率比QPSK高1.5倍但相位裕量窄了一半对载波同步的要求高了一个台阶。这就是标题里说的“从多普勒频移条件说起”——不是先研究调制再考虑频移而是一上来就得把频移当作约束条件来设计整个接收机。1.3 仿真里怎样才算“施加了多普勒频移条件”很多初学者的做法是在接收信号上加一个固定的相位偏移比如乘一个exp(j * π/8)。这算“相偏”不是“频偏”。多普勒的核心理特征是相位随时间连续累积所以仿真中必须使用上面的时变复指数形式。在MATLAB里生成一组频移信号非常简单fs Rs * sps; % 采样率 t (0:length(txSignal)-1). / fs; fd 500; % 多普勒频移目标值 rxDoppler txSignal .* exp(1j * 2 * pi * fd * t);这里有一个容易被忽略的细节t必须和采样点数严格对应并且要用列向量否则MATLAB广播机制会给你搞出一个矩阵来。另外多普勒频移是一个标称值真实信道中它可能是时变的比如低轨卫星过顶时fd从正到负连续变化本文先按固定频移处理最后再讨论时变扩展。2. 链路参数设计先定标称值再谈算法2.1 复基带仿真思路与参数表一颗频偏有多大不能只看fd必须联合符号速率一起看。符号速率越高同样fd对性能影响越小。为了让仿真结果有讨论价值我选择了一个“够低但是不算极端”的符号速率确保高频偏条件下BER崩溃现象明显。整个仿真链路参数如下参数取值说明调制方式8-PSK格雷映射每符号3比特符号速率 Rs100 kSps频偏敏感度适中采样率 fs1.6 MHz每符号16个采样点脉冲成形根升余弦RRC滚降系数0.35滤波器阶数10个符号延迟与截断权衡多普勒频移 fd0 Hz ~ 1 kHz对应0~450km/h2.4GHz信噪比Eb/N0 0~18 dB扫描BER仿真比特数30万保证低误码率统计稳定所有处理都在复基带完成。接收端链路顺序是RRC匹配滤波 → 符号采样理想定时或Gardner定时→ 频偏粗估计 → 频偏补偿 → DD-PLL细跟踪 → 解映射 → BER统计。2.2 为什么是根升余弦加0.35滚降很多教程里直接省略脉冲成形用矩形脉冲仿真也能跑通但一旦引入频偏和定时同步没有脉冲成形的系统会产生严重的符号间干扰ISI掩盖同步算法的真实效果。根升余弦滤波器是通信系统里的标准选择它的核心特性是发射端用根升余弦、接收端用同样的根升余弦匹配滤波后等效信道是升余弦滚降满足奈奎斯特第一准则理论上在最佳采样点不存在ISI。滚降系数0.35是工程中常见的折中——滚降越小频谱利用率越高但对定时误差越敏感滚降越大越鲁棒但占用带宽多。MATLAB里生成RRC滤波器用一行命令rrcFilter rcosdesign(0.35, 10, sps, sqrt);第二个参数10表示滤波器跨越10个符号周期sps是每个符号的采样点数。滤波器长度是10*sps1个采样点。这里的教训是滤波器有群延迟接收端必须补偿掉这个延迟再做后续处理很多新手仿真BER莫名其妙变差一大半是这个延迟对齐没处理干净。3. 8-PSK调制解调的MATLAB实现与星座验证3.1 格雷映射的实现细节8-PSK的映射有两种常用方式自然映射和格雷映射。自然映射下相邻星座点对应的3比特可能差好几位格雷映射下相邻相位点只差1个比特。由于噪声导致的最大概率错误是判到相邻星座点格雷映射可以把一次符号错误从“可能错3比特”降为“只错1比特”对BER是直接利好。MATLAB内置的pskmod本身可以指定gray选项但如果想理解映射关系可以手动构建constellation exp(1j * (0:7) * pi/4 1j*pi/8).; % 格雷映射表相邻星座点只差1比特 grayIdx [0 1 3 2 6 7 5 4]; mappingTable constellation(grayIdx 1);注意我加了pi/8的相位偏置让星座点整体旋转22.5度这样星座点不会落在I/Q轴正负方向上避免某些信道或同步算法下出现退化。如果直接使用pskmod需要这样调用symbols pskmod(bitIdx, 8, pi/8, gray);bitIdx是0~7的整数索引需要先把比特流每3位转成十进制。3.2 脉冲成形与匹配滤波发射端在符号序列后面跟上脉冲成形本质上是一个上采样加滤波的过程。upfirdn可以一步完成txShaped upfirdn(symbols, rrcFilter, sps);接收端匹配滤波就是再过一个同样的RRC滤波器然后再降采样到每符号1个采样点。但直接降采样到1倍符号率会让后续定时同步失去调整空间所以实践中常先降到2倍符号率再通过定时恢复找到最优采样点。这里先按理想定时处理rxShaped upfirdn(rxSignal, rrcFilter, 1, sps); rxSymbols rxShaped(delay1:end); % 先补偿滤波器延迟滤波器延迟对于sps16、跨10个符号的RRC是10*sps个采样点。加上upfirdn引入的边界效应最好去掉前后各一段数据再统计否则BER曲线的低信噪比端会异常。3.3 理想信道下的BER自检在加多普勒之前先验证一条干净链路是否正常。这是整个仿真工程最重要的自检步骤如果理想信道下BER都不对后面加入频移后根本没法判断问题是出在同步算法还是基础链路。8-PSK在高斯信道下的理论BER没有闭式表达式一般用SER公式近似SER ≈ 2 * Q(sqrt(2 * Es/N0) * sin(pi/8))在仿真时直接用蒙特卡洛统计BER在Eb/N0 12dB左右8-PSK的理论BER大约在10^-5量级。如果仿真得出的BER明显高于这个量级优先检查滤波器群延迟补偿是否准确符号采样点是否正好在最佳眼图张开处功率归一化是否正确发射符号能量是否为1常见的功率归一化陷阱是忘了对RRC滤波器做归一化导致等效SNR偏差好几个dB。最简单的方法是用rcosdesign生成的滤波器本来就是单位能量但如果手动设计了其他滤波器记得做rrcFilter rrcFilter / sqrt(sum(rrcFilter.^2))。4. 频偏崩溃实验不加同步先量化直观感受4.1 实验设置与代码在做同步算法之前先做一组“裸奔”实验理想定时、匹配滤波正常但完全不做频偏补偿只在接收端直接解调统计不同fd下的BER。这样能直观量化频偏的破坏力。fdList [0, 50, 200, 500, 1000]; for i 1:length(fdList) fd fdList(i); rx txShaped .* exp(1j * 2 * pi * fd * t); rx awgn(rx, snr, measured); rxFilt upfirdn(rx, rrcFilter, 1, sps); rxFilt rxFilt(delay1 : delay length(symbols)); rxDec rxFilt(1:sps:end); % 理想定时抽样 demodBits pskdemod(rxDec, 8, pi/8, gray); [~, ber(i)] biterr(bits, demodBits); end这段代码里rxFilt(delay1 : delay length(symbols))很关键它保证时序上每个抽样点正好对应原始符号位置。4.2 结果BER平台与星座旋转我跑出来的典型结果Eb/N012dB约3万符号大致如下fdHzfd/RsBER00约2e-5500.0005约3e-32000.002约2e-25000.005约6e-210000.01约1.3e-1即使频偏只有符号速率的0.05%BER已经从10^-5量级恶化到10^-3量级性能损失超过100倍。继续增大fdBER趋近于8-PSK在随机相位判决下的错误概率也就是接近0.5的水平。画星座图会看得更直观。没有频偏时8个相位点清晰分明加入500Hz频偏后星座点变成了一圈模糊的圆环每个符号采样时刻的相位都在旋转判决器基本是在“猜”。这说明了一个重要事实频偏导致的BER恶化并不会随着SNR提高而消失它会把BER曲线压成一个平台这就是常说的“错误地板”error floor。4.3 为什么会崩相位累积与判决边界要理解这个“错误地板”需要回到判决机制。每符号周期内频偏让相位旋转Δφ 2π * fd / Rs。当Δφ远小于22.5度时瞬时看起来每个符号都在正确区域附近但如果没有同步相位会随着符号序号不断累积。累积到某个时刻星座点整体旋转超过判决边界即使没有噪声也会开始连续出错。这个累积效应是积分性的和SNR无关所以SNR再高也无法消除。这就是为什么必须做载波同步——不是简单地把信号“转正”一次而是要让接收机的本地振荡器连续跟踪发射端的相位变化。频偏估计与补偿的本质就是估计出这个不断累加的相位斜率并在解调前把它乘回去。5. 载波同步粗频偏估计加DD-PLL双级结构5.1 M次方谱频偏估计的原理与实现载波同步最常见的做法是两级结构先做一个粗频率估计把大部分频偏补偿掉再用锁相环精确跟踪残余频偏和相位抖动。如果只做PLL不做粗估计环路带宽需要设得很宽才能捕获大频偏但宽带宽又引入大量噪声性能会打折扣。粗估计这里介绍M次方谱法原理非常巧妙。8-PSK星座有8个对称相位点把接收信号做8次方运算后数据相位被乘以8所有星座点都会旋转到同一个相位上——数据调制被“抹平”了。此时信号变成了一个频率为8*fd的单音剩余能量集中到一个谱线上。用FFT找到这个峰值的位置再除以8就得到频偏估计。MATLAB实现如下% rxSymbols: 匹配滤波后符号级采样每符号1点 xM rxSymbols .^ 8; Nfft 2^nextpow2(length(xM)); X fftshift(fft(xM, Nfft)); fAxis (-Nfft/2 : Nfft/2-1) / Nfft * Rs; % 注意此时采样率是符号速率Rs [~, peakIdx] max(abs(X)); fPeak fAxis(peakIdx); fdEst fPeak / 8;这里有几个关键细节值得圈出来。第一FFT的输入是符号率采样信号不是过采样信号。如果直接用过采样信号做FFT频谱范围变成[-fs/2, fs/2]峰值定位精度反而会因为点数分配更粗糙而下降。第二频率分辨率是Rs / Nfft。如果符号数太少比如只有1000个符号Nfft1024分辨率约为98Hz粗估计误差会非常大PLL需要花更长时间收敛。所以实际项目中符号数至少上万最好有5万个以上粗估计精度才能到Hz量级。第三8次方运算会把噪声也抬起来低SNR下峰值可能被噪声淹没。如果SNR低于5dB建议先做多个符号块的平均频谱再取峰值。5.2 二阶数字PLL的环路设计粗估计补偿后残留频偏通常在几十Hz量级此时用判决辅助锁相环Decision-Directed PLLDD-PLL做细跟踪。DD-PLL的核心思想是先用当前相位校正信号做硬判决得到参考符号再用接收符号与参考符号的相位差作为误差信号驱动环路滤波器校正NCO相位。一阶PLL只能消除固定相偏二阶PLL才能跟踪固定频偏。这里用标准二阶环路结构% 环路滤波器参数BL为噪声带宽zeta为阻尼系数 BL 500; % 噪声带宽约500Hz约为符号率的0.5% zeta 0.707; % 常规阻尼 theta BL / (Rs * (zeta 1/(4*zeta))); Kp 2 * zeta * theta; % 比例系数 Ki theta^2; % 积分系数 % DD-PLL 主循环 phaseAcc 0; freqAcc 0; pllOut zeros(size(rxComp)); constell exp(1j * (0:7) * pi/4 1j*pi/8).; % 接收判决星座 for n 1:length(rxComp) % 相位校正 corrected rxComp(n) * exp(-1j * phaseAcc); pllOut(n) corrected; % 硬判决 [~, idx] min(abs(corrected - constell)); decision constell(idx); % 鉴相器用相位误差而不是虚部误差 err angle(corrected * conj(decision)); % 环路滤波 freqAcc freqAcc Ki * err; phaseAcc phaseAcc freqAcc Kp * err; end这个环路有几个值得注意的参数选择。BL要折中BL太窄比如10HzPLL收敛慢且对于时变多普勒来不及跟踪BL太宽比如5000Hz噪声抑制差星座图散得厉害。对100kSps符号率500Hz是一个不错的起点。如果fd固定不变收敛后可以降低BL进一步压噪声如果fd时变需要保留一定带宽跟踪变化。鉴相器用angle()计算完整相位误差比直接用imag(r * conj(d))线性近似在相位误差大时更准确。代价是计算量略高但在仿真里无所谓。环路滤波器的参数推导基于连续时间二阶环的离散化这里不展开公式推导但可以记住结论噪声带宽BL、阻尼系数ζ和环路增益Kp、Ki的换算关系就是上面那三行。如果BER结果不理想优先检查BL是否过大其次是确认初始phaseAcc是否设成粗估计的残余相位。5.3 相位模糊仿真中最容易漏掉的一环这是整个同步链路里最坑的细节。在M次方谱估计和DD-PLL中8-PSK映射会有8重相位模糊。也就是说PLL锁定后系统可能收敛到任意一个45度倍数被“转到正确位置”的相位上但和真正的星座方向相比可能整整偏了一个或多个45度。原因在于8-PSK星座旋转45度后和原星座完全重合P
多普勒频移下的8-PSK MATLAB仿真与载波同步设计
做通信物理层仿真的人大概都有这个体会教科书里的调制解调链路跑起来很容易MATLAB里几行代码就能画出漂亮的BER曲线。但一旦往信道模型里加入“多普勒频移”四个字整条链路就开始表演花式翻车。我第一次在8-PSK仿真里加入1kHz频偏时星座图直接转成了陀螺误码率从10^-6量级一路崩到接近0.15当时第一反应以为是代码写错了。后来才意识到问题不在代码而在同步算法缺位。这篇文章就围绕这个场景完整走一遍MATLAB下的8-PSK调制解调、频移影响分析、载波同步和符号定时同步的设计思路适合正在做无线通信物理层仿真、或者准备把高阶PSK搬进移动场景的同学参考。1. 为什么8-PSK仿真首先绕不开多普勒频移1.1 多普勒频移在复基带里到底是什么先明确一个基本概念。多普勒频移不是信道里额外加的一段噪声而是收发两端存在相对径向运动时接收信号载波频率发生的偏移。它的物理公式是fd v * cosθ * fc / c其中v是相对运动速度fc是载波频率c是光速。举个具体例子载波2.4GHz终端以120km/h约33m/s接近基站fd大约是267Hz。如果符号速率只有100kSps这个频偏看似只有符号率的0.27%但问题在于它会让星座图持续旋转。在MATLAB里做复基带仿真时实际GHz量级的载波不需要显式建模。我们用复包络表示信号那么多普勒频移在接收端就等效为乘上一个时变的复指数rx(t) tx(t) * exp(j * 2π * fd * t)这个表达式直观上就说明频偏不是给星座图加噪声而是给每个采样点乘一个持续旋转的单位向量。相位会随时间累积这也是为什么频偏的杀伤力远比同等信噪比损失更隐蔽——你很难靠提高发射功率来对抗它。1.2 8-PSK的相位间隔对频偏有多敏感8-PSK每符号携带3比特8个星座点均匀分布在单位圆上相邻星座点的相位间隔是45度。解调时最常用的硬判决准则是把接收符号判到最近的那个星座点上所以判决边界正好落在相邻星座点的中间也就是22.5度。这里有一个非常关键的量化关系每个符号周期内频偏造成的相位旋转量是Δφ 2π * fd * Ts 2π * fd / Rs。当Δφ接近判决边界22.5度时判决几乎必然出错。换算成归一化频偏fd / Rs ≈ 22.5 / 360 6.25%这是理论上的绝对崩溃点。但实际工程中要保证BER退化不明显会让fd / Rs远小于这个值一般要求在一个数量级以上——比如1%以内对应每符号3.6度旋转就已经开始明显劣化0.1%以内才相对安全。QPSK为什么对频偏容忍度高一些因为它的相位间隔是90度判决边界45度同样绝对频偏下比8-PSK从容得多。8-PSK处于一个尴尬位置它的频谱效率比QPSK高1.5倍但相位裕量窄了一半对载波同步的要求高了一个台阶。这就是标题里说的“从多普勒频移条件说起”——不是先研究调制再考虑频移而是一上来就得把频移当作约束条件来设计整个接收机。1.3 仿真里怎样才算“施加了多普勒频移条件”很多初学者的做法是在接收信号上加一个固定的相位偏移比如乘一个exp(j * π/8)。这算“相偏”不是“频偏”。多普勒的核心理特征是相位随时间连续累积所以仿真中必须使用上面的时变复指数形式。在MATLAB里生成一组频移信号非常简单fs Rs * sps; % 采样率 t (0:length(txSignal)-1). / fs; fd 500; % 多普勒频移目标值 rxDoppler txSignal .* exp(1j * 2 * pi * fd * t);这里有一个容易被忽略的细节t必须和采样点数严格对应并且要用列向量否则MATLAB广播机制会给你搞出一个矩阵来。另外多普勒频移是一个标称值真实信道中它可能是时变的比如低轨卫星过顶时fd从正到负连续变化本文先按固定频移处理最后再讨论时变扩展。2. 链路参数设计先定标称值再谈算法2.1 复基带仿真思路与参数表一颗频偏有多大不能只看fd必须联合符号速率一起看。符号速率越高同样fd对性能影响越小。为了让仿真结果有讨论价值我选择了一个“够低但是不算极端”的符号速率确保高频偏条件下BER崩溃现象明显。整个仿真链路参数如下参数取值说明调制方式8-PSK格雷映射每符号3比特符号速率 Rs100 kSps频偏敏感度适中采样率 fs1.6 MHz每符号16个采样点脉冲成形根升余弦RRC滚降系数0.35滤波器阶数10个符号延迟与截断权衡多普勒频移 fd0 Hz ~ 1 kHz对应0~450km/h2.4GHz信噪比Eb/N0 0~18 dB扫描BER仿真比特数30万保证低误码率统计稳定所有处理都在复基带完成。接收端链路顺序是RRC匹配滤波 → 符号采样理想定时或Gardner定时→ 频偏粗估计 → 频偏补偿 → DD-PLL细跟踪 → 解映射 → BER统计。2.2 为什么是根升余弦加0.35滚降很多教程里直接省略脉冲成形用矩形脉冲仿真也能跑通但一旦引入频偏和定时同步没有脉冲成形的系统会产生严重的符号间干扰ISI掩盖同步算法的真实效果。根升余弦滤波器是通信系统里的标准选择它的核心特性是发射端用根升余弦、接收端用同样的根升余弦匹配滤波后等效信道是升余弦滚降满足奈奎斯特第一准则理论上在最佳采样点不存在ISI。滚降系数0.35是工程中常见的折中——滚降越小频谱利用率越高但对定时误差越敏感滚降越大越鲁棒但占用带宽多。MATLAB里生成RRC滤波器用一行命令rrcFilter rcosdesign(0.35, 10, sps, sqrt);第二个参数10表示滤波器跨越10个符号周期sps是每个符号的采样点数。滤波器长度是10*sps1个采样点。这里的教训是滤波器有群延迟接收端必须补偿掉这个延迟再做后续处理很多新手仿真BER莫名其妙变差一大半是这个延迟对齐没处理干净。3. 8-PSK调制解调的MATLAB实现与星座验证3.1 格雷映射的实现细节8-PSK的映射有两种常用方式自然映射和格雷映射。自然映射下相邻星座点对应的3比特可能差好几位格雷映射下相邻相位点只差1个比特。由于噪声导致的最大概率错误是判到相邻星座点格雷映射可以把一次符号错误从“可能错3比特”降为“只错1比特”对BER是直接利好。MATLAB内置的pskmod本身可以指定gray选项但如果想理解映射关系可以手动构建constellation exp(1j * (0:7) * pi/4 1j*pi/8).; % 格雷映射表相邻星座点只差1比特 grayIdx [0 1 3 2 6 7 5 4]; mappingTable constellation(grayIdx 1);注意我加了pi/8的相位偏置让星座点整体旋转22.5度这样星座点不会落在I/Q轴正负方向上避免某些信道或同步算法下出现退化。如果直接使用pskmod需要这样调用symbols pskmod(bitIdx, 8, pi/8, gray);bitIdx是0~7的整数索引需要先把比特流每3位转成十进制。3.2 脉冲成形与匹配滤波发射端在符号序列后面跟上脉冲成形本质上是一个上采样加滤波的过程。upfirdn可以一步完成txShaped upfirdn(symbols, rrcFilter, sps);接收端匹配滤波就是再过一个同样的RRC滤波器然后再降采样到每符号1个采样点。但直接降采样到1倍符号率会让后续定时同步失去调整空间所以实践中常先降到2倍符号率再通过定时恢复找到最优采样点。这里先按理想定时处理rxShaped upfirdn(rxSignal, rrcFilter, 1, sps); rxSymbols rxShaped(delay1:end); % 先补偿滤波器延迟滤波器延迟对于sps16、跨10个符号的RRC是10*sps个采样点。加上upfirdn引入的边界效应最好去掉前后各一段数据再统计否则BER曲线的低信噪比端会异常。3.3 理想信道下的BER自检在加多普勒之前先验证一条干净链路是否正常。这是整个仿真工程最重要的自检步骤如果理想信道下BER都不对后面加入频移后根本没法判断问题是出在同步算法还是基础链路。8-PSK在高斯信道下的理论BER没有闭式表达式一般用SER公式近似SER ≈ 2 * Q(sqrt(2 * Es/N0) * sin(pi/8))在仿真时直接用蒙特卡洛统计BER在Eb/N0 12dB左右8-PSK的理论BER大约在10^-5量级。如果仿真得出的BER明显高于这个量级优先检查滤波器群延迟补偿是否准确符号采样点是否正好在最佳眼图张开处功率归一化是否正确发射符号能量是否为1常见的功率归一化陷阱是忘了对RRC滤波器做归一化导致等效SNR偏差好几个dB。最简单的方法是用rcosdesign生成的滤波器本来就是单位能量但如果手动设计了其他滤波器记得做rrcFilter rrcFilter / sqrt(sum(rrcFilter.^2))。4. 频偏崩溃实验不加同步先量化直观感受4.1 实验设置与代码在做同步算法之前先做一组“裸奔”实验理想定时、匹配滤波正常但完全不做频偏补偿只在接收端直接解调统计不同fd下的BER。这样能直观量化频偏的破坏力。fdList [0, 50, 200, 500, 1000]; for i 1:length(fdList) fd fdList(i); rx txShaped .* exp(1j * 2 * pi * fd * t); rx awgn(rx, snr, measured); rxFilt upfirdn(rx, rrcFilter, 1, sps); rxFilt rxFilt(delay1 : delay length(symbols)); rxDec rxFilt(1:sps:end); % 理想定时抽样 demodBits pskdemod(rxDec, 8, pi/8, gray); [~, ber(i)] biterr(bits, demodBits); end这段代码里rxFilt(delay1 : delay length(symbols))很关键它保证时序上每个抽样点正好对应原始符号位置。4.2 结果BER平台与星座旋转我跑出来的典型结果Eb/N012dB约3万符号大致如下fdHzfd/RsBER00约2e-5500.0005约3e-32000.002约2e-25000.005约6e-210000.01约1.3e-1即使频偏只有符号速率的0.05%BER已经从10^-5量级恶化到10^-3量级性能损失超过100倍。继续增大fdBER趋近于8-PSK在随机相位判决下的错误概率也就是接近0.5的水平。画星座图会看得更直观。没有频偏时8个相位点清晰分明加入500Hz频偏后星座点变成了一圈模糊的圆环每个符号采样时刻的相位都在旋转判决器基本是在“猜”。这说明了一个重要事实频偏导致的BER恶化并不会随着SNR提高而消失它会把BER曲线压成一个平台这就是常说的“错误地板”error floor。4.3 为什么会崩相位累积与判决边界要理解这个“错误地板”需要回到判决机制。每符号周期内频偏让相位旋转Δφ 2π * fd / Rs。当Δφ远小于22.5度时瞬时看起来每个符号都在正确区域附近但如果没有同步相位会随着符号序号不断累积。累积到某个时刻星座点整体旋转超过判决边界即使没有噪声也会开始连续出错。这个累积效应是积分性的和SNR无关所以SNR再高也无法消除。这就是为什么必须做载波同步——不是简单地把信号“转正”一次而是要让接收机的本地振荡器连续跟踪发射端的相位变化。频偏估计与补偿的本质就是估计出这个不断累加的相位斜率并在解调前把它乘回去。5. 载波同步粗频偏估计加DD-PLL双级结构5.1 M次方谱频偏估计的原理与实现载波同步最常见的做法是两级结构先做一个粗频率估计把大部分频偏补偿掉再用锁相环精确跟踪残余频偏和相位抖动。如果只做PLL不做粗估计环路带宽需要设得很宽才能捕获大频偏但宽带宽又引入大量噪声性能会打折扣。粗估计这里介绍M次方谱法原理非常巧妙。8-PSK星座有8个对称相位点把接收信号做8次方运算后数据相位被乘以8所有星座点都会旋转到同一个相位上——数据调制被“抹平”了。此时信号变成了一个频率为8*fd的单音剩余能量集中到一个谱线上。用FFT找到这个峰值的位置再除以8就得到频偏估计。MATLAB实现如下% rxSymbols: 匹配滤波后符号级采样每符号1点 xM rxSymbols .^ 8; Nfft 2^nextpow2(length(xM)); X fftshift(fft(xM, Nfft)); fAxis (-Nfft/2 : Nfft/2-1) / Nfft * Rs; % 注意此时采样率是符号速率Rs [~, peakIdx] max(abs(X)); fPeak fAxis(peakIdx); fdEst fPeak / 8;这里有几个关键细节值得圈出来。第一FFT的输入是符号率采样信号不是过采样信号。如果直接用过采样信号做FFT频谱范围变成[-fs/2, fs/2]峰值定位精度反而会因为点数分配更粗糙而下降。第二频率分辨率是Rs / Nfft。如果符号数太少比如只有1000个符号Nfft1024分辨率约为98Hz粗估计误差会非常大PLL需要花更长时间收敛。所以实际项目中符号数至少上万最好有5万个以上粗估计精度才能到Hz量级。第三8次方运算会把噪声也抬起来低SNR下峰值可能被噪声淹没。如果SNR低于5dB建议先做多个符号块的平均频谱再取峰值。5.2 二阶数字PLL的环路设计粗估计补偿后残留频偏通常在几十Hz量级此时用判决辅助锁相环Decision-Directed PLLDD-PLL做细跟踪。DD-PLL的核心思想是先用当前相位校正信号做硬判决得到参考符号再用接收符号与参考符号的相位差作为误差信号驱动环路滤波器校正NCO相位。一阶PLL只能消除固定相偏二阶PLL才能跟踪固定频偏。这里用标准二阶环路结构% 环路滤波器参数BL为噪声带宽zeta为阻尼系数 BL 500; % 噪声带宽约500Hz约为符号率的0.5% zeta 0.707; % 常规阻尼 theta BL / (Rs * (zeta 1/(4*zeta))); Kp 2 * zeta * theta; % 比例系数 Ki theta^2; % 积分系数 % DD-PLL 主循环 phaseAcc 0; freqAcc 0; pllOut zeros(size(rxComp)); constell exp(1j * (0:7) * pi/4 1j*pi/8).; % 接收判决星座 for n 1:length(rxComp) % 相位校正 corrected rxComp(n) * exp(-1j * phaseAcc); pllOut(n) corrected; % 硬判决 [~, idx] min(abs(corrected - constell)); decision constell(idx); % 鉴相器用相位误差而不是虚部误差 err angle(corrected * conj(decision)); % 环路滤波 freqAcc freqAcc Ki * err; phaseAcc phaseAcc freqAcc Kp * err; end这个环路有几个值得注意的参数选择。BL要折中BL太窄比如10HzPLL收敛慢且对于时变多普勒来不及跟踪BL太宽比如5000Hz噪声抑制差星座图散得厉害。对100kSps符号率500Hz是一个不错的起点。如果fd固定不变收敛后可以降低BL进一步压噪声如果fd时变需要保留一定带宽跟踪变化。鉴相器用angle()计算完整相位误差比直接用imag(r * conj(d))线性近似在相位误差大时更准确。代价是计算量略高但在仿真里无所谓。环路滤波器的参数推导基于连续时间二阶环的离散化这里不展开公式推导但可以记住结论噪声带宽BL、阻尼系数ζ和环路增益Kp、Ki的换算关系就是上面那三行。如果BER结果不理想优先检查BL是否过大其次是确认初始phaseAcc是否设成粗估计的残余相位。5.3 相位模糊仿真中最容易漏掉的一环这是整个同步链路里最坑的细节。在M次方谱估计和DD-PLL中8-PSK映射会有8重相位模糊。也就是说PLL锁定后系统可能收敛到任意一个45度倍数被“转到正确位置”的相位上但和真正的星座方向相比可能整整偏了一个或多个45度。原因在于8-PSK星座旋转45度后和原星座完全重合P
