الانحدار التدريجي

الانحدار التدرجي في ثنائي الأبعاد

يُعدّ التدرج الهبوطي طريقةً للتحسين الرياضي غير المقيد . وهو خوارزمية تكرارية من الدرجة الأولى لتقليل دالة متعددة المتغيرات قابلة للتفاضل .

تعتمد الفكرة على اتخاذ خطوات متكررة في الاتجاه المعاكس لتدرج الدالة ( أو ما يقاربه) عند النقطة الحالية، لأن هذا هو اتجاه الانحدار الأسرع. وعلى العكس، فإن اتخاذ خطوات في اتجاه التدرج سيؤدي إلى مسار يُعظّم تلك الدالة؛ وتُعرف هذه العملية حينها باسم صعود التدرج . يجب عدم الخلط بين هبوط التدرج وخوارزميات البحث المحلي ، على الرغم من أن كليهما طريقتان تكراريتان للتحسين .

يُعدّ التدرج الهبوطي مفيدًا بشكل خاص في مجال التعلم الآلي والذكاء الاصطناعي لتقليل دالة التكلفة أو الخسارة. [ 1 ]

يُنسب الفضل في تطوير خوارزمية التدرج الانحداري عمومًا إلى أوغستين لويس كوشي ، الذي اقترحها لأول مرة عام 1847. [ 2 ] واقترح جاك هادامارد بشكل مستقل طريقة مماثلة عام 1907. [ 3 ] [ 4 ] ودرس هاسكل كاري خصائص تقاربها لمسائل التحسين غير الخطي لأول مرة عام 1944، [ 5 ] وأصبحت هذه الطريقة موضع دراسة واستخدام متزايدين في العقود اللاحقة. [ 6 ] [ 7 ]

يُعد التدرج العشوائي، وهو امتداد بسيط للتدرج الهبوطي ، بمثابة الخوارزمية الأساسية المستخدمة لتدريب معظم الشبكات العميقة اليوم.

وصف

رسم توضيحي لانحدار التدرج على سلسلة من مجموعات المستويات

يعتمد انحدار التدرج على ملاحظة أنه إذا كانت الدالة متعددة المتغيراتو(x){\displaystyle f(\mathbf {x} )}معرفة وقابلة للتفاضل في جوار نقطةأ{\displaystyle \mathbf {a} }، ثمو(x){\displaystyle f(\mathbf {x} )}يتناقص بشكل أسرع إذا انتقل المرء منأ{\displaystyle \mathbf {a} }في اتجاه التدرج السالب لـو{\displaystyle f}فيأ،أنا.هـ.،-و(أ){\displaystyle \mathbf {a} ,ie,-\nabla f(\mathbf {a} )}ويترتب على ذلك أنه إذا

أن+1=أن-ηو(أن){\displaystyle \mathbf {a} _{n+1}=\mathbf {a} _{n}-\eta \nabla f(\mathbf {a} _{n})}

لحجم خطوة صغير بما فيه الكفاية أو معدل تعلم صغير بما فيه الكفايةηR+{\displaystyle \eta \in \mathbb {R} _{+}}، ثم و(أن)و(أن+1){\displaystyle f(\mathbf {a_{n}} )\geq f(\mathbf {a_{n+1}} )}بعبارة أخرى، المصطلحηو(أ){\displaystyle \eta \nabla f(\mathbf {a} )}يتم طرحها منأ{\displaystyle \mathbf {a} }لأننا نريد التحرك عكس اتجاه التدرج، باتجاه الحد الأدنى المحلي. مع وضع هذه الملاحظة في الاعتبار، يبدأ المرء بتخمين.x0{\displaystyle \mathbf {x} _{0}}كحد أدنى محلي لـو{\displaystyle f}ويأخذ في الاعتبار التسلسلx0،x1،x2،...{\displaystyle \mathbf {x} _{0},\mathbf {x} _{1},\mathbf {x} _{2},\ldots }بحيث

xن+1=xن-ηنو(xن)، ن0.{\displaystyle \mathbf {x} _{n+1}=\mathbf {x} _{n}-\eta _{n}\nabla f(\mathbf {x} _{n}),\ n\geq 0.}

لدينا متتالية رتيبة

و(x0)و(x1)و(x2)،{\displaystyle f(\mathbf {x} _{0})\geq f(\mathbf {x} _{1})\geq f(\mathbf {x} _{2})\geq \cdots ,}

إذن التسلسل(xن){\displaystyle (\mathbf {x} _{n})}يتقارب إلى الحد الأدنى المحلي المطلوب. لاحظ أن قيمة حجم الخطوةη{\displaystyle \eta }يُسمح بتغييره في كل تكرار.

من الممكن ضمان التقارب إلى الحد الأدنى المحلي في ظل افتراضات معينة على الدالة.و{\displaystyle f}(على سبيل المثال،و{\displaystyle f}محدب وو{\displaystyle \nabla f}ليبشيتز ) وخيارات محددة منη{\displaystyle \eta }وتشمل هذه التسلسلات

ηن=|(xن-xن-1)[و(xن)-و(xن-1)]|و(xن)-و(xن-1)2{\displaystyle \eta _{n}={\frac {\left|\left(\mathbf {x} _{n}-\mathbf {x} _{n-1}\right)^{\top }\left[\nabla f(\mathbf {x} _{n})-\nabla f(\mathbf {x} _{n-1})\right]\right|}{\left\|\nabla f(\mathbf {x} _{n})-\nabla f(\mathbf {x} _{n-1})\right\|^{2}}}}

