离散系统的频率响应分析和零极点分布.docVIP

实验2离散系统的频率响应分析和零、极点分布 一、 实验目的 通过MATLAB仿真简单的离散时间系统，研究其时域特性，加深对离散系统的冲激 响应，频率响应分析和零、极点分布的概念的理解。 二、 基本原理 离散系统的吋域方程为 N M X - k) = Pkx(n - k) 其变换域分析方法如下： OO y[n] = x[n] * h[n] = ^x[m]h[n 一 m] 0 Y(严)=X(eja,)H(ejaf) 频域 H (严)==久+卩旷加+…+亿“旷 系统的频率响应为 D(严)心+ d严+... +厶八血 oo y[n] = x[n]*h[n] = ^x[m]h[n-m] o Y(z) = X(z)H(z) H(z)= P ⑵=A)+ PzT+???+PwZ「 系统的转移函数为 D⑵ 心+ “厂 +??? + % 2 Ph n(i-6z_) 恥)=—=K . Ydkz~l H(l —久. 分解因式 z 曰 ，其中h和称为零、极 点。 在MATLAB中，可以用函数[z, p, K]=tf2zp (num, den)求得有理分式形式的 系统转移函数的零、极点，用函数zplane (z, p)绘出零、极点分布图；也可 以用函数zplane (num, den)直接绘出有理分式形式的系统转移函数的零、极 点分布图。 另外，在MATLAB中，可以用函数[r, p, k]=residuez (num, den)完成部 分分式展开计算；可以用函数sos二zp2sos (z, p, K)完成将高阶系统分解为2 阶系统的串联。 三、实验内容及要求 一个LT1离散时间系统的输入输出差分方程为 y (n) -1. 6y (nT) +1. 28y (n-2) =0. 5x (n) +0. lx (n~l) 编程求出此系统的单位冲激响应序列，并画出其波形。 若输入序列 x(n)=5 (n) +2 5 (n-1)+3 6 (n-2) +4 5 (n~3) +5 5 (n-4),编程 求此系统输岀序列y(n),并画出其波形。 编程得到系统频响的幅度响应和相位响应，并画图。 编程得到系统的零极点分布图，分析系统的因果性和稳定性。 解答： (1) elf; N=40; num= [0. 5, 0. 1]; den=[l, -1. 6, 1. 28]; y=impz (num, den, N) stcm(y); xlabel (吋间信号n，); ylabel (信号幅度); titleC冲击响应); o o o o o o o6 4 2-2-4W sio80■5 10 15 o o o o o o o 6 4 2-2-4W si o 80 ■ 5 10 15 20 25 30 35 40 * “ 5翻爵n )Figure Ko. 1 匚1面冈 X 丄U Xu J Vi ew Insert Tools Window Help A Z / 妙 Q C 冲击响应 (2) elf; N 二 40; num=[0. 5, 0. 1]; den=[l, -1. 6, 1. 28];x=[l, 2, 3, 4, 5, zeros (1, N-l) J; y=fiIter (num, den, x) stem(y); xlabel (时间信号n); ylabel (信号幅度)； title(，输出波形); J^Figure Ho. 1 File Edit View Insert Tools Window Help AZ/ 更 Q C输出波形00T—10s暑 AZ/ 更 Q C 输出波形 00 T— o 80 00 6 400 悝10s暑 o 20 o 20 ■ o 40 - o 60 - o 00 8 ■ Illi 1 1 1 ■ ■ Q Q% Q o o ) o 1 1 1 ? 1丄 d — Illi 1 1 1 i ) ■ 10 15 20 25 30 35 40 45 时间信号n 函数 y-cov (x, h) 函数 y-cov (x, h)和 y=fi Iter (num, den, x)的区别 elf; N 二 40; num=[0. 5, 0. 1]; den=[l, -1. 6, 1. 28]; x=[l, 2, 3, 4, 5, zeros (1, N-l)]; h=impz (num, den, N); y=conv(x, h) subplot (2, 1, 1); stem(y); xlabel(，时间信号n); ylabel (信号幅度)； titleC输出波形)； y=f i 1 ter (num, den, x); subplot (2, 1, 2); stem(y);xlabelC时间信号n); yla