与系统响应、逆乙变换相关的 matlab函数 1. filter.m 本文件用来求离散系统的输出y(n)。 >若系统的h(n)已知,由yn)=x(n)*h(n),用conm 文件可求出y(n)。y=conv(x,h) >filter件是在A(z)、B(z)已知,但不知道hn)的 情况下求yn)的。 调用格式是:y= filter(b,a,x) x,ya和b都是向量。 H(z)≤ B(x)b(1)+6(2)x2+632x2+…+6(m+1)z A(x)1+a22+(3)2+…+a(n2+1)z
17 1．filter.m 本文件用来求离散系统的输出y(n) 。 ➢若系统的h(n)已知，由y(n)=x(n)*h(n)，用conv.m 文件可求出y(n) 。 y=conv (x, h) ➢filter文件是在A(z)、B(z)已知，但不知道h(n)的 情况下求y(n)的。 调用格式是: y=filter(b, a, x) x, y, a 和 b都是向量。 与系统响应、逆Z变换 相关的matlab 函数
与逆乙变换相关的 matlab函数 1. filter m X=[1,234:y=[346] z conV= conv(X,y)%x,y为输入和单位脉冲响应时输出 z conV=con(yX)%x,y为输入和单位脉冲响应时输出 z filter=fter(y1,x)%x为输入,y为FIR单位脉冲响应时输出 z filter= filter(X,1y)%y为输入,x为FIR单位脉冲响应时输 z conv=31023363424 z conV=31023363424 z filter =3 10 23 36 z filter =3 10 23 可见, conv(X,y)总是等于conv(yX)。而fter(x,1,y)却不 定等于 filter(y1,x),但是它们总是conv(Xy)截短的结 果,截短的长度等于 -length( filter的第三个参数) 18
18 1．filter.m x=[1,2,3,4]; y=[3,4,6] z_conv= conv(x,y) % x , y 为输入和单位脉冲响应时输出 z_conv_= conv(y, x) % x , y为输入和单位脉冲响应时输出 z_filter=filter(y,1,x) % x为输入, y为FIR单位脉冲响应时输出 z_filter_=filter(x,1,y) % y为输入, x为FIR单位脉冲响应时输出 与逆Z变换 相关的matlab 函数 可见，conv(x,y)总是等于conv(y,x)。而filter(x,1,y)却不 一定等于filter(y,1,x)，但是它们总是conv(x,y)截短的结 果，截短的长度等于length(filter的第三个参数) z_conv = 3 10 23 36 34 24 z_conv_=3 10 23 36 34 24 z_filter = 3 10 23 36 z_filter_= 3 10 23
与逆乙变换相关的 matlab函数 H(2)=B(2)b6)62)x21+03)x2++b%+1)z A(z)1+02(3)2+…+a(n2+1)z 2.impz.m 在A(z)、B(z)已知情况下,求系统的单 位抽样响应h(n)。调用格式是 h=impz(b,a,N)或 Th, t]=impz(b,a,N) N是所需的的长度。前者绘图时n从1开始, 而后者从0开始。 19
19 与逆Z变换 相关的matlab 函数 2.impz.m 在 A（z）、B（z）已知情况下, 求系统的单 位抽样响应 h(n)。调用格式是： h = impz(b, a, N) 或 [h,t]=impz(b,a,N) N是所需的的长度。前者绘图时n从1开始， 而后者从0开始
3. residue. m X(z)的有理分式分解成简单有理分式的和 因此可用来求逆Z变换。调用格式: rp, k]=residue(b, a 假如知道了向量,p和k,利用 residue. m还可反 过来求出多项式A(z)、B(z)。格式是 b, a]=residue(rp, k) X()=5(2) r(1) A(x)1-(l)z +…+n D7+K(1)++…
20 3. residuez.m 将X(z) 的有理分式分解成简单有理分式的和， 因此可用来求逆Z变换。调用格式： [r,p,k]= residuez(b,a) 假如知道了向量r, p和k，利用residuez.m还可反 过来求出多项式A(z)、B(z)。格式是 [b,a]= residuez(r,p,k)
4.频率响应函数: fretz.m 知A(z)、B(z),求系统的频率响应。基本的调用格 式是: [H, w]=freqz(b, a,N,'whole, Fs) N是频率轴的分点数,建议N为2的整次幂;w是返回 频率轴座标向量,绘图用;Fs是抽样频率,若Es= 1,频率轴给出归一化频率;’ Whole指定计算的 频率范围是从0~FS,缺省时是从0~FS/2 幅频响应:Hr=abs(H); 相频响应: h()=B() Hphase=angle(H) 解卷绕: Hphase=unwrap( Hphase
21 4.频率响应函数：freqz.m 已知A(z)、B(z), 求系统的频率响应。基本的调用格 式是： [H,w]=freqz(b,a,N,'whole',Fs) N是频率轴的分点数，建议N为2的整次幂；w是返回 频率轴座标向量，绘图用；Fs是抽样频率，若Fs＝ 1，频率轴给出归一化频率；’whole’指定计算的 频率范围是从0～FS,缺省时是从0～FS/2. ( ) ( ) ( ) B z H z A z = 幅频响应：Hr=abs(H); 相频响应： Hphase=angle(H); 解卷绕： Hphase=unwrap(Hphase);