GNSS四大类测距码生成方式（对标官方ICD、仿真实操要点）

先总述底层区分：

1. LFSR Gold码（双移位寄存器m序列异或）：GPS C/A、BDS B1I/B2a/B3I，传统主流，硬件/FPGA极易实现；
2. 单LFSR m序列：GPS L2C CM/CL，单寄存器直接输出，无第二路移位；
3. Weil数论码（无任何LFSR）：BDS B1C、GPS L1C，纯有限域模运算，新一代高精度信号；
4. 预存储查表码：Galileo E1 OS，离线生成全量码存入ROM，实时仅寻址读取。

前置通用概念：LFSR底层原理（所有移位寄存器类码共用）

LFSR=线性反馈移位寄存器，分斐波那契结构（导航标准ICD统一采用）

1. 寄存器：n个1bit存储单元，编号1（最左输入）~n（最右输出）；
2. 时钟：每1个码时钟脉冲，全体寄存器右移1位；
3. 反馈抽头：从指定寄存器位取出比特做模2异或(XOR，0⊕0=0、0⊕1=1、1⊕1=0)；
4. 反馈值送入第1位作为新输入；
5. 输出：寄存器第n位持续输出1bit码片；
6. m序列特性：n级LFSR最大周期L=2n−1L=2^n-1L=2n−1，自相关尖锐、互相关极低，适合测距。
   

一、第一类：双LFSR Gold码（GPS L1 C/A ）

基础固定参数

码速率：1.023 Mcps；码长1023 chip；周期1ms；10级G1、G2两套LFSR；32颗卫星靠G2抽头偏移区分PRN。

步骤1：全局时钟分频基准

卫星原子钟基准10.23MHz ÷10 =1.023MHz码时钟，同步驱动G1、G2寄存器，移位、输出严格同步。

步骤2：寄存器上电初始化（强制全1）

G1 = [1,1,1,1,1,1,1,1,1,1]（bit1~bit10）
G2 = [1,1,1,1,1,1,1,1,1,1]
⚠️易错：不能初始全0，全0会锁死输出恒0，无法生成序列。

步骤3：G1 m序列逐时钟生成（固定反馈多项式G1(x)=1+x3+x10G_1(x)=1+x^3+x^{10}G1​(x)=1+x3+x10）

每来1个码时钟，串行执行4步：

1. 取出bit3、bit10的值：fb=G1[3]⊕G1[10]fb = G1[3] \oplus G1[10]fb=G1[3]⊕G1[10]；
2. 全体寄存器右移1位：G1[2]→G1[1]、G1[3]→G1[2]……G1[10]→丢弃；
3. 反馈值fb填入G1[1]（最左端）；
4. 当前时刻G1输出码片：g1_out=G1[10]g1\_out = G1[10]g1_out=G1[10]。
   循环1023次完成一整周期m序列。

步骤4：G2 m序列生成（固定反馈G2(x)=1+x2+x3+x6+x8+x9+x10G_2(x)=1+x^2+x^3+x^6+x^8+x^9+x^{10}G2​(x)=1+x2+x3+x6+x8+x9+x10）

反馈抽头：bit2、3、6、8、9、10全部异或得到反馈fb；移位逻辑和G1完全一致，输出原始g2_raw=G2[10]。

步骤5：卫星专属相位偏移（核心区分不同PRN）

查表获取当前卫星两组抽头号S1、S2S_1、S_2S1​、S2​，取G2寄存器对应位置比特异或，得到延迟G2输出：
g2_delay=G2[S1]⊕G2[S2]g2\_delay = G2[S_1] \oplus G2[S_2]g2_delay=G2[S1​]⊕G2[S2​]
例：PRN1抽头(2,6)、PRN2抽头(3,7)，每颗卫星查表唯一。

步骤6：Gold码合成（模2加）

单时刻C/A码片：ca_chip=g1_out⊕g2_delayca\_chip = g1\_out \oplus g2\_delayca_chip=g1_out⊕g2_delay

步骤7：周期重置+电平映射（仿真必备）

