حول التوقف المبكر لطريقة الانحدار العشوائي المرآتي لمشاكل عكسية غير محددة
On early stopping of stochastic mirror descent method for ill-posed inverse problems

المجلة: Numerische Mathematik، المجلد: 157، العدد: 2
DOI: https://doi.org/10.1007/s00211-025-01458-7
تاريخ النشر: 2025-03-06
المؤلف: Jing Huang وآخرون
الموضوع الرئيسي: طرق عددية في المشاكل العكسية

نظرة عامة

في السنوات الأخيرة، ظهرت الخوارزميات العشوائية كأدوات فعالة لمعالجة المشاكل العكسية غير المحددة، خاصة بسبب قدرتها على اختيار مجموعات عشوائية من المعادلات في كل تكرار، مما يعزز قابلية التوسع والأداء في السياقات الكبيرة. ومع ذلك، تواجه هذه الخوارزميات تحديات بسبب الطبيعة غير المحددة للمشاكل والضوضاء في بيانات القياس، مما يؤدي إلى تذبذبات كبيرة وظاهرة شبه التقارب النموذجية لطرق التنظيم التكرارية. هذا يعقد تحقيق المخرجات ذات الخصائص التقريبية المرغوبة.

لمعالجة هذه التحديات، يقدم هذا البحث قاعدة إيقاف بعدية لطريقة الانحدار العشوائي، مستوحاة من مبدأ التباين، مصممة خصيصًا للمشاكل العكسية غير المحددة في فضاءات باناش. تضمن قاعدة الإيقاف المقترحة أن تنتهي الخوارزمية في عدد محدود من الخطوات، مع تقارب الحل الناتج تقريبًا إلى الحل الحقيقي مع انخفاض مستوى الضوضاء. يقدم المؤلفون محاكاة عددية توضح فعالية هذا النهج، مع تسليط الضوء على إمكانيته في تحسين موثوقية الحلول في سياق الخوارزميات العشوائية للمشاكل العكسية.

مقدمة

تناقش المقدمة المشاكل العكسية غير المحددة الممثلة بالنظام \( F_i(x) = y_i \) لـ \( i = 1, \ldots, N \)، حيث \( F_i \) هي مشغلين يربطون من فضاء باناش المشترك \( X \) إلى فضاءات باناش الملساء بشكل موحد \( Y_i \). يشير المؤلفون إلى أن الحلول الدقيقة غالبًا ما تكون غير قابلة للتحقيق بسبب وجود الضوضاء في البيانات، مما يؤدي إلى ضرورة استخدام تقنيات التنظيم. يبرزون طريقة الانحدار العشوائي كنهج قابل للتطبيق لحل هذه المشاكل، خاصة في سياق البيانات الم noisy، ويؤكدون على أهمية طرق التحسين العشوائية التي تتعامل بكفاءة مع مجموعات البيانات الكبيرة.

يهدف البحث إلى تعزيز طريقة الانحدار العشوائي من خلال تقديم قاعدة إيقاف بعدية لتحديد متى يجب إنهاء التكرارات. يتم تحفيز ذلك من خلال الملاحظة بأن الطريقة يمكن أن تحقق حدود متبقية مرضية، \( F_i(x_n) – y_i \leq \tau \delta_i \)، للمعلمات المختارة بشكل مناسب. يقترح المؤلفون استراتيجية للتحقق من هذه الحالة بشكل أقل تكرارًا لتقليل الحمل الحسابي مع الاستفادة من المعلومات التي تم الحصول عليها خلال التكرارات. يوضحون منهجيتهم، التي تجمع بين عناصر من نظرية فضاء باناش، التحليل المحدب، والاحتمالات، لتأسيس نتائج التقارب لنهجهم. هيكل البحث منظم في أقسام تغطي المواد الخلفية، تفاصيل الطريقة، تحليل التقارب، والمحاكاة العددية للتحقق من كفاءة الطريقة المقترحة.

نقاش

في هذا القسم، يناقش المؤلفون المفاهيم الأساسية من التحليل الوظيفي والمحدب التي تدعم المنهجية اللاحقة لحل الأنظمة غير المحددة باستخدام طريقة الانحدار العشوائي. يعرفون فضاء باناش \(X\) وثنائيه \(X^*\)، مع الخصائص ذات الصلة مثل السلاسة الموحدة ومفهوم المشتقات الفرعية. يقدم القسم مسافة بريغمان، التي تستخدم لقياس التباين بين النقاط في سياق الدوال المحدبة. كما يحدد المؤلفون الشروط التي يمكن بموجبها تحليل دالة صحيحة، شبه مستمرة من الأسفل، و\(p\)-محدبة \(R\)، مع التركيز بشكل خاص على نظير ليجاندر-فينشيل.

