Pattern recognition

2014-02-12T18:00:40Z

86.92.21.59: /* Supervised classification */

Hello! <br>My name is King and I'm a 24 years old girl from Austria.<br><br>my web blog :: [https://www.youtube.com/watch?v=uboG8ZUzlqI music making software download free]

Surface chemistry of paper

2013-06-30T13:26:32Z

86.92.236.252: /* Pigments, Binders, and Co-binders */

In [[signal processing]] and [[control theory]], the '''Bistritz criterion''' is a simple method to determine whether a [[discrete system|discrete]] [[LTI system theory|linear time invariant (LTI) system]] is [[BIBO stability|stable]] proposed in <ref name="[1]">Y. Bistritz (1984) [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=1457261&sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A31341%29 Zero location with respect to the unit circle of discrete-time linear system polynomials], Proc. of the IEEE, 72 (9): 1131–1142.</ref> see also.<ref name="[2]">Y. Bistritz (2002) [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=989164&sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A21319%29 Zero location of polynomials with respect to the unit circle unhampered by nonessential singularities], IEEE Trans. CAS I, 49(3): 305–314.</ref> Stability of a discrete LTI system requires that its [[characteristic polynomial]]s

:<math> D_n(z)=d_0+d_1 z+d_2 z^2+ \cdots + d_{n-1}z^{n-1} + d_n z^n </math>

(obtained from its difference equation, its dynamic matrix, or appearing as the denominator of its transfer function) is a [[stable polynomial]], where <math> D_n(z)</math> is said to be stable if all its zeros are inside the unit circle, viz.

:<math> | z_k | < 1 , k=1,\ldots,n </math>,