1. 累计生成1023个码片后，G1、G2强制重置为全1，开启下一个1ms周期；
2. 二进制电平转换：0→+1，1→-1，消除直流分量，用于扩频相乘。

延伸：北斗三号 B2a Gold码（BDS ICD B2a官方标准）

1. 寄存器级数：11级G1/G2，最大m序列周期211−1=20472^{11}-1=2047211−1=2047；
2. 反馈多项式固定，初始值统一规定；
3. 码长截短至2046 chip，码速率10.23Mcps，周期0.2ms；
4. 流程和C/A完全一致：双LFSR移位→G2抽头相位选择→模2加输出测距码；
5. 区分数据支路D、导频支路P：两套独立PRN抽头表，生成正交测距码，分别调制BOC(10,5)副载波。

二、第二类：单LFSR m序列（GPS L2C CM/CL）

基础参数

L2C码时钟511.5kHz，27级单LFSR，反馈多项式固定；

• CM（数据码）：码长10230 chip，周期20ms，每周期寄存器复位；
• CL（导频长码）：码长767250 chip，周期1.5s，超长相干积分抗干扰。

完整分步流程

1. 初始化27bit寄存器为PRN专属初始值（查表，非统一全1）；
2. 每个码时钟：多抽头比特异或得到反馈值；
3. 寄存器全体右移1位，反馈填入首位；
4. 寄存器末位直接输出码片，无第二路、无额外异或；
5. CM计数满10230、CL计数满767250，寄存器重置初始状态；
6. 电平映射0→+1、1→-1，时分复用CM/CL完成L2C扩频码。

关键区别Gold码：只有1套移位寄存器，不存在两路m序列模2加。

三、第三类：Weil数论码（BDS B1C / GPS L1C）

核心：基于有限域二次剩余（勒让德Legendre序列）纯模运算构造，无任何寄存器移位逻辑，和LFSR架构完全割裂，ICD明文区分。
以北斗B1C主码（10230chip）逐条拆解官方数学步骤：

步骤1：定义素数基底p=101，构造勒让德序列L(n)

