简介《信号与系统》实验报告是广东工业大学信息工程学院电子信息类学生的课程实践总结面向正在修读该课程、需要完成MATLAB仿真实验或整理实验文档的本科生。压缩包内为1个docx文档约1.68MB按实验一至实验四组织涵盖连续时间信号在MATLAB中的表示、卷积运算、连续与离散LTI系统的时域分析、周期信号的傅里叶级数及频谱分析、傅里叶变换性质、拉普拉斯变换与Z变换等内容。每章列出实验目的与实验内容配有MATLAB代码和结果波形图便于对照理论验证计算过程。文档已排版成册共90余页可作为实验报告的写作框架与参考范例帮助读者理解基本信号生成、卷积求解、采样量化影响及系统频率特性分析等关键环节节省从零搭建实验步骤的时间。目前已有1257人学习下载适合作为课程实验的辅助材料与复习提纲。1. 一份 92 页的《信号与系统》MATLAB 实验报告从画波形到算卷积的完整链路第一次翻这份《信号与系统》MATLAB 实验报告是在帮某高校电子信息类专业的同学定位一段画不出曲线的代码。报告从连续时间信号的表示一路排到 Z 域分析四组实验、十四个章节每节都是「公式、代码、波形、结论」四件套。真正值钱的地方不是答案而是它把课堂上那些抽象符号压成了能敲进命令行的东西阶跃怎么定义、卷积怎么近似、傅里叶级数频谱怎么画、零极点图怎么看。它适合三类人正在赶实验报告但想弄明白每一步为什么这么写的人、需要一份信号处理命令速查模板的人、以及想把时域频域这套东西重新捡起来的人。下面按我自己的复现顺序拆。2. 信号表示的基础把公式翻译成时间轴与匿名函数MATLAB 里没有真正意义上的连续信号plot画出来的永远是折线。这句话听起来像废话但它决定了后面每一个参数的取法。这一章先把「时间轴怎么切」「阶跃怎么定义」「复信号怎么拆四个子图」这三件事固定下来后面的时移、卷积、频谱全都建立在这三块砖上。2.1 时间步长的取法dt 不是随便填个 0.01报告里几乎每段代码都以t -2:0.01:4开头。0.01 这个数字看着随意实际上和信号最高频率绑在一起。以报告里的exp(-t).*sin(10*pi*t)为例角频率 10π频率 5 Hz一个周期 0.2 秒步长 0.01 意味着一个周期采 20 个点折线看起来已经足够光滑。换成sin(100*pi*t)50 Hz一个周期只剩 2 个点画出来是锯齿甚至因为混叠伪装成一个低频信号。% 时间轴三要素起点、步长、终点 dt 0.001; % 步长先按最高频率的 1/20 取画不平滑再减半 fmax 50; % 信号里最高的正弦分量频率单位 Hz t -2:dt:4; % 时间范围按想看到几个周期/几个转折点定 N numel(t); % 采样点数后面算频谱和卷积长度要用 fs 1/dt; % 等效采样率画频谱时必须知道这个值逻辑说明dt决定波形平滑度t的首尾决定视野两者互不影响。常见误区是把axis的范围当成时间轴范围来调——axis只裁剪显示窗口不改变参与运算的数据量。参数说明fmax这个变量很多人不写直接拍一个 0.01遇到高频信号就翻车养成先估fmax再定dt 1/(20*fmax)的习惯改起来有依据。2.2 单位阶跃与门函数匿名函数比符号工具箱稳报告里2.3.1(3)那段用了一个uCT(t)而同一页的其他代码又用(t 0)硬写风格不统一。这是个信号原作者中途换过写法。复现时最省事的方案是用匿名函数把阶跃和门函数一次定义好整个脚本复用。% 基本信号定义放在脚本最前面后面所有表达式都用它们 u (t) double(t 0); % 单位阶跃 u(t)t0 处取 1 g (t, a, b) u(t-a) - u(t-b); % 门函数区间 [a, b) 内为 1 t -2:0.01:4; y (2 - exp(-t)) .* u(t); % (2 - e^{-t})u(t) plot(t, y); axis([-2 4 -0.5 2.2]); % 顺序不能反先 plot 再 axis grid on; xlabel(t); ylabel(f(t));逻辑说明double(t 0)把逻辑数组转成 0/1 数值数组才能参与乘法。逐元素乘.*不能写成*因为u(t)返回的是和t同尺寸的数组矩阵乘会直接报维度错。参数说明门函数用u(t-a) - u(t-b)而不是ta tb好处是能直接代入时移后的变量比如g(2*t1, 0, 1)做尺度变换时不用改函数体。提示报告第 6 页2.3.1(2)的代码里有连续两行给y赋值前一行y t.*((t 0) - (t 1))是残留代码被后一行覆盖了。抄作业时留着不影响结果但调试时如果看到波形对不上先检查有没有这类覆盖行。2.3 复信号的四个子图real/imag/abs/angle 各管一段复信号f(t) 2 e^{j(π/4)t} e^{j(π/2)t}是报告里比较少见的「需要拆成四张图」的题。四个分量分别是实部、虚部、模、相角对应real、imag、abs、angle四个函数。看起来机械但angle有个必须知道的特性返回值落在 (-π, π]相位跨过 ±π 时会跳变 2π画出来是一条满是竖线的锯齿。t -5:0.01:10; y 2 exp(1j*(pi/4)*t) exp(1j*(pi/2)*t); subplot(2,2,1); plot(t, real(y)); title(实部); grid on; axis([-5 10 0.5 4.5]); subplot(2,2,2); plot(t, imag(y)); title(虚部); grid on; axis([-5 10 -2 2]); subplot(2,2,3); plot(t, abs(y)); title(模); grid on; axis([-5 10 0.5 4.5]); subplot(2,2,4); plot(t, unwrap(angle(y))); title(相角展开后); grid on;逻辑说明虚数单位用1j而不是i因为i极易被for i 1:N覆盖成循环变量一旦覆盖后面所有复数运算全错而且报错信息指向的位置往往离真正出问题的地方很远。unwrap把相位跳变补回去适合看整体趋势如果要严格对照理论值直接看angle的原始输出更实在。参数说明subplot(2,2,k)的编号是按行优先先左后右、先上后下不要按「四个图随意排」的思路填。验证这一步有个便宜办法实部理论上是2 cos(πt/4) cos(πt/2)模的理论范围是[0, 4]三个单位向量同相时最大。用一行max(abs(y))和min(abs(y))对一下比盯着图看靠谱。3. 时移、反褶与卷积三个脚本级运算的边界条件理论上的f(2-t)、f(2t1)是纸面变换写成代码只改一个字符但改错的时候波形不会报错只会安安静静地画错。这一章把报告的第三章和第四章合并处理先用函数句柄把分段信号封装好再谈反褶为什么依赖时间轴的对称性最后落到卷积的数值近似——这是整份报告里最容易被「抄个大概」抄坏的一节。3.1 分段信号封装成函数文件而不是拼字符串报告里的functl.m是整份文档中唯一一个独立函数文件用四个区间的差分阶跃拼出一条折线。这种写法的好处是可以在任意自变量上调用functl(t-1)、functl(2-t)、functl(2*t1)都是同一份代码。function y functl(t) % 分段折线信号转折点位于 t -2, -1, 0, 1, 2 % 对应区间取值-1→0斜线、1常数、2常数、1→0斜线 y (t1) .* double(t-2 t-1) ... 1 .* double(t-1 t 0) ... 2 .* double(t 0 t 1) ... (-t2) .* double(t 1 t 2); end逻辑说明每个区间用一个逻辑数组做「开关」与对应表达式逐元素相乘四段相加得到完整波形。这和报告里用(t-2)-(t-1)做差分的写法等价但区间边界写得更直白改转折点时不容易错位。参数说明函数体内所有乘法必须是.*因为传入的t是数组如果只测试单个点用functl(0.5)也能跑。调用时的关键认知变换的是自变量不是函数值。functl(2-t)里的2-t在函数内部会被当作「新的时间轴」参与分段判断所以图形整体被压缩、反褶并平移而不是把画好的图左右翻一下。这一点想通了后面五个小问其实是同一个操作。表达式自变量替换几何效果f(t-1)t → t-1右移 1f(2-t)t → 2-t关于 t1 反褶f(2t1)t → 2t1压缩 2 倍后左移 0.5f(4-t/2)t → 4-t/2展宽 2 倍关于 t4 反褶3.2 奇偶分量分解fliplr 成立的前提条件报告里求偶分量和奇分量的代码只有三行y1 fliplr(y)、ye (yy1)/2、yo (y-y1)/2。这三行在特定条件下正确换一组时间轴就错而且错得很隐蔽——因为波形看起来还是像那么回事。fliplr做的是索引反转y1(i) y(N1-i)。它要等价于时间反褶f(-t)必须同时满足两个条件时间轴关于 0 对称且 0 落在某个采样点上。报告里t -3:0.01:3一共 601 个点索引 301 对应t 0fliplr后 301 号位置还是 301 号位置的元素恰好满足。% 错误示范时间轴不对称 t -3:0.01:2.99; % 601 个点但中点是 -0.005不是 0 f double(t0 t2).*(1-t); fr fliplr(f); % 整体偏移一个采样点奇分量在 t0 处不为零 % 稳妥写法按自变量重新采样 fr interp1(t, f, -t, linear, 0); % 在 -t 处插值取值越界补 0 fe (f fr) / 2; fo (f - fr) / 2;逻辑说明interp1的第四个参数是插值方法第五个是越界外推值。用-t作为查询点等价于「把时间轴翻过来重新采一遍」只要t是等间隔的就能对上不依赖对称性。参数说明如果-t的范围超出t超出部分按第五个参数补 0这正是反褶后信号为零的物理含义。验证奇偶分解是否正确看两个数就够max(abs(fo))应该等于max(abs(f))/2附近且奇分量在t0处的值必须为 0。这两条对不上说明反褶出问题了。3.3 连续卷积的数值近似conv 缺了 dt 就只是个求和报告第四章是连续时间信号的卷积运算但 MATLAB 的conv做的是离散求和。两者差一个步长因子这个因子在图上表现为幅度整体差 100 倍dt 0.01时并且conv输出的时间轴也不在原位需要手动推。dt 0.01; t -1:dt:4; u (t) double(t0); x u(t) - u(t-1); % 输入宽度 1 的门脉冲 h exp(-2*t) .* u(t); % 冲激响应单边指数衰减 y conv(x, h) * dt; % 求和 × 步长 ≈ 积分 ty (t(1)t(1)) (0:numel(y)-1)*dt; % 卷积结果的时间轴 plot(ty, y); grid on;逻辑说明conv(x,h)输出长度是numel(x)numel(h)-1起点是两个时间轴起点之和所以ty从-1 (-1) -2开始步长仍是dt。乘dt是把离散和Σ x[k]h[n-k]还原成积分∫ x(τ)h(t-τ)dτ的量纲修正不乘的话幅度和理论值差两个数量级。参数说明dt越小近似越准但conv的计算量按点数平方增长dt0.001时点数上千普通笔记本跑几十毫秒没问题再小就要考虑fft加速。另外时间轴范围要保证卷积结果「走完」——指数衰减信号尾部没衰减到零就截断会看到曲线在右端突然掉下来那是截断伪影不是物理现象。理论值可以用来对表。上面这组参数下0 t 1时解析解是(1 - e^{-2t})/2t 1时是(1 - e^{-2})e^{-2(t-1)}/2。取t 0.5代入得 0.3161从y里取对应位置的值对比误差应该在小数点后三位以内。误差偏大优先怀疑dt太大或者时间轴没对齐。4. 频域三件套傅里叶级数、频率特性与零极点怎么对上时域画波形是体力活频域才是真正容易出「图能画出来但结论是错的」的地方。报告后半段从周期矩形脉冲的频谱一路讲到连续系统的频率特性、拉普拉斯变换和零极点分布。这一章挑三个最常被问到的点级数系数怎么算才和fft对得上、频率特性曲线为什么有人画出来是反的、零极点图该用pzmap还是zplane。4.1 周期矩形脉冲手算系数和 fft 是两套坐标报告第六章的周期矩形脉冲周期T 2脉宽τ 1幅度 1。双边傅里叶级数系数是c_k (Aτ/T)·sinc(kτ/T)这里sinc用的是归一化定义sin(πx)/(πx)MATLAB 内置的sinc正好就是这个定义可以直接用。T 2; tau 1; A 1; k -20:20; c (A*tau/T) * sinc(k*tau/T); % 解析系数sinc sin(pi*x)/(pi*x) stem(k, abs(c), filled); grid on; title(周期矩形脉冲的双边幅度谱);用fft复现同一件事要注意三件事除以N得到的是 DFS 系数而不是连续谱要用fftshift把零频挪到中间采样总时长必须是信号周期的整数倍否则频谱泄漏包络会和解析值错开。fs 100; N 1000; % 总时长 10 s正好 5 个周期 t (0:N-1)/fs; x double(mod(t, T) tau); % 周期矩形脉冲占空比 0.5 X fftshift(fft(x)) / N; % 除以 N 才是级数系数 kf (-N/2 : N/2-1); stem(kf, abs(X), .); xlim([-20 20]); grid on;逻辑说明mod(t, T) tau是构造周期脉冲最简洁的写法等价于「每个周期内前 τ 秒为 1」。除以N之后X的纵轴和手算的|c_k|同量纲可以直接叠在一张图上比对。参数说明N/fs必须是T的整数倍这一条不满足时谱线会从离散谱线扩散成一片看起来像噪声。注意stem和plot不要混用在同一张频谱图上。频谱是离散谱线用stem连续谱包络用plot。混着画容易让人误以为采样点之间也有能量。4.2 频率特性手工算和 bode 应该是同一条线报告第八章讲连续时间 LTI 系统的频率特性。以H(s) 1/(s² 2s 1)为例令s jω得到H(jω) 1/((1-ω²) j2ω)。手写代码时最容易漏括号写成1/(1j*w).^2 ...MATLAB 会先算(1j*w).^2的逐元素倒数再加后面的项结果是一条完全不同的曲线而且不报错。w linspace(0, 20, 800); s 1j * w; H 1 ./ (s.^2 2*s 1); % 分母整体加括号逐元素除 subplot(2,1,1); plot(w, 20*log10(abs(H))); grid on; ylabel(幅度 (dB)); subplot(2,1,2); plot(w, angle(H)*180/pi); grid on; ylabel(相位 (度)); xlabel(\omega (rad/s));逻辑说明abs(H)是幅频取20*log10转 dBangle(H)是相频乘180/pi转角度。分母用./是逐元素除法因为w是向量。参数说明linspace(0, 20, 800)的 800 是画图点数只影响曲线平滑度不影响数值结果如果只看谐振峰附近把范围收到[0 5]会看得更清楚。用bode(tf([1],[1 2 1]))画一遍对照两条线应该完全重合。对不上时先查括号再查tf的系数顺序——tf(num, den)里den是按s的降幂排的写反了极点位置全变但曲线形状可能还挺像属于典型的「看起来对其实错」。4.3 零极点图连续域用 roots离散域用 zplane报告第十一章和第十四章分别是连续系统的零极点分析和 Z 域分析。两章的图画法完全不同混用会出错。% 连续域H(s) (s1)^2 / (s^23s2) b [1 2 1]; % 分子系数降幂 a [1 3 2]; % 分母系数降幂 zs roots(b); % 零点 ps roots(a); % 极点 plot(real(zs), imag(zs), o, real(ps), imag(ps), x); grid on; axis([-3 1 -2 2]); % 离散域H(z) (z - 0.5) / (z^2 - 1.2z 0.72) figure; zplane([1 -0.5], [1 -1.2 0.72]); % 直接吃系数向量自动画单位圆逻辑说明连续域的零极点就是多项式根roots足够离散域必须画单位圆作参考zplane会自动加上roots不会。参数说明zplane(b,a)的b、a也是降幂系数和tf一致如果手上是零极点增益形式用zpk(z,p,k)再转zplane更直接。判断稳定性有个口诀连续系统看极点实部是否全为负离散系统看极点模值是否全小于 1。写代码验证就是all(real(ps) 0)和all(abs(pz) 1)比盯着图数格子可靠。5. 常见问题排查五个「图不对」的具体归因与修法这一章记录的是我在复现过程中真实遇到过的五类问题。它们的共同点是MATLAB 不报错脚本跑完了图也出来了但图是错的。信号处理类代码的调试难度大部分集中在这里。5.1 图是空的、只有坐标轴uCT 未定义与路径问题现象运行报告第 7 页的2.3.1(3)命令行没有任何红色报错弹出的图里只有一对灰色坐标轴没有曲线。原因代码里调用了uCT(t)但整份报告没有附带uCT.m文件。MATLAB 遇到未定义函数会报Undefined function但如果脚本外面包了try-catch或者用的是旧版eval写法错误会被吞掉y保持空数组plot画空数组不报错。解决不要去找这个不存在的函数直接用匿名函数替代。在脚本开头加一行u (t) double(t0);把uCT(t)全部替换成u(t)。如果坚持要用函数文件确认它和主脚本在同一目录或者用addpath显式加进去然后用which uCT确认 MATLAB 找到的是哪一份。提示这类「函数找不到」的问题which 函数名和exist(函数名,file)两条命令能省掉大量猜测。前者告诉你找到的是哪个文件后者只告诉你有没有。5.2 反褶后波形整体偏移fliplr 的前提被破坏现象求奇偶分量时偶分量看起来正常奇分量在t0处不为零整条曲线像是被往左推了半个格。原因fliplr是索引反转不是时间反褶。当时间轴的采样中点不落在t0上时反转后的数组对应的实际时间是-t 2*t_mid相当于多平移了一个2*t_mid。t -3:0.01:2.99的中点t_mid -0.005偏移量 0.01正好是一个采样点。解决把时间轴改成严格对称的-3:0.01:3或者改用interp1(t, f, -t, linear, 0)。后者更通用代价是需要多一行代码和一次插值。改完用fo(find(t0))检查奇分量在原点是否为 0这是最快的验收标准。5.3 卷积幅度差两个数量级、峰值位置也对不上现象做完卷积曲线形状和理论一致但纵轴数值是理论值的 100 倍或者峰值位置比理论晚了 1 秒。原因两个独立错误。幅度问题是没有乘dtconv算的是离散和连续卷积是积分差一个步长因子dt0.01时正好差 100 倍。位置问题是时间轴起点算错ty的起点应该是两个输入时间轴起点之和很多人直接沿用第一个信号的t导致整体平移。解决y conv(x, h) * dt;和ty (t(1)t(1)) (0:numel(y)-1)*dt;两行必须成对出现。检查方法找一个理论值已知的点比如上文的t0.5处应为 0.3161用interp1(ty, y, 0.5)取值对比。5.4 相频曲线满屏竖线angle 的分支切割现象画angle(y)得到的相位曲线在若干位置突然从 π 跳到 -π看起来像噪声完全读不出趋势。原因angle的值域是 (-π, π]这是数学定义决定的不是 bug。跨分支时输出会回绕 2π。解决看整体趋势用unwrap(angle(y))需要严格对照理论相位时保留原始angle输出再单独画一条理论曲线叠加跳变位置就能对上。另一种做法是把相位转成弧度后用phase atan2(imag(y), real(y))效果和angle一样同样需要unwrap。5.5 axis 把波形裁掉了、subplot 里同一格画了两遍现象曲线顶部被削平或者四个子图里有两个长得一模一样。原因axis([xmin xmax ymin ymax])是硬裁剪不是自动缩放ymax给小了就直接切掉峰值。子图重复则是subplot(2,2,1)连续写了两次第二次覆盖第一次视觉上像是「某一格画了两条一样的线」。解决先plot再axis并且axis的ymax取max(y)*1.1留 10% 余量或者干脆用axis tight让 MATLAB 自己算。子图编号抄一遍检查subplot(2,2,1)到subplot(2,2,4)每个只出现一次中间用hold on叠加而不是重开一格。6. 把整份报告收成一套模板参数自检表与三个固定习惯复现到第四遍的时候我把这份报告里的重复劳动抽成了一张检查表。每次写完一段信号代码跑之前对着表过一遍能挡掉八成「图不对但没报错」的情况。表里的判据都是从上面几章的具体坑里提炼的不是通用规范。检查项命令或方法合格判据步长是否够密1/(20*max_freq)dt小于该值阶跃/门函数来源which u返回脚本内定义的匿名函数不是空时间轴对称性t(1)t(end)结果为 0反褶前必查卷积量纲max(y)/max(y_theory)比值接近 1不是 100 或 0.01相位是否回绕max(abs(diff(angle(y))))大于 π 就用unwrap显示窗口axis tight峰值不被裁切最后一章想说的其实是一个动作把「出图」和「验算」绑在一起不要等图画完再去想对不对。% 每次出图后立刻执行的验算片段 assert(abs(max(abs(y)) - 4) 0.05, 复信号模的最大值应为 4检查 1j 是否被覆盖); assert(abs(fo(t0)) 1e-9, 奇分量在原点应为 0检查 fliplr 的对称性前提); fprintf(点数 %d步长 %.4f等效采样率 %.1f Hz\n, numel(t), dt, 1/dt);这三行代码不值钱值钱的是把它变成肌肉记忆。assert失败时给出的那句话比调试半小时后恍然大悟要划算得多。参数说明容差 0.05 和 1e-9 是经验值前者对应波形幅度的舍入误差后者对应浮点意义上的零。还有一个习惯值得固定下来每个.m文件开头写清楚三件事——信号最高频率、时间轴范围、预期峰值。这三行注释在换题重跑的时候是后悔药。我自己有过一次把functl(t-1)改成functl(2-t)之后忘了改axis范围图看起来正常其实是把反褶后的波形裁掉了大半直到交稿前对着理论表才发现横轴范围根本对不上。从那以后我每次改完自变量表达式都强制把axis改成axis tight先看一眼全貌再手工收窗口。这个动作只需要两秒但它拦住的错误往往要两小时才能找回来。希望帮到你。本文还有配套的精品资源点击获取
简介《信号与系统》实验报告是广东工业大学信息工程学院电子信息类学生的课程实践总结面向正在修读该课程、需要完成MATLAB仿真实验或整理实验文档的本科生。压缩包内为1个docx文档约1.68MB按实验一至实验四组织涵盖连续时间信号在MATLAB中的表示、卷积运算、连续与离散LTI系统的时域分析、周期信号的傅里叶级数及频谱分析、傅里叶变换性质、拉普拉斯变换与Z变换等内容。每章列出实验目的与实验内容配有MATLAB代码和结果波形图便于对照理论验证计算过程。文档已排版成册共90余页可作为实验报告的写作框架与参考范例帮助读者理解基本信号生成、卷积求解、采样量化影响及系统频率特性分析等关键环节节省从零搭建实验步骤的时间。目前已有1257人学习下载适合作为课程实验的辅助材料与复习提纲。1. 一份 92 页的《信号与系统》MATLAB 实验报告从画波形到算卷积的完整链路第一次翻这份《信号与系统》MATLAB 实验报告是在帮某高校电子信息类专业的同学定位一段画不出曲线的代码。报告从连续时间信号的表示一路排到 Z 域分析四组实验、十四个章节每节都是「公式、代码、波形、结论」四件套。真正值钱的地方不是答案而是它把课堂上那些抽象符号压成了能敲进命令行的东西阶跃怎么定义、卷积怎么近似、傅里叶级数频谱怎么画、零极点图怎么看。它适合三类人正在赶实验报告但想弄明白每一步为什么这么写的人、需要一份信号处理命令速查模板的人、以及想把时域频域这套东西重新捡起来的人。下面按我自己的复现顺序拆。2. 信号表示的基础把公式翻译成时间轴与匿名函数MATLAB 里没有真正意义上的连续信号plot画出来的永远是折线。这句话听起来像废话但它决定了后面每一个参数的取法。这一章先把「时间轴怎么切」「阶跃怎么定义」「复信号怎么拆四个子图」这三件事固定下来后面的时移、卷积、频谱全都建立在这三块砖上。2.1 时间步长的取法dt 不是随便填个 0.01报告里几乎每段代码都以t -2:0.01:4开头。0.01 这个数字看着随意实际上和信号最高频率绑在一起。以报告里的exp(-t).*sin(10*pi*t)为例角频率 10π频率 5 Hz一个周期 0.2 秒步长 0.01 意味着一个周期采 20 个点折线看起来已经足够光滑。换成sin(100*pi*t)50 Hz一个周期只剩 2 个点画出来是锯齿甚至因为混叠伪装成一个低频信号。% 时间轴三要素起点、步长、终点 dt 0.001; % 步长先按最高频率的 1/20 取画不平滑再减半 fmax 50; % 信号里最高的正弦分量频率单位 Hz t -2:dt:4; % 时间范围按想看到几个周期/几个转折点定 N numel(t); % 采样点数后面算频谱和卷积长度要用 fs 1/dt; % 等效采样率画频谱时必须知道这个值逻辑说明dt决定波形平滑度t的首尾决定视野两者互不影响。常见误区是把axis的范围当成时间轴范围来调——axis只裁剪显示窗口不改变参与运算的数据量。参数说明fmax这个变量很多人不写直接拍一个 0.01遇到高频信号就翻车养成先估fmax再定dt 1/(20*fmax)的习惯改起来有依据。2.2 单位阶跃与门函数匿名函数比符号工具箱稳报告里2.3.1(3)那段用了一个uCT(t)而同一页的其他代码又用(t 0)硬写风格不统一。这是个信号原作者中途换过写法。复现时最省事的方案是用匿名函数把阶跃和门函数一次定义好整个脚本复用。% 基本信号定义放在脚本最前面后面所有表达式都用它们 u (t) double(t 0); % 单位阶跃 u(t)t0 处取 1 g (t, a, b) u(t-a) - u(t-b); % 门函数区间 [a, b) 内为 1 t -2:0.01:4; y (2 - exp(-t)) .* u(t); % (2 - e^{-t})u(t) plot(t, y); axis([-2 4 -0.5 2.2]); % 顺序不能反先 plot 再 axis grid on; xlabel(t); ylabel(f(t));逻辑说明double(t 0)把逻辑数组转成 0/1 数值数组才能参与乘法。逐元素乘.*不能写成*因为u(t)返回的是和t同尺寸的数组矩阵乘会直接报维度错。参数说明门函数用u(t-a) - u(t-b)而不是ta tb好处是能直接代入时移后的变量比如g(2*t1, 0, 1)做尺度变换时不用改函数体。提示报告第 6 页2.3.1(2)的代码里有连续两行给y赋值前一行y t.*((t 0) - (t 1))是残留代码被后一行覆盖了。抄作业时留着不影响结果但调试时如果看到波形对不上先检查有没有这类覆盖行。2.3 复信号的四个子图real/imag/abs/angle 各管一段复信号f(t) 2 e^{j(π/4)t} e^{j(π/2)t}是报告里比较少见的「需要拆成四张图」的题。四个分量分别是实部、虚部、模、相角对应real、imag、abs、angle四个函数。看起来机械但angle有个必须知道的特性返回值落在 (-π, π]相位跨过 ±π 时会跳变 2π画出来是一条满是竖线的锯齿。t -5:0.01:10; y 2 exp(1j*(pi/4)*t) exp(1j*(pi/2)*t); subplot(2,2,1); plot(t, real(y)); title(实部); grid on; axis([-5 10 0.5 4.5]); subplot(2,2,2); plot(t, imag(y)); title(虚部); grid on; axis([-5 10 -2 2]); subplot(2,2,3); plot(t, abs(y)); title(模); grid on; axis([-5 10 0.5 4.5]); subplot(2,2,4); plot(t, unwrap(angle(y))); title(相角展开后); grid on;逻辑说明虚数单位用1j而不是i因为i极易被for i 1:N覆盖成循环变量一旦覆盖后面所有复数运算全错而且报错信息指向的位置往往离真正出问题的地方很远。unwrap把相位跳变补回去适合看整体趋势如果要严格对照理论值直接看angle的原始输出更实在。参数说明subplot(2,2,k)的编号是按行优先先左后右、先上后下不要按「四个图随意排」的思路填。验证这一步有个便宜办法实部理论上是2 cos(πt/4) cos(πt/2)模的理论范围是[0, 4]三个单位向量同相时最大。用一行max(abs(y))和min(abs(y))对一下比盯着图看靠谱。3. 时移、反褶与卷积三个脚本级运算的边界条件理论上的f(2-t)、f(2t1)是纸面变换写成代码只改一个字符但改错的时候波形不会报错只会安安静静地画错。这一章把报告的第三章和第四章合并处理先用函数句柄把分段信号封装好再谈反褶为什么依赖时间轴的对称性最后落到卷积的数值近似——这是整份报告里最容易被「抄个大概」抄坏的一节。3.1 分段信号封装成函数文件而不是拼字符串报告里的functl.m是整份文档中唯一一个独立函数文件用四个区间的差分阶跃拼出一条折线。这种写法的好处是可以在任意自变量上调用functl(t-1)、functl(2-t)、functl(2*t1)都是同一份代码。function y functl(t) % 分段折线信号转折点位于 t -2, -1, 0, 1, 2 % 对应区间取值-1→0斜线、1常数、2常数、1→0斜线 y (t1) .* double(t-2 t-1) ... 1 .* double(t-1 t 0) ... 2 .* double(t 0 t 1) ... (-t2) .* double(t 1 t 2); end逻辑说明每个区间用一个逻辑数组做「开关」与对应表达式逐元素相乘四段相加得到完整波形。这和报告里用(t-2)-(t-1)做差分的写法等价但区间边界写得更直白改转折点时不容易错位。参数说明函数体内所有乘法必须是.*因为传入的t是数组如果只测试单个点用functl(0.5)也能跑。调用时的关键认知变换的是自变量不是函数值。functl(2-t)里的2-t在函数内部会被当作「新的时间轴」参与分段判断所以图形整体被压缩、反褶并平移而不是把画好的图左右翻一下。这一点想通了后面五个小问其实是同一个操作。表达式自变量替换几何效果f(t-1)t → t-1右移 1f(2-t)t → 2-t关于 t1 反褶f(2t1)t → 2t1压缩 2 倍后左移 0.5f(4-t/2)t → 4-t/2展宽 2 倍关于 t4 反褶3.2 奇偶分量分解fliplr 成立的前提条件报告里求偶分量和奇分量的代码只有三行y1 fliplr(y)、ye (yy1)/2、yo (y-y1)/2。这三行在特定条件下正确换一组时间轴就错而且错得很隐蔽——因为波形看起来还是像那么回事。fliplr做的是索引反转y1(i) y(N1-i)。它要等价于时间反褶f(-t)必须同时满足两个条件时间轴关于 0 对称且 0 落在某个采样点上。报告里t -3:0.01:3一共 601 个点索引 301 对应t 0fliplr后 301 号位置还是 301 号位置的元素恰好满足。% 错误示范时间轴不对称 t -3:0.01:2.99; % 601 个点但中点是 -0.005不是 0 f double(t0 t2).*(1-t); fr fliplr(f); % 整体偏移一个采样点奇分量在 t0 处不为零 % 稳妥写法按自变量重新采样 fr interp1(t, f, -t, linear, 0); % 在 -t 处插值取值越界补 0 fe (f fr) / 2; fo (f - fr) / 2;逻辑说明interp1的第四个参数是插值方法第五个是越界外推值。用-t作为查询点等价于「把时间轴翻过来重新采一遍」只要t是等间隔的就能对上不依赖对称性。参数说明如果-t的范围超出t超出部分按第五个参数补 0这正是反褶后信号为零的物理含义。验证奇偶分解是否正确看两个数就够max(abs(fo))应该等于max(abs(f))/2附近且奇分量在t0处的值必须为 0。这两条对不上说明反褶出问题了。3.3 连续卷积的数值近似conv 缺了 dt 就只是个求和报告第四章是连续时间信号的卷积运算但 MATLAB 的conv做的是离散求和。两者差一个步长因子这个因子在图上表现为幅度整体差 100 倍dt 0.01时并且conv输出的时间轴也不在原位需要手动推。dt 0.01; t -1:dt:4; u (t) double(t0); x u(t) - u(t-1); % 输入宽度 1 的门脉冲 h exp(-2*t) .* u(t); % 冲激响应单边指数衰减 y conv(x, h) * dt; % 求和 × 步长 ≈ 积分 ty (t(1)t(1)) (0:numel(y)-1)*dt; % 卷积结果的时间轴 plot(ty, y); grid on;逻辑说明conv(x,h)输出长度是numel(x)numel(h)-1起点是两个时间轴起点之和所以ty从-1 (-1) -2开始步长仍是dt。乘dt是把离散和Σ x[k]h[n-k]还原成积分∫ x(τ)h(t-τ)dτ的量纲修正不乘的话幅度和理论值差两个数量级。参数说明dt越小近似越准但conv的计算量按点数平方增长dt0.001时点数上千普通笔记本跑几十毫秒没问题再小就要考虑fft加速。另外时间轴范围要保证卷积结果「走完」——指数衰减信号尾部没衰减到零就截断会看到曲线在右端突然掉下来那是截断伪影不是物理现象。理论值可以用来对表。上面这组参数下0 t 1时解析解是(1 - e^{-2t})/2t 1时是(1 - e^{-2})e^{-2(t-1)}/2。取t 0.5代入得 0.3161从y里取对应位置的值对比误差应该在小数点后三位以内。误差偏大优先怀疑dt太大或者时间轴没对齐。4. 频域三件套傅里叶级数、频率特性与零极点怎么对上时域画波形是体力活频域才是真正容易出「图能画出来但结论是错的」的地方。报告后半段从周期矩形脉冲的频谱一路讲到连续系统的频率特性、拉普拉斯变换和零极点分布。这一章挑三个最常被问到的点级数系数怎么算才和fft对得上、频率特性曲线为什么有人画出来是反的、零极点图该用pzmap还是zplane。4.1 周期矩形脉冲手算系数和 fft 是两套坐标报告第六章的周期矩形脉冲周期T 2脉宽τ 1幅度 1。双边傅里叶级数系数是c_k (Aτ/T)·sinc(kτ/T)这里sinc用的是归一化定义sin(πx)/(πx)MATLAB 内置的sinc正好就是这个定义可以直接用。T 2; tau 1; A 1; k -20:20; c (A*tau/T) * sinc(k*tau/T); % 解析系数sinc sin(pi*x)/(pi*x) stem(k, abs(c), filled); grid on; title(周期矩形脉冲的双边幅度谱);用fft复现同一件事要注意三件事除以N得到的是 DFS 系数而不是连续谱要用fftshift把零频挪到中间采样总时长必须是信号周期的整数倍否则频谱泄漏包络会和解析值错开。fs 100; N 1000; % 总时长 10 s正好 5 个周期 t (0:N-1)/fs; x double(mod(t, T) tau); % 周期矩形脉冲占空比 0.5 X fftshift(fft(x)) / N; % 除以 N 才是级数系数 kf (-N/2 : N/2-1); stem(kf, abs(X), .); xlim([-20 20]); grid on;逻辑说明mod(t, T) tau是构造周期脉冲最简洁的写法等价于「每个周期内前 τ 秒为 1」。除以N之后X的纵轴和手算的|c_k|同量纲可以直接叠在一张图上比对。参数说明N/fs必须是T的整数倍这一条不满足时谱线会从离散谱线扩散成一片看起来像噪声。注意stem和plot不要混用在同一张频谱图上。频谱是离散谱线用stem连续谱包络用plot。混着画容易让人误以为采样点之间也有能量。4.2 频率特性手工算和 bode 应该是同一条线报告第八章讲连续时间 LTI 系统的频率特性。以H(s) 1/(s² 2s 1)为例令s jω得到H(jω) 1/((1-ω²) j2ω)。手写代码时最容易漏括号写成1/(1j*w).^2 ...MATLAB 会先算(1j*w).^2的逐元素倒数再加后面的项结果是一条完全不同的曲线而且不报错。w linspace(0, 20, 800); s 1j * w; H 1 ./ (s.^2 2*s 1); % 分母整体加括号逐元素除 subplot(2,1,1); plot(w, 20*log10(abs(H))); grid on; ylabel(幅度 (dB)); subplot(2,1,2); plot(w, angle(H)*180/pi); grid on; ylabel(相位 (度)); xlabel(\omega (rad/s));逻辑说明abs(H)是幅频取20*log10转 dBangle(H)是相频乘180/pi转角度。分母用./是逐元素除法因为w是向量。参数说明linspace(0, 20, 800)的 800 是画图点数只影响曲线平滑度不影响数值结果如果只看谐振峰附近把范围收到[0 5]会看得更清楚。用bode(tf([1],[1 2 1]))画一遍对照两条线应该完全重合。对不上时先查括号再查tf的系数顺序——tf(num, den)里den是按s的降幂排的写反了极点位置全变但曲线形状可能还挺像属于典型的「看起来对其实错」。4.3 零极点图连续域用 roots离散域用 zplane报告第十一章和第十四章分别是连续系统的零极点分析和 Z 域分析。两章的图画法完全不同混用会出错。% 连续域H(s) (s1)^2 / (s^23s2) b [1 2 1]; % 分子系数降幂 a [1 3 2]; % 分母系数降幂 zs roots(b); % 零点 ps roots(a); % 极点 plot(real(zs), imag(zs), o, real(ps), imag(ps), x); grid on; axis([-3 1 -2 2]); % 离散域H(z) (z - 0.5) / (z^2 - 1.2z 0.72) figure; zplane([1 -0.5], [1 -1.2 0.72]); % 直接吃系数向量自动画单位圆逻辑说明连续域的零极点就是多项式根roots足够离散域必须画单位圆作参考zplane会自动加上roots不会。参数说明zplane(b,a)的b、a也是降幂系数和tf一致如果手上是零极点增益形式用zpk(z,p,k)再转zplane更直接。判断稳定性有个口诀连续系统看极点实部是否全为负离散系统看极点模值是否全小于 1。写代码验证就是all(real(ps) 0)和all(abs(pz) 1)比盯着图数格子可靠。5. 常见问题排查五个「图不对」的具体归因与修法这一章记录的是我在复现过程中真实遇到过的五类问题。它们的共同点是MATLAB 不报错脚本跑完了图也出来了但图是错的。信号处理类代码的调试难度大部分集中在这里。5.1 图是空的、只有坐标轴uCT 未定义与路径问题现象运行报告第 7 页的2.3.1(3)命令行没有任何红色报错弹出的图里只有一对灰色坐标轴没有曲线。原因代码里调用了uCT(t)但整份报告没有附带uCT.m文件。MATLAB 遇到未定义函数会报Undefined function但如果脚本外面包了try-catch或者用的是旧版eval写法错误会被吞掉y保持空数组plot画空数组不报错。解决不要去找这个不存在的函数直接用匿名函数替代。在脚本开头加一行u (t) double(t0);把uCT(t)全部替换成u(t)。如果坚持要用函数文件确认它和主脚本在同一目录或者用addpath显式加进去然后用which uCT确认 MATLAB 找到的是哪一份。提示这类「函数找不到」的问题which 函数名和exist(函数名,file)两条命令能省掉大量猜测。前者告诉你找到的是哪个文件后者只告诉你有没有。5.2 反褶后波形整体偏移fliplr 的前提被破坏现象求奇偶分量时偶分量看起来正常奇分量在t0处不为零整条曲线像是被往左推了半个格。原因fliplr是索引反转不是时间反褶。当时间轴的采样中点不落在t0上时反转后的数组对应的实际时间是-t 2*t_mid相当于多平移了一个2*t_mid。t -3:0.01:2.99的中点t_mid -0.005偏移量 0.01正好是一个采样点。解决把时间轴改成严格对称的-3:0.01:3或者改用interp1(t, f, -t, linear, 0)。后者更通用代价是需要多一行代码和一次插值。改完用fo(find(t0))检查奇分量在原点是否为 0这是最快的验收标准。5.3 卷积幅度差两个数量级、峰值位置也对不上现象做完卷积曲线形状和理论一致但纵轴数值是理论值的 100 倍或者峰值位置比理论晚了 1 秒。原因两个独立错误。幅度问题是没有乘dtconv算的是离散和连续卷积是积分差一个步长因子dt0.01时正好差 100 倍。位置问题是时间轴起点算错ty的起点应该是两个输入时间轴起点之和很多人直接沿用第一个信号的t导致整体平移。解决y conv(x, h) * dt;和ty (t(1)t(1)) (0:numel(y)-1)*dt;两行必须成对出现。检查方法找一个理论值已知的点比如上文的t0.5处应为 0.3161用interp1(ty, y, 0.5)取值对比。5.4 相频曲线满屏竖线angle 的分支切割现象画angle(y)得到的相位曲线在若干位置突然从 π 跳到 -π看起来像噪声完全读不出趋势。原因angle的值域是 (-π, π]这是数学定义决定的不是 bug。跨分支时输出会回绕 2π。解决看整体趋势用unwrap(angle(y))需要严格对照理论相位时保留原始angle输出再单独画一条理论曲线叠加跳变位置就能对上。另一种做法是把相位转成弧度后用phase atan2(imag(y), real(y))效果和angle一样同样需要unwrap。5.5 axis 把波形裁掉了、subplot 里同一格画了两遍现象曲线顶部被削平或者四个子图里有两个长得一模一样。原因axis([xmin xmax ymin ymax])是硬裁剪不是自动缩放ymax给小了就直接切掉峰值。子图重复则是subplot(2,2,1)连续写了两次第二次覆盖第一次视觉上像是「某一格画了两条一样的线」。解决先plot再axis并且axis的ymax取max(y)*1.1留 10% 余量或者干脆用axis tight让 MATLAB 自己算。子图编号抄一遍检查subplot(2,2,1)到subplot(2,2,4)每个只出现一次中间用hold on叠加而不是重开一格。6. 把整份报告收成一套模板参数自检表与三个固定习惯复现到第四遍的时候我把这份报告里的重复劳动抽成了一张检查表。每次写完一段信号代码跑之前对着表过一遍能挡掉八成「图不对但没报错」的情况。表里的判据都是从上面几章的具体坑里提炼的不是通用规范。检查项命令或方法合格判据步长是否够密1/(20*max_freq)dt小于该值阶跃/门函数来源which u返回脚本内定义的匿名函数不是空时间轴对称性t(1)t(end)结果为 0反褶前必查卷积量纲max(y)/max(y_theory)比值接近 1不是 100 或 0.01相位是否回绕max(abs(diff(angle(y))))大于 π 就用unwrap显示窗口axis tight峰值不被裁切最后一章想说的其实是一个动作把「出图」和「验算」绑在一起不要等图画完再去想对不对。% 每次出图后立刻执行的验算片段 assert(abs(max(abs(y)) - 4) 0.05, 复信号模的最大值应为 4检查 1j 是否被覆盖); assert(abs(fo(t0)) 1e-9, 奇分量在原点应为 0检查 fliplr 的对称性前提); fprintf(点数 %d步长 %.4f等效采样率 %.1f Hz\n, numel(t), dt, 1/dt);这三行代码不值钱值钱的是把它变成肌肉记忆。assert失败时给出的那句话比调试半小时后恍然大悟要划算得多。参数说明容差 0.05 和 1e-9 是经验值前者对应波形幅度的舍入误差后者对应浮点意义上的零。还有一个习惯值得固定下来每个.m文件开头写清楚三件事——信号最高频率、时间轴范围、预期峰值。这三行注释在换题重跑的时候是后悔药。我自己有过一次把functl(t-1)改成functl(2-t)之后忘了改axis范围图看起来正常其实是把反褶后的波形裁掉了大半直到交稿前对着理论表才发现横轴范围根本对不上。从那以后我每次改完自变量表达式都强制把axis改成axis tight先看一眼全貌再手工收窗口。这个动作只需要两秒但它拦住的错误往往要两小时才能找回来。希望帮到你。本文还有配套的精品资源点击获取
