简介一份面向通信工程学习者与科研初学者的MATLAB数字调制解调仿真源码包聚焦BPSK、QPSK、16QAM与64QAM四种主流数字调制技术代码围绕二进制数据映射、基带调制、加性高斯白噪声信道模拟、判决解调及误码率统计展开便于快速理解数字调制核心流程。压缩包共3个文件包含2个M脚本与1个txt说明文档整体约4KB结构精简适合直接运行与二次修改。已有352人学习浏览可作为通信原理课程设计、专业课实验或个人仿真练习的参考资料。通过阅读modulation.m可掌握不同阶数下的符号映射与星座点生成方式demodulation.m则演示了带噪信号下的判决门限设计在相同信噪比条件下运行仿真能直观对比四种调制方式的误码率差异深入体会频谱效率与抗噪声性能之间的权衡关系是数字通信入门与进阶的实用工具。1. 从mod-demod.rar说起一套能直接跑的调制解调框架该看什么如果你手头有一份名为 mod-demod.rar 的压缩包里面按目录堆着 BPSK、QPSK、16QAM、64QAM 的调制解调脚本那么它大概率是某个通信系统仿真或 SDR 实验的基础代码。这四种调制方式覆盖了从低阶到高阶的典型台阶BPSK 每符号 1 比特QPSK 每符号 2 比特16QAM 每符号 4 比特64QAM 达到 6 比特。它们共享同一套“映射成符号、送入信道、解映射成比特”的链路但星座点密度、判决边界和抗噪能力完全不同。这里先不纠结某个具体 rar 包里的文件怎么用重点说它背后那套最常见的从业方案星座图和格雷映射的理论基准MATLAB 里能直接跑通的最小仿真代码Simulink 与 AD9361 做硬件在环时的参数配合以及最后怎么用眼图、误码率曲线和自适应切换来验证链路没有问题。通信算法工程师、FPGA 验证人员和做软件无线电的学生都可以对应找到自己关心的段落。2. 调制解调的理论基准星座图、符号映射与判决边界2.1 先从BPSK和QPSK的相位域说起BPSK 把每个比特映射到两个相位之一通常取 0 和 π对应复平面上的 1 和 -1。解调时只需要看实部的正负虚部理论上不承载信息所以它对幅度失真和相位噪声的容忍度最高。QPSK 则是把两个比特合成一个符号映射到四个相位点上每符号携带 2 比特。从映射方式看QPSK 可以拆成 I、Q 两个正交分量各自的一个 BPSK因此判决也按 I、Q 独立进行。这种分解在实现定点和硬件逻辑时很有用两个支路可以并行处理判决边界就是 I0 和 Q0 两条坐标轴。星座图上的相邻点如何分配比特直接决定了误码性能。最常见的是格雷映射它保证任意相邻星座点之间的比特序列只差 1 位。比如 QPSK 的四个点按 00、01、11、10 的顺序沿相位圆周排列那么一个相位误判到相邻点只会引起 1 比特错误。如果不做格雷映射而用自然二进制 00、01、10、11 排列误判到邻近相位时可能错 2 比特。在高信噪比下这个差异会被放大误码率曲线会明显抬高。MATLAB 的 qammod、pskmod 默认支持 gray 选项这个参数在手动实现映射表时最容易漏掉。2.2 16QAM与64QAM的幅度相位联合映射16QAM 把 4 比特映射到 4×4 的复平面格点上I 路和 Q 路各取 -3、-1、1、3 四个电平。与 PSK 家族只靠相位区分不同QAM 的星座点同时依赖幅度和相位所以判决边界不再是简单的轴线而是一组间距为 2 的水平和垂直间线。64QAM 把点数扩展到 8×8电平扩展为 -7、-5、-3、-1、1、3、5、7每符号 6 比特。映射密度越大符号携带的信息越多但星座点之间的欧氏距离越小相同信噪比下误判概率更高。实际操作中64QAM 的判决器需要比 QPSK 多算 8 倍的距离或者做分级判决对定点实现时的数值精度也更敏感。实现 QAM 映射时能量归一化是一个绕不开的坑。如果直接用 -3、-1、1、3 去构造 16QAM 星座平均符号能量是 10而 64QAM 用 -7 到 7 的奇数电平平均符号能量是 42。若不归一化直接把两种调制放在同一个信道模型里比较误码率曲线会失真。常见做法是把构造好的星座点除以星座平均功率的开方使符号平均功率为 1。这样 BPSK 的 ±1、QPSK 的四点也都可以统一到同一功率基准上仿真时不同调制阶数之间的信噪比才可比。硬件上归一化因子会折算成基带信号乘以一个固定系数放在脉冲成型滤波器前面。2.3 调制阶数与信噪比门限的关系AWGN 信道下M-QAM 的理论误比特率由星座结构和格雷映射共同决定。工程上不需要每次都从 erfc 积分推导直接用 MATLAB 的 berawgn 函数给出参考曲线即可。下表列出四种调制方式在误码率达到 10^-4 时所需的典型每比特信噪比 Eb/N0数值基于相干解调和格雷映射实际系统通常还要留 0.5 到 1 dB 的实现余量。调制方式每符号比特数BER1e-4所需Eb/N0频谱效率bps/HzBPSK18.4 dB1QPSK28.4 dB216QAM412.2 dB464QAM616.5 dB6注意 BPSK 和 QPSK 的 Eb/N0 门限相同是基于每比特能量比较的必然结果。但从符号级看QPSK 的 Es/N0 比 BPSK 高 3 dB因为每个符号里有两个比特的能量。这张表的另一个用途是划定自适应调制切换的门限后面第 5 章会用到。需要提醒的是表格里的门限是理想相干解调下的值若链路中加入信道编码工作点会整体左移但相对关系不变。3. 在MATLAB里复现mod-demod的核心模块从波形生成到误码率统计3.1 最小可运行代码用qammod/qamdemod完成全链路先给出一段能直接复制到 MATLAB 里跑的最小完整代码。它实现了从随机比特生成、符号映射、加噪到解调判决的完整闭环通过修改 M 的值即可切换 BPSK、QPSK、16QAM 和 64QAM。% mod_demod_demo.m % 最小BPSK/QPSK/16QAM/64QAM基带误码率仿真 M 16; % 调制阶数2 / 4 / 16 / 64 k log2(M); % 每符号比特数 numBits 1e6; % 仿真总比特数 EbN0dB 12; % 每比特信噪比单位dB % 生成随机比特并按k的整数倍预先对齐 dataBits randi([0 1], numBits, 1); if mod(numBits, k) ~ 0 dataBits [dataBits; zeros(k - mod(numBits, k), 1)]; end numSymbols length(dataBits) / k; % 比特串转符号索引再用qammod做格雷映射 symbolIdx bi2de(reshape(dataBits, k, numSymbols)., left-msb); modSym qammod(symbolIdx, M, gray, InputType, integer); % 符号级信噪比 比特级信噪比 10*log10(k) EsN0dB EbN0dB 10*log10(k); % 加高斯白噪声awgn会先测量信号功率再按EsN0添加 noisySym awgn(modSym, EsN0dB, measured); % 解调并恢复比特流 rxIdx qamdemod(noisySym, M, gray, OutputType, integer); rxBits de2bi(rxIdx, k, left-msb).; rxBits rxBits(:); % 统计误比特数并计算误码率 numBitErr biterr(dataBits, rxBits); ber numBitErr / numel(dataBits); fprintf(M%d, Eb/N0%.1f dB, BER%.2e\n, M, EbN0dB, ber);逻辑说明代码先把比特流按 k 个一组转换成 0 到 M-1 的整数符号索引qammod在 gray 模式下会把符号索引映射到格雷码位置。加噪声用awgn这里必须先保证modSym的平均功率为 1因为awgn的 measured 选项是按信号实际功率来计算噪声方差的。EbN0dB 转 EsN0dB 必须加 10*log10(k)否则低阶调制和高阶调制之间的比较会错位。参数说明InputType,integer表明输入是整数索引而不是比特向量这个参数在较老版本的 MATLAB 中不能和 gray 一起使用需要先查一下你的版本是否支持。改 M 为 2、4、64 即可跑另外三种调制注意 64QAM 在高信噪比下需要更多比特数才能让误码率曲线平滑建议把 numBits 调大到 1e7 量级。这段代码没有考虑脉冲成型、定时同步和载波同步因此只适合验证调制解调本身的数学逻辑。3.2 手动实现符号映射与判决理解内部逻辑封装好的函数使用起来方便但不利于理解判决边界。以 QPSK 为例人为构造映射表和解调过程可以看清星座判决的本质。% qpsk_manual.m % 手动QPSK映射与最小距离判决 % 顺序映射仅用于演示判决逻辑实际工程请用格雷映射 mapTable [11j, -11j, -1-1j, 1-1j]; dataBits randi([0 1], 1e6, 1); symbolIdx bi2de(reshape(dataBits, 2, [])., left-msb); txSym mapTable(symbolIdx 1).; EsN0dB 10; rxSym awgn(txSym, EsN0dB, measured); % 每个接收符号与4个参考星座点求距离 dist abs(rxSym - mapTable); % 结果为 numSymbols x 4 [~, rxIdx] min(dist, [], 2); % 取距离最小的列索引 rxIdx rxIdx - 1; rxBits de2bi(rxIdx, 2, left-msb).; rxBits rxBits(:); ber sum(rxBits ~ dataBits) / length(dataBits);逻辑说明dist的每一行是当前接收符号到四个星座点的欧氏距离min(dist, [], 2)返回每行最小值所在的列也就是判决结果。这种穷举距离的方式在星座点少时很直观但 64QAM 有 64 个点每个符号要算 64 次距离仿真速度会明显下降。实际系统中会用 I/Q 分量的范围比较来代替全距离计算比如 16QAM 的 I 路大于 0 且大于 2 时判为 3这就是把复平面划分成判决域。参数说明rxSym - mapTable利用了 MATLAB 的隐式扩展rxSym 是 numSymbols×1 的列向量mapTable 是 1×4 的行向量两者相减得到 numSymbols×4 的距离矩阵。如果 MATLAB 版本较老需要改成abs(repmat(rxSym,1,4) - repmat(mapTable,numSymbols,1))。这里映射表用了自然的 00、01、10、11 顺序不是格雷映射所以误码率会比 qammod 带 gray 时差这里只用于理解原理。3.3 加性高斯白噪声信道下的仿真循环与参数设置单点误码率不足以验证链路需要在一组 Eb/N0 下循环仿真得到 BER 曲线。常见做法是把 3.1 节的逻辑封装成一个函数再写一个驱动脚本循环调用。这里给出驱动脚本的核心部分% ber_curve_driver.m M 4; % 可切换为2/16/64 EbN0dBSet 0:2:12; berSim zeros(size(EbN0dBSet)); for i 1:length(EbN0dBSet) berSim(i) run_ber_sim(M, EbN0dBSet(i), 2e5); end % 画仿真点 semilogy(EbN0dBSet, berSim, o-, LineWidth, 1.5); hold on; % 叠加理论曲线BPSK/QPSK用psk16QAM/64QAM用qam if M 2 || M 4 theo berawgn(EbN0dBSet, psk, M, nondiff); else theo berawgn(EbN0dBSet, qam, M); end semilogy(EbN0dBSet, theo, r-, LineWidth, 1); grid on; legend(仿真BER, 理论BER, Location, southwest);逻辑说明run_ber_sim是一个自定义函数内部实现与 3.1 节几乎一致只是把 fprintf 换成了返回结果。每个信噪比点至少执行 2×10^5 比特是为了让误码率在 1e-3 以下时有足够的统计样本64QAM 在低信噪比时误码率较高可以适当减少比特数但在高信噪比时误码率低必须增加比特数否则画出来的曲线会在底部出现台阶状抖动。参数说明berawgn的第四个输入 nondiff 表示差分编码关闭即采用绝对相位解调。如果仿真链路里用了差分编码理论曲线也要换成 diff二者不能混用。这段代码里 M4 时仿真链路没有差分编码所以用 nondiff 是对的。还有一个细节是semilogy的纵轴取对数误码率为 0 的点画不出来因此每个信噪比点要保证至少错 5 个比特以上否则该点会从图中消失。3.4 常见坑格雷映射、归一化因子、相位模糊第一个坑是格雷映射和自然映射混用。有的工程函数是自己写的映射表写成了自然二进制顺序却调用了带 gray 的 qammod 做解调导致星座点与比特对应关系错乱。判别方法很简单在无噪声条件下解调如果误码率不为 0映射顺序一定错了。另一个检查方法是用constellation(qammod(0:M-1,M,gray))打印星座图肉眼看相邻点是否只有 1 位变化。第二个坑是星座能量未归一化。用 ±1、±3 直接生成 16QAM 符号平均符号能量是 10但很多教材的误码率公式默认信号能量为 1。仿真时用mean(abs(txSym).^2)检查一下发射符号的功率。如果不是 1要么在生成星座后整体除以 sqrt(mean(abs(constellation).^2))要么在 awgn 中把信噪比设置为EsN0dB 10*log10(Es)其中 Es 是实际平均符号能量。更推荐前者因为后续接滤波器、均衡器时功率归一化的信号更好处理。第三个坑是相位模糊。AWGN 信道不会引入固定的相位旋转所以仿真链路很容易忽略载波同步问题。但一旦接入 AD9361 这类射频前端接收端本振与发射端本振存在未知相位差星座图会整体旋转。BPSK 会出现 0/π 模糊QPSK 会出现 π/2 整数倍的模糊。软件无线电工程里常用导频符号估计相位偏差或者在调制前做差分编码。仿真阶段可以暂时不管但心里要清楚这不是调制解调本身能解决的。提示如果只是验证调制解调算法建议在仿真初期把相位模糊排除在外先聚焦映射与判决逻辑到了硬件在环阶段再专门处理载波同步。4. 结合Simulink与AD9361从仿真到硬件在环的调制解调实现4.1 Simulink中调制解调模块的选型与连接Simulink 的 Communications Toolbox 里BPSK 和 QPSK 可以直接用 M-PSK Modulator Baseband 模块16QAM 和 64QAM 则用 Rectangular QAM Modulator Baseband 模块。两者名字里带“Baseband”表示工作于基带复数域输出是 IQ 符号流。调制模块的配置界面有几个参数直接影响后续链路行为M-ary number 设置为 16 或 64Constellation ordering 选 GrayNormalization method 建议选 Average Power 1这样 AWGN Channel 模块的 SNR 设置就是 Es/N0Phase offset 对 QPSK 可选 pi/4QAM 设为 0。模块连接上有一个容易漏掉的环节Rectangular QAM Modulator 的输入是整数符号而不是比特向量。因此 Bernoulli Binary Generator 输出的比特流需要先经过 Bit to Integer Converter按每 k 比特转成一个整数再进入调制器。解调链路对应地要用 Integer to Bit Converter 恢复比特。如果不加这个转换Simulink 会报维度不匹配错误而且这种错误在编译阶段才会暴露排查起来比较费时间。建议在模型里加一个 Display 模块直接观察调制器输出符号的实部虚部确认转换器工作正常。4.2 AD9361实现BPSK调制解调出数据的寄存器配置要点把 Simulink 仿真搬到 AD9361 板卡上是很多软件无线电项目的目标链路。AD9361 是一颗宽带 RF 收发器负责射频和模拟基带部分数字基带的调制解调仍然在 FPGA 或处理器里完成。因此“AD9361 实现 BPSK 调制解调出数据”这个问题的核心不是调制本身而是配置 AD9361 让它对 IQ 数据流的处理与你的调制参数匹配。常见做法是使用 ADI 官方的 AD9361 MATLAB/Simulink 支持包它提供了 Simulink 模型里可直接调用的 AD9361 Tx/Rx 模块。寄存器配置的关键点如下IQ 数据位宽与幅度缩放AD9361 的 ADC/DAC 是 12 位调制器输出的浮点符号要乘以一个缩放因子映射到整数范围建议留出约 30% 的峰均比余量也就是把信号峰值控制在满量程的 70% 左右防止削波失真又保证量化噪声不占主导。采样率配置方面AD9361 的采样率由 FPGA 时钟和芯片内部的 BBPLL 分频决定BPSK 符号速率如果是 1 Msps过采样倍数 4那么 IQ 采样率设为 4 Msps。配置时先查支持包里的采样率计算函数不要手动去改某个寄存器字段因为 BBPLL 和 FIR 滤波器级联关系很容易配错。增益控制也是 BPSK 调试阶段的重点。建议先用手动增益把 RX 增益固定在一个中等值比如 30 到 40 dB观察接收星座图。BPSK 的 1 和 -1 两个点如果幅度不一致或旋转再调整。链路稳定后再切换到 AGC让 AD9361 根据 RSSI 自动调节增益。注意 AD9361 的 AGC 有多种模式默认的 slow attack 模式适合语音类连续信号而 BPSK 突发信号要选择 fast attack否则 AGC 跟不上前导码的功率变化。4.3 采样率、滚降系数与匹配滤波器的配合无论是仿真还是硬件在环脉冲成型和匹配滤波都是 QPSK、16QAM、64QAM 链路不可跳过的一环。发送端用根升余弦滤波器做脉冲成型接收端再用一个相同的根升余弦滤波器做匹配滤波两级级联后等效为升余弦响应在采样点上满足奈奎斯特无码间干扰条件。这个结论在教科书里叫匹配滤波准则在工程里的直接含义是发射滤波器和接收滤波器必须使用相同的滚降系数和相同的过采样倍数否则星座点会散开。Simulink 里常用 Raised Cosine Transmit Filter 和 Raised Cosine Receive Filter 两个模块。发射滤波器的关键参数是滚降系数和过采样倍数。滚降系数一般在 0.2 到 0.35 之间取值越小带外辐射越少但脉冲长度越长受定时误差影响越大。工程里广播系统常用 0.2WLAN 系统常用 0.25 到 0.35具体看频谱约束。过采样倍数至少为 2推荐 4 或 8过采样倍数太小接收端采样相位不准确星座点发散太大FPGA 或处理器的时钟负担增加。参数推荐值说明滚降系数0.22 / 0.25 / 0.35频谱效率与定时敏感度的折衷过采样倍数4 或 8至少2推荐4起步滤波器符号跨度10 ~ 20截断误差与计算量的折衷调试时先固定过采样倍数为 4滚降系数 0.25观察接收端 QPSK 星座图的收敛情况。如果星座点发散得像一个环形先检查滤波器的采样率是否匹配如果星座点收拢但带明显的画圆轨迹说明定时同步没有锁定需要插入 Gardner 定时恢复模块。这些在 Simulink 的 Communications Toolbox 里都有现成模块但要注意它们默认参数与你的调制阶数是否一致比如 Gardner 模块内部的环路带宽在不同符号速率下需要重新调。5. 验证与调试技巧怎样判断你的mod-demod链路真的对了5.1 用眼图和星座图做第一轮体检当链路接好以后第一件事不是看误码率而是看接收符号的星座图。用 MATLAB 的scatterplot(rxSym)可以快速检查理想情况下 QPSK 应该只有四个点云16QAM 应该是十六个边界清晰的点簇。点云若以原点为中心旋转说明载波相位未修正点云半径分布比理论大很多说明定时同步或匹配滤波器配置有问题。眼图用eyediagram(rxSym, sps)看打开幅度最大的“眼睛”越方正说明码间干扰越小。硬件在环时把 AD9361 收到的 IQ 数据直接丢进这两个函数比盯着误码率数字更直观。5.2 误码率曲线与理论曲线对比的判定标准仿真链路的最终验证方式是把 BER 曲线叠加到理论曲线上比较。在误码率高于 1e-3 的区域仿真点与理论曲线的偏差应在 0.3 dB 以内低于 1e-3 后受比特数限制会有起伏可以接受在 0.5 dB 左右。如果偏差稳定在 1 dB 以上先查发射功率是否真的为 1再查噪声设置用的是 EsN0 还是 EbN0。还有一个常见的错误是 awgn 里用 measured但信号经过了没有归一化的星座映射导致实际信噪比与设定值不符。写成下面这样可以在仿真前主动做一次自检% check_power.m avgPower mean(abs(modSym).^2); if abs(avgPower - 1) 1e-6 modSym modSym / sqrt(avgPower); end这段代码的逻辑很简单先量出符号平均功率如果偏离 1 就归一化保证后续 awgn 和误码率比较都建立在同一功率基准上。注意这个检查要在脉冲成型之前做如果滤波器改变了信号功率需要在滤波器输出后再重新校准一次。5.3 一个实用的自适应调制切换策略把第 2.3 节的门限表用于自适应调制还需要加迟滞避免信道 SNR 在临界点附近抖动时反复切换。用估计 SNR 与两个阈值比较升阶阈值比理论门限高 1 dB降阶阈值比理论门限低 1 dB。% adaptive_mod_switch.m % estSNR由接收端估计单位dB if estSNR 17.5 M_next 64; elseif estSNR 13.2 M_next 16; elseif estSNR 9.4 M_next 4; elseif estSNR 5 M_next 2; else % SNR过低保持BPSK或进入重传 M_next 2; end参数说明阈值从第 2.3 节的 Eb/N0 门限上做了调整。16QAM 的理论门限 12.2 dB降阶阈值取 11.5 dB升阶阈值取 13.2 dB中间留了 1.7 dB 迟滞。64QAM 的门限 16.5 dB升到 17.5 dB 后才切掉到 16 dB 以下就退回 16QAM。实际信道还受衰落影响可以在 SNR 估计上再做时间平均比如对估计值做一阶低通滤波权重取 0.1 比 0.9这样切换频率就完全可控了。本文还有配套的精品资源点击获取