简介：面向通信原理课程设计场景的完整2FSK调制解调系统仿真文档，围绕MATLAB7.0平台完成设计与验证。内容按课程设计报告规范组织，包含设计任务与要求、方案选择论证、调制与解调原理、系统框图、源程序代码、仿真图形及性能分析：采用开关法产生2FSK信号，经加性高斯白噪声信道后以相干解调还原，并绘制信噪比与误码率关系曲线，便于对照理论结果评估系统性能。压缩包共1个文件，为doc格式文档，体积579KB，目录结构完整，可直接用Word打开阅读或按需修改其中参数。已有853人学习下载。对通信工程、电子信息类学生而言，可借此掌握2FSK系统建模思路、MATLAB仿真流程及课程设计报告写作规范，也能将其中代码与波形分析迁移到同类型数字调制系统设计中。
1. 2FSK调制解调：为什么相位不连续反而能成为经典教学案例
很多人在第一次接触2FSK时，以为它不过是“用两个频率代表0和1”这么简单。真正动手做MATLAB仿真时才会发现，坑全在解调端：明明加的是高斯白噪声，眼图却糊成一团；带通滤波器参数差一点，判决点就整体偏移；甚至载波频率不是某个特定倍数时，波形在一个码元周期内都凑不完整。这个课程设计恰好把2FSK最容易被忽略的两个底层特性暴露出来了——键控法产生的相位不连续，以及相干解调对本地载波同步的隐性依赖。搞清楚这两点，才算真正理解2FSK调制解调，而不是只会调用dmod和ddemod。
2FSK在数字通信里的位置很特殊：它实现最简单，但性能分析最完整。误码率公式可以通过相干和非相干两条路径推导出来，而且两条路径的差异恰好说明了“利用相位信息”到底带来了多少增益。这篇内容就从键控法的波形生成开始，拆到相干解调的判决机制，再落到MATLAB的仿真参数设置和误码率统计上。
2. 键控法生成2FSK：相位不连续是特征，不是缺陷
2.1 调频法和键控法的本质差异
2FSK信号有两种典型生成方式。调频法（VCO法）让载波频率随基带信号连续变化，对应的时间表达式是：
s(t) = cos(2π[f1 + k·m(t)]t + θ0)其中m(t)是二进制基带信号，k是频偏系数。这种方式产生的信号，在相邻码元切换时刻，载波相位是连续的。而本设计采用键控法（开关法），用二进制数据直接控制开关，在f1和f2两个独立频率源之间切换。由于两个振荡器各自独立运行，切换瞬间相位通常不连续，会产生相位跳变。
这个差异不是仿真的小细节，而是直接影响频谱结构。相位连续的2FSK，其功率谱能量更集中，旁瓣衰减更快；相位不连续的2FSK，因为存在相位跳变，频谱拖尾更严重，占用带宽更大。但键控法实现极其简单，在数字系统里就是“选通”，所以工程上反而用得更普遍。
2.2 MATLAB中dmod的参数语义
源程序里调制核心只有一行：
y = dmod(x, Fc, Fd, Fs, 'fsk', M, df);Fc=10是载频，Fd=1是码速率，Fs=40是系统采样频率，所以每个码元周期内采样点数N = Fs/Fd = 40。M=2表示二进制，df=10是两个载频的频偏。注意这里的df不是第二个载波的绝对频率，而是相对于中心载频的偏移量。所以实际两个载波频率是Fc - df/2和Fc + df/2，即5Hz和15Hz。
N = Fs/Fd这个比值是整个仿真链路的基石：它决定了每个比特对应多少个采样点，也决定了后续升余弦滤波器的群延迟设计。如果N太小，波形在码元周期内不足一个完整周期；如果N不是整数，滤波器设计和误码率统计都会出现奇怪的对不齐现象。常见的实现思路是先定码速率，再按整数倍关系选采样频率，避免分数采样点。
2.3 为什么f1和f2要“为pi的倍数”
源程序的调试记录里写了一句关键经验：振荡器频率f1和f2需要为pi的倍数，才能在一个时间单位内出现完整波形。这句话初看有点玄学，其实背后是离散采样点的周期对齐问题。
设载波频率为f，采样频率为Fs，那么一个载波周期占Fs/f个采样点。要让码元周期内的波形看起来“完整”，码元采样数N必须是Fs/f的整数倍。在Fs=40, N=40的条件下，若f=5Hz，一个载波周期占8个采样点，40个采样点正好5个完整周期；若f=15Hz，一个周期占40/15 ≈ 2.67个采样点，40个采样点相当于15个完整周期。后者的采样点在周期内的相位分布不均匀，画出来的波形就有明显的“毛刺感”。
更准确的说法是：采样频率必须满足奈奎斯特条件，同时采样点在每个载波周期内均匀分布才能得到光滑波形。当f2=15Hz而Fs=40时，单个周期内只采到约2.67个点，波形轮廓能看出来，但细致程度不够。把载频降为Fc=10, df=4左右，两个载波为8Hz和12Hz，采样点在周期内分布更均匀，波形更干净。课程设计里“调试多次才发现频率要取合适值”的血泪经验，本质就是离散采样对周期对齐的敏感性。
2.4 开关法2FSK的频谱特征
用键控法生成的2FSK信号，可以看作两个独立ASK信号的叠加：
s(t) = s1(t)·cos(2πf1t + θ1) + s2(t)·cos(2πf2t + θ2)其中s1(t)和s2(t)分别是二进制序列中“1”和“0”对应的基带脉冲。这个表达式意味着2FSK的功率谱是两套ASK谱的叠加，每套ASK谱的主瓣宽度取决于码速率Fd，而两个主瓣之间的间距取决于df。若df太小，两套谱重叠严重，接收端带通滤波很难分离；若df太大，带宽浪费严重。
典型的工程经验是取df ≥ Fd，即两个载频间距至少等于码速率。本设计中df=10, Fd=1，间距远大于码速率，两套谱完全分离，接收端即使滤波器做得很粗糙也能分开。这也是键控法的一个天然优势：两个频率源相互独立，载频间距想设多大就设多大。
3. 相干解调机制：从带通滤波到抽样判决的完整链路
3.1 相干解调为什么“相干”
2FSK的相干解调，核心思路是先让接收信号分别通过中心频率为f1和f2的带通滤波器，滤出各自对应的ASK分量，然后分别乘以同频同相的本地载波cos(2πf1t)和cos(2πf2t)，再低通滤除二倍频分量，最后在两个支路的低频输出之间做比较判决。
以f1支路为例，假设接收信号为r(t) = s1(t)cos(2πf1t + θ) + n(t)，乘以本地载波后得到：
r(t)·cos(2πf1t) = s1(t)·cos(2πf1t + θ)·cos(2πf1t) + n(t)·cos(2πf1t)积化和差后，第一项变成：
0.5·s1(t)·[cos(θ) + cos(4πf1t + θ)]低通滤波器滤除4πf1t分量后，剩下0.5·s1(t)·cos(θ)。注意这里出现了cos(θ)——这就是“相干”的代价：本地载波必须与发送端载波同频同相。若相位差θ = π/2，输出幅度直接归零，判决器什么都判不出来。这是相干解调和包络检波最本质的区别：包络检波不需要知道相位，相干解调必须同步。
3.2 为什么相干解调误码率更低
非相干解调的误码率公式是：
Pe_noncoh = 0.5 · exp(-SNR/2)相干解调的误码率公式是：
Pe_coh = 0.5 · erfc(sqrt(SNR/2))在高信噪比下，erfc函数按指数衰减且系数更优，所以相干解调的误码率始终低于非相干解调。物理直觉是：相干解调利用了载波相位这个额外信息维度，等于在判决前多做了一次“方向校正”，把噪声投影到信号正交方向上的分量也利用起来用于判决；而非相干解调只取包络幅度，把相位信息完全丢弃。
本设计在SNR=60dB的高信噪比下仿真，两种解调方式都接近零误码，但symerr统计出来的误码数通常是相干为0、非相干偶尔有零星错误。这也是为什么系统性能分析部分只给出了理论上相干优于非相干的结论，并明确说“本系统采用相干解调”。
3.3 ddemod中相干与非相干的调用差异
MATLAB中两种解调方式的调用差异在一个参数上：
z1 = ddemod(y, Fc, Fd, Fs, 'fsk/eye', M, df); % 相干 z2 = ddemod(y, Fc, Fd, Fs, 'fsk/eye/noncoh', M, df); % 非相干'fsk/eye'表示输出眼图数据的相干频移键控解调，'fsk/eye/noncoh'则是非相干版本。注意ddemod内部的实现并不真正做载波同步，而是假设接收端已知载波频率和相位，直接生成本地参考载波。这在仿真中是合理的简化，但也在提醒：实际系统里相干解调需要额外的载波恢复电路，比如PLL或者Costas环，这部分复杂度在仿真里被隐藏了。
3.4 抽样判决的位置选择
相干解调的最后一步是抽样判决。低通滤波器的输出在码元周期内并非恒定值，由于升余弦滤波器的群延迟，码元边界处存在过渡带。判决时刻应选在码元中心——也就是每个比特周期中间那个采样点——因为那里信号幅度最稳定，码间串扰最小。
源程序中dmod和ddemod默认在最佳时刻抽样，但在自己写解调器时会遇到这个问题。常见的实现思路是：先做匹配滤波，再以N = Fs/Fd为间隔分段，每段取中间点作为判决样本。Delay=3的升余弦滤波器引入约3个码元的延迟，判决时刻也要相应后移，否则会误判开头的几个比特。
4. 加入高斯白噪声：SNR修正项的来龙去脉
4.1 awgn的dB参数到底代表什么
源程序里加噪声那行是很多初学者最容易抄错、但又最值得深究的：
y = awgn(y, SNR - 10*log10(0.5) - 10*log10(N), 'measured', [], 'dB');SNR = SNRpBit/log2(M) = 60，因为M=2，所以SNRpBit和SNR数值相同。awgn函数的第二个参数要求输入“信噪比的dB值”，但这个信噪比的定义是信号功率对噪声功率的比值，且信号功率是按“每个采样点”还是“每个码元”计算，直接影响结果。
10*log10(0.5)这一项修正的是2FSK信号的占空比因子。2FSK信号在任一时刻只有一个载波在发送（另一个支路静默），平均功率是单载波的一半，即 −3dB。10*log10(N)修正的是码元能量和采样点能量的关系：每个码元包含N=40个采样点，码元能量是单采样点能量的40倍，即约16dB。把这两项从名义SNR里扣掉后，awgn内部按“每采样点信噪比”加噪声，最终等效到码元级时，实际比特信噪比才是60dB。
4.2 为什么选'measured'选项
awgn的'measured'选项会先测量输入信号的功率，再按目标SNR计算噪声方差。这么做的原因是：调制后的信号功率并非整数1，dmod输出的幅值和采样点数都会影响实际功率。若不用'measured'，而直接指定信号功率值，一旦载波频率和采样频率的比例不是整数，信号实际功率与估计值偏差可达百分之几十，最后算出来的误码率曲线就会整体偏移。
源程序中还设置了randn('state', seed(2))，即固定噪声生成器的随机种子。这是为了让仿真可复现——同一份seed=[12345 54321]下，每次运行生成的噪声序列完全相同。调试时固定种子极其重要，否则每次跑出来的眼图都不一样，没法判断是滤波器参数问题还是噪声随机性造成的。
4.3 SNR与Eb/N0的换算关系
通信原理教材里讨论误码率时用的信噪比通常是Eb/N0（每比特能量与噪声功率谱密度之比），而awgn给的往往是SNR。两者换算涉及三个量：码速率Fd、系统带宽B和过采样率N。
Eb/N0(dB) = SNR - 10*log10(B/Fd)在基带仿真的默认设置下，系统带宽近似等于采样频率Fs，所以B/Fd = N。这正是修正项10*log10(N)的来源。如果不做这个换算，直接用名义SNR画误码率曲线，理论曲线和仿真点会对不上。很多人的仿真相对于理论曲线差几个dB，原因就是这里少扣了一项。
5. 常见问题与排查：载频选择、滤波器失真和误码率统计
5.1 载频与采样率不匹配导致波形异常
现象：调制后的波形在码元边界处出现不规则的“凹陷”或“毛刺”，看起来频率切换点附近波形不完整；把t轴放大多个周期后，能明显看到有些码元周期内的载波周期数不是整数。
原因：采样频率Fs与载频f1/f2不是整数倍关系，导致每个码元周期内在同一相位位置反复出现采样点疏密变化。Fc=10, Fs=40时，Fs/f1=8（整数），Fs/f2=40/15≈2.67（非整数），所以f2支路的波形明显更粗糙。
解决：要么让载频为Fs的整数分之一（如取f2=20Hz或f2=10Hz），要么调整Fs到能同时整除两个载频的值。在Fc=10, df=10的配置下，把Fs提高为80，则Fs/f2=80/15≈5.33，仍非整数。更稳妥的做法是重新选取载频：如Fc=12, df=4，两个载波为10Hz和14Hz，Fs=40时40/10=4、40/14≈2.86，仍然不理想。实际课程设计中把df改成8得到载波2Hz和18Hz，或者改用Fc=10, df=10配合Fs=60，让两个载频都落在采样点均匀分布的位置上。
5.2 带通滤波器失真过大
现象：解调支路经过带通滤波器后，波形幅度明显衰减，且码元切换点出现严重的拖尾振荡，低通输出在判决时刻的幅度低于阈值，误判增多。
原因：带通滤波器的通带带宽设置过窄，落在通带边缘的频谱分量被大幅衰减；或者滤波器阶数过高，群延迟在通带内不恒定，造成波形畸变。Simulink中的模拟带通滤波器对f1和f2的间隔敏感，df越大越容易调。
解决：按“通带带宽≥2倍码速率”的经验设置上下限频率。对Fd=1Hz的码速率，带通滤波器通带至少要为 ±2Hz，即6dB带宽不小于4Hz。若使用FIR设计，则要关注滤波器阶数与群延迟的折中：阶数高则过渡带窄，但延迟大，判决时刻需要同步调整。
5.3 抽样判决实现不了
现象：低通滤波器输出波形已经恢复了基带信号的轮廓，但不知道用什么方法在时间轴上把码元值“抠”出来，直接对连续波形做阈值判决，结果在码元边界处不停翻转。
原因：低通输出的波形是连续的模拟信号，码元边界处存在过渡带和残余噪声，直接全局阈值判决会把过渡带里的随机波动当成码元变化。
解决：最可靠的实现是“先定时、再判决”。用Fd=1和Fs=40的约束，把信号每40个采样点分成一段，每段只取中点（第20个采样点）作为判决依据。这个位置的信号幅度最接近理想值，码间串扰最小。实现时用reshape(y(1:end-N), N, [])把流式信号变成矩阵，每列是一个码元，再取中间行做判决即可。若滤波器群延迟大，则要先丢弃前Delay*N个采样点再分段。
5.4 误码率统计模块报错
现象：Simulink里的误码率计算模块持续报错，换成代码实现后才正常。报错多发生在仿真时间与码元个数不匹配时。
原因：Simulink的误码率模块要求两个输入信号（发送序列和接收序列）具有完全相同的采样时间和帧格式。加入滤波器后，接收序列产生了群延迟，信号长度也变了，模块内部的对齐逻辑就会失效。
解决：改用MATLAB代码做误码率统计，即源程序中的symerr(x, z1)。但symerr也要求两个序列等长。ddemod的输出长度取决于输入长度，加入噪声后的y长度不变，所以这里能直接对齐。若自己实现滤波和解调导致长度变化，需要在对比前手动截掉滤波器延迟部分，并调整序列头部的偏移量。这就是源程序里注释Delay=3的含义：滤波器滞后3个码元，判决结果要和原始序列的前numSymb-3个码元对齐比较。
6. 设计取舍：2FSK为什么还在课设里屹立不倒
2FSK在频谱效率上几乎是最差的调制方式之一——两个独立载波各占一段频谱，总带宽约为2·df + 2·Fd，比同信息速率的2ASK翻倍。工程上还有QPSK、8PSK等高阶方案能用更小带宽传相同速率，为什么课程设计仍坚持用2FSK？
这个取舍背后是对“概念完整性”的重视。2FSK的调制端暴露了相位连续性和离散采样的矛盾，解调端暴露了相干与非相干对相位信息的依赖差异，噪声分析则暴露了SNR、Eb/N0和采样率之间的换算关系。这些概念在QPSK里同样存在，但被星座图和I/Q支路掩盖了。2FSK用最简单的方式把数字调制最底层的三个支撑点——载波、采样、判决——摊开放在桌面上。把这套机制吃透，后面学QPSK、MSK、OFDM时，遇到的无非是“多了几个载波”或“相位不再跳变”的变化。
从这个角度看，键控法产生的相位不连续反而比调频法的连续相位更有教学价值。它逼着你思考：相位跳变带来了什么代价（频谱扩展），又带来了什么便利（不需要VCO的线性度设计）。一次课程设计做完，能把这笔账算清楚，才算没白做。
我自己做这类仿真时，现在固定会先检查三件事：采样率与载频的整数倍关系、SNR修正项是否按码元能量折算、判决时刻是否落在码元中心。每一条都是从血泪经验里长出来的。从那以后，我每次拿到一个dmod/ddemod的仿真任务，都强制走一遍“波形→频谱→眼图→误码率”的完整链路，中间任何一个环节对不上，就先停下来查参数而不是继续调后面的模块。希望这些梳理过的底层细节，也能让你少走几步弯路。
本文还有配套的精品资源，点击获取
