基于MATLAB的BiAudio仿真电台：从傅里叶变换到双路AM调制解调全链路解析

简介：面向信号与系统课程设计的Matlab仿真电台项目，以Biaudio为主题，综合运用音频读取、播放控制、界面交互等知识，适合高校相关专业学生作为课设参考或二次开发基础。压缩包共33个文件，约34.91MB，核心包括MATLAB工程与脚本文件（.prj/.m）、演示音频（.mp3）、界面与素材图片（.png/.jpg）、辅助说明文档（.txt/.md）以及可直接运行的exe工具等，目录结构清晰，便于理清代码与资源关系。已有165人参与学习下载。内容涵盖完整项目源码、figure界面文件、编译脚本和音频素材，可直接运行观察效果，也可按需修改信号处理与播放逻辑，帮助理解“信号与系统”在实际电台仿真中的应用，是课设选题、代码研读与功能扩展的实用资料。 每年做信号与系统课设，最怕的就是题目听起来高大上，翻开课本却不知道从哪行代码开始写。BiAudio仿真电台这个题，表面是一个"仿真广播电台"的GUI演示项目，实际上它把傅里叶变换、调制解调、滤波器设计、系统响应这几章的核心知识全部串成了一条完整链路。做完这个项目，你对调幅广播为什么会有"串台"、收音机里沙沙声是哪来的、滤波器阶数到底影响什么，都会有非常具体的体感。

我把这个项目从选题到答辩完整拆一遍，包含方案取舍、参数推导、关键代码、踩坑记录和答辩追问准备。如果你是信号与系统课设里选了类似题目的学弟学妹，或者想用MATLAB做一套完整的信号链路练手，这篇可以帮你省下大量试错时间。

1. 这个课设到底做出来了个什么东西

1.1 题目的本质不是"做电台"，而是"用代码重演无线电广播"

大家看到"仿真电台"四个字，第一反应往往是做个漂亮界面，放两首歌进去，看起来像电台就行。但课程设计的评分老师真正想看的是：你能不能把"语音信号如何变成可发射的已调信号、已调信号如何在噪声信道里传输、接收端又如何把它从噪声和干扰中恢复出来"这一整套物理过程，用MATLAB完整地实现并演示出来。

BiAudio这个名字我自己的理解是"双音频通路"。也就是说，这个仿真电台系统里同时存在两路独立的音频节目源——比如A路是一段人声播报，B路是一段背景音乐。两路信号分别经过调制、混合、加噪、传输，在接收端再分别解调恢复。这样做的好处是：单路AM调制的课设代码网上到处都是，但双路同时传输、互不干扰，就需要你真正理解"频分复用"和"滤波器频率选择性"这两个概念，评分档次完全不同。

1.2 完成的系统从外到内是什么样

最终交付物是一个MATLAB App Designer应用程序，打开后是一个模拟电台控制面板。左侧是参数设置区，包括两路音频文件的选择、载波频率、信道信噪比滑块；中间四个坐标轴分别显示发射端频谱、接收端频谱、A路解调波形、B路解调波形；底部是"开始发射""播放A路""播放B路""停止"四个按钮。

点击"开始发射"后，系统按顺序完成：读取音频->限制基带带宽->两路AM调制->叠加合成->添加噪声->接收端带通分离->包络检波->低通滤波恢复音频。整个过程的信号波形和频谱会实时刷新在界面上，你甚至可以拖动信噪比滑块，亲眼看着接收波形从干净慢慢变成毛刺丛生，再点击播放听一听"滋滋"的噪声对语音的影响。

这里面每一块背后都对应课本上的一个知识点：频谱图对应傅里叶变换，带通滤波器对应频率选择性系统，调制解调对应卷积定理和频移性质。做完这套，课设报告里"原理分析"部分基本不用发愁没东西写。

2. 方案设计：为什么选了AM加频分复用加包络检波

2.1 调制方式对比之后的选择

信号与系统课设里能用的调制方式不少，AM、DSB、SSB、FM都能实现。我最终选了常规AM而非FM或DSB，原因很现实：常规AM的接收端可以用包络检波，也就是只做一次取绝对值和低通滤波，代码简单、原理直观，还能在报告里画经典的"包络线恢复原信号"示意图。

FM虽然抗噪声性能更好，但频率调制的瞬时频率和原信号成正比，接收端无论用鉴频器还是锁相环，对课程设计这个体量来说都有点重，调试时间会成倍增加。DSB又因为没有载波分量，接收端必须用相干解调，涉及载波同步问题。做课程设计最怕的就是"所有时间花在同步上"，所以常规AM配包络检波是性价比最高的选择。

要注意的是，常规AM的调制指数必须严格控制。调幅系数大于1时会出现"过调幅"，包络检波出来的波形和原始信号不再成正比，听起来就是严重的破音。这个坑我在后面调试部分会专门讲。

