数字信号处理实验报告全攻略：从采样混叠到滤波器设计

简介：这是一份面向数字信号处理课程学习者的DSP实验配套资料，覆盖循环操作、双操作数乘法、并行运算、小数/长字/浮点运算以及卷积、相关运算等核心实验，适合需要对照指导书完成上机实验、撰写实验报告或复习DSP基础操作的高校学生。压缩包共21个文件，以docx实验报告、doc指导书为主体，另含C源码、工程配置及编译输出文件等，约8.12MB，可满足从阅读实验说明到查看实际工程结果的多层次需求。已有599人学习下载，内容包含分步骤的详细实验报告，也提供dsp实验指导书和可参考的工程文件，便于理解实验原理、操作流程及结果分析。从循环结构到浮点运算，逐步覆盖DSP常见计算场景；卷积与相关实验还给出直接法与快速法对比思路，有助于理解不同实现方式的适用条件，对系统掌握数字信号处理基本运算有较好的参考价值。

1. DSP 实验报告包到底是什么：拆开 RAR 之前你需要知道的

这个压缩包的名字非常直白：一份数字信号处理课程的实验指导书，加上一份最终定稿的实验报告。它的目标用户通常是大三或者研一第一次接触 DSP 实验的学生，也可能是补做实验、重修或者整理归档的人。你下载它的原因大概率只有两个：不知道怎么开始写实验报告，或者想知道一份能过的报告到底长什么样。先说结论：这类包的价值不在那几页指导书，而在报告里数据、图表和结论的对应关系——能不能让实验步骤被复核、结果被复现、评审挑不出毛病。

接下来的拆解逻辑按一份常见实验包来走：指导书怎么读、实验怎么跑、报告怎么写、哪些坎容易翻车。即便你手里的包来自不同的课程平台，方法论是通用的。

2. 先读指导书再动手：把实验包拆成一份施工清单

2.1 解压后的文件清单怎么看：指导书、报告模板和可能的参考工程

拿到压缩包，第一步不是双击报告文件，而是先建立文件清单。常见命名规律是“实验指导书.pdf”“实验报告模板.doc(x)”，偶尔还有“参考代码”或“仿真工程”子文件夹。如果包只有两个文档，那说明是纯文档版，实验得靠自己在 MATLAB、Simulink 或其它环境里搭。如果还夹着 .m 文件、.slx 模型文件或 CCS 工程文件，参考价值会翻倍——至少能看到老师期望的仿真链路长什么样。

我一般会把指导书目录抄一遍，按实验编号列一张表。每个实验记录三件事：实验目的里的核心概念、实验步骤里的关键参数（采样率、滤波器阶数、截止频率这类）、以及最后要求提交的图表清单。这张表是后面的施工图。别嫌这一步慢，后面所有“不知道写什么”的卡壳，几乎都能追溯到这一步没做。

常见实验可按下表归类，拿到指导书后再逐项填充。这张表填完，你就清楚这个包到底覆盖了哪些实验，以及报告里最少要准备多少张图。填表过程中如果发现某个参数指导书里没给数值，比如只写“选择合适采样率”，这就是需要你自己补的实验条件——补多少，直接决定后面实验能不能跑通。

2.2 数字信号处理实验绕不开的四类板块：从采样混叠到滤波器设计

不同学校的 DSP 实验各有侧重，但翻来覆去逃不出四块：信号时域分析与采样混叠、离散傅里叶变换与频谱分析、数字滤波器设计（FIR 与 IIR）、卷积与系统响应。这是课程体系决定的。

采样混叠实验的核心是 fs 与 f0 的关系，判断标准是采样定理：fs 必须大于两倍最高频率，否则高频信号会折返进低频区域。FFT 实验的核心是点数 N、窗函数和频率分辨率 Δf = fs/N。这里有两个容易被忽略的连带问题：N 太小时谱线栅格粗，峰值可能偏离真实频率；加窗后主瓣变宽，两个频率相近的信号可能合成一个峰。滤波器设计里 FIR 重点看窗函数法和等波纹法，IIR 重点看巴特沃斯和切比雪夫的双线性变换，核心指标是通带波纹、阻带衰减和过渡带宽。卷积与系统响应实验相对简单，验证线性时不变系统的性质，但报告里要区分线性卷积和圆周卷积的区别，这是最常藏思考题的地方。

