p9_6.m

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

% Program P9_6
% Overflow Limit Cycles in Second-Order IIR Filter
alpha = input('alpha1 and alpha2 values = ');
yi1 = 0.75;yi2 = -0.75; 
for n = 1:41
    y(n) = - alpha(1)*yi1 - alpha(2)*yi2;
	y(n) = a2dR(y(n),3);
    yi2 = yi1; yi1 = y(n);
end
k = 0:40;
stem(k,y)
ylabel('Amplitude'); xlabel('Time index n')