数字信号处理20104200王大元实验4.doc

实验4 离散系统的变换域分析 实验目的 熟悉对离散系统的频率响应分析方法； 加深对零、极点分布的概念理解。 实验原理 离散系统的是时域方程为 其变换域分析方法如下： 频域： 系统的频率响应为： Z域： 系统的转移函数为： 分解因式： 其中 和 称为零、极点。 三、实验要求： 编程实现系统参数输入，绘出幅度频率响应和相位曲线和零、极点分布图。 求系统 num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336]; [z,p,k]=tf2zp(num,den); disp(零点);disp(z); disp(极点);disp(p); disp(增益系数);disp(k); figure(1) zplane(num,den) figure(2) freqz(num,den,128) 4.12 求滤波器输出的前十个采样值； num=[0.3 -0.25 0.1]; den=[1 0 0]; [z,p,k]=tf2zp(num,den); m=abs(p); disp(零点);disp(z); disp(极点);disp(p); [z,p,k]=tf2zp(num,den); m=abs(p); disp(零点);disp(z); disp(极点);disp(p); subplot(2,2,1); zplane(num,den); title(零 极点图); k=256; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,2); subplot(2,2,3); plot(w/pi,abs(h));grid title(幅度谱) xlabel(\omega/\pi);ylabel(幅值) subplot(2,2,4); plot(w/pi,angle(h));grid title(相位谱) xlabel(\omega/\pi);ylabel(弧度) clear; B=[0.3 -0.25 0.1]; A=1; n=0:9; X=[2*(n==0)-(n==2)]; Y=filter(A,B,X); stem(n,Y,.); title(滤波器输出采样值) 6.22 对于传输函数 求脉冲响应和阶跃响应 num=[5]; den=[1 -1 0.09]; subplot(2,1,1) impz(num,den,40); [h,t]=impz(num,den,40); title(脉冲响应) step=ones(1,100); subplot(2,1,2) y=conv(h,step); n=0:138; stem(n,y,filled) title(阶跃响应) 7.26 滤波器的传输函数为： H 不计算 ，画出幅度响应曲线。 k=256; num=[0,0,1,1]; den=[1,-1,0,0]; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1); plot(w/pi,real(h));grid title(实部); xlabel(\omega/\pi);ylabel(幅度) subplot(2,2,2); plot(w/pi,imag(h));grid title(虚部); xlabel(\omega/\pi);ylabel(Amplitude) subplot(2,2,3); plot(w/pi,abs(h));grid title(幅度谱); xlabel(\omega/\pi);ylabel(幅值) subplot(2,2,4); plot(w/pi,angle(h));grid title(相位谱); xlabel(\omega/\pi);ylabel(弧度) 7.28 由七项滑动平均滤波器的零极点推断滤波器形状 解：由七项滑动滤波器特点，写出其对应传输函数 num=[1 1 1 1 1 1 1]; den=[7 0 0 0 0 0 0]; [z,p,k]=tf2zp(num,den); disp(零点);disp(z); disp(极点);disp(z); disp(增益系数);disp(k); figure(1) zplane(num,den) figure(2) freqz(num,den,128)