معیار پایداری راث-هرویتز

در نظریه سامانه کنترل، معیار پایداری راث-هرویتز یک آزمون ریاضی است که شرط لازم و کافی برای پایداری سیستم کنترل تغییرناپذیر با زمان خطی (LTI) است. محک راث یک الگوریتم بازگشتی کارآمد است که ادوارد جان راث ریاضیدان انگلیسی در سال ۱۸۷۶ برای تعیین اینکه آیا همه ریشههای چندجمله‌ای مشخصه یک سامانه خطی دارای قسمت‌های حقیقی منفی هستند، پیشنهاد داد. ریاضیدان آلمانی، آدولف هورویتز، در سال ۱۸۹۵ به‌طور مستقل پیشنهاد کرد که ضرایب چندجمله‌ای را در یک ماتریس مربعی، به نام ماتریس هورویتز، مرتب کند و نشان داد که چندجمله‌ای پایدار است، اگر فقط اگر دنباله دترمینان‌های زیرماتریس‌های اصلی آن مثبت باشد. این دو روش معادل هستند، با محک روث روش کارآمدتری برای محاسبهٔ دترمینان‌های هرویتز نسبت به محاسبه مستقیم آنها ارائه می‌دهد. چندجمله‌ای ارضاء کننده معیار راث-هرویتز چندجمله‌ای هرویتز نامیده می‌شود.

استفاده از الگوریتم اقلیدس

این معیار مربوط به قضیه راث-هورویتز است. از بیان آن قضیه، ما داریم $p-q=w(+\mathrm{\infty })-w(-\mathrm{\infty })$ که:

• $p$ تعداد ریشه‌های چندجمله‌ای $f(z)$ با قسمت حقیقی منفی است؛
• $q$ تعداد ریشه‌های چندجمله‌ای $f(z)$ با قسمت حقیقی مثبت (طبق قضیه ، $f$ قرار است هیچ ریشه ای روی خط موهومی نباشد) است.
• w(x) تعداد تغییرات زنجیره استورم تعمیم یافته است که از$P_0(y)$ و$P_1(y)$ بدست می‌آید (توسط تقسیمات اقلیدوسی متوالی) که$f(iy)=P_0(y)+i{P}_1(y)$ برای Y حقیقی است.

با توجه به قضیه اساسی جبر، هر چند جمله ای درجه n باید دارای n ریشه در صفحه مختلط باشد (یعنی برای یک ƒ بدون ریشه در خط موهومی، p + q = n) بنابراین، این شرط را داریم که ƒ یک چندجمله‌ای پایدار (هرویتز) است اگر و فقط اگر p - q = n باشد(اثبات در زیر آورده شده‌است) با استفاده از قضیه راث-هرویتز، می‌توان شرط p و q را با زنجیره استورم تعمیم یافته جایگزین کرد، که به نوبه خود شرط ضرایب ƒ را خواهد داد.

مثال مرتبه-بالاتر

وقتی بدست آوردن ریشه‌های چندجمله‌ای مشخصه مرتبه بالاتر دشوار است می‌توان از روش جدولی برای تعیین پایداری استفاده کرد. برای یک چندجمله‌ای درجه n ام

• $D(s)=a_n{s}^n+a_{n-1}{s}^{n-1}+\cdots +a_{1}s+a_0$

جدول دارای n + 1 ردیف و ساختار زیر است:

$a_n$	$a_{n-2}$	$a_{n-4}$	$\dots$
$a_{n-1}$	$a_{n-3}$	$a_{n-5}$	$\dots$
$b_1$	$b_2$	$b_3$	$\dots$
$c_1$	$c_2$	$c_3$	$\dots$
$⋮$	$⋮$	$⋮$	$\ddots$

که در آن عناصر $b_i$ و $c_i$ به شرح زیر محاسبه می‌شود:

• $b_i=\frac{a_{n-1}\times a_{n-2i}-a_n\times a_{n-(2i+1)}}{a_{n-1}}.$
• $c_i=\frac{b_1\times a_{n-(2i+1)}-a_{n-1}\times b_{i+1}}{b_1}.$

پس از اتمام، تعداد تغییرات علامت در ستون اول تعداد ریشه‌های غیر منفی خواهد بود.