这四个板块不是孤立的：采样混叠是时域与频域的第一道连接，FFT 是观察频域的放大镜，滤波器是频域设计的工程落地，卷积则是时域分析的底层工具。写报告时如果能把这条主线理出来，复盘提问时能把前后实验串起来，远比孤立背诵每个实验目的更有说服力。每一块都有固定套路，指导书给的参数是出发点而不是最终答案。看到“取采样率为 8000Hz”照做可以，但你要明白为什么取这个值——语音信号带宽通常在 4kHz 以内，8kHz 采样恰好留足奈奎斯特余量。若被问“为什么不是 16kHz”，你答“为了和后续实验保持一致、降低计算量”，和答“不知道，指导书写的”，得分完全不同。

2.3 报告骨架：实验目的、原理、步骤、结果、分析缺一不可

一份能过审的实验报告，结构上其实非常死板。我见过太多人把“实验原理”一节写成教材目录，把“结果分析”写成“图一为…图二为…”两句话带过，这就是丢分点。按下面这个骨架走，至少不跑偏。

实验目的写两到三句话，说清楚这一个实验验证了什么概念，比如“验证采样定理中混叠产生的条件，观察不同采样率下同一信号的频谱变化”。实验原理推到可算的地步：关键等式必须具备，比如 FIR 滤波器的差分方程 y[n] = ∑ h[k]x[n-k]，或者 IIR 的传递函数表达式；大段推导可以省，但要保证从公式能推出实验里的参数。实验步骤按指导书顺序写，每一步注明在软件里做了什么操作，配上截图。实验结果放波形、频谱、幅频响应图，坐标轴必须标单位。结果分析部分，每张图回答三个问题：现象是什么、为什么出现、参数变了会怎样。

结论段把每个实验的核心发现压缩成一两句话，直接回应指导书的思考题。这里有个习惯值得培养：结论里只写你从数据里读出来的东西，不写教材原话。比如“当 fs=2kHz 时，3kHz 信号混叠为 1kHz，实测波形周期约 1ms，与计算值一致”，这是合格的结论；“本实验验证了采样定理”这十个字是没用的废话。

最后提醒一点：实验报告往往有页眉、页脚和字体要求，指导书里通常写了但容易被忽略。把“格式要求”单独列为施工清单最后一行，写完报告后逐项对照，不然可能因字体不对被退回重交——这种冤枉事在本科阶段遇过不止一次。

3. 从指导书到最终报告：每个实验的代码、参数和结果写法

3.1 第一个实验的最小命令：用 MATLAB 验证混叠现象

直接看代码。以下是最常见的混叠验证脚本，参数按语音实验常用的 8kHz 采样率来设。

% 混叠验证实验 fs = 8000; % 采样率 8kHz，语音常用 t = 0:1/fs:0.05; % 50ms 时长，正好 400 个采样点 f0 = 3000; % 原始信号频率 3kHz x = sin(2*pi*f0*t); % 3kHz 正弦波 plot(t, x); xlabel('时间/s'); ylabel('幅度'); title('fs=8kHz, f0=3kHz, 原始波形'); % 降采样到 2kHz 观察混叠 fs2 = 2000; t2 = 0:1/fs2:0.05; x2 = sin(2*pi*f0*t2); % 用 2kHz 采样率直接采样 3kHz 信号 figure; plot(t2, x2); xlabel('时间/s'); ylabel('幅度'); title('fs=2kHz 时直接采样 3kHz 信号，观察混叠');

这段代码的逻辑是：用足够高的采样率先把原始信号画出来，再用低采样率对同样频率的信号采样，观察时域波形是否“变了样”。3kHz 对 2kHz 采样率来说超过了奈奎斯特频率（fs/2 = 1kHz），所以采样后的频率会折返到 1kHz 上，你看到波形周期明显变长就是混叠的直接证据。

参数说明：时长取 50ms，能覆盖至少 150 个完整周期，如果时长太短波形稀疏，混叠效果不明显；f0 和 fs2 的比值决定折返频率，改成 3.5kHz 会折到 0.5kHz，波形会变得非常慢。跑完这个实验，报告要截两张图：一张高采样率下的原始波形，一张低采样率下的混叠波形。分析段写一句话加一个解释：观察到的混叠频率等于 |f0 - n*fs2| 落在 [0, fs2/2] 区间内的那个值，把 n 找出来代入计算。这个计算过程是评分时最想看到的，光贴图不推导等于白做。

