MATLAB 实现帧同步头捕获算法基于 PN 码 帧同步完整仿真链路1、前言上一篇文章我们介绍了最简单的QPSK调制和解调的基本流程现在我们来考虑一个问题对于一个TDMA系统发射端不是时时刻刻都在发送信息而是根据网控的时隙资源分配在对应的时隙发送自己的数据每个突发帧之间存在突发保护那么问题来了我们怎么找到每个突发帧呢数字通信系统接收端三大同步 —— 载波同步、位同步、帧同步帧同步依靠同步头前导码完成帧起始位置捕获是数据正确解调的前提。本文选用经典PN码作为同步头搭建完整 MATLAB 仿真链路完整实现同步头插入、AWGN 信道传输、滑动相关捕获、峰值判决、帧定位功能分析不同 SNR 下同步捕获性能代码全注释、适合通信原理课程学习、工程基带仿真参考。2、PN 码m 序列帧同步基础原理2.1、m 序列 PN 码生成原理PN码定义线性反馈移位寄存器 LFSR 生成最长线性伪随机序列。其特点为自相关函数具有尖锐主峰错位相关值极小适合做同步头。2.2、 PN 同步头滑动相关捕获流程发送端每帧数据起始位置插入一段固定长度 m 序列 PN 码作为前导同步头接收端本地预存相同 m 序列 PN 码模板滑动窗口互相关接收数据流逐比特滑动与本地 PN 码做滑动相乘累加阈值判决相关输出超过预设门限时判定捕获到同步头记录帧起始位置分割有效载荷数据。2.3 、同步性能评价指标捕获峰值PN 码完全对齐时的相关输出最大值漏警概率存在同步头但噪声压制导致峰值低于门限未捕获虚警概率纯噪声 / 业务数据段产生虚假峰值误判帧头捕获时延数据流输入到成功锁定 PN 同步头所需比特数SNR 适应性低信噪比下 PN 码相比短码捕获稳定性更强。3、仿真系统链路图下图为上篇文章中整个QPSK通信的大体流程今天讲的帧同步的位置在如下所示首先看发射端同步头和源发射数据一块组帧构成一个突发帧在进行调制发送。接着看接收端接收机在进行匹配滤波后进行同步头捕获找到有效突发帧在进行定时同步找到最佳采样点进行后续解调。4、matlab核心代码下面代码为PN码产生移位运算产生128位随机数。%% 同步头PN码生成%% 生成127位7阶m序列补1个0得到128点PN% 7阶本原多项式 x^7x1反馈tap[7,1]reg7[1000000];% 初始寄存器全1pn127zeros(1,127);fori1:127fbxor(reg7(7),reg7(1));% 反馈pn127(i)reg7(7);reg7[fb,reg7(1:6)];endpn128[pn127,0];% 补0至128长度滑动相关寻找最大峰代码如下xrx_match;x1x(1:4:end);%4倍采样选择13样点进入相关峰捕获x2x(3:4:end);x1_Idouble(real(x1)0)*2-1;%硬相关只取信号的符号位进行峰值计算x1_Qdouble(imag(x1)0)*2-1;x2_Idouble(real(x2)0)*2-1;x2_Qdouble(imag(x2)0)*2-1;pn128_useflip(pn128)*2-1;%滑动相关捕获本质是对应相乘再相加我这里用卷积代替capture_th1conv(x1_I,pn128_use).^2conv(x1_Q,pn128_use).^2;capture_th2conv(x2_I,pn128_use).^2conv(x2_Q,pn128_use).^2;% [M,I] max(capture_th1);% [M2,I2] max(capture_th2);th10000;%根据符号位数匹配度设定相关峰阈值fork1:length(capture_th1)%滑动查找相关峰if(capture_th1(k)th)yx(4*k1:end);%4倍采样所以位置坐标需要*4break;elseif(capture_th2(k)th)yx(4*k1:end);break;elseendendrx_capturey;5、仿真结果图相关峰捕获图高信噪比下峰值图如下低信噪比峰值图如下仿真结果分析信噪比降低会导致误码率上升同步头的符号匹配数会减少导致相关峰的峰值降低我们需要根据工程需要设定合适的阈值或者基于实时实时噪声功率动态调整判决阈值替代固定门限。全套代码见下面链接https://download.csdn.net/download/wobiaosha2753/93261141
MATLAB 实现帧同步头捕获算法基于 PN 码 帧同步完整仿真链路1、前言上一篇文章我们介绍了最简单的QPSK调制和解调的基本流程现在我们来考虑一个问题对于一个TDMA系统发射端不是时时刻刻都在发送信息而是根据网控的时隙资源分配在对应的时隙发送自己的数据每个突发帧之间存在突发保护那么问题来了我们怎么找到每个突发帧呢数字通信系统接收端三大同步 —— 载波同步、位同步、帧同步帧同步依靠同步头前导码完成帧起始位置捕获是数据正确解调的前提。本文选用经典PN码作为同步头搭建完整 MATLAB 仿真链路完整实现同步头插入、AWGN 信道传输、滑动相关捕获、峰值判决、帧定位功能分析不同 SNR 下同步捕获性能代码全注释、适合通信原理课程学习、工程基带仿真参考。2、PN 码m 序列帧同步基础原理2.1、m 序列 PN 码生成原理PN码定义线性反馈移位寄存器 LFSR 生成最长线性伪随机序列。其特点为自相关函数具有尖锐主峰错位相关值极小适合做同步头。2.2、 PN 同步头滑动相关捕获流程发送端每帧数据起始位置插入一段固定长度 m 序列 PN 码作为前导同步头接收端本地预存相同 m 序列 PN 码模板滑动窗口互相关接收数据流逐比特滑动与本地 PN 码做滑动相乘累加阈值判决相关输出超过预设门限时判定捕获到同步头记录帧起始位置分割有效载荷数据。2.3 、同步性能评价指标捕获峰值PN 码完全对齐时的相关输出最大值漏警概率存在同步头但噪声压制导致峰值低于门限未捕获虚警概率纯噪声 / 业务数据段产生虚假峰值误判帧头捕获时延数据流输入到成功锁定 PN 同步头所需比特数SNR 适应性低信噪比下 PN 码相比短码捕获稳定性更强。3、仿真系统链路图下图为上篇文章中整个QPSK通信的大体流程今天讲的帧同步的位置在如下所示首先看发射端同步头和源发射数据一块组帧构成一个突发帧在进行调制发送。接着看接收端接收机在进行匹配滤波后进行同步头捕获找到有效突发帧在进行定时同步找到最佳采样点进行后续解调。4、matlab核心代码下面代码为PN码产生移位运算产生128位随机数。%% 同步头PN码生成%% 生成127位7阶m序列补1个0得到128点PN% 7阶本原多项式 x^7x1反馈tap[7,1]reg7[1000000];% 初始寄存器全1pn127zeros(1,127);fori1:127fbxor(reg7(7),reg7(1));% 反馈pn127(i)reg7(7);reg7[fb,reg7(1:6)];endpn128[pn127,0];% 补0至128长度滑动相关寻找最大峰代码如下xrx_match;x1x(1:4:end);%4倍采样选择13样点进入相关峰捕获x2x(3:4:end);x1_Idouble(real(x1)0)*2-1;%硬相关只取信号的符号位进行峰值计算x1_Qdouble(imag(x1)0)*2-1;x2_Idouble(real(x2)0)*2-1;x2_Qdouble(imag(x2)0)*2-1;pn128_useflip(pn128)*2-1;%滑动相关捕获本质是对应相乘再相加我这里用卷积代替capture_th1conv(x1_I,pn128_use).^2conv(x1_Q,pn128_use).^2;capture_th2conv(x2_I,pn128_use).^2conv(x2_Q,pn128_use).^2;% [M,I] max(capture_th1);% [M2,I2] max(capture_th2);th10000;%根据符号位数匹配度设定相关峰阈值fork1:length(capture_th1)%滑动查找相关峰if(capture_th1(k)th)yx(4*k1:end);%4倍采样所以位置坐标需要*4break;elseif(capture_th2(k)th)yx(4*k1:end);break;elseendendrx_capturey;5、仿真结果图相关峰捕获图高信噪比下峰值图如下低信噪比峰值图如下仿真结果分析信噪比降低会导致误码率上升同步头的符号匹配数会减少导致相关峰的峰值降低我们需要根据工程需要设定合适的阈值或者基于实时实时噪声功率动态调整判决阈值替代固定门限。全套代码见下面链接https://download.csdn.net/download/wobiaosha2753/93261141