كما هو الحال في طريقة بارزيلاي-بورواين ، [ 8 ] [ 9 ] أو سلسلةηن{\displaystyle \eta _{n}}تحقيق شروط وولف (التي يمكن إيجادها باستخدام البحث الخطي ). عندما تكون الدالةو{\displaystyle f}إذا كانت محدبة ، فإن جميع القيم الصغرى المحلية هي أيضًا قيم صغرى عالمية، لذلك في هذه الحالة يمكن أن يتقارب انحدار التدرج إلى الحل العالمي.

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

تشبيه لفهم انحدار التدرج

الضباب في الجبال

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

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

اختيار حجم الخطوة واتجاه النزول

منذ استخدام حجم الخطوةη{\displaystyle \eta }إن صغر حجمها سيؤدي إلى إبطاء التقارب، وη{\displaystyle \eta }سيؤدي الحجم الكبير جدًا إلى تجاوز الحد والتباعد، لذا فإن إيجاد إعداد جيد لـη{\displaystyle \eta }تُعدّ هذه مشكلة عملية هامة. وقد دعا فيليب وولف أيضًا إلى استخدام "خيارات ذكية لاتجاه الانحدار" في الممارسة العملية. [ 10 ] في حين أن استخدام اتجاه ينحرف عن اتجاه الانحدار الأشد قد يبدو غير بديهي، إلا أن الفكرة تكمن في إمكانية تعويض الانحدار الأقل من خلال الحفاظ على الانحدار لمسافة أطول بكثير.

للتفكير في هذا الأمر رياضياً، ضع في اعتبارك اتجاهاً ما.صن{\displaystyle \mathbf {p} _{n}}وحجم الخطوةηن{\displaystyle \eta _{n}}ولننظر في التحديث الأكثر عمومية:

أن+1=أن-ηنصن{\displaystyle \mathbf {a} _{n+1}=\mathbf {a} _{n}-\eta _{n}\,\mathbf {p} _{n}}.

إيجاد إعدادات جيدة لـصن{\displaystyle \mathbf {p} _{n}}وηن{\displaystyle \eta _{n}}يتطلب الأمر بعض التفكير. أولاً، نريد أن يكون اتجاه التحديث نحو الأسفل. رياضياً، لنفترضθن{\displaystyle \theta _{n}}يشير إلى الزاوية بين-و(أن){\displaystyle -\nabla f(\mathbf {a_{n}} )}وصن{\displaystyle \mathbf {p} _{n}}وهذا يتطلب أنكوسθن>0.{\displaystyle \cos \theta _{n}>0.}لشرح الأمر بتفصيل أكبر، نحتاج إلى مزيد من المعلومات حول دالة الهدف التي نعمل على تحسينها. وذلك بافتراض ضعيف إلى حد ما أنو{\displaystyle f}إذا كانت قابلة للتفاضل باستمرار، فيمكننا إثبات ما يلي: [ 11 ]

تشير هذه المتباينة إلى أن المقدار الذي يمكننا من خلاله التأكد من الدالةو{\displaystyle f}يعتمد انخفاض قيمة على المفاضلة بين الحدين الموجودين بين قوسين مربعين. يقيس الحد الأول الزاوية بين اتجاه الانحدار والميل السالب. أما الحد الثاني فيقيس سرعة تغير الميل على طول اتجاه الانحدار.

من حيث المبدأ ، يمكن تحسين المتباينة ( 1 ) علىصن{\displaystyle \mathbf {p} _{n}}وηن{\displaystyle \eta _{n}}لاختيار حجم الخطوة الأمثل واتجاهها. تكمن المشكلة في أن حساب الحد الثاني بين قوسين مربعين يتطلب حسابو(أن-تηنصن){\displaystyle \nabla f(\mathbf {a} _{n}-t\eta _{n}\mathbf {p} _{n})}وتُعدّ عمليات حساب التدرج الإضافية مكلفة وغير مرغوب فيها عمومًا. ومن بين الحلول الممكنة لهذه المشكلة ما يلي:

  • تخلَّ عن فوائد اتجاه الهبوط الذكي عن طريق تحديدصن=و(أن){\displaystyle \mathbf {p} _{n}=\nabla f(\mathbf {a_{n}} )}واستخدم البحث الخطي للعثور على حجم خطوة مناسبγن{\displaystyle \gamma _{n}}مثل تلك التي تُحقق شروط وولف . يُعد البحث الخطي بالتراجع طريقةً أكثر اقتصادية لاختيار معدلات التعلم ، وهي طريقة تتمتع بضمانات نظرية جيدة ونتائج تجريبية موثوقة. تجدر الإشارة إلى أنه ليس من الضروري اختيار صن{\displaystyle \mathbf {p} _{n}}ليكون هو التدرج؛ أي اتجاه له حاصل ضرب داخلي موجب مع التدرج سيؤدي إلى انخفاض قيمة الدالة (لقيمة صغيرة بما فيه الكفاية لـηن{\displaystyle \eta _{n}}).
  • بافتراض أنو{\displaystyle f}إذا كانت قابلة للتفاضل مرتين، فاستخدم مصفوفة هيسيان الخاصة بها2و{\displaystyle \nabla ^{2}f}لتقديرو(أن-تηنصن)-و(أن)2تηن2و(أن)صن.{\displaystyle \|\nabla f(\mathbf {a} _{n}-t\eta _{n}\mathbf {p} _{n})-\nabla f(\mathbf {a} _{n})\|_{2}\approx \|t\eta _{n}\nabla ^{2}f(\mathbf {a} _{n})\mathbf {p} _{n}\|.}ثم اخترصن{\displaystyle \mathbf {p} _{n}}وηن{\displaystyle \eta _{n}}عن طريق تحسين عدم المساواة ( 1 ).
  • بافتراض أنو{\displaystyle \nabla f}إذا كانت دالة ليبشيتز ، فاستخدم ثابت ليبشيتز الخاص بها.ل{\displaystyle L}للتقييدو(أن-تηنصن)-و(أن)2لتηنصن.{\displaystyle \|\nabla f(\mathbf {a} _{n}-t\eta _{n}\mathbf {p} _{n})-\nabla f(\mathbf {a} _{n})\|_{2}\leq Lt\eta _{n}\|\mathbf {p} _{n}\|.}ثم اخترصن{\displaystyle \mathbf {p} _{n}}وηن{\displaystyle \eta _{n}}عن طريق تحسين عدم المساواة ( 1 ).
  • قم ببناء نموذج مخصص لـالأعلىت[0،1]و(أن-تηنصن)-و(أن)2و(أن)2{\displaystyle \max _{t\in [0,1]}{\frac {\|\nabla f(\mathbf {a} _{n}-t\eta _{n}\mathbf {p} _{n})-\nabla f(\mathbf {a} _{n})\|_{2}}{\|\nabla f(\mathbf {a} _{n})\|_{2}}}}لو{\displaystyle f}ثم اخترصن{\displaystyle \mathbf {p} _{n}}وηن{\displaystyle \eta _{n}}عن طريق تحسين عدم المساواة ( 1 ).
  • في ظل افتراضات أقوى بشأن الوظيفةو{\displaystyle f}مثل التحدب ، قد تكون هناك تقنيات أكثر تقدماً ممكنة.

