یک الگوریتم برای کنترل بهينه يك كلاس از سيستم‌هاي خطي متغير با زمان تأخیری با رویکرد کاهش زمان محاسبات و افزایش سرعت در مسائل مهندسی
الموضوعات :مهدی یوسفی طبری 1 , زهرا رحمانی 2 , علی وحیدیان کامیاد 3 , سیدجلیل ساداتی 4
1 - دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل
2 - دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل
3 - دانشکده علوم ریاضی، دانشگاه فردوسی مشهد
4 - دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل
الکلمات المفتاحية: کنترل بهینه, سیستمهای تأخیری, اصل حداکثر پونتریاگین,
ملخص المقالة :
سیستمهای تأخیر زمانی در چند دهه اخیر بسیار مورد توجه قرار گرفته‌اند و بسیاری از آنها در سیستم‌ها و شاخه‌های مختلف علوم مانند مهندسی، شیمی، فیزیک و مدل‌های بیماری ظاهر می‌شوند. وجود تأخیر، تحلیل و کنترل چنین سیستم‌هایی را بسیار پیچیده‌تر می‌کند. استفاده از اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی منجر به یک مسأله مقدار مرزی می‌شود که شامل هر دو شرایط تأخیر و تقدم است. در این مقاله، یک مسأله کنترل بهینه با تأخیر زمانی را در نظر می‌گیریم. در بخش اول ابتدا با استفاده از اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی، شرایط بهینه لازم را برای این مسأله به‌دست می‌آوریم و سپس الگوریتمی جدید برای حل عددی این مسأله ارائه می‌گردد که بر پایه یک تقریب برای مشتقات و درون‌یابی خطی برای جملات تأخیر است. نهایتاً معادلات حاصل به یک مسأله برنامه‌ریزی خطی تبدیل می‌شوند که می‌توان آنها را به‌صورت عددی حل نمود. کارايي روش پيشنهادی با شبیه‌سازی عددی مورد ارزیابی قرار می‌گیرد.
[1] M. R. Hestenes, Calculus of Variations and Optimal Control Theory, Wiley, 1966.
[2] A. E. Bryson, "Optimal control-1950 to 1985," IEEE Control Systems Magazine, vol. 16, no. 3, pp. 26-33, Jun. 1996.
[3] M. Jamshidi and C. Wang, "A computational algorithm for large-scale nonlinear time-delay systems," IEEE Trans. on Systems, Man, Cybernetics, vol. 14, no. 1, pp. 2-9, Jan.-Feb. 1984.
[4] M. Malek-Zavarei and M. Jamshidi, Time-Delay Systems: Analysis, Optimization and Applications, Elsevier Science Inc., 1987.
[5] H. J. Sussmann and J. C. Willems, "300 years of optimal control: from the brachystochrone to the maximum principle," IEEE Control Systems Magazine, vol. 17, no. 3, pp. 32-44, 1997.
[6] G. Kharatishdi, "The maximum principle in the theory of optimal processes with time lags," Dokl. Akad. Nauk SSSR, vol. 136, no. 1, pp. 39-43, 1961.
[7] D. Eller, J. Aggarwal, and H. Banks, "Optimal control of linear time-delay systems," IEEE Trans. on Automatic Control, vol. 14, no. 6, pp. 678-687, Dec. 1969.
[8] K. Palanisamy and R. G. Prasada, "Optimal control of linear systems with delays in state and control via Walsh functions," IEE Proc. D-Control Theory and Applications, vol. 130, no. 6, pp. 300-312, Nov. 1983.
[9] K. Inoue, H. Akashi, K. Ogino, and Y. Sawaragi, "Sensitivity approaches to optimization of linear systems with time delay," Automatica, vol. 7, no. 6, pp. 671-679, Nov. 1971.
[10] J. Banas and A. Vacroux, "Optimal piecewise constant control of continuous time systems with time-varying delay," Automatica, vol. 6, no. 6, pp. 809-811, Nov. 1970.
[11] H. R. Marzban and M. Razzaghi, "Optimal control of linear delay systems via hybrid of block-pulse and Legendre polynomials," J. of the Franklin Institute, vol. 341, no. 3, pp. 279-293, May 2004.
[12] L. Y. Lee, "Numerical solution of time‐delayed optimal control problems with terminal inequality constraints," Optimal Control Applications and Methods, vol. 14, no. 3, pp. 203-210, Jul./Sept. 1993.
[13] ها. چهکندی نژاد، م. فرشاد و ر هاونگی، "ارائه یک روش جدید به منظور تخمین برخط تأخیر زمانی در سیستمهای SISO-LTI با تأخیر زمانی متغیر با زمان و نامعلوم در ورودی کنترلی،" نشریه مهندسی برق و مهندسی کامپیوتر ایران، الف- مهندسی برق، سال 18، شماره 1، صص. 44-36، بهار 1399.
[14] H. Marzban and S. Hoseini, "Solution of linear optimal control problems with time delay using a composite Chebyshev finite difference method," Optimal Control Applications and Methods, vol. 34, no. 3, pp. 253-274, May/Jun. 2013.
[15] X. T. Wang, "Numerical solutions of optimal control for time delay systems by hybrid of block-pulse functions and Legendre polynomials," Applied Mathematics, vol. 184, no. 2, pp. 849-856, 15 Jan. 2007.
[16] A. Jajarmi and M. Hajipour, "An efficient finite difference method for the time‐delay optimal control problems with time‐varying delay," Asian J. of Control, vol. 19, no. 2, pp. 554-563, Mar. 2017.
[17] I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Elsevier, 1998.
[18] O. P. Agrawal and D. Baleanu, "A Hamiltonian formulation and a direct numerical scheme for fractional optimal control problems," J. of Vibration and Control, vol. 13, no. 9-10, pp. 1269-1281, 2007.
[19] R. L. Burden, J. D. Faires, and A. Reynolds, Numerical Analysis Prindle, Weber & Schmidt, pp. 28-34, 1985.
[20] H. Banks and J. Burns, "Hereditary control problems: numerical methods based on averaging approximations," SIAM J. on Control Optimization, vol. 16, no. 2, pp. 169-208, 1978.
[21] A. Jajarmi and D. Baleanu, "Suboptimal control of fractional-order dynamic systems with delay argument," J. of Vibration and Control, vol. 24, no. 12, pp. 2430-2446, 2017.
[22] L. Moradi, F. Mohammadi, and D. Baleanu, "A direct numerical solution of time-delay fractional optimal control problems by using Chelyshkov wavelets," J. of Vibration and Control, vol. 25, no. 2, pp. 310-324, 2019.
[23] S. Sabermahani, Y. Ordokhani, and S. A. Yousefi, "Fractional-order Lagrange polynomials: an application for solving delay fractional optimal control problems," Trans. of the Institute of Measurement and Control, vol. 41, no. 11, pp. 2997-3009, 2019.
[24] H. R. Marzban and F. Malakoutikhah, "Solution of delay fractional optimal control problems using a hybrid of block-pulse functions and orthonormal Taylor polynomials," J. of the Franklin Institute, vol. 356, no. 15, pp. 8182-8215, Oct. 2019.
[25] N. Haddadi, Y. Ordokhani, and M. Razzaghi, "Optimal control of delay systems by using a hybrid functions approximation," J. of Optimization Theory and Applications, vol. 153, no. 2, pp. 338-356, 12 Oct. 2012.
[26] K. Palanisamy, K. Balachandran, and R. Ramasamy, "Optimal control of linear time-varying delay systems via single-term Walsh series," IEE Proc. D-Control Theory and Applications, vol. 135, no. 4, pp. 332-332, Jul. 1988.
[27] S. Dadebo and R. Luus, "Optimal control of time‐delay systems by dynamic programming," Optimal Control Applications and Methods, vol. 13, no. 1, pp. 29-41, Jan./Mar. 1992.
[28] J. R. Ockendon and A. B. Tayler, "The dynamics of a current collection system for an electric locomotive," Proc. of the Royal Society A: Mathematical Physical Sciences, vol. 322, no. 1551, pp. 447-468, 4 May 1971. 1971.
[29] F. Ghomanjani, M. H. Farahi, and A. V. Kamyad, "Numerical solution of some linear optimal control systems with pantograph delays," IMA J. of Mathematical Control and Information, vol. 32, no. 2, pp. 225-243, Jun. 2015.
[30] N. Ghaderi and M. H. Farahi, "The numerical solution of some optimal control systems with constant and pantograph delays via bernstein polynomials," Iranian J. of Mathematical Sciences Informatics, vol. 15, no. 2, pp. 163-181, 2020.
[31] M. Fatehi, M. Vali, and M. Samavat, "State analysis and optimal control of linear time-invariant scale systems using the legendre wavelets," Canadian J. on Automation, Control & Intelligent Systems, vol. 3, no. 1, pp. 1-7, May 2012.
[32] M. Fatehi, M. Samavat, M. Vali, and F. Khaleghi, "State analysis and optimal control of linear time-invariant scaled systems using the Chebyshev wavelets," Contemporary Engineering Sciences, vol. 5, no. 2, pp. 91-105, 2012.
178 نشریه مهندسی برق و مهندسی کامپیوتر ایران، الف- مهندسی برق، سال 21، شماره 3، پاییز 1402
مقاله پژوهشی
یک الگوریتم برای کنترل بهينه يك كلاس از سيستم‌هاي خطي
متغير با زمان تأخیری با رویکرد کاهش زمان محاسبات
و افزایش سرعت در مسائل مهندسی
مهدی یوسفی طبری، زهرا رحمانی، علی وحیدیان کامیاد و سید جلیل ساداتی
چکیده: سیستمهای تأخیر زمانی در چند دهه اخیر بسیار مورد توجه قرار گرفته‌اند و بسیاری از آنها در سیستم‌ها و شاخه‌های مختلف علوم مانند مهندسی، شیمی، فیزیک و مدل‌های بیماری ظاهر می‌شوند. وجود تأخیر، تحلیل و کنترل چنین سیستم‌هایی را بسیار پیچیده‌تر می‌کند. استفاده از اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی منجر به یک مسأله مقدار مرزی می‌شود که شامل هر دو شرایط تأخیر و تقدم است. در این مقاله، یک مسأله کنترل بهینه با تأخیر زمانی را در نظر می‌گیریم. در بخش اول ابتدا با استفاده از اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی، شرایط بهینه لازم را برای این مسأله به‌دست می‌آوریم و سپس الگوریتمی جدید برای حل عددی این مسأله ارائه می‌گردد که بر پایه یک تقریب برای مشتقات و درون‌یابی خطی برای جملات تأخیر است. نهایتاً معادلات حاصل به یک مسأله برنامه‌ریزی خطی تبدیل می‌شوند که می‌توان آنها را به‌صورت عددی حل نمود. کارايي روش پيشنهادی با شبیه‌سازی عددی مورد ارزیابی قرار می‌گیرد.
کلیدواژه: کنترل بهینه، سیستمهای تأخیری، اصل حداکثر پونتریاگین.
1- مقدمه
مسایل کنترل بهینه به‌طور طبیعی در حوزه‌های مختلف علوم، مهندسی و ریاضیات به وجود می‌آیند. کارهای قابل توجهی در زمینه کنترل بهینه سیستمهای دینامیکی مرتبه صحیح انجام گردیده [1] و [2] و کنترل سیستمهای با تأخیر زمانی یکی از زمینه‌های مهم تحقیقاتی بوده است. تأخیر، غالباً در سیستمهای بیولوژیکی، شیمیایی، الکترونیکی و حمل‌ونقل رخ می‌دهد [3]. سیستمهای تأخیر زمانی دسته بسیار مهمی از سیستم‌ها هستند که کنترل و بهینه‌سازی آنها مورد توجه بسیاری از محققین است [4]. عمومي‌شدن تئوري کنترل بهينه کلاسیک به سال 1950 همراه با کاربردهاي آن در صنايع هوافضا در جريان جنگ سرد برمی‌گردد [5]. استفاده از اصل ماکسیمم پونتریاگین برای بهینه‌سازی سیستمهای کنترل با تأخیر زمانی- همان طور که در [6] مشخص شد- منجر به یک مسأله مقدار مرزی دونقطه‌ای می‌شود که شامل هر دو شرط تأخیر و تقدم است و راه حل دقیق آن به جز در موارد بسیار خاص، بسیار دشوار است. در [7] روشی پیشنهاد شده که شامل حل مجموعه‌ای از معادلات دیفرانسیل با مشتقات جزئی است که به‌طور کلی از نظر محاسباتی زمان‌بر است. مرجع [8] الگوریتم جدیدی با استفاده از توابع والش ارائه می‌دهد که در آن از ماتریس‌های عملیاتی والش برای تقریب جملات تأخیر و تقدم استفاده می‌شود. در [9] یک رهیافت حساسیت برای به‌دست‌آوردن کنترل زیربهینه برای سیستمهای خطی با تأخیر کوچک در حالت پیشنهاد شده است. آنها کنترل را در سری مک‌لارن در تأخیر گسترش دادند و ضرایب سری را از حل مسائل ساده مقدار مرزی دونقطه‌ای به‌دست آوردند. در [10] شرایط لازم برای کنترل بهینه سیستمهای تأخیری قطعه‌ای ثابت و پیوسته با زمان ارائه شده است. با استفاده از هامیلتونین معمولی، اثر تأخیر در تابع هزینه، گنجانده و یک مسأله مقدار نهایی ثابت برای یک سیستم فرموله شده است. در [11] روشی برای یافتن کنترل بهینه یک سیستم تأخیری خطی متغیر با زمان با معیار عملکرد مرتبه دوم مورد بحث قرار گرفته و ویژگی‌های توابع ترکیبی که از توابع بلوک پالس به اضافه چندجمله‌ای لژاندر تشکیل شده‌اند، ارائه گردیده است. در [12] از متغیرهای کمک حالت و تقریب پده2 برای تبدیل مسأله به یک مسأله بدون تأخیر زمانی استفاده شده که مسأله را می‌توان با یک الگوریتم بهینه‌سازی توسعه‌یافته حل نمود. با جايگزين‌كردن عامل تأخير زماني به وسيله يك تقريب پده و با استفاده از يك ترفند در [13]، ابتدا پارامتر تأخير زماني به يك متغير حالت از سيستم، تبديل و سپس از تئوري فيلتر كالمن براي طراحي يك تخمينگر حالت استفاده مي‌شود. در [14] یک روش تفاضل محدود ترکیبی چبیشف برای یافتن راه حل کنترل بهینه سیستمهای تأخیر زمانی با معیار عملکرد مرتبه دوم، معرفی و اعمال شده که این روش، توسعه‌ای از طرح تفاضل محدود چبیشف است. در [15] با استفاده از ویژگی‌های عملیاتی توابع بلوک پالس عمومی و چندجمله‌ای لژاندر، مسأله کنترل بهینه برای سیستمهای تأخیر زمانی به یک سیستم معادلات جبری تبدیل می‌شود. در [16] با استفاده از اصل ماکسیمم پونتریاگین، مسأله کنترل بهینه تأخیر زمانی اولیه، ابتدا به یک مسأله مقدار مرزی دونقطه‌ای تبدیل می‌شود که شامل هر تأخیر و تقدم است. سپس با استفاده از فرمول تفاضل محدود مرتبه دوم و درون‌یابی به ترتیب برای مشتقات مرتبه اول و عبارات تأخیری به سیستم معادلات جبری خطی تبدیل می‌شود.
وجود تأخیر زمانی در سیستمهای تأخیری باعث می‌شود که آنالیز و کنترل آنها بیش از پیش پیچیده گردد؛ بنابراین ارائه تکنیک‌های جدید تحلیلی و عددی برای حل مسأله کنترل بهینه سیستمهای دینامیکی تأخیردار، يکي از زمينه‌هاي تحقيقاتي فعال و مشکل در تئوري کنترل است. همچنین مقالات ارائه‌شده در زمینه مسأله کنترل بهینه سیستمهای تأخیری عموماً راه حل‌هایی را با مشکلات محاسباتی جهت حل معادلات نهایی به صورت عددی معرفی کرده‌اند.
هدف اصلي این تحقیق، ارائه یک روش طراحي کنترل بهينه براي دسته‌اي از سيستم‌هاي تأخیری است. اگرچه این مسأله توسط محققین مورد بررسی قرار گرفته و روش‌هایی برای حل آن ارائه شده است، اما همچنان تحقیقاتی برای ارائه راه حل‌های مؤثر و کارآمد برای حصول دقت بیشتر، کاستن از پیچیدگی و حجم محاسبات و افزایش سرعت برای رسیدن به پاسخ بهینه با توجه به تنوع سیستمهای کنترلی و کاربردهای آن در حال انجام است.
تا کنون گزارشی برای حل‌کردن مسأله کنترل بهینه تأخیری بر اساس روش پیشنهادی در سیستمهای تأخیری ارائه نشده است. به همین منظور هدف و نوآوری این پژوهش، ارائه الگوریتمی جدید بر مبنای تقریب مشتق و درون‌یابی خطی جهت حل مسأله کنترل بهینه تأخیری برای دسته‌ای از این سیستم‌ها است. برای این منظور، ابتدا شرایط لازم بهینگی برای مسأله کنترل بهینه تأخیری را بیان نموده و سپس یک الگوریتم جدید برای حل عددی این مسأله ارائه می‌گردد. اساس این الگوریتم بر پایه تقریب برای مشتقات و درون‌یابی خطی برای جملات تأخیر است.
بر اساس این روش، مسأله کنترل بهینه تأخیری به یک دستگاه معادلات جبری تبدیل می‌شود که می‌توان آنها را به صورت عددی حل نمود. کارايي روش پيشنهادي با حل چندین مثال عددی مورد ارزیابی قرار می‌گیرد و نشان داده می‌شود. بنابراین نوآوری مقاله حاضر به شرح زیر است:
- به‌دست‌آوردن معادلات اولر- لاگرانژ در سیستمهای با دینامیک تأخیری جهت برقراری شرایط لازم بهینگی
- استفاده از تقریب مشتق برای تبدیل مشتقات به معادلات خطی
- گسسته‌سازی جملات تأخیر با درون‌یابی خطی
- تبدیل مسأله به یک مسأله برنامه‌ریزی خطی
2- بیان مسأله
در یک مسأله کنترل بهینه تأخیری 3(DOCP)، معیار عملکرد و یا دینامیک سیستم و یا هر دو، حداقل شامل یک جمله مشتق و همچنین تأخیر زمانی هستند. در این مقاله، منظور از حل مسأله کنترل بهینه تأخیری، به‌دست‌آوردن کنترل بهینه برای سیستم‌هایی است که دینامیک سیستم شامل مشتقات مرتبه صحیح بوده و همچنین متغیر حالت دارای تأخیر زمانی ثابت است که به صورت زیر تعریف می‌شود. هدف در مسأله کنترل بهینه تأخیری، پیداکردن کنترل بهینه 
برای یک سیستم دینامیکی تأخیری به قسمی است که معیار عملکرد زیر حداقل گردد
(1)
در این مسأله، دینامیک سیستم و حالت اولیه به صورت رابطه زیر در نظر گرفته می‌شود
(2)
که در این رابطه 
و 
به ترتیب بردار حالت و کنترل، 
و 
ماتریس‌های نیمه‌معین مثبت4، 
ماتریس معین مثبت5، 
معرف مشتق مرتبه اول 
، 
و 
ماتریس‌های پیوسته معلوم با ابعاد مناسب، 
تأخیر زمانی ثابت در حالت سیستم، 
ماتریس ثابت و 
تابع حالت اولیه پیوسته و معلوم هستند. همچنین ساختار سیستم بدون اغتشاش و عدم قطعیت در نظر گرفته شده است.
با توجه به اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی، شرایط بهینه لازم را می‌توان به صورت زیر نوشت [6]
(3)
که در آن داریم
(4)
(5)
که 
و قانون کنترل بهینه به‌صورت زیر است
(6)
همان طور که نشان داده شد، (3) يك مسأله مقدار مرزی شامل جملات تأخیر و تقدم بوده که راه حل دقیق این مسأله بسیار دشوار است. برای حل این مسأله، یک طرح عددی کارآمد جهت پیداکردن متغیرهای حالت و کمک حالت در بخش بعدی ارائه خواهد شد.
3- روش پیشنهادی
اساس این روش بر پایه یک تقریب برای مشتقات و درون‌یابی خطی برای جملات تأخیر است. بر اساس این روش، مسأله کنترل بهینه تأخیری به یک دستگاه معادلات جبری تبدیل می‌شود که می‌توان آن را به‌صورت عددی حل نمود.
3-1 تقریب مشتقات
در بخش دوم از این تحقیق، جهت حل عددی معادلات حاصل، ابتدا بازه زمانی 
را به 
زیربازه که اندازه هر کدام با 
نمایش داده می‌شود تقسیم می‌نماییم. زمان در گره 
ام با 
مشخص گردیده و گره‌ها به‌صورت 
نمادگذاری می‌شوند. از تقریب پسرو و پیشرو مشتق مرتبه اول بسط تیلور برای حل عددی مسأله استفاده می‌نماییم
(7)
با استفاده از تقریب (7) برای مشتقات، معادله سیستم (3) را برای حل عددی آن بازنویسی می‌کنیم
(8)
تعمیم روش
در این بخش از تقریب اختلاف مرکزی مشتق مرتبه اول بسط تیلور برای حل عددی مسأله استفاده نموده‌ایم. می‌توان از تقریب گرانوالد- لتنیکف برای مشتقات کسری برای حل مسائل کنترل بهینه سیستمهای کسری تأخیری استفاده کرد که به اختصار به آن اشاره می‌گردد.
تبصره 1
فرض کنید 
تابع وابسته به زمان، 
مرتبه کسری مشتق و 
باشد. آنگاه مشتق کسری چپ و راست ریمان- لیوویل به‌صورت زیر تعریف می‌شوند [17]:
مشتق کسری ریمان- لیوویل چپ6
(9)
مشتق کسری ریمان- لیوویل راست7
(10)
همچنین اگر 
یک عدد صحیح باشد می‌توان نوشت [17]
(11)
(12)
به منظور حل عددی معادلات حاصل برای سیستمهای تأخیری مرتبه کسری، ابتدا بازه زمانی را به زیربازه که اندازه هر کدام
با نمایش داده می‌شود، تقسیم و زمان را در گره ام با مشخص می‌کنیم. گره‌ها به صورت نمادگذاری می‌شوند. تقریب گرانوالد- لتنیکف برای مشتقات به صورت زیر است [17] و [18]
(13)
(14)
که 
و 
تقریب عددی برای و 
در گره 
ام است و داریم
(15)
با استفاده از تقریب (13) و (14) برای مشتقات، معادله سیستم (3) را برای حل عددی آن بازنویسی می‌کنیم
(16)
3-2 تقریب تأخیر
از درون‌یابی خطی [19] برای تقریب جملات تأخیر استفاده خواهد شد. تابع درون‌ياب خطي 
بين دو نقطه 
و 
برای تابع 
در بازه 
به صورت زير است
(17)
در (16) می‌خواهیم که برای تابع 
و 
در بازه 
و 
و 
برای جملات شامل تأخیر، درون‌یابی خطی انجام دهیم. تابع درون‌یاب به شکل زیر خواهد بود
(18)
(19)
برای حل عددی (8) و تبدیل آن به معادلات خطی به کمک (18) و (19)، جمله شامل تأخیر و تقدم را به کمک درون‌یابی خطی تقریب می‌زنیم
(20)
(21)
(22)
(23)
از (22) و (23) و دستگاه (8)، فرم خطی یک دستگاه برای حل جملات شامل دستگاه (8) با معادله به‌دست می‌آید که به‌راحتی قابل حل است
(24)
بر اساس روش پیشنهادی، مسأله کنترل بهینه تأخیری (8) به کمک تقریب درون‌یابی خطی برای جملات شامل تأخیر به یک دستگاه جبری خطی (24) تبدیل شده که به راحتی قابل حل خواهد بود.
ویژگیهاي مهم روش پيشنهادي را می‌توان به صورت زير در نظر گرفت:
الف) استفاده از روش غیرمستقیم اصل پونتریاگین در حل مسأله کنترل بهینه
ب) استفاده از درون‌یابی خطی و تقریب مشتق برای جملات تأخیر و مشتقات
جدول 1: زمان پردازش سپری‌شده (در ثانیه) و تابع هدف 
برای مقادیر مختلف 
و 
در مثال 1.
CPU time | J | h | N |
00079/0 | 3055/1 | 1/0 | 20 |
00141/0 | 1703/1 | 05/0 | 40 |
01253/0 | 0727/1 | 01/0 | 200 |
05477/0 | 0611/1 | 005/0 | 400 |
23118/2 | 0519/1 | 001/0 | 2000 |
34670/12 | 0508/1 | 0005/0 | 4000 |
جدول 3: مقایسه مقادیر تابع هدف برای مقادیر مختلف 
و در مثال 1.
h | 0005/0 | 001/0 | 005/0 | 01/0 | 05/0 | 1/0 |
روش پیشنهادی | 0508/1 | 0519/1 | 0611/1 | 0727/1 | 1703/1 | 3055/1 |
[21] | - | 0513/1 | 0563/1 | 0629/1 | 1217/1 | 2018/1 |
ج) تبدیل مسأله کنترل بهینه تأخیری به مسأله برنامه‌ریزی خطی
د) کارایی، سرعت بالا و در عین حال سادگی روش پیشنهادی
4- مثال‌های عددی
در اين بخش با استفاده از مثال‌های عددی به بررسی قابليت‌های روش پيشنهادی از جمله سادگی و دقت بالای آن در حل مسأله کنترل بهينه خطی تأخیری می‌پردازيم. برای مثال‌های ارائه‌شده، نتايج عددی به‌دست‌آمده از روش پيشنهادی با نتايج موجود در سایر مقالات مقايسه شده است.
مثال 1: مسأله کنترل بهینه زير را با وجود تأخیر در حالت در نظر بگيريد
(25)
به‌طوری که سیستم تأخیری
(26)
با استفاده از روش پیشنهادی برای حل این مسأله در این مثال داریم: 
، 
، 
، 
و 
. بنابراین شرایط بهینه‌سازی‌شده با (8) به‌صورت زیر تبدیل می‌شوند
(27)
با استفاده از روابط بخش 3-2 ابتدا مسأله را گسسته نموده و سپس برای حل عددی (27) و تبدیل آن به معادلات خطی به کمک (20) و (21) جمله شامل تأخیر را به کمک درون‌یابی خطی تقریب می‌زنیم
(28)
جدول 2: مقایسه شاخص عملکرد با روشهای مختلف برای مثال 1 به ازای 
.
مقدار تابع هدف | روش حل مسأله |
6419/1 | [20] |
6497/1 | [8] |
6488/1 | [21] |
6478/1 | [22] |
6497/1 | Present method |
جدول 4: مقایسه زمان پردازش سپری‌شده (در ثانیه) برای مقادیر مختلف و در مثال 1.
h | 0005/0 | 001/0 | 005/0 | 01/0 | 05/0 | 1/0 |
روش پیشنهادی | 3467/12 | 2311/2 | 0547/0 | 0125/0 | 0014/0 | 0007/0 |
[21] | - | 0513/1 | 0719/0 | 0328/0 | 0126/0 | 0100/0 |
(29)
(30)
(31)
از (30) و (31) و دستگاه (27)، فرم خطی یک دستگاه برای حل جملات شامل دستگاه معادلات (27) با استفاده از (24) با تعداد معادله به‌دست می‌آید که به‌راحتی قابل حل است
(32)
روش پیشنهادی توصیف‌شده در بخش 3 به ازای و 
برای حل مسأله کنترل بهینه تأخیری (27) استفاده شده است. برای در جدول 1، شاخص عملکرد 
به دست آمده و مقادیر ارائه شده است. در جدول 2، شاخص عملکرد برای با نتایج به‌دست‌آمده در سایر مقالات مقایسه گردیده است. از جدول‌های 2 و 3 می‌توان نتیجه گرفت که نتایج حاصل از روش پیشنهادی این پژوهش با سایر روشهای حل واقعی ارائه‌شده در [8] و [20] تا [22] نزدیکی خوبی دارند. علاوه بر این، توابع حالت و کنترل به ازای مقادیر مختلف 
در شکل‌های 1 و 2 رسم شده است. شکل‌های به‌دست‌آمده از نتایج شبیه‌سازی این روش با سایر مراجع در دقتهای بالا انطباق خوبی با یکدیگر دارند. از جدول 4 می‌توان دریافت که روش پیشنهادی، زمان پردازش کمتری نسبت به [21] دارد. نتایج به‌دست‌آمده نشان می‌دهند که روش برنامه‌ریزی خطی با ترکیب تقریب مشتق و درون‌یابی خطی، کارآمد و در عین حال ساده است؛ حجم و مدت زمان پردازش را کاهش می‌دهد. نتایج ارائه‌شده در جداول 1 و 2 نشان می‌دهند که نتایج عددی، همگرایی خوبی دارند، هنگامی که تعداد بازه‌های درون‌یابی افزایش می‌یابد. علاوه بر این، زمان پردازش سپری‌شده (در ثانیه) و تابع هدف برای مقادیر مختلف در جدول 1 گزارش شده است.
تبصره 2
در [20] از زیرفضاهای 
متشکل از توابعی استفاده شده که به‌صورت تکه‌ای در بازه تأخیر 
ثابت هستند. این طرح، تقریب
(الف)
(ب)
[1] این مقاله در تاریخ 27 آبان ماه 1401 دریافت و در تاریخ 18 اردیبهشت ماه 1402 بازنگری شد.
مهدی یوسفی طبری، دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل، بابل، ايران، (email: mehdi.yt61@gmail.com).
زهرا رحمانی (نویسنده مسئول)، دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل، بابل، ايران، (email: zrahmani@nit.ac.ir).
علی وحیدیان کامیاد، دانشکده علوم ریاضی، دانشگاه فردوسی مشهد، مشهد، ايران، (email: vahidian@um.ac.ir).
سید جلیل ساداتی، دانشکده مهندسی برق، دانشگاه صنعتی نوشیروانی بابل، بابل، ايران، (email: j.sadati@nit.ac.ir).
[2] . Pade
[3] . Delay Optimal Control Problem
[4] . Positive Semi-Definite
[5] . Positive Definite
[6] . Left Riemann-Liouville Fractional Derivative
[7] . Right Riemann-Liouville Fractional Derivative
شکل 1: توابع حالت و کنترل مثال 1 برای و مقادیر مختلف .
(الف)
(ب)
شکل 2: توابع حالت و کنترل مثال 1 برای و مقادیر مختلف .
متوسط شناخته می‌شود. در [8] نویسندگان از توابع والش برای حل مسأله DFOC استفاده کرده‌اند. در [21] برای سیستم‌هایی که دینامیک آنها شامل مشتقات کسری ریمان- لیوویل است، ابتدا معادلات اولر- لاگرانژ در سیستمهای با دینامیک کسری تأخیری جهت برقراری شرایط لازم بهینگی بررسی و تشریح شده و سپس با استفاده از یک الگوریتم بر پایه یک تقریب برای مشتقات کسری و یک تقریب گسسته‌سازی برای جملات تأخیر، مسأله به یک مسأله برنامه‌ریزی خطی تبدیل می‌شود. در [22] با استفاده از موجک‌های چلیشکوف متعامد و ماتریس‌های عملیاتی، روش عددی مستقیم کارآمدی برای به‌دست‌آوردن راه حل تقریبی مسائل کنترل بهینه کسری تأخیر زمانی پیشنهاد شده است.
مثال 2: مسأله کنترل بهینه زير را با وجود تأخیر در حالت در نظر بگيريد
(33)
به‌طوری که سیستم تأخیری
(34)
در این مثال داریم: ، ، 
، و . بنابراین شرایط بهینه‌سازی‌شده با (16) به صورت زیر تبدیل می‌شوند
(35)
با استفاده از روابط بخش 3-2 ابتدا مسأله را گسسته نموده و سپس برای حل عددی (35) و تبدیل آن به معادلات خطی به کمک (20) و (21)، جمله شامل تأخیر را به کمک درون‌یابی خطی تقریب می‌زنیم
(36)
(37)
(الف)
(ب)
شکل 3: توابع حالت و کنترل مثال 2 برای مقادیر مختلف .
جدول 5: مقایسه شاخص عملکرد با روشهای مختلف برای مثال 2.
مقدار تابع هدف | روش حل مسأله |
2677/6 | [27] |
0079/6 | [26] |
9533/4 | [23] |
7407/4 | [25] |
7967/4 | [16] |
7985/4 | [24] |
7967/4 | [22] |
7978/4 | Present method |
(38)
(39)
از (38) و (39) و دستگاه (35)، فرم خطی یک دستگاه برای حل جملات شامل دستگاه (35) با استفاده از (24) با تعداد معادله به‌دست می‌آید که به‌راحتی قابل حل است
(40)
جدول 5 مقایسه‌ای را بین شاخص عملکرد به‌دست‌آمده و مقادیر ارائه‌شده در سایر منابع ارائه می‌دهد. از نتایج ارائه‌شده می‌توان نتیجه گرفت که نتایج حاصل از روش پیشنهادی این پژوهش با سایر روشهای ارائه‌شده در [16] و [22] تا [25] مطابقت خوبی دارند. به‌علاوه به‌راحتی می‌توان دریافت که مقادیر گزارش‌شده در [26] و [27] با نتایج عددی دیگر مطابقت خوبی ندارند.
توابع حالت و کنترل برای تکرارهای مختلف الگوریتم ارائه‌شده در شکل 3 رسم گردیده است. شکل‌های به‌دست‌آمده از نتایج شبیه‌سازی این روش با سایر مراجع ارائه‌شده در دقتهای بالا، انطباق خوبی با یکدیگر دارند. نتایج جدول 6 نشان می‌دهند هنگامی که تعداد بازه‌های درون‌یابی افزایش می‌یابد نتایج عددی، همگرایی خوبی دارند. علاوه بر این، زمان پردازش سپری‌شده (در ثانیه) و تابع هدف در جدول 6 ذکر شده است.
جدول 6: شاخص عملکرد و زمان پردازش سپری‌شده (در ثانیه) برای مقادیر مختلف
در مثال 2.
CPU time | J | h | N |
002160/0 | 86231/5 | 1/0 | 20 |
003354/0 | 30239/5 | 05/0 | 40 |
012735/0 | 89390/4 | 01/0 | 200 |
037865/0 | 84511/4 | 005/0 | 400 |
34243/2 | 80642/4 | 001/0 | 2000 |
8636/14 | 8016/4 | 0005/0 | 4000 |
جدول 7: مقایسه زمان پردازش سپری‌شده (در ثانیه) برای مقادیر مختلف در مثال 2.
h | 004/0 | 04/0 | 4/0 |
Present method | 064450/0 | 001461/0 | 000918/0 |
[16] | 354545/0 | 002062/0 | 000863/0 |
از جدول 7 می‌توان دریافت روش پیشنهادی زمان پردازش کمتری نسبت به [16] دارد. نتایج به‌دست‌آمده نشان می‌دهند که روش برنامه‌ریزی خطی با تقریب مشتق و درون‌یابی خطی، کارآمد و در عین حال ساده است و حجم و مدت زمان پردازش را کاهش می‌دهد.
تبصره 3
در [26] از توابع والش برای حل مثال 2 استفاده شده است. مرجع [27] روشی مبتنی بر برنامه‌ريزي پوياي تکرارشونده و [23] روشی بر اساس چندجمله‌ای لاگرانژ مرتبه کسری (FLP) و روش همبستگی برای حل مسائل کنترل بهینه کسری تأخیری ارائه کرده است. در [25] راه حلی بر اساس ویژگی‌های توابع ترکیبی ارائه شده که شامل توابع بلوک پالس و چندجمله‌ای برنولی است. در [16] از فرمول اختلاف محدود مرتبه دوم و چندجمله‌ای درون‌یابی aHermite برای مشتقات مرتبه اول برای حل مسأله استفاده شده است. در [24] روش حل بر اساس خواص توابع ترکیبی متشکل از توابع بلوک پالس و چندجمله‌ای برنولی ارائه شده است. مرجع [22] با استفاده از موجک‌های چلیشکوف متعامد و ماتریس‌های عملیاتی، یک روش عددی مستقیم برای حل تقریبی پیشنهاد داده است.
مثال 3: پانتوگراف وسیله‌ای است که جریان الکتریکی را از خطوط هوایی قطارهای الکتریکی یا تراموا جمع‌آوری می‌کند. معادله پانتوگراف از [28] نشأت می‌گیرد که سیستم جمع‌آوری برق بالای سر قطارها را
(الف)
(ب)
شکل 4: توابع حالت و کنترل مثال 3 برای مقادیر مختلف .
جدول 8: شاخص عملکرد و زمان پردازش سپری‌شده (در ثانیه) برای مقادیر مختلف
در مثال 3.
CPU time | J | h | N |
0089/0 | 2896/0 | 1/0 | 20 |
0170/0 | 1831/0 | 02/0 | 100 |
4436/0 | 1669/0 | 005/0 | 400 |
30/245 | 1624/0 | 0005/0 | 4000 |
مدل‌سازی و طراحی مجدد کردند تا اطمینان حاصل شود که تماس در سراسر آن برقرار است. در سال‌های اخیر، تحقیقات زیادی روی این معادله متمرکز شده که در اینجا از یک نمونه از این کارها برای حل به روش ارائه‌شده در این مقاله استفاده می‌نماییم. سیستمهای کنترل بهینه با تأخیر پانتوگراف زیر را در نظر بگیرید [29]
(41)
به‌طوری که
(42)
توجه: روش پیشنهادی علاوه بر سیستمهای ارائه‌شده، برای سیستمهای با تأخیر پانتوگراف قابل تعمیم است. در این مثال از این نوع سیستم‌ها استفاده شده است.
در این مثال داریم: ، 
، و ؛ بنابراین شرایط بهینه‌سازی‌شده با (16) به‌صورت زیر تبدیل می‌شود
(43)
با استفاده از روابط بخش 3، ابتدا مسأله را گسسته نموده و سپس برای حل عددی (43) و تبدیل آن به معادلات خطی به کمک (20) و (21)، جمله شامل تأخیر را به کمک درون‌یابی خطی تقریب می‌زنیم
(44)
(45)
(46)
(47)
از (46) و (47) و دستگاه (43)، فرم خطی یک دستگاه برای حل جملات شامل دستگاه (43) با استفاده از (24) با تعداد معادله به‌دست می‌آید که به‌راحتی قابل حل است
(48)
برای يافتن يک قانون کنترل بهينه، الگوريتم پيشنهادی در بخش 3 را به کار می‌گيريم. نتايج شبيه‌سازی الگوريتم پيشنهادی به ازای مقادیر مختلف در جدول 8 گزارش شده است. اين نتايج شامل مقدار عددی شاخص عملكرد و زمان صرف‌شده توسط پردازنده کامپيوتر در تكرارهای مختلف الگوريتم می‌باشد. همان گونه که در جدول 8 نشان داده شده، با افزایش دقت پاسخ و زمان پردازش افزایش یافته است.
در جدول 9، نتايج حاصل از روش پيشنهادی با پاسخ به‌دست‌آمده از روش ارائه‌شده در [29] تا [32] مقايسه شده است. همان گونه که مشاهده می‌گردد مقدار عددی شاخص عملكرد، نزديكی بسيار خوبی با پاسخ [29] تا [32] دارد و با افزایش دقت پاسخ نسبت به مراجع فوق افزایش می‌یابد. با توجه به نتایج ارائه‌شده در جداول 8 و 9، علاوه بر اینکه زمان
جدول 9: مقایسه شاخص عملکرد با روشهای مختلف برای مثال 3.
مقدار تابع هدف | روش حل مسأله |
173918/0 | [30] |
173968/0 | [29] |
1739290 | [31] |
1739460 | [32] |
1624/0 | Present method |
محاسبات روش پيشنهادی کمتر از روش [29] تا [32] بوده، دقت پاسخ نیز افزایش یافته است. علاوه بر این در شکل 4، توابع حالت و کنترل بر حسب زمان به ازای مقادیر مختلف رسم شده است. با بررسی انجام‌شده بین نتایج شبیه‌سازی در شکل 4 با شکل‌های [29] تا [32] مشاهده می‌شود که مطابقت خوبی با هم دارند. بنابراين روش پیشنهادی در مقايسه با نگرش [29] تا [32] از دقت کافی برخوردار بوده و در عين حال از لحاظ محاسباتی بسيار ساده‌تر است و باعث افزایش سرعت رسیدن به پاسخ و کاهش زمان پردازش می‌گردد.
تبصره 4
نویسندگان [30] یک روش عددی مبتنی بر چندجمله‌ای برنشتاین را برای حل سیستمهای کنترل بهینه با تأخیرهای ثابت ارائه داده‌اند. مرجع [29]، منحنی‌های Bezier را برای حل مسأله کنترل بهینه با تأخیر پانتوگراف پیشنهاد داده است. در [31] از یک روش عددی جایگزین برای حل مسأله کنترل بهینه درجه دوم خطی با استفاده از موجک‌های لژاندر استفاده شده است. در [32] روش عددی جایگزینی برای حل مسأله کنترل بهینه درجه دوم خطی با استفاده از موجک چبیشف ارائه شده است.
5- نتیجه‌گیری
در این مقاله، یک مسأله کنترل بهینه با تأخیر زمانی در نظر گرفته شده است. با استفاده از اصل حداکثر پونتریاگین برای مسائل کنترل بهینه با تأخیر زمانی، شرایط بهینه لازم را برای این مسأله به‌دست می‌آوریم و سپس الگوریتمی جدید برای حل عددی آن ارائه می‌دهیم. این الگوریتم بر پایه تقریب مشتق و درون‌یابی خطی برای جملات تأخیر است.
با استفاده از این روش، مسأله مورد بررسی به یک مسأله برنامه‌ریزی خطی تبدیل شده که حل آن بسیار ساده‌تر از مسأله اصلی است. نتایج به‌دست‌آمده نشان می‌دهند که روش برنامه‌ریزی خطی مبتنی بر ترکیب تقریب مشتق و گسسته‌سازی درون‌یابی خطی، علاوه بر اینکه کارآمد، مؤثر و دقیق است در عین حال ساده نیز می‌باشد و برخلاف سایر روشهای ارائه‌شده از پیچیدگی و حجم کمتری از محاسبات برخوردار است؛ به همین دلیل زمان پردازش کمتری دارد. با وجود این در دقتهای بسیار بالا و افزایش تکرار الگوریتم پیشنهادی، زمان محاسبات افزایش می‌یابد و لازم است در حل مسائلی که محدودیت زمان اهمیت دارد با توجه به دقت پاسخ، محدودیت‌هایی برای تکرار الگوریتم در نظر گرفت. در بخش تعمیم مسأله که در بخش 3 گزارش شده است، حل مسائل کنترل بهینه تأخیری با وجود مشتقات کسری که در تعدادی از روش‌ها امکان‌پذیر نیست به سهولت ممکن می‌شود و یکی از ویژگی‌های مهم این روش است. در انتها با ذکر مثال، اثربخشی و کارایی روش پیشنهادی را نشان داده‌ایم.
در کارهای آینده سعی خواهیم کرد روش جدید ارائه‌شده در این مقاله را در مسائل کاربردی پیچیده‌تری استفاده نماییم. همچنین در نظر داریم این روش را برای حل مسائل کنترل بهینه کسری تأخیری که در قسمت تعمیم مسأله در بخش 3 به آن اشاره شد و تأخیر ترکیبی و سیستمهای غیرخطی گسترش دهیم.
مراجع
[1] M. R. Hestenes, Calculus of Variations and Optimal Control Theory, Wiley, 1966.
[2] A. E. Bryson, "Optimal control-1950 to 1985," IEEE Control Systems Magazine, vol. 16, no. 3, pp. 26-33, Jun. 1996.
[3] M. Jamshidi and C. Wang, "A computational algorithm for large-scale nonlinear time-delay systems," IEEE Trans. on Systems, Man, Cybernetics, vol. 14, no. 1, pp. 2-9, Jan.-Feb. 1984.
[4] M. Malek-Zavarei and M. Jamshidi, Time-Delay Systems: Analysis, Optimization and Applications, Elsevier Science Inc., 1987.
[5] H. J. Sussmann and J. C. Willems, "300 years of optimal control: from the brachystochrone to the maximum principle," IEEE Control Systems Magazine, vol. 17, no. 3, pp. 32-44, 1997.
[6] G. Kharatishdi, "The maximum principle in the theory of optimal processes with time lags," Dokl. Akad. Nauk SSSR, vol. 136, no. 1, pp. 39-43, 1961.
[7] D. Eller, J. Aggarwal, and H. Banks, "Optimal control of linear time-delay systems," IEEE Trans. on Automatic Control, vol. 14, no. 6, pp. 678-687, Dec. 1969.
[8] K. Palanisamy and R. G. Prasada, "Optimal control of linear systems with delays in state and control via Walsh functions," IEE Proc. D-Control Theory and Applications, vol. 130, no. 6, pp. 300-312, Nov. 1983.
[9] K. Inoue, H. Akashi, K. Ogino, and Y. Sawaragi, "Sensitivity approaches to optimization of linear systems with time delay," Automatica, vol. 7, no. 6, pp. 671-679, Nov. 1971.
[10] J. Banas and A. Vacroux, "Optimal piecewise constant control of continuous time systems with time-varying delay," Automatica,
vol. 6, no. 6, pp. 809-811, Nov. 1970.
[11] H. R. Marzban and M. Razzaghi, "Optimal control of linear delay systems via hybrid of block-pulse and Legendre polynomials," J. of the Franklin Institute, vol. 341, no. 3, pp. 279-293, May 2004.
[12] L. Y. Lee, "Numerical solution of time‐delayed optimal control problems with terminal inequality constraints," Optimal Control Applications and Methods, vol. 14, no. 3, pp. 203-210, Jul./Sept. 1993.
[13] ها. چهکندی نژاد، م. فرشاد و ر هاونگی، "ارائه یک روش جدید به منظور تخمین برخط تأخیر زمانی در سیستمهای SISO-LTI با تأخیر زمانی متغیر با زمان و نامعلوم در ورودی کنترلی،" نشریه مهندسی برق و مهندسی کامپیوتر ایران، الف- مهندسی برق، سال 18، شماره 1، صص. 44-36، بهار 1399.
[14] H. Marzban and S. Hoseini, "Solution of linear optimal control problems with time delay using a composite Chebyshev finite difference method," Optimal Control Applications and Methods,
vol. 34, no. 3, pp. 253-274, May/Jun. 2013.
[15] X. T. Wang, "Numerical solutions of optimal control for time delay systems by hybrid of block-pulse functions and Legendre polynomials," Applied Mathematics, vol. 184, no. 2, pp. 849-856, 15 Jan. 2007.
[16] A. Jajarmi and M. Hajipour, "An efficient finite difference method for the time‐delay optimal control problems with time‐varying delay," Asian J. of Control, vol. 19, no. 2, pp. 554-563, Mar. 2017.
[17] I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Elsevier, 1998.
[18] O. P. Agrawal and D. Baleanu, "A Hamiltonian formulation and a direct numerical scheme for fractional optimal control problems," J. of Vibration and Control, vol. 13, no. 9-10, pp. 1269-1281, 2007.
[19] R. L. Burden, J. D. Faires, and A. Reynolds, Numerical Analysis Prindle, Weber & Schmidt, pp. 28-34, 1985.
[20] H. Banks and J. Burns, "Hereditary control problems: numerical methods based on averaging approximations," SIAM J. on Control Optimization, vol. 16, no. 2, pp. 169-208, 1978.
[21] A. Jajarmi and D. Baleanu, "Suboptimal control of fractional-order dynamic systems with delay argument," J. of Vibration and Control, vol. 24, no. 12, pp. 2430-2446, 2017.
[22] L. Moradi, F. Mohammadi, and D. Baleanu, "A direct numerical solution of time-delay fractional optimal control problems by using Chelyshkov wavelets," J. of Vibration and Control, vol. 25, no. 2, pp. 310-324, 2019.
[23] S. Sabermahani, Y. Ordokhani, and S. A. Yousefi, "Fractional-order Lagrange polynomials: an application for solving delay fractional optimal control problems," Trans. of the Institute of Measurement and Control, vol. 41, no. 11, pp. 2997-3009, 2019.
[24] H. R. Marzban and F. Malakoutikhah, "Solution of delay fractional optimal control problems using a hybrid of block-pulse functions
and orthonormal Taylor polynomials," J. of the Franklin Institute, vol. 356, no. 15, pp. 8182-8215, Oct. 2019.
[25] N. Haddadi, Y. Ordokhani, and M. Razzaghi, "Optimal control of delay systems by using a hybrid functions approximation," J. of Optimization Theory and Applications, vol. 153, no. 2, pp. 338-356, 12 Oct. 2012.
[26] K. Palanisamy, K. Balachandran, and R. Ramasamy, "Optimal control of linear time-varying delay systems via single-term Walsh series," IEE Proc. D-Control Theory and Applications, vol. 135, no. 4, pp. 332-332, Jul. 1988.
[27] S. Dadebo and R. Luus, "Optimal control of time‐delay systems by dynamic programming," Optimal Control Applications and Methods, vol. 13, no. 1, pp. 29-41, Jan./Mar. 1992.
[28] J. R. Ockendon and A. B. Tayler, "The dynamics of a current collection system for an electric locomotive," Proc. of the Royal Society A: Mathematical Physical Sciences, vol. 322, no. 1551, pp. 447-468, 4 May 1971. 1971.
[29] F. Ghomanjani, M. H. Farahi, and A. V. Kamyad, "Numerical solution of some linear optimal control systems with pantograph delays," IMA J. of Mathematical Control and Information, vol. 32, no. 2, pp. 225-243, Jun. 2015.
[30] N. Ghaderi and M. H. Farahi, "The numerical solution of some optimal control systems with constant and pantograph delays via bernstein polynomials," Iranian J. of Mathematical Sciences Informatics, vol. 15, no. 2, pp. 163-181, 2020.
[31] M. Fatehi, M. Vali, and M. Samavat, "State analysis and optimal control of linear time-invariant scale systems using the legendre wavelets," Canadian J. on Automation, Control & Intelligent Systems, vol. 3, no. 1, pp. 1-7, May 2012.
[32] M. Fatehi, M. Samavat, M. Vali, and F. Khaleghi, "State analysis and optimal control of linear time-invariant scaled systems using the Chebyshev wavelets," Contemporary Engineering Sciences, vol. 5, no. 2, pp. 91-105, 2012.
مهدی یوسفی طبری تحصيلات خود را در مقاطع كارشناسي و كارشناسي ارشد مهندسی برق به ترتيب در سالهاي 1383 و 1391 در دانشگاه آزاد اسلامی به پايان رسانده است. وی در حال حاضر دانشجوی دوره دكتري مهندسي برق در دانشکده مهندسی برق و کامپیوتر دانشگاه صنعتی نوشیروانی بابل است. زمينه‌هاي تحقيقاتي مورد علاقه ايشان عبارتند از: کنترل بهینه، سیستمهای مرتبه کسری و سیستمهای تاخیری.
زهرا رحمانی در سال 1377 مدرک كارشناسي مهندسي برق خود را از دانشگاه صتعتي شريف تهران و بهترتیب در سال 1379 مدرك كارشناسي ارشد و در سال 1386 دکتری مهندسي برق خود را از دانشگاه علم و صنعت ایران دريافت نمود. از سال 1387 بهعنوان عضو هيأت علمي گروه مهندسی برق و کامپیوتر دانشگاه صنعتی نوشیروانی بابل مشغول به فعاليت گرديد و اينك نيز دانشیار اين دانشكده ميباشد. زمينه‌هاي تحقیقاتی مورد علاقه ایشان متنوع بوده و شامل موضوعاتي مانند روشهای هوشمند، سیستمهای غیرخطی، پیچیده و ترکیبی است.
علی وحیدیان کامیاد تحصيلات خود را در مقاطع كارشناسي و كارشناسي ارشد ریاضیات کاربردی به ترتيب در سالهاي 1348، 1352 از دانشگاه فردوسی مشهد و موسسه ریاضیات ایران به پايان رسانیده است. وی مدرک دکتری رشته ریاضیات کاربردی- کنترل خود را از دانشگاه لیدز انگلستان در سال 1367 اخذ نمود. ایشان استاد تمام دانشکده علوم ریاضی دانشگاه فردوسی می‌‌باشند و زمينه‌هاي تحقيقاتي مورد علاقه ايشان عبارتند از: ریاضیات کاربردی، دینامیک کسری، کنترل بهینه، کاربردهای تئوری فازی، تئوری اندازهگیری و سیستمهای کنترلی است.
سید جلیل ساداتی رستمی تحصيلات خود را در مقاطع كارشناسي و كارشناسي ارشد و دکتری مهندسی برق گرایش کنترل بهترتيب در سالهاي 1381، 1384 و 1391 از دانشگاه صنعتی خواجه نصیرالدین طوسی و فردوسی مشهد و دانشگاه مازندران به پايان رسانده است و هماكنون استادیار دانشكده مهندسي برق و كامپيوتر دانشگاه صنعتي نوشیروانی بابل ميباشد. زمينه‌هاي تحقيقاتي مورد علاقه ايشان عبارتند از: کنترل سیستمهای مرتبه کسری، کنترل سیستمهای دارای تاخیر، کنترل غیرخطی، کنترل مدل پیشبین، کنترل یادگیر تکرارشونده و یادگیری تقویتی.
