program_9_8.m

来自「这是国外一本经典教材附带的程序」· M 代码 · 共 16 行

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