简介这是一套基于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码二次扩展 endD1电文的特别之处在于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的同时保持高信息率那是进阶玩法。采样率I/Q量化时长文件大小20Mspsint810s400MB20Mspsint1610s800MB62Mspsint810s1.24GB62Mspsint1610s2.48GB文件大小 采样率 × 量化字节数 × 2I/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秒所有子帧的周内秒和解出的星历参数都对不上。排查时第一步找PreambleB1I的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米量级才算闭环。第四步如果有条件和商用模拟器或者真实射频前端采集的数据做对比看频谱包络和相关峰形状这是最后一道后悔药。四步全过这个模拟器的输出才值得拿去做跟踪和定位的算法调优。我的习惯是每次改场景参数都要把第一步和第三步重跑一遍再交给环路不然你永远说不清是算法不行还是信号源不可靠整个过程就成了黑匣子调试。希望这些经验能帮你在北斗信号模拟器这条路上少翻几次车。本文还有配套的精品资源点击获取 
长河君 · 2026/10/27 1:32:58 · 阅读 1.2k 基于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码二次扩展 endD1电文的特别之处在于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的同时保持高信息率那是进阶玩法。采样率I/Q量化时长文件大小20Mspsint810s400MB20Mspsint1610s800MB62Mspsint810s1.24GB62Mspsint1610s2.48GB文件大小 采样率 × 量化字节数 × 2I/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秒所有子帧的周内秒和解出的星历参数都对不上。排查时第一步找PreambleB1I的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米量级才算闭环。第四步如果有条件和商用模拟器或者真实射频前端采集的数据做对比看频谱包络和相关峰形状这是最后一道后悔药。四步全过这个模拟器的输出才值得拿去做跟踪和定位的算法调优。我的习惯是每次改场景参数都要把第一步和第三步重跑一遍再交给环路不然你永远说不清是算法不行还是信号源不可靠整个过程就成了黑匣子调试。希望这些经验能帮你在北斗信号模拟器这条路上少翻几次车。本文还有配套的精品资源点击获取