n为序列索引，n∈[0,100]n \in [0,100]n∈[0,100]
L(n)={0,n mod 101=01,∃k:k2≡n(mod101)(二次剩余)−1,非二次剩余 L(n)= \begin{cases} 0,\quad n \bmod 101 = 0 \\ 1,\quad \exists k: k^2 \equiv n \pmod{101} \quad (\text{二次剩余}) \\ -1,\quad \text{非二次剩余} \end{cases}L(n)=⎩⎨⎧​0,nmod101=01,∃k:k2≡n(mod101)(二次剩余)−1,非二次剩余​
计算方式：遍历k=1~100求k2mod 101k^2 \mod 101k2mod101，标记剩余集合。

步骤2：生成基础Weil短序列W(k;ω)

ω=卫星专属偏移参数（每PRN唯一查表）
W(k;ω)=L(k)⊕L((k+ω) mod 101)W(k;\omega) = L(k) \oplus L\left((k+\omega) \bmod 101\right)W(k;ω)=L(k)⊕L((k+ω)mod101)
长度固定101点，⊕为模2运算（-1等效1、0等效0）。

步骤3：二维矩阵交织扩展至1010点基础序列

将101点W序列填入10行×101列矩阵，行间循环移位交织，打散周期性自相关旁瓣，得到长度1010的基础Weil序列。

步骤4：级联拼接生成母序列（长度10243 chip）

1010×10=10100，再补143点完成10243超长母序列WallW_{all}Wall​。

步骤5：循环截断得到卫星专属10230主码

每颗卫星给定截断起点ppp，循环截取10230个连续码片：
c(n;ω,p)=Wall((n+p−1) mod 10243;ω),n=0,1,...,10229c(n;\omega,p) = W_{all}\left((n+p-1) \bmod 10243;\omega\right),\quad n=0,1,...,10229c(n;ω,p)=Wall​((n+p−1)mod10243;ω),n=0,1,...,10229

步骤6：数据/导频支路分离后处理

1. 导频P支路：直接使用截断Weil主码，无导航电文调制；
2. 数据D支路：主码叠加固定相位偏移序列，与导频正交；
3. 电平统一映射±1，送入QMBOC副载波相乘调制。

核心特征（仿真重点区分）

全程只有模乘、模加、矩阵切片、循环截取；不存在移位寄存器、反馈抽头、异或移位循环，不能复用LFSR代码，必须单独实现数论运算模块。

四、第四类：Memory预存储查表码（Galileo E1 OS E1-B数据/E1-C导频）

完全放弃实时生成算法，离线预计算全量码固化硬件ROM，接收机/仿真器仅寻址读取，无实时运算。

完整工程流程

离线预处理

1. 根据每颗卫星PRN，离线生成完整主码：E1-B/E1-C主码长4092chip；
2. 二进制码压缩为十六进制，存入ROM码表数组；
3. 导频E1-C额外预存25位Neumann-Hoffman辅码CS25。

在线实时读取（仿真/FPGA运行阶段）

1. 码时钟驱动地址计数器自增0~4091，循环溢出归零；
2. 按当前地址直接读取ROM中预存码片bit，作为实时测距码；
3. 导频通道二次编码：每4092个主码周期，乘以1位辅码符号，25个主码周期完成1次辅码循环，相干积分周期拉长至100ms；
4. 电平映射±1，BOC(1,1)调制输出。

优势&局限

优势：实时无运算、资源占用极低、时序稳定；
局限：码长不能随意修改，更换PRN需重新烧录码表，无法实时动态生成。

五、所有测距码通用后处理

无论LFSR Gold/m序列、Weil、查表码，得到0/1二进制序列后统一处理：

1. 电平映射

二进制0 → +1、1 → -1
原因：二进制0/1存在直流分量，和载波相乘会引入载波泄漏；±1双极性序列无直流，频谱干净。

2. 码片上采样/插值

码速率固定（如C/A 1.023Mcps），仿真IQ基带采样率更高（2.6MHz），每个码片插值N个采样点：

• C/A码：fs=2.6MHz，单码片采样点数=2.6/1.023≈2.54点；
• B2a码速率10.23Mcps，采样26MHz，单码片约2.5采样点；
  插值方式：零阶保持（工程仿真最常用，运算量最小）。

3. 副载波叠加（BPSK无副载波，BOC/QMBOC必须乘）

• GPS C/A、BDS B1I：BPSK调制，无副载波，测距码直接乘载波；
• BDS B2a(BOC(10,5))、B1C(QMBOC)、Galileo E1：测距码 × 方波副载波，分裂频谱、提升测距精度、降低和GPS L1干扰。

六、四大码生成架构完整对比

七、仿真实操关键易错点

1. LFSR初值坑：Gold码C/A统一全1；L2C CM/CL是PRN专属初值，不能全1；全0寄存器会锁死序列；
2. 复数IQ采样适配：你仿真用2.6MHz复数基带采样，不受实数奈奎斯特2倍带宽限制；硬件ADC实采样必须≥2倍最高频率；
3. Weil码绝对不能套LFSR逻辑：很多新手写移位寄存器生成B1C，完全违反北斗ICD标准，生成码相关性完全错误，捕获失效；
4. B2a与C/A区分：同为Gold码，但11级寄存器、码长2046、码速率10.23M，不能直接复用C/A的10级LFSR代码；
5. 查表码不能实时生成：Galileo E1如果在线实时计算Gold码会和官方预存码存在微小相位差异，接收机无法匹配捕获。

八、总结

1. 传统老旧民用信号（GPS L1、北斗二号B1I/B3I、北斗三号B2a）全部基于双LFSR Gold码，依靠移位寄存器+异或组合生成；
2. GPS L2C属于简化单LFSR m序列，无两路寄存器；
3. 新一代高精度现代化信号B1C/GPS L1C彻底抛弃移位寄存器，采用纯数论Weil码，底层数学逻辑完全独立；
4. 伽利略E1采用查表存储方案，完全规避实时序列生成；