p31_16.m

eda实例DSP上机程序

%%%%%%%%%%%%%%%%%% 数字信号处理第一章习题P31第16小题程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear
B1=1;
A1=[1,-0.6,0.08];
B2=[2,0,-1];
A2=[1,-0.7,0.1];
%%%%系统1
xn=[1,zeros(1,30)];
hn1=filter(B1,A1,xn);
n=0:length(hn1)-1;
subplot(2,2,1),stem(n,hn1,'o');
title('(a)系统1的单位脉冲响应');
xlabel('n');
ylabel('h(n)');
xn=ones(1,30);
sn1=filter(B1,A1,xn);
n=0:length(sn1)-1;
subplot(2,2,2),stem(n,sn1,'o');
title('(b)系统1的单位阶跃响应');
xlabel('n');
ylabel('s(n)');
%%%%系统2
xn=[1,zeros(1,30)];
hn2=filter(B2,A2,xn);
n=0:length(hn2)-1;
subplot(2,2,3),stem(n,hn2,'o');
title('(c)系统2的单位脉冲响应');
xlabel('n');
ylabel('h(n)');
xn=ones(1,30);
sn2=filter(B2,A2,xn);
n=0:length(sn2)-1;
subplot(2,2,4),stem(n,sn2,'o');
title('(d)系统2的单位阶跃响应');
xlabel('n');
ylabel('s(n)');