हैमिल्टनियन (नियंत्रण सिद्धांत)

हैमिल्टनियन फलन (गणित) है जिसका उपयोग गतिशील प्रणाली के इष्टतम नियंत्रण की समस्या को हल करने के लिए किया जाता है। इसे उस समस्या के लैग्रेंज गुणक के तात्कालिक वृद्धि के रूप में समझा जा सकता है जिसे निश्चित समय अवधि में अनुकूलित किया जाना है। हेमिल्टनियन यांत्रिकी से प्रेरित, लेकिन उससे अलग, इष्टतम नियंत्रण सिद्धांत के हैमिल्टनियन को लेव पोंट्रीगिन ने अपने पोंट्रीगिन के न्यूनतम सिद्धांत के भागों के रूप में विकसित किया था। पोंट्रीगिन ने सिद्ध किया कि इष्टतम नियंत्रण समस्या को हल करने के लिए आवश्यक शर्त यह है कि हैमिल्टन को अनुकूलित करने के लिए नियंत्रण को चुना जाना चाहिए।

समस्या कथन और हैमिल्टनियन की परिभाषा

की गतिशील प्रणाली पर विचार करें $n$ प्रथम-क्रम अंतर समीकरण

${\dot {\mathbf {x} }}(t)=\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)$

जहाँ $\mathbf {x} (t)=\left[x_{1}(t),x_{2}(t),\ldots ,x_{n}(t)\right]^{\mathsf {T}}$ स्थिति चर के वेक्टर को दर्शाता है, और $\mathbf {u} (t)=\left[u_{1}(t),u_{2}(t),\ldots ,u_{r}(t)\right]^{\mathsf {T}}$ नियंत्रण चर का वेक्टर एक बार प्रारंभिक शर्तें $\mathbf {x} (t_{0})=\mathbf {x} _{0}$ और नियंत्रित करता है $\mathbf {u} (t)$ निर्दिष्ट हैं, अंतर समीकरणों का समाधान, जिसे प्रक्षेपवक्र कहा जाता है $\mathbf {x} (t;\mathbf {x} _{0},t_{0})$, पाया जा सकता है। इष्टतम नियंत्रण की समस्या चुनना है $\mathbf {u} (t)$ (कुछ सेट से ${\mathcal {U}}\subseteq \mathbb {R} ^{r}$) जिससे $\mathbf {x} (t)$ प्रारंभिक समय के बीच निश्चित हानि फलन को अधिकतम या कम करता है $t=t_{0}$ और टर्मिनल समय $t=t_{1}$ (जहाँ $t_{1}$ अनंत हो सकता है)। विशेष रूप से, लक्ष्य प्रदर्शन सूचकांक का अनुकूलन करना है $I(\mathbf {x} (t),\mathbf {u} (t),t)$ समय के प्रत्येक बिंदु पर,

$\max _{\mathbf {u} (t)}J=\int _{t_{0}}^{t_{1}}I[\mathbf {x} (t),\mathbf {u} (t),t]\,\mathrm {d} t$

स्थिति चर की गति के उपरोक्त समीकरणों के अधीन। समाधान विधि में सहायक कार्य को परिभाषित करना सम्मिलित है जिसे नियंत्रण हैमिल्टन के रूप में जाना जाता है

$H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)\equiv I(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)$

जो ऑब्जेक्टिव फलन और स्टेट समीकरण को स्टैटिक ऑप्टिमाइज़ेशन प्रॉब्लम में लैग्रेंज मल्टीप्लायर की तरह जोड़ता है, केवल यह कि मल्टीप्लायर $\mathbf {\lambda } (t)$, कॉस्टेट वेरिएबल्स के रूप में संदर्भित, स्थिरांक के अतिरिक्त समय के कार्य हैं।