۰٫۷۵	۱٫۵	۰	۰
-3	۶	۰	۰
۳	۰	۰	۰
۶	۰	۰	۰

در ستون اول، دو تغییر علامت وجود دارد (۰٫۷۵ → ۳−، و ۳ → ۳−)، بنابراین دو ریشه غیرمنفی وجود دارد که سیستم ناپایدار است.

معادله مشخصه سیستم سروو داده‌شده‌است توسط:

• $b_0{s}^4+b_1{s}^3+b_2{s}^2+b_{3}s+b_4=0$

$b_0$	$b_2$	$b_4$	۰
$b_1$	$b_3$	۰	۰
$\frac{b_1{b}_2-b_0{b}_3}{b_1}$	$\frac{b_1{b}_4-b_0\times 0}{b_1}$ = $b_4$	۰	۰
$\frac{(b_1{b}_2-b_0{b}_3)b_3-b_1^2{b}_4}{b_1{b}_2-b_0{b}_3}$	۰	۰	۰
$b_4$	۰	۰	۰

برای پایداری، تمام عناصر ستون اول آرایه راث باید مثبت باشند؛ بنابراین شرایطی که برای پایداری سیستم داده شده باید ارضاء شود به شرح زیر فراهم شود:

$b_1>0,b_1{b}_2-b_0{b}_3>0,(b_1{b}_2-b_0{b}_3)b_3-b_1^2{b}_4>0,b_4>0$

می‌بینیم که اگر

$(b_1{b}_2-b_0{b}_3)b_3-b_1^2{b}_4\ge 0$

سپس

$b_1{b}_2-b_0{b}_3>0$

ارضاء شده‌است.

• $s^4+6s^3+11s^2+6s+200=0$

جدول زیر را داریم:

۱	۱۱	۲۰۰	۰
6 1	6 1	۰	۰
10 1	200 20	۰	۰
-۱۹	۰	۰	۰
۲۰	۰	۰	۰

دو تغییر علامت وجود دارد این سیستم ناپایدار است، زیرا دارای دو قطبِ نیم-صفحه-راست و دو قطبِ نیم-صفحه-چپ است. از آنجا که یک ردیف صفر در جدول راث ظاهر نمی‌شود، سیستم نمی‌تواند قطب‌های jω داشته باشد.

• $s^6+2s^5+8s^4+12s^3+20s^2+16s+16=0.$

جدول زیر را داریم:

۱	۸	۲۰	۱۶
۲	۱۲	۱۶	۰
۲	۱۲	۱۶	۰
۰	۰	۰	۰

در چنین حالتی چندجمله‌ای کمکی $A(s)=2s^4+12s^2+16.$ است که باز هم برابر با صفر است. مرحله بعدی تفکیک معادله فوق است که چند جمله ای $B(s)=8s^3+24s^1.$ را ارائه می‌دهد. ضرایب ردیف حاوی صفر، اکنون "۸" و "۲۴" می‌شوند. روند آرایه راث با استفاده از این مقادیر که دارای دو نقطه در محور موهومی هستند پیش می‌رود. این دو نقطه در محور موهومی عامل اصلی پایداری حاشیه‌ای هستند.

جستارهای وابسته

• مهندسی کنترل
• مشتق از آرایه روت
• معیار پایداری نایکوئیست
• قضیه روث-هورویتس
• مکان ریشه
• تابع انتقال
• معیار لینارد-شیپر (گونه دیگری که به محاسبات کمتری احتیاج دارد)
• قضیه خاریتونوف (گونه دیگری برای ضرایب مجهول محدودشده در فواصل زمانی)
• معیار پایداری ژوری (قیاسی برای سیستم‌های LTI زمان گسسته)
• معیار پایداری بیستریتز (قیاسی برای سیستم‌های LTI زمان گسسته)الگو:Div col end

منابع

• Felix Gantmacher (J.L. Brenner translator) (1959) Applications of the Theory of Matrices, pp 177–80, New York: Interscience.
• الگو:Cite journal
• الگو:Cite book
• الگو:Cite book
• الگو:Cite journal

پیوند به بیرون

• فایل‌اجرایی متلب که محک راث-هرویتز را اجرا می‌کند
• اجرای آنلاین معیار راث-هرویتز