p74_31.m

eda实例DSP上机程序

%%%%%%%%%%%%%%%%%% 数字信号处理第二章习题P74第31小题程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear
A=[2,-2.98,0.17,2.3418,-1.5147];
B=[0,0,1,5,-50];
%%%%%用极点分布判断系统是否稳定
subplot(2,1,1);
zplane(B,A);
p=roots(A)
pm=abs(p)
if max(pm<1)
    disp('系统因果稳定');
else
    disp('系统非因果稳定');
end
%%%%%%%%画出u(n)的系统输出波形进行判断
un=ones(1,700);
sn=filter(B,A,un);
n=0:length(sn)-1;
subplot(2,1,2);
plot(n,sn);
xlabel('n');
ylabel('s(n)');