الگوریتم اولیه- دوگانه توزیع شده با پارامترهای متغیر و ساختار مشارکتی افزایشی دوجهته
محورهای موضوعی : مهندسی برق و کامپیوتر
1 - دانشگاه صنعتی ارومیه
کلید واژه: الگوریتم پارامتر متغیر, بازسازی توزیع‌شده, حسگری فشرده, مد افزایشی دوجهته,
چکیده مقاله :
به دلیل شرایط خاص شبکههای حسگری بیسیم از نقطه‌نظرهایی نظیر محدودیت انرژی، تسریع سرعت همگرایی الگوریتم‌های این حوزه اهمیت پیدا می‌کند. این امر در مورد حسگری فشرده توزیع‌شده که فاز بازسازی پیچیده‌ای دارد، ضروری‌تر به نظر می‌رسد. بر همین اساس در این مقاله، الگوریتم بازسازی حسگری فشرده توزیع‌شده‌ای ارائه می‌شود که امکان بازسازی با نرخ همگرایی بهبودیافته‌تری را میسر می‌سازد. الگوریتم پیشنهادی، یک الگوریتم اولیه- دوگانه توزیع‌شده در یک ساختار افزایشی دوجهته است که در آن پارامترها با زمان تغییر می‌کنند. تغییرات پارامترها به‌صورت ضابطه‌مند و برای آن دسته از مسائل بهینه‌سازی محدبی انجام می‌گیرد که در آنها توابعی که بیان‌کننده قید مسئله و مدل‌کننده مشارکت بین گره‌ها هستند، قویاً محدب می‌باشند. شیوه پیشنهادی با شبیه‌سازی‌هایی تضمین شده که نشان از عملکرد بالای الگوریتم پیشنهادی به لحاظ سرعت همگرایی، حتی در شرایط سختگیرانهتری نظیر تعداد اندک اندازه‌گیری‌ها و یا درجه تنکی پایین‌تر دارد.
Special conditions of wireless sensor networks, such as energy limitation, make it essential to accelerate the convergence of algorithms in this field, especially in the distributed compressive sensing (DCS) scenarios, which have a complex reconstruction phase. This paper presents a DCS reconstruction algorithm that provides a higher convergence rate. The proposed algorithm is a distributed primal-dual algorithm in a bidirectional incremental cooperation mode where the parameters change with time. The parameters are changed systematically in the convex optimization problems in which the constraint and cooperation functions are strongly convex. The proposed method is supported by simulations, which show the higher performance of the proposed algorithm in terms of convergence rate, even in stricter conditions such as the small number of measurements or the lower degree of sparsity.
[1] ح. نظری، م. رئیس‌دانایی و م. سپه‌وند، "مکان‌یابی بر اساس تفاضل توان سیگنال دریافتی با به‌کارگیری بهینه‌سازی محدب در شبکه حسگر بیسیم،" نشريه مهندسی برق و مهندسی کامپيوتر ايران، ب- مهندسي كامپيوتر، سال 17، شماره 4، صص. 310-306، زمستان 1398.
[2] ف. پدیداران مقدم و ح. مقصودی، "مسیریابی بهبودیافته برای توازن بار در شبکه حسگر بیسیم در بستر اینترنت اشیا بر پایه الگوریتم کلونی مورچگان چندگانه،" فصلنامه فناوری اطلاعات و ارتباطات ایران، سال 14، شماره 52/51، صص. 255-220، تابستان 1401.
[3] ق. آذرنیا، م. ع. طینتی و ت. یوسفی رضایی، "الگوریتم توزیع‌شده و مشارکتی به‌منظور بازسازی سیگنالهای تنک در شبکههای حسگری بیسیم با توپولوژی افزایشی دوجهته،" فصلنامه علمی پردازش علائم و داده‌ها، جلد 18، شماره 3، صص. 76-65، زمستان 1400.
[4] ق. آذرنیا و ع. ع. شریفی، "الگوریتم طول متغیر کسری نفوذی با قابلیت پیشبرد مشارکت مبتنی بر گرادیان برای افزایش کارایی شبکههای تطبیقی با لینک‌های نویزی،" مجله علمی پژوهشی رایانش نرم و فناوری اطلاعات، جلد 10، شماره 4، صص. 14-1، زمستان 1400.
[5] S. El Khediri, "Wireless sensor networks: a survey, categorization, main issues, and future orientations for clustering protocols," Computing, vol. 104, no. 8, pp. 1775-1837, Aug. 2022.
[6] G. Azarnia and A. A. Sharifi, "Steady-state analysis of distributed incremental variable fractional tap-length LMS adaptive networks," Wireless Networks, vol. 27, no. 7, pp. 4603-4614, 2021.
[7] G. Azarnia, "Diffusion fractional tap-length algorithm with adaptive error width and step-size," Circuits, Systems, and Signal Processing, vol. 41, no. 1, pp. 321-345, Jan. 2022.
[8] G. Azarnia, M. A. Tinati, A. A. Sharifi, and H. Shiri, "Incremental and diffusion compressive sensing strategies over distributed networks," Digital Signal Processing, vol. 101, Article ID: 102732, Jun. 2020.
[9] C. Li, G. Li, and P. K. Varshney, "Distributed detection of sparse signals with censoring sensors in clustered sensor networks," Information Fusion, vol. 83/84, pp. 1-18, Jul. 2022.
[10] G. Azarnia and A. A. Sharifi, "Fully cooperative and distributed focal underdetermined system solver compressive sensing recovery algorithm for wireless sensor networks," International J. of Communication Systems, vol. 35, no. 9, Article ID: e5126, Jun. 2022.
[11] O. Karpis, et al., "Compressed sensing-a way to spare energy in WSN for UAV," IFAC-PapersOnLine, vol. 55, no. 4, pp. 170-176, 2022.
[12] A. Salim, W. Osamy, A. M. Khedr, A. Aziz, and M. Abdel-Mageed, "A secure data gathering scheme based on properties of primes and compressive sensing for IoT-based WSNs," IEEE Sensors J., vol. 21, no. 4, pp. 5553-5571, 15 Feb. 2020.
[13] Y. Xu, G. Sun, T. Geng, and B. Zheng, "Compressive sparse data gathering with low-rank and total variation in wireless sensor networks," IEEE Access, vol. 7, pp. 155242-155250, 2019.
[14] C. Wang, X. Shen, H. Wang, and H. Mei, "Energy-efficient collection scheme based on compressive sensing in underwater wireless sensor networks for environment monitoring over fading channels," Digital Signal Processing, vol. 127, Article ID: 103530, Jul. 2022.
[15] M. Al Mazaideh and J. Levendovszky, "A multi-hop routing algorithm for WSNs based on compressive sensing and multiple objective genetic algorithm," J. of Communications and Networks, vol. 23, no. 2, pp. 138-147, Apr. 2021.
[16] J. Chen, et al., "Factor graphs for support identification in compressive sensing aided wireless sensor networks," IEEE Sensors J., vol. 21, no. 23, pp. 27195-27207, 1 Dec. 2021.
[17] R. Torkamani, H. Zayyani, and R. A. Sadeghzadeh, "Model-based decentralized Bayesian algorithm for distributed compressed sensing," Signal Processing: Image Communication, vol. 95, Article ID: 116212, Jul. 2021.
[18] F. Amini, Y. Hedayati, and H. Zanddizari, "Exploiting the inter-correlation of structural vibration signals for data loss recovery: a distributed compressive sensing-based approach," Mechanical Systems and Signal Processing, vol. 152, Article ID: 107473, May 2021.
[19] G. Azarnia, "Distribution agnostic Bayesian compressive sensing with incremental support estimation," Multidimensional Systems and Signal Processing, vol. 33, no. 2, pp. 327-340, Jun. 2022.
[20] A. Chambolle and T. Pock, "A first-order primal-dual algorithm for convex problems with applications to imaging," J. of Mathematical Imaging and Vision, vol. 40, no. 1, pp. 120-145, May 2011.
نشریه مهندسی برق و مهندسی کامپیوتر ایران، ب- مهندسی کامپیوتر، سال 21، شماره 2، تابستان 1402 129
مقاله پژوهشی
الگوریتم اولیه- دوگانه توزیع‌شده با پارامترهای متغیر
و ساختار مشارکتی افزایشی دوجهته
قنبر آذرنیا
چکیده: به دلیل شرایط خاص شبکههای حسگری بیسیم از نقطه‌نظرهایی نظیر محدودیت انرژی، تسریع سرعت همگرایی الگوریتم‌های این حوزه اهمیت پیدا می‌کند. این امر در مورد حسگری فشرده توزیع‌شده که فاز بازسازی پیچیده‌ای دارد، ضروری‌تر به نظر می‌رسد. بر همین اساس در این مقاله، الگوریتم بازسازی حسگری فشرده توزیع‌شده‌ای ارائه می‌شود که امکان بازسازی با نرخ همگرایی بهبودیافته‌تری را میسر می‌سازد. الگوریتم پیشنهادی، یک الگوریتم اولیه- دوگانه توزیع‌شده در یک ساختار افزایشی دوجهته است که در آن پارامترها با زمان تغییر می‌کنند. تغییرات پارامترها به‌صورت ضابطه‌مند و برای آن دسته از مسائل بهینه‌سازی محدبی انجام می‌گیرد که در آنها توابعی که بیان‌کننده قید مسئله و مدل‌کننده مشارکت بین گره‌ها هستند، قویاً محدب می‌باشند. شیوه پیشنهادی با شبیه‌سازی‌هایی تضمین شده که نشان از عملکرد بالای الگوریتم پیشنهادی به لحاظ سرعت همگرایی، حتی در شرایط سختگیرانهتری نظیر تعداد اندک اندازه‌گیری‌ها و یا درجه تنکی پایین‌تر دارد.
کلیدواژه: الگوریتم پارامتر متغیر، بازسازی توزیع‌شده، حسگری فشرده، مد افزایشی دوجهته.
1- مقدمه
محدودیت‌های توان و پهنای باند در شبکههای حسگری بیسیم، طراحی‌ها، پیکربندی‌ها و الگوریتم‌های مرتبط با چنین شبکه‌هایی را به چالش کشیده است [1] تا [5]. برای رفع این چالش، پردازش‌های توزیع‌شده و مشارکتی [6] و [7] در قالب توپولوژی‌های مختلف [8] در شبکههای حسگری بیسیم مورد توجه قرار گرفته است.
در برخی از کاربردهای شبکههای حسگری بیسیم، سیگنال‌ها اغلب نوعی ساختار را شامل می‌شوند که پردازش و نمایش هوشمند را برای آنها ممکن می‌سازد. ساختاری که اخیراً در این قبیل شبکه‌ها مورد توجه قرار گرفته است، تنکی2 سیگنال مورد پردازش می‌باشد [9]. سیگنال تنک، سیگنالی است که اکثر درایه‌های آن صفر باشد. در یک بیان دیگر، یک سیگنال 
تنک نامیده می‌شود، اگر حداکثر 
درایه غیرصفر را شامل شود. این ویژگی از سیگنال به استفاده توزیع‌شده از تکنیک حسگری فشرده3 در شبکههای حسگری بیسیم انجامیده است [10] و [11]. از جمله کاربردهای حسگری فشرده در شبکههای حسگری بیسیم می‌توان به کاربرد آن در جمع‌آوری داده و مسیریابی اشاره کرد [12] تا [15].
در [16] حسگری فشرده در شبکههای حسگری بیسیم از نقطه‌نظر تشخیص ساپورت (مکان ورودی‌های غیرصفر در یک سیگنال تنک) مورد توجه قرار گرفته است.
در [17] یک الگوریتم حسگری فشرده توزیع‌شده بیزین مبتنی بر مدل پیشنهاد شده که از هر دو همبستگی درون سیگنالی و بین سیگنالی برای بازیابی توأم چندین سیگنال تنک استفاده می‌کند.
در [18] بازیابی داده‌های ازدست‌رفته چندین سیگنال با نسبت تلفات متغیر بررسی شده و برای این منظور، [18] همبستگی بین سیگنالی را مورد توجه قرار داده و از حسگری فشرده توزیع‌شده بهره گرفته است.
در [19] یک روش بیزین برای بازیابی توزیع‌شده سیگنالهای تنک ارائه گردیده که نسبت به توزیع ساپورت، وفقی بوده و نیازی به معلوم‌بودن توزیع ساپورت ندارد. در شیوه پیشنهادی [19] در شبکه‌ای از حسگرها، بازسازی یک سیگنال تنک به‌طور مشارکتی و بر اساس مد همکاری افزایشی انجام می‌گیرد. در این روش، هر گره معیار تشخیص ساپورت غالب را از گره قبلی خود دریافت، مقدار معیار تشخیص ساپورت غالب خود را به آن اضافه و نتیجه را به گره بعدی ارسال می‌کند. در آخرین گره حلقه، داده‌های کل شبکه برای تصمیم‌گیری در مورد ساپورت غالب مورد استفاده قرار می‌گیرد که این بهره‌گیری از دایورسیتی مکانی در فرایند تشخیص ساپورت عملکرد برتری نتیجه می‌دهد.
در [3] یک الگوریتم بازسازی حسگری فشرده توزیع‌شده با مشارکت همه گره‌های شبکه ارائه گردیده که در آن، عملکرد حالت دائم بهتر و پیچیدگی پایین محاسباتی، مهم‌ترین مشخصه این الگوریتم به شمار می‌رود. الگوریتم پیشنهادی [3] یک الگوریتم بازگشتی اولیه- دوگانه است که با بهینه‌سازی ضرایب لاگرانژ حاصل شده است. این الگوریتم در چهارچوب مد مشارکتی افزایشی دوجهته ارائه گردیده که در آن هر حسگر در هر تکرار تنها با دو گره همسایه خود به تبادل داده می‌پردازد. لذا در مقایسه با ساختار مشارکتی نفوذی [10] که در آن هر حسگر با تمام گره‌های همسایه خود به تبادل داده می‌پردازد، شیوه پیشنهادی [3] به میزان قابل ملاحظه‌ای مصرف توان پایین‌تری دارد. شیوه پیشنهادی [3] یک چهارچوب توزیع‌شده جامع برای بازسازی سیگنالهای تنک در شبکههای حسگری بیسیم است که می‌توان آن را در قالب مسائل بهینه‌سازی محدب متفاوتی پیاده کرد. مزیت مسائل محدب نسبت به مسائل غیرمحدب این است که در مسائل بهینه‌سازی محدب، یک بهینه مطلق را می‌توان در اکثر موارد با دقت خوب در مدت زمانی معقول و مستقل از مقداردهی اولیه محاسبه کرد.
نقطه ضعف الگوریتم ارائه‌شده در [3] وجود پارامترهایی است که نحوه انتخاب آنها تأثیر مستقیمی بر همگرایی و سرعت همگرایی این الگوریتم دارد. در اینجا قصد داریم تا نحوه انتخاب این پارامترها را در راستای تسریع‌بخشیدن به سرعت همگرایی این الگوریتم فرمول‌بندی کنیم. در این مقاله، تحلیلی ریاضی پیرامون شیوه ارائه‌شده در [3] را ارائه خواهیم داد. در نتیجه تحلیل انجام‌گرفته، اولاً شرایط لازم برای همگرایی این الگوریتم استخراج خواهد شد و ثانیاً روابطی بازگشتی برای پارامترهای اساسی این الگوریتم ارائه خواهیم داد که انتخاب پارامترها بر اساس این روابط، همگرایی الگوریتم را تسریع خواهد نمود. شیوه پارامتر متغیر پیشنهادی می‌تواند در قالب تمام آن دسته از مسائل بهینه‌سازی پیاده شود که در آنها توابعی که بیان‌کننده قید مسئله و مدل‌کننده مشارکت بین گره‌ها هستند، قویاً محدب4 باشند.
مطالعات ما نشان دادند که در [20]، الگوریتمی برای تسریع همگرایی مسائل بهینه‌سازی محدب اولیه- دوگانه و برای مسائلی که درجه‌ای از همواری را شامل می‌شوند، ارائه گردیده است. اما شیوه [20] با روش پیشنهادی ما در این مقاله تفاوت‌های اساسی دارد. نخست اینکه الگوریتم پیشنهادی [20] شیوه توزیع‌شده‌ای نبوده و برای بازسازی یک سیگنال به‌صورت غیرمشارکتی در یک گره منفرد می‌تواند کاربرد داشته باشد؛ حال آنکه شیوه پیشنهادی ما در قالب شبکههای حسگری بیان شده است. همین امر، تحلیل‌های صورت‌گرفته در این مقاله را پرچالش‌تر کرده و باید در تحلیل‌های ریاضی انجام‌گرفته، اثر مشارکت بین گره‌ها را نیز منظور کنیم. همچنین در نتیجه مشارکت بین گره‌ها پارامتر دیگری وارد الگوریتم پیشنهادی [3] شده که استخراج قیدهایی برای این پارامتر به‌منظور تضمین همگرایی الگوریتم و یافتن رابطه‌ای برای این پارامتر در کنار دیگر پارامترها برای تسریع همگرایی الگوریتم، خود چالش دیگری است.
2- مسئله حسگری فشرده
2-1 بیان مسئله
از جمله مسائل عملی پرکاربرد در پردازش سیگنال، مسئله بازسازی سیگنال 
از داده مشاهده‌شده 
است که در آن ماتریس حسگری 
فرایند اندازه‌گیری خطی را مدل می‌کند. مادامی که 
باشد، این مسئله یک مسئله مرسوم در جبر کلاسیک خطی بوده و با روش‌های کلاسیک قابل حل است. اما زمانی که تعداد مشاهدات 
کمتر از طول 
سیگنال 
باشد، بدون اطلاعات اضافی‌تر، بازسازی از مشاهدات ممکن نخواهد بود. یکی از فرضیاتی که بازسازی سیگنال مطلوب را در چنین شرایطی ممکن می‌سازد، تنکی این سیگنال است به این معنا که اکثر درایه‌های صفر باشد. تحت فرض تنکی، الگوریتم 
تنک‌ترین برداری را نتیجه خواهد داد که با داده مشاهده‌شده 
سازگار باشد. اما از آنجایی که این الگوریتم در حالت کلی یک مسئله NP-hard است، الگوریتم‌های کارایی با قابلیت پیاده‌سازی در حسگری فشرده ارائه گردیده‌اند. این الگوریتم‌ها را می‌توان در چند دسته جای داد: شیوه‌های بهینه‌سازی محدب [10]، روش‌های حریصانه، روش‌های مبتنی بر آستانه [8]، شیوه‌های بهینه‌سازی غیرمحدب، الگوریتم‌های ترکیبی یا زیرخطی و روش‌های آماری [17].
وقتی مسئله حسگری فشرده در چهارچوب شبکههای سنسوری بیسیم مطرح می‌شود، الگوریتم‌های بازسازی سازگاریافته برای این چهارچوب، علاوه بر اینکه باید به‌طور منطقی سریع باشند، باید قابلیت پیاده‌سازی مشارکتی و توزیعی را دارا باشند؛ موردی که بر چالش طرح الگوریتم‌های بازسازی در این چهارچوب می‌افزاید. یکی از چنین الگوریتم‌های کارایی در [3] ارائه شده که در ادامه به معرفی این الگوریتم خواهیم پرداخت.
2-2 الگوریتم اولیه- دوگانه توزیعی با توپولوژی افزایشی دوجهته
یک شبکه حسگری بیسیم شامل 
حسگر را در نظر می‌گیریم. مأموریت این شبکه آن است که هر حسگر با داشتن بردار اندازه‌گیری 
و ماتریس اندازه‌گیری 
سیگنال تنک را به نحوی بازسازی کند که 
با 
برقرار باشد. برای این منظور، [3] مسئله بهینه‌سازی زیر را در نظر گرفته است
(1)
که 
و 
توابع محدب، نیمه‌پیوسته پایین‌تر5 و حقیقی‌مقدار تعمیم‌یافته6 هستند. همچنین متغیر محلی 
بهعنوان تخمینی برای سیگنال در گره 
در نظر گرفته شده و فرض گردیده که 
است. مسئله (1) را می‌توان به‌صورت زیر بازنویسی کرد
(2)
لاگرانژین مسئله (2) عبارت است از
(3)
که در آن 
. مقدار تابع لاگرانژ در نقاط بهینه آن، همان مقدار بهینه مسئله saddle-point زیر است
(4)
در (4)، 
و 
به ترتیب مزدوج محدب7 توابع 
و 
هستند. مزدوج محدب یک تابع نوعی مثل به‌صورت زیر تعریف می‌شود
(5)
در [3] یک رویکرد بازگشتی برای حل مسئله (4) به‌صورت زیر آمده است
(6)
در این الگوریتم 
نگاشت پروگزیمال8 تابع 
بوده و به‌صورت زیر تعریف می‌شود
(7)
همچنین نگاشت پروگزیمال و می‌تواند مستقیماً توسط (7) و یا با استفاده از اتحاد Moreau به‌صورت زیر حساب شود
(8)
در الگوریتم پیشنهادی [3] نحوه انتخاب پارامترهای 
، 
، 
و 
تأثیر مستقیمی بر همگرایی و سرعت همگرایی این الگوریتم دارد. به همین دلیل در ادامه الگوریتمی با پارامترهای متغیر با زمان 
، 
، 
و 
ارائه خواهیم کرد که انتخاب پارامترها بر اساس این روابط همگرایی الگوریتم را تسریع خواهد بخشید. در پیشنهاد این الگوریتم چنین فرض کرده‌ایم که و قویاً محدب باشند.
3- الگوریتم پارامتر متغیر پیشنهادی
در این بخش به ارائه الگوریتم پارامتر متغیری خواهیم پرداخت که می‌تواند عملکرد الگوریتم [3] را ارتقا ببخشد. به‌منظور ارائه الگوریتم پیشنهادی، نخست به ارائه لم و قضیه‌ای خواهیم پرداخت که پایه‌های منطقی برای ارائه الگوریتم پیشنهادی را بنا می‌کنند.
لم 1: فرض کنیم که و هر دو قویاً محدب با پارامترهای تحدب قوی 
و 
باشند. در این صورت با درنظرگرفتن پارامترهای متغیر با زمان ، ، و ، تکرارهای (6) رابطه زیر را در نقطه بهینه 
برآورده خواهند کرد
(9)
اثبات: به پیوست الف مراجعه شود.
قضیه 1: فرض کنیم که و هر دو قویاً محدب با پارامترهای تحدب قوی و باشند. انتخاب‌های زیر را انجام می‌دهیم
(10)
و قیدهای زیر را اعمال می‌کنیم
(11)
در این صورت تکرارهای 
از الگوریتم (6) کراندار خواهند ماند.
اثبات: به پیوست ب مراجعه شود.
قضیه 1، برهان لازم برای ارائه الگوریتم پیشنهادی را فراهم می‌آورد به این معنا که می‌توان پارامترهای ، ، و را مطابق (10) متغیر با زمان انتخاب کرد، اما هنوز همگرایی الگوریتم (6) را از دست نداد، مشروط بر اینکه قیدهای (11) رعایت شوند. بر این اساس، الگوریتم پارامتر متغیر پیشنهادی به‌صورت زیر ارائه می‌شود
(12)
4- ارزیابی و شبیه‌سازی الگوریتم پیشنهادی
در این بخش به ارزیابی عملکرد الگوریتم پیشنهادی می‌پردازیم. برای این منظور مدل تنکی زیر را در نظر می‌گیریم
(13)
اولین عبارت (13) تنکی سیگنال، دومین عبارت قید مسئله بهینه‌سازی و سومین عبارت آن مشارکت بین هر گره و همسایه آن گره را مدل می‌کند. این مسئله معادل مسئله (1) و (2) است که در آن 
،
و 
می‌باشد.
مزدوج محدب و عبارتند از
(14)
و
(15)
همان طور که ملاحظه می‌شود، هر دو تابع و قویاً محدب بوده و پارامتر تحدبشان به ترتیب عبارت است از 
و 
. حال به محاسبه نگاشت‌های پروگزیمال لازم در (12) می‌پردازیم. نگاشت پروگزیمال عبارت خواهد بود از 
با درایه 
ام 
که در آن 
عملگر آستانه نرم مختلط9 بوده و به‌صورت زیر تعریف می‌شود
(16)
که در آن 
است. نگاشت پروگزیمال توابع و
نیز به ترتیب عبارت هستند از 
و 
.
با توجه به نگاشت پروگزیمال ، عبارت متناظر با به‌هنگام‌رسانی 
در (12) به‌صورت 
خواهد بود. با توجه به وجود 
در این عبارت، اینکه گره ام در هر تکرار به ماتریس اندازه‌گیری گره همسایه خود نیاز دارد، شاید بهعنوان یک مشکل تلقی شود، اما با تعریف متغیر جدید 
این مبادله ماتریس مطابق الگوریتم زیر به تبادل یک بردار کاهش می‌یابد
(17)
در ادامه به شبیه‌سازی این الگوریتم پرداخته شده و آن را با شیوه پارامتر ثابت ارائه‌شده در [3] مقایسه می‌کنیم. در شبیه‌سازی‌ها این الگوریتم را
شكل 1: NMSE بر حسب تکرار برای 
.
الگوریتم بازسازی تنک توزیع‌شده افزایشی دوجهته پارامتر متغیر10 یا به اختصار VPDB-ISPARE و نقطه مقابل پارامتر ثابت آن را که نتیجه حل (13) بر طبق شیوه ارائه‌شده در [3] است، الگوریتم بازسازی تنک توزیع‌شده افزایشی دوجهته11 یا به اختصار DIBI-SPARE می‌نامیم. برای ارزیابی الگوریتم‌ها از خطای میانگین مربعات نسبی استفاده خواهیم کرد که به‌صورت زیر تعریف می‌شود
(18)
که در آن 
نرم فربنیوس ماتریس 
است. ماتریس از بار کنار هم قرار گرفتن بردار تنک ایجاد شده و 
ماتریسی است که ستون ام آن، بردار بازسازی‌شده توسط گره است. برای تحقق میانگین موجود در معیار ارزیابی (18) نتایج 100 آزمایش مستقل از هم متوسط‌گیری شده‌اند. در هر آزمایش، درایه‌های ماتریس اندازه‌گیری 
با توزیع گوسی استاندارد، ایجاد و سپس ستون‌های این ماتریس نرمالیزه شده‌اند. به‌علاوه در هر آزمایش ساپورت بردار تنک به‌طور تصادفی و با توزیع یکنواخت انتخاب شده که مؤلفه‌های غیرصفر آن از توزیع گوسی استاندارد تبعیت می‌کنند. همچنین مقداردهی اولیه شبیه‌سازی‌ها به‌صورت 
، 
، 
و 
است. برای پارامتر مقدار اولیه 
انتخاب شده است. مقدار اولیه پارامترهای و طوری انتخاب شده‌اند که قید (11) برآورده گردد. مقادیر این پارامترها در الگوریتم [3] نیز برابر همین مقادیر اولیه انتخاب شده‌اند. در الگوریتم [3] برای مقدار واحد انتخاب شده و مقادیر پارامترهای دیگر نیز عبارت هستند از 
، 
، 
و 
.
شکل 1 منحنی NMSE بر حسب تکرار را برای 
نشان می‌دهد که در آن 
بوده و 
و 
قرار داده شده‌اند. همان طور که مشاهده می‌شود، VPDB-ISPARE حدود 20 تکرار سریع‌تر از DIBI-SPARE همگرا می‌شود که به‌وضوح نشان از عملکرد همگرایی بالای الگوریتم پیشنهادی دارد.
تغییرات و بر حسب تکرار و برای گره 
در شکل 2 نشان داده شده است. منحنی تغییرات شبیه بوده و لذا در این شکل
شكل 2: منحنی (a) 
و (b) 
بر حسب تکرار و برای گره 
.
شكل 3: NMSE بر حسب تکرار برای 
.
شكل 4: NMSE بر حسب تکرار برای 
.
نمایش داده نشده است. همان طور که مشاهده می‌شود روند نزولی و روند صعودی داشته و این شیوه تغییرات به خوبی توانسته که مطابق شکل 1 موجب تسریع سرعت همگرایی شود.
شکل‌های 3 و 4 منحنی NMSE بر حسب تکرار را به ترتیب برای 
و 
نشان می‌دهند که در آن بوده و 
و 
قرار داده شده‌اند. همان طور که مشاهده می‌شود، الگوریتم VPDB-ISPARE برای حدود 60 تکرار و برای حدود 190 تکرار سریع‌تر از DIBI-SPARE همگرا می‌شود. به عبارتی از شکل‌های 3 و 4 می‌توان نتیجه گرفت که اولاً
شكل 5: NMSE بر حسب تکرار برای 
.
شكل 6: NMSE بر حسب تکرار برای 
.
الگوریتم پیشنهادی در مقادیر خیلی کوچک‌تر نیز همگرا شده و ثانیاً اگرچه با کاهش سرعت همگرایی کاهش پیدا کرده است، اما تأثیر این کاهش سرعت بر الگوریتم پیشنهادی به مراتب کمتر از الگوریتم DIBI-SPARE بوده است.
در شبیه‌سازی‌های قبل، پارامتر تنکی برابر در نظر گرفته شده است. برای بررسی عملکرد الگوریتم پیشنهادی بر حسب پارامترهای تنکی مختلف، منحنی NMSE بر حسب تکرار برای 
و 
به ترتیب در شکل‌های 5 و 6 ترسیم گردیده که در آنها بوده و 
و 
قرار داده شده‌اند. همان طور که از این شکل‌ها مشاهده می‌شود، حتی در شرایط سختگیرانه‌تر برای پارامتر تنکی، الگوریتم پیشنهادی عملکرد به‌مراتب بهتری از خود به نمایش می‌گذارد.
5- نتیجه
در اين مقاله، یک الگوریتم توزیع‌شده اولیه- دوگانه با پارامترهای متغیر و در یک ساختار افزایشی دوجهته که در آن هر گره می‌تواند در یک ساختار چرخه‌ای از گره ماقبل خود دیتا دریافت کرده و نیز به آن دیتا ارسال کند، ارائه شده است. انتخاب پارامترها به‌صورت متغیر با زمان و در یک ساختار نظام‌مند این امکان را در الگوریتم پیشنهادی فراهم آورد تا با سرعت بالاتری به حالت دائم خود همگرا شود. ضابطه‌های استخراج‌شده برای پارامترهای الگوریتم پیشنهادی در شرایط سختگیرانهتری نظیر پارامتر تنکی بزرگ‌تر یا تعداد اندازه‌گیری پایین‌تر نیز عملکرد همگرایی سریع‌تری را نتیجه می‌دهند.
[1] این مقاله در تاریخ 30 شهریور ماه 1401 دریافت و در تاریخ 17 دی ماه 1401 بازنگری شد.
قنبر آذرنیا (نويسنده مسئول)، دانشکده فنی و مهندسی خوی، دانشگاه صنعتی ارومیه، ارومیه، ایران، (email: g.azarnia@uut.ac.ir).
[2] . Sparsity
[3] . Compressive Sensing
[4] . Strongly Convex
[5] . Lower Semicontinuous
[6] . Extended Real-Valued
[7] . Convex Conjugate
[8] . Proximal Mapping
[9] . Complex Soft Thresholding
[10] . Variable Parameter Distributed Bidirectional Incremental Sparse Reconstruction
[11] . Distributed Bidirectional Incremental Sparse Reconstruction
(پ- 5)
(پ- 6)
پیوست
پیوست الف- اثبات لم 1
با درنظرگرفتن پارامترهای متغیر با زمان ، ، و ، تکرارهای الگوریتم (6) را می‌توان به‌صورت زیر بیان کرد
(پ- 1)
که در آن 
به زیردیفرانسیل1 تابع اشاره دارد. در (پ- 1) از این واقعیت استفاده شده که برای تابع محدب داریم 
، اگر و تنها اگر 
. با بهره‌گیری از مفهوم زیردیفرانسیل و با این فرض که و هر دو قویاً محدب هستند از (پ- 1) داریم
(پ- 2)
که در آن 
پارامترهای تحدب قوی هستند. با جمع طرفین سه نامساوی در (پ- 2) خواهیم داشت
(پ- 3)
با توجه به 
داریم
(پ- 4)
با جمع طرفین رابطه فوق به ازای 
، (پ- 5) را داریم. در (پ- 5) عبارت 
را جایگذاری کرده و 
را نقطه بهینه در نظر می‌گیریم تا (پ- 6) نتيجه گردد. در از این واقعیت استفاده شده که در نقطه بهینه
[1] . Subdifferential
(پ- 7)
(پ- 12)
(پ- 13)
(پ- 16)
تفاضل دو عبارت اول سمت راست نامساوی (پ- 5) با توجه به (4) نامنفی خواهد بود. رابطه (پ- 6) را می‌توان به‌صورت (پ- 7) بازنویسی کرد که در آن از این واقعیت استفاده شده که 
، 
و 
. با توجه به 
و 
برای هر 
داریم
(پ- 8)
به‌طور مشابه داریم
(پ- 9)
با جایگذاری (پ- 8) و (پ- 9) در (پ- 7) و با اندکی محاسبات، (9) نتیجه خواهد شد.
پیوست ب- اثبات قضیه 1
از (10) داریم: 
، 
، 
و 
. از طرفی از قیدهای (11) و (10) همواره خواهیم داشت
(پ- 10)
و از این رابطه می‌توان نتیجه گرفت که
(پ- 11)
لذا با توجه به لم 1، (9) را می‌توان به‌صورت (پ- 12) نوشت. طرفین (پ- 12) را از 
تا 
جمع می‌بندیم. با فرض 
، (پ- 13) را خواهیم داشت و از طرفی با توجه به
(پ- 14)
و
(پ- 15)
می‌توان (پ- 13) را به‌صورت (پ- 16) نوشت. با توجه به (پ- 10) تمام عبارت‌های سمت راست نامساوی (پ- 16) مثبت می‌باشد و لذا با انتخاب پارامترهای متغیر با زمان مطابق (10) و با اعمال قید (11) تکرارهای کراندار می‌مانند.
مراجع
[1] ح. نظری، م. رئیس‌دانایی و م. سپه‌وند، "مکان‌یابی بر اساس تفاضل توان سیگنال دریافتی با به‌کارگیری بهینه‌سازی محدب در شبکه حسگر بیسیم،" نشريه مهندسی برق و مهندسی کامپيوتر ايران، ب- مهندسي كامپيوتر، سال 17، شماره 4، صص. 310-306، زمستان 1398.
[2] ف. پدیداران مقدم و ح. مقصودی، "مسیریابی بهبودیافته برای توازن بار در شبکه حسگر بیسیم در بستر اینترنت اشیا بر پایه الگوریتم کلونی مورچگان چندگانه،" فصلنامه فناوری اطلاعات و ارتباطات ایران، سال 14، شماره 52/51، صص. 255-220، تابستان 1401.
[3] ق. آذرنیا، م. ع. طینتی و ت. یوسفی رضایی، "الگوریتم توزیع‌شده و مشارکتی به‌منظور بازسازی سیگنالهای تنک در شبکههای حسگری بیسیم با توپولوژی افزایشی دوجهته،" فصلنامه علمی پردازش علائم و داده‌ها، جلد 18، شماره 3، صص. 76-65، زمستان 1400.
[4] ق. آذرنیا و ع. ع. شریفی، "الگوریتم طول متغیر کسری نفوذی با قابلیت پیشبرد مشارکت مبتنی بر گرادیان برای افزایش کارایی شبکههای تطبیقی با لینک‌های نویزی،" مجله علمی پژوهشی رایانش نرم و فناوری اطلاعات، جلد 10، شماره 4، صص. 14-1، زمستان 1400.
[5] S. El Khediri, "Wireless sensor networks: a survey, categorization, main issues, and future orientations for clustering protocols," Computing, vol. 104, no. 8, pp. 1775-1837, Aug. 2022.
[6] G. Azarnia and A. A. Sharifi, "Steady-state analysis of distributed incremental variable fractional tap-length LMS adaptive networks," Wireless Networks, vol. 27, no. 7, pp. 4603-4614, 2021.
[7] G. Azarnia, "Diffusion fractional tap-length algorithm with adaptive error width and step-size," Circuits, Systems, and Signal Processing, vol. 41, no. 1, pp. 321-345, Jan. 2022.
[8] G. Azarnia, M. A. Tinati, A. A. Sharifi, and H. Shiri, "Incremental and diffusion compressive sensing strategies over distributed networks," Digital Signal Processing, vol. 101, Article ID: 102732, Jun. 2020.
[9] C. Li, G. Li, and P. K. Varshney, "Distributed detection of sparse signals with censoring sensors in clustered sensor networks," Information Fusion, vol. 83/84, pp. 1-18, Jul. 2022.
[10] G. Azarnia and A. A. Sharifi, "Fully cooperative and distributed focal underdetermined system solver compressive sensing recovery algorithm for wireless sensor networks," International J. of Communication Systems, vol. 35, no. 9, Article ID: e5126, Jun. 2022.
[11] O. Karpis, et al., "Compressed sensing-a way to spare energy in WSN for UAV," IFAC-PapersOnLine, vol. 55, no. 4,
pp. 170-176, 2022.
[12] A. Salim, W. Osamy, A. M. Khedr, A. Aziz, and M. Abdel-Mageed, "A secure data gathering scheme based on properties of primes and compressive sensing for IoT-based WSNs," IEEE Sensors J., vol. 21, no. 4, pp. 5553-5571, 15 Feb. 2020.
[13] Y. Xu, G. Sun, T. Geng, and B. Zheng, "Compressive sparse data gathering with low-rank and total variation in wireless sensor networks," IEEE Access, vol. 7, pp. 155242-155250, 2019.
[14] C. Wang, X. Shen, H. Wang, and H. Mei, "Energy-efficient collection scheme based on compressive sensing in underwater wireless sensor networks for environment monitoring over fading channels," Digital Signal Processing, vol. 127, Article ID: 103530, Jul. 2022.
[15] M. Al Mazaideh and J. Levendovszky, "A multi-hop routing algorithm for WSNs based on compressive sensing and multiple objective genetic algorithm," J. of Communications and Networks, vol. 23, no. 2, pp. 138-147, Apr. 2021.
[16] J. Chen, et al., "Factor graphs for support identification in compressive sensing aided wireless sensor networks," IEEE Sensors J., vol. 21, no. 23, pp. 27195-27207, 1 Dec. 2021.
[17] R. Torkamani, H. Zayyani, and R. A. Sadeghzadeh, "Model-based decentralized Bayesian algorithm for distributed compressed sensing," Signal Processing: Image Communication, vol. 95, Article ID: 116212, Jul. 2021.
[18] F. Amini, Y. Hedayati, and H. Zanddizari, "Exploiting the inter-correlation of structural vibration signals for data loss recovery:
a distributed compressive sensing-based approach," Mechanical Systems and Signal Processing, vol. 152, Article ID: 107473, May 2021.
[19] G. Azarnia, "Distribution agnostic Bayesian compressive sensing with incremental support estimation," Multidimensional Systems and Signal Processing, vol. 33, no. 2, pp. 327-340, Jun. 2022.
[20] A. Chambolle and T. Pock, "A first-order primal-dual algorithm for convex problems with applications to imaging," J. of Mathematical Imaging and Vision, vol. 40, no. 1, pp. 120-145, May 2011.
قنبر آذرنیا تحصیلات خود را در مقاطع کارشناسی و کارشناسی ارشد در رشته مهندسی برق مخابرات سیستم در دانشگاه تبريز به پايان رساند؛ سپس بهعنوان دانش آموخته برتر تحصیلات خود را در مقطع دکترای تخصصی رشته مهندسی برق مخابرات سیستم ادامه داده و در سال 1396 موفق به کسب درجه دکترای تخصصی از دانشگاه تبريز شد. ايشان هماکنون استاديار دانشکده فنی و مهندسی خوی دانشگاه صنعتی ارومیه هستند.
زمینههای پژوهشی مورد علاقه ايشان عبارتند از: پردازش توزيعشده در شبکههای حسگری بیسیم، تخمین کانال، پردازش سیگنالهای پزشکی، حسگری فشرده، فیلتر وفقی و تبديل موجک.