数字信号处理第二章上机题作业

数字信号处理作业实验题报告

第一章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)')

运行结果：

实验结论：

该系统稳定。