简介基于MATLAB的Simulink仿真工具完成的一篇通信系统课程设计报告面向通信工程、电子信息类专业学生以及需要开展通信系统仿真实验的读者为从系统建模、模块设计到性能测试的完整仿真过程提供了可直接参考的方案。报告围绕QPSK正交移相键控通信系统展开系统介绍了信源模块、信源编码包括抽样、量化与DPCM差分脉冲编码调制、QPSK调制与解调、加性高斯白噪声信道、误码率测试与功率谱分析等核心环节并包含系统组成结构框图、各模块原理说明、方案论证、仿真结果分析与设计总结等完整章节。压缩包共1个PDF文件大小仅1.18MB内容完整便于随时查阅。报告章节结构清晰、步骤完整既可作为通信系统课程设计报告写作模板也能帮助初学者掌握Simulink环境下数字通信系统建模、调试与性能验证方法对完成课程设计或期末仿真作业具有直接参考价值。已有227人学习下载适合通信工程及相关专业学生参考使用。1. 通信系统设计不能只靠公式Simulink 把链路变成看得见的模型做通信系统的人都有这种经验算法推演到公式层已经很完整但码本映射、成形滤波、信道噪声叠在一起理论曲线和实际链路差几个 dB 是常有的事。用 MATLAB 写脚本做数值验证可以逐级跑通可一旦想观察信号波形、眼图、频谱或者换一种调制方式重跑整条链路脚本就得拆了重写。Simulink 的优势在于把发射机、信道、接收机拖成框图信号在哪一级变成什么样一目了然参数改动不用动核心逻辑。这篇文章从理论底子讲到工程落地的参数设置、模型组织和验证方法适合正在做通信系统仿真验证的工程师和学生。2. 从信源到信宿Simulink 通信仿真的建模基础2.1 为什么先建基带等效模型射频仿真的代价高得吓人。如果载频是 2.4 GHz符号速率 1 Mbps直接仿真一个完整的载波周期需要步长小于 0.4 ns仿真 1 秒数据就要跑几十亿步普通 PC 根本扛不住。通信系统仿真的标准做法是搬到基带把调制信号分解为 I/Q 两路复基带信号把载波从数学上约掉仿真步长只需覆盖符号周期通常是符号周期的几分之一。Simulink 的 Communications Toolbox 里很多调制模块默认输出的就是复数基带信号比如 BPSK Modulator Baseband、Rectangular QAM Modulator Baseband它们的输出端口虽然看着是标量或向量实际上内部表示的是一个复基带序列。这样仿真时间和内存占用都从射频量级降到基带量级数字信号处理部分的时钟关系也和真实硬件更接近。2.2 信源、信道与信宿的模块化映射Simulink 通信仿真链路说到底就是三个环节信源产生比特信道叠加损伤信宿恢复信息。信源侧常用 Bernoulli Binary Generator 生成随机 0/1 序列或者用 Random Integer Generator 生成 M 进制符号信道侧用 AWGN Channel 对应高斯白噪声用 Multipath Rayleigh Fading Channel 对应多径衰落信宿侧先做匹配滤波、同步、解调再用 Error Rate Calculation 模块对比收发比特得到实时误码率。把这几个模块按信号流向接起来再连一个 Display 或 Scope一条可观察的链路就跑起来了。表 2-1 列出了最常见的模块及其用途设计新链路时可以直接根据这张表找模块。表 2-1 通信系统设计时常用的 Simulink 模块模块典型用途所属模块库Bernoulli Binary Generator生成等概率 0/1 比特流Communications Toolbox / SourcesRandom Integer Generator生成均匀分布的 M 进制符号Communications Toolbox / SourcesBPSK Modulator Baseband基带 BPSK 调制Communications Toolbox / ModulationAWGN Channel加性高斯白噪声信道Communications Toolbox / ChannelsMultipath Rayleigh Fading Channel多径衰落信道模型Communications Toolbox / ChannelsError Rate Calculation接收/发送比特对比输出 BERCommunications Toolbox / Comm Sinks这里有个细节调制器既接受比特序列也接受整数符号模式由输入端口的类型决定。在 Comm Toolbox 中调制模块的输入可以配置为比特或整数这会影响后面的符号映射。如果解调端和调制端设置不一致误码率会直接飙升哪怕信道是理想的。2.3 仿真粒度符号级还是比特级新手最常踩的坑是把符号和比特混在一起。BPSK 一个符号一个比特没有歧义到了 QPSK一个符号 2 个比特在 Simulink 里调制器输入是比特流输出是符号流噪声模块作用在符号域解调器输出又是比特流。若用脚本计算理论曲线必须清楚你说的是符号差错率SER还是比特差错率BER。对于 QPSK若差分编码或 Gray 映射不同BER 和 SER 的关系也不同。所以我在每个模型里都先写一段初始化脚本把所有调制阶数、符号速率、采样率、每符号采样数定义成 MATLAB 工作区变量然后在 Simulink 模块参数里引用这些变量名而不是直接填数字。% 通信链路初始化脚本各模块参数引用这些变量 M 4; % 调制阶数QPSK k log2(M); % 每符号比特数 2 Rs 1e6; % 符号速率 1 Msps sps 20; % 每符号采样数用于离散仿真 Fs Rs * sps; % 采样率 20 Msps EbNo 8; % 目标比特信噪比单位 dB Nbits 1e6; % 仿真总比特数这段脚本没什么魔法但把参数集中放有两个好处一是改调制方式时只需要改 M下面所有换算自动生效二是后面做参数扫描时外层循环可以直接修改EbNo、Nbits不必动模型。sps这个变量特别重要它决定了数据流在 Simulink 中的采样网格密度。如果sps设成 1就是符号率采样适合理想匹配滤波下的蒙特卡洛仿真如果要做成形滤波或观察眼图sps必须大于等于 2工程上通常取 8 到 32。采样率过低会产生频谱混叠过高会增加仿真时间所以取一个适中值是模型设计的第一道平衡。3. 搭一个能跑的最小通信系统BPSK AWGN 仿真链路3.1 最小链路模块选择现在搭一条最经典的链路BPSK 调制经过 AWGN 信道BPSK 解调统计误码率。这个链路虽然简单但包含了通信仿真所有核心环节后续扩展成 QPSK、8PSK 甚至 OFDM 都只是替换调制器和增加同步模块的问题。在 Simulink 中新建一个模型从模块库中拖出以下模块Bernoulli Binary Generator信源BPSK Modulator Baseband调制AWGN Channel信道BPSK Demodulator Baseband解调Error Rate Calculation误码率统计Display实时显示连接顺序为信源 → 调制器 → AWGN 信道 → 解调器 → 误码率统计。同时误码率统计模块还需要第二路输入即原始的发送比特所以从信源直接引一条线到 Error Rate Calculation 的 Tx 端口。注意这里存在一个时序问题发送比特走过调制、信道、解调后到达误码率模块的时间与直接从信源到达误码率模块的时间差必须用延迟模块对齐否则误码率模块会拿不同时刻的比特做比较。常见做法是在信源直通路径上增加一个 Delay 模块延迟值设为调制器和解调器引入的总延迟。3.2 关键参数与典型取值参数设置不是在模块对话框里瞎填的每个参数都和前面的初始化脚本绑定。表 3-1 是这一条链路里几个关键参数的设置建议。表 3-1 BPSK AWGN 链路关键参数模块参数设置值说明Bernoulli Binary Generator - Probability of a zero0.5等概率比特流BER 统计有意义Bernoulli Binary Generator - Sample time1/Rs每个比特占一个符号周期若做成形滤波则改为 1/FsBPSK Modulator Baseband - Output typeBit输出符号流每符号一个样本AWGN Channel - ModeEb/N0直接给出比特信噪比省去手动换算AWGN Channel - Eb/N0 (dB)EbNo引用工作区变量便于扫描Error Rate Calculation - Output dataPort将 BER 输出到端口供 MATLAB 读取Error Rate Calculation - Stop simulationNbits仿真多少比特后自动停止注意 AWGN Channel 的 Mode 下拉框里有Eb/N0、Es/N0、SNR三种。Eb/N0和Es/N0的区别是后者基于符号能量对于格雷映射的 QPSKEs/N0 Eb/N0 10*log10(k)。在 BPSK 中 k1两者一样。如果后面换 QPSK但 AWGN 模块仍选SNR模式还要额外考虑过采样因子sps的影响容易算错。所以我的习惯是能选Eb/N0就不选SNR。3.3 用脚本驱动仿真和读取结果模型搭好后直接点 Run 也能看到 Display 上的实时 BER但要想批量跑不同EbNo点还是得用脚本。下面这段脚本演示了如何从 MATLAB 命令行启动仿真并读取结果。% 设置仿真停止时间让误码率模块内部统计完再停 set_param(bpsk_awgn/Error Rate Calculation, StopSim, on); set_param(bpsk_awgn/Error Rate Calculation, StopSimTime, Nbits/Rs); % 运行仿真 open_system(bpsk_awgn); simOut sim(bpsk_awgn); % 从仿真输出中读取误码率统计 ber simOut.get(BER); % 假设在模型中定义了 To Workspace 模块输出 BER fprintf(模拟 Eb/N0 %g dB 时 BER %g\n, EbNo, ber(1));这里sim函数的第一个参数是模型名StopTime实际不用显式设置因为 Error Rate Calculation 模块内部已经配置了停止条件。读取结果前最好在模型中加一个 To Workspace 模块连接到 Error Rate Calculation 的输出端口变量名填BER否则要从 scope 历史里想办法导数据很麻烦。上面脚本里的simOut.get(BER)只是其中一种方式如果 To Workspace 保存格式默认是timeseries那么ber(1)提取的是最后一次统计的误码率值。Nbits/Rs这个表达式算的是仿真物理时间同时要确认模型中没有其他连续模块造成额外的时间偏移。4. 把仿真跑稳采样率、Eb/N0 与发散排查4.1 采样率与步长该匹配到多细通信系统仿真经常出现“跑起来不报错但结果明显不对”的情况。最常见的根源是采样率和步长的关系没有理顺。如果整个链路都是离散模块Simulink 会用模块指定的采样时间自动计算每个信号的离散速率。比如前面参数里把 Bernoulli Binary Generator 的采样周期设为1/Rs那后续每个符号一个样本整个链路的离散积分步长就是符号周期。这种设置适合做纯蒙特卡洛 BER 仿真因为不需要恢复波形只需要符号判决。但如果你后面要从接收信号中画眼图、星座图或者要做匹配滤波就需要过采样。过采样时把信源采样时间改成1/Fs同时调制器输出端每符号会产生sps个样本这时要检查 AWGN Channel 的输入采样率是否匹配。通常 AWGN 模块会自动继承输入端的采样率不会报错但它内部对噪声功率的计算依赖sps如果sps不对噪声方差就错最终 BER 曲线会整体偏移。所以调模型时第一件事就是把每个模块的采样时间显示出来。在 Simulink 模型界面用快捷键CtrlJ可以切换显示采样时间颜色红色表示固定步长离散任务颜色不同代表不同的采样速率。如果看到一条链路里相邻两个模块的采样时间颜色不一样又没有插入 interpolate 模块那多半就是采样率不匹配。4.2 Eb/N0 与噪声功率的换算细节理论公式里加性高斯白噪声的单边功率谱密度是 N0信号每符号能量是 Es每比特能量是 Eb Es / k。用伏特表示时噪声的标准差和信号幅度有关具体值取决于采样率与符号率的比值sps。在 Simulink 的 AWGN Channel 模块里如果选Eb/N0模式模块内部会利用输入信号的功率自动计算所需噪声方差理论上是省心的。但要注意模块的输入信号功率计算的是实际样本序列的平均功率。如果信号序列中含有直流分量或没有做归一化那么计算出的 SNR 会和理论设定不一致。结果就是仿真 BER 曲线整体向低信噪比方向偏移看起来比理论还好这显然是错的。我遇到这种问题时先检查信源输出的信号功率是否归一化比如 BPSK 调制器输出的 ±1 序列理想功率是 1若把映射改成 0/1功率就是 0.5噪声功率计算就会出现 3 dB 偏差。表 4-1 常见调制阶数下 k 值与 Es/N0 和 Eb/N0 的关系调制方式Mk log2(M)Es/N0 - Eb/N0 (dB)BPSK210QPSK423.018PSK834.7716QAM1646.02需要手动换算时直接在 MATLAB 里写EsNo EbNo 10*log10(k)即可。4.3 仿真发散先分清代数环还是采样冲突通信仿真里说的“发散”往往不是数值无穷大而是波形乱跳、BER 为零或接近 1。最常见的罪魁祸首是代数环和采样时间冲突。代数环多出现在把统计量的输出又回接到信源或调制使能端的闭环模型里。比如你想用实时误码率控制信源是否继续发送直接把 Error Rate Calculation 的输出接到 Bernoulli Binary Generator 的使能端口形成一个环。由于模块的输出依赖输入的当前值Simulink 在每个步长里都必须解一个代数方程解不出来就会报 “Algebraic loop detected” 或者无限迭代。排查方法很简单把这类回接路径断开改用条件子系统或者在 MATLAB 工作区里用逻辑控制仿真进程。采样时间冲突则更隐蔽。比如一个离散模块的采样时间是1/Rs后面接了一个连续求解器模块Simulink 可能会报“Invalid setting for input port”。如果不报错也会造成输出信号的采样混乱。这时用 Model Advisor 的 “Check for modeling errors” 或直接看 Sample Time Legend很快能定位。总的原则是纯离散链路尽量全部使用离散模块求解器设为固定步长比如discrete求解器。下面是一个典型的排查命令用于检查模型的离散求解器状态set_param(bpsk_awgn, Solver, FixedStepDiscrete); set_param(bpsk_awgn, FixedStep, 1/Fs);把求解器固定成离散以后Simulink 会忽略连续模块并提示警告这恰恰能逼你发现哪些模块是不必要连续采样的。很多仿真发散其实只是初始步长太小或太大换成离散求解器后问题自然消失。5. 从单模型到交付模型引用、外部模式与代码生成5.1 用模型引用拆分发射机、信道和接收机一个复杂的通信系统模型如果全放在一张画布上几百个模块连起来后期维护就是灾难。我一般把发射机、信道部分、接收机分别建成独立模型再用 Model Reference 在顶层模型里引用。这样每个子模型可以单独打开、单独仿真、单独做单元测试。在 MATLAB 中用Model Reference的方式加载子模型代码层面可以用set_param指定。% 顶层模型 top_model 中引用子模型 tx_model set_param(top_model/Transmitter, ModelName, tx_model); set_param(top_model/Transmitter, SimulationMode, Normal);这里Transmitter是顶层模型中 Model 模块的实例名tx_model是子模型文件名。需要注意子模型的输入输出端口名称必须和顶层信号线名称匹配否则会报端口找不到。引用模式下子模型的参数仍然可以在顶层用set_param调节吗实际上可以的通过Simulink.BlockDiagram的接口但更常见的做法是在子模型里把需要调用的参数定义为Simulink.Parameter或直接引用工作区变量这样工作量最小。5.2 外部模式一边跑仿真一边改参数做联合调试时比如接收机算法和信道模型同时迭代你不想每次改一个系数就重新从头仿真。Simulink 的外部模式External Mode允许你在模型运行过程中通过宿主的界面修改参数参数会实时下传到运行中的模型。这个模式原本是给实时硬件用的但对普通软件仿真也能用。打开步骤在 Simulink 窗口中菜单 Simulation Mode External然后点击运行。要注意的是外部模式要求信号接口支持数据交换模型会额外注入一组通信通道仿真速度会有一些下降但好处是可以动态调节例如环路增益、滤波器系数而不中断数据流。如果你想在外部模式下也改变采样率那就没戏了采样结构是编译期固定的只能在仿真前改。5.3 从 Simulink 模型生成 C 代码链路在 Simulink 里验证完最终要移植到嵌入式平台或作为算法交付给 C 工程不能每次都在 MATLAB 环境里跑。Simulink Coder 可以把模型生成独立 C 代码。对于通信系统我常用配置是解决器固定步长、目标选ert.tlcEmbedded Real-Time Target、代码包装为可重入函数。关键设置如下表。表 5-1 通信系统模型生成代码的关键配置配置项推荐值说明Solver - TypeFixed-step代码生成不支持变步长Solver - Solverdiscrete离散求解器生成代码最简单Code Generation - System target fileert.tlc嵌入式实时目标Code Generation - LanguageC可选 C/CCode Generation - Pack dataCharacter arrays方便嵌入网络字节流生成命令一行搞定slbuild(bpsk_awgn);slbuild会在模型所在目录生成bpsk_awgn_ert_rtw文件夹里面是完整的 C 源码和 Makefile。生成代码前一定要用 Model Advisor 的 Embedded Coder 检查项过一遍比如检查是否有连续时间模块、是否有动态内存分配。通信链路里的AWGN Channel模块在代码生成时可能需要替换成手写的随机数生成器因为默认的 TLC 不一定支持所有通信库模块。如果你发现某个模块不支持生成代码常见做法是将那个模块替换为 MATLAB Function 模块自己写等价的系统函数保证算法可移植。6. 用理论曲线验证仿真误码率结果6.1 用理论曲线做基准在批量跑完 BER 后仿真数据本身不会告诉你对错只有拿理论曲线对比才能发现系统性偏差。对 BPSK理论误码率是Q(sqrt(2*Eb/N0))在 MATLAB 中可以用berawgn直接拿到不同调制方式的理论值。下面这段代码是最常用的对比脚本EbNo_vec 0:1:10; ber_theory berawgn(EbNo_vec, psk, 2, nondiff); % 假设仿真结果存成了 ber_sim 向量长度与 EbNo_vec 一致 figure; semilogy(EbNo_vec, ber_theory, b-, LineWidth, 1.5); hold on; semilogy(EbNo_vec, ber_sim, ro, MarkerFaceColor, r); grid on; xlabel(Eb/N0 (dB)); ylabel(Bit Error Rate); legend(理论, 仿真);如果仿真点出现在理论曲线上方 0.5 dB 以内基本可以认为模型是正确的如果整体偏到下方通常说明噪声功率被低估检查 AWGN 模块的载波功率归一化如果偏到上方则可能是同步或延迟对齐出了问题。另一种验证思路是用bertool图形界面它内置多种调制理论和仿真配置能自动帮你画出对比曲线。6.2 蒙特卡洛次数怎么够验证时不要只画最终 BER还要监控仿真过程中错误比特数。蒙特卡洛仿真要得到稳定估计至少需要观察到 100 个错误比特如果目标 BER 是 1e-4那么仿真比特数就要达到 1e6 量级。一个取巧的做法是在 Error Rate Calculation 模块输出中还带有符号数、误比特数直接判断误比特数是否大于某个阈值再停止仿真避免白白跑满时长。通常把最小错误次数设为 50这样 BER 的相对置信区间大概在 ±30% 以内要更精确就设成 500。这个阈值写在初始化脚本里再通过set_param注入到模块参数中整个流程就能自动完成。最后一个小技巧如果你的链路包含有色散、多径或者非线性器件单靠 AWGN 理论曲线做基准是不够的这时可以用两个不同实现方式互相验证——比如同一个信道用 Simulink 模块搭另一个用 System object 在 MATLAB 里搭两者跑同样比特流和信道种子比较输出的原始软比特序列是否一致。只要软比特误差在浮点精度范围内就说明 Simulink 模型和参考实现完全等价。这个方法比单纯看 BER 更灵敏任何时序错位、极性反转都会直接暴露出来。本文还有配套的精品资源点击获取