2.2 两路信号共存靠FDM而不是TDM

两路音频如何互不干扰地在同一个信道里传输，有两种直观方案：频分复用（FDM）和时分复用（TDM）。我选择了FDM，也就是给两路信号分配不同的载波频率，让它们在频谱上错开。

FDM的物理意义非常明确：A路调到8kHz载波上，B路调到16kHz载波上，只要两个载波离得足够远，两路信号的频谱就不会重叠。接收端用两个不同中心频率的带通滤波器，就能像分拣机一样把混在一起的信号重新分开。这个过程和"收音机调台"的直觉完全一致，答辩演示时观众一看频谱图就明白。

TDM也可以做，但那是数字通信的典型思路，需要精确的时隙同步。模拟仿真里做TDM，还得自己设计收发端的同步机制，属于给自己挖坑。课程设计的本质是验证课堂理论，FDM和AM组合起来，用到的恰恰是傅里叶变换频移性质，最贴合"信号与系统"这门课的主线。

2.3 解调策略如何兼顾"好写"与"好讲"

解调方案我最终实现了包络检波作为主流程，同时把相干解调留作扩展模块。包络检波的数学原理很简单：AM信号 s(t)=(1 + mx(t))cos(2πfct)，取绝对值后，载波的正负半周全被翻到正半轴，再经过一个截止频率低于载波、高于基带信号最高频率的低通滤波器，就能得到(1 + m*x(t))这个缓变包络。去直流后，原信号x(t)就回来了。

相干解调则需要接收端本地产生一个与发射端同频同相的载波，把接收信号乘上去，再用低通滤波取出基带分量。这个方案多了"载波同步"这一层麻烦，因为仿真里可以用同一个载波发生器，所以实现不难，但在答辩时容易被追问"实际系统中载波不同步怎么办"。我建议基础版用包络检波，把相干解调作为报告里的"改进与扩展"，既展示你理解了两种解调方式，又不至于把复杂度堆高。

3. 信号链路各模块的实现与参数推导

3.1 先把系统参数整整齐齐定下来

做信号链路仿真，第一步不是写代码，而是算参数。参数设计得离谱，后面所有波形都会乱掉。我最终用的参数如下表：

最关键的推导是载波频率的选取。AM调制后每路信号占用双边带，带宽等于2B=6000Hz。A路用8kHz载波，频谱占据5kHz到11kHz；B路用16kHz载波，频谱占据13kHz到19kHz。两路之间从11kHz到13kHz留了2kHz保护间隔，这样接收端带通滤波器不需要做到极端的陡降，就能把两路分开。

采样率48kHz则保证了19kHz的上边带远小于24kHz的奈奎斯特频率。如果你直接把原始44.1kHz采样率的音频读进来不重采样，16kHz载波导出的19kHz分量就非常贴近奈奎斯特边界，滤波器很容易失真。所以代码里统一做了resample到48kHz。

3.2 发射端代码的结构与关键点

发射端代码的核心思路是："先限带，再调幅，最后叠加"。限带操作很多人会忽略，直接拿完整频响的音频去做调制，结果就是两路信号边带互相重叠，怎么滤波都分不开。

% 读取音频并统一到48kHz采样率 [x1, fs_orig] = audioread('voice.wav'); x1 = resample(x1, fs, fs_orig); x1 = x1 / max(abs(x1)); % 幅值归一化 % 限制基带带宽，只保留0~3kHz lpFilt = designfilt('lowpassfir', 'FilterOrder', 64, ... 'CutoffFrequency', 3000, 'SampleRate', fs); x1 = filter(lpFilt, x1); % 生成A路AM已调信号 t = (0:length(x1)-1).' / fs; m = 0.6; s1 = (1 + m * x1) .* cos(2*pi*fc1*t);

第二路B完全同理。之后把两路已调信号相加，就得到发射端的复合信号。

这里有个容易被忽视的细节：载波信号cos(2pifc*t)的长度必须和音频信号对齐。很多初学者会习惯性地用linspace生成时间轴，一旦音频长度不是偶数，t和x1就无法逐元素相乘，报错"矩阵维度不一致"。用(0:length(x1)-1).'/fs这种方式生成列向量时间轴是最稳妥的。

在频域观察发射端复合信号时，你会看到三个清晰的谱峰族：8kHz附近是A路双边带，16kHz附近是B路双边带，两者之间有明显空隙。这个频谱图就是答辩报告里最有说服力的"中间产物"截图。

3.3 信道模块：噪声功率怎么算才正确

信道模型我选择了加性高斯白噪声（AWGN）。实现很简单，但"信噪比如何转化为噪声幅度"这一步很多人搞错。