3.2 FFT 频谱分析的参数怎么设：点数、窗函数和坐标轴陷阱

FFT 实验里最容易翻车的地方不是 FFT 本身，而是你画出频谱后坐标轴没标对、谱线没找对。先看最小可运行的代码。

% FFT 频谱分析实验 fs = 8000; N = 1024; % FFT 点数 t = (0:N-1)/fs; % 时长 0.128s f1 = 1000; f2 = 1500; % 双频信号 x = sin(2*pi*f1*t) + 0.5*sin(2*pi*f2*t); X = fft(x, N); f_axis = (0:N-1)/N*fs; % 完整频率轴，0 到 fs % 单边谱处理：只取前 N/2 点，幅值乘以 2 X_single = X(1:N/2); f_single = f_axis(1:N/2); mag = 2*abs(X_single)/N; plot(f_single, 20*log10(mag + eps)); xlabel('频率/Hz'); ylabel('幅值/dB'); title('双频信号的幅度谱（1024点 FFT）'); grid on;

逻辑说明：fft(x, N) 算的是 N 点离散傅里叶变换，输出是复数序列，包含正负频率。幅值谱的物理意义是每个频率分量的幅度，除以 N 是归一化，乘以 2 是因为单边谱只取一半。参数说明：频率分辨率 Δf = fs/N = 8000/1024 ≈ 7.8Hz，这意味着两个频率相差小于 7.8Hz 时谱峰可能分不开；如果要分辨 5Hz 以内的频率差，就得把 N 加到 4096 甚至更大。加窗的问题在于，直接对截断的时域信号做 FFT，矩形窗的旁瓣会污染频谱——如果看到 1kHz 谱峰旁边出现一串逐渐衰减的“裙边”，那不是信号本身有这些频率，而是泄漏，换成汉明窗能压下去一部分，但主瓣会变宽。

报告里必须画两版图：一版是直接用矩形窗（不额外加窗）的结果，另一版是加汉明窗后的结果，然后对比说明窗函数对频谱泄漏的影响。坐标轴单位一个用 Hz，一个用 dB 或线性幅度，不要混着用。FFT 参数的设置在没搞懂之前都像是玄学，其实就是 N、fs 和窗函数三个变量在起作用。

注意：FFT 结果图的横轴单位是 Hz，不是示波器上的“格”。如果指导书要求的是归一化频率，换算公式是 f_norm = f / (fs/2)，报告里两种轴都要标注。

3.3 FIR 滤波器实验：设计指标、差分方程和幅频曲线一个都不能少

滤波器实验常常是课程考核的重头戏，因为参数多、可调空间大，很容易看出学生到底懂不懂。下面给一个窗函数法设计 FIR 低通滤波器的完整例子。

% FIR 低通滤波器设计（汉明窗） fs = 8000; fc = 1500; % 截止频率 1.5kHz N = 33; % 滤波器长度 33 h = fir1(N-1, fc/(fs/2), 'low', hamming(N)); freqz(h, 1, 1024, fs); % 画幅频和相频响应 % 滤波前后对比 t = (0:1000)/fs; x = sin(2*pi*800*t) + sin(2*pi*3000*t); % 800Hz 通过，3kHz 被滤除 y = filter(h, 1, x); figure; subplot(2,1,1); plot(t, x); title('滤波前（含 800Hz 和 3kHz）'); subplot(2,1,2); plot(t, y); title('滤波后（保留 800Hz）');

逻辑说明：fir1 第一个参数是滤波器长度减一，第二个参数是归一化截止频率，即 fc 除以奈奎斯特频率（fs/2），所以 1500/4000 = 0.375。通带范围由这个归一化值决定，阻带衰减依赖窗函数和阶数。参数怎么调：阶数 N 越大过渡带越窄，但计算复杂度上升，相位延迟变大；汉明窗阻带衰减大约 53dB，如果实验要求 60dB，就得换凯泽窗，把 beta 参数调到 6 以上。注意 filter 输出会有群延迟 (N-1)/2 个采样点，对比滤波前后的时域波形时，必须把输入信号也做同样延迟，否则图形上会看到明显错位。

报告部分要把设计指标写成表：通带截止频率、阻带截止频率、通带波纹、阻带衰减，每个指标对应实际算出来的值。然后贴幅频响应曲线，在图上把 -3dB 点标出来，说明实际截止频率和设计值是否吻合。这个“标点”的动作特别加分，因为大多数学生的图里只有一条光溜溜的曲线。

