简介这份《2FSK调制与解调系统设计与仿真》课程设计文档面向通信工程及相关专业学生完整展示了基于MATLAB7.0实现二进制频移键控系统设计的过程。内容涵盖信源随机比特流生成、开关法调制、高斯白噪声信道、相干解调链路建模并配有主要程序代码、仿真波形与眼图以及SNR-误码率性能分析思路可直接作为通信原理课程设计、实验报告的撰写模板或调试参考。资源包共1个文件为doc文档包体大小579KB适合作为轻量参考材料重点在于通过阅读代码与仿真流程理解2FSK系统的完整实现方法。该资源已有853人学习下载说明其在同类课程设计中的实用性与复用价值较高。读者可借助其中的设计任务、方案论证、源程序和结果分析四部分快速掌握从系统建模到性能评价的关键步骤。1. 2FSK 课程设计在做什么一次把调制、信道、解调全链路跑通如果你正在为通讯系统原理的课程设计发愁2FSK 这个题目十有八九是逃不掉的。它不像 QPSK 那样要处理相位模糊也不像 QAM 要考虑幅度衰落本质就是「用两个频率分别代表 0 和 1」看起来简单但真要把调制、信道加噪、解调、误码率统计这一整条链路在 MATLAB 里跑通并且能在答辩时讲清楚每个参数为什么这么设很多同学要花掉两周。这篇文章按我平时辅导课设的习惯把 2FSK 调制与解调系统设计与仿真拆成可复现的步骤先写一个手撸的调制脚本再做解调再搭 Simulink 模型最后用 Monte Carlo 曲线验证。适合大三通信工程学生、以及想把 FSK 仿真快速捡起来的从业者。2. 从数学式到 MATLAB 波形手写 2FSK 调制器与参数选择2.1 2FSK 的数学定义与频移指数 h为什么频偏不能随便设2FSK 的时域表达式可以写成发送比特 1 时输出一段频率为 f1 的余弦发送比特 0 时输出频率为 f0 的余弦。也就是说信号在时间轴上被切成一码元一码元的片段每一片的频率跟着比特走。这里面最关键的设计参数不是载频本身而是频偏 Δf f1 - f0 与码元速率 Rb 的比值工程上叫频移指数 h。h 的取值直接决定了频谱形状和解调难度。h 远大于 1 时两个频率在频谱上分得很开带通滤波器容易把两路分开代价是占用带宽大h 接近 0.5 时两个频率几乎正交频谱利用率高但滤波器分离的难度也大这就是后来 MSK 的思路。课程设计里最常见的做法是让 h 取 1 到 3 之间既有明显的双峰频谱又不会让滤波器设计太尴尬。我用的一组典型参数是Rb 1000 bpsf0 2000 Hzf1 4000 HzΔf 2000 Hzh 2采样率 fs 20000 Hz。这样每个码元有 20 个采样点频谱图上两根谱线落在 2 kHz 和 4 kHz视觉上非常清晰。还需要注意一个容易被忽略的点2FSK 有两种产生方式。一种是键控法直接用一个开关电路在两个独立振荡器之间切换相位不连续另一种是模拟调频法让一个振荡器的瞬时频率跟着基带信号连续变化相位连续。课程设计报告里两种都可以写但仿真时如果直接用cos(2*pi*f0*t)和cos(2*pi*f1*t)拼接波形在码元边界会有相位跳变FFT 之后频谱会多出一圈「裙边」。我一般推荐用连续相位方式生成波形也就是下面要写的cumsum法物理上更接近调频器写进报告也更有说法。2.2 手写连续相位调制脚本repelem 与 cumsum 的配合直接贴一段我在课设里常用的调制脚本注释里写了每个参数的作用% 参数定义 fs 20000; % 采样率 20 kHz Rb 1000; % 码元速率 1 kbps sps fs / Rb; % 每码元采样点数 20 N 200; % 码元数量总共仿真 0.2 秒 f0 2000; % 符号 0 的载频 f1 4000; % 符号 1 的载频 df f1 - f0; % 频偏 2000 Hz频移指数 h df / Rb 2 rng(1); % 固定随机种子保证结果可复现 bits randi([0 1], 1, N); % 随机生成 0/1 比特序列 n 0 : N*sps - 1; % 总采样点序号 % 连续相位 2FSK 调制先映射频率再积分出相位 freq_map f0 (f1 - f0) * bits; % 每个码元一个频率值 freq_seq repelem(freq_map, sps); % 扩展成每个采样点的瞬时频率 phase 2 * pi * cumsum(freq_seq) / fs; % 相位是频率的积分 tx cos(phase); % 输出通带波形这段代码的核心是repelem和cumsum的配合。repelem把每个码元对应的频率值重复 sps 次得到一个长度 N*sps 的瞬时频率序列cumsum对这个频率序列做累加再乘上 2π/fs就得到了连续的相位轨迹。这样生成的正弦波在任何时刻的相位都是连续的不会像键控法那样在码元边界出现跳变。h 2意味着两个载频相差 2 倍码速率这个值偏大带宽会比较宽但课程设计里展示频谱图时反而好看。如果你想在报告里顺便提一句「频谱效率」可以对比 h1 和 h2 的带宽但主线参数就用上面这一组因为它坑最少。没有 MATLAB 的同学可以直接用 GNU Octave 跑这段脚本repelem、cumsum、randi这些函数在 Octave 里都有几乎不用改。2.3 用 fskmod 对照验证输出是复基带不是通带波形手写代码最容易写错的是参数映射关系所以我会用通信工具箱的fskmod做一次交叉验证。fskmod的调用格式是y fskmod(x, M, freqsep, nsamp, fs)其中 x 是取值 0 到 M-1 的整数符号序列freqsep是相邻频率间隔单位是 Hznsamp是每个符号的采样点数。% 用 fskmod 生成复基带 2FSK 信号和手写版本做频谱对比 x_int bits; % fskmod 需要整数符号 0/1 y_base fskmod(x_int, 2, df, sps, fs); % 查看频谱确认两个频偏位置 figure; pwelch(y_base, [], [], [], fs); title(fskmod 复基带信号频谱);这里有个非常常见的误解fskmod的输出是复基带等效信号频率是相对于 0 的偏移不是真实的 2 kHz 和 4 kHz 通带波形。要在通带里显示得自己乘一个载波取实部。课程设计里如果只想验证误码率用复基带就够了但如果你要画「像课本里那样的 2FSK 波形」还得用手写代码生成的tx。两种方式各有用途fskmod省事、和理论误码率曲线对得上手写代码能讲清楚原理答辩时老师问「相位为什么连续」你能答上来。我一般是两个都做频谱图、波形图用通带信号BER 曲线用基带信号各取所长。3. 解调不是只有一种非相干、相干、过零检测的代码与误码率对比3.1 三种解调方式的原理与误码率公式2FSK 的解调方案比调制更值得写进报告因为它有三种常见路线每种都对应不同的实现成本和抗噪性能。下表是我在课设报告里直接用过的对比解调方式实现思路是否需要载波同步AWGN 下理论误码率非相干解调两路带通滤波 包络检波比较包络大小不需要0.5 * exp(-Eb/(2*N0))相干解调两路乘以本地载波 低通滤波比较输出需要且要同频同相0.5 * erfc(sqrt(Eb/(2*N0)))过零检测统计每个码元内过零次数次数多判为 f1不需要无闭式表达实测抗噪最差三种解调方式在低信噪比下的差异很明显。相干解调比非相干在相同 Eb/N0 下大约好 1 到 2 dB但因为要恢复本地载波实际系统里代价不小非相干是课设最稳妥的选择过零检测实现最简单但高噪声环境下过零次数会被噪声严重干扰通常只作为原理演示。写报告时把这条表格放进去基本能覆盖「方案选择依据」这部分。3.2 非相干解调代码分路滤波、包络检波与判决时刻非相干解调的关键是两路带通滤波器要能把 2 kHz 和 4 kHz 的能量分开。我建议滤波器带宽取 f ± Rb也就是每路带宽 1 kHz。下面的代码用butter设计带通再用filtfilt做零相位滤波这一步能帮你少踩一个延迟补偿的坑。% 两路带通滤波器中心频率分别在 f0 和 f1 [b0, a0] butter(4, [f0-Rb, f0Rb] / (fs/2), bandpass); [b1, a1] butter(4, [f1-Rb, f1Rb] / (fs/2), bandpass); % 零相位滤波filtfilt 不会引入群延迟判决时刻不用补偿 y0 filtfilt(b0, a0, tx); y1 filtfilt(b1, a1, tx); % 包络检波取绝对值再用滑动平均平滑 env0 abs(y0); env1 abs(y1); win round(sps / 2); % 滑动窗长 10 个采样点 env0 movmean(env0, win); env1 movmean(env1, win); % 在每码元中点采样并判决 sam_idx round(sps/2 : sps : N*sps); dec_bit double(env1(sam_idx) env0(sam_idx)); ber sum(dec_bit ~ bits) / N;判决时刻选在码元中点而不是末尾是因为码元边界附近有滤波器的瞬态响应包络幅度不稳定。采样点偏移几个采样点通常对结果影响不大但如果你在 Simulink 里做判决时刻和滤波器延迟一定要对齐这点第 5 章会专门讲。movmean的窗长取 sps/2 是我试过比较稳的值窗太小包络还有纹波窗太大码元间会互相串扰。这里的dec_bit是 0/1 序列直接和原始bits对比就能算出一组特定信噪比下的误码率。3.3 相干解调与过零检测本地载波假设与阈值推导相干解调的脚本在结构上和非相干很像只是把「带通 包络」换成了「乘法器 低通」% 相干解调假设本地载波与发送端同频同相 t_axis (0 : N*sps - 1) / fs; lo0 cos(2 * pi * f0 * t_axis); lo1 cos(2 * pi * f1 * t_axis); mix0 tx .* lo0; % 混频产生直流分量和 2 倍频分量 mix1 tx .* lo1; % 低通滤波器只保留直流附近的基带分量 [b_lp, a_lp] butter(4, Rb / (fs/2), low); z0 filtfilt(b_lp, a_lp, mix0); z1 filtfilt(b_lp, a_lp, mix1); sam_idx round(sps/2 : sps : N*sps); dec_bit double(z1(sam_idx) z0(sam_idx));这段代码的隐含假设是本地载波和发送载波完全同步。实际通信系统里这个同步要靠锁相环或导频来实现课程设计里可以直接「白嫖」一个理想载波但写报告时要说明这是理想假设。如果本地载波差了一个相位 φ输出会乘上 cos(φ)φ 接近 90 度时信号直接消失这就是为什么真实系统不敢轻易用相干解调。过零检测的思路则完全不同% 过零检测统计每个码元内信号穿过零点的次数 zc_event diff(sign(tx) ~ 0); % 过零位置标记 zc_count zeros(1, N); for k 1:N seg zc_event((k-1)*sps1 : k*sps); zc_count(k) sum(seg); % 每个码元的过零次数 end % 阈值推导频率 f 的正弦每秒过零 2f 次每码元 2f/Rb 次 th (f0 f1) / Rb; % 平均过零次数 6 dec_bit double(zc_count th);这个阈值 6 是这么来的f0 2 kHz 时每秒过零 4000 次每码元 4 次f1 4 kHz 时每秒过零 8000 次每码元 8 次。取中值 6过零次数大于 6 判为 1。这方法在无噪声时很漂亮但信噪比一低噪声会让信号反复穿越零点统计次数严重偏大。所以过零检测只适合做原理演示不适合作为课设的主打方案——但它出现在报告里会显得你做得全面。4. 用 Simulink 搭 2FSK 系统模块选型、SNR 设置与运行提速4.1 基础模块搭链路Switch 选频与两路带通滤波Simulink 里搭 2FSK 有两条路直接用通信工具箱的 M-FSK Modulator Baseband 和 M-FSK Demodulator Baseband两个模块一拖就完事或者用基础模块从零搭。课设答辩时老师往往会瞄一眼模型里有没有「黑匣子」所以我推荐用基础模块搭通带链路再用通信工具箱模块做交叉验证。基础链路的信号流向是伯努利二进制信源产生 0/1 → 送给 Switch 作为控制端 → Switch 在频率 f0 和 f1 两路正弦之间切换 → 加 AWGN 噪声 → 两路带通滤波器 → Abs 取包络 → 滑动平均平滑 → 码元中点采样判决。这里最关键的是 Switch 模块的参数Criteria 选u2 ThresholdThreshold 设为 0.5因为信源输出恰好是 0 和 1。如果你把信源设置成了 ±1那 Switch 的阈值就要改成 0这是很多人第一个翻车点。两路正弦波的 Sine Wave 模块一定要把 Sample time 设为1/fs也就是 0.00005 秒而不是默认的 0。默认的 0 表示连续时间信号会拖垮后面要说的求解器。带通滤波器用 Analog Filter Design 模块Design method 选 ButterworthFilter type 选 Bandpass上下截止频率分别是 f ± Rb。两路滤波器的参数必须严格对称否则解调判决会出现固定偏移。4.2 AWGN 信道 SNR 与 Eb/N0 的换算Simulink 里填哪个数这是整个 Simulink 仿真里最容易让人懵的地方。AWGN Channel 模块里让你填的 SNR 是「每个采样点的信号功率与噪声功率之比」不是论文里常说的 Eb/N0。两者差一个每码元采样点数 sps。如果直接在模块里填 10 dB实际等效的 Eb/N0 会变成 10 10*log10(sps) dB也就是说你以为是 10 dB实际等效是 23 dBBER 曲线当然对不上理论值。换算公式按复基带模型写是SNR_dB EbN0_dB - 10log10(sps)。例如 sps 20 时10log10(20) ≈ 13 dB所以想模拟 Eb/N0 8 dB采样点 SNR 应该填 -5 dB。通带实信号模型还要再考虑 3 dB 的单边带噪声问题工程上更稳的做法是先按这个公式填跑一两个点对比误码率理论值如果偏差固定就做一个查表校准——这比在报告里硬推公式省力还显得你有工程意识。如果你用的是通信工具箱的 M-FSK Modulator Baseband 模块AWGN Channel 和 Error Rate Calculation 之间会有模块内部的 Eb/N0 换算不用你操心但这也意味着模型的可解释性变弱。我的建议是基础模块链路用来截图和展示原理通信工具箱链路用来快速出误码率数据两者互补。4.3 求解器与仿真时长设置让仿真在 1 分钟内跑完Simulink 模型跑得慢九成是求解器设置问题。只要链路里有连续时间信号比如某个 Sine Wave 的 Sample time 留了 0Simulink 就会用变步长连续求解器而通信系统里信号频率高、仿真时长又长步长会被压得很小跑 1 秒仿真可能要几分钟。解决方法是把所有的 Source 模块都改成离散采样时间然后在 Solver 里选Fixed-step求解器类型选Discrete (no continuous states)步长填1/fs。仿真时长一般填 0.2 到 0.5 秒就够了。200 个码元 × 0.001 秒 0.2 秒足够统计出误码率的一个点想画 BER 曲线就循环改 SNR 多次运行每次重新初始化随机种子。如果老师要求看更长时间的波形可以单独拉一个 Scope 看前 0.05 秒不必全程保留示波器数据——Scope 模块默认会存所有采样点存 10 万个点界面就开始卡了。顺便说一句如果你后续想用 Multisim 或 Spice 这类电路级工具搭 2FSK逻辑完全一样只是带通滤波器要换成 RLC 谐振回路包络检波换成二极管和电容这会是个很扎实的加分项。5. 2FSK 仿真避坑指南5 个最容易翻车的细节与排查方法5.1 频谱图上两条谱线糊成一片频偏与 FFT 分辨率两个原因现象用 pwelch 或 FFT 看生成信号2 kHz 和 4 kHz 处看不到两个清晰的峰整个频谱像一团糊状。原因有两种。第一种是频移指数 h 太小比如你设了 Rb 1000、f0 1000、f1 1200Δf 200 Hz两路带通滤波器根本分不开频谱自然糊。第二种是 FFT 点数太少或没加窗频率分辨率不够。解决把 h 调到 1 以上用我前面推荐的那组参数频谱分析时用pwelch(tx, hann(1024), [], [], fs)加窗能显著压平旁瓣。先改参数再看窗别一上来就怀疑代码写错。5.2 解调波形整体延迟一个码元开头误码一片滤波器延迟没补偿现象无噪声情况下解调出错的比特都集中在开头或者判决结果比发送比特整体晚了一个码元位置。原因是带通滤波器是因果系统输出相对输入有群延迟如果你用的是filter而不是filtfilt在码元边界采样必然取到上一码元的残余。解决有两个办法一是把filter换成filtfilt零相位滤波不引入延迟课程设计仿真里这是最省事的手段二是如果非要用filter用grpdelay(b, a, 1)算出群延迟点数把判决索引整体后移。我在课设里一开始没注意这个前 5 个码元全错后来一查是滤波器延迟了 8 个采样点把采样的起点往后挪 8 个点就正常了。5.3 BER 曲线永远离理论值差几个 dBSNR 换算少了符号采样数现象Monte Carlo 仿真的误码率曲线和理论曲线形状一致但整体右移偏乐观或左移偏悲观固定几个 dB。原因几乎都是 SNR 和 Eb/N0 换算错误最常见的是忘了除以 sps。比如你用awgn(tx, 10)以为加了 10 dB 噪声然后和理论0.5*exp(-Eb/(2*N0))对比时误把 10 dB 当 Eb/N0。解决统一按SNR_dB EbN0_dB - 10*log10(sps)换算并且做一次自检——把 SNR 设到 30 dB 以上误码率必须降到 0如果不是 0说明解调链路里还有 bug别急着跑曲线。5.4 Simulink 仿真慢到怀疑人生连续模块与离散步进混用现象模型结构没问题波形也能出来但跑一次要几分钟改一个参数又要重跑。原因是部分模块用了连续时间Sample time 0Simulink 被迫用变步长求解器而通信信号频率高步长被压得极小。解决把所有 Source 模块的 Sample time 改成1/fsSolver 选Fixed-stepDiscrete (no continuous states)步长设为1/fs。改完之后0.2 秒的仿真基本秒出。另外Scope 的 History 不要勾选全部数据点否则数据量大了照样卡。5.5 过零检测在低信噪比下直接罢工噪声过零与阈值失效现象信噪比低时过零检测判决结果几乎全是 1或者误码率接近 0.5。原因是噪声会让信号在真实过零点之外反复穿越零电平过零次数统计严重膨胀固定阈值 6 已经完全失去区分度。解决过零检测只在高信噪比下用比如 Eb/N0 12 dB 时做演示如果一定要在低信噪比下用先对接收信号做带通滤波再数过零可以把一部分噪声毛刺滤掉但改善有限。这个现象本身很适合写进报告的分析部分为什么实际系统不常用过零检测答案就在这。6. 验证与进阶Monte Carlo 误码率曲线、眼图和 MSK 扩展方向6.1 Monte Carlo 误码率脚本包含理论与仿真两条曲线验证 2FSK 仿真是否正确的最终手段是把误码率曲线跑出来和理论公式对比。下面的脚本用fskmod/fskdemod在复基带上做 Monte Carlo跑完直接画图EbN0_dB 0:2:12; Nbits 5000; % 每个信噪比点仿真的比特数 ber_sim zeros(size(EbN0_dB)); for i 1:length(EbN0_dB) EbN0 10^(EbN0_dB(i)/10); SNR_dB EbN0_dB(i) - 10*log10(sps); bits randi([0 1], 1, Nbits); tx_base fskmod(bits, 2, df, sps, fs); rx_base awgn(tx_base, SNR_dB, measured); rx_bit fskdemod(rx_base, 2, df, sps, fs); ber_sim(i) sum(rx_bit ~ bits) / Nbits; end % 理论曲线相干 2FSK 和 非相干 2FSK EbN0_lin 10.^(EbN0_dB/10); ber_coh 0.5 * erfc(sqrt(EbN0_lin / 2)); ber_non 0.5 * exp(-EbN0_lin / 2); semilogy(EbN0_dB, ber_sim, o, EbN0_dB, ber_coh, -, EbN0_dB, ber_non, --); legend(仿真, 相干理论, 非相干理论);fskdemod在 MATLAB 中实现的是最大似然检测结果会更贴近相干理论曲线非相干曲线作为下界参考。每个信噪比点 5000 个比特在 12 dB 时大约能统计到 10 个误码曲线毛刺属于正常现象想要更光滑就把 Nbits 加到 20000代价是低信噪比点跑得慢。这条脚本可以直接拍进报告附录答辩老师问起「为什么仿真点落在两条理论曲线之间」你就能解释最大似然检测比简单包络检波增益了约 1 dB。6.2 眼图与频谱图答辩时最能撑场面的两张图除了 BER 曲线课设报告里一定要放眼图和频谱图。眼图用eyediagram(rx_base, sps)生成它能直观展示码间串扰和噪声容限信噪比高时眼睛张开大信噪比低时眼皮变厚。频谱图用pwelch(tx, hann(1024), [], [], fs)能看到 2 kHz 和 4 kHz 两根谱线以及带外滚降。画图时记得加grid on把f0和f1标出来截图放在报告里会显得你确实理解了信号的频域特征。我习惯把两张图放在一个 Figure 里左右对比无噪声和有噪声各一组这样答辩时不用多解释图自己会说话。6.3 从 2FSK 到 MSK/GFSK一个能说清楚原理的扩展方向如果课设要求里写着「可选做提高部分」优先做 MSK 和 GFSK而不是 QPSK。MSK 就是把频移指数 h 从 2 降到 0.5实现最小频移的正交 FSKMATLAB 里fskmod(bits, 2, Rb/2, sps, fs)就能生成GFSK 是在 FSK 之前加一个高斯滤波器把相位轨迹磨平蓝牙物理层用的就是 GFSK。扩展部分不必重新搭一整条链路只要在现有脚本里把df换掉重新生成眼图和频谱再对比 2FSK、MSK、GFSK 的带宽这份报告就不只是「完成了课程设计」而是有了对比实验的深度。我第一次做 2FSK 课设时频谱图一出来就懵了两道谱线糊在一起调了整整一个晚上才发现是 h 设成了 0.2频偏小于码速率滤波器和频谱一个都不对。后来凡是做数字调制仿真我都会先把 h、判决时刻、SNR 换算这三件事写在纸上再动笔。这套流程陪我过了很多次课设和项目预研希望对你有同样的帮助。本文还有配套的精品资源点击获取