where <math> D_n(z)=d_n \prod_{k=1}^n (z-z_k) </math>. The test determines whether <math> D_n(z) </math> is stable algebraically (i.e. without numerical determination of the zeros) . The method also solves the full zero location (ZL) problem. Namely, it can count the number of inside the unit-circle (IUC) zeros <math>(~ | z_k| < 1 ~) </math>, on the unit-circle zeros (UC) zeros <math>(~ | z_k| = 1 ~)</math> and outside the unit-circle (OUC) zeros <math> (~ | z_k| > 1 ~) </math> for any real or complex polynomial.<ref name="[1]"/><ref name="[2]"/>
The Bistritz test is the discrete equivalent of [[Routh–Hurwitz stability criterion|Routh]] criterion used to test stability of continuous LTI systems. This title was introduced for it in <ref>E. I. Jury and M. Mansour (1985), [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=1457477&sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A31348%29%26pageNumber%3D2 On the terminology relationship between continuous and discrete systems criteria], Proc. of the IEEE, 73(4):884.</ref> soon after its presentation. It has been also recognized to be more efficient than previously available stability tests for discrete systems like the Schur-Cohn and the [[Jury stability criterion|Jury test]].<ref>K. Premaratne, and E. I. Jury (1993) [http://www.sciencedirect.com/science/article/pii/001600329390028S On the Bistritz tabular form and its relationship with the Schur–Cohn minors and inner determinants], J. of the Franklin Institute, 30(1):165-182.</ref>
In the following, the focus is only on how to test stability of a real polynomial. However, as long as the basic recursion needed to test stability remains valid, ZL rules are also brought. For more general needs the references <ref name="[1]"/><ref name="[2]"/> should be consulted.

==Algorithm==

Consider <math> D_n(z)</math> as above and assume <math> D_n(1) \neq 0 </math>. (If <math> D_n(1)=0 </math> the polynomial is not stable.) Define its reciprocal polynomial

:<math>D^\sharp_n(z)=z^n D_n(1/z)=d_n+d_{n-1}z+d_{n-2} z^2+\cdots+d_{n-1}z^{n-1} + d_0 z^n </math>.

The algorithm assigns to <math>D_n(z)</math> a sequence of [[symmetric polynomial]]s

:<math>T_m(z)=T^\sharp_m(z), m=n,n-1, \ldots , 0 </math>

created by a three-term polynomial recursion. Write out the polynomials by their coefficients,

:<math>T_m(z)=\sum_{k=1}^m t_{m,k} z^k </math>,

symmetry means that

:<math>T_m(z)=t_{m,0}+t_{m,1} z + \cdots + t_{m,1} z^{m-1}+t_{m,0} z^m </math>,

so that it is enough to calculate for each polynomial only about half of the coefficients. The recursion begins with two initial polynomials driven from the sum and difference of the tested polynomial and its reciprocal, then each subsequent polynomial of reduced degree is produced from the last two known polynomials.

Initiation:
: <math> T_n(z)=D_n(z)+D^\sharp_{n}(z) \quad, \quad T_{n-1}(z)=\frac{ D_n(z)-D^\sharp_{n}(z) }{z-1}</math>

Recursion: For <math> m=n-1,\ldots,1 </math> do:

: <math> \delta_{m+1}=\frac{T_{m+1}(0)}{T_m(0)} </math>

: <math> T_{m-1}(z)=\frac{\delta_{m+1}(1+z) T_m(z) - T_{m+1}(z) }{z}</math>

==Stability condition==

The successful completion of the sequence with the above recursion requires
<math> T_m(0) \neq 0 , \quad m=n-1, \ldots ,1 </math>. The expansion of these conditions into
<math> T_m(0) \neq 0 , \quad m=n, \ldots , 0 </math>
are called normal conditions.

Normal conditions are necessary for stability. This means that, the tested polynomial can be declared as not stable as soon as a <math> T_m(0) = t_{m,0}=t_{m,m}=0 </math> is observed. It also follows that the above recursion is broad enough for testing stability because the polynomial can be declared as not stable before a division by zero is encountered.

'''Theorem'''. If the sequence is not normal then <math> D_n(z) </math> is not stable.
If normal conditions hold then the complete sequence of symmetric polynomials is well defined. Let

: <math> \nu = Var \{ T_n (1), T_{n-1}(1) , \ldots , T_1(1) , t_{0,0} \} </math>

denote the count of the number of sign variations in the indicated sequence. Then
<math> D_n(z) </math> is stable if and only if <math> \nu =0 </math>.
More generally, if normal condition hold than <math> D_n(z) </math> has no UC zeros,
<math> \nu </math> OUC zeros and <math> n- \nu </math> IUC zeros.

Violation of various necessary conditions for stability may be used advantageously as early indications that the polynomial is not stable (has at least one UC or OUC zero). The polynomial can be declared not stable as soon as a <math> T_m (0)=0 </math>, or a <math> \delta_{m} < 0 </math>, or a change of sign in the sequence of <math> T_m (1) </math>'s is observed.

==Example==
Consider the polynomial <math> D_3(z)=2+Kz -22 z^2 +24 z^3</math>, where <math> K </math> is a real parameter.

Q1:For what values of <math> K </math> the polynomial is stable?

Construct the sequence:

: <math> T_3(z)=26+(K-22)z +(K-22)z^2 +26 z^3</math>

: <math> T_2(z)=22-Kz +22z^2 </math>

: <math> T_1(z)=\frac{24(22-K)}{11} (1+z) </math>

: <math> T_0(z)=44+k </math>

Use their values at z=1 to form

: <math> \operatorname{Var}(8+2K, 44-K , 48(22-K)/11, 44+k) \, </math>

All the entries in the sequence are positive for -4 < K < 22 (and for no ''K'' are they all negative). Therefore D(z) is stable for −4 < ''K'' < 22.

Q2: Find ZL for K = 33 Var { 71, 11, -48, 11 }=2 => 2 OUC, 1 IUC zeros.

Q3: Find ZL for K = -11 Var{ -14, 55, 144, 33 }=1 => 1 OUC, 2 IUC zeros.

== Comments ==

(1) The test bears a remarkable similarity to the [[Routh–Hurwitz stability criterion|Routh]] test. This is best observed when the Routh test is arranged appropriately into a corresponding three-term polynomial recursion.

(2) The Bistritz test uses three-term polynomial recursion that propagates polynomials with symmetry as opposed to previously available classical tests for discrete systems that propagate polynomials with no particular structure using a two-term recursion. It stimulated the discovery of more algorithms in the area of digital signal processing (e.g. solving the [[linear prediction]] problem) and discrete systems (e.g. testing stability of higher dimensional systems) collectively called "immittance" or "split" algorithms that adopted this technique to more efficient counterparts to also other classical so called "scattering" algorithms, see <ref>P. Delsarte and E. Genin (1986) [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=1164830&sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A26198%29%26pageNumber%3D2 The split Levinson algorithm] IEEE Trans. ASSP 34(3):470-478.</ref><ref>Y. Bistritz, H. Lev-Ari and T. Kailath (1989) [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=30994&sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A1333%29%26pageNumber%3D2 Immittance-domain Levinson algorithms] IEEE Trans. IT, 35(3):675-682.</ref> and the chapter on linear prediction in the textbook.<ref>S. J. Orfanidis, [http://www.ece.rutgers.edu/~orfanidi/osp2e/osp2e.pdf Optimum signal processing: An introduction], 2nd Edition, Macmillan, 1988.</ref> The Bistritz test forms the "immittance" counterpart of the "scattering" type classical tests of Schur–Cohn and [[Jury stability criterion|Jury]].

==References==
<references />

{{DEFAULTSORT:Linear System}}
[[Category:Systems theory]]
[[Category:Stability theory]]
[[Category:Digital signal processing]]

formulasearchengine - User contributions [en]

Pattern recognition

Surface chemistry of paper