program_12_8.m

来自「数字信号处理—基于计算机的方法教材上的所有例题源程序」· M 代码 · 共 14 行

% Program 12_8
% Illustration of Overflow Limit Cycles
%
a = input('alpha_1 and alpha_2 values = ');
yi1 = -0.625; yi2 = -0.125;
for n = 1:41;
    y(n) = - a(1)*yi1 - a(2)*yi2;
    y(n) = a2dR(y(n),3);
    yi2 = yi1; yi1 = y(n);
end
k = 0:40;
stem(k,y)
xlabel('Time index n');ylabel('Amplitude');
title(['\alpha_1 = ' num2str(a(1)),'\alpha_2 = ' num2str(a(2))]);