基于MATLAB的北斗B1I信号模拟器：从星座生成到中频IQ输出

简介：这是一套基于MATLAB的北斗卫星导航信号模拟器源码包，面向卫星导航接收机开发、信号处理算法验证及BDS教学研究人员，可用于模拟B1C/B1I/B2a/B3I等频点的导航信号生成与特性分析。压缩包共13个文件，包含11个MATLAB脚本、1张信号功率谱图像和1份说明文档，整体仅144KB，脚本涵盖扩频码生成、主码/子码合成、BOC调制等核心模块，图像用于展示信号频谱特性，文档帮助快速上手。已有199人下载学习，适合具备一定MATLAB和卫星导航基础、希望深入理解BDS信号结构的读者。源码实现紧扣BDS信号模拟的关键技术，包括坐标系统与轨道参数处理、多路径与传播损耗仿真、接收机前端模型以及性能评估流程，并配有可视化输出。相比仅阅读理论，读者可直接运行代码观察不同频点信号波形与功率谱，结合注释与说明文档实现参数调整和二次开发，是衔接理论与实践的高性价比资源。

1. 拿到“北斗信号模拟器_MATLAB_下载.zip”之前，先弄清它能干什么

你在调北斗B1I接收机的捕获模块，楼顶天线收下来的信号里到底有没有目标卫星，噪声里扒了半天也说不清；想复现一个弱信号场景，又总不能真等卫星过境。这个zip里的东西，就是在电脑里把北斗星座“发射”成一个文件：按你指定的时间、位置、卫星号，生成一段数字中频IQ数据，让接收机算法把它当成从天线+射频前端送进来的信号。它能解决的核心问题是，把捕获、跟踪、电文解调、伪距定位这些环节放到完全可控的条件下调，不受天气、时段和干扰影响。适合做接收机算法验证、实验室教学、课程设计的人——尤其是没有硬件模拟器、又想认真跑通北斗链路的学生和工程师。

2. 北斗信号模拟器的生成链路：从星座到中频IQ的六步

拿到包先别急着点运行。模拟器不是一把噪声撒出去就完事，一条可信的生成链路至少要经过星历、位置、伪距、电文、扩频、合路六个环节。每多懂一个环节，后面调参和排错就少走一段弯路。

2.1 星历与星座：模拟器的“时间基准”

第一步永远是确定“此刻天上哪颗星在什么位置”。星历来源常见有两种：广播星历和精密星历。做接收机算法验证我一般用广播星历，精度是米级，足够把捕获跟踪链路跑通；更关键的是，接收机最终从电文里解出来的就是广播星历参数，模拟器用它生成、接收机用它解算，两端对得上，验证闭环才成立。精密星历事后才能拿到，适合做高精度定位研究，不适合日常调试。

拿到的RINEX导航文件里是一组开普勒参数。这里有个特别容易踩的坑：文件里的半长轴是平方根格式（sqrtA），单位是米开方，很多人直接拿来当半长轴用，后面伪距误差会到几十公里量级。另外模拟器内部的时间基准要统一用北斗时（BDT）。BDT和GPS时差14秒，这14秒里卫星已经在轨道上跑了大约55公里，伪距直接崩掉。

% 解析RINEX 3.03北斗导航电文，按ICD-B1I字段对应 % rinexFile: RINEX文件路径; prn: 目标卫星号 eph = extractBDSBroadcastEph(rinexFile, prn); % eph.sqrtA -> 半长轴平方根, 使用时要平方 % eph.e -> 偏心率 % eph.i0 -> 轨道倾角 % eph.OMEGA0 -> 升交点赤经 % eph.M0 -> 平近点角 t_bdt = currentBDT(); % 统一用北斗时，别用GPS时 [posX, posY, posZ] = calcSatPosEcef(eph, t_bdt);

这段代码的逻辑是标准的“星历转位置”：拿到开普勒根数后，解开普勒方程求平近点角到真近点角，再经过一系列坐标旋转得到地心地固系（ECEF）坐标。参数上，sqrtA平方后才是半长轴；时间不能用系统UTC直接代入，要先转成BDT。

2.2 伪距、多普勒与电文：让信号“动”起来

