matlab判断IIR传递函数稳定性

num = input('输入分子： ');
den = input('输入分母： ');
N = max(length(num),length(den));
x = 1; y0 = 0; S = 0;zi = zeros(1,N-1);
for n = 1:1000
    [y] = filter(num,den,x);
    if abs(y) < 0.000001, break, end
end

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