与系统响应、逆Z变换相关的matlab函数 1.filter.m 本文件用来求离散系的输出y(n)。 >若系统的h(n)已知，由y(n)=x(n)*h(n),用conv.m 文件可求出yn)。y-conv(&,h) >filter文件是在A(z)、B(z)已知，但不知道h(n)的 情沉下求yn)的。 调用格式是：y=-filter(b,a,x) x,y,a和b都是向量。 H(x)= B{之) h(1)+(2)x1+b(3)x2+…+(+1)z” A(z) 1十a(2)z+a(3)z2+…+a(n.+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 函数
与逆Z变换相关的matlab函数 1.filter.m ×=[1,2,3,4];y=[3,4,6] Z_conV=conV(X,y)oX,y为输入和单位脉冲响应时输出 Z_conV_=conv(y,X)%x,y为输入和单位脉冲响应时输出 z_fter=-filter(y,1,X)%x为输入，y为FIR单位脉冲响应时输出 z_filter_.=filter(X,1,y)%y为输入，x为FIR单位脉冲响应时输H z conv=31023363424 Z_c0nV_=31023363424 z filter =3 10 23 36 z_filter=3 10 23 可见，conv(Xy)总是等于conv(y,x)。而filter(x,1,y)却不 一定等于filter(y,1,x),但是它们总是conv(x,y)截短的结 果，截短的长度等于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
与逆Z变换相关的matlab函数 H(x)= B{之) (1)+(2)x1+(3)x2+…+(%十1)z A(z) 1十a(2)z+a(3)z+…+a(na+1)z". 2.impz.m 在A(z)、B(z)已知情况下，求系统的单 位抽样响应h(n)。调用格式是： h=impz(b,a,N)或 [h,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.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). X()= B(z) r(1) r(n) A1-02+m+-2+1)+2z1+m 20
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.频率响应函数：freqz.m 已知A(z)、B(z),求系统的频率响应。基本的调用格 式是： [H,w]=freqz (b,a,N,'whole',Fs) N是频率轴的分点数，建议N为2的整次幂；w是返回 频率轴座标向量，绘图用；s是抽样频率，若Es= 1,频率轴给出归一化频率；’who1e指定计算的 频率范围是从0~FS,缺省时是从0~FS/2. 幅频响应：Hr=abs(H); H(z)= B() 相频响应： Hphase=angle(H); A( 解卷绕： Hphase=unwrap(Hphase); 21
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);