program_6_7.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

%Program_6_7
% Stability Test
%
num = input('Type in the numerator vector = ');
den = input('Type in the denominator vector = ');
N = max(length(num),length(den));
x = 1; y0 = 0; S = 0;zi = zeros(1,N-1);
for n = 1:1000
    [y,zf] = filter(num,den,x,zi);
    if abs(y) < 0.000001, break, end
    x = 0;
    S = S + abs(y);
    y0 = y;zi = zf;
end

if n < 1000
    disp('Stable Transfer Function');
else
    disp('Unstable Transfer Function');
end