重庆大学通信学院数字信号处理实验第四次实验报告.docVIP

图1 第一题的频率响应 图2 第二题的频率响应 图3 第三题的频率响应 图4 第四题的频率响应 第五题： 零、极点分布图 结果单位抽样响应： y = 1.0000 0.2000 -0.4000 -0.2000 1.0000 0 0 用filter函数： b=[1 0.2 -0.4 -0.2 1]; a=1; %注意此处分母的表示 x=zeros(1,10); x(1)=1;%单位抽样脉冲 y=filter(b,a,x) stem(y); title(单位抽样响应) 用impz函数： b=[1 0.2 -0.4 -0.2 1]; a=1; y=impz(b,a); 单位阶跃响应： y = 1.0000 1.2000 0.8000 0.6000 1.6000 1.6000 1.6000 1.6000 1.6000 1.6000 程序： b=[1 0.2 -0.4 -0.2 1]; a=1; x=ones(1,20); x(1)=1;%单位阶跃信号 y=filter(b,a,x) 第六题： 零、极点分布图 B=[0.3 -0.5 0]; A=[1 -0.7 0.8]; zplane(B,A) %零、极点分布图 B=[0.3 -0.5]; A=[1 -0.7 0.8]; subplot(2,2,1);zplane(B,A) title(系统函数极零点分布图) %最小相移网络： B=[1 -0.6]; A=[1 -0.7 0.8]; subplot(2,2,2) zplane(B,A); title(最小相移网络零点分布图) %全通网络： B=[-0.6 1]; A=[1 -0.6]; subplot(2,2,3); zplane(B,A); title(全通网络零点分布图) %修正网络： B=[-0.6 1]; A=[1 -0.7 0.8]; subplot(2,2,4); zplane(B,A); title(修正网络零点分布图) 对比区别： B=[0.3 -0.5]; A=[1 -0.7 0.8]; w=0:0.01:6*pi; h=freqz(B,A,w) subplot(2,2,1); plot(abs(h)); title(系统函数极频率响应) subplot(2,2,2); plot(angle(h)); title(系统函数极频率相位响应) B=[-0.6 1]; A=[1 -0.7 0.8]; w=0:0.01:6*pi; [h,w]=freqz(B,A,w) subplot(2,2,3); plot(abs(h)); title(修正网络频率响应) subplot(2,2,4); plot(angle(h)); title(修正网络频率相位响应) 第七题： 说明grpdelay的使用方法，计算1和5中的群延时 b=[0.2 0 -0.3]; a=[1 -0.4 0.8]; w=0:0.01:3*pi; grpdelay(b,a,w); 第五题的群延时 b=[1 0.2 -0.4 -0.2 1]; a=[1]; w=0:0.01:3*pi; grpdelay(b,a,w); 群延时函数的使用方法： help grpdelay GRPDELAY Group delay of a digital filter. [Gd,W] = GRPDELAY(B,A,N) returns length N vectors Gd and W containing the group delay and the frequencies (in radians) at which it is evaluated. Group delay is -d{angle(w)}/dw. The frequency response is evaluated at N points equally spaced around the upper half of the unit circle. For an FIR filter where N is a power of two, the computation i