2FSK调制解调MATLAB仿真：从课设代码到误码率性能分析

简介这份资源是面向通信工程、电子信息类专业学生的课程设计参考资料围绕2FSK调制与解调系统的设计与仿真展开适合正在学习通信系统原理、需要完成MATLAB课程设计或想加深数字调制理解的学习者。资源包内共1个doc文档约579KB内容涵盖设计任务与要求、方案论证、2FSK调制与解调原理、源程序代码及仿真结果分析完整呈现了从信源产生、键控法调制、加性高斯白噪声信道到相干解调与误码率性能评估的全过程。读者可从中获取课程设计的完整框架与实现思路包括二进制随机比特流生成、发送与接收滤波器设置、相干解调器结构以及SNR与误码率曲线的绘制方法并借助MATLAB7.0源程序复现仿真波形对照分析系统性能。目前已有853人学习下载适合作为通信原理课程设计的实操范本与排错参考。1. 从一份课设代码说起2FSK 调制解调到底在仿真什么很多人第一次接触数字通信仿真都是从 2FSK 开始的。原因很直接它足够简单又足够完整——信源、调制、信道、解调、误码率统计一条链路全走一遍。这份《2FSK 调制与解调系统设计与仿真》就是典型的通信原理课程设计产物基于 MATLAB 7.0 实现核心目标是用代码把二进制频移键控的调制、加噪、相干/非相干解调、误码率对比全部跑出来。它适合两类人正在做通信课设、需要一份能直接跑通的参考实现以及想用 MATLAB 快速验证调制解调算法、不想从零搭链路的开发者。整份代码不到 200 行但把 dmod/ddemod、awgn、symerr、rcosine 这些通信工具箱函数串成了一条完整链路拿来改参数、换调制方式都很方便。2. 系统架构与模块拆解一条 2FSK 链路是怎么串起来的2.1 整体链路与模块划分这份设计的结构非常清晰按信号流向可以拆成六个模块信源产生、2FSK 调制、高斯信道加噪、解调相干/非相干、误码率统计、滤波与眼图观察。每个模块在代码里对应一段独立逻辑模块之间通过信号变量传递没有复杂的类封装全是脚本式顺序执行。这种写法在课设里很常见好处是调试直观坏处是变量全在 workspace 里改一处容易影响后面。从原理层面看2FSK 信号可以理解为两个不同载频的 ASK 信号叠加。本设计采用开关法产生二进制“1”对应载频 f1“0”对应载频 f2两个频率源独立所以相邻码元之间相位不一定连续。这一点和调频法产生的连续相位 2FSK 有本质区别在后续解调时会影响相干载波的相位对齐。解调部分设计了两套方案做对比相干解调和非相干解调。相干解调需要本地恢复与发送端同频同相的载波分别与 f1、f2 带通滤波后的信号相乘再低通滤波、抽样判决非相干解调则用包络检波不依赖相位信息。代码里通过ddemod的不同参数来切换这两种方式后面会详细讲。2.2 关键参数与初始化任何仿真链路参数设置都是第一步也是最容易翻车的地方。这份代码的参数集中在开头几行Fc 10; % 载频 Fs 40; % 系统采样频率 Fd 1; % 码速率 N Fs/Fd; % 每个码元的采样点数 df 10; % 频率间隔 numSymb 25; % 仿真信息码元个数 M 2; % 进制数2FSK 固定为 2 SNRpBit 60; % 每比特信噪比(dB) SNR SNRpBit/log2(M); % 转换为符号信噪比 seed [12345 54321]; % 随机种子保证可复现 numPlot 15; % 绘图显示的码元数这里有几个参数需要重点理解。Fs/Fd 40意味着每个码元采 40 个点采样率足够高才能画出平滑的调制波形。df 10是 f1 和 f2 的频率间隔它决定了两个载频差多少直接影响解调时带通滤波器能不能分开。SNRpBit 60这个值偏大实际跑出来的误码率会很低课设里为了看到“干净”的解调结果通常会这么设但做性能曲线时一般会扫一个 SNR 范围。seed数组控制随机数状态randsrc和randn分别用它来固定信源和噪声这样每次运行结果一致方便对比调试。这个习惯在实际工程里也很重要不然每次跑出来的误码率都在跳根本没法判断改动是否有效。2.3 信源产生与调制实现信源部分用randsrc产生 0 到 M-1 的随机整数序列M2 时就是二进制比特流x randsrc(numSymb, 1, [0:M-1]); % 产生 25 个二进制随机码 figure(1) stem([0:numPlot-1], x(1:numPlot), bx); title(二进制随机序列); xlabel(Time); ylabel(Amplitude);randsrc(numSymb, 1, [0:M-1])生成一个 25 行 1 列的列向量每个元素等概率取 0 或 1。stem画杆图只显示前 15 个码元避免图太密看不清。这里用单极性矩形波表示基带信号符合课设要求。调制用通信工具箱的dmod函数一行搞定y dmod(x, Fc, Fd, Fs, fsk, M, df); numModPlot numPlot * Fs; t [0:numModPlot-1]./Fs; figure(2) plot(t, y(1:length(t)), b-); axis([min(t) max(t) -1.5 1.5]); title(调制后的信号);dmod的参数依次是输入符号、载频、码速率、采样率、调制方式、进制数、频率间隔。fsk指定频移键控M2就是 2FSKdf10控制两个频率的间隔。返回的y是已调信号长度是numSymb * Fs。画图时只取前numPlot * Fs个点对应 15 个码元的波形这样能看到频率切换的细节。2.4 高斯信道加噪信道部分用awgn函数加入高斯白噪声randn(state, seed(2)); y awgn(y, SNR - 10*log10(0.5) - 10*log10(N), measured, [], dB); figure(3) plot(t, y(1:length(t)), b-); axis([min(t) max(t) -1.5 1.5]); title(加入高斯白噪声后的已调信号);这里有个细节值得注意awgn的第二个参数是信噪比但代码里做了一串修正——SNR - 10*log10(0.5) - 10*log10(N)。原因是awgn默认按信号功率和噪声功率比来算而通信原理里的 SNR 通常定义在符号能量和噪声功率谱密度上两者之间差了一个采样点数 N 和调制因子。这个修正是课设里常见的“血泪经验”不修正的话实际信噪比会和预期差很多误码率曲线对不上理论值。randn(state, seed(2))固定噪声生成器的状态保证每次加噪结果一致。measured表示让awgn先测量信号功率再按指定 SNR 加噪比默认的假设功率更准确。3. 解调与性能分析相干和非相干到底差在哪3.1 相干解调实现与眼图相干解调用ddemod的fsk/eye模式figure(4) z1 ddemod(y, Fc, Fd, Fs, fsk/eye, M, df); title(相干解调后的信号的眼图); figure(5) stem([0:numPlot-1], x(1:numPlot), bx); hold on; stem([0:numPlot-1], z1(1:numPlot), ro); hold off; axis([0 numPlot -0.5 1.5]); title(相干解调后的信号原序列比较); legend(原输入二进制随机序列, 相干解调后的信号);ddemod的参数和dmod对应fsk/eye表示解调并绘制眼图。返回的z1是解调后的符号序列。图 5 把原始序列x和解调序列z1画在一起对比蓝色杆是原信号红色圈是解调结果重合就说明解调正确。眼图是判断解调质量的重要工具。眼张开越大说明码间干扰越小、判决余量越大。如果眼图闭合说明滤波或采样时刻有问题解调误码率会很高。课设里通常要求观察眼图并分析这一步不能省。3.2 非相干解调与对比非相干解调只需要在ddemod里加noncoh标志figure(6) z2 ddemod(y, Fc, Fd, Fs, fsk/eye/noncoh, M, df); title(非相干解调后的信号的眼图); figure(7) stem([0:numPlot-1], x(1:numPlot), bx); hold on; stem([0:numPlot-1], z2(1:numPlot), ro); hold off; axis([0 numPlot -0.5 1.5]); title(非相干解调后的信号); legend(原输入二进制随机序列, 非相干解调后的信号);fsk/eye/noncoh告诉ddemod用非相干方式解调内部走包络检波路径。非相干解调不需要恢复载波相位实现简单但在相同信噪比下误码率比相干解调高。代码里把两种解调结果都画出来方便直观对比。3.3 误码率统计与 SNR-误码率曲线误码率统计用symerr函数[errorSym, ratioSym] symerr(x, z1); figure(8) simbasebandex([0:1:5]); title(相干解调后误码率统计); [errorSym, ratioSym] symerr(x, z2); figure(9) simbasebandex([0:1:5]); title(非相干解调后误码率统计);symerr(x, z1)比较原始序列和解调序列返回错误符号数errorSym和误码率ratioSym。simbasebandex是通信工具箱自带的示例函数用来画基带仿真的误码率曲线这里直接调用它来展示不同 SNR 下的性能。从理论上看2FSK 相干解调的误码率在信噪比很大时近似为0.5 * erfc(sqrt(SNR/2))非相干解调近似为0.5 * exp(-SNR/2)。相同信噪比下相干解调的误码率更低性能更好。但相干解调需要载波同步实现复杂度更高。课设里把两种方式都跑一遍就是为了验证这个结论。3.4 升余弦滤波与带外噪声抑制最后一段代码用升余弦滤波器对加噪信号做滤波观察滤波前后的波形对比Delay 3; R 0.5; PropD 0; [yf, tf] rcosine(Fd, Fs, fir, R, Delay); [yo2, to2] rcosflt(y, Fd, Fs, filter, yf); t [0:numModPlot-1]./Fs; figure(10) plot(t, y(1:length(t)), r-); hold on; plot(to2, yo2, b-); hold off; axis([0 30 -1.5 1.5]); xlabel(Time); ylabel(Amplitude); legend(加入高斯白噪声后的已调信号, 经过升余弦滤波器后的已调信号); title(升余弦滤波前后波形比较); eyediagram(yo2, N); title(加入高斯白噪声后的已调信号的眼图);rcosine生成升余弦滤波器的系数R 0.5是滚降系数Delay 3是滤波器延迟。rcosflt对信号做滤波滤除带外噪声。图 10 里红色是滤波前蓝色是滤波后可以看到滤波后波形更平滑毛刺明显减少。最后用eyediagram画滤波后的眼图眼张开程度比滤波前更好。注意升余弦滤波器的延迟会导致输出信号比输入滞后画图对比时时间轴要对齐否则会误以为波形变形。4. 避坑与常见问题课设里最容易翻车的几个点4.1 振荡器频率必须是 pi 的倍数现象调制模块跑出来的波形不是完整的正弦波频率切换处出现断裂或畸变。原因MATLAB 里正弦波用sin(2*pi*f*t)生成如果 f 不是整数或不是 pi 的倍数关系在一个时间单位内可能采不到完整周期导致波形不完整。解决把 f1、f2 设置为与采样率成整数倍关系的值或者直接用dmod内部生成的载波不要自己手写正弦。课设里用dmod就避开了这个问题。4.2 带通滤波器上下限频率设置不当现象Simulink 仿真中带通滤波器输出失真严重解调误码率异常高。原因带通滤波器的中心频率应该对准 f1 或 f2带宽要能覆盖调制信号的频谱主瓣。上下限设得太窄会滤掉有用信号设得太宽会引入相邻频率的干扰。解决反复调试上下限频率一般取f1 ± df/2和f2 ± df/2作为通带边界。课设里经过多次调整才把失真降到可接受范围。4.3 抽样判决时刻不对齐现象解调出来的序列和原始序列错位误码率接近 0.5。原因抽样判决需要在每个码元的中间时刻采样如果采样时刻偏到码元边界判决结果会随机跳变。解决根据滤波器延迟和码元周期计算最佳采样时刻通常取Delay N/2附近。课设里这部分是请教后解决的说明抽样时刻的确定需要结合具体滤波器参数。4.4 awgn 信噪比修正遗漏现象误码率曲线和理论值对不上实际性能比预期差很多。原因awgn默认按信号功率比算 SNR而通信原理里的 SNR 定义在符号能量和噪声功率谱密度上两者差一个10*log10(N)和调制因子。解决在调用awgn时手动修正如代码里的SNR - 10*log10(0.5) - 10*log10(N)。不修正的话采样率越高偏差越大。4.5 Simulink 误码率模块报错现象Simulink 里计算误码率的模块始终报错找不到原因。原因Simulink 的误码率模块对输入信号维度、数据类型和采样率有严格要求配置稍有不一致就会报错。解决课设里最终放弃 Simulink 模块改用 MATLAB 代码symerr实现误码率统计。这也是个实用建议Simulink 适合搭链路看波形但精细的误码率统计用代码更可控。5. 二次开发与改造从课设代码到可复用的仿真脚本这份代码虽然是为课设写的但骨架是完整的稍加改造就能变成自己的仿真工具。我一般会从三个方向动手。第一个方向是参数化。把Fc、Fs、Fd、df、SNRpBit这些硬编码的值抽成函数参数封装成一个sim_2fsk(Fc, Fs, Fd, df, SNR_dB)函数返回误码率。这样扫 SNR 范围就变成一行循环SNR_range 0:2:20; ber_coh zeros(size(SNR_range)); ber_noncoh zeros(size(SNR_range)); for k 1:length(SNR_range) [ber_coh(k), ber_noncoh(k)] sim_2fsk(10, 40, 1, 10, SNR_range(k)); end semilogy(SNR_range, ber_coh, b-o, SNR_range, ber_noncoh, r-s); xlabel(SNR (dB)); ylabel(BER); legend(相干解调, 非相干解调); grid on;这样跑出来的曲线可以直接和理论公式对比验证仿真是否正确。注意semilogy用对数纵轴因为误码率通常跨好几个数量级。第二个方向是替换调制方式。dmod和ddemod支持多种调制把fsk换成psk或qam改一下M和df参数就能对比不同调制方式的性能。比如 2PSK 的误码率理论上是0.5*erfc(sqrt(SNR))比 2FSK 相干解调好 3dB跑一遍就能验证。第三个方向是加蒙特卡洛仿真。课设里numSymb 25太少误码率统计不准确。做性能曲线时至少要把码元数加到 10000 以上或者跑多次取平均。我一般会设numSymb 1e5然后循环 10 次取平均误码率这样曲线才平滑。改造方向关键改动适用场景参数化封装抽成函数返回误码率扫 SNR、对比不同参数替换调制方式改 dmod/ddemod 的调制参数对比 FSK/PSK/QAM蒙特卡洛增大 numSymb多次循环精确误码率曲线还有一个容易忽略的点rcosine和rcosflt的延迟补偿。如果要在滤波后做误码率统计必须把滤波延迟去掉否则解调序列和原始序列错位误码率会虚高。我一般会在rcosflt后面加一段截断把前Delay*Fs/Fd个点去掉再比较。从那以后我每次做通信仿真都强制先跑一遍无噪声的理想链路确认调制解调本身没问题再加噪声、加滤波、加延迟。这个顺序能省掉大量排查时间——如果理想链路都错了后面加什么都是白搭。希望帮到你。本文还有配套的精品资源点击获取