लक्ष्य इष्टतम नियंत्रण नीति कार्य खोजना है $\mathbf {u} ^{\ast }(t)$ और, इसके साथ, स्थिति चर का इष्टतम प्रक्षेपवक्र $\mathbf {x} ^{\ast }(t)$, जो पोंट्रीगिन के अधिकतम सिद्धांत के अनुसार वे तर्क हैं जो हैमिल्टनियन को अधिकतम करते हैं,

$H(\mathbf {x} ^{\ast }(t),\mathbf {u} ^{\ast }(t),\mathbf {\lambda } (t),t)\geq H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)$ सभी के लिए $\mathbf {u} (t)\in {\mathcal {U}}$

अधिकतम के लिए प्रथम-क्रम आवश्यक शर्तें किसके द्वारा दी गई हैं

${\frac {\partial H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)}{\partial \mathbf {u} }}=0$ जो अधिकतम सिद्धांत है,

${\frac {\partial H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)}{\partial \mathbf {\lambda } }}={\dot {\mathbf {x} }}$ जो स्थिति संक्रमण फलन उत्पन्न करता है $\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)={\dot {\mathbf {x} }}$,

${\frac {\partial H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)}{\partial \mathbf {x} }}=-{\dot {\mathbf {\lambda } }}(t)$ जो उत्पन्न करता है ${\dot {\mathbf {\lambda } }}(t)=-\left[I_{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} _{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)\right]$

जिनमें से बाद वाले को कॉस्टेट समीकरण कहा जाता है। साथ में, स्थिति और कॉस्टेट समीकरण हैमिल्टनियन गतिशील प्रणाली का वर्णन करते हैं (भौतिक विज्ञान में हैमिल्टनियन प्रणाली से फिर से समान लेकिन अलग), जिसके समाधान में दो-बिंदु सीमा मूल्य समस्या सम्मिलित है, यह देखते हुए कि $2n$ समय में दो अलग-अलग बिंदुओं को सम्मिलित करने वाली सीमा की स्थिति, प्रारंभिक समय ( $n$ स्थिति चर के लिए अंतर समीकरण), और टर्मिनल समय ( $n$ कॉस्टेट वेरिएबल्स के लिए डिफरेंशियल समीकरण; जब तक कोई अंतिम कार्य निर्दिष्ट नहीं किया जाता है, सीमा की स्थिति होती है $\mathbf {\lambda } (t_{1})=0$, या $\lim _{t_{1}\to \infty }\mathbf {\lambda } (t_{1})=0$ अनंत समय क्षितिज के लिए)।

अधिकतम के लिए पर्याप्त शर्त हैमिल्टनियन की अवतलता है जिसका मूल्यांकन समाधान में किया गया है, अर्थात

$H_{\mathbf {uu} }(\mathbf {x} ^{\ast }(t),\mathbf {u} ^{\ast }(t),\mathbf {\lambda } (t),t)\leq 0$

जहाँ $\mathbf {u} ^{\ast }(t)$ इष्टतम नियंत्रण है, और $\mathbf {x} ^{\ast }(t)$ स्थिति चर के लिए इष्टतम प्रक्षेपवक्र का परिणाम है। वैकल्पिक रूप से,ओल्वी एल मंगसेरियन के परिणामस्वरूप, कार्य करने पर आवश्यक शर्तें पर्याप्त हैं $I(\mathbf {x} (t),\mathbf {u} (t),t)$ और $\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)$ दोनों अवतल हैं $\mathbf {x} (t)$ और $\mathbf {u} (t)$.

लाग्रंगियन से व्युत्पत्ति

विवश अनुकूलन समस्या जैसा कि ऊपर कहा गया है, विशेष रूप से लाग्रंगियन अभिव्यक्ति का सुझाव देती है

$L=\int _{t_{0}}^{t_{1}}I(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\left[\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)-{\dot {\mathbf {x} }}(t)\right]\,\mathrm {d} t$

जहाँ $\mathbf {\lambda } (t)$ स्थिर अनुकूलन समस्या में लैग्रेंज गुणक की तुलना करता है लेकिन अब, जैसा कि ऊपर उल्लेख किया गया है, समय का कार्य है। मिटाने के लिए ${\dot {\mathbf {x} }}(t)$, दाईं ओर के अंतिम शब्द को भागों द्वारा एकीकरण का उपयोग करके फिर से लिखा जा सकता है, जैसे कि