4. 报告翻车避坑：五个高频问题的现象、原因和解决

这一章按实验报告里最常见的问题来写，每条都按“现象→原因→解决”三步走，写报告前对照一遍，能救回不少冤枉分。

4.1 坑一：频谱图左右完全对称，你却没做单边谱处理

现象：FFT 结果画出来，频谱图以 fs/2 为轴左右对称，两个相同高度的峰，频率轴从 0 到 fs。报告交了，被批“频谱分析不规范”。

原因：fft 输出包含正负频率，实信号的频谱天然共轭对称，直接画全谱就是两边都有。一半信息是冗余的，而且用线性幅度坐标时，负频率部分的峰会让主峰幅度看起来只有真实值的一半。

解决：用单边谱，只取前 N/2 点，幅值乘 2（直流分量不乘）。前面代码里已经写了，关键是报告里的图也要成对出现——一版说明完整谱为什么对称，一版展示单边谱的正确读法。如果你用的是 Python，numpy.fft.rfft 直接返回单边谱，能省一步。

4.2 坑二：滤波器实际截止频率和设计值对不上

现象：设计一个 fc=1500Hz 的低通滤波器，freqz 图里 -3dB 点在 1400Hz 左右，你觉得是程序错了，硬把数据改成 1500Hz。

原因：窗函数设计法里，实际截止频率取决于窗的主瓣宽度，窗函数会改变滤波器过渡带位置，尤其是阶数较短（比如 N=15）时偏差更大。这不是 bug，是算法特性。

解决：把 freqz 结果放大，用光标读取 -3dB 对应频率，在报告里写“实测截止频率约 1430Hz，与设计值 1500Hz 偏差约 4.7%，原因是汉明窗主瓣有限宽度导致的过渡带偏移”。这样写就是分析，不是纠错。千万别把代码里的 fc 直接抄到报告里当实测值，一旦被问细节就穿帮。

4.3 坑三：Simulink 仿真图里信号根本没混叠，你却写了混叠分析

现象：有的同学按指导书在 Simulink 里搭了一个采样加示波器模型，示波器显示波形看起来“还挺平滑”，报告里写“该信号无混叠”。但理论计算明明应该混叠。

原因：示波器的显示采样率太高了。Simulink 示波器默认用变步长求解器内部时间步长去显示波形，把连续信号显示得很细腻，采样器输出的阶梯波形被平滑掉，肉眼看不出混叠后的包络。

解决：在示波器设置里把 Decimation 或采集间隔调成和采样率一致，每 N 个点取一个显示；更稳妥的做法是把采样器输出接一个零阶保持器再进示波器。报告里如果用了 Simulink 截图，一定要检查示波器左侧时间轴刻度，确认显示步长等于 1/fs，否则截图里的波形和计算值对不上。

4.4 坑四：滤波前后的时域波形错位，你说“滤波后信号变小了”

现象：filter 滤波后，用 plot(t, y) 直接把波形画在原始输入同一张图上，发现 y 比 x 整体向右平移了多个采样点，起始处还有一段幅度很小。你以为是滤波器衰减太大。

原因：FIR 滤波器是因果系统，输出相对输入有固定群延迟 (N-1)/2 个采样点。滤波器初始状态为零，所以输出从零开始爬升，这些都不是衰减造成的。

解决：画对比图时用 delay = (N-1)/2，把滤波输出往前移 delay 个点，或把输入信号也往后移 delay 个点。分析段写“系统群延迟为 16 个采样点，对应 2ms，符合预期”，直接证明你理解因果性和滤波器的群延迟特性。

delay = (N-1)/2; figure; plot(t, x, 'b'); hold on; plot(t(1:end-delay), y(delay+1:end), 'r'); legend('原始信号', '滤波后（对齐群延迟）');

4.5 坑五：图表编号混乱，参数表和图片对不上

现象：报告里图 3 引用的是图 5 的内容，或者代码里 fs=8000，报告里写“采样率 10kHz”。

原因：写报告时复制了之前的模板，改图但没改文字，数字和图表脱节。这是典型的复制粘贴事故。

