《 数 字 信 号 处 理 (第 三 版 ） 》第 10章 实 验 程 序 注 意 ： 每 次 实 验 前 ， 请 新 建 以 下 8个 m 文 件 并 依 实 验 需 要 （ 已注 明 ） 移 动 至 matlab当 前 工 作 目 录 下 ， 便 于 主 程 序 调 用 。各程 序 均已 验 证， 请 直接复制即可。 文 件 名 ： tstem.m（ 实 验 一 、二需要） 程序： function tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名 称（ 字符串） n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn'); axis([0,n(end),min(xn),1.2*max(xn)]); 文 件 名 ： tplot.m（ 实 验 一 、四需要） 程序： function tplot(xn,T,yn) %时 域 序 列 连 续 曲 线 绘 图 函 数 %xn:信 号 数 据 序 列 ,yn:绘 图 信 号 的 纵 坐 标 名 称 （ 字 符 串 ） %T为 采 样 间 隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文 件 名 ： myplot.m（ 实 验 一 、四需要） 程序 ： %(1)myplot;计算时 域 离散系统损耗函 数 并绘 制曲 线 图 。 function myplot(B,A) %B为 系统函 数 分子多项式系数 向量 %A为 系统函 数 分母多项式系数 向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅 度 (dB)') axis([0,1,-80,5]);title('损 耗 函数曲线'); 文 件 名 ： mstem.m（ 实 验 一 、三需要） 程序 ： function mstem(Xk) %mstem(Xk)绘 制 频 域 采 样 序 列 向 量 Xk的 幅 频 特 性 图 M=length(Xk); k=0:M-1;wk=2*k/M;%产 生 M 点 DFT对 应 的 采 样 点 频 率 （ 关 于pi归 一 化 值 ） stem(wk,abs(Xk),'.');box on;%绘 制 M 点 DFT的 幅 频 特 性 图 xlabel('w/\pi');ylabel('幅度'); axis([0,2,0,1.2*max(abs(Xk))]); 文 件 名 ： mpplot.m（ 实 验 一 需 要 ） 程 序 ： %(2)mpplot;计 算 时 域 离 散 系 统 损 耗 函 数 和 相 频 特 性 函 数 ， 并 绘 制 曲线 图 。 function mpplot(B,A,Rs) %mpplot(B,A,Rs) %时 域 离 散 系 统 损 ...
