p31_17.m

eda实例DSP上机程序

%%%%%%%%%%%%%%%%%% 数字信号处理第一章习题P31第17小题程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear
B=0.866;
A=[1,-0.8,0.64];
%%%%系统单位脉冲响应
xn=[1,zeros(1,48)];
hn=filter(B,A,xn);
n=0:length(hn)-1;
subplot(2,1,1),stem(n,hn,'o');
title('(a)系统的单位脉冲响应');
xlabel('n');
ylabel('h(n)');
%%%%系统单位阶跃响应
xn=ones(1,100);
sn=filter(B,A,xn);
n=0:length(sn)-1;
subplot(2,1,2),stem(n,sn,'o');axis([0,30,0,2]);
title('(b)系统的单位阶跃响应');
xlabel('n');
ylabel('s(n)');