p9_5.m

来自「《数字信号处理实验指导书（MATLAB 版）》 著者： （美）米特拉著 」· M 代码 · 共 12 行

% Program P9_5
% Granular Limit Cycles in First-Order IIR Filter
alpha = input('Type in the value of alpha = ');
yi = 0;x = 0.411; 
for n = 1:21
	y(n) = a2dR(alpha*yi,5) + x;
    yi = y(n); x = 0;
end
k = 0:20;
stem(k,y)
ylabel('Amplitude'); xlabel('Time index n')