解决：写完初稿后，用脚本提取报告里所有图片标题、坐标轴标注和数值，与参数表逐项核对。最简单的方法是把实验代码里的参数都集中在脚本开头，报告里的数值全部从运行结果复制，不要手动敲。手动敲数字是血泪教训重灾区——一不小心把 3kHz 敲成 3MHz，评审一眼就能看出你没跑过实验。

5. 让报告从“能交”变成“能过审”：平台适配、数据加工和思考题答法

5.1 平台差异：MATLAB、Simulink、CCS 三种实验环境的报告侧重点

不同学校 DSP 实验平台不一样，报告侧重点完全不同。拿到指导书后第一件事是确认平台。

MATLAB 为主的课程，报告重点是代码逻辑和频谱图；Simulink 为主的课程，报告重点是模型结构图和每个模块的参数设置，截图必须把模块参数对话框打开再截，光截模型图看不出设了什么；CCS 这类面向开发板的课程，重点则是代码、内存分配和实际运行时波形。如果指导书里出现“新建工程”“链接命令文件”“查看寄存器”这些词，说明是嵌入式方向的实验，报告里要有工程配置截图，并单独留一个调试章节，把连接失败、数据溢出这类过程写进去。这类报告最值钱的部分是你怎样发现问题、怎样定位到程序哪一行——只贴成功截图、不写调试过程的“幸运报告”得分通常不高。

5.2 实验数据怎么二次加工：用表格呈现参数对比，用局部放大图强调细节

我见过的好报告都有共性：不只贴原始图表，还把关键信息提炼成表。比如设计滤波器后，把设计指标、实测值和理论值放在同一张表里，三列一一对应，偏差写在第四列。看到这种表，评审基本不会再怀疑你对数据的理解。

图表处理上，频谱图一定要有局部放大。原始频谱图范围太长，细节都在某个窄频段里，不放大不一定能看到。用 xlim([900, 1100]) 截一段，把主瓣宽度和旁瓣幅度标出来，再在分析段写“主瓣宽度约 23Hz，第一旁瓣比主瓣低约 55dB”。这种定量描述是拉满分的动作。

5.3 思考题的加分答法：不止回答，还要反推验证

指导书最后的思考题别空着，也别只写一行“不会”。这类题常见两种：一是理论验证型，如“什么条件下混叠可以避免”，你就代入公式算一个具体边界值，再把实验观察到的现象用公式解释一遍；二是方案对比型，如“比较矩形窗和汉明窗的优劣”，你就按阻带衰减、过渡带宽和主瓣宽度三个指标列表对比，用实验数据引证。

更进阶的做法是：在思考题答案末尾加一行“如果在实际系统中要用这个滤波器，我会选 XX，原因是 YY”——提前用工程判断收尾。这一句话能让你从“完成了作业”变成“具备初步工程思维”。多写几个点没坏处，只要每个点都能对应实验数据，就不是注水。

6. 提交前自查清单：五个关键核对点和一个命名习惯

这里不写大道理，只给你一张提交前必过的自查表。

第一，检查压缩包命名是否符合课程要求的格式，例如“学号_姓名_实验报告.rar”。有的课程要求压缩包内不放指导书原文件，只放报告和代码，看清楚再打包，顺序反了会被直接退回。

第二，报告里的图都标了图号、图题、坐标轴单位，没有一张是黑白模糊的截图。尤其检查示波器截图，如果波形边缘出现锯齿，说明分辨率不足，重新截。

第三，所有数值从代码运行结果复制，报告里没有手动敲的 magic number。如果有手动填写的参数，在备注里标“由运行结果读取”，并保留对应图。

第四，实验结果分析和图里的现象互相印证。挑三张核心图，对着图说一句“为什么会这样”，说不出来就回去补分析。

第五，如果参考了别人的代码或开源 demo，在报告末尾注明出处并说明做了哪些修改。这是学术规范问题，不是可有可无。

另外说一个我的习惯：从第一次跑通实验开始，把每一版截图命名为“fig3_fs8k_f0_3k_混叠.png”这种带参数的名字，不要叫“截图1”。写报告时直接从图库按文件名挑图，永远不会有“哪张图对应哪个实验”的混乱。我已经因为图名太抽象吃过亏——交报告前发现图库里的“新建位图”根本分不清是哪组参数跑的，只能重跑一遍全部实验，整个周末搭进去。

希望帮到你：把这套流程走一遍，指导书读透、实验跑真、数据写清，剩下的交给评审，问题不大。

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