$-\int _{t_{0}}^{t_{1}}\mathbf {\lambda } ^{\mathsf {T}}(t){\dot {\mathbf {x} }}(t)\,\mathrm {d} t=-\mathbf {\lambda } ^{\mathsf {T}}(t_{1})\mathbf {x} (t_{1})+\mathbf {\lambda } ^{\mathsf {T}}(t_{0})\mathbf {x} (t_{0})+\int _{t_{0}}^{t_{1}}{\dot {\mathbf {\lambda } }}^{\mathsf {T}}(t)\mathbf {x} (t)\,\mathrm {d} t$

जिसे देने के लिए लाग्रंगियन अभिव्यक्ति में वापस प्रतिस्थापित किया जा सकता है

$L=\int _{t_{0}}^{t_{1}}\left[I(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)+{\dot {\mathbf {\lambda } }}^{\mathsf {T}}(t)\mathbf {x} (t)\right]\,\mathrm {d} t-\mathbf {\lambda } ^{\mathsf {T}}(t_{1})\mathbf {x} (t_{1})+\mathbf {\lambda } ^{\mathsf {T}}(t_{0})\mathbf {x} (t_{0})$

इष्टतम के लिए प्रथम-क्रम की स्थिति प्राप्त करने के लिए, मान लें कि समाधान मिल गया है और लाग्रंगियन अधिकतम हो गया है। फिर किसी तरह की गड़बड़ी $\mathbf {x} (t)$ या $\mathbf {u} (t)$ लाग्रंगियन के मूल्य में गिरावट का कारण होना चाहिए। विशेष रूप से, का कुल व्युत्पन्न $L$ का अनुसरण करता है

$\mathrm {d} L=\int _{t_{0}}^{t_{1}}\left[\left(I_{\mathbf {u} }(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} _{\mathbf {u} }(\mathbf {x} (t),\mathbf {u} (t),t)\right)\mathrm {d} \mathbf {u} (t)+\left(I_{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} _{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)+{\dot {\mathbf {\lambda } }}(t)\right)\mathrm {d} \mathbf {x} (t)\right]\mathrm {d} t-\mathbf {\lambda } ^{\mathsf {T}}(t_{1})\mathrm {d} \mathbf {x} (t_{1})+\mathbf {\lambda } ^{\mathsf {T}}(t_{0})\mathrm {d} \mathbf {x} (t_{0})\leq 0$

इस अभिव्यक्ति के लिए शून्य के बराबर होने के लिए निम्नलिखित इष्टतमता शर्तों की आवश्यकता होती है:

${\begin{aligned}I_{\mathbf {u} }(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} _{\mathbf {u} }(\mathbf {x} (t),\mathbf {u} (t),t)&=0\\I_{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} _{\mathbf {x} }(\mathbf {x} (t),\mathbf {u} (t),t)+{\dot {\mathbf {\lambda } }}(t)&=0\end{aligned}}$

यदि दोनों प्रारंभिक मान $\mathbf {x} (t_{0})$ और टर्मिनल मान $\mathbf {x} (t_{1})$ निश्चित हैं, अर्थात् $\mathrm {d} \mathbf {x} (t_{0})=\mathrm {d} \mathbf {x} (t_{1})=0$, कोई शर्त नहीं है $\mathbf {\lambda } (t_{0})$ और $\mathbf {\lambda } (t_{1})$ आवश्यक है। यदि टर्मिनल मूल्य मुक्त है, जैसा कि अधिकांशतः होता है, अतिरिक्त शर्त $\mathbf {\lambda } (t_{1})=0$ श्रेष्ठता के लिए आवश्यक है। उत्तरार्द्ध को निश्चित क्षितिज समस्या के लिए ट्रांसवर्सलिटी स्थिति कहा जाता है।

यह देखा जा सकता है कि आवश्यक शर्तें हैमिल्टनियन के लिए ऊपर बताई गई शर्तों के समान हैं। इस प्रकार हैमिल्टनियन को पहले क्रम की आवश्यक शर्तों को उत्पन्न करने के लिए उपकरण के रूप में समझा जा सकता है।

असतत समय में हैमिल्टनियन

जब समस्या असतत समय में तैयार की जाती है, तो हैमिल्टनियन को इस प्रकार परिभाषित किया जाता है:

$H(x_{t},u_{t},\lambda _{t+1},t)=\lambda _{t+1}^{\top }f(x_{t},u_{t},t)+I(x_{t},u_{t},t)\,$

और कॉस्टेट समीकरण हैं

$\lambda _{t}={\frac {\partial H}{\partial x_{t}}}$

(ध्यान दें कि समय पर असतत समय हैमिल्टनियन $t$ समय पर कॉस्टेट वैरिएबल सम्मिलित है $t+1.$ यह छोटा विवरण आवश्यक है जिससे जब हम इसके संबंध में अंतर करें $x$ हमें शब्द सम्मिलित है $\lambda (t+1)$ कॉस्टेट समीकरणों के दाहिने हाथ की ओर। यहां गलत सम्मेलन का उपयोग करने से गलत परिणाम हो सकते हैं, यानी कॉस्टेट समीकरण जो पीछे की ओर अंतर समीकरण नहीं है)।

हैमिल्टनियन का समय के साथ व्यवहार

पोंट्रीगिन के अधिकतम सिद्धांत से, हैमिल्टनियन के लिए विशेष शर्तें प्राप्त की जा सकती हैं। जब आखिरी बार $t_{1}$ निश्चित है और हैमिल्टन समय पर स्पष्ट रूप से निर्भर नहीं करता है $\left({\tfrac {\partial H}{\partial t}}=0\right)$, तब:

$H(x^{*}(t),u^{*}(t),\lambda ^{*}(t))=\mathrm {constant} \,$

या यदि टर्मिनल समय निःशुल्क है, तो:

$H(x^{*}(t),u^{*}(t),\lambda ^{*}(t))=0.\,$

इसके अतिरिक्त, यदि टर्मिनल समय अनंत तक जाता है, तो हैमिल्टनियन पर पारलौकिक स्थिति स्थिति प्रयुक्त होती है।

$\lim _{t\to \infty }H(t)=0$

यांत्रिकी के हैमिल्टनियन की तुलना में नियंत्रण का हैमिल्टनियन

विलियम रोवन हैमिल्टन ने प्रणाली के यांत्रिकी का वर्णन करने के लिए हैमिल्टनियन यांत्रिकी को परिभाषित किया। यह तीन चरों का कार्य है:

${\mathcal {H}}={\mathcal {H}}(p,q,t)=\langle p,{\dot {q}}\rangle -L(q,{\dot {q}},t)$

जहाँ $L$ लाग्रंगियन यांत्रिकी है, जिसका चरमोत्कर्ष गतिकी को निर्धारित करता है (ऊपर परिभाषित लाग्रंगियन नहीं), $q$ स्थिति चर है और ${\dot {q}}$ इसका काल व्युत्पन्न है।

$p$ तथाकथित संयुग्म गति है, द्वारा परिभाषित

$p={\frac {\partial L}{\partial {\dot {q}}}}$

हैमिल्टन ने तब सिस्टम की गतिशीलता का वर्णन करने के लिए अपने समीकरण तैयार किए

${\frac {d}{dt}}p(t)=-{\frac {\partial }{\partial q}}{\mathcal {H}}$

${\frac {d}{dt}}q(t)=~~{\frac {\partial }{\partial p}}{\mathcal {H}}$

नियंत्रण सिद्धांत का हैमिल्टन प्रणाली की गतिशीलता का वर्णन नहीं करता है, लेकिन नियंत्रण चर के संबंध में कुछ स्केलर फलन (लैग्रैंगियन) को चरम पर पहुंचाने की स्थिति $u$. जैसा कि सामान्य रूप से परिभाषित किया गया है, यह 4 चरों का कार्य है

