Authors
Department of Electrical Engineering, Imam Khomeini International University, Qazvin, Iran
Abstract
Keywords
1- مقدمه[1]
تقریب سیستمهای مرتبه بالا با مدلهای مرتبه پایین یکی از مسائل مهم در تئوری سیستم است. استفاده از مدلهای مرتبه کاهشی تحلیل، شبیهسازی و طراحی سیستم کنترلی را آسانتر میکند. همچنین بُعد پایین سیستم و کنترلکننده، پیادهسازی کنترلکننده، تعمیر و نگهداری آن را ارزان و ساده میکند.
دیویسون، کاهش مرتبة سیستم را نخستینبار در سال 1966 پیشنهاد کرد [1]. چن و همکاران در سال 1968، با استفاده از بسط نردبانی تابع تبدیل سیستم ابعاد وسیع رهیافتی مبتنی بر بسطهای حوزه فرکانسی برای کاهش مرتبه ارائه کردند [2]. ویلسون در سال 1970، با کمینهسازی انتگرال مربع خطا که در آن، خطا اختلاف میان پاسخ ضربة مدل ابعاد وسیع و مدل کاهش مرتبه یافتة آن بود، نخستین روش مبتنی بر بهینهسازی را در کاهش مرتبة سیستمهای SISO و MIMO ارائه کرد [3,4]. شاماش در سال 1974، با استفاده از تقریب پید[1] برای کاهش مرتبة سیستمهای پایدار استفاده کرد [5].
ال - عطار و ویدیاساگار در اواخر دهه 70، با کمینهکردن نرم و روشی برای کاهش مرتبة سیستمهای ابعاد وسیع ارائه کردند [6]. آبیناتا و همکاران [7] در اوایل دهه1980، با کمینهکردن معادلة خطا به فرم بستهای از جواب دست یافت. تا پیش از سال 1981 روشهای کاهش مرتبه عموماً بر پایة مفهوم مقادیر ویژه استوار بودند. در سال 1981 مور با ارائه روشی مبتنی بر مدل متعادلشده، انقلابی در کاهش مرتبه ایجاد کرد [8]. در این روش که به برش متعادل[2] معروف است، از مفهوم مدل متعادل و مفاهیم گرامیان کنترلپذیری و رؤیتپذیری برای تعیین مودهای کم اهمیتتر استفاده شده است. پس از مور، گلوور با ارائه مفهومی مشابه به نام تقریب بهینه هانکل[3] به کاهش مرتبة سیستمهای ابعاد وسیع پرداخت [9].
اندرسون و لیو در اواخر دهه 80 روشی بسیار دقیق به نام انحراف تکین[4] ارائه کردند که در فرکانسهای پایین، عملکرد بسیار مطلوبی داشت [10].
فلدمن و فروند در سال 1995 روشی بر پایة تقریب پید و مبتنی بر روش لانسوز[5] ارائه کردند که بسیاری از مشکلات محاسبة ممانها[6] را حل کرد [11]. یکی دیگر از روشهایی که در کاهش مرتبه بسیار به آن توجه شد، استفاده از الگوریتمهای تصویرسازی زیرفضای کریلوف[7] است [12]. در سال 1996 جمشیدی در [13] روشهای کاهش مرتبة سیستمهای ابعاد وسیع را بررسی و دستهبندی کرد. در اوایل قرن اخیر، بهبودهای بسیاری برای روش برش متعادل ارائه شدند [14-15].
در دهه اخیر، استفاده از الگوریتمهای تکاملی در کاهش مرتبه بسیار درخور توجه قرار گرفته است. چن و همکاران از الگوریتم تفاضلی[8] برای کاهش مرتبة سیستم استفاده کردند [16]. پارمار و همکاران در سال 2007 از الگوریتم انبوه ذرات[9] برای کاهش مرتبة سیستمها استفاده کردند [17]. استفاده از چندجملهایهای متعامد و الگوریتمهای تکاملی از دیگر روشهای کاهش مرتبة مدل است. در این روشها سیستم مرتبه کامل با توابع متعامد بسط داده شد و با کمینهسازی تابع هدف با الگوریتمهای هوشمند مدل مرتبه کاهشی تعیین شد [18-20]. در سالهای اخیر به کاهش مرتبة سیستمهای غیرخطی و سیستمهای دوخطی بسیار توجه شده است. در این روشها با تعمیم روشهای کاهش مرتبة سیستمهای خطی از قبیل برش متعادل [21-23]، بهینهسازی نرم H2 [24-25]، استفاده از بسط سری ولترا[10] [26] و روشهای مبتنی بر زیرفضای کریلوف [27-28]، کاهش مرتبة مدل بررسی شده است.
در این مقاله، با استفاده از کمینهسازی تابع معیار چندهدفه مجموع ITSE و نرم H∞ سیستم، کاهش مرتبة سیستمهای ناپایدار بررسی شده است. بهطور کلی در کاهش مرتبة سیستمهای ناپایدار دو رهیافت وجود دارد؛ رهیافت تجزیه و رهیافت انتقال [13]. در روش انتقال، مودهای سیستم ناپایدار بهگونهای انتقال داده میشوند که سیستم پایدار شود. سپس کاهش مرتبة سیستم پایدار بررسی میشود و در آخر، با عکس انتقال لحاظ شده، مدل ناپایدار کاهش مرتبه یافته به دست میآید؛ درحالیکه در روش تجزیه، مودهای پایدار و ناپایدار از یکدیگر تفکیک میشوند و سپس بخش پایدار کاهش مرتبه مییابد. با توجه به اینکه در روش انتقال مکان مودهای ناپایدار تغییر میکند، این روش نسبت به روش تجزیه محبوبیت کمتری دارد. در صورتی که در روش تجزیه، مودهای ناپایدار تغییر نمیکنند. در این مقاله، پس از تفکیک مودهای پایدار و ناپایدار با روش تجزیه، از مفهوم مقادیر تکین هانکل[11] برای تعیین تعداد مودهای کم اهمیتتر بخش پایدار استفاده شده است. در ادامه، ساختار ثابت مناسبی برای مدل مرتبه کاهشی در نظر گرفته میشود که پارامترهای صورت و مخرج این مدل مرتبه کاهشی مجهولاند. با کمینهسازی تابع مجموع وزنی انتگرال حاصلضرب زمان در مجذور خطا و اختلاف نرم H∞، پارامترهای مجهول سیستم مرتبه کاهشی تعیین میشوند که در آن، خطا اختلاف پاسخ پله بخش پایدار سیستم اصلی و سیستم مرتبه کاهشی است. این کمینهسازی با کمک الگوریتم خفاش[12] انجام میشود. برای ارضای شرط پایداری از محک روث[13] استفاده شده است که این شرط بهصورت قید در مسئلۀ بهینهسازی بیان میشود و درنتیجه، مسئلۀ بهینهسازی به مسئلۀ بهینهسازی مقید تبدیل میشود. در آخر، بخش ناپایدار تفکیکشده به بخش مرتبه کاهشی افزوده میشود. برای نشاندادن کارایی و دقت روش پیشنهادی، دو سیستم ناپایدار با استفاده از روش پیشنهادی کاهش مرتبه داده شده و با برخی روشهای کاهش مرتبه مقایسه شدهاند. از نوآوریهای این مقاله، استفاده از ترکیب دو معیار مبتنی بر حوزة فرکانس و حوزة زمان، یعنی معیارهای ITSEو H∞، بهعنوان معیار بهینهسازی چندهدفه است. این معیار پیشنهادی قابلیت اثرگذاری در هر دو حوزة زمان و فرکانس را دارد. همچنین، استفاده از الگوریتم خفاش در کاهش مرتبة مدل بهعنوان الگوریتم فرا ابتکاری از دیگر نوآوریهای روش پیشنهادی است که یکی از مهمترین ویژگیهای آن، سرعت بالای همگرایی آن نسبت به دیگر الگوریتمها است. در بخش دوم الگوریتم خفاش، معرفی و در بخش سوم، روش پیشنهادی ارائه شده است. در بخش چهارم شبیهسازی و نتایج ارائه میشوند و درنهایت، در بخش پنجم جمعبندی صورت میگیرد.
2- الگوریتم خفاش
الگوریتم خفاش براساس رفتار شنوایی خفاشها در مکانیابی طعمههایشان توسعه داده شده است. نخستینبار یانگ، الگوریتم بهینهسازی خفاش را معرفی کرد [29]. خفاشها یکسری پالسهای مافوق صوت را منتشر میکنند و به بازتاب این پالسها از اشیای اطراف گوش میدهند. پهنای باند موجهای مافوق صوت منتشرشده با توجه به نوع اشیا و گونهها تغییر میکند و با استفاده از هارمونیکها افزایش مییابد. این موجها با تأخیر زمانی و سطوح مختلفی از صوت بازتاب داده میشوند و باعث میشوند خفاش طعمة مدنظرش را مکان یابی و شکار کند. در سالهای اخیر، بهینهسازی با الگوریتم خفاش در علوم مهندسی بسیار شایان توجه قرار گرفته است. تنظیم بهینة پایدارساز سیستم قدرت غیرمتمرکز برای سیستم چند ماشین [30]، برنامهریزی شبکة انتقال [31]، طراحی بهینة پارامترهای کنترل پیشبین مدل برای کنترل فرکانس بار [32] ازجمله کاربردهای اخیر الگوریتم خفاشاند. خلاصهای از مراحل بهینهسازی الگوریتم خفاش در ادامه ارائه شده است [33].
گام 1. تمام خفاشها برای سنجش فاصله و دستهبندی میان طعمه و مانع از شنوایی استفاده میکنند.
گام 2. هر خفاش با سرعت در موقعیت ، با فرکانس ثابت ، طول موج متغیر و بلندی صوت برای شکار طعمه پرواز میکند. خفاش فرکانس پالس منتشرشدة خود را در محدودة تنظیم میکند و نرخ انتشار پالس را در محدوده براساس نزدیکی به هدف تغییر میدهد.
گام 3. فرکانس، بلندی صوت و نرخ انتشار پالس هر خفاش تغییر میکند.
گام 4. بلندی صوت خفاشها از مقدار بزرگ به حداقل مقدار ثابت تغییر میکند.
در طول فرایند بهینهسازی، موقعیت و سرعت بهصورت (1) تا (3) بهروز میشوند:
(1) |
|
(2) |
|
(3) |
موقعیت و سرعتهای بهترتیب موقعیت و سرعت هر خفاش iام در تکرار t و بردار تصادفی بهدستآمده از تابع توزیع یکنواخت است. همچین، بهترین مکان فعلی است که در هر تکرار پس از مقایسه با موقعیت خفاشها انتخاب میشود. در هر تکرار در جستجوی محلی، بهترین جواب، انتخاب و موقعیت جدید هر خفاش بهطور محلی با گام تصادفی بهصورت (4) بهروز میشود:
(4) |
عددی تصادفی و متوسط بلندی صوت تمام خفاشها در هر تکرار است. همچنین، بلندی صوت و نرخ پالس ارسالی r در هر تکرار بهصورت (5) تا (6) بهروز میشود:
(5) |
|
(6) |
بلندی صوت صفر نشان میدهد خفاش طعمه را یافته و بهطور موقت ارسال پالس را متوقف کرده است. هنگامی که تعداد تکرار به بینهایت میل کند، بلندی صوت به صفر میرسد و است. فلوچارت الگوریتم خفاش در شکل (1) نشان داده شده است [34].
شکل (1): فلوچارت الگوریتم خفاش
3- روش پیشنهادی
فرض کنید یک سیستم ابعاد وسیع ناپایدار باشد. ابتدا این سیستم ناپایدار به دو بخش پایدار و ناپایدار تفکیک شده است:
(7) |
اکنون بخش پایدار که از درجة n است را بهصورت (8) در نظر بگیرید:
(8) |
و ثابتهای معلوماند.
فرض کنید مدل مرتبه کاهشی پایدار بهصورت (9) باشد:
(9) |
و ثابتهای مجهول و r مرتبة مدل کاهش یافتهاند که r<n است. برای تعیین r از مفهوم مقادیر تکین هانکل استفاده شده است.
اگر گرامیان کنترلپذیری و رؤیتپذیری بهترتیب با P و Q نمایش داده شوند، با جذر مقادیر ویژه حاصلضرب P و Q، یعنی مقادیر تکین هانکل محاسبه میشوند. با استفاده از این تعریف میتوان آن دسته از متغیرهای حالت که اثر ضعیفی بر کنترلپذیری و رؤیتپذیری دارند را مشخص کرد و متعاقباً مرتبة سیستم کاهش یافته را به دست آورد.
هدف، تعیین پارامترهای مجهول سیستم مرتبه کاهشی پایدار است؛ بهگونهایکه سیستم مشخصههای سیستم را داشته باشد.
برای دستیابی به مدل مرتبه کاهشی با روش پیشنهادی، ابتدا بردار پاسخ پله سیستم ، تعیین شده است. سپس همانند رابطة (9) ساختاری ثابت با ضرایب صورت و مخرج مجهول، بخش پایدار مدل مرتبه کاهشی در نظر گرفته شده است. این پارامترهای مجهول با استفاده از الگوریتم خفاش تعیین خواهند شد. هدف بهینهسازی، یافتن بهترین پارامترها برای است. اگر و بهترتیب پاسخ پله سیستم و باشند، انتگرال حاصلضرب زمان در مجذور خطا بهصورت (10) تعریف میشود:
(10) |
همچنین، نرم H∞ بهصورت (11) تعریف میشود:
(11) |
که و میباشند.
با کمینهکردن تابع برازش مجموع وزنی معیار ITSE و نرم H∞ پارامترهای مدل مرتبه کاهشی تعیین خواهند شد:
(12) |
و ضرایب وزنی تابع برازشاند که بهصورت تجربی تعیین میشوند. با توجه به تأثیر بیشتر معیار ITSE نسبت به نرم H∞ ، و انتخاب شدهاند.
رابطة (12) معیاری ترکیبی است که در هر دو حوزة زمان و فرکانس اثرگذار است. با کمینهسازی معیار (12) خطا در هر دو حوزة زمان و فرکانس شایان توجه قرار خواهد گرفت.
بهدلیل اینکه بخش پایدار سیستم ابعاد وسیع ناپایدار باید کاهش داده شود، مدل مرتبه کاهشی نیز باید پایدار باشد؛ بنابراین، روش پیشنهادی باید پایداری مدل مرتبه کاهشی را نیز تضمین کند. به همین منظور، از محک روث برای تعیین شرط پایداری بهصورت زیر استفاده شده است [19]:
مخرج مدل مرتبه کاهشی (9) را بهصورت زیر میتوان نشان داد [35]:
(13) |
ضرایب دو سطر اول آرایة روث بههمراه عناصر ستون اول را بهصورت میتوان زیر نمایش داد:
(14) |
k برای r زوج، 1 و برای r فرد صفر است.
با مقایسة درایههای سطر اول با و سطر دوم با ، رابطة زیر به دست میآید:
(15) |
با قراردادن روابط بالا در مخرج مدل مرتبه کاهشی، رابطة (15) به دست میآید؛ بنابراین، شرط لازم و کافی برای اینکه تمام قطبهای مدل مرتبه کاهشی اکیداً در نیمصفحة چپ باشند، عبارت است از:
(16) |
و در نتیجه:
(17) |
بنابراین، برای اینکه مدل مرتبه کاهشی پایدار باشد، پارامترهای مدل مرتبه کاهشی، با کمینهکردن (12) و رعایت شرط (17) تعیین میشوند که مسئلۀ بهینهسازی به مسئلۀ بهینهسازی مقید تبدیل میشود. به عبارت دیگر، مدل مرتبه کاهشی با کمینهکردن تابع برازش زیر به دست میآید:
(18) |
با کمینهکردن تابع برازش مقید (18) با استفاده از الگوریتم خفاش، پارامترهای بهگونهای تعیین میشوند که مشخصههای آن با مشخصههای نزدیک به هم باشند.
پس از تعیین پارامترهای ، بخش ناپایدار را بهصورت (19) به بخش پایدار کاهش مرتبه یافته اضافه میکنیم تا مدل مرتبه کاهشی نهایی به دست آید:
(19) |
گفتنی است در کاهش مرتبة سیستمهای ناپایدار بهطور کلی دو رهیافت وجود دارد؛ رهیافت تجزیه و رهیافت انتقال. در روش انتقال، مودهای سیستم ناپایدار بهگونهای انتقال داده میشوند که سیستم پایدار شود. سپس کاهش مرتبة سیستم پایدار بررسی میشود و در آخر، با عکس انتقال لحاظشده، مدل ناپایدار کاهش مرتبه یافته به دست میآید؛ با این حال، با توجه به اینکه در روش انتقال مکان مودهای ناپایدار تغییر میکند، این روش نسبت به روش تجزیه از محبوبیت کمتری برخوردار است؛ در صورتی که در روش تجزیه مودهای ناپایدار تغییر نمیکنند.
روش پیشنهادی در گامهای زیر خلاصه میشود:
گام اول: تفکیک سیستم ابعاد وسیع ناپایدار به دو بخش پایدار و ناپایدار؛
گام دوم: در نظر گرفتن بخش پایدار برای کاهش مرتبه؛
گام سوم: انتخاب مرتبة مناسب برای مدل مرتبه کاهشی پایدار با استفاده از مفهوم مقادیر تکین هانکل؛
گام چهارم: انتخاب ساختار ثابت مناسبی همانند (9) برای مدل مرتبه کاهشی که در آن و پارامترهای مجهولاند و باید تعیین شوند؛
گام پنجم: استفاده از الگوریتم خفاش برای تعیین پارامترهای مجهول مدل مرتبه کاهشی پایدار با کمینهسازی تابع برازش مقید (18)؛
گام ششم: افزودن بخش ناپایدار به بخش کاهش مرتبه پایدار برای دستیابی به مدل مرتبه کاهشی ناپایدار.
4- شبیهسازی و نتایج
در این بخش برای تأیید عملکرد و توانایی روش پیشنهادی، سه سیستم ناپایدار ابعاد وسیع کاهش مرتبه داده شدهاند. بهمنظور فهم بهتر و آسانتر روش پیشنهادی، مثال اول بهصورت مرحله به مرحله توضیح داده شده است.
مثال 1: سیستم ناپایدار از مرتبه 5 را در نظر بگیرید که در [36] ارائه شده است:
(20) |
گام 1: تفکیک سیستم ابعاد وسیع ناپایدار (20) به دو بخش پایدار و ناپایدار.
(21) |
گام 2: در نظر گرفتن بخش پایدار رابطة (20) بهعنوان مدلی که قصد داریم مرتبة آن را کاهش دهیم.
(22) |
گام 3: رسم مقادیر تکین بخش پایدار سیستم ابعاد وسیع برای تعیین مرتبة مناسب سیستم مرتبه کاهشی.
شکل (2): مقادیر تکین هانکل بخش پایدار سیستم نمونة 1
همانطور که مشاهده میشود مودهای 2 و 3 نسبت به مود 1 انرژی بسیار کمتری دارند و درنتیجه، از مودهای 2 و 3 میتوان صرفنظر کرد؛ بنابراین، مرتبة مدل کاهش مرتبه یافته 1 است.
گام 4: انتخاب ساختار مناسبی با ضرایب مجهول با توجه به مرتبة بهدستآمده.
(23) |
گام 5: کمینهسازی تابع برازش (18) با استفاده از الگوریتم خفاش برای تعیین پارامترهای مجهول رابطة (23).
(24) |
گام 6: افزودن بخش ناپایدار به بخش پایدار کاهش مرتبه یافته برای دستیابی به مدل کاهش مرتبه کل.
(25) |
شکل (3): پاسخ پله بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و برخی
روشهای کلاسیک کاهش مرتبه
برای نشاندادن دقت و کارایی روش پیشنهادی، پاسخ پله بخش پایدار و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی، با دیگر روشهای رایج کاهش مرتبه از قبیل برش متعادل (BT) [37] و روش تقریب هانکل(HSV) [37] مقایسه شده است. با توجه به شکلهای (3) و (4) مشاهده میشود مشخصههای مدل بهدستآمده با روش پیشنهادی نسبت به روشهای کلاسیک کاهش مرتبه به سیستم اصلی بسیار نزدیکتر است. در جدول (1) برخی مشخصههای مهم سیستم از قبیل مقدار حالت ماندگار، زمان صعود، زمان نشست و بیشینة فراجهش مقایسه شدهاند. همچنین، معیار انتگرال حاصلضرب زمان در مربع خطا و نرم بررسی شدند که تابع هدف در کمینهسازی در نظر گرفته شدهاند.
شکل (4): نمودار اندازة خطا برای بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و روشهای برش متعادل و روش تقریب هانکل برای مثال 1
جدول (1): مقایسة مشخصههای بخش پایدار سیستم ابعاد وسیع و مدل کاهش مرتبه یافتة آن برای مثال 1
|
مقدار حالت ماندگار به ورودی پله واحد |
بیشینه فراجهش (%) |
زمان صعود (ثانیه) |
زمان نشست (ثانیه) |
انتگرال حاصلضرب زمان در مربع خطا (ITSE) |
نرم |
سیستم اصلی |
86/1 |
0 |
9902/2 |
1970/5 |
- |
8571/1 |
روش پیشنهادی |
86/1 |
0 |
8142/2 |
0110/5 |
7004/23 |
8571/1 |
برش متعادل |
64/1 |
0 |
7554/1 |
1258/3 |
8/1142 |
6403/1 |
هانکل |
59/1 |
0 |
2662/2 |
0353/4 |
7/1752 |
5930/1 |
با توجه به شکلهای (3) و (4) و همچنین، جدول (1) مشاهده میشود روش پیشنهادی نسبت به دیگر روشهای کاهش مرتبه، عملکرد مشابهتر و مشخصههای نزدیکتری به بخش پایدار سیستم اصلی دارد.
در ادامه، با افزدون بخش ناپایدار به بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و متعاقباً روشهای کلاسیک کاهش مرتبه از قبیل BT و HSV، پاسخ پله سیستم ناپایدار مقایسه شده است.
شایان ذکر است که در [36] مرتبة سیستم ناپایدار به 4 کاهش داده شده است؛ درحالیکه در روش پیشنهادی مرتبة مدل کاهش مرتبه یافته 3 شد و از همین رو در مقایسة عملکرد روشهای کاهش مرتبه، روش پیشنهادشده در [36] بررسی نشد.
مثال 2: تابع انتقال ناپایدار و ناکمینة فاز از مرتبة 4 ارائهشده در [16] را در نظر بگیرید:
(26) |
شکل (5): پاسخ پله سیستم ناپایدار ابعاد وسیع و مدلهای کاهش مرتبه یافته برای مثال 1
با توجه به روش پیشنهادی باید سیستم ناپایدار (26) را به دو بخش پایدار و ناپایدار تفکیک کنیم.
(27) |
سپس مقادیر تکین هانکل بخش پایدار سیستم ابعاد وسیع (26) را رسم میکنیم تا ابعاد بخش مرتبه کاهشی پایدار تعیین شود.
شکل (6): مقادیر تکین بخش پایدار سیستم مثال 2
با توجه به شکل (6) مشاهده میشود مود 3، انرژی کمی دارد و از آن چشمپوشی میشود؛ بنابراین، مرتبه مدل پایدار کاهشی 2 انتخاب میشود.
با انتخاب ساختار ثابت مناسب با ضرایب مجهول و اعمال گامهای 4 تا 5 به مدل پایدار مرتبه کاهشی (28) دست مییابیم:
(28) |
پاسخ پله بخش پایدار سیستم ابعاد وسیع و مدل مرتبه کاهشی متناظر بهدستآمده با روش پیشنهادی، روشهای BT و HSV و روش پیشنهادی [16] در شکل (7) نشان داده شدهاند. همچنین در شکل (8)، نمودار اندازة خطا برحسب زمان نشان داده شده است که با توجه به آن، بهوضوح مشاهده میشود روش پیشنهادی نسبت به دیگر روشهای کاهش مرتبه کلاسیک خطای بسیار کمتری دارد.
شکل (7): پاسخ پله بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و برخی روشهای کلاسیک کاهش مرتبه برای مثال 2
شکل (8): نمودار اندازة خطا برای بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و روشهای برش متعادل (BT)، روش تقریب هانکل (HSV) و روش [16] برای مثال 2
در جدول (2)، مشخصههای بخش پایدار سیستم اصلی و مدلهای کاهش مرتبه یافتة متناظر با بخش پایدار مقایسه شدهاند. در ادامه، بخش ناپایدار تفکیکشده در رابطة (27)، به بخش پایدار کاهش مرتبه یافته افزوده میشود تا مدل کاهش مرتبه یافتة نهایی تعیین شود.
(29) |
با توجه به جدول (2) و شکلهای (7-8) مشاهده میشود روش پیشنهادی نسبت به دیگر روشهای کاهش مرتبه، دقت و عملکرد بهتری دارد.
جدول (2): مقایسة مشخصههای بخش پایدار سیستم ابعاد وسیع و مدل کاهش مرتبه یافتة آن برای مثال 2
|
مقدار حالت ماندگار به ورودی پله واحد |
بیشینه فراجهش (%) |
زمان صعود (ثانیه) |
زمان نشست (ثانیه) |
انتگرال حاصلضرب زمان در مربع خطا (ITSE) |
نرم |
سیستم اصلی |
07/6 |
0 |
2775/0 |
5211/0 |
- |
0701/6 |
روش پیشنهادی |
07/6 |
0 |
1816/0 |
7715/0 |
1262/85 |
0701/6 |
برش متعادل |
13/7 |
0 |
1313/1 |
5141/2 |
29894 |
1314/7 |
هانکل |
13/7 |
0 |
5002/0 |
0238/1 |
31264 |
1314/7 |
روش [16] |
09/7 |
0 |
2052/0 |
4044/0 |
30077 |
0942/7 |
مثال 3: سیستم ناپایدار از مرتبة 7 ارائهشده توسط یانگ در [38] را در نظر بگیرید.
(30) |
با تفکیک سیستم (30) به دو بخش پایدار و ناپایدار، بخش پایدار و ناپایدار بهترتیب بهصورت (31) و (32) خواهند بود.
(31) |
|
(32) |
در این بخش، مقادیر تکین هانکل رسم شده است تا مرتبة سیستم مرتبه کاهشی تعیین شود.
شکل (9): مقادیر تکین بخش پایدار سیستم مثال 2
با توجه به شکل (9) مشاهده میشود مرتبة مناسب بخش پایدار مرتبه کاهشی 2 است. با کمینهسازی تابع هدف ارائهشده در (18) با استفاده از الگوریتم خفاش بهمدل مرتبه کاهشی از رابطة (31) دست مییابیم.
(33) |
برای بررسی دقت و کارایی روش پیشنهادی، در شکل (10) پاسخ پله بخش پایدار بهدستآمده با روش پیشنهادی و برخی روشهای کلاسیک از قبیل BT و HSV و روش [38] مقایسه شدهاند. برای بررسی بیشتر نیز نمودار اندازة خطا برحسب زمان در شکل (11) نشان داده شده است.
شکل (10): پاسخ پله بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و برخی روشهای
کلاسیک کاهش مرتبه برای مثال 3
شکل (11): نمودار اندازة خطا برای بخش پایدار سیستم اصلی و بخش پایدار کاهش مرتبه یافته با روش پیشنهادی و روشهای
برش متعادل (BT)، روش تقریب هانکل (HSV) و روش [38] برای مثال 3
جدول (3): مقایسة مشخصههای بخش پایدار سیستم ابعاد وسیع و مدل کاهش مرتبه یافتة آن برای مثال 3
|
مقدار حالت ماندگار به ورودی پله واحد |
بیشینه فراجهش (%) |
زمان صعود (ثانیه) |
زمان نشست (ثانیه) |
انتگرال حاصلضرب زمان در مربع خطا (ITSE) |
نرم |
سیستم اصلی |
00241/0- |
4371/171 |
1068/0 |
9207/8 |
- |
0246/0 |
روش پیشنهادی |
00235/0- |
3890/112 |
1267/0 |
6088/11 |
0026/0 |
0248/0 |
برش متعادل |
00329/0- |
8276/86 |
1550/0 |
4591/8 |
0074/0 |
0240/0 |
هانکل |
00351/0- |
8559/84 |
1584/0 |
9240/7 |
0114/0 |
0239/0 |
روش [38] |
0032/0- |
6138/87 |
1560/0 |
6577/8 |
0062/0 |
0232/0 |
در جدول (3)، مشخصههای اصلی مدلهای مرتبه کاهشی بهدستآمده با بخش پایدار سیستم ابعاد وسیع مقایسه شدهاند. با توجه به شکلهای (10) و (11) و جدول (3) بهوضوح مشاهده میشود نتایج بهدستآمده با روش پیشنهادی نسبت به دیگر روشهای کاهش مرتبه، دقت و کارایی بهتری دارند. در ادامه، با افزودن بخش ناپایدار به بخش پایدار مرتبه کاهشی، مدل نهایی سیستم مرتبه کاهشی به دست میآید. برای بررسی کارایی و دقت مدل مرتبه کاهشی نهایی، در شکل (12) پاسخ پله سیستم ابعاد وسیع و مدلهای نهایی بهدستآمده با روش پیشنهادی و برخی روشهای کاهش مرتبه نشان داده شدهاند. نتایج نشان میدهند روش پیشنهادی بیشترین تشابه با مدل ناپایدار اصلی را دارد.
شکل (12): پاسخ پله سیستم ناپایدار ابعاد وسیع و مدلهای کاهش مرتبه یافته برای مثال 3
5- نتیجهگیری
در این مقاله، روشی مبتنی بر کمینهسازی تابع دو هدفه مجموع وزنی معیارهای ITSE و نرم H∞ با استفاده از الگوریتم خفاش، برای کاهش مرتبة سیستمهای ابعاد وسیع ناپایدار ارائه شد. پس از تفکیک سیستم ابعاد وسیع ناپایدار به دو بخش پایدار و ناپایدار، با بررسی مقادیر تکین هانکل بخش پایدار، مرتبة بخش پایدار مدل مرتبه کاهشی تعیین شده است. با انتخاب ساختاری ثابت با ضرایب مجهول و کمینهسازی تابع برازش با استفاده از الگوریتم خفاش، با توجه به محک روث بهعنوان قید، مدل مرتبه کاهشی بخش پایدار تعیین شده است. در آخر، با افزودن بخش ناپایدار به بخش کاهش مرتبه یافته، مدل کاهش مرتبه یافتة نهایی به دست آمده است. توانمندی و کارایی روش پیشنهادی با مقایسة نتایج بهدستآمده با برخی روشهای رایج کاهش مرتبه از قبیل برش متعادل و روش تقریب هانکل نشان داده شده است.
[1]تاریخ ارسال مقاله: 15/11/1396
تاریخ پذیرش مقاله: 18/08/1398
نام نویسندۀ مسئول: حسن زرآبادیپور
نشانی نویسندۀ مسئول: ایران، قزوین، دانشگاه بینالمللی امام خمینی قزوین، گروه برق