数字信号处理作业实验题报告
第一章16.(1)
实验目的：
求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。

实验要求：
运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。

源程序：
B1=1;A1=[1, -0.6, 0.08];
ys=2;
%设差分方程
xn=[1, zeros(1, 20)];
%xn=单位脉冲序列，长度N=31
xi=filtic(B1, A1, ys);
hn1=filter(B1, A1, xn, xi);
%求系统输出信号hn1
n=0:length(hn1)-1;
subplot(2, 1, 1);stem(n, hn1, '.')
title('单位脉冲响应');
xlabel('n');ylabel('h(n)')
xn=ones(1, 20);
sn1=filter(B1, A1, xn, xi);
%求系统输出信号sn1
n=0:length(sn1)-1;
Subplot(2, 1, 2);
stem(n, sn1, '.')
title('单位阶跃响应');
xlabel('n');
ylabel('s(n)')
运行结果：
实验分析：
单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。

第二章31题
实验目的：
用matlab判断系统是否稳定。

实验要求：
用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。

源程序：
A=[2, -2.98, 0.17, 2.3418, -1.5147];
B=[0, 0, 1, 5, -50];
subplot(2,1,1);
zplane(B,A);
%求H(z)的极点
p=roots(A);
%求H(z)的模
pm=abs(p);
if max(pm)<1 disp('系统因果稳定'),
else,disp('系统因果不稳定'),end
un=ones(1,800);
sn=filter(B, A, un);
n=0:length(sn)-1;
subplot(2, 1, 2);plot(n, sn)
xlabel('n');ylabel('s(n)')
运行结果：
实验结论：
该系统稳定。