信噪比SNR的定义是信号功率与噪声功率之比，单位是dB。所以要先算出复合信号s_total的平均功率，再反推噪声功率：

signal_power = mean(s_total.^2); snr_linear = 10^(snr_db/10); noise_power = signal_power / snr_linear; noise = sqrt(noise_power) * randn(size(s_total)); r = s_total + noise;

注意噪声标准差是sqrt(noise_power)，不是noise_power本身。我曾经看到有人直接把噪声功率当成标准差乘进去，结果信噪比即使设30dB，接收端也全是雪花噪声。这个错误在波形图上非常容易露馅：30dB应该基本无感，你的图上却啥都看不清，那多半就是噪声幅度差了数量级。

更精细的做法是分别计算有用信号和噪声的功率谱密度，但对于课设来说，时域平均功率法已经把概念讲清楚了。

3.4 接收端：带通分离、包络检波、低通恢复

接收端是整个系统里滤波器发挥作用最集中的地方。

% 带通提取A路 bpfA = designfilt('bandpassfir', 'FilterOrder', 256, ... 'CutoffFrequency1', 5000, 'CutoffFrequency2', 11000, ... 'SampleRate', fs); yA = filter(bpfA, r); % 包络检波：取绝对值后再低通 envA = abs(yA); envA = filter(lpFilt, envA); % 低通截止3kHz envA = envA - mean(envA); % 去直流，恢复零均值原信号

这段代码每一行都有它的物理含义。带通滤波器的作用是"频率选择性接收"，只让A路所在频段的信号通过，B路和大部分噪声被压掉。取绝对值实现了非线性整流，这是包络检波的核心。低通滤波器则负责把整流后残留的载波频率成分（8kHz附近的能量）滤除，只留下包络信号。去直流是因为AM的包络是(1+m*x(t))，其中的直流分量不包含音频信息，必须去掉才能得到x(t)。

这里滤波器阶数的选择要讲一下。256阶FIR听上去很大，但在48kHz采样率下，5kHz和11kHz之间的过渡带约20kHz带宽，256阶的过渡带宽度大概在几个kHz量级，足以分隔。如果你用32阶，频率选择性太弱，A路带通滤波器会把B路的19kHz边上带残骸也放进来一部分，听感就是"串台"。

3.5 解调效果怎么量化评价

课设报告里光贴波形图还不够，最好有个量化指标。我用两个指标评价解调质量：一是恢复信号与原始信号的波形相关系数，二是解调前后的信噪比改善量。

R = corrcoef(x1, envA_normalized);

相关系数越接近1，说明波形保真度越高。实测在无噪声环境下，包络检波的相关系数能做到0.95以上。噪声加到10dB时，相关系数会跌到0.8左右，听感上已经有明显背景噪声但语音仍可辩别。

这个"计算相关系数"的操作在答辩时非常加分，因为它表明你不是只会画图，而是会用统计量评估系统性能。

4. 从命令行到GUI：用App Designer搭一个能交互的电台操作台

4.1 界面布局和数据流设计

命令行脚本跑通信号链路之后，下一步是把它搬进App Designer。我的布局是经典的三段式：左边参数区，中间显示区，底部操作区。控件不多，但数据流要理清楚。

在App Designer里，我建议把所有中间信号放在一个名叫app.data的结构体里，字段包括fs、x1、x2、s_total、r、env1、env2、corr等。这样回调函数之间传数据非常方便，不会在多个handles之间来回倒腾。startupFcn里预加载两个默认音频文件，界面一打开就有东西可以演示，不用每次手动选文件。

4.2 按钮回调与播放音频的一个大坑

"开始发射"按钮的回调函数，本质就是把命令行脚本按顺序执行一遍，然后把结果填到app.data和各个坐标轴上。这个过程不复杂，但播放音频时容易踩坑。

MATLAB的play函数播放声卡对象是异步的，如果你在按钮回调里连续创建两个audioplayer对象播放A路和B路，后创建的对象可能中断前一个。更隐蔽的问题是audioplayer对象是句柄类，如果只是创建后不保存引用，它可能被MATLAB的垃圾回收机制提前销毁，表现就是点击播放按钮后声音刚响起就戛然而止。

解决办法是把player对象存成App的属性：

app.playerA = audioplayer(app.data.env1, app.data.fs); play(app.playerA);

并且播放前先stop已有的player对象。我在这个坑上耗了两个小时，最后发现就是对象生命周期的问题。

4.3 用Timer实现频谱"实时刷新"

如果每次拖动滑块都触发一次完整的重算加绘图，整个界面会卡得没法用。我采用的做法是：滑块回调里只更新app.data的SNR值并重新加噪解调，频谱图的刷新由一个定时器timer控制。

