数字信号处理 实验报告1 时域离散系统的频率响应

1、数字信号处理实验报告专 业：.班 级：学 号：姓 名：指导教师：课程名称数字信号处理班级实验日期姓名学号实验成绩实验名称实验一时域离散系统的频率响应实验 目的 要求1、理解差分方程的求解方法；2、掌握系统函数零、极点的求解方法3、掌握用matlab求解系统的频率响应、单位冲激响应的方法。实验环境Matlab软件实 验 内 容1、用MATLAB计算差分方程y(n) + 0,7y伽-1)- 0 45y(n - 2) - 0.6-3)=0.8x(w) - 0 44x( -1) + 0.36x(n - 2) + 0.02x(-3)当输入序列x(n)=u(n)-u(n-5时，该差分方程输出结果的前41个

2、样点0 402、用MATLAB求以下系统函数的零、极点，并画出零、极点在Z平面上的 分布图皿、1-0.3z-2 - 0.3z-3 - 0.2z”Hz)=.- i04-1 + 0. Iz-1 + 0,2z-2 + 0. 2z-3 + 0.5z“3、用MATLAB计算差分方程y + 0.7y(忽-1) - 0.45yg-2)- 0.6y(冏-3)=0.8x()- 0.44x(n - 1) + 0,36式% -2) + 0.02K伽 一 3)所对应的离散LTI系统的频率响应，并画出其实部、虚部、幅度、相位图。 再求出并绘出该系统的单位冲激响应。1、求解差分方程clear all;N=41;a=0.8

3、 -0.44 0.36 0.02;实验程序清单b=l 0.7 -0.45 -0.6;实验程序清单x=l 1111 zeros(l,N-5);k=0:l:N-l;y=filter(a,b,x);stem(k,y) xlabel(H);ylabelC 幅度) grid on2、系统函数的零、极点 clear all;num=l -0.1 -0.3 -0.3 -0.2;den=l 0.1 0.2 0.2 0.5;z,p,k=tf2zp(num,den);dispC 零点);disp ;dispC极点);disp(p);zplane(num,den)grid ontitle,零、极点分布图)xlabe

4、l，实部)ylabe(虚部)grid on3、离散LTI系统的频率响应、单位冲激响应 clear all;k=256;num=0. 8 -0. 44 0. 36 0. 02;den=l 0. 7 -0. 45 -0. 6;w=0:pi/k:pi;h=freqz (num, den, w);subplot (2, 2, 1);plot(w/pi, real (h);gridtitle (实部)xlabel C omega/pi，) ; ylabel C 幅度) subplot (2, 2, 2);plot (w/pi, imag(h);gridtitle (虚部)xlabelC omega/pi

5、J);ylabel( Amplitude5) subplot (2, 2, 3);plot (w/pi, abs (h);gridtitle(幅度谱)xlabel C omega/piJ) ; ylabel (幅值) subplot (2, 2, 4);plot(w/pi, angle (h);gridtitleC相位谱)xlabel C omega/pi，) ; ylabel C 弧度)figure (2);dimpulse (num, den) ; title (单位冲激响应); grid on;1、求解差分方程程序运行后的输出图像如下:2、系统函数的零、极点实 验 结 果实 验 结 果实 验 结 果(1)系统函数的零、极点为:实 验 结 果(2)系统函数的零、极点在Z平面上的分布图如下：3、离散LTI系统的频率响应、单位冲激响应(1)该离散LTI系统的频率响应的实部、虚部、幅度、相位图如下:(2)该离散LTI系统的单位冲激响应如下:1、简述离散LT