姓名: 学号: 实验日期:
实验题目:因果离散线性系统的时域分析
实验目的:实现由差分方程构成的数字滤波器:
实验内容:
1.用两种不同的程序计算数字滤波器的单位抽样响应 ,
给定差分方程:y(n)-y(n-1)+0.9y(n-2)=x(n)
a.编制文件himpz.m,实现数字滤波器的单位抽样响应
b. 编制函数文件hfilter.m,实现数字滤波器的单位抽样响应
2.给定 计算此数字滤波器的单位阶跃响应 ,并判断系统的稳定性。用stem(n,y)画出相应的图形。参考程序gfilter.m
实验地点:4305机房
实验结果:
a=[1,-1,0.9];b=1;
impz(b,a)
title('单位抽样响应')
a=[1,-1,0.9];b=1;
% 求h(n)
x=impuls (0,-10,100);n=[-10:100];
h=filter(b,a,x);
subplot(2,1,1);stem(n,h)
axis([-10,100,-1.1,1.1])
title(‘单位冲激响应’);xlabel(‘n’);ylabel(‘h(n)’)
%求单位阶跃响应
a=[1,-1,0.9];b=1;
x=[zeros(1,10),ones(1,101)];
n=-10:100;
s=filter(b,a,x);
subplot(2,1,2);stem(n,s)
axis([-10,100,-.5,2.5])
title(‘单位阶跃响应’);xlabel(‘n’);ylabel(‘s(n)’)
%判断系统稳定
sum(abs(h))
z=roots(a);
magz=abs(z)
结果分析:
总结: