مضاعف لاغرانج

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

ملخص وتبرير

تتمثل الفكرة الأساسية في تحويل مسألة مقيدة إلى صيغة تسمح بتطبيق اختبار المشتقة الخاص بمسألة غير مقيدة. وتؤدي العلاقة بين تدرج الدالة وتدرجات القيود بشكل طبيعي إلى إعادة صياغة المسألة الأصلية، والمعروفة بدالة لاغرانج أو دالة لاغرانج. [ 2 ] في الحالة العامة، تُعرَّف دالة لاغرانج على النحو التالي:

ل(x،λ)و(x)+λ،ز(x){\displaystyle {\mathcal {L}}(x,\lambda )\equiv f(x)+\langle \lambda ,g(x)\rangle }

للوظائفو،ز{\displaystyle f,g}؛ التدوين،{\displaystyle \langle \cdot ,\cdot \rangle }يشير إلى الضرب الداخلي . القيمةλ{\displaystyle \lambda }يُطلق عليه اسم مُضاعِف لاغرانج .

في الحالات البسيطة، حيث يُعرَّف الضرب الداخلي بأنه الضرب النقطي ، يكون لاغرانجيان

ل(x،λ)و(x)+λز(x){\displaystyle {\mathcal {L}}(x,\lambda )\equiv f(x)+\lambda \cdot g(x)}

يمكن تلخيص الطريقة على النحو التالي: لإيجاد القيمة العظمى أو الصغرى لدالة ماو{\displaystyle f}مع مراعاة شرط المساواةز(x)=0{\displaystyle g(x)=0}أوجد النقاط الثابتة لـل{\displaystyle {\mathcal {L}}}باعتبارها دالة لـx{\displaystyle x}ومضاعف لاغرانجλ {\displaystyle \lambda ~}وهذا يعني أن جميع المشتقات الجزئية يجب أن تساوي صفرًا، بما في ذلك المشتقة الجزئية بالنسبة إلىλ {\displaystyle \lambda ~}[ 3 ]

لx=0{\displaystyle {\frac {\partial {\mathcal {L}}}{\partial x}}=0} و  ل λ=0 ؛{\displaystyle {\frac {\ \partial {\mathcal {L}}\ }{\partial \lambda }}=0\ ;}

أو ما يعادل ذلك

و(x)x+λز(x)x=0{\displaystyle {\frac {\partial f(x)}{\partial x}}+\lambda \cdot {\frac {\partial g(x)}{\partial x}}=0} و ز(x)=0 .{\displaystyle g(x)=0~.}

إن الحل المقابل للتحسين المقيد الأصلي هو دائمًا نقطة سرجية لدالة لاغرانج، [ 4 ] [ 5 ] والتي يمكن تحديدها من بين النقاط الثابتة من خلال تحديد مصفوفة هيسيان ذات الحدود . [ 6 ]

تتمثل الميزة الكبرى لهذه الطريقة في أنها تسمح بحل مسائل التحسين دون الحاجة إلى تحديد المعاملات بشكل صريح بدلالة القيود. ونتيجة لذلك، تُستخدم طريقة مُضاعِفات لاغرانج على نطاق واسع لحل مسائل التحسين المقيدة المعقدة. علاوة على ذلك، تُعمَّم طريقة مُضاعِفات لاغرانج باستخدام شروط كاروش-كون-تاكر ، التي يمكنها أيضًا مراعاة قيود المتباينات من الشكل التالي:ح(x)ج{\displaystyle h(\mathbf {x} )\leq c}لثابت معينج{\displaystyle c}.

إفادة

ما يلي يُعرف باسم نظرية مضاعف لاغرانج. [ 7 ]

يتركو:RنR{\displaystyle f\colon \mathbb {R} ^{n}\to \mathbb {R} }لتكن دالة الهدف ولتكنز:RنRج{\displaystyle g\colon \mathbb {R} ^{n}\to \mathbb {R} ^{c}}لتكن دالة القيود، وكلاهما ينتمي إلىج1{\displaystyle C^{1}}(أي أن تكون مشتقاتها الأولى متصلة). لننظر في مسألة التحسين المقيد التالية:

أقصى و(x)رهناً بما يلي: ز(x)=0{\displaystyle {\begin{aligned}&{\text{maximize }}f(x)\\&{\text{subject to: }}g(x)=0\end{aligned}}}

يتركx{\displaystyle x_{\star }}ليكن حلاً أمثل لمسألة التحسين المذكورة أعلاه بحيث يكون، بالنسبة لمصفوفة المشتقات الجزئية[دز(x)]ج،ك= زج xك{\displaystyle \left[\operatorname {D} g(x_{\star })\right]_{j,k}={\frac {\ \partial g_{j}\ }{\partial x_{k}}}}،رتبة(دز(x))=جن{\displaystyle \operatorname {rank} (\operatorname {D} g(x_{\star }))=c\leq n}إذن، يوجد مُضاعِف لاغرانج فريدλRج{\displaystyle \lambda _{\star }\in \mathbb {R} ^{c}}بحيثدو(x)=λتيدز(x) .{\displaystyle \operatorname {D} f(x_{\star })=\lambda _{\star }^{\mathsf {T}}\operatorname {D} g(x_{\star })~.}(في هذه المعادلة،λ{\displaystyle \lambda _{\star }}هو متجه عمودي، لذا فإن منقولهλتي{\displaystyle \lambda _{\star }^{\mathsf {T}}}هو متجه صف. بدلاً من ذلك، يمكننا إعادة تعريف مُضاعِف لاغرانج مباشرةً كمتجه صف، وبالتالي تجنب عملية النقل.

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

قيد واحد

الشكل 1: يوضح المنحنى الأحمر القيد g ( x , y ) = c . أما المنحنيات الزرقاء فهي خطوط كفاف f ( x , y ) . النقطة التي يلامس فيها القيد الأحمر خط كفاف أزرق بشكل مماس هي القيمة القصوى لـ f ( x , y ) على طول القيد، لأن d1 > d2 .

في حالة وجود قيد واحد فقط ومتغيرين للاختيار فقط (كما هو موضح في الشكل 1)، ضع في اعتبارك مسألة التحسينأقصىx،yو(x،y)رهناً بـز(x،y)=0.{\displaystyle {\begin{aligned}{\underset {x,y}{\text{maximize}}}\quad &f(x,y)\\{\text{subject to}}\quad &g(x,y)=0.\end{aligned}}} (أحيانًا يتم عرض الثابت الإضافي بشكل منفصل بدلاً من تضمينه فيز{\displaystyle g}وفي هذه الحالة، يُكتب القيدز(x،y)=ج،{\displaystyle g(x,y)=c,}كما هو موضح في الشكل 1.) نفترض أن كليهماو{\displaystyle f}وز{\displaystyle g}لدينا مشتقات جزئية أولى متصلة . نقدم متغيرًا جديدًا (λ{\displaystyle \lambda }) يُسمى مُضاعِف لاغرانج (أو مُضاعِف لاغرانج غير المُحدَّد ) ودراسة دالة لاغرانج (أو لاغرانجيان أو تعبير لاغرانجيان ) المُعرَّفة بواسطة ل(x،y،λ)=و(x،y)+λز(x،y)،{\displaystyle {\mathcal {L}}(x,y,\lambda )=f(x,y)+\lambda \cdot g(x,y),} حيثλ{\displaystyle \lambda }يمكن إضافة المصطلح أو طرحه. إذاو(x0،y0){\displaystyle f(x_{0},y_{0})}الحد الأقصى هوو(x،y){\displaystyle f(x,y)}بالنسبة للمسألة المقيدة الأصلية وز(x0،y0)0،{\displaystyle \nabla g(x_{0},y_{0})\neq 0,}ثم يوجدλ0{\displaystyle \lambda _{0}}بحيث (x0،y0،λ0{\displaystyle x_{0},y_{0},\lambda _{0}}) هي نقطة ثابتة لدالة لاغرانج (النقاط الثابتة هي تلك النقاط التي تكون عندها المشتقات الجزئية الأولى لـل{\displaystyle {\mathcal {L}}}(صفر). الافتراضز0{\displaystyle \nabla g\neq 0}يُطلق على هذه العملية اسم "تأهيل القيد". مع ذلك، لا تُؤدي جميع النقاط الثابتة إلى حل للمسألة الأصلية، إذ تُقدّم طريقة مُضاعفات لاغرانج شرطًا ضروريًا فقط للأمثلية في المسائل المُقيدة. [ 9 ] [ 10 ] [ 11 ] [ 12 ] [ 13 ] توجد أيضًا شروط كافية للحد الأدنى أو الأقصى ، ولكن إذا استوفى حل مُرشّح مُعيّن الشروط الكافية، فإنه يُضمن فقط أن هذا الحل هو الأفضل محليًا - أي أنه أفضل من أي نقاط مجاورة مسموح بها. يُمكن إيجاد الحل الأمثل العالمي بمقارنة قيم دالة الهدف الأصلية عند النقاط التي تُحقق الشروط الضرورية والشروط الكافية محليًا.

تعتمد طريقة مُضاعِفات لاغرانج على الحدس القائل بأنه عند القيمة العظمى، لا يمكن أن تكون الدالة f ( x , y ) متزايدة في اتجاه أي نقطة مجاورة لها أيضًا g = 0. فلو كانت كذلك، لأمكننا التحرك على طول g = 0 للوصول إلى قيمة أعلى، مما يعني أن نقطة البداية لم تكن في الواقع القيمة العظمى. من هذا المنظور، تُشابه هذه الطريقة تمامًا اختبار ما إذا كانت مشتقة دالة غير مقيدة تساوي صفرًا ، أي أننا نتحقق من أن المشتقة الاتجاهية تساوي صفرًا في أي اتجاه ذي صلة (ممكن).

يمكننا تصور خطوط الكفاف لـ f المعطاة بواسطة f ( x , y ) = d لقيم مختلفة من d ، وخطوط الكفاف لـ g المعطاة بواسطة g ( x , y ) = c .

لنفترض أننا نسير على طول خط الكفاف حيث g = c . نحن مهتمون بإيجاد النقاط التي لا تتغير فيها f تقريبًا أثناء سيرنا، لأن هذه النقاط قد تكون نقاطًا عظمى.

هناك طريقتان يمكن أن يحدث بهما هذا:

  1. يمكننا لمس خط كفاف الدالة f ، لأن f ، بحسب التعريف ، لا تتغير أثناء سيرنا على طول خطوط كفافها. وهذا يعني أن المماسات لخطوط كفاف الدالتين f و g متوازية هنا.
  2. لقد وصلنا إلى جزء "مستوي" من f ، مما يعني أن f لا تتغير في أي اتجاه.

للتحقق من الاحتمال الأول (حيث نلامس خطًا كفافيًا للدالة f )، لاحظ أنه بما أن تدرج الدالة عمودي على الخطوط الكفافية، فإن المماسات للخطوط الكفافية للدالتين f و g تكون متوازية إذا وفقط إذا كان تدرج الدالتين f و g متوازيًا. وبالتالي، نريد النقاط ( x , y ) حيث g ( x , y ) = c و x،yو=λx،yز،{\displaystyle \nabla _{x,y}f=\lambda \,\nabla _{x,y}g,} بالنسبة للبعضλ{\displaystyle \lambda }أين x،yو=(وx،وy)،x،yز=(زx،زy)\displaystyle \nabla _{x,y}f=\left(\frac {\partial f}{\partial x},\frac {\partial f}{\partial y}\right),\qquad \nabla _{x,y}g=\left(\frac {\partial g}{\partial x},\frac {\partial g}{\partial y}\right)} هي التدرجات المعنية. الثابتλ{\displaystyle \lambda }هذا ضروري لأنه على الرغم من أن متجهي التدرج متوازيان، إلا أن مقداريهما لا يكونان متساويين عمومًا. يُسمى هذا الثابت مُضاعِف لاغرانج. (في بعض الاصطلاحات)λ{\displaystyle \lambda }(يسبقها علامة ناقص).

لاحظ أن هذه الطريقة تحل أيضًا الاحتمال الثاني، وهو أن تكون الدالة f مستوية: إذا كانت f مستوية، فإن ميلها يساوي صفرًا، وبالتالي فإنλ=0{\displaystyle \lambda =0}هو حل بغض النظر عنx،yز{\displaystyle \nabla _{x,y}g}.

ولدمج هذه الشروط في معادلة واحدة، نقدم دالة مساعدة ل(x،y،λ)و(x،y)+λز(x،y)،{\displaystyle {\mathcal {L}}(x,y,\lambda )\equiv f(x,y)+\lambda \cdot g(x,y)\,,} وحل x،y،λل(x،y،λ)=0 .{\displaystyle \nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )=0~.}لاحظ أن هذا يعني حل ثلاث معادلات بثلاثة مجاهيل. هذه هي طريقة معاملات لاغرانج.

لاحظ أن λل(x،y،λ)=0 {\displaystyle \ \nabla _{\lambda }{\mathcal {L}}(x,y,\lambda )=0\ }يشير إلى ز(x،y)=0 ،{\displaystyle \ g(x,y)=0\ ,}كمشتق جزئي لـل{\displaystyle {\mathcal {L}}}بالنسبة إلىλ{\displaystyle \lambda }يكون ز(x،y) .{\displaystyle \ g(x,y)~.}

باختصار x،y،λل(x،y،λ)=0{x،yو(x،y)=-λx،yز(x،y)ز(x،y)=0{\displaystyle \nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )=0\iff {\begin{cases}\nabla _{x,y}f(x,y)=-\lambda \,\nabla _{x,y}g(x,y)\\g(x,y)=0\end{cases}}}يمكن تعميم هذه الطريقة بسهولة على الدوال الموجودة علىن{\displaystyle n}المتغيرات x1،...،xن،λل(x1،...،xن،λ)=0{\displaystyle \nabla _{x_{1},\dots ,x_{n},\lambda }{\mathcal {L}}(x_{1},\dots ,x_{n},\lambda )=0} وهو ما يعادل حل n + 1 معادلة في n + 1 مجهول.

تمثل القيم القصوى المقيدة للدالة f نقاطًا حرجة في دالة لاغرانج.ل{\displaystyle {\mathcal {L}}}لكنها ليست بالضرورة قيمًا قصوى محلية لـل{\displaystyle {\mathcal {L}}}(انظر المثال  2 أدناه).

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

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

قيود متعددة

الشكل 2: قطع مكافئ مقيد على طول خطين متقاطعين.
الشكل 3: خريطة كونتورية للشكل 2.

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

على وجه التحديد، لنفترض أن لدينام{\displaystyle M}القيود ويسيرون على طول مجموعة النقاط التي تحققزأنا(x)=0،أنا=1،...،م.{\displaystyle g_{i}(\mathbf {x} )=0,i=1,\dots ,M\,.}كل نقطةx{\displaystyle \mathbf {x} }على محيط دالة قيد معينةزأنا{\displaystyle g_{i}}يحتوي على فضاء من الاتجاهات المسموح بها: فضاء المتجهات العمودية علىزأنا(x).{\displaystyle \nabla g_{i}(\mathbf {x} )\,.}وبالتالي، فإن مجموعة الاتجاهات المسموح بها وفقًا لجميع القيود هي فضاء الاتجاهات العمودية على تدرجات جميع القيود. نرمز إلى فضاء الحركات المسموح بها هذا بـ أ {\displaystyle \ A\ }ونرمز إلى مدى تدرجات القيود بـS.{\displaystyle S\,.}ثمأ=S،{\displaystyle A=S^{\perp }\,,}فضاء المتجهات العمودية على كل عنصر منS.{\displaystyle S\,.}

ما زلنا مهتمين بإيجاد نقاط حيثو{\displaystyle f}لا يتغير أثناء سيرنا، لأن هذه النقاط قد تكون نقاطًا قصوى (مقيدة). لذلك نسعى إلىx{\displaystyle \mathbf {x} }بحيث يكون أي اتجاه مسموح به للحركة بعيدًا عنx{\displaystyle \mathbf {x} }عمودي علىو(x){\displaystyle \nabla f(\mathbf {x} )}(وإلا فقد نزيدو{\displaystyle f}عن طريق التحرك في ذلك الاتجاه المسموح به). بعبارة أخرى،و(x)أ=S.{\displaystyle \nabla f(\mathbf {x} )\in A^{\perp }=S\,.}وبالتالي توجد كميات قياسيةλ1،λ2، ...،λم{\displaystyle \lambda _{1},\lambda _{2},\ \dots ,\lambda _{M}}بحيث و(x)=ك=1مλكزك(x)و(x)-ك=1مλكزك(x)=0 .{\displaystyle \nabla f(\mathbf {x} )=\sum _{k=1}^{M}\lambda _{k}\,\nabla g_{k}(\mathbf {x} )\quad \iff \quad \nabla f(\mathbf {x} )-\sum _{k=1}^{M}{\lambda _{k}\nabla g_{k}(\mathbf {x} )}=0~.}

هذه الكميات القياسية هي معاملات لاغرانج. لدينا الآنم{\displaystyle M}واحد منهم، واحد لكل قيد.

كما في السابق، نقدم دالة مساعدة ل(x1،...،xن،λ1،...،λم)=و(x1،...،xن)-ك=1مλكزك(x1،...،xن) {\displaystyle {\mathcal {L}}\left(x_{1},\ldots ,x_{n},\lambda _{1},\ldots ,\lambda _{M}\right)=f\left(x_{1},\ldots ,x_{n}\right)-\sum \limits _{k=1}^{M}{\lambda _{k}g_{k}\left(x_{1},\ldots ,x_{n}\right)}\ } وحل x1،...،xن،λ1،...،λمل(x1،...،xن،λ1،...،λم)=0{و(x)-ك=1مλكزك(x)=0ز1(x)==زم(x)=0{\displaystyle \nabla _{x_{1},\ldots ,x_{n},\lambda _{1},\ldots ,\lambda _{M}}{\mathcal {L}}(x_{1},\ldots ,x_{n},\lambda _{1},\ldots ,\lambda _{M})=0\iff {\begin{cases}\nabla f(\mathbf {x} )-\sum _{k=1}^{M}{\lambda _{k}\,\nabla g_{k}(\mathbf {x} )}=0\\g_{1}(\mathbf {x} )=\cdots =g_{M}(\mathbf {x} )=0\end{cases}}} وهو ما يعادل حلن+م{\displaystyle n+M}المعادلات في ن+م {\displaystyle \ n+M\ }مجهول.

الافتراض المتعلق بتأهيل القيود عند وجود قيود متعددة هو أن تدرجات القيود عند النقطة ذات الصلة مستقلة خطيًا.

صياغة حديثة عبر مشعبات قابلة للتفاضل

يمكن تعميم مشكلة إيجاد القيم العظمى والصغرى المحلية الخاضعة للقيود لتشمل إيجاد القيم العظمى والصغرى المحلية على مشعب قابل للتفاضل م .{\displaystyle \ M~.}[ 14 ] فيما يلي، ليس من الضروري أنم{\displaystyle M}قد يكون فضاءً إقليديًا، أو حتى متعدد الشعب الريماني . جميع مظاهر التدرج  {\displaystyle \ \nabla \ }(الذي يعتمد على اختيار المقياس الريماني) يمكن استبداله بالمشتق الخارجي د{\displaystyle \ \operatorname {d} }.

قيد واحد

يترك م {\displaystyle \ M\ }ليكن متعدد الأبعاد أملسًا م .{\displaystyle \ m~.}لنفترض أننا نرغب في إيجاد النقاط الثابتة x {\displaystyle \ x\ }دالة سلسة و:مR {\displaystyle \ f:M\to \mathbb {R} \ }عندما يقتصر على الشعبة الفرعية شمال {\displaystyle \ N\ }محدد بواسطة ز(x)=0 ،{\displaystyle \ g(x)=0\ ,}أين ز:مR {\displaystyle \ g:M\to \mathbb {R} \ }هي دالة سلسة يكون فيها الصفر قيمة منتظمة .

يترك دو {\displaystyle \ \operatorname {d} f\ }و دز {\displaystyle \ \operatorname {d} g\ }تكون المشتقات الخارجية لـ و {\displaystyle \ f\ }و ز {\displaystyle \ g\ }. استقرار القيد و|شمال {\displaystyle \ f|_{N}\ }في xشمال {\displaystyle \ x\in N\ }وسائل د(و|شمال)x=0 .{\displaystyle \ \operatorname {d} (f|_{N})_{x}=0~.}وبعبارة أخرى، النواة كير(دوx) {\displaystyle \ \ker(\operatorname {d} f_{x})\ }يتضمن تيxشمال=كير(دزx) .{\displaystyle \ T_{x}N=\ker(\operatorname {d} g_{x})~.}بعبارة أخرى، دوx {\displaystyle \ \operatorname {d} f_{x}\ }و دزx {\displaystyle \ \operatorname {d} g_{x}\ }هي أشكال تناسبية من الدرجة الأولى. ولذلك، من الضروري والكافي أن يكون النظام التالي من 12م(م-1) {\displaystyle \ {\tfrac {1}{2}}m(m-1)\ }المعادلة صحيحة: دوxدزx=0Λ2(تيx*م){\displaystyle \operatorname {d} f_{x}\wedge \operatorname {d} g_{x}=0\in \Lambda ^{2}(T_{x}^{\ast }M)} أين  {\displaystyle \ \wedge \ }يشير إلى المنتج الخارجي . النقاط الثابتة x {\displaystyle \ x\ }هي حلول نظام المعادلات المذكور أعلاه بالإضافة إلى القيد ز(x)=0 .{\displaystyle \ g(x)=0~.}لاحظ أن 12م(م-1) {\displaystyle \ {\tfrac {1}{2}}m(m-1)\ }المعادلات ليست مستقلة، لأن الطرف الأيسر من المعادلة ينتمي إلى مجموعة فرعية من Λ2(تيx*م) {\displaystyle \ \Lambda ^{2}(T_{x}^{\ast }M)\ }يتكون من عناصر قابلة للتحلل .

في هذه الصيغة، ليس من الضروري إيجاد مُضاعِف لاغرانج صراحةً، وهو عدد λ {\displaystyle \ \lambda \ }بحيث دوx=λدزx .{\displaystyle \ \operatorname {d} f_{x}=\lambda \cdot \operatorname {d} g_{x}~.}

قيود متعددة

يترك م {\displaystyle \ M\ }و و {\displaystyle \ f\ }يكون الأمر كما هو موضح في القسم أعلاه فيما يتعلق بحالة القيد الواحد. بدلاً من الدالةز{\displaystyle g}كما هو موضح هناك، لنفترض الآن دالة سلسة جي:مRص(ص>1) ،{\displaystyle \ G:M\to \mathbb {R} ^{p}(p>1)\ ,}مع وظائف المكونات زأنا:مR ،{\displaystyle \ g_{i}:M\to \mathbb {R} \ ,}والتي0Rص{\displaystyle 0\in \mathbb {R} ^{p}}هي قيمة منتظمة . ليكنشمال{\displaystyle N}ليكن التشعب الفرعي لـ م {\displaystyle \ M\ }محدد بواسطة جي(x)=0 .{\displaystyle \ G(x)=0~.}

 x {\displaystyle \ x\ }هي نقطة ثابتة لـو|شمال{\displaystyle f|_{N}}إذا وفقط إذا كير(دوx) {\displaystyle \ \ker(\operatorname {d} f_{x})\ }يتضمن كير(دجيx) .{\displaystyle \ \ker(\operatorname {d} G_{x})~.}للتسهيل، دعنا لx=دوx {\displaystyle \ L_{x}=\operatorname {d} f_{x}\ }و كx=دجيx ،{\displaystyle \ K_{x}=\operatorname {d} G_{x}\ ,}أين دجي{\displaystyle \ \operatorname {d} G}يشير إلى خريطة المماس أو جاكوبيان تيمتيRص {\displaystyle \ TM\to T\mathbb {R} ^{p}~}( تيxRص{\displaystyle \ T_{x}\mathbb {R} ^{p}}يمكن تحديدها بشكل أساسي بـ Rص{\displaystyle \ \mathbb {R} ^{p}}). الفضاء الجزئيكير(كx){\displaystyle \ker(K_{x})}أبعاده أصغر من أبعادكير(لx){\displaystyle \ker(L_{x})}، أي خافت(كير(لx))=ن-1 {\displaystyle \ \dim(\ker(L_{x}))=n-1\ }و خافت(كير(كx))=ن-ص .{\displaystyle \ \dim(\ker(K_{x}))=n-p~.}كير(كx){\displaystyle \ker(K_{x})}ينتمي إلى كير(لx) {\displaystyle \ \ker(L_{x})\ }إذا وفقط إذالxتيx*م{\displaystyle L_{x}\in T_{x}^{\ast }M}ينتمي إلى صورة كx*:Rص*تيx*م .{\displaystyle \ K_{x}^{\ast }:\mathbb {R} ^{p\ast }\to T_{x}^{\ast }M~.}من الناحية الحسابية، الشرط هو أنلx{\displaystyle L_{x}}ينتمي إلى فضاء الصفوف لمصفوفة كx ،{\displaystyle \ K_{x}\ ,}أو ما يعادله، فضاء أعمدة مصفوفةكx*{\displaystyle K_{x}^{\ast }}(المنقول). إذا ωxΛص(تيx*م) {\displaystyle \ \omega _{x}\in \Lambda ^{p}(T_{x}^{\ast }M)\ }يرمز إلى الضرب الخارجي لأعمدة مصفوفة كx* ،{\displaystyle \ K_{x}^{\ast }\ ,}الشرط الثابت لـ و|شمال {\displaystyle \ f|_{N}\ }في x {\displaystyle \ x\ }يصبح لxωx=0Λص+1(تيx*م){\displaystyle L_{x}\wedge \omega _{x}=0\in \Lambda ^{p+1}\left(T_{x}^{\ast }M\right)} ومرة أخرى، في هذه الصيغة ليس من الضروري إيجاد معاملات لاغرانج بشكل صريح، الأرقام λ1،...،λص {\displaystyle \ \lambda _{1},\ldots ,\lambda _{p}\ }بحيث  دوx=أنا=1صλأناد(زأنا)x .{\displaystyle \ \operatorname {d} f_{x}=\sum _{i=1}^{p}\lambda _{i}\operatorname {d} (g_{i})_{x}~.}

تفسير معاملات لاغرانج

في هذا القسم، نقوم بتعديل معادلات القيد من الشكل التالي:زأنا(x)=0{\displaystyle g_{i}({\bf {x}})=0}إلى النموذج زأنا(x)=جأنا ،{\displaystyle \ g_{i}({\bf {x}})=c_{i}\ ,}حيث جأنا {\displaystyle \ c_{i}\ }هي m ثوابت حقيقية تُعتبر وسائط إضافية للتعبير اللاغرانجيل{\displaystyle {\mathcal {L}}}.

غالبًا ما تُفسَّر معاملات لاغرانج على أنها كمية معينة ذات أهمية. على سبيل المثال، من خلال تحديد خط الكفاف للقيد، أي إذا كان تعبير لاغرانج هو ل(x1،x2،...؛λ1،λ2،...؛ج1،ج2،...)=و(x1،x2،...)+λ1(ج1-ز1(x1،x2،...))+λ2(ج2-ز2(x1،x2،...))+{\displaystyle {\begin{aligned}&{\mathcal {L}}(x_{1},x_{2},\ldots ;\lambda _{1},\lambda _{2},\ldots ;c_{1},c_{2},\ldots )\\[4pt]={}&f(x_{1},x_{2},\ldots )+\lambda _{1}(c_{1}-g_{1}(x_{1},x_{2},\ldots ))+\lambda _{2}(c_{2}-g_{2}(x_{1},x_{2},\dots ))+\cdots \end{aligned}}} ثم  لجك=λك .{\displaystyle \ {\frac {\partial {\mathcal {L}}}{\partial c_{k}}}=\lambda _{k}~.}

إذن، λk هو معدل تغير الكمية المراد تحسينها كدالة لمعامل القيد. على سبيل المثال، في ميكانيكا لاغرانج، تُشتق معادلات الحركة بإيجاد النقاط الثابتة للفعل ، وهو التكامل الزمني للفرق بين الطاقة الحركية والطاقة الكامنة. وبالتالي، يمكن تفسير القوة المؤثرة على جسيم نتيجة جهد قياسي، F = −∇V ، على أنها مُضاعِف لاغرانج يُحدد التغير في الفعل (انتقال الطاقة الكامنة إلى الطاقة الحركية) تبعًا لتغير في مسار الجسيم المُقيد. في نظرية التحكم، يُصاغ هذا بدلاً من ذلك كمعادلات الحالة المرافقة .

علاوة على ذلك، وبحسب نظرية الغلاف، فإن القيمة المثلى لمضاعف لاغرانج تُفسَّر على أنها التأثير الهامشي لثابت القيد المقابل على القيمة المثلى الممكنة لدالة الهدف الأصلية: إذا رمزنا للقيم عند القيمة المثلى بعلامة النجمة ({\displaystyle \star })، عندئذٍ يمكن إثبات أن  دو( x1(ج1،ج2،...)، x2(ج1،ج2،...)، ... ) دجك=λك .{\displaystyle {\frac {\ \operatorname {d} f\left(\ x_{1\star }(c_{1},c_{2},\dots ),\ x_{2\star }(c_{1},c_{2},\dots ),\ \dots \ \right)\ }{\operatorname {d} c_{k}}}=\lambda _{\star k}~.}

على سبيل المثال، في علم الاقتصاد، يتم حساب الربح الأمثل للاعب في ظل مجموعة محدودة من الإجراءات، حيث يمثل مُضاعِف لاغرانج التغير في القيمة المثلى لدالة الهدف (الربح) نتيجة لتخفيف قيد معين (مثلاً من خلال تغيير في الدخل)؛ في مثل هذا السياق λك {\displaystyle \ \lambda _{\star k}\ }يمثل التكلفة الحدية للقيد، ويشار إليه باسم السعر الظلي . [ 15 ]

شروط كافية

يمكن تحديد الشروط الكافية لوجود قيمة عظمى أو صغرى محلية مقيدة بدلالة سلسلة من المحددات الرئيسية (محددات المصفوفات الفرعية المحاذية للأعلى الأيسر) لمصفوفة هيسيان ذات الحدود للمشتقات الثانية للتعبير اللاغرانجي. [ 6 ] [ 16 ]

أمثلة

المثال 1

توضيح لمسألة التحسين المقيد 1 

لنفترض أننا نرغب في تحقيق أقصى قدر من و(x،y)=x+y {\displaystyle \ f(x,y)=x+y\ }رهناً بالقيود x2+y2=1 .{\displaystyle \ x^{2}+y^{2}=1~.}مجموعة الحلول الممكنة هي دائرة الوحدة، ومجموعات مستوى الدالة f هي خطوط قطرية (بميل -1)، لذا يمكننا أن نرى بيانيًا أن القيمة القصوى تحدث عند (12،12) ،{\displaystyle \ \left({\tfrac {1}{\sqrt {2}}},{\tfrac {1}{\sqrt {2}}}\right)\ ,}وأن الحد الأدنى يحدث عند (-12،-12) .{\displaystyle \ \left(-{\tfrac {1}{\sqrt {2}}},-{\tfrac {1}{\sqrt {2}}}\right)~.}

بالنسبة لطريقة معاملات لاغرانج، فإن القيد هو ز(x،y)=x2+y2-1=0 ،{\displaystyle g(x,y)=x^{2}+y^{2}-1=0\ ,} ومن ثم دالة لاغرانج، ل(x،y،λ)=و(x،y)+λز(x،y)=x+y+λ(x2+y2-1) ،{\displaystyle {\begin{aligned}{\mathcal {L}}(x,y,\lambda )&=f(x,y)+\lambda \cdot g(x,y)\\[4pt]&=x+y+\lambda (x^{2}+y^{2}-1)\ ,\end{aligned}}} هي دالة مكافئة لـ و(x،y) {\displaystyle \ f(x,y)\ }متى ز(x،y) {\displaystyle \ g(x,y)\ }تم ضبطه على 0 .

والآن يمكننا حساب التدرج: x،y،λل(x،y،λ)=(لx،لy،لλ)=(1+2λx،1+2λy،x2+y2-1) ،{\displaystyle {\begin{aligned}\nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )&=\left({\frac {\partial {\mathcal {L}}}{\partial x}},{\frac {\partial {\mathcal {L}}}{\partial y}},{\frac {\partial {\mathcal {L}}}{\partial \lambda }}\right)\\[4pt]&=\left(1+2\lambda x,1+2\lambda y,x^{2}+y^{2}-1\right)\ \color {gray}{,}\end{aligned}}} وبالتالي: x،y،λل(x،y،λ)=0{1+2λx=01+2λy=0x2+y2-1=0{\displaystyle \nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )=0\quad \Leftrightarrow \quad {\begin{cases}1+2\lambda x=0\\1+2\lambda y=0\\x^{2}+y^{2}-1=0\end{cases}}}

لاحظ أن المعادلة الأخيرة هي القيد الأصلي.

المعادلتان الأوليان تعطيان x=y=-12λ،λ0 .{\displaystyle x=y=-{\frac {1}{2\lambda }},\qquad \lambda \neq 0~.} بالتعويض في المعادلة الأخيرة نحصل على: 14λ2+14λ2-1=0 ،{\displaystyle {\frac {1}{4\lambda ^{2}}}+{\frac {1}{4\lambda ^{2}}}-1=0\ ,} لذا λ=±12  ،{\displaystyle \lambda =\pm {\frac {1}{\sqrt {2\ }}}\ ,} مما يعني أن النقاط الثابتة لـل{\displaystyle {\mathcal {L}}}نكون (2 2،2 2،-12 )،(-2 2،-2 2،12 ) .{\displaystyle \left({\tfrac {\sqrt {2\ }}{2}},{\tfrac {\sqrt {2\ }}{2}},-{\tfrac {1}{\sqrt {2\ }}}\right),\qquad \left(-{\tfrac {\sqrt {2\ }}{2}},-{\tfrac {\sqrt {2\ }}{2}},{\tfrac {1}{\sqrt {2\ }}}\right)~.}

تقييم دالة الهدف f عند هذه النقاط ينتج عنه و(2 2،2 2)=2  ،و(-2 2،-2 2)=-2  .{\displaystyle f\left({\tfrac {\sqrt {2\ }}{2}},{\tfrac {\sqrt {2\ }}{2}}\right)={\sqrt {2\ }}\ ,\qquad f\left(-{\tfrac {\sqrt {2\ }}{2}},-{\tfrac {\sqrt {2\ }}{2}}\right)=-{\sqrt {2\ }}~.}

وبالتالي فإن الحد الأقصى المقيد هو 2  {\displaystyle \ {\sqrt {2\ }}\ }والحد الأدنى المقيد هو-2{\displaystyle -{\sqrt {2}}}.

المثال 2

توضيح لمسألة التحسين المقيد 2 

الآن نعدل دالة الهدف في المثال 1 بحيث نخفض  و(x،y)=(x+y)2 {\displaystyle \ f(x,y)=(x+y)^{2}\ }بدلاً من و(x،y)=x+y ،{\displaystyle \ f(x,y)=x+y\ ,}مرة أخرى على طول الدائرة ز(x،y)=x2+y2-1=0 .{\displaystyle \ g(x,y)=x^{2}+y^{2}-1=0~.}الآن مجموعات المستويات منو{\displaystyle f}لا تزال هذه الخطوط ذات ميل -1، والنقاط على الدائرة المماسية لهذه المستويات هي مرة أخرى (2/2،2/2) {\displaystyle \ ({\sqrt {2}}/2,{\sqrt {2}}/2)\ }و (-2/2،-2/2) .{\displaystyle \ (-{\sqrt {2}}/2,-{\sqrt {2}}/2)~.}تمثل نقاط التماس هذه القيم القصوى لـ و .{\displaystyle \ f~.}

من ناحية أخرى، تحدث القيم الدنيا على مجموعة المستويات لـ و=0 {\displaystyle \ f=0\ }(منذ بنائها) و {\displaystyle \ f\ }لا يمكن أن تأخذ قيمًا سالبة عند (2/2،-2/2) {\displaystyle \ ({\sqrt {2}}/2,-{\sqrt {2}}/2)\ }و (-2/2،2/2) ،{\displaystyle \ (-{\sqrt {2}}/2,{\sqrt {2}}/2)\ ,}حيث منحنيات المستوى لـ و {\displaystyle \ f\ }لا تكون مماسية للقيد. الشرط هو أن x،y،λ(و(x،y)+λز(x،y))=0 {\displaystyle \ \nabla _{x,y,\lambda }\left(f(x,y)+\lambda \cdot g(x,y)\right)=0\ }يُحدد البرنامج جميع النقاط الأربع بشكل صحيح على أنها نقاط قصوى؛ وتتميز النقاط الدنيا بـ λ=0 {\displaystyle \ \lambda =0\ }والقيم القصوى بواسطة λ=-2 .{\displaystyle \ \lambda =-2~.}

المثال 3

توضيح لمسألة التحسين المقيد 3 . 

يتناول هذا المثال حسابات أكثر تعقيدًا، ولكنه لا يزال مشكلة ذات قيد واحد.

لنفترض أن المرء يريد إيجاد القيم القصوى لـ و(x،y)=x2y{\displaystyle f(x,y)=x^{2}y} بشرط أن x {\displaystyle \ x\ }- و y {\displaystyle \ y\ }تقع الإحداثيات على الدائرة المحيطة بنقطة الأصل بنصف قطر 3  .{\displaystyle \ {\sqrt {3\ }}~.}أي، رهناً بالقيد ز(x،y)=x2+y2-3=0 .{\displaystyle g(x,y)=x^{2}+y^{2}-3=0~.}

بما أن هناك قيدًا واحدًا فقط، فهناك مُضاعِف واحد، على سبيل المثال λ .{\displaystyle \ \lambda ~.}

القيد ز(x،y) {\displaystyle \ g(x,y)\ }يساوي الصفر تمامًا على الدائرة التي نصف قطرها 3  .{\displaystyle \ {\sqrt {3\ }}~.}أي مضاعف لـ ز(x،y) {\displaystyle \ g(x,y)\ }يمكن إضافته إلى ز(x،y) {\displaystyle \ g(x,y)\ }مغادرة ز(x،y) {\displaystyle \ g(x,y)\ }دون تغيير في المنطقة محل الاهتمام (على الدائرة حيث يتم استيفاء القيد الأصلي).

يؤدي تطبيق طريقة مضاعف لاغرانج العادية إلى ل(x،y،λ)=و(x،y)+λز(x،y)=x2y+λ(x2+y2-3) ،{\displaystyle {\begin{aligned}{\mathcal {L}}(x,y,\lambda )&=f(x,y)+\lambda \cdot g(x,y)\\&=x^{2}y+\lambda (x^{2}+y^{2}-3)\ ,\end{aligned}}} ومنها يمكن حساب التدرج: x،y،λل(x،y،λ)=(لx،لy،لλ)=(2xy+2λx،x2+2λy،x2+y2-3) .{\displaystyle {\begin{aligned}\nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )&=\left({\frac {\partial {\mathcal {L}}}{\partial x}},{\frac {\partial {\mathcal {L}}}{\partial y}},{\frac {\partial {\mathcal {L}}}{\partial \lambda }}\right)\\&=\left(2xy+2\lambda x,x^{2}+2\lambda y,x^{2}+y^{2}-3\right)~.\end{aligned}}} وبالتالي: x،y،λل(x،y،λ)=0{2xy+2λx=0x2+2λy=0x2+y2-3=0{x(y+λ)=0(أنا)x2=-2λy(ii)x2+y2=3(ثالثاً){\displaystyle \nabla _{x,y,\lambda }{\mathcal {L}}(x,y,\lambda )=0\quad \iff \quad {\begin{cases}2xy+2\lambda x=0\\x^{2}+2\lambda y=0\\x^{2}+y^{2}-3=0\end{cases}}\quad \iff \quad {\begin{cases}x(y+\lambda )=0&{\text{(i)}}\\x^{2}=-2\lambda y&{\text{(ii)}}\\x^{2}+y^{2}=3&{\text{(iii)}}\end{cases}}} (iii) هو القيد الأصلي فقط. (i) يستلزم x=0 {\displaystyle \ x=0\ }أو λ=-y .{\displaystyle \ \lambda =-y~.}لوx=0{\displaystyle x=0}ثم y=±3  {\displaystyle \ y=\pm {\sqrt {3\ }}\ }بواسطة (ثالثاً) وبالتالي  λ=0 {\displaystyle \ \lambda =0\ }من (ii). إذا λ=-y ،{\displaystyle \ \lambda =-y\ ,}وباستبدال هذا في (ii) ينتج x2=2y2 .{\displaystyle \ x^{2}=2y^{2}~.}بتعويض هذا في (iii) وحل المعادلة لإيجاد y {\displaystyle \ y\ }أعطِ y=±1 .{\displaystyle \ y=\pm 1~.}وبالتالي، توجد ست نقاط حرجة ل :{\displaystyle \ {\mathcal {L}}\ :} (2 ،1،-1)؛(-2 ،1،-1)؛(2 ،-1،1)؛(-2 ،-1،1)؛(0،3 ،0)؛(0،-3 ،0) .{\displaystyle ({\sqrt {2\ }},1,-1);\quad (-{\sqrt {2\ }},1,-1);\quad ({\sqrt {2\ }},-1,1);\quad (-{\sqrt {2\ }},-1,1);\quad (0,{\sqrt {3\ }},0);\quad (0,-{\sqrt {3\ }},0)~.}

بتقييم الهدف عند هذه النقاط، يجد المرء أن و(±2 ،1)=2؛و(±2 ،-1)=-2؛و(0،±3 )=0 .{\displaystyle f(\pm {\sqrt {2\ }},1)=2;\quad f(\pm {\sqrt {2\ }},-1)=-2;\quad f(0,\pm {\sqrt {3\ }})=0~.}

وبالتالي، تصل دالة الهدف إلى القيمة القصوى العالمية (مع مراعاة القيود) عند (±2 ،1 ){\displaystyle \ (\pm {\sqrt {2\ }},1\ )}والحد الأدنى العالمي عند (±2 ،-1) .{\displaystyle \ (\pm {\sqrt {2\ }},-1)~.}النقطة (0،3 ) {\displaystyle \ (0,{\sqrt {3\ }})\ }هو الحد الأدنى المحلي لـ و {\displaystyle \ f\ }و (0،-3 ) {\displaystyle \ (0,-{\sqrt {3\ }})\ }هي قيمة عظمى محلية لـ و ،{\displaystyle \ f\ ,}كما يمكن تحديده من خلال النظر في مصفوفة هيسيان لـ ل(x،y،0) .{\displaystyle \ {\mathcal {L}}(x,y,0)~.}

لاحظ أنه بينما (2 ،1،-1) {\displaystyle \ ({\sqrt {2\ }},1,-1)\ }تُعد نقطة حاسمة من ل ،{\displaystyle \ {\mathcal {L}}\ ,}إنها ليست حالة قصوى محلية لـ ل .{\displaystyle \ {\mathcal {L}}~.}لدينا ل(2 +ε،1،-1+دلتا)=2+دلتا(ε2+(22 )ε) .{\displaystyle {\mathcal {L}}\left({\sqrt {2\ }}+\varepsilon ,1,-1+\delta \right)=2+\delta \left(\varepsilon ^{2}+\left(2{\sqrt {2\ }}\right)\varepsilon \right)~.}

بالنظر إلى أي حي من (2 ،1،-1) ،{\displaystyle \ ({\sqrt {2\ }},1,-1)\ ,}يمكن للمرء أن يختار جانبًا إيجابيًا صغيرًا ε {\displaystyle \ \varepsilon \ }وصغير دلتا {\displaystyle \ \delta \ }من أي من العلامتين للحصول على ل{\displaystyle \ {\mathcal {L}}}القيم الأكبر والأصغر من 2 .{\displaystyle \ 2~.}ويمكن ملاحظة ذلك أيضًا من خلال مصفوفة هيسيان الخاصة بـ ل {\displaystyle \ {\mathcal {L}}\ }يتم تقييمها عند هذه النقطة (أو في الواقع عند أي من النقاط الحرجة) وهي مصفوفة غير محددة . كل نقطة من النقاط الحرجة لـ ل {\displaystyle \ {\mathcal {L}}\ }هي نقطة سرجية ل .{\displaystyle \ {\mathcal {L}}~.}[ 4 ]

المثال 4 - الإنتروبيا

لنفترض أننا نرغب في إيجاد التوزيع الاحتمالي المنفصل على النقاط {ص1،ص2،...،صن} {\displaystyle \ \{p_{1},p_{2},\ldots ,p_{n}\}\ }مع أقصى إنتروبيا معلوماتية . وهذا يعني أننا نرغب في إيجاد أقل توزيع احتمالي منظم على النقاط {ص1،ص2،،صن} .{\displaystyle \ \{p_{1},p_{2},\cdots ,p_{n}\}~.}بمعنى آخر، نرغب في تعظيم معادلة إنتروبيا شانون : و(ص1،ص2،...،صن)=-ج=1نصجسجل2صج .{\displaystyle f(p_{1},p_{2},\ldots ,p_{n})=-\sum _{j=1}^{n}p_{j}\log _{2}p_{j}~.}

لكي يكون هذا توزيعًا احتماليًا، يجب أن يكون مجموع الاحتمالات صأنا {\displaystyle \ p_{i}\ }عند كل نقطة xأنا {\displaystyle \ x_{i}\ }يجب أن يساوي 1، لذا فإن قيدنا هو: ز(ص1،ص2،...،صن)=ج=1نصج=1 .{\displaystyle g(p_{1},p_{2},\ldots ,p_{n})=\sum _{j=1}^{n}p_{j}=1~.}

نستخدم معاملات لاغرانج لإيجاد نقطة أقصى إنتروبيا، ص* ،{\displaystyle \ {\vec {p}}^{\,*}\ ,}عبر جميع التوزيعات الاحتمالية المنفصلة ص {\displaystyle \ {\vec {p}}\ }على {x1،x2،...،xن} .{\displaystyle \ \{x_{1},x_{2},\ldots ,x_{n}\}~.}نحن نشترط ما يلي: ص(و+λ(ز-1))|ص=ص*=0 ،{\displaystyle \left.{\frac {\partial }{\partial {\vec {p}}}}(f+\lambda (g-1))\right|_{{\vec {p}}={\vec {p}}^{\,*}}=0\ ,} مما ينتج عنه نظام من n معادلة، ك=1، ...،ن ،{\displaystyle \ k=1,\ \ldots ,n\ ,}بحيث: صك{-(ج=1نصجسجل2صج)+λ(ج=1نصج-1)}|صك=صك=0 .{\displaystyle \left.{\frac {\partial }{\partial p_{k}}}\left\{-\left(\sum _{j=1}^{n}p_{j}\log _{2}p_{j}\right)+\lambda \left(\sum _{j=1}^{n}p_{j}-1\right)\right\}\right|_{p_{k}=p_{\star k}}=0~.}

بإجراء عملية التفاضل لهذه المعادلات n ، نحصل على -(1ln2+سجل2صك)+λ=0 .{\displaystyle -\left({\frac {1}{\ln 2}}+\log _{2}p_{\star k}\right)+\lambda =0~.}

هذا يدل على أن كل صك {\displaystyle \ p_{\star k}\ }متساويتان (لأنهما تعتمدان على λ فقط). باستخدام القيد جصج=1 ،{\displaystyle \sum _{j}p_{j}=1\ ,} نجد صك=1ن .{\displaystyle p_{\star k}={\frac {1}{n}}~.}

وبالتالي، فإن التوزيع المنتظم هو التوزيع ذو أكبر إنتروبيا، من بين التوزيعات على n نقطة.

المثال 5 - التحسين العددي

تتسبب مضاعفات لاغرانج في حدوث النقاط الحرجة عند نقاط السرج (المثال 5 ). 
يمكن استخدام مقدار التدرج لإجبار النقاط الحرجة على الحدوث عند الحد الأدنى المحلي (المثال 5 ). 

تقع النقاط الحرجة لدوال لاغرانج عند نقاط السرج ، وليس عند القيم العظمى (أو الصغرى) المحلية. [ 4 ] [ 17 ] لسوء الحظ، صُممت العديد من تقنيات التحسين العددي، مثل تسلق التلال ، وانحدار التدرج ، وبعض طرق شبه نيوتن ، وغيرها، لإيجاد القيم العظمى (أو الصغرى) المحلية وليس نقاط السرج. لهذا السبب، يجب إما تعديل الصيغة لضمان أنها مسألة تصغير (على سبيل المثال، عن طريق إيجاد القيمة القصوى لمربع تدرج دالة لاغرانج كما هو موضح أدناه)، أو استخدام تقنية تحسين تجد النقاط الثابتة (مثل طريقة نيوتن بدون بحث خطي عن القيم القصوى ) وليس بالضرورة القيم القصوى.

كمثال بسيط، لننظر في مشكلة إيجاد قيمة x التي تقلل من و(x)=x2 ،{\displaystyle \ f(x)=x^{2}\ ,}مقيدة بحيث x2=1 .{\displaystyle \ x^{2}=1~.}(هذه المشكلة غير نمطية إلى حد ما لأنه لا يوجد سوى قيمتين تحققان هذا القيد، لكنها مفيدة لأغراض التوضيح لأنه يمكن تصور الدالة غير المقيدة المقابلة في ثلاثة أبعاد.)

باستخدام معاملات لاغرانج، يمكن تحويل هذه المشكلة إلى مشكلة تحسين غير مقيدة: ل(x،λ)=x2+λ(x2-1) .{\displaystyle {\mathcal {L}}(x,\lambda )=x^{2}+\lambda (x^{2}-1)~.}

تحدث النقطتان الحرجتان عند نقاط السرج حيث x = 1 و x = −1 .

لحل هذه المشكلة باستخدام تقنية التحسين العددي، يجب أولاً تحويلها بحيث تقع النقاط الحرجة عند القيم الدنيا المحلية. ويتم ذلك بحساب مقدار تدرج مسألة التحسين غير المقيدة.

أولاً، نقوم بحساب المشتقة الجزئية للمسألة غير المقيدة بالنسبة لكل متغير: لx=2x+2xλلλ=x2-1 .{\displaystyle {\begin{aligned}&{\frac {\partial {\mathcal {L}}}{\partial x}}=2x+2x\lambda \\[5pt]&{\frac {\partial {\mathcal {L}}}{\partial \lambda }}=x^{2}-1~.\end{aligned}}}

إذا لم تكن الدالة المستهدفة قابلة للتفاضل بسهولة، فيمكن تقريب التفاضل بالنسبة لكل متغير على النحو التالي:  ل xل(x+ε،λ)-ل(x،λ)ε، ل λل(x،λ+ε)-ل(x،λ)ε،{\displaystyle {\begin{aligned}{\frac {\ \partial {\mathcal {L}}\ }{\partial x}}\approx {\frac {{\mathcal {L}}(x+\varepsilon ,\lambda )-{\mathcal {L}}(x,\lambda )}{\varepsilon }},\\[5pt]{\frac {\ \partial {\mathcal {L}}\ }{\partial \lambda }}\approx {\frac {{\mathcal {L}}(x,\lambda +\varepsilon )-{\mathcal {L}}(x,\lambda )}{\varepsilon }},\end{aligned}}} أينε{\displaystyle \varepsilon }قيمة صغيرة.

بعد ذلك، نحسب مقدار التدرج، وهو الجذر التربيعي لمجموع مربعات المشتقات الجزئية: ح(x،λ)=(2x+2xλ)2+(x2-1)2 ( ل(x+ε،λ)-ل(x،λ) ε)2+( ل(x،λ+ε)-ل(x،λ) ε)2  .{\displaystyle {\begin{aligned}h(x,\lambda )&={\sqrt {(2x+2x\lambda )^{2}+(x^{2}-1)^{2}\ }}\\[4pt]&\approx {\sqrt {\left({\frac {\ {\mathcal {L}}(x+\varepsilon ,\lambda )-{\mathcal {L}}(x,\lambda )\ }{\varepsilon }}\right)^{2}+\left({\frac {\ {\mathcal {L}}(x,\lambda +\varepsilon )-{\mathcal {L}}(x,\lambda )\ }{\varepsilon }}\right)^{2}\ }}~.\end{aligned}}}

(بما أن المقدار دائمًا غير سالب، فإن التحسين على أساس مربع المقدار يكافئ التحسين على أساس المقدار. وبالتالي، يمكن حذف "الجذر التربيعي" من هذه المعادلات دون أي فرق متوقع في نتائج التحسين.)

تحدث النقاط الحرجة للدالة h عند x = 1 و x = −1 ، تمامًا كما فيل .{\displaystyle {\mathcal {L}}~.}بخلاف النقاط الحرجة فيل،{\displaystyle {\mathcal {L}}\,,}ومع ذلك، فإن النقاط الحرجة في h تحدث عند الحد الأدنى المحلي، لذلك يمكن استخدام تقنيات التحسين العددي لإيجادها.

التطبيقات

ميكانيكا لاغرانج

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

عند استخدام معاملات لاغرانج، يجب حل معادلات القيد بالتزامن مع معادلات أويلر-لاغرانج. وبالتالي، تصبح المعادلات نظامًا من المعادلات التفاضلية الجبرية (بدلاً من نظام المعادلات التفاضلية العادية ). [ 19 ]

تُعدّ طريقة مُضاعِفات لاغرانج مفيدةً عندما يصعب كتابة دالة لاغرانج بدلالة مجموعة من الإحداثيات المعممة المستقلة خطيًا. على سبيل المثال، تُستخدم في خوارزميات نمذجة الأنظمة الديناميكية البرمجية، أو في نمذجة الأنظمة ذات السلاسل الحركية المغلقة. [ 20 ] كما أنها مفيدة لفرض قيود غير هولونومية. [ 18 ] [ 20 ]

بالنظر إلى مجموعة من معادلات القيود الهولونوميةوج(q،ت)=0{\displaystyle f_{j}(\mathbf {q} ,t)=0}ويمكن كتابة معادلات أويلر-لاغرانج مع مضاعفات لاغرانج على النحو التالي [ 18 ] [ 19 ]

ددتلq˙أنا-لqأنا+ج=1جλجوجqأنا-τأنا،قيد=τأنا{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\frac {\partial L}{\partial {\dot {q}}_{i}}}-{\frac {\partial L}{\partial q_{i}}}+\underbrace {\sum _{j=1}^{C}\lambda _{j}{\frac {\partial f_{j}}{\partial q_{i}}}} _{-\tau _{i,{\text{constraint}}}}=\tau _{i}}

معنىτأنا،قيد{\displaystyle \tau _{i,{\text{constraint}}}}يمكن تفسير ذلك بنقله إلى الجانب الآخر من المعادلة ودمجه في حد القوة المعممةτأنا{\displaystyle \tau _{i}}في هذا التفسير، يمتلك النظامج{\displaystyle C}عدد من درجات الحرية الإضافية، ولا توجد قيود مفروضة إضافية، ولكن قوى القيدτأنا،قيد{\displaystyle \tau _{i,{\text{constraint}}}}تصادف أن تكون لها القيم المناسبة بحيث تتحقق القيود. [ 18 ] [ 19 ]

نظرية التحكم

في نظرية التحكم الأمثل ، يتم تفسير مضاعفات لاغرانج على أنها متغيرات الحالة المرافقة ، ويتم إعادة صياغة مضاعفات لاغرانج على أنها تقليل الهاميلتوني ، في مبدأ بونترياجين الأقصى .

البرمجة غير الخطية

تتضمن طريقة مُضاعِف لاغرانج عدة تعميمات. في البرمجة غير الخطية، توجد عدة قواعد للمُضاعِف، مثل قاعدة كاراثيودوري-جون للمُضاعِف وقاعدة المُضاعِف المحدب، لقيود المتباينات. [ 21 ]

الاقتصاد

في العديد من نماذج الاقتصاد الرياضي ، مثل نماذج التوازن العام ، يُعبَّر عن سلوك المستهلك بتعظيم المنفعة ، وعن سلوك الشركة بتعظيم الربح ، مع خضوع كلا الكيانين لقيود مثل قيود الميزانية وقيود الإنتاج . ويتم عادةً تحديد الحل الأمثل من خلال تعظيم دالة معينة، حيث تُفرض القيود باستخدام مُضاعِفات لاغرانج. [ 22 ] [ 23 ] [ 24 ] [ 25 ]

أنظمة الطاقة

تُستخدم الطرق القائمة على مضاعفات لاغرانج في أنظمة الطاقة ، على سبيل المثال في تحديد مواقع موارد الطاقة الموزعة (DER) وتخفيف الأحمال. [ 26 ]

التعلم الآمن المعزز

تُطبَّق طريقة مُضاعِفات لاغرانج على عمليات اتخاذ القرار ماركوف المقيدة . [ 27 ] وهي تُنتج بشكل طبيعي خوارزميات ثنائية أولية قائمة على التدرج في التعلم المعزز الآمن. [ 28 ]

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

انظر أيضاً

مراجع

  1. هوفمان، لورانس د.؛ برادلي، جيرالد ل. (2004). حساب التفاضل والتكامل للأعمال والاقتصاد والعلوم الاجتماعية والحياتية (  الطبعة الثامنة). ماكجرو هيل للتعليم العالي. الصفحات 575-588 . ISBN  0-07-242432-X.
  2. بيفيس، برايان؛ دوبس، إيان م. (1990). "التحسين الساكن" . نظرية التحسين والاستقرار للتحليل الاقتصادي . نيويورك: مطبعة جامعة كامبريدج. ص 40. ISBN  0-521-33605-8.
  3. بروتر، موراي هـموري، تشارلز ب. الابن (1985). حساب التفاضل والتكامل المتوسط ​​( الطبعة الثانية). نيويورك، نيويورك: سبرينغر. ص 267. ISBN   0-387-96058-9.
  4. 1 2 3 والش، جي آر (1975). "خاصية نقطة السرج لدالة لاغرانج" . طرق التحسين . نيويورك، نيويورك: جون وايلي وأولاده. ص 39-44 . ISBN  0-471-91922-5.
  5. كالمان، دان (2009). " التسوية باستخدام لاغرانج: منظور بديل للتحسين المقيد". مجلة الرياضيات . 82 (3): 186-196 . doi : 10.1080/0025570X.2009.11953617 . JSTOR 27765899. S2CID 121070192 .  
  6. 1 2 سيلبربرغ، يوجين؛ سوين، وينغ (2001). بنية الاقتصاد: تحليل رياضي ( الطبعة الثالثة). بوسطن: إيروين ماكجرو هيل. ص 134-141 . ISBN   0-07-234352-4.
  7. دي لا فوينتي ، أنجيل (2000). الأساليب والنماذج الرياضية للاقتصاديين . كامبريدج: مطبعة جامعة كامبريدج. ص 285. doi : 10.1017/CBO9780511810756 . ISBN  978-0-521-58512-5.
  8. لونبرغر، ديفيد ج. (1969). التحسين باستخدام طرق الفضاء المتجهي . نيويورك: جون وايلي وأولاده. ص 188-189 . 
  9. بيرتسيكاس، ديمتري ب. (1999). البرمجة غير الخطية ( الطبعة الثانية). كامبريدج، ماساتشوستس: أثينا ساينتيفيك. ISBN  1-886529-00-0.
  10. فابنيارسكي، آي بي (2001) [1994]، "مضاعفات لاغرانج" ، موسوعة الرياضيات ، دار نشر إي إم إس.
  11. لاسدون، ليون س. (2002) [1970]. نظرية التحسين للأنظمة الكبيرة (طبعة مُعاد طباعتها). مينولا، نيويورك، نيويورك: دوفر. ISBN   0-486-41999-1MR 1888251 . 
  12. ^ هيريارت أوروتي، جان بابتيست؛ ليمارشال، كلود (1993). “الفصل الثاني عشر: الازدواجية المجردة للممارسين”. التحليل المحدب وخوارزميات التقليل . Grundlehren der Mathematischen Wissenschaften [المبادئ الأساسية للعلوم الرياضية]. المجلد. 306. برلين، DE: Springer-Verlag. ص 136 – 193 (والتعليقات الببليوغرافية ص 334 – 335). رقم ISBN    3-540-56852-2MR 1295240. المجلد الثاني : النظرية المتقدمة وطرق الحزم.  
  13. ^ ليمارشال، كلود (15-19 مايو 2000). "استرخاء لاغرانج". في جونجر، مايكل؛ نادف، دينيس (محرران). التحسين التوافقي الحسابي: أوراق من مدرسة الربيع التي عقدت في شلوس داغستوهل . أقيمت مدرسة الربيع في شلوس داغستوهل في الفترة من 15 إلى 19 مايو 2000 . ملاحظات محاضرة في علوم الكمبيوتر. المجلد. 2241. برلين، DE: Springer-Verlag (نشرت عام 2001). ص 112 – 156. دوى : 10.1007 / 3-540-45586-8_4 . رقم ISBN   3-540-42877-1MR 1900016 . S2CID 9048698 .​  
  14. لافونتين، جاك (2015). مقدمة في المشعبات التفاضلية . سبرينغر. ص 70. ISBN  978-3-319-20735-3.
  15. ديكسيت، أفيناش ك. (1990). "أسعار الظل" . الأمثلية في النظرية الاقتصادية ( الطبعة الثانية). نيويورك: مطبعة جامعة أكسفورد. ص 40-54 . ISBN   0-19-877210-6.
  16. ↑ تشيانغ ، ألفا سي. (1984). الأساليب الأساسية للاقتصاد الرياضي ( الطبعة الثالثة). ماكجرو هيل. ص 386. ISBN   0-07-010813-7.
  17. هيث، مايكل ت. (2005). الحوسبة العلمية: دراسة تمهيدية . ماكجرو هيل. ص 203. ISBN  978-0-07-124489-3.
  18. 1 2 3 4 غولدشتاين، هربرت ؛ بول، تشارلز ب. الابن ؛ سافكو، جون ل. (2002). الميكانيكا الكلاسيكية ( الطبعة الثالثة). سان فرانسيسكو، كاليفورنيا: أديسون ويسلي. الصفحات 45-51 . ISBN   0-201-65702-3.
  19. 1 2 3 برينان، ك. إي.؛ كامبل، س. ل.؛ بيتزولد، ل. ر. (1995). الحل العددي لمسائل القيمة الابتدائية في المعادلات التفاضلية الجبرية . جمعية الرياضيات الصناعية والتطبيقية. ص 4-5 . doi : 10.1137/1.9781611971224 . 
  20. 1 2 روي، فيذرستون (2008). خوارزميات ديناميكيات الأجسام الصلبة . سبرينغر نيويورك. ص 42-45 ، 141-143 . doi : 10.1007/978-1-4899-7560-7 . 
  21. بوركيو، بروس هـ. (1980). "قواعد المضاعف الحديثة" . المجلة الرياضية الأمريكية الشهرية . 87 (6): 433-452 . doi : 10.2307/2320250 . JSTOR 2320250 . 
  22. كامين، إم آي ؛ شوارتز، إن إل (1991). التحسين الديناميكي: حساب التفاضل والتكامل والتحكم الأمثل في الاقتصاد والإدارة ( الطبعة الثانية). نيويورك: إلسيفير. ISBN  0-444-01609-0.
  23. غلوتزل، إرهارد؛ غلوتزل، فلورنتين؛ ريخترز، أوليفر (2019). "من التحسين المقيد إلى الديناميكيات المقيدة: توسيع أوجه التشابه بين الاقتصاد والميكانيكا". مجلة التفاعل والتنسيق الاقتصادي . 14 (3): 623-642 . doi : 10.1007/s11403-019-00252-7 . hdl : 10419/171974 .
  24. باكسلي، جون ف.؛ مورهاوس، جون س. (1984). "مسائل مضاعف لاغرانج في الاقتصاد". المجلة الرياضية الأمريكية الشهرية . 91 (7): 404-412 . doi : 10.1080/00029890.1984.11971446 ..
  25. جانوفا، جيتكا (2011). "تطبيقات منهجية الميكانيكا المقيدة في الاقتصاد". المجلة الأوروبية للفيزياء . 32 (6): 1443-1463 . arXiv : 1106.3455 . Bibcode : 2011EJPh...32.1443J . doi : 10.1088/0143-0807/32/6/001 .
  26. غوتام، موكيش؛ بهوسال، نارايان؛ بنيدريس، محمد (2020). منهج قائم على الحساسية لتخفيف الأحمال التكيفي عند انخفاض التردد . مؤتمر IEEE تكساس للطاقة والكهرباء 2020 (TPEC). معهد مهندسي الإلكترونيات والكهرباء . الصفحات 1-5 . doi : 10.1109/TPEC48276.2020.9042569 . 
  27. ألتمان، إيتان (2021). عمليات اتخاذ القرار ماركوف المقيدة . روتليدج .
  28. دينغ، دونغشنغ؛ تشانغ، كايتشينغ؛ يوفانوفيتش، ميهايلو؛ باسار، تامر (2020). طريقة التدرج الطبيعي للسياسة الأولية-الثنائية لعمليات اتخاذ القرار ماركوف المقيدة . التطورات في أنظمة معالجة المعلومات العصبية.

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

  • بيفيس، برايان؛ دوبس، إيان م. (1990). "التحسين الساكن" . نظرية التحسين والاستقرار للتحليل الاقتصادي . نيويورك، نيويورك: مطبعة جامعة كامبريدج. ص 32-72 . ISBN  0-521-33605-8.
  • بيرتسيكاس، ديمتري ب. (1982). التحسين المقيد وطرق مضاعف لاغرانج . نيويورك، نيويورك: أكاديميك برس. ISBN 0-12-093480-9.
  • بيفريدج، جوردون إس جي؛ شيشتر، روبرت إس. (1970). "المضاعفات اللاغرانجية" . التحسين: النظرية والتطبيق . نيويورك، نيويورك: ماكجرو هيل. ص 244-259 . ISBN  0-07-005128-3.
  • بينجر، برايان ر.؛ هوفمان، إليزابيث (1998). "التحسين المقيد". الاقتصاد الجزئي مع حساب التفاضل والتكامل (  الطبعة الثانية). ريدينغ: أديسون-ويسلي. ص 56-91 . ISBN  0-321-01225-9.
  • كارتر، مايكل (2001). "قيود المساواة" . أسس الاقتصاد الرياضي . كامبريدج، ماساتشوستس: مطبعة معهد ماساتشوستس للتكنولوجيا. ص 516-549 . ISBN  0-262-53192-5.
  • هيستينز، ماغنوس ر. (1966). "القيم الصغرى للدوال الخاضعة لقيود المساواة". حساب التفاضل والتكامل ونظرية التحكم الأمثل . نيويورك، نيويورك: وايلي. ص 29-34 . 
  • وايلي، سي. راي؛ باريت، لويس سي. (1995). "القيم القصوى للتكاملات تحت القيود". الرياضيات الهندسية المتقدمة (  الطبعة السادسة). نيويورك، نيويورك: ماكجرو هيل. الصفحات 1096-1103 . ISBN  0-07-072206-4.

معرض

نصوص إضافية وتطبيقات تفاعلية