عادةً، باتباع إحدى الوصفات المذكورة أعلاه، يمكن ضمان التقارب إلى الحد الأدنى المحلي. عندما تكون الدالةو{\displaystyle f}إذا كانت محدبة ، فإن جميع القيم الصغرى المحلية هي أيضًا قيم صغرى عالمية، لذلك في هذه الحالة يمكن أن يتقارب انحدار التدرج إلى الحل العالمي.

حل نظام خطي

تم تطبيق خوارزمية الانحدار الأسرع على مرشح وينر [ 12 ]

يمكن استخدام خوارزمية التدرج الهبوطي لحل نظام من المعادلات الخطية

أx-ب=0{\displaystyle \mathbf {A} \mathbf {x} -\mathbf {b} =0}

تمت إعادة صياغتها كمسألة تقليل تربيعي. إذا كانت مصفوفة النظامأ{\displaystyle \mathbf {A} }إذا كانت دالة الهدف حقيقية ومتناظرة وموجبة التحديد ، فإنها تُعرَّف على أنها دالة تربيعية، مع تقليل قيمة

و(x)=xأx-2xب،{\displaystyle f(\mathbf {x} )=\mathbf {x} ^{\top }\mathbf {A} \mathbf {x} -2\mathbf {x} ^{\top }\mathbf {b} ,}

لهذا السبب.

و(x)=2(أx-ب).{\displaystyle \nabla f(\mathbf {x} )=2(\mathbf {A} \mathbf {x} -\mathbf {b} ).}

بالنسبة لمصفوفة حقيقية عامةأ{\displaystyle \mathbf {A} }تُعرّف طريقة المربعات الصغرى الخطية

و(x)=أx-ب2.{\displaystyle f(\mathbf {x} )=\left\|\mathbf {A} \mathbf {x} -\mathbf {b} \right\|^{2}.}

في طريقة المربعات الصغرى الخطية التقليدية للأرقام الحقيقيةأ{\displaystyle \mathbf {A} }وب{\displaystyle \mathbf {b} }يتم استخدام المعيار الإقليدي، وفي هذه الحالة

و(x)=2أ(أx-ب).{\displaystyle \nabla f(\mathbf {x} )=2\mathbf {A} ^{\top }(\mathbf {A} \mathbf {x} -\mathbf {b} ).}

تقليل البحث الخطي ، وإيجاد حجم الخطوة الأمثل محليًاη{\displaystyle \eta }في كل تكرار، يمكن إجراؤها تحليليًا للدوال التربيعية، وصيغ صريحة للحل الأمثل محليًاη{\displaystyle \eta }معروفة. [ 6 ] [ 13 ]

على سبيل المثال، بالنسبة للمصفوفة الحقيقية المتناظرة والمصفوفة الموجبة المحددةأ{\displaystyle \mathbf {A} }، يمكن أن تكون الخوارزمية البسيطة كما يلي، [ 6 ]

كرر في الحلقة:ر:=ب-أxη:=رر/رأرx:=x+ηرلو رر إذا كانت صغيرة بما يكفي، فقم بالخروج من الحلقةنهاية حلقة التكراريعود x ونتيجة لذلك{\displaystyle {\begin{aligned}&{\text{repeat in the loop:}}\\&\qquad \mathbf {r} :=\mathbf {b} -\mathbf {Ax} \\&\qquad \eta  :={\mathbf {r} ^{\top }\mathbf {r} }/{\mathbf {r} ^{\top }\mathbf {Ar} }\\&\qquad \mathbf {x}  :=\mathbf {x} +\eta \mathbf {r} \\&\qquad {\hbox{إذا كانت }}\mathbf {r} ^{\top }\mathbf {r} {\text{ صغيرة بما يكفي، فاخرج من الحلقة}}\\&{\text{أكمل حلقة التكرار}}\\&{\text{أرجع }}\mathbf {x} {\text{ كنتيجة}}\end{aligned}}}

