مقدمه
شبکه های عصبی چند لایه پیش خور1 به طور وسیعی د ر زمینه های متنوعی از قبیل طبقه بندی الگوها، پردازش تصاویر، تقریب توابع و ...
مورد استفاده قرار گرفته است.
الگوریتم یادگیری پس انتشار خطا2، یکی از رایج ترین الگوریتم ها جهت آموزش شبکه های عصبی چند لایه پیش خور می باشد.
این الگوریتم، تقریبی از الگوریتم بیشترین تنزل3 می باشد و در چارچوب یادگیری عملکردی 4 قرار می گیرد.
عمومیت یافتن الگوریتمBP ، بخاطر سادگی و کاربردهای موفقیت آمیزش در حل مسائل فنی- مهندسی می باشد.
علیرغم، موفقیت های کلی الگوریتم BP در یادگیری شبکه های عصبی چند لایه پیش خور هنوز، چندین مشکل اصلی وجود دارد:
- الگوریتم پس انتشار خطا، ممکن است به نقاط مینیمم محلی در فضای پارامتر، همگرا شود.
بنابراین زمانی که الگوریتم BP همگرا می شود، نمی توان مطمئن شد که به یک جواب بهینه رسیده باشیم.
- سرعت همگرایی الگوریتم BP، خیلی آهسته است.
از این گذشته، همگرایی الگوریتم BP، به انتخاب مقادیر اولیه وزنهای شبکه، بردارهای بایاس و پارامترها موجود در الگوریتم، مانند نرخ یادگیری، وابسته است.
در این گزارش، با هدف بهبود الگوریتم BP، تکنیک های مختلفی ارائه شده است.
نتایج شبیه سازیهای انجام شده نیز نشان می دهد، الگوریتم های پیشنهادی نسبت به الگوریتم استاندارد BP، از سرعت همگرایی بالاتری برخوردار هستند.
خلاصه ای از الگوریتم BP
از قانون یادگیری پس انتشار خطا (BP)، برای آموزش شبکه های عصبی چند لایه پیش خور که عموماً شبکه های چند لایه پرسپترون 5 (MLP) هم نامیده می شود، استفاده می شود، استفاده می کنند.
به عبارتی توپولوژی شبکه های MLP، با قانون یادگیری پس انتشار خطا تکمیل می شود.
این قانون تقریبی از الگوریتم بیشترین نزول (S.D) است و در چارچوب یادگیری عملکردی قرار می گیرد.
بطور خلاصه، فرایند پس انتشار خطا از دو مسیر اصلی تشکیل می شود.
مسیر رفت6 و مسیر برگشت 7 .
در مسیر رفت، یک الگوی آموزشی به شبکه اعمال می شود و تأثیرات آن از طریق لایه های میانی به لایه خروجی انتشار می یابد تا اینکه
_________________________________
1.
Multi-Layer Feedforward Neural Networks
2.
Back-Propagation Algorithm
3.
Steepest Descent (S.D)
4.
Performance Learning
5.
Multi Layer Perceptron
6.
Forward Path
7.
Backward Path
نهایتاً خروجی واقعی شبکه MLP، به دست می آید.
در این مسیر، پارامترهای شبکه (ماتریس های وزن و بردارهای بایاس)، ثابت و بدون تغییر در نظر گرفته می شوند.
در مسیر برگشت، برعکس مسیر رفت، پارامترهای شبکه MLP تغییر و تنظیم می گردند.
این تنظیمات بر اساس قانون یادگیری اصلاح خطا1 انجام می گیرد.
سیگنال خطا، رد لایه خروجی شبکه تشکیل می گردد.
بردار خطا برابر با اختلاف بین پاسخ مطلوب و پاسخ واقعی شبکه می باشد.
مقدار خطا، پس از محاسبه، در مسیر برگشت از لایه خروجی و از طریق لایه های شبکه به سمت پاسخ مطلوب حرکت کند.
در شبکه های MLP، هر نرون دارای یک تابع تحریک غیر خطی است که از ویژگی مشتق پذیری برخوردار است.
در این حالت، ارتباط بین پارامترهای شبکه و سیگنال خطا، کاملاً پیچیده و و غیر خطی می باشد، بنابراین مشتقات جزئی نسبت به پارامترهای شبکه به راحتی قابل محاسبه نیستند.
جهت محاسبه مشتقات از قانون زنجیره ای2 معمول در جبر استفاده می شود.
فرمول بندی الگوریتم BP
الگوریتم یادگیری BP، بر اساس الگوریتم تقریبی SD است.
تنظیم پارامترهای شبکه، مطابق با سیگنالهای خطا که بر اساس ارائه هر الگو به شبکه محاسبه می شود، صورت می گیرد.
الگوریتم بیشترین تنزل با معادلات زیر توصیف می شود:
(1)
(2)
به طوری WLji و bLj، پارامترهای نرون j ام در لایه iام است.
α، نرخ یادگیری2 و F، میانگین مربعات خطا می باشد.
(3)
(4)
(5)
به طوریکه SLj(k)، حساسیت رفتار شبکه در لایه L ام است.
_________________________________
1.
Error-Correctting Learning Rule
2.
Chain Rule
3.
Learning Rate
معایب الگوریتم استاندارد پس انتشار خطا1 (SBP)
الگوریتم BP، با فراهم آوردن روشی از نظر محاسباتی کارا، رنسانسی در شبکه های عصبی ایجاد نموده زیرا شبکه های MLP، با قانون یادگیری BP، بیشترین کاربرد را در حل مسائل فنی- مهندسی دارند.
با وجود، موفقیت های کلی این الگوریتم در یادگیری شبکه های عصبی چند لایه پیش خود، هنوز مشکلات اساسی نیز وجود دارد:
- اولاً سرعت همگرایی الگوریتم BP آهسته است.
همانطور که می دانیم، تغییرات ایجاد شده در پارامترهای شبکه (ماتریس های وزن و بردارهای بایاس)، پس از هر مرحله تکرار الگوریتم BP، به اندازه ، است، به طوریکه F، شاخص اجرایی، x پارامترهای شبکه و α، طول قدم یادگیری است.
از این، هر قدر طول قدم یادگیری، α، کوچکتر انتخاب گردد، تغییرات ایجاد شده در پارامترهای شبکه، پس از هر مرحله تکرار الگوریتم BP، کوچکتر خواهد بود، که این خود منجر به هموار گشتن مسیر حرت پارامترها به سمت مقادیر بهینه در فضای پارامترها می گردد.
این مسئله موجب کندتر گشتن الگوریتم BP می گردد.
بر عکس با افزایش طول قدم α، اگرچه نرخ یادگیری و سرعت یادگیری الگوریتم BP افزایش می یابد، لیکن تغییرات فاحشی در پارامترهای شکه از هر تکراربه تکرار بعد ایجاد می گردد، که گاهی اوقات موجب ناپایداری و نوسانی شدن شبکه می شود که به اصطلاح می گویند پارامترهای شبکه واگرا شده اند:
در شکل زیر، منحنی یادگیری شبکه برای جدا سازیالگوها در مسأله XOR، به ازای مقادیر مختلف نرخ یادگیری، نشان داده شده است.
به ازای مقادیر کوچک، α، شبکه کند اما هموار، یاد نمی گیرد الگوهای XOR را از هم جدا نماید، ددر صورتی که به ازای 9/0= α شبکه واگرا می شود.
_________________________________
1.
Standard Back-Propagation Algorithm
شکل (1).
منحنی یادگیری شبکه برای نرخ های یادگیری مختلف در مسأله XOR
- ثانیاً احتمالاً به دام افتادن شبکه در نقاط مینیمم محلی وجود دارد.
در شبکه های MLP، میانگین مجوز خطا، در حالت کلی خیلی پیچیده است و از تعداد زیادی نقطه اکسترمم در فضای پارامترهای شبکه برخوردار می باشد.
بنابراین الگوریتم پس انتشار خطا با شروع از روی یک سری شرایط اولیه پارامترهای شبکه، به نقطه مینیمم سراسری و با شروع از یک مجموعه شرایط اولیه دیگر به تقاط مینیمم محلی در فضای پارامترها همگرا می گردد، بنابراین زمانی که الگوریتم BP همگرا می شود، نمی توان مطمئن شد که به یک جواب بهینه رسیده باشیم.
- ثالثاً: همگرایی الگوریتم BP، به یقین مقادیر اولیه پارامترهای شبکه عصبی MLP وابسته است، بطوری که یک انتخاب خوب می تواند کمک بزرگی در همگرایی سریعتر الگوریتم BP فراهم آورد.
برعکس انتخاب اولیه نادرست پارامترهای شبکه MLP، منجر به گیر افتادن شبکه در نقاط مینیمم محلی در فضای برداری پارامترهای شبکه می گردد که این خود منجر به این می شود که شبکه خیلی زودتر از معمول به موضعی بیفتد که منحنی یادگیری شبکه برای تعداد بزرگی از دفعات تکرار، تغییر نکند.
به عنوان مثال، فرض می کنیم مقدار اولیه پارامترهای شبکه خیلی بزرگ باشند، در حالی که می دانیم توابع تبدیل نرونها مخصوصاً در لایه های میانی از نوع زیگموئید هستند.
در این حالت برای نرون i ام، اندازه ورودی تابع تبدیل (ni) خیلی بزرگ می باشد و خروجی نرون (ai) به مقدار 1± میل می کند.
لذا مشتق بردار خروجی شبکه، a ، خیلی کوچک می باشد.
فرض کنیم که باید مقدار واقعی ai، 1 باشد یعنی ti = 1، لیکن به خاطر انتخاب بر مقادیر اولیه، ai = -1 گردد.
در این حالت خطای حداکثر را داریم در حالی که چون ai ≈ 0 می باشد تغییرات ناچیزی در پارامترهای متناظر با نرون i ام داریم.
این چیزی است که بیانگر رسیدن زودتر از معمول نرونها به حد اشباع خود می باشند، جایی که پاسخ واقعی با پاسخ شبکه کاملاً فرق دارد و زمان زیادی طول خواهد کشید که نرون از این حالت خارج شود.
از این رو با پیشرفت پروسه یادگیری، پارامترهای منتسب به نرورنهایی که به مرز اشباع نرسیده اند، سریعتر تنظیم می شوند، چرا که سیگنال خطار گرادیانهای محلی از مقدار از اندازه بزرگتری برخوردار می باشند.
این عمل منجر به کاهش در مجموع مربعات خطای لحظه ای می گردد و اگر در این مرحله، نرونهای به حد اشباع رسیده تغییری در وضعیت تحریکشان رخ ندهد، شبکه برای مدتی طولانی از یک شکل هموار منحنی خطا برخوردار خواهدبود.
بهبود الگوریتم استاندارد پس انتشار خطا (SBP)
- الگوریتم BP از نوع دسته ای1 (BBP)
الگوریتم استاندارد BP، بر اساس فرم الگو به الگو است، بدین ترتیب که پارامترهای شبکه پس از ارائه هریک از الگوهای یادگیری که عموماً بطور تصادفی انتخاب می شوند، تنظیم می گردند، اما در الگوریتم BBP، تنظیم پارامترهای شبکه پس از اعمال تمامی ورودی ها صورت می پذیرد.
پردازش دسته ای موجب می شود که گرادیانهای محلی به گرادیان محلی واقعی نزدیکتر باشند و نهایتاً الگوریتم BP به الگوریتم بیشترین نزول نزدیکتر گردد که این خود موجب می شود همگرایی الگوریتم BP افزایش یابد.
در شکل زیر مسئله XOR با متد الگوریتم BP به فرم دسته ای پیاده شده است.
به راحتی می توان دید که الگوریتم BBP از سرعت همگرایی بالاتری به الگوریتم SBP برخوردار است.
شکل (2).
رفتار شبکه با الگوریتم BBP در مسأله XOR ( )
رفتار شبکه با الگوریتم SBP (0)
_________________________________
1.
Batching Back-Propagation Algorithm
- روش ممنتم 1 برای الگوریتم BP (MBP)
همانطور که مشاهده شد، اگر نرخ یادگیری α، کوچک انتخاب شود، متد BP که در واقع همان تقریب الگوریتم SD است، بسیار کند می گردد.
و اگر α، بزرگتر انتخاب شود، شبکه نوسانی خواهد بود.
یک راه ساده و مؤثر که عموماً جهت افزایش و بهبود نرخ یادگیری، استفاده می شود- جایی که خطر ناپایداری و نوسانی شدن شبکه جلوگیری می گردد- افزودن یک جمله ممنتم در الگوریتم تقریبی SD می باشد، یعنی به هر پارامتر از شبکه MLP، یک مقدار اینرسی یا اندازه حرکت اضافه می شود تا اینکه پارامتر مورد نظر در مسیری تمایل به تغییر داشته باشد که کاهش تابع انرژی احساس شود.
الگوریتم یادگیری MBP با معادلات زیر قابل توصیف است:
(6)
(7)
جایی که (1و0) ، ترم ممنتم را نشان می دهد و عموماً با نرخ یادیگیری به صورت زیر رابطه دارد:
(8)
معادلات فوق، ترم های اصلاحی پارامترهای شبکه را از فیلتر پایین گذر عبور می دهند واین یعنی تغییرات با فرکانس بالا(نوسانات شدید) فیلترمی شوند.
شکل (3)، مسأله XOR را که به وسیله الگوریتم MBP، پیاده شده است، نشان می دهد.
شبیه سازی نشان می دهد که الگوریتم MBP نسبت به الگوریتم استاندارد BP، از سرعت همگرایی بالاتری برخوردار است.
شکل(3): رفتار شبکه با الگوریتم MBP درمسأله XOR ()
رفتار شبکه با الگوریتم SBP(.)
_______________________________
1.
Momentum
- نرخ یادگیری متغیر1 (VLR)
درالگوریتم BP استاندارد، نرخ یادگیری درطول فرآیند یادگیری ثابت نگه داشته می شود.
عملکرد الگوریتم به انتخاب مناسب نرخ یادگیری خیلی حساس می باشد.
اگر نرخ یادگیری خیلی بزرگ انتخاب شود ممکن است الگوریتم نوسان کرده وناپایدار شود.
اگر نرخ یادگیری خیلی کوچک باشد زمان زیادی طول خواهد کشید تا الگوریتم همگرا شود.
انتخاب نرخ یادگیری اپتیمم قبل از یادگیری، عملی نبوده ودرحقیقت نرخ یادگیری اپتیمم به هنگام پروسه آموزش، همچنان که الگوریتم برروی سطح خطا حرکت می کنددائماً تغییر می کند.
اگر اجازه دهیم نرخ یادگیری بهنگام پروسه آموزش تغییر کند عملکرد الگوریتم BP استاندارد را می توان بهبود بخشید.
نرخ یادگیری تطبیقی سعی می کند که نرخ یادگیری را تا آنجایی که ممکن است و سیستم ناپایدار نشده است، افزایش دهد.
اگر اجازه دهیم نرخ یادگیری بهنگام پروسه آموزش تغییر کند عملکرد الگوریتم BP استاندارد را می توان بهبود بخشید.
نرخ یادگیری تطبیقی سعی می کند که نرخ یادگیری را تا آنجایی که ممکن است و سیستم ناپایدار نشده است، افزایش دهد.
نرخ یادگیری تطبیقی نیاز به تغییراتی در الگوریتم BP استاندارد دارد.
مراحل الگوریتم VLR، به طور خلاصه در زیر بیان شده است: _______________________________ 1.
Variable Learning Rate معمولاً مقادیر 1.04، 0.7 و 1.05 به ترتیب برای ضرائب 1 inc- perf- max، 2 dec- Lr و 3 inc- Lr در نظر گرفته می شود.
در این روش نرخ یادگیری به اندازه ای افزایش می یابد که موجب افزایش زیاد در خطاهای آموزش نگردد.
بنابرین، یک نرخ یادگیری نزدیک به بهینه بدست می آید.
الگوریتم VLR بر روی مسأله XOR پیاده شده است.
شبیه سازی نسان می دهد، که سرعت یادگیری این الگوریتم را از الگوریتم SBP، بیشتر است.
نتایج شبیه سازی در شکل (4) نشان داده شده است.
شکل (4).
– رفتار شبکه با الگوریتم VLR برای مسأله XOR ( ـــ ) رفتار شبکه با الگوریتم SBP (0 ـــ) - تغییرات نرخ یادگیری (α) در کل فرآیند یادگیری برای مسأله XOR در این، خلاصه ای از نتایج مقالات مورد بررسی، جهت بهبود الگوریتم یادگیری پس انتشار خطا، ذکر می گردد: 1- الگوریتم BP، با وجود آنکه کاربردهای فراوانی دارد، متأسفانه نرخ همگرایی آن بسیار پائین است.
اگر از نرخ یادگیری کوچک استفاده شود، این مسأله می تواند سبب نرخ پائین همگرایی.
به منظور جلوگیری از این پدیده نامطلوب، یک راه حل، استفاده از نرخ یادگیری می باشد.
این راه حل بویژه در مواقعی که در نقطه ای از سطح قرار داریم که شیب تندی دارد، نامطلوب بوده و می تواند باعث واگرایی شبکه شود.
_________________________________ 1.
Maximum Performance Increase 2.
Learning Rate decrease 3.
Learning Rate increase موقعیکه شکل سطح خطا از سطح خطای درجه دوم خیلی فاصله داشته باشد، در اینصورت سطح خطا، شامل مناطق با شیب تند زیادی خواهد بود.
در این صورت الگوریتم BP با نرخ یادگیری ثابت دارای راندمان پائینی می باشد، دلیل این مسأله این است که به منظور جلوگیری از نوسان در مناطقی که سطح خطا دارای شیب زیادی است بایستی نرخ یادگیری کوچک انتخا شود، در ننیجه بردار وزن، موقعی که در مناطق مسطح قرار داریم به دلیل کوچک بودن گرادیان، خیلی کند حرکت خواهد کرد.
بنابراین نیاز به یک الگوریتم یادگیری کارآمد می باشد تا بتواند بطور پویا نرخ یادگیری را تغییر دهد.
- الگوریتم پس انتشار خطای تطبیقی1 (ABP) در الگوریتم پس انتشار خطای تطبیقی، نرخ یادگیری، بطور، اتوماتیک و بر اساس خطای آموزش تنظیم می گردد.
]1[ ایده اصلی از الگوریتم پس انتشار خطای تطبیقی آن است که: - اگر علامتهای گرادیان ، در طول دو تکرار متوالی مخالف هم باشد، دلالت بر این امر دارد که خطای جدید نسبت به خطای قبلی، افزایش یافته است و تکرار جاری با ارزش نیست.
لذا نرخ یادگیری را بایستی کاهش داد.
- اگر علامتهای امتدادگرادیان در طول دو تکرار متوالی، یکسان باشند، دلیل بر این امر است که نرخ تنزل آهسته است و نرخ یادگیری را بایستی افزایش داد.
ایده فوق را می توان با معادلات زیر، توصیف کرد: (8) به طوری که را به صورت زیر، تعریف می شود: (9) مشخص است که زمانی که λ>0 ، نرخ یادگیری افزایش می یابد.
الگوریتم BP تطبیقی با معادلات زیر بیان می شود: (10) شکل (5)، رفتار شبکه را با قانون یادگیری فوق، جهت جداسازی الگئهای XOR، نشان می دهد.
Adaptive Back-Propagation Algorithm شکل (5).
منحنی یادگیری الگوریتم BP تطبیقی برای XOR - الگوریتم پس انتشار خطا با نرخ یادگیری و ضریب ممنتم تطبیقی1 (BPALM) در این الگوریتم نرخ یادگیری و ضریب ممنتم در هر سیکل به طور تطبیقی تنظیم می شوند تا به همگرایی الگوریتم BP استاندارد بهبود بخشیده شود]2[.
در ابتدا، فاکتور نسبی er (k)، به صورت زیر تعریف می شود: (11) اصلاح و تنظیمات نرخ یادگیری و ضریب ممنتم به صورت زیر انجام می گیرد: Case 1: for er (k) α (k+1) = α (k) [1+ue –er(k)] ; uε (0,1) (12a) δ (k+1) = δ (k) [1+ve –er(k)] ; vε (0,1) (12b) Case 2: for er (k) ≥ 0 α (k+1) = α (k) [1-ue –er(k)] ; uε (0,1) (13a) δ (k+1) = δ (k) [1-ve –er(k)] ; vε (0,1) (13b) _________________________________ 1.
Back- Propagation with Adaptive Learning rate and Momentum term شکل زیر، عملکرد شبکه را با قانون یادگیری BPALM، جهت جداسازی الگوهای XOR، نشان می دهد.
شکل(6).
– منحنی یادگیری الگوریتم BPALM در مسأله XOR - تغییرات نرخ یادگیری - تغییرات ضریب ممنتم - تغییرات علامت1 بر اساس این الگوریتم، اگر علامت مشتق شاخص اجرایی نسبت به پارامترهای شبکه دردوتکرار متوالی تغییر نکند، نرخ یادگیری افزایش می یابد و در غیر اینصورت، نرخ یادگیری کاهش می یابد.
[6] و این یعنی: (14 a) αji(k) = αji(k-1).u if (14 b) αji(k) = αji(k-1).u if آبگاه، اصلاح وزن بر اساس معادله زیر، انجام می گیرد: (15) معموماً مقادیر 1.1-1.3 برای u و 0.7-0.9 برای d، به کار برده می شود.
Sign Changes بر اساس این الگوریتم، اگر علامت مشتق شاخص اجرایی، نسبت به پارامترهای شبکه در دوتکرارمتوالی تغییر نکند و نرخ یادگیری از مقدار حداکثری، کمتر باشد، آنگاه افزودن یک مقدار ثابت به نرخ یادگیری، نرخ یادگیری افزتیش می یابد.
اگر علامت مشتق شاخص اجرایی نسبت به پارامترهای شبکه دردوتکرارمتوالی تغییر کند، نرخ یادگیری با ضرب شدن در مقدار کوچکتر، کاهش می یلبد و در غیر اینصورت تغییری در نرخ یادگیری نخواهیم داشت.
[6] عبارت فوق را می توان با روابط زیر توصیف کرد: (16 a) αji(k) = αji(k-1)+u if (16 b) αji(k) = αji(k-1).d if (16 c) αji(k) = αji(k-1) if شکل زیر، منحنی یادگیری الگوریتم فوق را برای جداسازی الگوهای XOR، نشان می دهد.
شکل (7).
منحنی یادگیری الگوریتم Delta Bar Delta Rule در مسأله XOR - الگوریتم یادگیری Super SAB این الگوریتم، مانند الگوریتم قبل می باشد با این تفاوت که اگر علامت تغییرات مشتق شاخص اجرایی نسبت به پارامترهای شبکه دردوتکرار متوالی تغییر نکند و نرخ یادگیری از مقدار حداکثری کمتر باشد، نرخ یادگیری با ضرب شدن در یک پارامتر ثابتی، افزایش می یابد [6].
(17 a) αji(k) = αji(k-1)+u if (17 b) αji(k) = αji(k-1).d if (17 c) αji(k) = αji(k-1) if مقادیر پیشنهادی بر u، 1.05 و d، 0.5 می باشد.
α max، عددی بین 1,0 است.
از الگوریتم مطرح شده، برای جداسازی الگوها در مسأله XOR، استفاده شده است.
در شکل (8)، منحنی یادگیری این الگوریتم را نشان می دهد.
شکل (8).
منحنی یادگیری الگوریتم Super SAB برای مسأله XOR 2- الگوریتم پس انتشار خطا با سه ترم در این الگوریتم، ترم جدیدی به نام ضریب تناسبی1 (PF)، به الگوریتم استاندارد پس انتشار خطا اضافه شده است.
[3] افزایش سرعت یادگیری، همراه با حفظ سادگی اگوریتم BP، هدف اصی این مقاله است.
الگوریتم مطرح شده را می توان با معاله زیر توصیف کرد: (18) بطوری که C، ضریب تناسبی (PF) می باشد.
e(W(k))، اختلاف بین خروجی مطلوب و خروجی واقعی در تکرار k ام است.
الگوریتم فوق، دراای سه ترم است: 1- ترمی که متناسب با مشتق F(W(k)) است.
2- ترمی که متناسب با مقادیر قبلی تغییرات وزن هاست.
3- ترمی که متناسب با e(W(k)) است.
بنابراین، مشاهده می شود که این سه ترم، همانند ضرایب کنترل کننده PID، که به طور متداول در کاربردهای کنترلی استفاده می شوند، عمل می کنند.
آنالیز همگرایی در این قسمت، الگوریتم BP با سه ترم، آنالیزو تحلیل می شود و نشان داده می شود که نقاط مینیمم محلی تابع حداقل مربعات خطا، تنها نقاط پایدار مجانبیمحلی الگوریتم هستند.
معادله (18) را می توان به صورت زیر، بازنویسی کرد: (19) فرض کنید x2(k) = W(k) -W(k-1) , x1(k) = W(k)، آنگاه یک تحقق متغیر حالت معادله (19)، چنین است: (20a) (20b) لم (1): m = (m1,m2)، یک نقطه تعادل از سیستم یا معادلات (20b) , (20a) است، اگر _________________________________ 1.
Proproortional Factor اثبات: اگر m = (m1,m2)، یک نقطه تعادل است، آنگاه x1(k) = m1 , x2(k) = m2 و x1(k+1) - x1(k) = 0 x2(k+1) - x2(k) = 0 با جایگذاری در معادلات (20b) , (20a)، خواهیم داشت: (21a) (21b) با کم کردن معادله (21b) از معادله (21a)، نتیجه می شود: x2(k) = 0 m2 = 0 و با جایگذاری x2(k) = 0 در معادله (21a) و یا (21b)، عبارت زیر به دست می آید: (22) تذکر(1): از آنجاییکه m =(m1,m2)،س یک نقطه تعادل معادلات (20b), (20a) است و e(x1(k)) =0 ، آنگاه F(x1(k)) =0.
ویژگی پایداری محلی، حول نقطه تعادل (m1,m2) می تواند با استفاده از آنالیزهای سیگنال کوچک، امتحان شود.
فرض کنید سیگنال های آشفته به صورت زیر تعریف شوند: آنگاه معادلات حالت به فرم زیر تبدیل می شوند: (23a) (23b) با خطی سازی حول محور تعادل m، معادلات (23) , (23a) به صورت زیر خواهند بود: (24a) (24b) بطوریکه Q، سایز بردار وزن است.
آنگاه : (25) معادله ماتریسی بالا، می تواند به صورت زیر نوشته شود: (26) واضح است که سیستم گسسته- زمان به فرم (26)، پایدار مجانبی است اگر θ دارای مقادیر ویژه مجزای (iΨ) باشد و هر مقدار ویژه (iΨ) در شرط رو به رو صدق کند: i| لم (2): اگر λ یک مقدار ویژه تابع E باشد، جایی که .
آنگاه مقادیر ویژه متناظر با ماتریسθ، توسط جوابهای معادله درجه دوم زیر بدست می آیند: (27) اثبات: ماتریس θ برای هر D, A 0, ≠ δ، معکوس پذیر است.فرض کنید Ψ، یک مقدار ویژهθ باشد و آن غیر صفر است (زیراθ ناویژه1 است).
فرض کنید که Z =(x,y)، بردار ویژه غیر صفر متناظر با مقدار ویژه Ψ است.
آنگاه: (28) Ψz= θz که منجر می شود به: x - αAx + CDx + δy = Ψx (29) و -αAx + CDx + δy = Ψy (30) با کم کردن معادله (30) از معادله (29)، خواهیم داشت: (31) با جایگذاری معادله (31) در معادله (30)، خواهیم داشت: (32) با انتخاب و جایگذاری در معادله (32)، رابطه زیر بدست می آید: (33) _________________________________ 1.
Non - singular از آنجایی که ، اسکالر و غیر صفر است، اگر بردارx در این معادله صدق کند، x یک بردار ویژه ماتریس E است.
آنگاه: Ex = λx (34) با مقایسه معادلات (33), (34)، معادله زیر بدست خواهد آمد: (35) با مرتب کردن معادله (35)، خواهیم داشت: ) + δ =0 (36) - Ψ (1+δ-αλc 2 Ψ الف.
تست پایداری جوری: تست جوری، برای آنالیز معادلات از هرمرتبه، به کار می رود.
معادله (36)، یک معادله درجه دوم است و می توان آن را به صورت f(z)=a2 z2 + a1 z + a0 نوشت.
شرایط لازم و کافی برای هرچند جمله ای ، بطوری که ریشه ای خارج و روی دایراه واحد نداشته باشد، آن است: | a0 | f (1) > 0 (38) (-1)2 f(-1) > 0 (39) با بکار بردن تست جوری در معادله (36)، نتیجه می شود: | δ | (1+δ) > ((1+δ) - αλc) (41) (1+δ) > (-(1+δ) + αλc) (42) از آنجایی که، ضریب ممنتم مثبت است: 0 نامعادله (41)، منجر خواهد شد به: αλc > 0 (44) و نا معادله (42)، به فرم زیر تبدیل خواهد شد: (45) با استفاده از نامعادلات (43) و (45)، نتیجه می شود: 0 مقادیر α و c، بایستی مثبت باشند تا سیستم یاد بگیرد.
بنابراین آن شرط کافی است که همه مقادیر ویژه E، مثبت باشند تا نتیجه بگیریم که m، یک نقطه مینیمم محلی است.
بنابراین همه نقاط مینیمم محلی از F(0)، به طور مجانبی پایدار محلی1 هستند.
همچنین از آنجایی که F(0)، یک نقطه مینیمم محلی دارد، ماتریس هسیان2 A، ماتریسی مثبت معین3 است.
نشان داده شده است که اگر ضرایب الگوریتم BP دارای سه ترم، در شرایط (46), (44), (43) صدق کنند، پایداری سیستم تضمین می شود و سیستم به یک نقطه مینیمم محلی، همگرا خواهد شد.
اگر مقادیر ویژه ماتریس E، نسبتاً بزرگ باشندف ممکن است شرط (46)، نقض شود ولی در اکثر موارد E محدود است، بنابراین اگر α و c، به اندازه کافی کوچک باشند، تمامی نقاط، مینیمم محلی پایدار هستند.
ب.
شرط پایداری برای ماتریس D نامعادله (46)، مقادیر نرخ یادگیری α و ضریب تناسبی c و نیز مقادیر ویژه E را محدود می کنند.
هدف از این بخش آن است که بر روی شرایط کافی پایداری بحث کنیمس به طوری که، رابطه بین پارامترهای α وc را بدست آوریم.
بنابراین باید به بحث روی ماتریس D بپردازیم، به طوری که ماتریس E، مثبت معین باشد: (47) تئوری (1): E، مثبت معین است، اگر و فقط اگر، ماتریس بالا مثلثی، حقیقی و ناویژه H، به گونه ای یافت شود که E = HTH.
دو حالت را در نظر می گیریم: حالت اول: می پذیریم که ماتریس سd، ماتریس منفی نیمه معین4 است، آن گاه ماتریس -D ، مثبت نیمه معین5 خواهد بود.
و مقادیر ویژه STBS، مثبت می باشند.
بنابراین مقادیر ویژه معادله (50)، کمتر از نیستند، یعنی: (51) _________________________________ 1.
Locally Asymptotically Stable 2.
Hessian 3.
Positive Definite Matrix 4.
Negative Semidefinite Matrix 5.
Positive Semidefinite Matrix بنابراین STBS، مثبت معین است، از اینرو: (52) det (A) = det (S-TS-1) = (DET (S))-2 بنابراین: (53) معادله فوق را می توان به صورت زیر، مرتب کرد: (54) حالت دوم: فرض می کنیم که D، ماتریسی مثبت نیمه معین است.
تئوری (3): فرض کنید که ماتریس D، مثبت نیمه معین است و A، مثبت معین می باشدآنگاه ماتریس ، مثبت معین است و .
اثبات: بر اساس اثبات قبلی، (55) برای آنکه E، مثبت معین باشد، بایستی طرف راست معادله (55)، مثبت باشد، از اینرو: (56) از الگوریتم مطرح شده برای جداسازی الگوها در مسأله XOR استفاده شده است.
شبیه سازی، نشان می دهد که سرعت همگرایی این الگوریتم، از الگوریتم استانداردپس انتشار خطا بیشتر است.
شکل (9).
منحنی یادگیری الگوریتم BP دارای سه ترم، در مسأله XOR 3.
از رابطه (5)، به خوبی مشخص است که بردارهای حساسیت رفتار شبکه (δjL(K))، شامل ترم مشتق تابع تحریک نرون می باشد با فرض آنکه تابع تحریک نورونها، در لایه های میانی از نوع تابع زیگموئیدی لگاریتمی می باشند این ترم، ajL(K)(1- ajL(K) است.
بردارهای حساسیت از لایه آخر به لایه های میانی و سپس به لایه ورودی، برگشت داده می شوند تا پارامترهای شبکه را تصحیح کنند.
زمانی که خروجی ajL(K) به مقادیر 1 و یا 0، میل می کند.
مشتق تابع تحریک که دارای فاکتور ajL(K)(1- ajL(K) می باشد، سبب می شود که بردار حساسیت δjL(K) ، بسیار کوچک شود، بنابراین تغییرات ناچیز در پارامتر متناظر با نرون j ام خواهیم داشت.
لذا، فرایند یادگیری و اصلاح پارامترها در الگوریتم BP، بسیار آرام خواهد بود و یا حتی متوقف خواهد شد، حتی اگر مقادیر پارامترهای شبکه، از مقادیر بهینه اشان، فاصله ای زیاد داشته باشند.
بنابراین اصلاح مشتق تابع تحریک f، ضروری است تا بدین ترتیب عملکرد الگوریتم BP را بهبود بخشیم.
در یکی از مقالات مورد بررسی [4]، ایده فوق در نظر گرفته شده است و جهت افزایش سیگنال حساسیت بر روی مشتقات جزئی یا شیب تابع تحریک، بهبودهایی صورت گرفته است.
بردار حساسیت در معادله (5)، به صورت زیر تغییر کرده است: (57) (58) بطوریکه، s یک عدد حقیقی بزرگتر از یک است.
(S ≥ 1) زمانیکه S =1 است، الگوریتم همان الگوریتم استاندارد BP می باشد.
برای S >1، بردار حساسیت، زمانی که خروجی به مقدار اشاتباه میل می کند، زیاد خواهد شد.
بنابراین نرخ همگرایی فرایند یادگیری می تواند افزایش یابد.
با این وجود، مقدار S، نباید از حد خاصی بیشتر شود زیرا در غیر اینصورت ترم به یک میل می کند و بردار حساسیت δjL(K)، ناپایدار خواهد شد.
شکل زیر، شبیه سازی که بر اساس این الگوریتم روی مسءله XOR انجام شده است، نشان می دهد.
به خوبی مشخص است که سرعت همگرایی این الگوریتم به مراتب بیشتر زا الگوریتم BP می باشد.
شکل (10).
منحنی یادگیری الگورتم GBP در مسأله XOR به ازای S=2 و S=1 - الگوریتم پس انتشار خطای بهبود پذیر1 (Rprop) این الگوریتم، جهت اصلاح مشکل فوق ارائه شده است.
بر اساس این الگوریتم، تنها از علامت مشتق تابع تحریک، جهت اصلاح پارامترهای شبکه استفاده می شود.
اندازه مشتق تابع تحریک، هیچ اثری بر تنظیم پارامترهای شبکه ندارد [6], [5].
میزان تغییرات در پارامترهای شبکه، توسط فاکتور delt-inc، افزوده می شود، زمانی که علامت مشتق شاخص اجرایی، نسبت به پارامترهای شبکه دردوتکرار متوالی، تغییر نکند.
و زمانی که مشتق خاص اجرایی دردوتکرار متوالی هم علامت نباشند، تغییرات در پارامترهای شبکه توسط فاکتور delt-dec، کاهش می یابد.
الگوریتم Rprop، در زیر خلاصه می گردد.
مقادیر پیشنهادی برای پارامترها عبارتند از: Δmax = 50, Δmin = 0.000001, delt- inc = 1.2, delt-dec = 0.5 .
Resilient Back-Propagation Algorithm نتیجه گیری از قانون یادگیری پس انتشار خطا (BP) برای آموزش شبکه ها عصبی چند لایه پیش خور استفاده می شود.
با وجود کاربردهای فراوان این الگوریتم یادگیری، هنوز مشکلاتی نیز وجود دارد: سرعت همگرایی الگوریتم BP، پائین است و ممکن است شبکه به آسانی به نقاط مینیمم محلی همگرا شود.
از طرفی، انتخاب نرخ یادگیری، تأثر بسزایی در سرعت همگرایی آموزش شبکه عصبی دارند.
در این گزارش، الگوریتم های جدیدی، جهت بهبود الگوریتم BP، ارائه شده است.
برخی از این روش ها بر مبنای نرخ یادگیری تطبیقی می باشند.
بدین صورت که نرخ یادگیری به هنگام پروسه آموزش تغییر می کند تا عملکرد در الگوریتم BP استاندارد بهبود بخشیده شود، نرخ یادگیری تطبیقی سعی می کند که نرخ یادگیری را تا آنجایی که ممکن است و سیستم ناپایدار نشده است، افزایش دهد.
الگوریتم دیگری که جهت بهبود سرعت همگرایی الگوریتم BP، ارائه شده است، الگوریتم BP با سه ترم است.
در این الگوریتم، ترم جدیدی به نام ضریب تناسبی (PE)، علاوه بر دوترم نرخ یادگیری و ضریب ممنتم، به الگوریتم BP، اضافه شده است.
الگوریتم مطرح شده، همانند کنترل کننده PID، عمل می کند.
همچنین در این گزارش، به بررسی و آنالیز پایداری الگوریتم مطرح شده، پرداخته شده است.
آنالیز پایداری به دلیل آن است که، شرایطی را که باید پارامترهای یادگیری در آن صدق کنند، تا الگوریتم پایدار بماند، بدست آوریم.
در آخر نیز، الگوریتم هایی ارائه شده است، یکی از این الگوریتم ها، الگوریتم پس انتشار خطای بهبود پذیر ( Rprop) است.
در این الگوریتم، تنها از علامت مشتق تابع تحریک نسب به پارامترهای شبکه، جهت تنظیم پارامترهای شبکه استفاده می شود و اندازه مشتق تابع تحریک، هیچ تأثیری بر تنظیم پارامترهای شبکه ندارد.
در الگوریتم دیگر، جهت افزایش سیگنال حساسیت، اصلاحاتی بر روی مشتق تابع تحریک نرونها، انجام گرفته است.
شبیه سازی انجام شده نشان می دهد، که سرعت همگرایی الگوریتم های مطرح شده نسبت به الگوریتم استاندارد BP، بیشتر است.
از طرفی الگوریتم های مطرح شده، محاسبات پیچیده ای ندارد و از همان سادگی الگوریتم BP، برخوردار می باشند.
مراجع [1] Wen Jin-Wei, Zhao Jia-Li, Luo Si-Wei and Han Zhen, ''The Improvements of BP Neural Network Learninig Algorithm'', IEEE, 2000.
[2] Chien-Cheng Yu, and Bin-Da Liu, ''A BACKPROPAGATION ALGORITHM WITH ADAPTIVE LEARNING RATE AND MOMENTOM COEFFI CIENT'', IEEE, 2000 [3] Yahya H.Zweiri, James F.Whidborne, Kaspar Althoefer and Lakmal D.Seneviratne, ''A New Three-Term Backpropagation Algorithm with Convergence Analysis'', Proceeding of the IEEE International Conference on Robotics 8Automation, 2002.
[4] S.C.Ng, S.H.Leung and A.Luk, ''A Generalized Back-Propagtion Algorithm for Faster Convergence'', IEEE, 1996.
[5] Riedmiller.M and H.Braun, ''A Direct Adaptive Method for Faster Backpropagation Learning: The RPROP Algorithm'', proceeding of the IEEE International Conference on Neural Networks, 1993.
[6] Z.Zainuddin, and Y.Abu Hassan, ''Improving the Convergence of the Backpropagation Algorithm USING Local Adaptive Techniques'', Tnternational Journal of Computational Intelligence, Vol.
1, No.
3, 2004.
[7] Zaiyong Tang and Gary J.KOEHLER, ''A Convergent Neural Network Learning Algorithm'', IEEE, 1992.