În matematică , în special în algebra liniară , criteriul Routh - Hurwitz determină numărul de rădăcini cu parte reală pozitivă și negativă a unui polinom pornind de la coeficienții săi, îmbunătățind criteriul Descartes . Este util, de exemplu, pentru a determina stabilitatea unui sistem dinamic liniar și invariant în timp cu o singură intrare și o singură ieșire (SISO).
Derivare
Criteriul este legat de teorema Routh-Hurwitz : {\ displaystyle ab = w (+ \ infty) -w (- \ infty)} , unde este:
- {\ displaystyle a} este numărul rădăcinilor în afară de partea negativă reală a polinomului {\ displaystyle f (z)} ;
- {\ displaystyle b} este numărul rădăcinilor în afară de partea reală pozitivă a polinomului {\ displaystyle f (z)} ;
- {\ displaystyle w (x)} este numărul de variații ale secvenței Sturm obținute din {\ displaystyle P_ {0} (x)} Și {\ displaystyle P_ {1} (x)} (pentru iterațiile ulterioare ale algoritmului lui Euclid ) unde{\ displaystyle f (ix) = P_ {0} (x) + iP_ {1} (x)} pentru un număr real {\ displaystyle x} . Prin teorema fundamentală a algebrei , orice polinom de grad {\ displaystyle n} trebuie avut {\ displaystyle n} rădăcini complexe. Prin urmare, avem condiția ca. {\ displaystyle f} este un polinom stabil (Hurwitz) dacă și numai dacă {\ displaystyle ab = n} . Prin urmare, putem înlocui condiția de pe {\ displaystyle a} Și {\ displaystyle b} cu o condiție pe secvența Sturm, care la rândul nostru ne va da o condiție pe coeficienții lui {\ displaystyle f} .
Descriere
Să se dea polinomul finit {\ displaystyle p (x) = a_ {n} x ^ {n} + a_ {n-1} x ^ {n-1} + \ ldots + a_ {0}} în care se presupune {\ displaystyle a_ {n}> 0.} Matricea Routh este construită:
- {\ displaystyle \ mathbf {R} (p) = {\ begin {bmatrix} a_ {n} & a_ {n-2} & a_ {n-4} & a_ {n-6} & \ ldots \\ a_ { n-1} & a_ {n-3} & a_ {n-5} & \ ldots \\ b_ {n-1} & b_ {n-2} & \ ldots & \\ c_ {n-2} & c_ {n-3} && \ end {bmatrix}}}
unde elementele {\ displaystyle b_ {i}} iar cele ulterioare sunt legate de coeficienți. Fiecare element corespunde raportului dintre determinantul matricei compus din elementele celor două rânduri superioare, în prima coloană și în coloana care urmează celei a elementului, și primul coeficient (schimbat în semn) al rândului imediat deasupra elementul fiind calculat:
- {\ displaystyle b_ {n-1} = {\ frac {\ begin {vmatrix} a_ {n} & a_ {n-2} \\ a_ {n-1} și a_ {n-3} \ end {vmatrix} } {-a_ {n-1}}}}
- {\ displaystyle b_ {n-2} = {\ frac {\ begin {vmatrix} a_ {n} & a_ {n-4} \\ a_ {n-1} și a_ {n-5} \ end {vmatrix} } {-a_ {n-1}}}}
- {\ displaystyle c_ {n-2} = {\ frac {\ begin {vmatrix} a_ {n-1} & a_ {n-3} \\ b_ {n-1} & b_ {n-2} \ end { vmatrix}} {- b_ {n-1}}}}
- {\ displaystyle c_ {n-3} = {\ frac {\ begin {vmatrix} a_ {n-1} & a_ {n-5} \\ b_ {n-1} & b_ {n-3} \ end { vmatrix}} {- b_ {n-1}}}}
- adică mai general
- {\ displaystyle k_ {i, j} = {\ frac {\ begin {vmatrix} k_ {i-2,1} & k_ {i-2, j + 1} \\ k_ {i-1,1} & k_ {i -1, j + 1} \ end {vmatrix}} {- k_ {i-1,1}}}}
Construcția se termină imediat ce există o singură matrice pătrată cu determinant nul, adică doi coeficienți consecutivi cu o singură valoare fiecare. De fapt, acolo unde nu sunt prezente, elementele matricilor trebuie considerate nule.
Fiecare variație (permanență) a semnului coeficienților primei coloane corespunde unei rădăcini cu o parte reală pozitivă (negativă).
Prezența de zerouri pe prima coloană
În cazul în care un termen din prima coloană este nul, există patru metode diferite.
Prima metodă
Prin substituirea simbolului {\ displaystyle \ varepsilon} să reprezinte un număr foarte mic în valoare absolută, având tendința de a {\ displaystyle 0 ^ {+}} oa {\ displaystyle 0 ^ {-}} .
Dacă celelalte numere din prima coloană sunt toate pozitive atunci {\ displaystyle \ varepsilon \ to 0 ^ {+}} .
Dacă celelalte numere din prima coloană sunt negative, atunci {\ displaystyle \ varepsilon \ to 0 ^ {-}} .
În caz contrar, trebuie luate în considerare ambele cazuri {\ displaystyle \ varepsilon \ to 0 ^ {+}} Și {\ displaystyle \ varepsilon \ to 0 ^ {-}} .
De exemplu:
- {\ displaystyle \ mathbf {R} (x ^ {5} + 2x ^ {4} + 2x ^ {3} + 4x ^ {2} + 3x) = {\ begin {bmatrix} 1 & 2 & 3 \\ 2 & 4 & 0 \\ 0 & 3 & 0 \ \\ end {bmatrix}}}
devine prin înlocuirea acestuia cu {\ displaystyle \ varepsilon \ to 0 ^ {+}}
- {\ displaystyle \ mathbf {R} (x ^ {5} + 2x ^ {4} + 2x ^ {3} + 4x ^ {2} + 3x) = {\ begin {bmatrix} 1 & 2 & 3 \\ 2 & 4 & 0 \\\ varepsilon & 3 & 0 \\ {\ frac {4 \ varepsilon -6} {\ varepsilon}} & 0 & 0 \\ 3 & 0 & 0 \\\ end {bmatrix}}}
se vede clar că {\ displaystyle \ varepsilon> 0} dar {\ displaystyle {\ frac {4 \ varepsilon -6} {\ varepsilon}} <0}
Această metodă este strict justificată numai atunci când polinomul nu are zerouri pe axa imaginară; din acest motiv, în unele cazuri, poate da naștere la erori (așa cum se poate vedea analizând {\ displaystyle \ mathbf {R} (x ^ {6} + x ^ {5} + 3x ^ {4} + 3x ^ {3} + 3x ^ {2} + 2x + 1)} ).
A doua metodă
Polinomul dat poate fi înmulțit cu un binom cu zero negativ (adăugând astfel un zero negativ la polinom, permițând analiza semnelor celorlalte zerouri): fiind {\ displaystyle p (x)} polinomul original, trecem la studiul polinomului {\ displaystyle p (x) \ cdot (xc)}
De exemplu:
- {\ displaystyle \ mathbf {R} (x ^ {4} + 2x ^ {3} + 3x ^ {2} + 6x + 4) = {\ begin {bmatrix} 1 & 3 & 4 \\ 2 & 6 & 0 \\ 0 & 4 & 0 \\\ end {bmatrix}}}
deci putem adăuga, de exemplu, un zero în -1:
- {\ displaystyle p (x) \ cdot (x + 1) = x ^ {5} + 3x ^ {4} + 5x ^ {3} + 9x ^ {2} + 10x + 4}
- {\ displaystyle \ mathbf {R} (x ^ {5} + 3x ^ {4} + 5x ^ {3} + 9x ^ {2} + 10x + 4) = {\ begin {bmatrix} 1 & 5 & 10 \ \ 3 & 9 & 4 \\ 2 & {\ frac {26} {3}} & 0 \\ - 4 & 4 & 0 \\ {\ frac {32} {3}} & 0 & 0 \\ 4 & 0 și 0 \\\ end {bmatrix}}}
A treia metodă
Se aplică și în prezența mai multor zerouri consecutive pe aceeași linie. Acesta constă în înlocuirea liniei în cauză cu șirul de numere obținut prin adăugarea la element {\ displaystyle i} -alea din linie elementul loc {\ displaystyle i + j} în același rând înmulțit cu {\ displaystyle (-1) ^ {j}} , fiind {\ displaystyle j} numărul primelor elemente nule.
De exemplu:
- {\ displaystyle \ mathbf {R} (x ^ {5} + 5x ^ {3} + 10x + 4) = {\ begin {bmatrix} 1 & 5 & 10 \\ 0 & 0 & 4 \\ 0 & 0 & 0 \\\ end {bmatrix}}}
În locul primei prime {\ displaystyle 4} (care este primul element non-nul al liniei) și înmulțiți-l cu {\ displaystyle (-1) ^ {2}} fiind {\ displaystyle 2} zerouri consecutive înainte {\ displaystyle 4} .
- {\ displaystyle \ mathbf {R} (x ^ {5} + 5x ^ {3} + 10x + 4) = {\ begin {bmatrix} 1 & 5 & 10 \\ 4 & 0 & 4 \\ 5 & 9 & 0 \\ - {\ frac {36} {5}} & 4 & 0 \\ {\ frac {106} {9}} & 0 & 0 \\ 4 & 0 & 0 \\\ end {bmatrix}}}
A patra metodă
Se poate întâmpla ca toți termenii unui rând să fie nul doar dacă rândul este de ordin neobișnuit; de fapt, cele două linii precedente trebuie să fie proporționale și, prin urmare, trebuie să aibă același număr de elemente (rețineți că atunci când treceți de la o linie impară la una pară, numărul de elemente nu se schimbă). În această împrejurare se poate concluziona că polinomul considerat este produsul a două polinoame: primul va avea zerouri care au o parte reală caracterizată prin variațiile semnului elementelor primei coloane a tabelului construite până acum ( zerouri ale {\ displaystyle p_ {1} (x)} în afară de realul pozitiv există atâtea variații de semn care au apărut în prima coloană a tabelului construită până în acel moment); al doilea polinom {\ displaystyle p_ {2} (x)} (care se numește ecuație auxiliară) este de grad egal cu indicele rândului care precede rândul care este anulat, are doar puteri de grad egal și coeficienții săi sunt în ordinea de la cel al gradului maxim la cel al gradului, i coeficienții rândului care precede cel care a fost anulat.
De exemplu:
- {\ displaystyle \ mathbf {R} (x ^ {5} + x ^ {4} + 3x ^ {3} + 3x ^ {2} + x + 1) = {\ begin {bmatrix} 1 & 3 & 1 \ \ 1 & 3 & 1 \\ 0 & 0 & 0 \\\ end {bmatrix}}}
- {\ displaystyle x ^ {5} + x ^ {4} + 3x ^ {3} + 3x ^ {2} + x + 1 = p_ {1} (x) p_ {2} (x)}
{\ displaystyle p_ {1} (x)} este gradul I și are un zero negativ, {\ displaystyle p_ {2} (x)} de grad {\ displaystyle 4} numai cu puteri de ordine uniformă:
- {\ displaystyle p_ {2} (x) = x ^ {4} + 3x ^ {2} +1}
Pentru {\ displaystyle p_ {2} (x)} putem calcula zerourile, dar dacă gradul este prea mare ar putea exista dificultăți. Apoi, construcția mesei poate fi reluată în alt mod. Se derivă {\ displaystyle p_ {2} (x)}
- {\ displaystyle p_ {2} '(x) = 4x ^ {3} + 6x}
iar coeficienții rândului nul (în acest caz al treilea) sunt înlocuiți cu acești noi coeficienți:
- {\ displaystyle \ mathbf {R} (4x ^ {3} + 6x) = {\ begin {bmatrix} 1 & 3 & 1 \\ 1 & 3 & 1 \\ 4 & 6 & 0 \\ {\ frac {3 } {2}} & 1 & 0 \\ {\ frac {10} {3}} & 0 & 0 \\ 1 & 0 & 0 \\\ end {bmatrix}}}
Observați că zerourile polinomului {\ displaystyle p_ {2} (x)} , cu puteri de numai grad egal, au o simetrie dublă față de fiecare axă a planului complex. Acest lucru asigură că, dacă nu există zero, altele decât partea reală pozitivă, toate sunt pe axa imaginară. Mai exact, a doua parte a tabelului trebuie interpretată după cum urmează: sunt numărate doar variațiile de semn corespunzătoare rădăcinilor cu partea reală pozitivă, pe care le indicăm cu {\ displaystyle np} ; sunt {\ displaystyle nn = np} rădăcinile cu partea reală negativă (dată simetria dublă numită și simetria cadranului), dacă ecuația auxiliară este de grad {\ displaystyle n} , apoi cele rămase {\ displaystyle n-np-nn} rădăcinile se află pe axa imaginară. În exemplul nostru, ecuația auxiliară este gradul {\ displaystyle 4} , asa de {\ displaystyle n = 4} . Variațiile sunt {\ displaystyle np = 0} , deci există {\ displaystyle nn = 0} rădăcini în afară de negativul real și deci {\ displaystyle 4-0-0 = 4} rădăcini pe axa imaginară.
Bibliografie
- A.Hurwitz, Math. Ann., Vol. 46, 1895, p. 273-284.
Elemente conexe
linkuri externe