لتجنب الضرب فيأ{\displaystyle \mathbf {A} }نلاحظ مرتين في كل تكرار أنx:=x+ηر{\displaystyle \mathbf {x} :=\mathbf {x} +\eta \mathbf {r} } يعنير:=ر-ηأر{\displaystyle \mathbf {r} :=\mathbf {r} -\eta \mathbf {Ar} } ، مما يعطي الخوارزمية التقليدية، [ 14 ]

ر:=ب-أxكرر في الحلقة:η:=رر/رأرx:=x+ηرلو رر إذا كانت صغيرة بما يكفي، فقم بالخروج من الحلقةر:=ر-ηأرنهاية حلقة التكراريعود x ونتيجة لذلك{\displaystyle {\begin{aligned}&\mathbf {r} :=\mathbf {b} -\mathbf {Ax} \\&{\text{كرر في الحلقة:}}\\&\qquad \eta  :={\mathbf {r} ^{\top }\mathbf {r} }/{\mathbf {r} ^{\top }\mathbf {Ar} }\\&\qquad \mathbf {x}  :=\mathbf {x} +\eta \mathbf {r} \\&\qquad {\hbox{إذا كانت }}\mathbf {r} ^{\top }\mathbf {r} {\text{ صغيرة بما يكفي، فاخرج من الحلقة}}\\&\qquad \mathbf {r}  :=\mathbf {r} -\eta \mathbf {Ar} \\&{\text{نهاية حلقة التكرار}}\\&{\text{إرجاع }}\mathbf {x} {\text{ كنتيجة}}\end{aligned}}}
مسار التقارب لطريقة الانحدار الأسرع لـ A = [[2, 2], [2, 3]]

