www.gusucode.com > signal 工具箱matlab源码程序 > signal/@dfilt/@df2t/thislimitcycle.m

    function [y,zi,overflows] = thislimitcycle(Hd,x)
%THISLIMITCYCLE   

%   Author(s): V. Pellissier
%   Copyright 2005 The MathWorks, Inc.

q = Hd.filterquantizer;

stateq = quantizer([q.StateWordLength q.StateFracLength]);
% The values are distributed over the range of the state format
zi = randquant(stateq, size(Hd.HiddenStates)); 

num = Hd.privNum;
den = Hd.privDen;

% Reach steady state
N = impzlength(Hd);
[y,zii] = df2tfilter(q,num,den,x(1:N,:),zi);

% Look for limitcycles in steady state only
[y,zf,overflows] = df2tfilter(q,num,den,x(N+1:end),zii,'limitcycle');

Hd.HiddenStates = zf;

% [EOF]