app.timer = timer('Period', 0.2, 'TimerFcn', @(~,~) app.updateDisplay()); start(app.timer);

每200毫秒刷新一次显示，既流畅又不至于让CPU持续满载。这里要格外小心的是App Designer关闭时，一定要在closeRequestFcn里delete这个timer，否则关掉界面后MATLAB命令窗口会不断报错"无效的定时器对象"。

我还加了"关闭界面自动清理"的逻辑，这一步虽然小，但能让项目显得非常完整，有的老师会专门检查程序退出是否干净。

5. 调试阶段踩过的坑和答辩实战准备

5.1 排查"串台"问题的完整链路

系统第一次调通后，我播放解调后的A路，竟然隐隐约约能听到B路音乐的声音。这就是经典的邻道干扰。我没有直接去调滤波器系数，而是按"信号越界排查"的思路一步步来。

第一步，把发射端复合信号的频谱放大看，确认两路信号有没有在频域重叠。如果B路音频没有做3kHz限带，它的边带会远远超出设计范围，直接插到A路的频段里。限带滤波器一加，这个源头就堵住了。

第二步，看带通滤波器的幅频响应。A路滤波器5kHz到11kHz通带，在13kHz处应该已经有明显衰减。如果阶数太低或者过渡带太宽，B路的13kHz边带会残留在输出里。我最初用32阶FIR，13kHz处只衰减了10dB左右，B路声音自然清晰可闻；换成256阶之后，带外衰减超过40dB，串台彻底消失。

第三步，验证包络检波后的低通滤波是否足够抑制残余载波。如果低通截止频率选高了，比如5kHz，那载波分量会和音频混在一起，表现为"高频吱吱声"而不是串台。

这个排查过程本身就是报告里很好的"问题分析与解决"章节素材：你展示了从频域定位、到滤波器参数调整、再到效果验证的完整闭环。

5.2 过调幅失真：看上去很美的波形其实已经废了

另一个让我印象深刻的坑是调幅系数设置不当。我为了图省事直接设m=1，结果在信号幅度较大的时刻，出现了过调幅。包络检波后的波形和原始信号差得离谱，播放出来就是"扑扑"的爆破声。

原理上，(1 + mx(t))这个因子在mx(t)小于-1时会变负，AM信号在这里发生相位反转，包络不再是原信号的线性映射。解决方式有两个：一是把m设为0.6到0.8，留足余量；二是在发射端加一个保护判断：

if max(abs(m * x1)) >= 1 error('调制系数过大，请减小m'); end

课程设计的实际意义往往就在这种细节里：不是代码跑通就完事，而是要让每个参数都工作在合理区间。

5.3 答辩高频追问和参考应答思路

做完项目到答辩环节，老师最常问的问题基本聚焦在"为什么这么选"和"如果那样会怎样"两个方向。我把高频问题整理一下：

• "为什么用AM不用FM？"应答核心：AM接收端可以用包络检波，实现简单，演示直观；FM需要鉴频器或锁相环，抗噪声虽好但复杂度高，作为课程设计选择AM能突出信号与系统课程的核心知识点。
• "为什么频分复用而不是时分复用？"应答核心：FDM在频域直接区分，和AM、滤波器这些本课程知识点衔接紧密；TDM需要严格同步，更多用于数字通信。
• "滤波器为什么用FIR不用IIR？"应答核心：FIR具有严格线性相位特性，信号波形经过滤波后不会产生相位失真，这对音频解调非常重要；IIR阶数低效率高但相位非线性。
• "系统总带宽怎么算？"应答核心：每路AM双边带带宽2B为6kHz，两路是12kHz，再加2kHz保护间隔，总带宽14kHz。如果能当场写出这个计算，基本稳了。
• "信道噪声再大一点会怎样？"这时候直接拖动SNR滑块到5dB，展示接收波形和播放效果，比任何口头解释都有说服力。

5.4 代码结构的一些个人习惯

最后分享几个让项目整体质量提升的小习惯。第一，脚本文件按功能命名：main_app.m、transmitter.m、channel.m、receiver.m、evaluate.m，每个模块职责单一，报告里贴代码时也清晰。第二，所有参数集中在文件头部或者一个params.m脚本里，修改参数不需要满文件搜索。第三，中间关键节点都用save存一次.mat文件，方便调试时直接加载某个阶段数据，不用每次都从头跑一遍完整链路。

这些习惯不是课程设计的硬性要求，但对后续找工作面试时展示项目经历非常有帮助：当你被问到"项目里代码怎么组织的"，能讲出模块划分和参数管理思路，和只回答"写了几个脚本文件"完全是两个水平。做完这个BiAudio仿真电台，你会发现再回头看书本上那些公式，脑子里会浮现出真实的波形和声音。这就是课设最大的价值。

本文还有配套的精品资源，点击获取