نادرًا ما تُستخدم هذه الطريقة لحل المعادلات الخطية، حيث تُعد طريقة التدرج المترافق من أكثر البدائل شيوعًا. ويتناسب عدد تكرارات انحدار التدرج عادةً مع رقم الحالة الطيفية.κ(أ){\displaystyle \kappa (\mathbf {A} )}مصفوفة النظامأ{\displaystyle \mathbf {A} }(نسبة القيم الذاتية القصوى إلى القيم الذاتية الدنيا لـأأ{\displaystyle \mathbf {A} ^{\top }\mathbf {A} }بينما يُحدد تقارب طريقة التدرج المترافق عادةً بالجذر التربيعي لرقم الحالة، أي أنها أسرع بكثير. يمكن لكلا الطريقتين الاستفادة من التكييف المسبق ، حيث قد يتطلب هبوط التدرج افتراضات أقل على المُكيِّف المسبق. [ 14 ]

السلوك الهندسي والتعامد المتبقي

في الانحدار الأشد المطبق على حلأx=ب{\displaystyle \mathbf {Ax} =\mathbf {b} }، أينأ{\displaystyle \mathbf {A} }إذا كانت متناظرة وموجبة التحديد، فإن متجهات البواقيرك=ب-أxك{\displaystyle \mathbf {r} _{k}=\mathbf {b} -\mathbf {A} \mathbf {x} _{k}}متعامدة عبر التكرارات:

رك+1،رك=0.{\displaystyle \langle \mathbf {r} _{k+1},\mathbf {r} _{k}\rangle =0.}

لأن كل خطوة تُتخذ في الاتجاه الأكثر انحدارًا، فإن خطوات الانحدار الأشد تتناوب بين اتجاهات محاذية للمحاور القصوى لمجموعات المستويات الممتدة. عندماκ(أ){\displaystyle \kappa (\mathbf {A} )}عندما يكون حجمها كبيرًا، ينتج عن ذلك مسار متعرج مميز. سوء التكييف لـأ{\displaystyle \mathbf {A} }يُعد السبب الرئيسي للتقارب البطيء، كما أن تعامد البواقي المتتالية يعزز هذا التناوب.

كما هو موضح في الصورة على اليمين، فإن الانحدار الأسرع يتقارب ببطء بسبب ارتفاع رقم الحالة لـأ{\displaystyle \mathbf {A} }وتجبر خاصية تعامد البواقي كل اتجاه جديد على تصحيح التجاوز الناتج عن الخطوة السابقة. والنتيجة هي مسار متعرج نحو الحل. يُعد هذا القصور أحد أسباب تفضيل طرق التدرج المترافق أو طرق التكييف المسبق. [ 15 ]

حل نظام غير خطي

يمكن استخدام خوارزمية التدرج الهبوطي لحل نظام من المعادلات غير الخطية . فيما يلي مثال يوضح كيفية استخدام هذه الخوارزمية لإيجاد قيم ثلاثة متغيرات مجهولة، x1 و x2 و x3 . يوضح هذا المثال دورة واحدة من خوارزمية التدرج الهبوطي.

لنفترض نظام المعادلات غير الخطي

رسم متحرك يوضح أول 83 تكرارًا لخوارزمية التدرج الهبوطي المطبقة على هذا المثال. الأسطح هي أسطح متساوية القيمة لـو(x(ن)){\displaystyle f(\mathbf {x} ^{(n)})}حسب التقديرات الحاليةx(ن){\displaystyle \mathbf {x} ^{(n)}}وتشير الأسهم إلى اتجاه الهبوط. ونظرًا لصغر حجم الخطوة وثباته، يكون التقارب بطيئًا.
{3x1-كوس(x2x3)-32=04x12-625x22+2x2-1=0خبرة(-x1x2)+20x3+10π-33=0{\displaystyle {\begin{cases}3x_{1}-\cos(x_{2}x_{3})-{\tfrac {3}{2}}=0\\4x_{1}^{2}-625x_{2}^{2}+2x_{2}-1=0\\\exp(-x_{1}x_{2})+20x_{3}+{\tfrac {10\pi -3}{3}}=0\end{cases}}}

لنقدم الدالة المرتبطة

جي(x)=[3x1-كوس(x2x3)-324x12-625x22+2x2-1خبرة(-x1x2)+20x3+10π-33]،{\displaystyle G(\mathbf {x} )={\begin{bmatrix}3x_{1}-\cos(x_{2}x_{3})-{\tfrac {3}{2}}\\4x_{1}^{2}-625x_{2}^{2}+2x_{2}-1\\\exp(-x_{1}x_{2})+20x_{3}+{\tfrac {10\pi -3}{3}}\\\end{bmatrix}},}

أين

x=[x1x2x3].{\displaystyle \mathbf {x} ={\begin{bmatrix}x_{1}\\x_{2}\\x_{3}\\\end{bmatrix}}.}

يمكن الآن تعريف دالة الهدف

و(x)=12جي(x)جي(x)=12[(3x1-كوس(x2x3)-32)2+(4x12-625x22+2x2-1)2+(خبرة(-x1x2)+20x3+10π-33)2]،{\displaystyle {\begin{aligned}f(\mathbf {x} )&={\frac {1}{2}}G^{\top }(\mathbf {x} )G(\mathbf {x} )\\&={\frac {1}{2}}\left[\left(3x_{1}-\cos(x_{2}x_{3})-{\frac {3}{2}}\right)^{2}+\left(4x_{1}^{2}-625x_{2}^{2}+2x_{2}-1\right)^{2}+\right.\\&{}\qquad \left.\left(\exp(-x_{1}x_{2})+20x_{3}+{\frac {10\pi -3}{3}}\right)^{2}\right],\end{aligned}}}

والتي سنحاول تقليلها. كتخمين أولي، دعونا نستخدم

x(0)=0=[000].{\displaystyle \mathbf {x} ^{(0)}=\mathbf {0} ={\begin{bmatrix}0\\0\\0\\\end{bmatrix}}.}

نحن نعلم ذلك

x(1)=0-η0و(0)=0-η0ججي(0)جي(0)،{\displaystyle \mathbf {x} ^{(1)}=\mathbf {0} -\eta _{0}\nabla f(\mathbf {0} )=\mathbf {0} -\eta _{0}J_{G}(\mathbf {0} )^{\top }G(\mathbf {0} ),}

حيث مصفوفة جاكوبيججي{\displaystyle J_{G}}يُعطى بواسطة

ججي(x)=[3الخطيئة(x2x3)x3الخطيئة(x2x3)x28x1-1250x2+20-x2خبرة(-x1x2)-x1خبرة(-x1x2)20].{\displaystyle J_{G}(\mathbf {x} )={\begin{bmatrix}3&\sin(x_{2}x_{3})x_{3}&\sin(x_{2}x_{3})x_{2}\\8x_{1}&-1250x_{2}+2&0\\-x_{2}\exp {(-x_{1}x_{2})}&-x_{1}\exp(-x_{1}x_{2})&20\\\end{bmatrix}}.}

نحسب:

ججي(0)=[3000200020]،جي(0)=[-2.5-110.472].{\displaystyle J_{G}(\mathbf {0} )={\begin{bmatrix}3&0&0\\0&2&0\\0&0&20\end{bmatrix}},\qquad G(\mathbf {0} )={\begin{bmatrix}-2.5\\-1\\10.472\end{bmatrix}}.}

هكذا

x(1)=0-η0[-7.5-2209.44]،{\displaystyle \mathbf {x} ^{(1)}=\mathbf {0} -\eta _{0}{\begin{bmatrix}-7.5\\-2\\209.44\end{bmatrix}},}

و

و(0)=0.5((-2.5)2+(-1)2+(10.472)2)=58.456.{\displaystyle f(\mathbf {0} )=0.5\left((-2.5)^{2}+(-1)^{2}+(10.472)^{2}\right)=58.456.}

والآن، مناسبη0{\displaystyle \eta _{0}}يجب إيجادها بحيث

و(x(1))و(x(0))=و(0).{\displaystyle f\left(\mathbf {x} ^{(1)}\right)\leq f\left(\mathbf {x} ^{(0)}\right)=f(\mathbf {0} ).}

يمكن القيام بذلك باستخدام أي من خوارزميات البحث الخطي المتنوعة . وقد يقوم المرء أيضاً بالتخمين ببساطة.η0=0.001،{\displaystyle \eta _{0}=0.001,}مما يعطي

x(1)=[0.00750.002-0.20944].{\displaystyle \mathbf {x} ^{(1)}={\begin{bmatrix}0.0075\\0.002\\-0.20944\\\end{bmatrix}}.}

بتقييم دالة الهدف عند هذه القيمة، ينتج

و(x(1))=0.5((-2.48)2+(-1.00)2+(6.28)2)=23.306.{\displaystyle f\left(\mathbf {x} ^{(1)}\right)=0.5\left((-2.48)^{2}+(-1.00)^{2}+(6.28)^{2}\right)=23.306.}

الانخفاض منو(0)=58.456{\displaystyle f(\mathbf {0} )=58.456}إلى قيمة الخطوة التالية لـ

و(x(1))=23.306{\displaystyle f\left(\mathbf {x} ^{(1)}\right)=23.306}

يمثل ذلك انخفاضًا كبيرًا في دالة الهدف. وستؤدي الخطوات اللاحقة إلى تقليل قيمتها أكثر حتى يتم التوصل إلى حل تقريبي للنظام.

تعليقات

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

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

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

يمكن تعديل خوارزمية التدرج الهبوطي باستخدام الزخم [ 16 ] ( نيستيروف ، بولياك [ 17 ] ، وفرانك-وولف [ 18 ] ) ومعاملات الكرة الثقيلة (المتوسطات المتحركة الأسية [ 19 ] والزخم الموجب-السالب [ 20 ] ). ومن أبرز الأمثلة على هذه الخوارزميات المُحسِّنة: آدم، وديفغراد، ويوجي، وآدا بيليف، وغيرها.

قد تكون الطرق القائمة على طريقة نيوتن وعكس مصفوفة هيسيان باستخدام تقنيات التدرج المترافق بدائل أفضل. [ 21 ] [ 22 ] عمومًا، تتقارب هذه الطرق في عدد أقل من التكرارات، ولكن تكلفة كل تكرار أعلى. ومن الأمثلة على ذلك طريقة BFGS التي تتضمن حساب مصفوفة في كل خطوة يُضرب بها متجه التدرج للانتقال إلى اتجاه "أفضل"، بالإضافة إلى خوارزمية بحث خطي أكثر تطورًا ، لإيجاد القيمة "المثلى".η.{\displaystyle \eta .}بالنسبة للمشاكل الكبيرة للغاية، حيث تهيمن مشاكل ذاكرة الكمبيوتر، يجب استخدام طريقة ذات ذاكرة محدودة مثل L-BFGS بدلاً من BFGS أو طريقة الانحدار الأسرع.

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

يمكن اعتبار خوارزمية التدرج الهبوطي بمثابة تطبيق لطريقة أويلر لحل المعادلات التفاضلية العادية.x(ت)=-و(x(ت)){\displaystyle x'(t)=-\nabla f(x(t))}إلى تدفق متدرج . وبدورها، يمكن اشتقاق هذه المعادلة كوحدة تحكم مثلى [ 23 ] لنظام التحكمx(ت)=u(ت){\displaystyle x'(t)=u(t)}معu(ت){\displaystyle u(t)}تم تقديمها في نموذج التقييمu(ت)=-و(x(ت)){\displaystyle u(t)=-\nabla f(x(t))}.

التعديلات

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

طرق التدرج السريع

اقترح يوري نيستيروف [ 24 ] تعديلًا بسيطًا يُتيح تقاربًا أسرع للمسائل المحدبة، وقد تم تعميمه لاحقًا. بالنسبة للمسائل الملساء غير المقيدة، تُسمى هذه الطريقة طريقة التدرج السريع (FGM) أو طريقة التدرج المُسرّع (AGM). تحديدًا، إذا كانت الدالة قابلة للتفاضلو{\displaystyle f}محدب وو{\displaystyle \nabla f}هي ليبشيتز ، ولا يُفترض ذلكو{\displaystyle f}إذا كانت الدالة محدبة بشدة ، فإن الخطأ في قيمة الهدف الناتج في كل خطوةك{\displaystyle k}سيتم تحديد قيمة باستخدام طريقة التدرج الهبوطي بواسطةيا(ك-1){\textstyle {\mathcal {O}}\left({k^{-1}}\right)}باستخدام تقنية تسريع نيستروف، ينخفض ​​الخطأ عنديا(ك-2){\textstyle {\mathcal {O}}\left({k^{-2}}\right)}[ 25 ] [ 26 ] من المعروف أن المعدليا(ك-2){\displaystyle {\mathcal {O}}\left({k^{-2}}\right)}يُعدّ خفض دالة التكلفة الأمثل لطرق التحسين من الدرجة الأولى. ومع ذلك، توجد إمكانية لتحسين الخوارزمية بتقليل العامل الثابت. تُقلّل طريقة التدرج المُحسَّن (OGM) [ 27 ] هذا العامل الثابت إلى النصف، وهي طريقة مثلى من الدرجة الأولى للمسائل واسعة النطاق. [ 28 ]

بالنسبة للمسائل المقيدة أو غير الملساء، يُطلق على طريقة التدرج التقريبي السريع (FPGM) الخاصة بـ Nesterov اسم طريقة التدرج التقريبي السريع (FPGM)، وهي تسريع لطريقة التدرج التقريبي .

طريقة الزخم أو طريقة الكرة الثقيلة

في محاولة لكسر نمط التذبذب في خوارزمية التدرج الهبوطي، تستخدم طريقة الزخم أو طريقة الكرة الثقيلة حد الزخم قياسًا على كرة ثقيلة تنزلق على سطح قيم الدالة المراد تصغيرها، [ 6 ] أو على حركة كتلة في ديناميكا نيوتن عبر وسط لزج في مجال قوة محافظ . [ 29 ] تحتفظ خوارزمية التدرج الهبوطي مع الزخم بتحديث الحل في كل تكرار، وتحدد التحديث التالي كمزيج خطي من التدرج والتحديث السابق. بالنسبة للتصغير التربيعي غير المقيد، يكون حد معدل التقارب النظري لطريقة الكرة الثقيلة متقاربًا تقريبًا مع حد معدل التقارب لطريقة التدرج المترافق الأمثل . [ 6 ]

تُستخدم هذه التقنية في خوارزمية التدرج العشوائي، وكامتداد لخوارزميات الانتشار العكسي المستخدمة لتدريب الشبكات العصبية الاصطناعية . [ 30 ] [ 31 ] في اتجاه التحديث، تُضيف خوارزمية التدرج العشوائي خاصية عشوائية. ويمكن استخدام الأوزان لحساب المشتقات.

الإضافات

يمكن توسيع خوارزمية التدرج الهبوطي لتشمل التعامل مع القيود عن طريق تضمين إسقاط على مجموعة القيود. لا تكون هذه الطريقة مجدية إلا عندما يكون حساب الإسقاط فعالاً على الحاسوب. في ظل افتراضات مناسبة، تتقارب هذه الطريقة. تُعد هذه الطريقة حالة خاصة من خوارزمية التقديم-التراجع للتضمينات الرتيبة (والتي تشمل البرمجة المحدبة والمتباينات التغيرية ). [ 32 ]

يُعدّ التدرج الهبوطي حالة خاصة من التدرج الهبوطي المرآوي باستخدام المسافة الإقليدية المربعة كمعامل تباعد بريغمان المعطى . [ 33 ]

الخصائص النظرية

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

أمثلة

انظر أيضاً

مراجع

  1. 1 2 بويد، ستيفن؛ فاندنبيرغ، ليفين (2004-03-08). التحسين المحدب . مطبعة جامعة كامبريدج. doi : 10.1017/cbo9780511804441 . ISBN 978-0-521-83378-3.
  2. ليمارشال، سي. (1 يناير 2012). "كوشي وطريقة التدرج" (ملف PDF) . في: غروتشل، م. (محرر). قصص التحسين . سلسلة دوكيومنتا ماثيماتيكا. المجلد 6 ( الطبعة الأولى). دار نشر EMS. الصفحات 251-254 . doi : 10.4171/dms/6/27 . ISBN    978-3-936609-58-5أُرشف من النسخة الأصلية (PDF) بتاريخ 29-12-2018 . تم الاطلاع عليه بتاريخ 26-01-2020 .
  3. ^ هادامارد، جاك (1908). "ذاكرة حول مشكلة التحليل المتعلقة بتوازن اللوحات المرنة المغلفة". مذكرات مقدمة من علماء متنوعين من أكاديمية العلوم بالمعهد الفرنسي . 33 .
  4. كوران، ر. (1943). "الأساليب التباينية لحل مسائل التوازن والاهتزازات" . نشرة الجمعية الرياضية الأمريكية . 49 (1): 1-23 . doi : 10.1090/S0002-9904-1943-07818-4 .
  5. كاري، هاسكل ب. (1944). "طريقة الانحدار الأسرع لمسائل التصغير غير الخطية" . مجلة الرياضيات التطبيقية الفصلية 2 (3): 258-261 . doi : 10.1090/qam/10667 .
  6. 1 2 3 4 5 بولياك، بوريس (1987). مقدمة في التحسين .
  7. 1 2 أكيلوف، جي بي؛ كانتوروفيتش، إل في (1982). التحليل الوظيفي (الطبعة الثانية ). مطبعة بيرغامون. رقم ISBN  0-08-023036-9.
  8. بارزيلاي، جوناثان؛ بورواين، جوناثان م. (1988). "طرق التدرج ذات حجم الخطوة بنقطتين". مجلة IMA للتحليل العددي . 8 (1): 141-148 . doi : 10.1093/imanum/8.1.141 .
  9. فليتشر، ر. (2005). "حول طريقة بارزيلاي-بورواين". في: تشي، ل.؛ تيو، ك.؛ يانغ، ش. (محررون). التحسين والتحكم مع التطبيقات . التحسين التطبيقي. المجلد 96. بوسطن: سبرينغر. الصفحات 235-256 . ISBN   0-387-24254-6.
  10. وولف، فيليب (أبريل 1969). "شروط التقارب لطرق الصعود". مجلة SIAM . 11 (2): 226-235 . doi : 10.1137/1011036 .
  11. بيرنشتاين، جيريمي؛ فهدات، أراش؛ يو، ييسونغ؛ ليو، مينغ يو (2020-06-12). "حول المسافة بين شبكتين عصبيتين واستقرار التعلم". arXiv : 2002.03432 [ cs.LG ].
  12. هايكين، سيمون س. نظرية المرشح التكيفي. بيرسون للتعليم الهند، 2008. - ص 108-142، 217-242
  13. سعد، يوسف (2003). الطرق التكرارية للأنظمة الخطية المتفرقة ( الطبعة الثانية ). فيلادلفيا، بنسلفانيا: جمعية الرياضيات الصناعية والتطبيقية. ص 195. ISBN   978-0-89871-534-7.
  14. 1 2 بووميستر، هنريكوس؛ دوغيرتي، أندرو؛ كنيازيف، أندرو ف. (2015). "التكييف المسبق غير المتماثل لطرق التدرج المترافق والانحدار الأشد" . وقائع علوم الحاسوب . 51 : 276-285 . arXiv : 1212.6680 . doi : 10.1016/j.procs.2015.05.241 .
  15. هولمز، م. (2023). مقدمة في الحوسبة العلمية وتحليل البيانات، الطبعة الثانية . سبرينغر. ISBN 978-3-031-22429-4.
  16. عبد القادر، روسلان؛ لياخوف، بافيل؛ ناغورنوف، نيكولاي (يناير 2023). "دراسة استقصائية لخوارزميات التحسين في الشبكات العصبية الحديثة" . الرياضيات . 11 (11): 2466. doi : 10.3390/math11112466 . ISSN 2227-7390 . 
  17. دياكونيكولاس، يلينا؛ جوردان، مايكل آي. (يناير 2021). "طرق الزخم المعممة: منظور هاميلتوني" . مجلة SIAM للتحسين . 31 (1): 915-944 . arXiv : 1906.00436 . doi : 10.1137/20M1322716 . ISSN 1052-6234 . 
  18. Meyer, Gerard G. L. (November 1974). "Accelerated Frank–Wolfe Algorithms". SIAM Journal on Control. 12 (4): 655–663. doi:10.1137/0312050. ISSN 0036-1402.
  19. Kingma, Diederik P.; Ba, Jimmy (2017-01-29), Adam: A Method for Stochastic Optimization, arXiv:1412.6980
  20. Xie, Zeke; Yuan, Li; Zhu, Zhanxing; Sugiyama, Masashi (2021-07-01). "Positive-Negative Momentum: Manipulating Stochastic Gradient Noise to Improve Generalization". Proceedings of the 38th International Conference on Machine Learning. PMLR: 11448–11458. arXiv:2103.17182.
  21. Press, W. H.; Teukolsky, S. A.; Vetterling, W. T.; Flannery, B. P. (1992). Numerical Recipes in C: The Art of Scientific Computing (2nd ed.). New York: Cambridge University Press. ISBN 0-521-43108-5.
  22. Strutz, T. (2016). Data Fitting and Uncertainty: A Practical Introduction to Weighted Least Squares and Beyond (2nd ed.). Springer Vieweg. ISBN 978-3-658-11455-8.
  23. Ross, I.M. (July 2019). "An optimal control theory for nonlinear optimization". Journal of Computational and Applied Mathematics. 354: 39–51. doi:10.1016/j.cam.2018.12.044. S2CID 127649426.
  24. Nesterov, Yurii (2004). Introductory Lectures on Convex Optimization: A Basic Course. Springer. ISBN 1-4020-7553-7.
  25. Vandenberghe, Lieven (2019). "Fast Gradient Methods"(PDF). Lecture notes for EE236C at UCLA.
  26. Walkington, Noel J. (2023). "Nesterov's Method for Convex Optimization". SIAM Review. 65 (2): 539–562. doi:10.1137/21M1390037. ISSN 0036-1445.
  27. كيم ، د.؛ فيسلر، ج. أ. (2016). "طرق مُحسَّنة من الدرجة الأولى للتقليل المحدب السلس" . البرمجة الرياضية . 151 ( 1-2 ): 81-107 . arXiv : 1406.5468 . doi : 10.1007/s10107-015-0949-3 . PMC 5067109. PMID 27765996. S2CID 207055414 .   
  28. دروري، يوئيل (2017). "التعقيد الدقيق القائم على المعلومات للتقليل المحدب السلس". مجلة التعقيد . 39 : 1-16 . arXiv : 1606.01424 . doi : 10.1016/j.jco.2016.11.001 . S2CID 205861966 . 
  29. تشيان، نينغ (يناير 1999). " حول مصطلح الزخم في خوارزميات تعلم التدرج الهبوطي". الشبكات العصبية . 12 (1): 145-151 . CiteSeerX 10.1.1.57.5612 . doi : 10.1016/S0893-6080(98)00116-6 . PMID 12662723. S2CID 2783597 .   
  30. "الزخم وتكيف معدل التعلم" . جامعة ويلاميت . تم الاطلاع عليه بتاريخ 17 أكتوبر 2014 .
  31. جيفري هينتون ؛ نيتيش سريفاستافا؛ كيفن سويرسكي. "طريقة الزخم" . كورسيرا . تم الاطلاع عليه بتاريخ 2 أكتوبر 2018 .جزء من سلسلة محاضرات لدورة Coursera عبر الإنترنت بعنوان " الشبكات العصبية للتعلم الآلي" مؤرشفة بتاريخ 31-12-2016 في Wayback Machine .
  32. كومبيتس، ب. ل.؛ بيسكيه، ج.-س. (2011). "طرق التقسيم التقريبي في معالجة الإشارات". في: بوشكه، هـ. هـ.؛ بوراشيك، ر . س.؛ كومبيتس، ب. ل.؛ إلسر، ف.؛ لوك، د. ر.؛ وولكوفيتش، هـ. (محررون). خوارزميات النقطة الثابتة للمسائل العكسية في العلوم والهندسة . نيويورك: سبرينغر. ص 185-212 . arXiv : 0912.3522 . ISBN  978-1-4419-9568-1.
  33. "خوارزمية الهبوط المرآوي" .
  34. 1 2 بوبيك، سيباستيان (2015). "التحسين المحدب: الخوارزميات والتعقيد". arXiv : 1405.4980 [ math.OC ].

للمزيد من القراءة

  • بويد، ستيفن ؛ فاندنبيرغ، ليفين (2004). "التقليل غير المقيد" (ملف PDF) . التحسين المحدب . نيويورك: مطبعة جامعة كامبريدج. الصفحات 457-520 . ISBN  0-521-83378-7.
  • تشونغ، إدوين كيه بي؛ زاك، ستانيسلاف إتش. (2013). "طرق التدرج" . مقدمة في التحسين (  الطبعة الرابعة). هوبوكين: وايلي. ص 131-160 . ISBN  978-1-118-27901-4.
  • هيملبلاو، ديفيد م. (1972). "إجراءات التصغير غير المقيدة باستخدام المشتقات". البرمجة غير الخطية التطبيقية . نيويورك: ماكجرو هيل. ص 63-132 . ISBN  0-07-028921-2.