简介基于Matlab与Simulink的CDMA仿真系统源码包面向电子信息工程、计算机、数学等专业的学习者用于理解码分多址通信原理、序列生成及扩频解扩流程。资源内同时包含可运行源码与配套仿真数据适合课程设计、毕业设计或课题研究作为参考资料。包内共18个文件含8个.m脚本如mseq.m、goldseq.m、spread.m、despread.m等分别覆盖m序列与Gold序列生成、扩频、解扩及系统主流程6个.mat数据文件用于保存仿真输入输出便于复现实验另有Simulink模型相关文件.slxc、.xml及说明文件整体体积约97KB。代码结构清晰重点演示CDMA基带信号处理的关键环节可直接在Matlab环境中加载学习。该资源已有365人学习属轻量紧凑的示例包适合具备一定Matlab基础、希望自行调试代码并扩展功能的读者对照代码快速梳通CDMA系统从序列产生到解调输出的完整实现路径。1. CDMA仿真系统把码分多址原理变成能跑的模型做 CDMA 仿真系统最怕的不是原理看不懂而是原理写进 Simulink 之后模型跑出来的误码率怎么调都不对。这套基于 Matlab 和 Simulink 的 CDMA 仿真系统正好把最难啃的这块拆开了发射链路、PN 序列扩频、信道模型、接收端解扩判决、误码率统计一条龙拉通源码和数据都齐全。适合正在做通信课程设计、本科毕设或者刚入职想补 CDMA 物理层基础的人。我按一线工程师拆项目的习惯带你把每个模块的作用、参数、坑全部过一遍看完你能直接照着复现也能改成自己的实验参数跑出新结果。2. 系统拆分发射链路、信道模型与接收端的 Simulink 搭建2.1 从原理到模型动手前先想清楚的三个问题拿到这套资源别急着打开模型就跑。我每次接手一个通信仿真项目都会先花十分钟画一张数据流草图明确三个问题速率怎么对齐、同步怎么实现、误码率在哪一级统计。CDMA 的核心是扩频信息比特先乘以高速 PN 序列把窄带信号扩展成宽带信号接收端再用同步的 PN 序列做相关解扩把宽带信号重新收拢成窄带信号同时获得处理增益。这个过程在 Simulink 里落地第一个要解决的问题是采样率。信息速率假设是 Rb扩频因子是 SF那么码片速率 Rc SF × Rb。这意味着扩频模块前后的信号速率不一致Simulink 模型里必须用明确的样本时间Sample Time区分比特周期和码片周期否则数据流会乱。第二个问题是同步。实际 CDMA 接收机要做 PN 码捕获和跟踪但仿真模型里最常见做法是让发送端和接收端共用同一套时钟也就是用同一个 PN 序列生成器接收端直接“白嫖”发送端的码相位。这样做的好处是省掉同步模块把精力集中在扩频、解扩和误码率统计上代价是如果真想研究同步算法需要手动把接收端 PN 序列延迟几个码片观察相关峰变化。第三个问题是误码率统计口径。误码率应该统计判决之后的信息比特不是扩频后的码片。很多新手把 Error Rate Calculation 模块接在 BPSK 解调之后、解扩之前统计的是码片错误率这和系统真实误码率是两码事数值看起来会莫名其妙地高。2.2 发射链路PN 序列生成、扩频与 BPSK 调制发射链路在 Simulink 里主要由三块拼起来随机信息源、PN 序列生成器、扩频与调制模块。这套系统用的是 m 序列PN Sequence Generator 模块里需要配置多项式系数和寄存器长度。% 用脚本配置PN序列生成器的本原多项式 % 特征多项式: x^5 x^2 1对应多项式系数 [1 0 0 1 0 1] poly [1 0 0 1 0 1]; % 寄存器初始状态不能全0否则输出恒为0 init_state [1 0 0 0 0]; % 在Simulink模型里通过set_param配置PN模块 set_param(cdma_model/PN_Sequence, Polynomial, mat2str(poly)); set_param(cdma_model/PN_Sequence, InitialState, mat2str(init_state));这段代码的核心是 m 序列的两个硬约束多项式系数必须包含最高次项和常数项寄存器初态不能全零。常见做法是选 n5 或 n7 的 m 序列周期分别是 31 和 127。扩频因子 SF 通常取等于 m 序列周期比如 SF31对应的处理增益是 10×log10(31) ≈ 15 dB这个增益会直接体现在误码率曲线上——同样 Eb/N0 下扩频系统的误码率会比未扩频系统明显低。扩频操作在 Simulink 里的实现方式有两种。一种是用 Repetition重复模块把每个信息比特复制 SF 次再与 PN 序列逐码片相乘另一种是直接用 Gain 模块配合 PN 序列做乘法。这套系统看起来采用的是第二种思路PN 序列输出的 ±1 序列直接与 BPSK 映射后的数据相乘相当于每个比特和 SF 个码片按位相乘完成扩频。BPSK 调制用的是 BPSK Modulator Baseband 模块输入 0/1 比特映射为 1/-1。这里有一个容易忽略的细节BPSK 解调器输出是软值带幅度信息如果后续直接接解扩和判决必须搞清楚软值是正负号还带大小。如果是带大小的软值在 Rayleigh 信道下幅度起伏会直接影响判决需要做归一化。2.3 信道与接收AWGN、瑞利衰落与解扩判决信道模型是仿真系统里最值得替换的部分。系统自带两个信道配置AWGN 和瑞利平坦衰落。AWGN 信道用 AWGN Channel 模块实现参数主要是信噪比 SNR 或者 Eb/N0瑞利信道用 Multipath Rayleigh Fading Channel 模块需要配置多径延迟和增益。接收端的核心是解扩。解扩在数学上就是本地 PN 序列与接收信号的逐点相乘再在一个比特周期内累加最后做符号判决。Simulink 实现时一般用 Product 模块乘上同一个 PN 序列再用 Discrete-Time Integrator 做比特周期内的累加最后用 Sign 模块判为 ±1对应回 0/1 比特。% 配置AWGN信道的信噪比单位dB set_param(cdma_model/AWGN_Channel, SNR, 10); % 配置瑞利信道参数 % 多径延迟 [0 2e-6]对应两个径相差2微秒 set_param(cdma_model/Rayleigh_Channel, ... DelayVector, [0 2e-6], ... GainVector, [0 -3]);这里参数的含义分别是SNR 是接收端信号功率与噪声功率比值单位 dBRayleigh 信道的 DelayVector 是每条径相对第一条径的延迟秒数GainVector 是每条径的平均增益 dB。常见做法是让第一条径增益 0 dB后面径逐渐衰减。如果仿真瑞利信道却不加解扩误码率会非常难看——这正是 CDMA 多径分集要解决的问题。系统里自带一个简单的双径合并等增益合并把两条径的解扩结果直接相加再判决这已经能体现 RAKE 接收的基本思想。3. 在 Simulink 里跑通仿真脚本驱动、参数设置与误码率统计3.1 用脚本驱动仿真一套参数跑完整个实验这套系统的一个友好之处在于模型配置和仿真运行都可以用脚本控制。我一般会把所有可变参数集中在脚本开头模拟一个“参数面板”改参数不用进 Simulink 界面里翻模块。% 系统参数面板 SF 31; % 扩频因子等于m序列周期 Rb 1000; % 信息速率单位bps Rc SF * Rb; % 码片速率 numBits 5000; % 仿真总比特数建议至少5000 EbN0_dB 0:2:12; % Eb/N0扫描范围 % 计算采样率与仿真时长 samplePerChip 1; % 每个码片一个采样点 Ts 1 / (Rc * samplePerChip); % 采样时间 simTime numBits * SF * Ts; % 总仿真秒数 % 写入模型参数 set_param(cdma_model/PN_Sequence, SampleTime, num2str(Ts)); set_param(cdma_model/AWGN_Channel, SNR, num2str(EbN0_dB(1))); set_param(cdma_model, StopTime, num2str(simTime));这段脚本的作用是把系统参数和模型解耦想换码速率、换扩频因子、换仿真时长只改脚本第一段就行。注意 samplePerChip 这个变量它决定每个码片采几个点默认 1 就够了但如果后面接脉冲成型滤波器或频谱分析建议改成 4 或者 8否则频谱会有混叠。仿真停止时间必须按总比特数计算。当初我犯过一个低级错误固定 StopTime 为 1 秒结果信息速率一改统计的比特数就变了误码率曲线波动大得没法看。正确的做法是让仿真时长和比特数挂钩保证每次实验统计的符号数量一致。3.2 误码率统计从误码到曲线的完整链路Simulink 里统计误码率最直接的工具是 Error Rate Calculation 模块把它接在判决输出和原始信息比特之间模块输出三个数误码率、误码个数、对比的总比特数。用 To Workspace 模块把结果存到 Matlab 工作区仿真结束后统一处理。% 仿真运行结束后读取误码率结果 sim(cdma_model); ber BER_Out(1); % BER_Out是To Workspace输出的数组 numErrors BER_Out(2); % 对比理论BPSK误码率 berTheory 0.5 * erfc(sqrt(10^(EbN0_dB(1)/10))); % 打印对比结果 fprintf(Eb/N0 %.1f dB, 仿真BER %.5f, 理论BER %.5f\n, ... EbN0_dB(1), ber, berTheory);误码率统计最关键的是对齐比较的两路信号。Error Rate Calculation 模块的 Tx 端口和 Rx 端口必须延迟对齐否则统计出的误码率不是系统性能而是“信号没对齐造成的虚假误码”。系统中所有模块都带 Sample Time一旦前后速率不匹配模块输出会产生固有延迟需要在 Rx 路径上加 Delay 模块补偿。经验值是看仿真输出的误码个数如果误码率稳定在 0.5 附近说明对齐彻底错了信号完全错位。另一个细节是误码率统计的冷启动问题。仿真前几个比特因为模块初始状态未稳定误码率会偏高。常见做法是在脚本里设置一个“丢弃前 100 比特”的逻辑或者干脆把仿真比特数设大一些让前几百个比特的影响占比小到可以忽略。3.3 参数表扩频因子、信噪比与仿真时长的设定原则把系统里最常用的参数列成一张表方便对照修改参数项取值范围设定原则扩频因子 SF31、63、127越大处理增益越高但带宽占用越大仿真越慢信息速率 Rb100 ~ 10000 bps仿真模型里速率不影响结果但影响仿真速度每码片采样数1、4、8只观察误码率用 1做频谱或眼图建议 4 以上Eb/N0-2 ~ 14 dB根据期望误码率下限设定低于 -2 dB 误码率基本 0.5仿真总比特数5000 ~ 50000目标误码率越低需要统计的比特数越多多径延迟0 ~ 100 码片时长延迟超过码片周期后RAKE 合并效果才明显这套参数表几乎决定了仿真结果可信度。比如你目标误码率是 1e-4按统计学要求至少要统计 100 个错误比特才有意义也就是至少跑 100 万比特的仿真。虽然是 Simulink 模型100 万比特可能要跑十几分钟但低于这个量级跑出的 1e-4 完全不可信属于“样本太少纯属玄学”。4. 避坑与排查CDMA 仿真中最常见的五个翻车点4.1 误码率先高后低解扩时刻没对齐现象误码率曲线在低信噪比时接近 0.5信噪比升高后突然下落但始终比理论曲线差很多。原因接收端的本地 PN 序列和接收信号的码相位没有对齐解扩变成了“错位相乘”扩频增益完全没有发挥信号被自己扩频码的旁瓣干扰压制。解决检查接收端 PN 序列生成器的 InitialState 是否和发送端完全一致再检查 PN 序列模块的输出是否有额外延迟。在接收端 Path 上加一个 Delay 模块用仿真前 100 个比特做一次手动对齐让误码率先降到正常水平再调其他参数。这是 CDMA 仿真里最“玄学”的问题十有八九是这里。4.2 频谱显示全是毛刺采样率与码片速率关系搞错了现象用 Spectrum Scope 看发射信号频谱本该是主瓣清晰的扩频谱结果全是尖刺主瓣形状都看不出来。原因频谱分析仪需要一个合适的采样率如果每个码片只采 1 个点频谱会周期性延拓出现明显混叠。解决把 samplePerChip 改成 4 或 8也就是让码片速率乘上过采样倍数后再进入滤波器频谱就干净了。注意改了这个参数所有模块的采样时间都要跟着变不能只在频谱仪前面加一个插值模块。4.3 多用户干扰为零扩频码正交性没用上现象加了第二个用户的信号之后误码率竟然和单用户完全一样没有任何恶化。原因两个用户可能用的是完全相同的 PN 序列或者同一序列的不同相位相乘之后被自己解扩时完全消除了自干扰为零这本质上是同一个信号叠加了两次不是多用户干扰。解决给两个用户分配两条不同相位的 m 序列或者直接用 Gold 码。Gold 码互相关性可控是 CDMA 多用户仿真的标准选择。在 Simulink 中可以用两个 PN Sequence Generator 模块配置不同 InitialState 来生成不同相位的序列。4.4 仿真跑不动仿真时长与采样率组合失控现象模型仿真速度越来越慢跑几分钟只走了总时长的 1%。或者仿真直接卡死。原因仿真时长是按秒设置的而采样率是码片速率乘上过采样倍数。SF127、每码片采样 8 点、跑 5 万比特时总采样点数是 127×8×50000 5080 万点Simulink 每个采样点都要执行几十个模块计算量直接爆炸。解决先用 SF31、每码片采样 1 点跑通流程确认无误后再往大参数调。调大参数后先用 1000 比特验证速度速度能接受再跑全量。血泪经验能用 5000 比特验证的结论不要等 5 万比特跑完才看结果。4.5 误码率曲线不随信噪比下降判决阈值与幅度归一化问题现象Eb/N0 从 0 dB 调到 12 dB误码率始终徘徊在 0.1 左右曲线是一条平线。原因BPSK 解调器输出的是带幅度的软值经过瑞利信道后幅度起伏大直接用固定阈值 0 判决没有做幅度归一化低信噪比时噪声压过信号但高信噪比时深衰落瞬间的符号照样判错。解决在解扩累加之后用一个 Gain 模块或者 Divide 模块按信号平均幅度做归一化把软值归一化到 ±1 附近再判决。瑞利信道下不归一化误码率存在“地板效应”加再大信噪比也没用。5. 进阶把仿真系统变成研究工具的三件事5.1 批量跑误码率曲线脚本化参数扫描单点跑仿真只能验证一个 Eb/N0 下的误码率做研究必须画一整条 BER 曲线。常见做法是把仿真放进 for 循环每次只改 AWGN 信道的 SNR把误码率存成数组最后统一绘图。% 批量扫描Eb/N0绘制BER曲线 EbN0_dB 0:2:12; ber_sim zeros(size(EbN0_dB)); for k 1:length(EbN0_dB) set_param(cdma_model/AWGN_Channel, SNR, num2str(EbN0_dB(k))); sim(cdma_model); ber_sim(k) BER_Out(1); fprintf(Eb/N0 %.1f dB, BER %.5f\n, EbN0_dB(k), ber_sim(k)); end % 理论曲线对比 ber_theory 0.5 * erfc(sqrt(10.^(EbN0_dB/10))); semilogy(EbN0_dB, ber_sim, bo-, EbN0_dB, ber_theory, r*-); grid on; xlabel(Eb/N0 (dB)); ylabel(误码率); legend(仿真, 理论);这段脚本跑起来慢是正常的实时看打印输出注意趋势即可。理论上 AWGN 信道下 BPSK 的误码率曲线和 CDMA 仿真曲线应该基本重合差别在于 CDMA 的处理增益。如果画的曲线比理论低好几个数量级先别高兴检查是不是统计口径错把码片误码当成比特误码。5.2 可视化升级眼图、星座图与误码分布误码率只是一维指标真正定位问题还是要看眼图和星座图。把前端的信号接一个 Constellation Diagram 模块可以实时观察 BPSK 星座点的聚集情况。AWGN 信道下星座点应该是两团球簇瑞利信道下两团球簇会“糊”成一片这正是深衰落的直观表现。眼图更能说明问题。在 BPSK 调制之前、扩频之后加一个发射滤波器比如升余弦滤波器再在接收端加匹配滤波器眼图就会呈现出教科书里的形状。这套系统如果不加滤波眼图张不开但频谱利用率低加了滤波眼图好看了符号间干扰又冒出来。我一般用眼图观测接收端解扩前的信号质量用来判断信道损伤是噪声主导还是失真主导。5.3 从单用户到多用户Gold 码与干扰控制做完单用户系统下一步就是加用户。多用户 CDMA 的实质是同时同频传输靠码分区分用户。用 Gold 码替代 m 序列两个用户用不同生成多项式或不同相位的 Gold 序列接收端用对应码解扩时其他用户的信号表现为类似噪声的干扰。多用户系统的误码率公式不再只是热噪声还包含多用户干扰项。当用户数增加时BER 曲线会抬升这是因为干扰功率随用户数线性增长。我在仿真里通常让用户数从 2 加到 6观察 BER 曲线之间的距离这比任何理论公式都直观。从那以后我每次接手新的 CDMA 仿真模型都会强制自己先跑一遍 0~12 dB 的单用户 AWGN 曲线确认解扩同步、误码统计口径都正常之后再往模型里加信道、加多用户、加 RAKE 合并。模型每加一个复杂度就重复一遍这条基线流程出问题时能快速定位是哪一步引入的。希望这套拆解思路帮到你。本文还有配套的精品资源点击获取 基于Matlab/Simulink的CDMA仿真系统拆解与误码率分析
简介基于Matlab与Simulink的CDMA仿真系统源码包面向电子信息工程、计算机、数学等专业的学习者用于理解码分多址通信原理、序列生成及扩频解扩流程。资源内同时包含可运行源码与配套仿真数据适合课程设计、毕业设计或课题研究作为参考资料。包内共18个文件含8个.m脚本如mseq.m、goldseq.m、spread.m、despread.m等分别覆盖m序列与Gold序列生成、扩频、解扩及系统主流程6个.mat数据文件用于保存仿真输入输出便于复现实验另有Simulink模型相关文件.slxc、.xml及说明文件整体体积约97KB。代码结构清晰重点演示CDMA基带信号处理的关键环节可直接在Matlab环境中加载学习。该资源已有365人学习属轻量紧凑的示例包适合具备一定Matlab基础、希望自行调试代码并扩展功能的读者对照代码快速梳通CDMA系统从序列产生到解调输出的完整实现路径。1. CDMA仿真系统把码分多址原理变成能跑的模型做 CDMA 仿真系统最怕的不是原理看不懂而是原理写进 Simulink 之后模型跑出来的误码率怎么调都不对。这套基于 Matlab 和 Simulink 的 CDMA 仿真系统正好把最难啃的这块拆开了发射链路、PN 序列扩频、信道模型、接收端解扩判决、误码率统计一条龙拉通源码和数据都齐全。适合正在做通信课程设计、本科毕设或者刚入职想补 CDMA 物理层基础的人。我按一线工程师拆项目的习惯带你把每个模块的作用、参数、坑全部过一遍看完你能直接照着复现也能改成自己的实验参数跑出新结果。2. 系统拆分发射链路、信道模型与接收端的 Simulink 搭建2.1 从原理到模型动手前先想清楚的三个问题拿到这套资源别急着打开模型就跑。我每次接手一个通信仿真项目都会先花十分钟画一张数据流草图明确三个问题速率怎么对齐、同步怎么实现、误码率在哪一级统计。CDMA 的核心是扩频信息比特先乘以高速 PN 序列把窄带信号扩展成宽带信号接收端再用同步的 PN 序列做相关解扩把宽带信号重新收拢成窄带信号同时获得处理增益。这个过程在 Simulink 里落地第一个要解决的问题是采样率。信息速率假设是 Rb扩频因子是 SF那么码片速率 Rc SF × Rb。这意味着扩频模块前后的信号速率不一致Simulink 模型里必须用明确的样本时间Sample Time区分比特周期和码片周期否则数据流会乱。第二个问题是同步。实际 CDMA 接收机要做 PN 码捕获和跟踪但仿真模型里最常见做法是让发送端和接收端共用同一套时钟也就是用同一个 PN 序列生成器接收端直接“白嫖”发送端的码相位。这样做的好处是省掉同步模块把精力集中在扩频、解扩和误码率统计上代价是如果真想研究同步算法需要手动把接收端 PN 序列延迟几个码片观察相关峰变化。第三个问题是误码率统计口径。误码率应该统计判决之后的信息比特不是扩频后的码片。很多新手把 Error Rate Calculation 模块接在 BPSK 解调之后、解扩之前统计的是码片错误率这和系统真实误码率是两码事数值看起来会莫名其妙地高。2.2 发射链路PN 序列生成、扩频与 BPSK 调制发射链路在 Simulink 里主要由三块拼起来随机信息源、PN 序列生成器、扩频与调制模块。这套系统用的是 m 序列PN Sequence Generator 模块里需要配置多项式系数和寄存器长度。% 用脚本配置PN序列生成器的本原多项式 % 特征多项式: x^5 x^2 1对应多项式系数 [1 0 0 1 0 1] poly [1 0 0 1 0 1]; % 寄存器初始状态不能全0否则输出恒为0 init_state [1 0 0 0 0]; % 在Simulink模型里通过set_param配置PN模块 set_param(cdma_model/PN_Sequence, Polynomial, mat2str(poly)); set_param(cdma_model/PN_Sequence, InitialState, mat2str(init_state));这段代码的核心是 m 序列的两个硬约束多项式系数必须包含最高次项和常数项寄存器初态不能全零。常见做法是选 n5 或 n7 的 m 序列周期分别是 31 和 127。扩频因子 SF 通常取等于 m 序列周期比如 SF31对应的处理增益是 10×log10(31) ≈ 15 dB这个增益会直接体现在误码率曲线上——同样 Eb/N0 下扩频系统的误码率会比未扩频系统明显低。扩频操作在 Simulink 里的实现方式有两种。一种是用 Repetition重复模块把每个信息比特复制 SF 次再与 PN 序列逐码片相乘另一种是直接用 Gain 模块配合 PN 序列做乘法。这套系统看起来采用的是第二种思路PN 序列输出的 ±1 序列直接与 BPSK 映射后的数据相乘相当于每个比特和 SF 个码片按位相乘完成扩频。BPSK 调制用的是 BPSK Modulator Baseband 模块输入 0/1 比特映射为 1/-1。这里有一个容易忽略的细节BPSK 解调器输出是软值带幅度信息如果后续直接接解扩和判决必须搞清楚软值是正负号还带大小。如果是带大小的软值在 Rayleigh 信道下幅度起伏会直接影响判决需要做归一化。2.3 信道与接收AWGN、瑞利衰落与解扩判决信道模型是仿真系统里最值得替换的部分。系统自带两个信道配置AWGN 和瑞利平坦衰落。AWGN 信道用 AWGN Channel 模块实现参数主要是信噪比 SNR 或者 Eb/N0瑞利信道用 Multipath Rayleigh Fading Channel 模块需要配置多径延迟和增益。接收端的核心是解扩。解扩在数学上就是本地 PN 序列与接收信号的逐点相乘再在一个比特周期内累加最后做符号判决。Simulink 实现时一般用 Product 模块乘上同一个 PN 序列再用 Discrete-Time Integrator 做比特周期内的累加最后用 Sign 模块判为 ±1对应回 0/1 比特。% 配置AWGN信道的信噪比单位dB set_param(cdma_model/AWGN_Channel, SNR, 10); % 配置瑞利信道参数 % 多径延迟 [0 2e-6]对应两个径相差2微秒 set_param(cdma_model/Rayleigh_Channel, ... DelayVector, [0 2e-6], ... GainVector, [0 -3]);这里参数的含义分别是SNR 是接收端信号功率与噪声功率比值单位 dBRayleigh 信道的 DelayVector 是每条径相对第一条径的延迟秒数GainVector 是每条径的平均增益 dB。常见做法是让第一条径增益 0 dB后面径逐渐衰减。如果仿真瑞利信道却不加解扩误码率会非常难看——这正是 CDMA 多径分集要解决的问题。系统里自带一个简单的双径合并等增益合并把两条径的解扩结果直接相加再判决这已经能体现 RAKE 接收的基本思想。3. 在 Simulink 里跑通仿真脚本驱动、参数设置与误码率统计3.1 用脚本驱动仿真一套参数跑完整个实验这套系统的一个友好之处在于模型配置和仿真运行都可以用脚本控制。我一般会把所有可变参数集中在脚本开头模拟一个“参数面板”改参数不用进 Simulink 界面里翻模块。% 系统参数面板 SF 31; % 扩频因子等于m序列周期 Rb 1000; % 信息速率单位bps Rc SF * Rb; % 码片速率 numBits 5000; % 仿真总比特数建议至少5000 EbN0_dB 0:2:12; % Eb/N0扫描范围 % 计算采样率与仿真时长 samplePerChip 1; % 每个码片一个采样点 Ts 1 / (Rc * samplePerChip); % 采样时间 simTime numBits * SF * Ts; % 总仿真秒数 % 写入模型参数 set_param(cdma_model/PN_Sequence, SampleTime, num2str(Ts)); set_param(cdma_model/AWGN_Channel, SNR, num2str(EbN0_dB(1))); set_param(cdma_model, StopTime, num2str(simTime));这段脚本的作用是把系统参数和模型解耦想换码速率、换扩频因子、换仿真时长只改脚本第一段就行。注意 samplePerChip 这个变量它决定每个码片采几个点默认 1 就够了但如果后面接脉冲成型滤波器或频谱分析建议改成 4 或者 8否则频谱会有混叠。仿真停止时间必须按总比特数计算。当初我犯过一个低级错误固定 StopTime 为 1 秒结果信息速率一改统计的比特数就变了误码率曲线波动大得没法看。正确的做法是让仿真时长和比特数挂钩保证每次实验统计的符号数量一致。3.2 误码率统计从误码到曲线的完整链路Simulink 里统计误码率最直接的工具是 Error Rate Calculation 模块把它接在判决输出和原始信息比特之间模块输出三个数误码率、误码个数、对比的总比特数。用 To Workspace 模块把结果存到 Matlab 工作区仿真结束后统一处理。% 仿真运行结束后读取误码率结果 sim(cdma_model); ber BER_Out(1); % BER_Out是To Workspace输出的数组 numErrors BER_Out(2); % 对比理论BPSK误码率 berTheory 0.5 * erfc(sqrt(10^(EbN0_dB(1)/10))); % 打印对比结果 fprintf(Eb/N0 %.1f dB, 仿真BER %.5f, 理论BER %.5f\n, ... EbN0_dB(1), ber, berTheory);误码率统计最关键的是对齐比较的两路信号。Error Rate Calculation 模块的 Tx 端口和 Rx 端口必须延迟对齐否则统计出的误码率不是系统性能而是“信号没对齐造成的虚假误码”。系统中所有模块都带 Sample Time一旦前后速率不匹配模块输出会产生固有延迟需要在 Rx 路径上加 Delay 模块补偿。经验值是看仿真输出的误码个数如果误码率稳定在 0.5 附近说明对齐彻底错了信号完全错位。另一个细节是误码率统计的冷启动问题。仿真前几个比特因为模块初始状态未稳定误码率会偏高。常见做法是在脚本里设置一个“丢弃前 100 比特”的逻辑或者干脆把仿真比特数设大一些让前几百个比特的影响占比小到可以忽略。3.3 参数表扩频因子、信噪比与仿真时长的设定原则把系统里最常用的参数列成一张表方便对照修改参数项取值范围设定原则扩频因子 SF31、63、127越大处理增益越高但带宽占用越大仿真越慢信息速率 Rb100 ~ 10000 bps仿真模型里速率不影响结果但影响仿真速度每码片采样数1、4、8只观察误码率用 1做频谱或眼图建议 4 以上Eb/N0-2 ~ 14 dB根据期望误码率下限设定低于 -2 dB 误码率基本 0.5仿真总比特数5000 ~ 50000目标误码率越低需要统计的比特数越多多径延迟0 ~ 100 码片时长延迟超过码片周期后RAKE 合并效果才明显这套参数表几乎决定了仿真结果可信度。比如你目标误码率是 1e-4按统计学要求至少要统计 100 个错误比特才有意义也就是至少跑 100 万比特的仿真。虽然是 Simulink 模型100 万比特可能要跑十几分钟但低于这个量级跑出的 1e-4 完全不可信属于“样本太少纯属玄学”。4. 避坑与排查CDMA 仿真中最常见的五个翻车点4.1 误码率先高后低解扩时刻没对齐现象误码率曲线在低信噪比时接近 0.5信噪比升高后突然下落但始终比理论曲线差很多。原因接收端的本地 PN 序列和接收信号的码相位没有对齐解扩变成了“错位相乘”扩频增益完全没有发挥信号被自己扩频码的旁瓣干扰压制。解决检查接收端 PN 序列生成器的 InitialState 是否和发送端完全一致再检查 PN 序列模块的输出是否有额外延迟。在接收端 Path 上加一个 Delay 模块用仿真前 100 个比特做一次手动对齐让误码率先降到正常水平再调其他参数。这是 CDMA 仿真里最“玄学”的问题十有八九是这里。4.2 频谱显示全是毛刺采样率与码片速率关系搞错了现象用 Spectrum Scope 看发射信号频谱本该是主瓣清晰的扩频谱结果全是尖刺主瓣形状都看不出来。原因频谱分析仪需要一个合适的采样率如果每个码片只采 1 个点频谱会周期性延拓出现明显混叠。解决把 samplePerChip 改成 4 或 8也就是让码片速率乘上过采样倍数后再进入滤波器频谱就干净了。注意改了这个参数所有模块的采样时间都要跟着变不能只在频谱仪前面加一个插值模块。4.3 多用户干扰为零扩频码正交性没用上现象加了第二个用户的信号之后误码率竟然和单用户完全一样没有任何恶化。原因两个用户可能用的是完全相同的 PN 序列或者同一序列的不同相位相乘之后被自己解扩时完全消除了自干扰为零这本质上是同一个信号叠加了两次不是多用户干扰。解决给两个用户分配两条不同相位的 m 序列或者直接用 Gold 码。Gold 码互相关性可控是 CDMA 多用户仿真的标准选择。在 Simulink 中可以用两个 PN Sequence Generator 模块配置不同 InitialState 来生成不同相位的序列。4.4 仿真跑不动仿真时长与采样率组合失控现象模型仿真速度越来越慢跑几分钟只走了总时长的 1%。或者仿真直接卡死。原因仿真时长是按秒设置的而采样率是码片速率乘上过采样倍数。SF127、每码片采样 8 点、跑 5 万比特时总采样点数是 127×8×50000 5080 万点Simulink 每个采样点都要执行几十个模块计算量直接爆炸。解决先用 SF31、每码片采样 1 点跑通流程确认无误后再往大参数调。调大参数后先用 1000 比特验证速度速度能接受再跑全量。血泪经验能用 5000 比特验证的结论不要等 5 万比特跑完才看结果。4.5 误码率曲线不随信噪比下降判决阈值与幅度归一化问题现象Eb/N0 从 0 dB 调到 12 dB误码率始终徘徊在 0.1 左右曲线是一条平线。原因BPSK 解调器输出的是带幅度的软值经过瑞利信道后幅度起伏大直接用固定阈值 0 判决没有做幅度归一化低信噪比时噪声压过信号但高信噪比时深衰落瞬间的符号照样判错。解决在解扩累加之后用一个 Gain 模块或者 Divide 模块按信号平均幅度做归一化把软值归一化到 ±1 附近再判决。瑞利信道下不归一化误码率存在“地板效应”加再大信噪比也没用。5. 进阶把仿真系统变成研究工具的三件事5.1 批量跑误码率曲线脚本化参数扫描单点跑仿真只能验证一个 Eb/N0 下的误码率做研究必须画一整条 BER 曲线。常见做法是把仿真放进 for 循环每次只改 AWGN 信道的 SNR把误码率存成数组最后统一绘图。% 批量扫描Eb/N0绘制BER曲线 EbN0_dB 0:2:12; ber_sim zeros(size(EbN0_dB)); for k 1:length(EbN0_dB) set_param(cdma_model/AWGN_Channel, SNR, num2str(EbN0_dB(k))); sim(cdma_model); ber_sim(k) BER_Out(1); fprintf(Eb/N0 %.1f dB, BER %.5f\n, EbN0_dB(k), ber_sim(k)); end % 理论曲线对比 ber_theory 0.5 * erfc(sqrt(10.^(EbN0_dB/10))); semilogy(EbN0_dB, ber_sim, bo-, EbN0_dB, ber_theory, r*-); grid on; xlabel(Eb/N0 (dB)); ylabel(误码率); legend(仿真, 理论);这段脚本跑起来慢是正常的实时看打印输出注意趋势即可。理论上 AWGN 信道下 BPSK 的误码率曲线和 CDMA 仿真曲线应该基本重合差别在于 CDMA 的处理增益。如果画的曲线比理论低好几个数量级先别高兴检查是不是统计口径错把码片误码当成比特误码。5.2 可视化升级眼图、星座图与误码分布误码率只是一维指标真正定位问题还是要看眼图和星座图。把前端的信号接一个 Constellation Diagram 模块可以实时观察 BPSK 星座点的聚集情况。AWGN 信道下星座点应该是两团球簇瑞利信道下两团球簇会“糊”成一片这正是深衰落的直观表现。眼图更能说明问题。在 BPSK 调制之前、扩频之后加一个发射滤波器比如升余弦滤波器再在接收端加匹配滤波器眼图就会呈现出教科书里的形状。这套系统如果不加滤波眼图张不开但频谱利用率低加了滤波眼图好看了符号间干扰又冒出来。我一般用眼图观测接收端解扩前的信号质量用来判断信道损伤是噪声主导还是失真主导。5.3 从单用户到多用户Gold 码与干扰控制做完单用户系统下一步就是加用户。多用户 CDMA 的实质是同时同频传输靠码分区分用户。用 Gold 码替代 m 序列两个用户用不同生成多项式或不同相位的 Gold 序列接收端用对应码解扩时其他用户的信号表现为类似噪声的干扰。多用户系统的误码率公式不再只是热噪声还包含多用户干扰项。当用户数增加时BER 曲线会抬升这是因为干扰功率随用户数线性增长。我在仿真里通常让用户数从 2 加到 6观察 BER 曲线之间的距离这比任何理论公式都直观。从那以后我每次接手新的 CDMA 仿真模型都会强制自己先跑一遍 0~12 dB 的单用户 AWGN 曲线确认解扩同步、误码统计口径都正常之后再往模型里加信道、加多用户、加 RAKE 合并。模型每加一个复杂度就重复一遍这条基线流程出问题时能快速定位是哪一步引入的。希望这套拆解思路帮到你。本文还有配套的精品资源点击获取