北斗B1I扩频码Matlab生成原理与完整实现

简介本资源是一份面向卫星导航信号处理初学者与MATLAB仿真实践者的北斗B1I路扩频码生成与基础接收链路仿真代码包聚焦北斗系统核心调制技术解决扩频码构造、BPSK调制、AWGN信道建模及粗同步捕获等关键环节的动手实现问题。压缩包共3个MATLAB源文件.m总大小仅7KB轻量紧凑涵盖BD2B1_AcquisitionB1I信号粗捕获模块、BD2B1_SoftReceiver软接收机主流程及配套GPS_Tracking参考函数便于理解北斗与GPS扩频体制差异并开展对比实验。已有575人学习下载代码严格依据北斗公开接口规范设计完整实现Gold长码生成、短码叠加、扩频调制、噪声注入与相关峰检测全流程可直接运行观察扩频序列自相关特性、捕获门限响应及不同SNR下的同步性能是开展北斗信号处理课程设计、毕业设计或算法验证的理想起点。 北斗B1I的扩频码是做接收机基带仿真绕不开的一个基础模块。很多刚接触GNSS方向的同学一开始都会去翻文档、找现成代码结果发现网上的实现五花八门生成的码要么跟资料对不上要么自己都没法验证对错。这篇文章用matlab把B1I扩频码的生成原理、完整代码和相关特性验证一次讲透提供一套可以直接跑的完整实现也会分享一些我实际调试时踩过的坑。适合正在做北斗信号仿真、课程设计或者接收机捕获跟踪算法验证的读者参考。1. 北斗B1I信号与扩频码基础1.1 B1I信号的关键参数北斗B1I信号是北斗二号、北斗三号系统公开服务信号中的一种工作在L频段。很多人第一次接触这个信号时容易被一堆参数搞得一头雾水。这里先列出一张核心参数表把后续生成扩频码需要关注的信息都放进去。参数名称参数值说明载波频率1561.098 MHzB1I信号的中心频率码速率2.046 Mcps扩频码的码片速率码长2046 chips一个完整扩频码周期的码片数调制方式BPSK(2)二进制相移键控测距码类型截短Gold码由两个m序列模2加生成主瓣带宽约4.092 MHz主瓣零点间带宽从这张表可以看到B1I信号的码速率是2.046 Mcps码长是2046个码片所以一个完整码周期的时间正好是1毫秒。这个1毫秒的周期跟北斗的导航电文速率、比特对齐方式都有直接关系接收机在做捕获的时候通常就是按1毫秒的相干积分长度来处理扩频码的。1.2 为什么导航信号选Gold码GPS的C/A码用的是Gold码北斗B1I也用了Gold码这背后是有明确工程考虑的。扩频码在导航信号里承担两个核心任务一是通过码分多址区分不同卫星二是作为测距的标尺。要测距码的自相关性要好主峰要尖锐要区分卫星不同卫星的码之间互相关性要尽量低。如果用单一m序列虽然自相关性很好但m序列的互相关特性不稳定不同相位偏移下会出现比较大的互相关峰值用户量一多就容易误捕。Gold码是由两个优选m序列模2加生成的码族它的巧妙之处在于码族内任意两个码的互相关值都被限制在一个已知的范围内不会出现离谱的尖峰。这就像给每颗卫星发了一个身份证号接收机靠这个号码来识别是哪颗卫星在发信号同时号码之间足够“不相似”不会认错人。北斗B1I选择Gold码还有一个细节就是码长取2046而不是标准的2047。标准的m序列周期是2的11次方减1等于2047B1I在生成后截短了一位变成2046码片。这样做一方面是为了码周期刚好等于1毫秒与整毫秒的帧结构对齐另一方面也说明B1I的码不是教科书上标准的完整周期Gold码在相关特性、频谱特性上会有一些细微变化这个在后面验证环节会提到。1.3 B1I扩频码的组成结构B1I的扩频码结构可以用一句话概括两个11级线性反馈移位寄存器分别生成m序列G1和G2再把G1与经过相位调整的G2逐位模2加得到某颗卫星的测距码。更具体地说G1和G2的生成多项式在北斗ICD文件里有明确规定。G1的生成多项式是1 x x^3 x^4 x^8 x^9 x^11G2的生成多项式是1 x x^2 x^3 x^4 x^5 x^8 x^9 x^11两个寄存器都是11级初始状态为全1。G1序列直接输出G2序列则根据不同卫星的PRN号做不同的相位偏移或抽头组合再和G1模2加就得到了该卫星独有的扩频码。这个结构跟GPS C/A码的生成方式非常像只是寄存器级数、生成多项式和码长不同。理解了这个框架再看matlab代码就会清晰很多无非是把这两个多项式翻译成代码再把PRN到G2相位的映射表整理好。2. matlab生成B1I扩频码的设计思路2.1 线性反馈移位寄存器实现原理在matlab里生成m序列最直接的方式就是模拟线性反馈移位寄存器LFSR的工作过程。一个11级的LFSR本质上是一个11位的状态存储器每个时钟周期做两件事输出当前状态中的某一位同时根据反馈抽头计算一个新的输入位放到寄存器的入口。用代码来表示就是维护一个长度为n的向量每次循环取向量最后一位作为输出然后根据参与反馈的位置做模2加也就是异或再把寄存器整体右移一位把异或结果放到最前面。这个“取模2加”的过程对应了生成多项式中各位的系数。这里有个容易混淆的点多项式里的x^i对应的是寄存器中的第几个位置。在我的实现里规定寄存器向量reg的第1位到第n位对应多项式中的x^1到x^n项输出取reg(n)。而生成多项式中除了常数项1之外的非零项就是反馈抽头的位置。比如G1多项式从x^1到x^11的非零项在1、3、4、8、9、11这些位置那么反馈计算就是把这些位置的寄存器状态取出来做异或。2.2 G1、G2生成多项式的matlab映射把上面的原理翻译成matlab代码就是下面这个核心函数。这个函数可以复用于G1和G2只需要传入不同的抽头向量。function seq lfsr_generate(taps, n, len) % LFSR生成m序列 % taps: 反馈抽头位置向量如[1 3 4 8 9 11] % n: 寄存器级数 % len: 输出序列长度 reg ones(1, n); seq zeros(1, len); for i 1:len seq(i) reg(n); fb mod(sum(reg(taps)), 2); reg [fb, reg(1:n-1)]; end end对于G1抽头向量是[1 3 4 8 9 11]对于G2抽头向量是[1 2 3 4 5 8 9 11]。调用方式如下n 11; len 2046; G1_taps [1 3 4 8 9 11]; G2_taps [1 2 3 4 5 8 9 11]; G1 lfsr_generate(G1_taps, n, len); G2 lfsr_generate(G2_taps, n, len);这一段代码跑完之后G1和G2就是两条长度为2046的m序列。注意初始状态全部为1输出取的是寄存器最高位这两点决定了序列的相位和极性。如果你的实现里初始状态或者输出位定义不一致生成的序列可能会跟别人的实现差一个反向或者差一个循环移位后面会专门讲这个问题。2.3 卫星号到G2相位的对应关系有了G1和G2的基础序列接下来就是最关键的一步如何根据PRN号得到该卫星的扩频码。北斗ICD中每颗卫星的PRN号都对应一个G2序列的相位配置。这个配置在工程实现上有两种常见形式。第一种是延迟移位方式。就是ICD直接给出该PRN对应的G2延迟码片数比如PRN 1对应延迟d个码片那么在软件里把G2序列做循环移位再跟G1异或即可。这种方式非常直观适合在matlab里实现。第二种是抽头组合方式。这种方法跟GPS C/A码的经典做法一致在G2寄存器运行过程中同时取两个抽头位置的状态做异或将这个异或结果作为与G1相加的G2相位序列。数学上这个序列等价于G2序列的某个延迟版本但实现上不需要先缓存完整的G2再移位而是可以边跑边输出更贴近硬件实现。下面代码用延迟移位方式实现PRN到G2延迟的映射表需要参考你手头的ICD版本。我在代码里放了一个示例配置正式开发时请务必核对ICD原文。% PRN号到G2延迟量映射表示例数据以实际ICD为准 g2_delay_table [ 1; % PRN 1 2; % PRN 2 3; % PRN 3 4; % PRN 4 5; % PRN 5 % ... 后续按ICD补齐 ]; function code beidou_b1i_code(prn) n 11; len 2046; G1_taps [1 3 4 8 9 11]; G2_taps [1 2 3 4 5 8 9 11]; G1 lfsr_generate(G1_taps, n, len); G2 lfsr_generate(G2_taps, n, len); delay g2_delay_table(prn); G2_shifted circshift(G2, [0, delay]); code mod(G1 G2_shifted, 2); end这个实现的优点是代码短、逻辑清楚替换查表数据即可。缺点是需要一次性生成完整的G2序列再移位对于2046码片的长度来说完全没有性能压力。如果你希望用更接近硬件的方式实现抽头组合可以在LFSR循环中同时取G2寄存器的两个位置做异或效果是等效的。3. 完整matlab实现与相关特性验证3.1 一套可直接运行的生成代码为了方便读者直接跑通流程我把前面的代码整合成一个完整的脚本。这个脚本会生成指定PRN的B1I扩频码并画出它的时域波形。%% 北斗B1I扩频码生成脚本 clear; clc; close all; n 11; len 2046; G1_taps [1 3 4 8 9 11]; G2_taps [1 2 3 4 5 p a hrefhttps://download.csdn.net/download/wouderw/87444338 stylecolor:#ec7500;font-size:14px; 本文还有配套的精品资源点击获取 /a img altmenu-r.4af5f7ec.gif srchttps://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif stylewidth:16px;margin-left:4px;vertical-align:text-bottom;cursor:text; /p