$H(q,u,p,t)=\langle p,{\dot {q}}\rangle -L(q,u,t)$

जहाँ $q$ स्थिति चर है और $u$ नियंत्रण चर है जिसके संबंध में हम चरम सीमा पर हैं।

अधिकतम के लिए संबद्ध शर्तें हैं

${\frac {dp}{dt}}=-{\frac {\partial H}{\partial q}}$

${\frac {dq}{dt}}=~~{\frac {\partial H}{\partial p}}$

${\frac {\partial H}{\partial u}}=0$

यह परिभाषा सस्मान और विलेम्स के लेख द्वारा दी गई परिभाषा से सहमत है। (पृष्ठ 39 देखें, समीकरण 14)। सुस्मान और विलेम्स दिखाते हैं कि हैमिल्टनियन नियंत्रण को गतिशीलता में कैसे प्रयोग किया जा सकता है उदा। ब्राचिस्टोक्रोन समस्या के लिए, लेकिन इस दृष्टिकोण पर कैराथोडोरी के पूर्व कार्य का उल्लेख न करें।

वर्तमान मूल्य और वर्तमान मूल्य हैमिल्टनियन

अर्थशास्त्र में, गतिशील अनुकूलन समस्याओं में उद्देश्य कार्य अधिकांशतः केवल घातीय छूट के माध्यम से सीधे समय पर निर्भर करता है, जैसे कि यह रूप लेता है

$I(\mathbf {x} (t),\mathbf {u} (t),t)=e^{-\rho t}\nu (\mathbf {x} (t),\mathbf {u} (t))$

जहाँ $\nu (\mathbf {x} (t),\mathbf {u} (t))$ तात्क्षणिक उपयोगिता फलन या परमानंद फलन के रूप में जाना जाता है। यह हैमिल्टनियन को फिर से परिभाषित करने की अनुमति देता है $H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)=e^{-\rho t}{\bar {H}}(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t))$ जहाँ

${\begin{aligned}{\bar {H}}(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t))\equiv &\,e^{\rho t}\left[I(\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\lambda } ^{\mathsf {T}}(t)\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)\right]\\=&\,\nu (\mathbf {x} (t),\mathbf {u} (t),t)+\mathbf {\mu } ^{\mathsf {T}}(t)\mathbf {f} (\mathbf {x} (t),\mathbf {u} (t),t)\end{aligned}}$

जिसे हैमिल्टनियन के वर्तमान मूल्य के विपरीत, वर्तमान मूल्य हैमिल्टनियन कहा जाता है $H(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t),t)$ पहले खंड में परिभाषित। विशेष रूप से कॉस्टेट वेरिएबल्स को फिर से परिभाषित किया गया है $\mathbf {\mu } (t)=e^{\rho t}\mathbf {\lambda } (t)$, जो संशोधित प्रथम-क्रम स्थितियों की ओर जाता है।

${\frac {\partial {\bar {H}}(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t))}{\partial \mathbf {u} }}=0$,

${\frac {\partial {\bar {H}}(\mathbf {x} (t),\mathbf {u} (t),\mathbf {\lambda } (t))}{\partial \mathbf {x} }}=-{\dot {\mathbf {\mu } }}(t)+\rho \mathbf {\mu } (t)$

जो उत्पाद नियम से तुरंत अनुसरण करता है। आर्थिक रूप से, $\mathbf {\mu } (t)$ पूंजीगत वस्तुओं के लिए वर्तमान-मूल्यवान छाया कीमतों $\mathbf {x} (t)$ का प्रतिनिधित्व करते हैं .

उदाहरण: रैमसे-कैस-कूपमन्स मॉडल

अर्थशास्त्र में, रैमसे-कैस-कूपमन्स मॉडल का उपयोग अर्थव्यवस्था के लिए इष्टतम बचत व्यवहार निर्धारित करने के लिए किया जाता है। उद्देश्य फलन $J(c)$ सामाजिक कल्याण कार्य है,

$J(c)=\int _{0}^{T}e^{-\rho t}u(c(t))dt$

इष्टतम उपभोग पथ के चुनाव द्वारा अधिकतम किया जाना $c(t)$. कार्यक्रम $u(c(t))$ उपभोग के प्रतिनिधि एजेंट उपयोगिता को इंगित करता है $c$ किसी भी समय पर। कारण $e^{-\rho t}$ छूट का प्रतिनिधित्व करता है। अधिकतमकरण समस्या पूंजी तीव्रता के लिए निम्नलिखित अंतर समीकरण के अधीन है, जो प्रति प्रभावी कार्यकर्ता पूंजी के समय के विकास का वर्णन करती है:

${\dot {k}}={\frac {\partial k}{\partial t}}=f(k(t))-(n+\delta )k(t)-c(t)$

जहाँ $c(t)$ अवधि टी खपत है, $k(t)$ प्रति कर्मचारी अवधि टी पूंजी है (के साथ $k(0)=k_{0}>0$), $f(k(t))$ अवधि टी उत्पादन है, $n$ जनसंख्या वृद्धि दर है, $\delta$ पूंजी मूल्यह्रास दर है, एजेंट भविष्य की उपयोगिता दर पर छूट देता है $\rho$, साथ $u'>0$ और $u''<0$.

यहाँ, $k(t)$ स्थिति चर है जो उपरोक्त समीकरण के अनुसार विकसित होता है, और $c(t)$ नियंत्रण चर है। हैमिल्टनियन बन जाता है

$H(k,c,\mu ,t)=e^{-\rho t}u(c(t))+\mu (t){\dot {k}}=e^{-\rho t}u(c(t))+\mu (t)[f(k(t))-(n+\delta )k(t)-c(t)]$

इष्टतम स्थिति हैं

${\frac {\partial H}{\partial c}}=0\Rightarrow e^{-\rho t}u'(c)=\mu (t)$

${\frac {\partial H}{\partial k}}=-{\frac {\partial \mu }{\partial t}}=-{\dot {\mu }}\Rightarrow \mu (t)[f'(k)-(n+\delta )]=-{\dot {\mu }}$

ट्रांसवर्सलिटी कंडीशन के अतिरिक्त $\mu (T)k(T)=0$. अगर हम जाने दें $u(c)=\log(c)$, फिर लॉगरिदमिक विभेदीकरण | लॉग-डिफरेंशियेटिंग पहली इष्टतम स्थिति के संबंध में $t$ पैदावार

$-\rho -{\frac {\dot {c}}{c(t)}}={\frac {\dot {\mu }}{\mu (t)}}$

इस समीकरण को दूसरी अनुकूलतम स्थिति में सम्मिलित करने से प्राप्त होता है

$\rho +{\frac {\dot {c}}{c(t)}}=f'(k)-(n+\delta )$

जिसे कीन्स-रैमसे नियम के रूप में जाना जाता है, जो हर अवधि में खपत के लिए शर्त देता है, जिसका पालन करने पर अधिकतम जीवनकाल उपयोगिता सुनिश्चित होती है।

अग्रिम पठन

• Léonard, Daniel; Long, Ngo Van (1992). "The Maximum Principle". Optimal Control Theory and Static Optimization in Economics. New York: Cambridge University Press. pp. 127–168. ISBN 0-521-33158-7.
• Takayama, Akira (1985). "Developments of Optimal Control Theory and Its Applications". Mathematical Economics (2nd ed.). New York: Cambridge University Press. pp. 600–719. ISBN 0-521-31498-4.
• Wulwick, Nancy (1995). "The Hamiltonian Formalism and Optimal Growth Theory". In Rima, I. H. (ed.). Measurement, Quantification, and Economic Analysis. London: Routledge. ISBN 978-0-415-08915-9.