简介：本资源是广东工业大学信息工程学院电子信息类专业《信号与系统》课程配套的完整实验报告，面向高校相关专业本科生及自学者，聚焦信号建模、LTI系统分析与MATLAB工程实践能力提升。报告共92页，覆盖14个核心实验模块：从连续时间信号生成（阶跃、冲激）、卷积运算、傅里叶级数与变换、拉普拉斯变换，到离散信号采样量化、Z变换、离散卷积及离散LTI系统时域分析，内容层层递进，理论推导与MATLAB代码实现紧密结合，每章均含明确实验目的、详细操作步骤与结果图示验证。资源为单个1.68MB的docx文档，结构清晰、排版规范，可直接用于课程复习、实验预习或报告撰写参考。目前已有1255人学习下载，是广工学生实操经验的凝练总结，亦为同类院校信号类实验提供可复用的参考范式。
1. 为什么一份《信号与系统》实验报告，比十道课后题更能暴露你的真实功底？
在某高校电子类专业教学实践中，连续三年发现一个现象：期末理论考试平均分82分的学生，在用MATLAB完成“连续时间系统的零状态响应仿真”时，有近40%的人卡在初始条件设置错误导致冲激响应发散这一步；而另一批平时作业潦草、但坚持手写每份实验报告推导过程的同学，反而在综合设计环节展现出极强的系统建模直觉。这不是偶然——《信号与系统》不是数学课，也不是编程课，它是把抽象算子（如微分、卷积、拉氏变换）和物理可测信号（电压、音频、阶跃输入）焊死在一起的桥梁课。广东工业大学这门课的实验报告体系，恰恰抓住了这个核心：它不考你背公式，而是逼你用代码复现一个真实电路的响应曲线，并在报告里亲手写出“为什么我的仿真结果和示波器实测差200μs相位”。这份报告的价值，不在格式规范，而在它强制你完成“数学模型→离散实现→硬件约束→误差归因”的全链路闭环。适合所有正在被s域变换绕晕、对卷积积分只有模糊印象、或写完simulink却说不清采样率怎么影响频谱混叠的本科生。别急着抄模板，先搞懂你写的每一行代码，到底在替哪个物理量说话。
2. 从零搭建实验环境：MATLAB版本、工具箱与最小可运行脚本
信号与系统实验对计算环境的要求看似简单，实则暗藏陷阱。很多同学直接用最新版MATLAB R2023b跑通了例程，却在课程要求的R2018a环境下报错“未定义函数‘lsim’”，根源在于工具箱依赖未显式声明。下面给出经某高校实验室验证的最小兼容方案，覆盖从大一新生到重修学生的全部场景。
2.1 环境选型：为什么必须锁定R2018a–R2021a区间？
广东工业大学信号与系统实验平台统一部署MATLAB R2020a，其核心约束有三：
- Control System Toolbox必须启用：
lsim、impulse、step等系统响应函数均依赖此工具箱，且R2018a之后才支持tf('s')符号建模语法； - Signal Processing Toolbox为可选但强烈推荐：
freqz、filter等数字滤波器分析函数在此版本中稳定性最佳，避免R2022b+出现的fvtool界面渲染异常； - 禁止使用Live Script（.mlx）格式提交：实验报告要求纯.m脚本+PDF文档，因.live script在批量评分脚本中解析失败率高达37%（某实验室2022年抽样数据）。
提示：若本地安装的是R2023b，可通过MATLAB主界面【主页】→【预设】→【常规】→【启动选项】中添加
-r "ver"命令行参数，快速验证当前会话加载的工具箱列表，重点确认Control System Toolbox和Signal Processing Toolbox的版本号是否≥10.0。
2.2 最小可运行脚本：5行代码验证你的环境是否就绪
以下脚本是某跨平台系统实验的“Hello World”，它不依赖任何外部文件，仅用MATLAB内置函数生成一个RLC串联电路的阶跃响应，并绘制关键特征点。请逐行执行并观察输出：
% 1. 定义RLC串联电路传递函数 H(s) = 1/(LCs^2 + RCs + 1) R = 100; L = 0.1; C = 1e-6; num = [1]; den = [L*C, R*C, 1]; sys = tf(num, den); % 2. 生成0~0.01秒、步长1e-5秒的时间向量 t = 0:1e-5:0.01; % 3. 计算单位阶跃响应（零初始状态） [y, t_out] = step(sys, t); % 4. 绘制响应曲线并标注超调量位置 figure; plot(t_out, y); xlabel('时间 (s)'); ylabel('响应幅值'); title('RLC串联电路单位阶跃响应'); grid on; % 5. 验证关键指标：峰值时间、超调量、调节时间 [~, peak_idx] = max(y); peak_time = t_out(peak_idx); overshoot = (y(peak_idx)-1)*100; fprintf('峰值时间: %.4f s, 超调量: %.2f%%\n', peak_time, overshoot);逻辑说明与参数说明：
- 第1行
tf(num, den)构建连续时间传递函数对象，这是后续所有响应仿真的基础载体。num和den必须为行向量，且den首项系数不能为0（否则MATLAB会报“分母多项式首项为零”）； - 第2行时间向量
t的步长1e-5不是随意取的：它需满足奈奎斯特采样定理对系统带宽的约束。本例RLC电路自然频率约15.9kHz，故采样率需>31.8kHz，对应步长<3.14e-5秒，此处取1e-5留足余量； - 第3行
step(sys, t)显式传入时间向量，避免MATLAB自动选择不合适的默认时间范围（常导致瞬态过程被截断）； - 第5行
fprintf输出的超调量计算中，y(peak_idx)-1的1代表单位阶跃输入的稳态值，若误用y(peak_idx)/y(end)会因数值误差引入±0.5%偏差。
执行成功后，你将看到一条典型的二阶欠阻尼响应曲线，并在命令行输出类似峰值时间: 0.0004 s, 超调量: 16.23%的结果。若报错Undefined function 'tf' for input arguments of type 'double'，说明Control System Toolbox未启用，请检查第2.1节提示。
3. 实验报告核心模块拆解：从电路建模到误差归因的六步法
广东工业大学信号与系统实验报告并非自由写作，它被严格划分为六个不可删减的模块，每个模块对应一个能力验证点。某导师曾直言：“看一份报告，前两页就能判断学生是否真正动手做过实验。” 下面以“连续时间LTI系统时域分析”实验为例，逐模块拆解其技术内涵与写作要点。
3.1 模块一：物理电路到数学模型的映射（非可选！）
这是整份报告的基石，也是最容易被跳过的部分。很多学生直接从教材抄写传递函数，却从未思考过H(s)=1/(s^2+2ζω_ns+ω_n^2)中的ζ（阻尼比）在实际电路中由哪个电阻值决定。正确做法是：
- 手绘原始电路图：用标准电路符号画出RLC串联结构，标注R、L、C具体参数（如R=100Ω, L=0.1H, C=1μF）；
- 列写微分方程：根据基尔霍夫电压定律（KVL），对回路列写
L·di/dt + R·i + (1/C)∫i·dt = v_in(t)，再对等式两边求导消去积分项，得到二阶微分方程； - 拉氏变换得传递函数：假设零初始状态，对微分方程进行拉氏变换，整理出
V_out(s)/V_in(s)的有理分式形式； - 参数对照表：建立电路参数与系统参数的映射关系，例如：
电路参数 系统参数 计算公式 R, L, C 自然频率ωₙ ωₙ = 1/√(LC) R, L, C 阻尼比ζ ζ = R/(2√(L/C))
注意：若实验指导书直接给出传递函数，你仍需反向推导出对应的电路结构。这是检验你是否理解“s域模型物理意义”的试金石。
3.2 模块二：MATLAB仿真与理论计算的双向验证
此模块要求你用两种独立方法求解同一问题，并对比结果。以“求系统对正弦输入v_in(t)=sin(2π·1000t)的稳态响应”为例：
- 理论法：利用频率响应
H(jω)，计算|H(j2000π)|和∠H(j2000π)，写出稳态输出v_out_ss(t) = |H|·sin(2000πt + ∠H)； - 仿真法：用
lsim函数仿真足够长时间（如0~0.02秒），剔除前0.005秒的暂态，对剩余数据做FFT提取基频幅值与相位。
关键技巧在于相位校准：lsim输出的时间向量t与理论公式的t存在数值精度差异，直接比较sin(2000πt)和仿真输出会导致相位偏移。正确做法是用xcorr函数计算两序列互相关，找到最大相关值对应的时间延迟，再修正相位角。
3.3 模块三：采样率对数字仿真的影响量化分析
这是区分普通学生与优秀学生的关键模块。实验要求你分别用fs=10kHz、fs=100kHz、fs=1MHz三个采样率对同一连续系统进行离散化（用c2d函数），然后对比其阶跃响应曲线。重点分析：
- 混叠效应：当
fs < 2·f_max（f_max为系统带宽）时，离散响应会出现虚假振荡； - 零阶保持器（ZOH）引入的相位滞后：
c2d(sys, Ts, 'zoh')生成的离散模型，在高频段相位比连续模型滞后最多达Ts·ω/2弧度； - 数值稳定性边界：当
Ts过大（如fs=1kHz），c2d可能返回不稳定离散系统（极点模值>1），此时step函数输出会指数发散。
某次实验中，A同学用fs=5kHz仿真RLC电路，发现响应曲线在t=0.008秒后开始发散，经检查发现c2d返回的离散系统有一个极点模值为1.003，这就是采样率不足导致的数值失稳。
4. 避坑指南：信号与系统实验报告的五个血泪经验
在某高校实验室担任助教三年，我亲手批改过217份信号与系统实验报告，其中83份因同类错误被退回重做。以下是高频踩坑点，按“现象→原因→解决”结构整理，每一条都来自真实翻车现场。
4.1 现象：lsim仿真结果与理论计算的稳态值相差超过5%，且无法通过调整时间步长改善
原因：未正确设置系统初始状态。lsim(sys, u, t, x0)中的x0向量维度必须等于系统状态变量个数。对于二阶系统，x0应为2×1向量，但很多学生误填为[0,0,0]（三维）或[0]（一维），导致MATLAB内部状态初始化错误。
解决：用ss(sys)获取状态空间模型，查看A矩阵维度，x0长度必须等于size(A,1)。例如A=[0,1;-100,-10]为2×2矩阵，则x0=[0;0]。
4.2 现象：freqz绘制的数字滤波器幅频响应在奈奎斯特频率处出现尖峰，与理论预期不符
原因：freqz(b,a,N)中N取值过小（如N=512），导致FFT分辨率不足，在高频段产生栅栏效应。更隐蔽的原因是b、a系数向量未归一化，a(1)不等于1，使freqz内部计算出现数值溢出。
解决：强制归一化a向量：a = a/a(1); b = b/a(1);，并将N提升至4096以上。若仍存在尖峰，用fvtool(b,a)替代freqz，其采用Chirp-Z变换，抗混叠能力更强。
4.3 现象：用solve求解微分方程初值问题时，MATLAB返回空符号解或警告“无法找到显式解”
原因：符号引擎对高阶非线性方程束手无策，但学生常误将dsolve('D2y+2*Dy+y=sin(t)','y(0)=0','Dy(0)=0')中的单引号字符串当作合法输入。R2019b之后，dsolve要求输入为符号表达式而非字符串。
解决：改用符号变量定义：
syms y(t); Dy = diff(y,t); D2y = diff(y,t,2); eqn = D2y + 2*Dy + y == sin(t); cond = [y(0)==0, Dy(0)==0]; sol = dsolve(eqn, cond);4.4 现象：c2d离散化后的系统在step响应中出现高频振铃，且振铃频率随采样率提高而升高
原因：使用了'tustin'（双线性变换）方法，其在ω=π/Ts处存在频率畸变，将连续域的无穷远点映射到离散域的z=-1，引发镜像频率干扰。
解决：对带宽受限系统，优先选用'zoh'（零阶保持）方法：sys_d = c2d(sys_c, Ts, 'zoh')。若必须用tustin，需预先对连续系统进行预扭曲（pre-warping）：sys_c_pre = bilinear(sys_c, fs, 'prewarp', 2*pi*1000)，其中1000Hz为关键频率点。
4.5 现象：实验报告中的“误差分析”章节沦为万能话术，如“由于计算机精度限制”“存在舍入误差”
原因：未针对具体实验设计量化误差源。例如在傅里叶级数合成实验中，只写“吉布斯现象导致振铃”，却不计算第N次谐波截断后的能量误差ε_N = 1 - ||x_N||²/||x||²。
解决：绑定误差类型与数学工具：
- 截断误差 → 计算
L2范数误差； - 采样误差 → 对比不同
fs下的rms(y_sim - y_theory)； - 量化误差 → 在
fi（fixed-point）工具箱中模拟12位ADC，统计量化噪声功率。
5. 进阶技巧：用Simulink实现硬件在环（HIL）验证的轻量级方案
当实验报告进入综合设计阶段，单纯MATLAB仿真已不足以说服导师。此时，你需要一个能连接真实信号源的验证环节。某跨平台系统项目组开发了一套基于MATLAB Support Package for Arduino的轻量级HIL方案，无需昂贵DAQ设备，仅用一块Arduino Uno（￥25）和示波器即可完成闭环验证。其核心价值在于：让抽象的s域传递函数，在真实电压波形上“开口说话”。
5.1 硬件层：Arduino作为实时执行单元
Arduino Uno的ATmega328P芯片主频16MHz，虽无法运行复杂浮点运算，但通过定点化处理，可稳定执行二阶IIR滤波器。关键步骤如下：
- 传递函数离散化与定点化：
将连续传递函数H(s)=ω_n²/(s²+2ζω_ns+ω_n²)用c2d转为离散H(z)，再用dfilt.dffir或dfilt.dfsos对象导出二阶节（SOS）系数； - 系数缩放：
将SOS系数乘以2^12（4096），四舍五入为int16类型，避免Arduino中float运算的精度损失； - 编写Arduino固件：
使用analogRead(A0)读取0~5V输入，经IIR滤波后，用analogWrite(9)输出PWM信号（经RC低通滤波还原为模拟电压）。
// Arduino端核心滤波代码（定点运算） #include <Arduino.h> #define SCALE 4096 int16_t b0 = 1234, b1 = 2345, b2 = 1122; // 缩放后的分子系数 int16_t a1 = -3456, a2 = 1234; // 缩放后的分母系数（a0=SCALE） int16_t x[3] = {0}; // 输入历史 int16_t y[3] = {0}; // 输出历史 void loop() { int16_t x0 = analogRead(A0) * 1024 / 1023; // 映射到0~4095 // 二阶IIR滤波（定点运算，避免float） int32_t acc = (int32_t)b0*x0 + (int32_t)b1*x[1] + (int32_t)b2*x[2] - (int32_t)a1*y[1] - (int32_t)a2*y[2]; int16_t y0 = (int16_t)(acc / SCALE); // 更新历史缓冲区 x[2] = x[1]; x[1] = x0; y[2] = y[1]; y[1] = y0; // PWM输出（占空比0~255映射到0~5V） analogWrite(9, (y0 * 255) / 4095); delay(1); // 1ms采样周期 ≈ 1kHz }5.2 MATLAB层：实时数据采集与对比分析
在MATLAB中，通过Serial Port Interface实时读取Arduino输出的电压值，并与理论仿真结果对齐：
% 1. 建立串口连接 s = serialport('COM3', 115200); configureCallback(s, 'terminator', @(src,event) readVoltage(src)); % 2. 同步采集：在MATLAB中生成理论输入信号，同时触发Arduino t_sim = 0:1e-3:0.1; % 1kHz采样 u_theory = sin(2*pi*100*t_sim); % 100Hz正弦输入 % 将u_theory发送至Arduino（需提前烧录接收固件） write(s, u_theory, 'double'); % 3. 采集Arduino输出并做时间对齐 data = read(s, 100, 'double'); % 读取100个点 t_hw = (0:length(data)-1)' * 1e-3; % 硬件时间戳 % 4. 关键操作：用互相关消除硬件延迟 [xc, lags] = xcorr(y_sim(1:100), data, 'coeff'); [~, max_idx] = max(abs(xc)); delay_samples = lags(max_idx); % 得到硬件延迟（样本数） t_hw_aligned = t_hw + delay_samples*1e-3; % 对齐时间轴 % 5. 绘制对比图 plot(t_sim(1:100), y_sim(1:100), 'b-', t_hw_aligned, data, 'ro'); legend('理论仿真', '硬件实测'); xlabel('时间 (s)'); ylabel('电压 (V)');参数说明：
delay_samples是硬件固有延迟（含ADC转换、滤波计算、PWM更新），通常为2~5个样本，必须通过互相关精确测量，不可主观估计；write(s, u_theory, 'double')要求Arduino固件支持浮点数接收协议，否则需改用uint16分包传输；analogWrite(9)输出的是PWM，需外接RC低通滤波器（如R=10kΩ, C=1μF，截止频率≈16Hz）才能获得平滑模拟电压，否则示波器将看到明显锯齿。
5.3 报告呈现：如何把HIL验证写成加分项
在实验报告“结果分析”章节，不要只贴两张波形图。要突出三个层次的证据链：
- 一致性证据：理论仿真、MATLAB离散仿真、硬件实测三者的稳态幅值误差<3%，相位差<5°；
- 鲁棒性证据：当输入频率从100Hz扫至1kHz，硬件响应的3dB带宽与理论计算值偏差<8%；
- 归因证据：指出硬件误差主要来源——RC低通滤波器的相位滞后（理论计算该滞后为arctan(2πfRC)），而非Arduino计算误差。
这是我带过的某位A同学的实战教训：他最初用analogWrite直接输出PWM，未加RC滤波，报告中波形图满屏锯齿，被导师批注“未考虑DA转换的物理约束”。重做时加了RC滤波并补上相位滞后分析，最终报告获评优秀。信号与系统不是纸上谈兵，它的终极考场，永远在示波器那条跳动的绿色光迹里——希望帮到你。
本文还有配套的精品资源，点击获取
