پیوندهای مفید
اخبار و اعلانات
آمار
| تعداد نشریات | 44 |
| تعداد شمارهها | 1,877 |
| تعداد مقالات | 15,278 |
| تعداد مشاهده مقاله | 43,731,473 |
| تعداد دریافت فایل اصل مقاله | 17,563,688 |
کاربرد الگوریتم کاوش باکتری دینامیکی تطبیقی برای پخش بار بهینه با در نظر گرفتن تأثیر نقطهی شیر و توان بادی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| هوش محاسباتی در مهندسی برق | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| مقاله 1، دوره 3، شماره 4، بهمن 1391، صفحه 1-14 اصل مقاله (222.69 K) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| نوع مقاله: مقاله پژوهشی فارسی | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| نویسندگان | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| احسان افضلان1؛ محمود جورابیان* 2 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 1- کارشناس ارشد مهندسی برق قدرت- | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 2استاد، دانشکده فنی و مهندسی- دانشگاه شهید چمران- اهواز- ایران | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| چکیده | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| این مطالعه به ارائه الگوریتم کاوش باکتری دینامیکی تطبیقی (DABFA) برای حل مسأله پخش بار اقتصادی (ED) با در نظر گرفتن تأثیر نقطه شیر و تلفات توان میپردازد. به علاوه توان باد نیز در فرمولسازی مسأله در نظر گرفته شده است. منابع انرژی تجدید پذیر و به خصوص انرژی بادی، اخیراً به علت ملاحظات محیطی و اقتصادی گوناگون جلب توجه زیادی نمودهاند. الگوریتم (BFA) اولیه، روش بهینه سازی تکاملی اخیراً پیشرفت یافتهای است که با الگوگیری از رفتار کاوشی باکتری Escherichia coli به دست آمده است. روش BFA ابتدایی به طرز موفقیت آمیزی برای حل مسائل بهینه سازی کوچک اجرا شده است، با این وصف، این روش مشخصه همگرایی ضعیفی برای مسائل مقید بزرگتر نشان میدهد. برای حل مشکل پیچیدگی و بعد زیاد فضای جستجوی مسأله، اصلاحاتی لازم برای بهبود عملکرد الگوریتم معرفی شدهاند. گام پایهای به گونهای تنظیم شده که رفتار دینامیکی غیر خطی داشته باشد تا بین جستجوی عمومی و محلی تعادل برقرار شود. معیار توقف روش BFA ابتدایی نیز اصلاح شده تا به صورت تطبیقی وابسته به پیشرفت حل مسأله باشد، به جای این که توقف مسأله وابسته به ماکزیمم تعیین شده تعداد تکرارها باشد. الگوریتم پیشنهادی با استفاده از تعدادی سیستم تست صحهگذاری شده است. نتایج به دست آمده با نتایج سایر الگوریتمهای اعمال شده برای حل مسأله با در نظر گرفتن اثر نقطه شیر و تلفات توان همراه با توان بادی مقایسه شدهاند. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| کلیدواژهها | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| الگوریتم کاوش باکتری دینامیکی تطبیقی (DABFA)؛ مسأله پخش بار اقتصادی (ED) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| اصل مقاله | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
بهینه سازی پخش بار اقتصادی دههها به عنوان یک مسأله بهرهبرداری اقتصادی سیستم قدرت موضوع مطالعات تحقیقاتی وسیعی بوده است. هدف مسأله، کمینهسازی تابع هزینه سوخت واحدهای تولیدی به روش تخصیص بار کل مابین این واحدهاست. تابع هزینه کمینهسازی شده دارای تعدادی قیود است که شامل تعادل بار، حدود بالا و پایین واحدهای تولیدی و نرخ شیبهاست [1]. در روشهای بهینهسازی کلاسیک، مسأله توسط یک تابع هدف مسطح مشتقپذیر مربعی یا تکهتکه مربعی تخمین زده میشود. با وجود این به علت اثرات نقطه شیر، مشخصههای بالاتر و ناپیوستگی هستند که به یک تابع هزینه سوخت غیر محدب و غیر مسطح منجر میشود، اثرات نقطه شیر با استفاده از دو رویکرد مختلف ارائه شدهاند [2]: در اولین رویکرد این اثرات به صورت قیود نامساوی که آنها را به صورت نقاط عملکرد ممنوعه نشان میدهند، فرمول بندی شدهاند [3،4]. رویکرد دوم که در این مقاله استفاده میشود، حاوی یک جمله سینوسی یکسو شده در تابع هدف اولیه است تا این اثرات را مدل کند [5]. به علت افزایش شدید نگرانیهای محیطی و اقتصادی، انرژی بادی به عنوان یک انرژی تجدید پذیر مطمئن جلب توجه زیادی نموده است [6]. در نظر گرفتن انرژی بادی در شبکههای قدرت چالشهای زیادی را در استراتژیهای بهرهبرداری و برنامهریزی ایجاد میکند. این حقیقت که توان بادی نه به سادگی قابل پیشبینی است و نه قابلیت توزیع دارد، علت اصلی نگرانیهای امنیت و قابلیت اطمینان در سیستمهای قدرت دارای انرژی بادی است. انرژی بادی اخیراً موضوع تحقیقاتی در حال تکاملی بوده است [7]. بسته به مشخصههای بادی محلی و پروفایل بار، سهم توان بادی تولید شده میتواند 10 تا 40% از تقاضای بار باشد. برخی مطالعات سهم توان بادی را تا سال 2030 حدود 20 تا 30% پیشبینی کردهاند [6-7]. روشهای بهینهسازی قطعی گوناگونی برای حل مسأله کلاسیک استفاده شدهاند، از جمله روش تکراری لامدا [1](LIM)، جستجوی گرادیان، برنامهریزی خطی، برنامهریزی مربعی، برنامهریزی دینامیکی و روشهای بر پایه نیوتون [1، 8]. با اینکه این روشهای سنتی بر پایه حسابان عملکرد مناسبی برای حل مسأله ED کلاسیک از خود نشان میدهند، ولی در حل مسائل ED با در نظر گرفتن اثرات نقطهی شیر موفقیت قابل قبولی به دست نمیآورند [9]. اخیراً روشهای مکاشفهای غیر کلاسیک برای حل این مسأله پیشنهاد شدهاند. آنها شامل برنامهریزی تکاملی [10]، الگوریتم ژنتیک [11]، جستجوی Tabu [12] و بهینه سازی جستجوی ذرات[2](PSO) [13] هستند. این روشهای غیر وابسته به مشتق، عملکرد خوبی در حل مسألهی ED صرفنظر از شکل غیر خطی و غیر مسطح مشخصههای ورودی- خروجی واحد تولیدی حرارتی از خود نشان میدهند [9]. روش الگوریتم کاوش باکتری (BFA) تکنیک بهینهسازی اخیراً معرفی شدهای است که از رفتار کاوشی باکتری Escherichia coli الگو گرفته است [14]. روش BFA به طرز موفقیتآمیزی برای حل مسایل بهینهسازی مختلف، از جمله بهینهسازی و کنترل توزیع شده [14]، پخش بار بهینه [15،16]، طراحی پایدارسازهای بهینه سیستم قدرت [17] و تخمین هارمونیک [18] اعمال شده است. مشخصاً یک روش BFA ترکیبی اخیراً برای حل مسأله ED اعمال شده و به نتایج خوبی منجر شده است [19]. با وجود این، نتایج شبیهسازی نشان دادند که روش BFA دچار ویژگیهای همگرایی ضعیف بوده، به زمان طولانی نیاز دارد. این رفتار ضعیف در محیطهای دینامیکی و فضاهای جستجوی با ابعاد بالا مرتبط با مسائل پیچیده بدتر نیز میشود [20،21]. در این مقاله الگوریتم کاوش باکتری تطبیقی دینامیکی (DABFA) ارائه شده و برای حل مسأله ED با در نظر گرفتن اثرات نقطه شیر و تلفات توان اعمال شده است. مسأله ED به گونهای فرمولبندی شده که در معادله تعادل توان خود توان بادی را نیز در نظر بگیرد. در شبیهسازیها نفوذ انرژی بادی به اندازه 10% کل تقاضای بار در نظر گرفته شده است. چهار حالت سیستمهای تست برای صحهسنجی عملکرد روش DABFA پیشنهادی و مقایسه آن با تکنیکهای معروف دیگر بهینهسازی استفاده شدهاند. نتایج بیانگر دقت، سرعت و توانایی بالای الگوریتم پیشنهادی در یافتن پاسخهای بهینه نسبت به سایر روشهای مقایسه شده است. در ادامه این مقاله به این صورت ترتیببندی شده است: بخش2 فرمولبندی مسأله ED را ارائه میکند، در بخش3 الگوریتم DABFA توضیح داده شده است، نتایج شبیهسازی در بخش4 نشان داده شده و نتیجهگیریها در بخش5 آمدهاند.
1- فرمولبندی مسأله پخش بارمسأله ED برای تعیین بارگیری بهینه تمامی واحدهای در مدار بهمنظور کمینهسازی تابع هزینه مشروط به قیود سیستم طراحی شده است [1].
2-1- تابع هدف هزینه سوخت هر واحد تولیدی معمولاً به صورت تقریبی به صورت یک تابع مربعی مرتبه دوم تقریب زده میشود [22].
که در آن Fi(Pgi) تابع هزینه سوخت واحد i، Pgi توان تولیدی واحدi،ai ، bi و ci ضرایب هزینه سوخت واحد i هستند. برای بهدست آوردن یک مدل دقیق تابع هزینه، اثرات نقطه شیر به صورت زیر در تابع سوخت مدل میشوند [23].
که در آن ei وfi ضرایب ثابت هزینه سوخت واحد i با در نظر گرفتن اثرات نقطه شیر بوده و Pgimin مینیمم توان تولیدی واحد i است. تابع هزینه کل مورد کمینهسازی میتواند به صورت برداری زیر نوشته شود [24]:
که در آن FT تابع هزینه کل، Pg بردار ستونی توان خروجی واحد تولیدی، a مجموع بردار ضریب ثابت، b بردار ستونی ضریب خطی، c ماتریس قطری ضریب مربعی و Ng تعداد واحدهای تولیدی است.
2-2- قیود تابع هزینه نشان داده شده با رابطه (3) دارای تعدادی قید شامل موارد زیر است:
که در آن PL کل تلفات توان حقیقی سیستم و PD تقاضای بار کل سیستم است.
که در آن Pgimax بیشینهی توان تولیدی واحد i است.
تلفات خط انتقال به صورت تابعی از توان حقیقی و ماتریس ضرایب B بیان میشوند [25].
که در آن پارامترهای Bij ضرایب اتلاف نامیده میشوند. عبارت (6) که به عنوان فرمول تلفات Korn شناخته میشود [23] را میتوان به صورت برداری زیر بیان کرد [24].
که در آن B ماتریس متقارن ضرایب اتلاف، B0 بردار ستونی ضرایب اتلاف و B00 ثابت اسکالر ضریب اتلاف است. 2-3- نفوذ توان بادی تعادل توان با در نظر گرفتن توان باد: رابطهی (4) به گونهای اصلاح میشود که توان تولید شده توسط منابع بادی PW از تقاضای کل بار کم شود.
در دسترس بودن توان بادی: توان باد PWدر رابطه (8) توسط مقدار توان در دسترس از پارک بادی Pav محدود میشود.
3- روش الگوریتم کاوش باکتری دینامیکی تطبیقی پیشنهادی روش DABFA پیشنهادی بر اساس روش BFA اولیه بنا شده است. در قسمت زیر روش BFA اصلی توضیح داده شده است.
3-1- الگوریتم کاوش باکتری روش BFA یک تکنیک بهینه سازی تکاملی است که از رفتار کاوشی باکتری E.coli الهام گرفته شده است. جنبههای بیولوژیک روشهای کاوش باکتری و رفتار حرکتی آنها، علاوه بر مکانیزمهای تصمیمگیری آنها را میتوان در [14] یافت. به عنوان روشی مکاشفهای ، روش BFA به منظور غلبه بر مشکلات بهینهسازی بدون گرادیان و سر و کار داشتن با توابع هدف پیچیده و مشتق ناپذیر طراحی شده است. جستجوی ابرفضا از طریق سه عملکرد اصلی اجرا میشود: chemotaxis ، تولید مثل و فعالیتهای حذف- پراکندگی [14]. فرایند chemotaxis از طریق شنا کردن و غلت خوردن صورت میگیرد. باکتری طول عمر خود را با گذار از یکی از این مودهای حرکتی به مود دیگر طی میکند. در روش BFA، غلت خوردن با یک طول واحد در جهتی تصادفی ، نمایانده میشود، که جهت حرکت پس از غلت خوردن را مشخص میکند. اندازه گام برداشته شده در جهت تصادفی با ثابت واحد طول حرکت نشان داده میشود. برای جمعیتی از باکتریها، مکان i امین باکتری در j امین گام، k امین گام تولید مثل و l امین رخداد حذف/ پراکندگی با نشان داده میشود. در این مکان تابع هزینه با نشان داده میشود که به نام تابع مواد غذایی نیز شناخته میشود. پس از یک غلت خوردن، مکان i امین باکتری توسط رابطه زیر نشان داده میشود:
وقتی در تابع هزینه بهتر (پایینتر) از باشد، گام دیگری به اندازه در همان جهت طی میشود. عمل شنا کردن تا زمانی که هزینه کمتری به دست آید و به تعداد بیشینه گام از قبل معین شده برسیم تکرار میشود. تابع هزینه هر باکتری در جمعیت متأثر است از نوع ایجاد ازدحام اجرا شده توسط ارتباط سلول به سلول آزاد شده توسط باکتریها. برای شکل دادن به الگوهای ازدحامی عملیات هجوم آوردن به صورت زیر بیان میشود:
که در آن و و و ضرایب نمایانگر مشخصههای سیگنالهای جذبی و دفعی آزاد شده توسط سلول بوده و مؤلفه m ام باکتری i ام است. مکان هر عضو جمعیت S باکتری بوده، به این صورت تعریف میشود:
که S اندازه جمعیت باکتریهاست. تابع (11) که نمایانگر تأثیر ارتباط سلول به سلول است، با تابع هزینه جمع میشود:
فرایند تولید مثل پس از طی ماکزیمم تعداد گامهای chemotactic، ، اجرا میشود، جمعیت به دو نیمه میشود به گونهای که نیمه کمتر سالم میمیرد و هر باکتری در آن نیمه دیگر سالمتر به دو باکتری تقسیم میشود که همان مکان را اختیار میکنند.
پس از طی گام تولید مثل فرایند حذف/ پراکندگی به تعداد بار اجرا میشود. در این عملیات، هر باکتری میتواند برای کاوش قسمتهای دیگر فضای جستجو حرکت کند. احتمال این که هر باکتری اتفاق حذف/ پراکندگی را تجربه کند، با نسبت از پیش تعیین شده تعیین میشود.
3-2- الگوریتم کاوش باکتری دینامیکی تطبیقی در مورد روش BFA ابتدایی طول واحد گام، پارامتری ثابت است که میتواند در مورد مسائل بهینه سازی کوچک جستجوی مناسبی را تضمین نماید، اما در مورد مسایل پیچیده مقیاس بالا با ابعاد زیاد عملکرد ضعیفی از خود نشان میدهد. پارامتر طول حرکت عامل کلیدی در کنترل قابلیت جستجوی الگوریتم BFA است. از این دیدگاه توازن بین اکتشاف و بهره برداری جستجو میتواند با تنظیم واحد طول حرکت به دست آید. به خاطر داشته باشید که در سایر الگوریتمهای تکاملی برای کنترل قابلیت جستجوی محلی و کلی الگوریتم چند تابع دینامیکی اجرا میشوند [26-27]. در این مقاله به جای تابع ثابت پله یک تابع دینامیکی غیر خطی کاهشی برای اجرای حرکت شنا پیشنهاد داده میشود. این تابع به صورت زیر بیان میشود:
که در آن j طول و بیشینه تعداد گامهای chemotactic بوده و پارامتری از قبل تعیین شده است. معیار توقف اجرای الگوریتم BFA اصلی بیشینه گامهای chemotactic، تعداد گامهای تولید مثل و تعداد وقایع حذف/ پراکندگی است. این معیار موجب افزایش نیازمندی محاسباتی الگوریتم در برخی موارد میشود. در این مقاله یک معیار توقف تطبیقی اعمال شده تا الگوریتم بیشینه تعداد تکرارها را بسته به پیشرفت تابع هزینه تنظیم نماید. عملیات chemotaxis هنگامی متوقف میشود که یا پیشرفتی در حل مسأله حاصل نشود یا هنگامی که به بیشینه تعداد گامهای chemotactic برسیم. تکنیک الگوریتم پیشنهادی به صورت زیر است: گام 1) مقداردهی اولیه پارامترهای زیر: : بعد فضای جستجو : تعداد باکتریها در جمعیت : تعداد گامهای chemotactic : طول شنا وقتی روی شیب (مماس) باشد : عداد گامهای تولید مثل : تعداد وقایع حذف/ پراکندگی : احتمال این که هر باکتری حذف یا پراکنده شود. : واحد اولیه طول حرکت : واحد طول حرکت در پایان گامهای chemotactic (j=Nc ) : مکان تصادفی اولیه هر باکتری گام 2) حلقه حذف/ پراکندگی، گام 3) حلقه تولید مثل، گام 4) حلقه chemotaxis، . به ازای گام chemotactic را برای هر باکتری به صورت زیر اجرا کن:
تا زمانی که (تا مقدار بسیار طولانی مسیر طی نشود). قرار بده: . اگر باشد، آنگاه قرار بده سپس گامی دیگر را در همان جهت برداشته، مقدار جدید را محاسبه کن.
آنگاه (یعنی عملیات chemotactic را پایان بده). گام 5) اگر به گام چهارم برو ( ). گام 6) تولید مثل. برای k و l داده شده سلامتی هر باکتری i را به صورت زیر برآورد کن:
سلامتی باکتری i معیار این است که این باکتری در ظول زندگی خود چقدر مواد غذایی به دست آورده است.
گام 7) اگر به گام 3 برو( ). گام 8) حذف/ پراکندگی: با احتمال به صورت تصادفی هر باکتری i را حذف و پراکنده کن؛ به گونهای که اندازه جمعیت ثابت بماند. گام 9) اگر به گام 2برو( )،در غیر این صورت متوقف شو. فلوچارت الگوریتم پیشنهادی در شکل(1) ارائه شده است.
شکل(1): فلوچارت الگوریتم پیشنهادی
3-3- کارباقیود در مسأله پخش بار اقتصادی قیود تساوی و نامساوی به ترتیب با روابط (4) و (5) بیان میشوند. نوع دوم قیود (قیود نامساوی) بیانگر حدود مرزی تولید توان است. وقتی یک باکتری خارج از مرزهای فضای جستجوی خود شنا کند، قید نامساوی نقض شده است. در این حالت مدخل متناظر در بردار موقعیت باکتری دوباره در آخرین موقعیت مورد تحقق خود مکانیابی میشود. این در واقع اصل رد حلهای غیر قابل تحقق و نگهداری حلهای قابل تحقق است [29]. این روال به منظور حفظ ویژگیهای تصادفی الگوریتم و نگهداشتن هر باکتری در حالت فعال به جای اعمال جریمه مرگ، حتی اگر خارج از مرزها حرکت کند، اجرا میشود. قیود تساوی یعنی معادلات توازن بار از طریق فاکتور جریمه با تابع هدف تعمیم داده میشوند. معادله به دست آمده به نام تابع ارزیابی خوانده میشود. در الگوریتم پیشنهادی تابع جریمه به کار گرفته شده به صورت یک تابع دینامیکی غیر خطی افزایشی فرمول بندی میشود [30]. تابع ارزیابی پیشنهادی در روش DABFA که در هر تکرار محاسبه خواهد شد، به صورت ریاضی زیر بیان میشود:
جایی که ، و مقادیری ثابت بوده و j اندیسی برای گام chemotactic است. فرمول بندی نشان میدهد که قسمت اول تابع جریمه با افزایش تکرار بزرگتر میشود. یک انتخاب معقول برای ثوابت ، و میتواند به ترتیب 5/0 ، 2 و 2 باشد [30]. در بسیاری از شرایط عملی انتخاب پارامترها وابسته به مسأله است.
4- نتایج شبیهسازی الگوریتم DABFA پیشنهادی با استفاده از چهار مطالعه موردی تست شده است. برخی از این موارد نمایانگر مسأله پخش بار اقتصادی با در نظر گرفتن تلفات سیستم بوده و سایر آنها در تابع هدف خود شامل اثر نقطه شیر هستند. توان بادی نیز در هر مطالعه موردی در نظر گرفته شده است. الگوریتم در نرمافزار MATLAB 7.8و روی یک کامپیوتر شخصی Intel Core 2 Duo 1.66 GHzاجرا شده است. در هر تست به تعداد 30 اجرای برنامه به صورت مستقل با مقادیر متفاوت و تصادفی حل اولیه برای هر اجرا انجام شده است. نتایج به دست آمده در هر مورد با نتایج روشهای دیگر مقایسه شدهاند. مقایسه انجام شده با رویکردهای مختلفی همانند روشهای قطعی و مکاشفهای اجرا شده است. در همگی موارد، الگوریتم DABFA پیشنهادی یا عملکرد بهتری از سایر روشها دارد یا نتایج قابل مقایسهای به دست میدهد. مسأله پخش بار اقتصادی در هر مورد برای هر سیستم یک بار بدون در نظر گرفتن توان بادی و سپس با در نظر گرفتن یک ژنراتور بادی معادل حل شده است. در هریک از چهار مطالعه موردی سهم توان بادی تولید شده به میزان 10 درصد از کل تقاضای بار فرض شده است.
4-1- سیستم تست اول این مطالعه موردی سیستم 30 باسه IEEE با 6 ژنراتور و تقاضای کل بار MW1800 است [24،32]. تلفات سیستم با استفاده از فرمول تلفات (6) و (7) با ماتریس ضرایب تلفات زیر در نظر گرفته شدهاند [41].
الگوریتم DABFA پیشنهادی ابتدا بدون در نظر گرفتن توان تولیدی بادی اعمال میشود و سپس یک واحد معادل بادی با سهم 10 درصد از بار در نظر گرفته میشود. نتایج به دست آمده در هر دو زیر بخش در جدولهای1 و2 نشان داده شدهاند. نتایج به دست آمده در زیر بخش اول با نتایج الگوریتم(SWT-NR)[3] که در مرجع [31] مطرح شده، روش برنامهریزی مربعی متوالی [4](SQP) و بهینهسازی با جستجوی الگو مطرح شده [5](PS) در [24] مقایسه شدهاند. نتایج و مقایسه آنها در جدول 1 ارائه شدهاند.
جدول (1): مقایسه و نتایج شبیهسازی سیستم تست اول
جدول (2): نتایج شبیهسازی سیستم تست اول برای 30 بار اجرای برنامه
4-2- سیستم تست دوم در این حالت سیستم متشکل از 20 واحد تولیدی و تقاضای کل بار MW2500 است. دادههای سیستم شامل ماتریس ضرایب B را میتوان در [32] یافت. دو زیر بخش در نظر گرفته شدهاند: در اولی هیچ توان بادی وجود ندارد و در دومی یک ژنراتور بادی معادل که 10 درصد بار را تأمین میکند، در نظر گرفته شده است. نتایج به دست آمده در زیر بخش اول با نتایج روش Hopfield[32]، روش LIM و رویکرد PS[24] مقایسه شدهاند. جدولهای 3 و 4 نتایج اعمال الگوریتم پیشنهادی را برای هر دو زیر بخش به همراه مقایسهها برای زیر بخش اول نشان میدهند. در مقایسه با روشهای دیگر برای سیستم بدون توان بادی تابع هزینه به دست آمده با روش DABFA بسیار بهتر از توابع هزینه به دست آمده با مدل LIM و Hopfield بوده، مقداری کمتر از تابع هزینه به دست آمده از رویکرد PS است.
4-3- سیستم تست سوم این سیستم تست شامل 13 واحد تولیدی با تقاضای بار MW1800 و با در نظر گرفتن اثرات نقطه شیر است. دادههای سیستم در [10] ارائه شده و در جدول 5 نیز نشان داده شدهاند. این جدول همچنین پخش بار بهینه را در هر دو حالت سیستم با و بدون در نظر گرفتن توان بادی نشان میدهد. آمارهای همگرایی برای هر دو زیر بخش در جدول 6 نشان داده شدهاند. نتایج الگوریتم پیشنهادی برای تولید بدون توان بادی با نتایج روشهای دیگر مقایسه شدهاند. سایر روشها شامل برنامهریزی تکاملی [6](EP)[32]، بهینه سازی جستجوی ذرات[7](PSO)، برنامهریزی تکاملی ترکیبی با روش برنامهریزی مربعی متوالی(HEP-SQP)، روش ترکیبی جستجوی ذرات با روش برنامهریزی مربعی متوالی (HPSO-SQP) [9] و روش ترکیبی تکامل دیفرانسیلی آشوبی[8](CED) با روش برنامهریزی مربعی متوالی (CED-SQP) [33] هستند. جدول 7 این مقایسه را نشان میدهد. الگوریتم DABFA پیشنهادی به مقدار زیادی عملکرد بهتری از سایر روشها دارد، زیرا به مینیمم هزینه ($/h) 82/17845 منتج شده که به معنی صرفهجویی سالانه ($) 816000 در مقایسه با کمترین هزینه به دست آمده از سایر روشهای جدول 7 است.
جدول (3): نتایج مدلهای Hopfield، LIM، PSو DABFA برای سیستم تست دوم
جدول (4): نتایج شبیهسازی سیستم تست دوم برای 30 بار اجرای برنامه
جدول (5): اطلاعات سیستم شامل 13 واحد تولیدی درسیستم تست سوم
جدول (6): نتایج شبیهسازی سیستم تست سوم برای 30 بار اجرای برنامه
جدول (7): مقایسه و نتایج شبیهسازی برای سیستم تست سوم
4-4- سیستم تست چهارم در این مطالعه موردی تعداد واحدهای تولیدی 40 و تقاضای کل بار MW10500 بوده و اثرات نقطه شیر هم در نظر گرفته شده است. دادههای سیستم در[10] قابل دسترسی هستند. الگوریتم DABFA در این سیستم یک بار بدون توان بادی و بار دیگر با در نظر گرفتن توان بادی آزمایش شده است. تولید بهینه هر واحد برای هر دو زیر بخش در جدول 8 نشان داده شده است. خلاصهای از رفتار همگرایی 30 اجرای الگوریتم برای دو زیر بخش در جدول 9 نشان داده شده است. نتایج برای حالت بدون توان بادی با روشهای فهرست شده در مورد تست سوم به علاوه روش PSO اصلاح شده (MPSO) ارائه شده در [34] مقایسه شدهاند. نتایج مقایسه در جدول 10 نشان داده شدهاند. کمترین هزینه به دست آمده با الگوریتم DABFA پیشنهادی ($/h) 1970/119898 بوده که کمتر از هزینههای به دست آمده با سایر روشهای جدول 10 است. مقدار سالانه کاهش هزینه با استفاده از روش پیشنهادی ($) 320/16151530 است. این مقدار با مقایسه نمودن هزینه روش پیشنهادی و کمترین هزینه به دست آمده از سایر روشهای جدول 10 محاسبه شده است.
جدول (8): حل بهینه سیستم آزمون چهارم با 40 واحد تولیدی
جدول (9): نتایج شبیهسازی سیستم آزمون چهارم برای 30 بار اجرای برنامه
جدول (10): مقایسه و نتایج شبیهسازی برای سیستم آزمون چهارم
این حالت (سیستم تست با 40 واحد تولیدی) نشان میدهد که الگوریتم DABFA پیشنهادی با وجود ابعاد بالای ابر فضای جستجو و ناهمواری تابع هدف مربعی ناشی از اثرات نقطه شیر، عملکرد مناسبی در مسأله پخش بار اقتصادی دارد.
5- نتیجهگیری در این مقاله یک الگوریتم DABFA برای پخش بار اقتصادی بهینه معرفی شد. تأثیر در نظر گرفتن توان بادی در سیستمهای تولیدی و روی پخش بار اقتصادی مطالعه گردید. فرمولبندی مسأله پخش بار اقتصادی با در نظر گرفتن اثرات نقطه شیر و تلفات توان ارائه شد. به علاوه، فرمولبندی مسأله توان بادی را نیز به منظور نشان دادن و بررسی مزایای اقتصادی قرار دادن توان بادی در شبکههای سیستم قدرت مورد نظر قرار داد. الگوریتم پیشنهادی یک تکنیک کاوش باکتری اصلاح شده است که یک روش دینامیکی غیر خطی را برای به روز کردن بردار جستجو و بهبود همگرایی الگوریتم اعمال میکند. همچنین، از یک معیار توقف تطبیقی برای کاهش بار محاسباتی الگوریتم استفاده میکند. نتایج شبیهسازی موثر و مقاوم بودن الگوریتم پیشنهادی را نشان دادهاند. مقایسه با سایر روشهای قطعی و مکاشفهای نشان میدهد که الگوریتم پیشنهادی در اکثر حالات به نتایج بسیار بهتری رسیده و در سایر حالات هم نتایج حاصله از روش پیشنهادی چندان بدتر نبوده است. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| مراجع | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
[1] WOOD, A.J., WOLLENBERG, B.F., ‘Power generation operation and control’ (Wiley, New York, USA, 1996, 2nd edn.)
[2] VICTOIRE, T.A.A., JEYAKUMAR, A.E., ‘Reserve constrained dynamic dispatch of units with valve-point effects’, IEEE Trans. Power Syst, Vol. 20, No. 3, pp. 1273-1282, 2005.
[3] LEE, F.N., BREIPOHL, A.M., ‘Reserve constrained economic dispatch with prohibited operating zones’, IEEE Trans.Power Syst, Vol. 8, No. 1, pp. 246-254, 1993.
[4] GAING, Z., ‘Particle swarm optimization to solving the economic dispatch considering the generator constraints’, IEEE Trans. Power Syst, Vol. 18, No. 3, pp. 1187–1195, 2003.
[5] ATTAVIRIYANUPAP, P., KITA, H., TANAKA, E., HASEGAWA, J., ‘A hybrid EP and SQP for dynamic economic dispatch with nonsmooth fuel cost function’, IEEE Power Eng. Rev, Vol. 22, No. 4, pp. 77-77, 2002.
[6] SMITH, J., THRESHER, R., ZAVADIL, R., ET, A.L., ‘A mighty wind’, IEEE Power Energy Mag, Vol. 7, No. 2, pp. 41-51, 2009.
[7] DEMEO, E.A., GRANT, W., MILLIGAN, M.R., SCHUERGER, M.J., ‘Wind plant integration’, IEEE Power Energy Mag, Vol. 3, No. 6 , pp. 38-46, 2005.
[8] ARAVINDHABABU, P., NAYAR, K.R., ‘Economic dispatch based on optimal lambda using radial basis function network’, Int. J. Electr.Power Energy Syst, Vol. 24, No. 7, pp. 551-556, 2002.
[9] VICTOIRE, T.A.A., JEYAKUMAR, A.E., ‘Hybrid PSO-SQP for economic dispatch with valve-point effect’, Electr.Power Syst. Res, Vol. 71, No. 1, pp. 51-59, 2004.
[10] SINHA, N., CHAKRABARTI, R., CHATTOPADHYAY, P.K., ‘Evolutionary programming techniques for economic load dispatch’, IEEE Trans. Evol. Comput, Vol. 7, No. 1, pp. 83–94, 2003.
[11] TIPPAYACHAI, J., ONGSAKUL, W., NGAMROO, I., ‘Parallel micro genetic algorithm for constrained economic dispatch’, IEEE Trans. Power Syst, Vol. 17, No. 3, pp. 790–797, 2002.
[12] LIN, W., CHENG, F., TSAY, M., ‘An improved Tabu search for economic dispatch with multiple minima’, IEEE Trans.Power Syst, Vol. 17, No. 1, pp. 108-112, 2002.
[13] SELVAKUMAR, A.I., THANUSHKODI, K., ‘A new particle swarm optimization solution to nonconvex economic dispatch problems’, IEEE Trans. Power Syst, Vol. 22, No. 1, pp. 42-51, 2007.
[14] PASSINO, K.M., ‘Biomimicry of bacterial foraging for distributed optimization and control’, IEEE Control Syst.Mag, Vol. 22, No. 3, pp. 52-67, 2002.
[15] TANG, W.J., WU, Q.H., SAUNDERS, J.R., ‘Bacterial foraging algorithm for dynamic environments’.IEEE Congress Evolutionary Computation (CEC 2006), pp. 1324-1330, 2006.
[16] TANG, W.J., LI, M.S., HE, S., WU, Q.H., SAUNDERS, J.R., ‘Optimal power flow with dynamic loads using bacterial foraging algorithm’.Int. Conf. on Power System Technology (PowerCon 2006), pp.1-5, 2006.
[17] DAS, T.K., VENAYAGAMOORTHY, G.K., ALIYU, U.O., ‘Bio-inspired algorithms for the design of multiple optimal power system stabilizers: SPPSO and BFA’, IEEE Trans. Ind. Appl, Vol. 44, No. 5, pp. 1445-1457, 2008.
[18] MISHRA, S., ‘A hybrid least square-fuzzy bacterial foraging strategy for harmonic estimation’, IEEE Trans. Evol. Comput, Vol. 9, No. 1, pp. 61-73, 2005.
[19] PANIGRAHI, B.K., PANDI, V.R., ‘Bacterial foraging optimisation: Nelder-Mead hybrid algorithm for economic load dispatch’, IET Gener.Transm. Distrib, Vol. 2, No. 4, pp. 556-565, 2008.
[20] CHEN, H., ZHU, Y., HU, K., ‘Self-adaptation in bacterial foraging optimization algorithm’.Third Int. Conf. on Intelligent System and Knowledge Engineering (ISKE 2008), Vol. 1, pp. 1026-1031, 2008.
[21] CHU, Y., MI, H., LIAO, H., JI, Z., WU, Q.H., ‘A fast bacterial swarming algorithm for high-dimensional function optimization’. IEEE Congress Evolutionary Computation (IEEE World Congress on Computational Intelligence)(CEC2008), pp.3135-3140,2008.
[22] EL-HAWARY, M.E., CHRISTENSEN, G.S., ‘Optimal economic operation of electric power systems’ (Academic Press, New York, USA, 1979).
[23] FLETCHER, R., ‘Practical methods of optimization’ (Wiley, New York, USA, 2002, 2nd edn.).
[24] ALHAJRI, M.F., EL-HAWARY, M.E., ‘Pattern search optimization applied to convex and non-convex economic dispatch’.IEEE Int. Conf. on Systems, Man and Cybernetics (ISIC 2007), pp. 2674-2678, 2007.
[25] BERGEN, R., VITTAL, V., ‘Power systems analysis’ (Prentice-Hall, Upper Saddle River, NJ, USA, 1999, 2nd edn.).
[26] SHI, Y., EBERHART, R., ‘A modified particle swarm optimizer’.Proc. IEEE Int. Conf. on Evolutionary Computation, 1998(IEEE World Congress on Computational Intelligence), pp. 69-73, 1998.
[27] NAKA, S., GENJI, T., YURA, T., FUKUYAMA, Y., ‘Practical distribution state estimation using hybrid particle swarm optimization’.2001 IEEE Power Eng. Society Winter Meeting, Vol. 2,pp. 815-820, 2001.
[28] PERAM, T., VEERAMACHANENIM, K., | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
آمار تعداد مشاهده مقاله: 3,190 تعداد دریافت فایل اصل مقاله: 874 |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||