简介这份资源面向通信原理学习者与数字信号处理入门者聚焦ASK、FSK、PSK、AM、PM、FM六种调制方式的仿真实现帮助读者在MATLAB环境中理解模拟与数字调制的原理差异、波形特征及解调思路。压缩包共5个文件均为.m脚本整体约2KB分别对应六种调制方式的仿真代码便于按需运行与对比观察。目前已有1595人学习下载说明其在通信课程实验与自学场景中具有一定参考价值。读者可借助这些脚本直观查看载波幅度、频率、相位随基带信号变化的规律理解二进制ASK、FSK、BPSK以及AM包络检波、PM相位检波、FM鉴频等关键环节并在此基础上修改参数、扩展多进制调制或加入噪声分析适合作为课程设计、实验报告与算法验证的起点素材。1. 从六种调制仿真说起ASK、FSK、PSK、AM、PM、FM 到底怎么在一套代码里跑通做通信仿真的人大多有过这样的经历课设题目写着“实现六种调制方式”打开 MATLAB 却不知道从哪下手网上搜到的代码要么只有一种调制要么参数对不上、波形出不来。这份资源就是冲着这个场景来的——它把 ASK、FSK、PSK、AM、PM、FM 六种调制方式放进同一套仿真框架里用统一的采样率、载波频率和符号速率跑通输出时域波形和频谱图。适合正在做通信原理课设、准备信号处理面试、或者想快速验证调制算法的人。你不需要从零推导公式但需要知道每种调制的参数怎么设、波形为什么长那样、频谱里的旁瓣意味着什么。接下来我按“先跑通再调参”的顺序拆一遍把能抄的代码和容易翻车的地方都摆出来。2. 六种调制的仿真框架统一参数与模块划分2.1 为什么要把六种调制塞进同一个框架单独写一种调制不难难的是六种放在一起还能对比。常见做法是每种调制写一个独立脚本结果采样率不一致、载波相位不统一画出来的频谱没法横向比较。这份资源采用“基带信号生成 → 载波调制 → 信道加噪 → 解调判决”的四段式结构六种调制共用同一套时间轴和采样参数。这样做的直接好处是你改一个fs或fc六种波形同步变化对比误码率时不会因为参数差异产生假象。具体到实现上基带信号统一用随机二进制序列生成符号速率Rb和采样率fs满足fs sps * Rb其中sps是每个符号的采样点数。载波频率fc一般取Rb的 4 到 10 倍避免频谱混叠。AM、PM、FM 属于模拟调制基带信号换成正弦波或三角波但采样框架不变。这种统一带来的边界是数字调制和模拟调制共用一套解调判决逻辑时模拟调制的“判决”其实是包络检波或鉴频代码里用不同函数分支处理。2.2 核心参数表与初始化代码下面这张表是六种调制共用的基础参数我按常见课设要求设了一组值你可以直接改参数符号典型值说明采样率fs200 kHz必须大于 2 倍最高频率符号速率Rb10 kbps数字调制的码元速率每符号采样点sps20fs / Rb载波频率fc40 kHz一般取 4~10 倍 Rb载波幅度Ac1归一化信噪比SNR10 dB加噪用仿真符号数N100总码元数初始化代码块如下% 统一参数初始化六种调制共用 fs 200e3; % 采样率 200 kHz Rb 10e3; % 符号速率 10 kbps sps fs / Rb; % 每符号采样点 20 fc 40e3; % 载波频率 40 kHz Ac 1; % 载波幅度 SNR_dB 10; % 信噪比 N 100; % 符号数 t (0 : N*sps - 1) / fs; % 统一时间轴 % 生成随机二进制基带信号 data randi([0 1], 1, N); baseband repelem(data, sps); % 每个符号重复 sps 次逻辑说明repelem把每个比特扩展成sps个采样点形成矩形脉冲。参数sps决定了波形的时间分辨率太小会导致边沿陡峭、频谱展宽太大会增加计算量。fc取 40 kHz 是为了在 200 kHz 采样率下每个载波周期有 5 个采样点满足奈奎斯特条件且波形平滑。如果你把Rb改成 20 kbpssps自动变成 10此时fc建议不低于 80 kHz否则一个载波周期采样点太少解调时相位误差会明显变大。2.3 六种调制的实现分支与关键差异数字调制 ASK、FSK、PSK 的差异在载波参数如何随基带变化。ASK 是幅度键控s_ask baseband .* cos(2*pi*fc*t)FSK 是频率键控需要两个频率f1和f2分别对应 0 和 1常见做法是f1 fc - dff2 fc dfdf取Rb/2保证正交PSK 是相位键控BPSK 用s_psk cos(2*pi*fc*t pi*baseband)这里baseband是 0/1 序列映射到 0 和 π 相位。模拟调制 AM、PM、FM 的基带换成连续正弦波m sin(2*pi*fm*t)fm取 1 kHz。AM 是s_am (1 m) .* cos(2*pi*fc*t)调制度设为 1PM 是s_pm cos(2*pi*fc*t kp*m)kp是相位灵敏度一般取 π/2FM 是s_fm cos(2*pi*fc*t kf*cumsum(m)/fs)kf是频率灵敏度取 2pi5e3 左右。注意 FM 里的cumsum(m)/fs是积分近似直接写cumsum会引入直流漂移这是很多人翻车的地方。% 数字调制示例BPSK s_psk cos(2*pi*fc*t pi*baseband); % 模拟调制示例FM注意积分用 cumsum 后除以 fs fm_tone sin(2*pi*1e3*t); % 1 kHz 基带 kf 2*pi*5e3; % 频率灵敏度 s_fm cos(2*pi*fc*t kf*cumsum(fm_tone)/fs);参数说明kf决定最大频偏kf/(2*pi)就是频偏值这里约 5 kHz。如果kf太大频谱会超出采样率范围产生混叠太小则频偏不明显解调后信噪比恶化。常见做法是先算卡森带宽BW ≈ 2*(df fm)确保fs 2*BW。3. 从波形到频谱解调与误码率验证的实操步骤3.1 加噪与解调六种方式不能共用一套判决加噪用awgn函数最省事但要注意awgn的 SNR 定义是信号功率与噪声功率之比而数字调制常用Eb/N0。换算关系是SNR Eb/N0 10*log10(Rb/fs)因为每个符号能量被摊到了sps个采样点上。如果你直接拿Eb/N0 10 dB去调awgn实际信噪比会低 13 dB 左右误码率曲线整体右移这是仿真里最常见的“玄学”偏差。解调分支ASK 用包络检波或相干解调相干解调乘cos(2*pi*fc*t)后过低通滤波再抽样判决FSK 用两个匹配滤波器分别匹配f1和f2比较输出能量大小PSK 相干解调乘载波后低通判决相位。AM 用包络检波加隔直PM 和 FM 用鉴相器和鉴频器实际仿真里 FM 常用diff(unwrap(angle(s_fm)))近似鉴频再除以2*pi得到瞬时频率。% 加噪注意 SNR 与 Eb/N0 的换算 EbN0_dB 10; SNR_dB EbN0_dB 10*log10(Rb/fs); % 换算后约 -3 dB rx awgn(s_psk, SNR_dB, measured); % BPSK 相干解调 demod_carrier cos(2*pi*fc*t); mixed rx .* demod_carrier; % 低通滤波截止频率取 Rb [b, a] butter(4, Rb/(fs/2)); filtered filtfilt(b, a, mixed); % 抽样判决 samples filtered(1:sps:end); decided samples 0;逻辑说明filtfilt做零相位滤波避免群延迟导致抽样点偏移。samples从每个符号的起始点抽取如果滤波器延迟没补偿抽到的可能是过渡点误码率会异常高。常见做法是先用grpdelay看滤波器延迟再调整抽样偏移。参数butter的阶数取 4太高会振铃太低则旁瓣抑制不够。3.2 频谱分析与误码率曲线绘制频谱用pwelch或fft都行我一般用pwelch因为窗函数和平均能压随机起伏。画的时候横轴用f/1e3转成 kHz纵轴用 dB。六种调制叠在同一张图里对比时ASK 和 AM 的频谱主瓣宽度约2*RbFSK 的带宽取决于频偏PSK 主瓣也是2*Rb但旁瓣衰减更快。FM 的带宽用卡森公式估算频谱上能看到明显的边带。误码率验证数字调制扫Eb/N0从 0 到 12 dB每个点跑 1000 个符号以上统计误码数。理论值 BPSK 是0.5*erfc(sqrt(EbN0))FSK 相干解调是0.5*erfc(sqrt(EbN0/2))ASK 相干解调同 BPSK 但非相干差 3 dB。仿真曲线和理论曲线对不上时先检查 SNR 换算再检查抽样时刻最后看滤波器带宽。% 误码率扫描示例 EbN0_range 0:2:12; ber zeros(size(EbN0_range)); for k 1:length(EbN0_range) snr EbN0_range(k) 10*log10(Rb/fs); err 0; for trial 1:100 data randi([0 1], 1, N); bb repelem(data, sps); tx cos(2*pi*fc*t pi*bb); rx awgn(tx, snr, measured); % ... 解调判决 ... err err sum(decided ~ data); end ber(k) err / (N * 100); end semilogy(EbN0_range, ber, o-);参数说明trial循环 100 次是为了平均随机性N取 100 时总符号数 10000误码率低到 1e-4 以下需要更多符号。常见做法是误码率低于 1e-3 时把N加到 1000 以上否则曲线尾部全是零看不出趋势。4. 避坑与排查六种调制仿真里最容易翻车的五个地方4.1 现象频谱出现镜像波形不对称原因采样率fs不够高或者载波频率fc超过fs/2。ASK 和 AM 的频谱关于fc对称如果fc接近fs/2上边带会被折叠到低频看起来像多了个杂散。解决确保fs 2*(fc Rb)数字调制里fc至少取4*Rb模拟调制里fc至少取10*fm。4.2 现象误码率曲线整体偏高和理论差 3 dB 以上原因awgn的 SNR 没换算成Eb/N0或者解调时抽样点没对准。先检查SNR_dB EbN0_dB 10*log10(Rb/fs)这行有没有写再检查samples的起始索引是不是1。如果滤波器有延迟抽样点要往后挪grpdelay个采样。解决用filtfilt代替filter或者手动补偿延迟。4.3 现象FM 解调后波形失真低频漂移原因cumsum积分没有除以fs或者没有去直流。cumsum(m)的幅度会随时间累积导致相位线性增长解调后出现斜坡。解决写成cumsum(m)/fs并在解调后加detrend或高通滤波。另外kf太大会导致相位卷绕unwrap之前先确认相位范围。4.4 现象PSK 星座图旋转判决全错原因载波相位不同步。仿真里发射和接收载波都用cos(2*pi*fc*t)但如果信道加了延迟或滤波器引入了相移接收相位会偏。解决在解调前用angle估计相位偏移并补偿或者改用差分 PSK。常见做法是加一个phase_offset变量手动扫 0 到 2π 看误码率变化。4.5 现象六种调制画在同一张图里坐标轴对不上原因每种调制的信号功率不一样AM 有直流分量FM 是恒包络直接叠图会互相遮挡。解决画频谱时统一归一化到 0 dB画时域波形时分subplot或者用不同颜色但限制ylim。另外注意 AM 的调制度如果大于 1 会过调包络出现反相解调后严重失真调制度设 0.8 到 1 之间比较稳。5. 进阶技巧用 App Designer 把六种调制做成可交互工具跑通脚本之后下一步自然是把参数暴露出来做成能拖滑块看波形的工具。MATLAB 的 App Designer 适合干这个但有几个细节值得说。第一六种调制切换用DropDown控件回调里用switch分支重新计算信号不要每次重建整个 UI。第二频谱刷新用drawnow limitrate否则拖动滑块时卡顿明显。第三误码率计算放在后台parfeval里跑避免阻塞界面。我自己的习惯是先在脚本里把六种调制的函数封装成modulate(type, params)和demodulate(type, rx, params)再在 App 里调用。这样脚本和 GUI 共用一套核心代码改参数不会两边不同步。参数面板里fs、fc、Rb、SNR各一个EditField加一个Run按钮触发计算。波形区用两个UIAxes一个画时域一个画频谱。误码率曲线单独一个窗口扫Eb/N0时显示进度条。% 封装后的调制函数签名 function tx modulate(type, data, params) switch type case ASK tx data .* cos(2*pi*params.fc*params.t); case FSK f1 params.fc - params.Rb/2; f2 params.fc params.Rb/2; tx cos(2*pi*(f1*(1-data) f2*data).*params.t); case PSK tx cos(2*pi*params.fc*params.t pi*data); case AM tx (1 data) .* cos(2*pi*params.fc*params.t); case PM tx cos(2*pi*params.fc*params.t pi/2*data); case FM tx cos(2*pi*params.fc*params.t ... 2*pi*5e3*cumsum(data)/params.fs); end end参数说明data对数字调制是 0/1 序列对模拟调制是连续波形。params结构体里存fs、fc、Rb、t。这样封装之后App 里只需要改params字段再调modulate就行。验证方法把type依次设成六种看时域波形是否符合预期——ASK 幅度跳变、FSK 频率跳变、PSK 相位跳变、AM 包络起伏、PM 相位随基带变、FM 频率随基带变。如果哪个波形不对先查data的维度是不是和t一致再查params有没有漏字段。从那以后我每次做调制仿真都强制先把六种调制的参数表列出来再写一行assert检查fs 2*(fc Rb)省得后面频谱出问题回头找半天。希望帮到你。本文还有配套的精品资源点击获取