有了卫星坐标，下一步是算用户到卫星的几何距离。模拟器里伪距一般是几何距离加卫星钟差，再按需叠加电离层、对流层和相对论修正；教学版模拟器通常只保留几何距离和钟差，够用。真正的坑在多普勒上：多普勒是距离变化率的负数，不是拿起始位置算一个固定值就完了。动态场景下，同一段信号里每个采样点的伪距都在变，如果只更新到毫秒级甚至更粗，接收机跟踪环路一拉就断。

导航电文的组织也直接决定模拟器能不能用于解调验证。B1I的D1电文是50bps，一个子帧300比特、历时6秒，帧结构依次是8位Preamble、子帧ID、周内秒和星历参数，再带上BCH校验。很多简化模拟器直接拿随机比特填电文，捕获测试没问题，但你要解星历做定位，这种包就废了。

% 生成一个B1I D1子帧，只列结构骨架 preamble = [1 0 0 0 1 0 1 1]; % 0x8B subId = de2bi(3, 4, 'left-msb'); % 子帧3 bits = [preamble, subId, howBits, ephBits]; % 拼接 % BCH(15,11)校验：按ICD要求在30位字边界上计算 words = bch15_11_encode(reshape(bits, 24, []), 'left-msb'); if enableNh nbi = reshape(repmat(words', 20, 1), [], 1); % NH码二次扩展 end

D1电文的特别之处在于，50bps的比特流先经过20位NH码做二次编码，每个电文比特被扩展成20个1ms码片，再和测距码相乘。接收端必须先把NH码边界找回来，把20个码片累加成1个电文比特，否则解出来全是噪声。D2电文是500bps、没有NH码，千万别拿D1的流程硬套。

2.3 测距码、NH码与BPSK调制：把数据变成射频前的形态

北斗B1I的测距码是Gold码，码长2046、码率2.046Mcps、周期1ms，由两个11级移位寄存器G1和G2生成，每颗卫星对应G2不同的抽头延迟。MATLAB里用通信工具箱的生成函数最省事，手写移位寄存器也不难，关键是码表和ICD要一致，错一个抽头相关峰出不来。

扩频之后是调制。基带信号乘上电文比特和测距码，就得到一个码率2.046Mcps的基带波形。接下来两条路：要么零中频，直接输出I/Q两路基带；要么搬到低中频，比如4.092MHz，更接近真实接收机的下变频输出。我一般选后者，毕竟真实天线下来都有中频，捕获器也按中频模型写。多颗卫星各自生成后按功率加权累加，就得到最终IQ。

% 生成单颗卫星基带信号的最小骨架 code = generateB1ICode(prn, 'g2Delay', delayTable(prn)); code = code * 2 - 1; % 0/1 -> +1/-1 % 电文比特nbi(已含NH码扩展)映射到码片 spread = reshape(repmat(code', 1, 20)', 1, [])'; baseband = spread .* nbi; % BPSK扩频 % 搬到中频，fs为采样率 t = (0:length(baseband)-1) / fs; iq = baseband .* exp(1j * 2 * pi * f_if * t);

逻辑说明：spread把一个码周期复制成20份，对应一个电文比特占20个码周期；baseband完成扩频；最后乘复指数完成数字上变频。参数上，fs要和码率有合理的倍数关系，比如20Msps下每个码片约9.78个采样点，相位分辨率和捕获精度都比较平衡。非要抠每个码片几个点，这类细节我会放到“参数设置”那一章展开。

3. 在MATLAB里跑通下载包：目录结构与最小复现步骤

压缩包下载下来，多数人的第一反应是直接双击主脚本，然后面对一片报错发愣。我的建议是，先花两分钟把目录结构看明白，再决定怎么跑。

3.1 拿到zip先做三件事：解压、看README、查路径

zip解压到没有中文和空格的路径是最基本的一条。MATLAB对中文路径总体能处理，但有些旧版脚本里硬编码了相对路径，中文目录名一掺和就找不到了。解压后别急着开MATLAB，先在命令行里看一遍目录结构，有没有README，哪个是主入口文件，函数文件放哪个子目录。

unzip 北斗信号模拟器_MATLAB_下载.zip -d bds_sim cd bds_sim ls -R | head -60

如果目录里有README或快速上手指南，先读完再动手。很多包的主脚本要依赖子文件夹里的函数，但MATLAB当前路径不对，会出现“未定义函数或变量”一类的报错，这其实是路径问题而不是代码问题。我习惯在浏览器里先在文件列表里点开主脚本看一眼，确认入口位置，再在MATLAB里addpath把根目录和src目录加进去。用which命令可以快速确定某个函数是否在路径上。

3.2 主脚本与最小参数集：先让默认配置跑出IQ文件

典型的目录布局长这样，虽然不同包有差异，但大体逃不出这个骨架：

bds_sim/ config/ % 参数配置 data/ % RINEX星历、输出IQ src/ % 信号生成函数 run_sim.m % 主入口

拿到主脚本先把参数区浏览一遍，不急着改。最小参数集就四个：采样率、信号时长、卫星号列表、输出文件名。我第一次跑通一个模拟器包，永远是用默认参数跑一遍，确认它能吐出一个IQ文件，然后立刻做频谱和相关峰验证，确认“货是真的”，再去调后面那些花活。上来就改参数，出了问题分不清是包的问题还是自己改的问题，这种排错体验最折磨人。

% run_sim.m 的典型调用方式 cfg.sampleRate = 20e6; % 采样率：20 MHz cfg.durationSec = 10; % 信号时长：10 秒 cfg.ifFreq = 4.092e6; % 中频：4.092 MHz cfg.prnList = [1 3 7 11]; % 卫星PRN列表 cfg.outputFile = 'b1i_20m_10s_int8.iq'; iq = bdsSignalGenerator(cfg); % 返回复数基带/中频 writeComplexIQ(cfg.outputFile, iq, 'int8'); % 落盘

输出写盘时有一个常见翻车点：直接对复数数组cast成int8会把小信号截成一片噪声，应该按IQ两路分别乘127取整，也就是幅度归一化后再量化。参数层面，采样率20MHz是老规矩，B1I码率2.046Mcps，一个码片大约9.8个采样点，够捕获使用，文件也控制在合理范围；10秒时长足够后续做电文解调和简易定位验证。

4. 模拟器关键参数怎么设：采样率、时长、卫星数与量化位宽

参数设置决定了模拟器输出能不能被接收机用、能用多久、磁盘放不放得下。这几个参数相互牵制，分开调容易顾此失彼。

4.1 采样率与中频：决定文件大小与捕获性能

B1I信号的主瓣带宽约4.092MHz。采样率低了，码片分辨率不够，捕获的相关峰变钝；采样率高了，文件大小上涨，但捕获精度提升有限。20MHz是很多接收机前端的典型采样率，兼顾精度和存储；62MHz这类高采样率常见于SDR前端，如果模拟器主要喂软件接收机，按实际硬件的Fs来输出就行，省得接收链路里再内插。

中频的选择直接影响后端处理。零中频最简单，但直流偏置直接落在信号带内，频谱上中心一大坨DC分量，抓不准；4.092MHz低中频最常用，好处是它是码率2.046Mcps的整数倍，采样率20MHz时信号落在第一奈奎斯特区内，后续数字下变频干净利落。也有模拟器支持带通采样，降低Fs的同时保持高信息率，那是进阶玩法。

文件大小 = 采样率 × 量化字节数 × 2（I/Q两路）× 秒数。设参数前先按这个公式算一笔，我吃过亏，默认62Msps加int16生成20秒直接写满硬盘。

% 参数配置：采样率、中频、量化 cfg.sampleRate = 20e6; % 码片采样点约9.78，够用 cfg.ifFreq = 4.092e6; % B1I码率整数倍 cfg.quantBits = 8; % int8: 文件减半

参数说明：中频建议贴4.092MHz的整数倍，量化不低于8bit，低于8bit信号会被削波，弱星直接跳不出来。

4.2 信号时长与动态场景：决定你要跑多久

时长这个问题没有统一答案，取决于你想验证什么。捕获测试给1到2秒就够，捕获器只需要几毫秒数据；解调电文至少要完整一个子帧，B1I子帧是6秒，我习惯给10秒，留出收敛余量；做定位解算测试要给30秒以上，因为星历参数在多个子帧里循环播发，时间太短凑不齐一套星历。

动态场景是另一个维度。模拟车、无人机运动时，伪距不是恒定值，每个采样点的伪距由初始距离、径向速度、径向加速度决定。只做常数伪距的模拟器在静态捕获测试里没问题，但跟踪环路一上，载波频率一直固定，环路以为信号频漂，锁不住。动态场景建议给伪距加一个二阶模型，粗糙但有效。

% 动态场景：伪距按时间二次多项式更新 cfg.dyn.p0 = 2.3e7; % 初始伪距，约2.3万公里 cfg.dyn.v = -120; % 径向速度，m/s，负号表示靠近 cfg.dyn.a = 5; % 径向加速度，m/s^2 for k = 1:nSamples t = k / fs; rho(k) = cfg.dyn.p0 + cfg.dyn.v * t + 0.5 * cfg.dyn.a * t^2; end

逻辑说明：每个采样点更新一次伪距，多普勒自然通过相位差分体现，不用单独算。参数上，径向速度决定多普勒偏移量，加速度决定多普勒变化率，这两项和接收机环路带宽要匹配，否则信号超出环路跟踪能力。长时长场景要用分段生成，比如10秒拆成10段，每段生成1秒写盘，避免内存一次装下整段IQ。

4.3 卫星数与功率控制：影响接收机冷启动速度

卫星数量直接决定接收机能不能完成定位。捕获验证用1颗卫星完全够；定位解算至少4颗；和真实环境对标用8到12颗。卫星增多后，多颗星信号累加，互相关干扰上升，弱星会被强星压制，这个在真实环境里存在，模拟器里同样存在，正好用来验证接收机在复杂环境下的性能。

功率控制的核心是CN0，单位dB-Hz。真实接收机天线下，强信号50dB-Hz左右，弱信号35dB-Hz也能跟踪。模拟器一般不给每颗星单独设幅度，而是给一个CN0表，生成时按CN0换算成幅度，再叠加到总信号上。

% 每颗星的CN0配置，单位dB-Hz，顺序和prnList对应 cfg.cn0 = [45 44 43 40]; % 换算成采样域幅度；这里假设噪声底已按int8量化归一 A = sqrt(2) * 10.^((cfg.cn0 - 30) / 20); signal = sum(A .* sigList, 2) + noise; % 多星合路加噪声

参数说明：A越大信号越强，CN0差10dB，幅度差约3.16倍。不要把所有星都设成同样功率，真实场景里仰角不同、CN0天然不同，全一样反而测不出接收机的灵敏度差异。弱信号测试把CN0压到30到35dB-Hz，捕获搜索步长要相应缩小到200Hz左右，否则弱星的多普勒bin一格扫过去，能量损失太多。

5. 避坑指南：北斗模拟器常见翻车点与排查

用这类zip包，最大的成本不在跑通，而在跑通之后信号能不能被接收机正经解调。下面几条是常见的踩坑记录，全是我自己栽过的跟头。

5.1 现象：解压后脚本报错“未定义函数”

刚解压就在MATLAB里双击主脚本，弹窗“未定义函数或变量”。原因基本是两个：MATLAB路径没切到解压目录，函数文件不在搜索路径上；缺工具箱，比如模拟器调用了通信工具箱的Gold序列生成函数或信号处理工具箱的滤波函数，本机没装。

排查用which 函数名看能不能定位，用ver看工具箱清单。确认缺工具箱后，要么安装对应工具箱，要么把调用改成自己手写的移位寄存器生成函数。我的习惯是拿到包先补齐工具箱再跑，MATLAB版本也别太老，较新的版本对RINEX读取和文件I/O支持都好一些，能少很多兼容性报错。

5.2 现象：IQ文件巨大，磁盘瞬间写满

模拟器默认参数可能配到62.5MHz采样率加int16量化，10秒就是2.5GB。很多包为了“通用”把默认值设得偏高，新手不看参数直接跑，一下就写满磁盘，主脚本还在那报“磁盘空间不足”。

跑之前按采样率、量化位宽、时长算一遍文件大小。2秒、20MHz、int8是80MB，适合先验证链路；确实需要长时长，就用分段写盘，一段一段生成追加写入，别让数组一直在内存涨着。20MHz对B1I足够，除非你要测接收机在高采样率下的行为，否则没必要上62MHz。

5.3 现象：捕获相关峰不明显或完全无峰

信号生成出来，丢进捕获器却找不到峰值。先别怀疑包坏了，九成是参数不匹配。第一个是多普勒搜索没覆盖，MEO卫星在静态用户下多普勒范围大约±5kHz，你搜索只设了±2kHz，峰值出不来；第二个是搜索步进太大，1kHz步进会让信号落在相邻bin的缝隙里，信号弱一点就被噪声淹没；第三个是本地码生成端和信号端采样率不一致，比如信号按20.48MHz生成，本地码却按20MHz生成，码相位根本对不齐。

排查时先把场景拉回最简单：单颗卫星、静态用户、零多普勒、不加噪声，能出峰再加多普勒、加卫星、加噪声。这种“由简到繁”的排除法比对着参数瞎猜高效得多。

5.4 现象：电文解调全错但捕获正常

捕获出峰说明码相位和多普勒是对的，但解电文全错，常见原因是NH码同步没做。D1电文每个比特被NH码扩展成20个1ms码片，接收机必须先找到NH码组的边界，把20个码片累加成1个电文比特，直接按50bps去解，出来的全是乱码。

另一个高频原因是时间基准混用。模拟器内部按北斗时生成，接收机解星历却用GPS时，两颗系统时间差14秒，所有子帧的周内秒和解出的星历参数都对不上。排查时第一步找Preamble，B1I的Preamble是0x8B，能稳定找到说明子帧边界对；第二步检查NH码边界；最后核对时间基准是不是BDT。还有一个隐蔽坑：星历过期。RINEX里的广播星历外推超过2小时，伪距误差会到公里级，解出来的定位自然错得离谱，换新鲜的星历文件就好。

5.5 现象：定位结果差几十公里

电文能解，伪距能算，但定位解算出来差几十公里——这问题不在接收机，在模拟器端。最常见的三个原因：星历参数单位没用对，半长轴忘了平方，星钟参数忘了乘光速；坐标框架混用，模拟器按CGCS2000算，接收机按WGS-84转经纬度，虽然两者米级一致，但如果你自己写转换时用的椭球参数不对，偏差会到几百米甚至更大；伪距里只放了几何距离，漏了卫星钟差和地球自转修正。

我的排查习惯是把伪距残差打出来。哪颗星残差特别大，就重点查那一颗星的电文参数，多半是那一颗的参数解析错了。模拟器端设的用户坐标也值得核一遍，ECEF坐标转经纬高时常见单位写错，米和千米之间差三个数量级，定位结果直接飞出天际。

6. 验证模拟器输出的四步技巧：从相关峰到定位解算

生成IQ文件之后，要习惯性地验证它“是真的”，别拿到手里就闷头调环路。我从四个层次查，每一层过了，模拟器输出的可信度才成立。

第一步，也是我最先做的，是捕获验证。取1ms数据做相关，本地码和信号做循环互相关，找峰值。

% 读IQ，取1ms数据做循环相关 iq = readComplexIQ('b1i_20m_10s_int8.iq', 20e6, 1); % 读1秒 seg = iq(1:20e6/1000); % 1ms [corr, lag] = xcorr(seg, localCodeB1I(prn, fs)); plot(lag, abs(corr)); % 主峰应明显高于旁瓣

这一步过了，码相位和多普勒基本没问题。第二步是验电文：找Preamble，查子帧ID是不是1到5，看周内秒是否递增，能对上就说明NH码和解调链路是通的。第三步是伪距单点定位：把模拟器已知的伪距喂给接收机的定位解算，算出来的用户位置和模拟器设定位置差在10米量级，才算闭环。第四步如果有条件，和商用模拟器或者真实射频前端采集的数据做对比，看频谱包络和相关峰形状，这是最后一道后悔药。

四步全过，这个模拟器的输出才值得拿去做跟踪和定位的算法调优。我的习惯是，每次改场景参数，都要把第一步和第三步重跑一遍再交给环路；不然你永远说不清是算法不行，还是信号源不可靠，整个过程就成了黑匣子调试。希望这些经验能帮你在北斗信号模拟器这条路上少翻几次车。

本文还有配套的精品资源，点击获取