يستعرض قسم المنهجية خوارزمية الانحدار العشوائي، موضحًا الافتراضات اللازمة لتنفيذها، بما في ذلك خصائص فضاءات باناش المعنية ومحدودية الدالة \(R\). يقدم المؤلفون معيار إيقاف استنادًا إلى مبدأ التباين، والذي يضمن أن تنتهي الخوارزمية بعد عدد محدود من التكرارات مع احتمال مرتفع. يقترحون استراتيجية عملية لتنفيذ الخوارزمية، والتي تتضمن التناوب بين خطوات الانحدار العشوائي وخطوة انحدار كاملة لتعزيز التقارب. يختتم القسم بوصف هيكل الخوارزمية، مع التأكيد على كفاءتها وتأثيرات التنظيم التي تسهم بها كل من الحلقات الداخلية والخارجية للعملية التكرارية.

Journal: Numerische Mathematik, Volume: 157, Issue: 2
DOI: https://doi.org/10.1007/s00211-025-01458-7
Publication Date: 2025-03-06
Author(s): Jing Huang et al.
Primary Topic: Numerical methods in inverse problems

Overview

In recent years, stochastic algorithms have emerged as effective tools for addressing ill-posed inverse problems, particularly due to their ability to select random subsets of equations at each iteration, which enhances scalability and performance in large-scale contexts. However, these algorithms are challenged by the ill-posed nature of the problems and noise in measurement data, leading to significant oscillations and a semi-convergence phenomenon typical of iterative regularization methods. This complicates the achievement of outputs with desirable approximation properties.

To address these challenges, this paper introduces an a posteriori stopping rule for the stochastic mirror descent method, inspired by the discrepancy principle, specifically tailored for ill-posed inverse problems in Banach spaces. The proposed stopping rule guarantees that the algorithm terminates in a finite number of steps, with the resulting solution converging almost surely to the true solution as the noise level diminishes. The authors provide numerical simulations that illustrate the effectiveness of this approach, highlighting its potential for improving the reliability of solutions in the context of stochastic algorithms for inverse problems.

Introduction

The introduction discusses ill-posed inverse problems represented by the system \( F_i(x) = y_i \) for \( i = 1, \ldots, N \), where \( F_i \) are operators mapping from a common Banach space \( X \) to uniformly smooth Banach spaces \( Y_i \). The authors note that exact solutions are often unattainable due to the presence of noise in the data, leading to the necessity for regularization techniques. They highlight the mirror descent method as a viable approach for solving these problems, particularly in the context of noisy data, and emphasize the importance of stochastic optimization methods that efficiently handle large datasets.

The paper aims to enhance the stochastic mirror descent method by introducing an a posteriori stopping rule to determine when to terminate iterations. This is motivated by the observation that the method can achieve satisfactory residual bounds, \( F_i(x_n) – y_i \leq \tau \delta_i \), for appropriately chosen parameters. The authors propose a strategy to check this condition less frequently to reduce computational load while still leveraging the information obtained during iterations. They outline their methodology, which combines elements from Banach space theory, convex analysis, and probability, to establish convergence results for their approach. The structure of the paper is organized into sections covering background material, method details, convergence analysis, and numerical simulations to validate the proposed method’s efficiency.

Discussion

In this section, the authors discuss foundational concepts from functional and convex analysis that underpin the subsequent methodology for solving ill-posed systems using the stochastic mirror descent method. They define a Banach space \(X\) and its dual \(X^*\), along with relevant properties such as uniform smoothness and the concept of subdifferentials. The section introduces the Bregman distance, which is utilized to measure the divergence between points in the context of convex functions. The authors also establish conditions under which a proper, lower semi-continuous, and \(p\)-convex function \(R\) can be analyzed, particularly focusing on its Legendre-Fenchel conjugate.

The methodology section outlines the stochastic mirror descent algorithm, detailing the assumptions necessary for its implementation, including the properties of the Banach spaces involved and the convexity of the function \(R\). The authors present a stopping criterion based on the discrepancy principle, which ensures that the algorithm terminates after a finite number of iterations with high probability. They propose a practical strategy for implementing the algorithm, which involves alternating between the stochastic mirror descent steps and a full mirror descent step to enhance convergence. The section concludes with a description of the algorithm’s structure, emphasizing its efficiency and the regularization effects contributed by both the inner and outer loops of the iterative process.