आव्यूह अवकल समीकरण: Difference between revisions

From Vigyanwiki
(minor changes)
(Work done)
Line 37: Line 37:
n = 2 स्थिति में (दो अवस्था चर के साथ), स्थिरता की स्थिति है कि परिवर्ती आव्यूह A के दो आइगेनमान ​​प्रत्येक में एक ऋणात्मक वास्तविक भाग होता है, इस स्थिति के समकक्ष होता है कि A का [[ ट्रेस (गणित) |ट्रेस]] ऋणात्मक होना चाहिए और इसकी सारणिक धनात्मक होनी चाहिए।
n = 2 स्थिति में (दो अवस्था चर के साथ), स्थिरता की स्थिति है कि परिवर्ती आव्यूह A के दो आइगेनमान ​​प्रत्येक में एक ऋणात्मक वास्तविक भाग होता है, इस स्थिति के समकक्ष होता है कि A का [[ ट्रेस (गणित) |ट्रेस]] ऋणात्मक होना चाहिए और इसकी सारणिक धनात्मक होनी चाहिए।


== आव्यूह रूप में समाधान ==
== आव्यूह रूप में हल ==
का औपचारिक समाधान <math>\mathbf{\dot{x}}(t)=\mathbf{A}[\mathbf{x}(t)-\mathbf{x}^*]</math> आव्यूह घातीय रूप है
<math>\mathbf{\dot{x}}(t)=\mathbf{A}[\mathbf{x}(t)-\mathbf{x}^*]</math> के औपचारिक हल में आव्यूह घातीय रूप निम्नलिखित है
:<math>\mathbf{x}(t)=\mathbf{x}^*+e^{\mathbf{A}t}[\mathbf{x}(0)-\mathbf{x}^*] ~,</math> कई तकनीकों का उपयोग करके मूल्यांकन किया गया।
:<math>\mathbf{x}(t)=\mathbf{x}^*+e^{\mathbf{A}t}[\mathbf{x}(0)-\mathbf{x}^*] ~,</math>
विभिन्न तकनीकों में से किसी का उपयोग करके मूल्यांकन किया गया।


=== कंप्यूटिंग के लिए पुत्जर एल्गोरिथम {{math|e<sup>'''A'''''t''</sup>}}===
=== {{math|e<sup>'''A'''''t''</sup>}} की गणना के लिए पुट्जर एल्गोरिथम===
eigenvalues ​​​​के साथ एक आव्यूह ए दिया <math>\lambda_1,\lambda_2,\dots,\lambda_n</math>,   
<math>\lambda_1,\lambda_2,\dots,\lambda_n</math> आइगेनमानों वाला आव्यूह A दिया गया है,   
:<math>e^{\mathbf{A}t} = \sum_{j=0}^{n-1}r_{j+1}{\left(t\right)}\mathbf{P}_{j}</math>
:<math>e^{\mathbf{A}t} = \sum_{j=0}^{n-1}r_{j+1}{\left(t\right)}\mathbf{P}_{j}</math>
कहां
जहाँ
:<math>\mathbf{P}_0 = \mathbf{I}</math>
:<math>\mathbf{P}_0 = \mathbf{I}</math>
:<math>\mathbf{P}_j = \prod_{k=1}^{j}\left(\mathbf{A}-\lambda_k \mathbf{I}\right)= \mathbf{P}_{j-1} \left(\mathbf{A}-\lambda_j \mathbf{I}\right),  \qquad j=1,2,\dots,n-1</math>
:<math>\mathbf{P}_j = \prod_{k=1}^{j}\left(\mathbf{A}-\lambda_k \mathbf{I}\right)= \mathbf{P}_{j-1} \left(\mathbf{A}-\lambda_j \mathbf{I}\right),  \qquad j=1,2,\dots,n-1</math>
Line 51: Line 52:
:<math>\dot{r}_{j} = \lambda_j r_j + r_{j-1}, \qquad  j=2,3,\dots,n</math>
:<math>\dot{r}_{j} = \lambda_j r_j + r_{j-1}, \qquad  j=2,3,\dots,n</math>
:<math>r_j{\left(0\right)}=0,  \qquad  j=2,3,\dots,n</math>
:<math>r_j{\left(0\right)}=0,  \qquad  j=2,3,\dots,n</math>
के लिए समीकरण <math>r_i (t)</math> साधारण प्रथम कोटि के विषम ODE हैं।
<math>r_i (t)</math> के समीकरण साधारण प्रथम कोटि के विषम ओडीई हैं।


ध्यान दें कि एल्गोरिदम की आवश्यकता नहीं है कि आव्यूह [[ विकर्ण मैट्रिक्स | विकर्ण आव्यूह]] हो और आमतौर पर उपयोग किए जाने वाले जॉर्डन कैनोलिक रूपों की जटिलताओं को छोड़ दें।
ध्यान दें कि एल्गोरिथम के लिए यह आवश्यक नहीं है कि आव्यूह A [[ विकर्ण मैट्रिक्स |विकर्णीय]] हो और सामान्य रूप से उपयोग किए जाने वाले जॉर्डन कैनोनिकल रूपों की जटिलताओं को बाईपास कर दिया जाता है।


== एक आव्यूह साधारण अवकल समीकरण == का विखंडित उदाहरण
== आव्यूह साधारण अवकल समीकरण का विखंडित उदाहरण ==
 
दो फलनों x(t) और y(t) में एक प्रथम-कोटि सघातीय आव्यूह साधारण अवकल समीकरण, जब आव्यूह रूप में लिया जाता है, तो इसका निम्नलिखित रूप प्राप्त होता है:
दो कार्यों x(t) और y(t) में एक प्रथम-क्रम सजातीय आव्यूह साधारण अवकल समीकरण, जब आव्यूह रूप से बाहर निकाला जाता है, तो इसका निम्न रूप होता है:


: <math>\frac{dx}{dt}=a_1x+b_1y,\quad\frac{dy}{dt}=a_2x+b_2y</math>
: <math>\frac{dx}{dt}=a_1x+b_1y,\quad\frac{dy}{dt}=a_2x+b_2y</math>
कहां <math>a_1</math>, <math>a_2</math>, <math>b_1</math>, और <math>b_2</math> कोई मनमाना स्केलर हो सकता है।
जहाँ <math>a_1</math>, <math>a_2</math>, <math>b_1</math>, और <math>b_2</math> कोई भी यादृच्छिक अदिश हो सकते हैं।


उच्च क्रम आव्यूह ओडीई के पास अधिक जटिल रूप हो सकता है।
उच्च कोटि आव्यूह ओडीई का अधिक जटिल रूप हो सकता है।


== विघटित आव्यूह साधारण अवकल समीकरणों को हल करना ==
== विघटित आव्यूह साधारण अवकल समीकरणों को हल करना ==


उपरोक्त समीकरणों को हल करने और इस विशेष क्रम और रूप के आवश्यक कार्यों को खोजने की प्रक्रिया में 3 मुख्य चरण होते हैं। इनमें से प्रत्येक चरण का संक्षिप्त विवरण नीचे सूचीबद्ध है:
उपरोक्त समीकरणों को हल करने और इस विशेष क्रम और रूप के आवश्यक फलनों को प्राप्त करने की प्रक्रिया में 3 मुख्य चरण होते हैं। इन चरणों में से प्रत्येक का संक्षिप्त विवरण नीचे सूचीबद्ध किया गया हैं:


*[[ eigenvalues ]] का पता लगाना
*[[ eigenvalues |आइगेनमानों]] को ज्ञात करना
* [[ egenvectors ]] ढूँढना
*[[ egenvectors |आइगेनसदिशों]] को ज्ञात करना
* आवश्यक कार्यों का पता लगाना
* आवश्यक फलनों को ज्ञात करना


इस तरह के सामान्य अवकल समीकरणों को हल करने में अंतिम, तीसरा, चरण आमतौर पर पिछले दो चरणों में गणना किए गए मानों को एक विशेष सामान्य रूप समीकरण में प्लगिंग के माध्यम से किया जाता है, जिसका उल्लेख इस लेख में बाद में किया गया है।
इस प्रकार के साधारण अवकल समीकरणों को हल करने में अंतिम, तीसरा, चरण सामान्य रूप से दो पूर्व चरणों में गणना किए गए मानों को एक विशेष सामान्य रूप समीकरण में प्लगिंग के माध्यम से किया जाता है, जिसका उल्लेख इस आलेख में बाद में किया गया है।


== एक आव्यूह ODE == का हल उदाहरण
== आव्यूह ओडीई का हल किया गया उदाहरण ==
{{see also|Matrix exponential#Linear differential equations}}
{{see also|आव्यूह चरघातांक#रैखिक अवकल समीकरण}}
ऊपर दिए गए तीन चरणों के अनुसार एक आव्यूह ODE को हल करने के लिए, प्रक्रिया में सरल मैट्रिसेस का उपयोग करते हुए, हम एक फ़ंक्शन पाते हैं {{mvar|x}} और एक समारोह {{mvar|y}} दोनों एकल स्वतंत्र चर के संदर्भ में {{mvar|t}}, पहले क्रम के निम्नलिखित सजातीय [[ रैखिक अंतर समीकरण | रैखिक अवकल समीकरण]] में,
प्रक्रिया में सरल आव्यूहों का उपयोग करते हुए, ऊपर दिए गए तीन चरणों के अनुसार आव्यूह ओडीई को हल करने के लिए, मान लीजिए, एक फलन {{mvar|x}} और एक फलन {{mvar|y}} दोनों को एक स्वतंत्र चर {{mvar|t}} के संदर्भ में प्रथम कोटि के निम्नलिखित सघातीय [[ रैखिक अंतर समीकरण |रैखिक अवकल समीकरण]] में प्राप्त होते हैं।
: <math>\frac{dx}{dt}=3x-4y,\quad\frac{dy}{dt}=4x-7y~.</math>
: <math>\frac{dx}{dt}=3x-4y,\quad\frac{dy}{dt}=4x-7y~.</math>
इस विशेष सामान्य अवकल समीकरण प्रणाली को हल करने के लिए, समाधान प्रक्रिया में किसी बिंदु पर, हमें दो प्रारंभिक स्थितियों (प्रारंभिक बिंदु पर दो राज्य चर के अनुरूप) के एक सेट की आवश्यकता होगी। इस मामले में, हम चुनते हैं {{math|1=''x''(0) = ''y''(0) = 1}}.
इस विशेष सामान्य अवकल समीकरण प्रणाली को हल करने के लिए, हल प्रक्रिया में किसी बिंदु पर, हमें दो प्रारंभिक मानों के एक समुच्चय की आवश्यकता होगी (शुरुआती बिंदु पर दो अवस्था चर के अनुरूप)इस स्थिति में, हम {{math|1=''x''(0) = ''y''(0) = 1}} चयनित करते हैं।


=== पहला कदम ===
=== प्रथम चरण ===
पहला चरण, जिसका पहले ही ऊपर उल्लेख किया जा चुका है, में A के eigenvalues ​​का पता लगा रहा है
प्रथम चरण, जिसका पहले ही ऊपर उल्लेख किया जा चुका है, में A का आइगेनमान ज्ञात करना है
: <math>\begin{bmatrix} x'\\y' \end{bmatrix} = \begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix}\begin{bmatrix} x\\y \end{bmatrix}~.</math>
: <math>\begin{bmatrix} x'\\y' \end{bmatrix} = \begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix}\begin{bmatrix} x\\y \end{bmatrix}~.</math>
उपरोक्त सदिशों में से किसी एक में देखे गए व्युत्पन्न संकेतन x' आदि को लाग्रेंज के संकेतन के रूप में जाना जाता है (पहली बार [[ जोसेफ लुइस लाग्रेंज ]] द्वारा प्रस्तुत किया गया। यह पिछले समीकरण में प्रयुक्त व्युत्पन्न संकेतन dx/dt के बराबर है, जिसे लाइबनिज के संकेतन के रूप में जाना जाता है, सम्मान करते हुए [[ गॉटफ्रीड लीबनिज ]] का नाम।)
उपरोक्त सदिशों में से किसी एक में देखे गए अवकलज अंकन x' आदि को लाग्रेंज के संकेतन के रूप में जाना जाता है (सबसे पहले [[ जोसेफ लुइस लाग्रेंज |जोसेफ लुइस लाग्रेंज]] द्वारा प्रस्तुत किया गया। यह अवकलज संकेतन dx/dt के बराबर है जिसका उपयोग पूर्व समीकरण में किया गया था, जिसे लेबनिज़ के संकेतन के रूप में जाना जाता है, [[ गॉटफ्रीड लीबनिज |गॉटफ्राइड लाइबनिज़]] के नाम का सम्मान करते हुए।)


एक बार ऊपर प्रदर्शित [[ मैट्रिक्स (गणित) | आव्यूह (गणित)]] फॉर्म 'ए' में दो चर के गुणांक लिखे जाने के बाद, आइगेनवैल्यू का मूल्यांकन किया जा सकता है। उस अंत तक, एक आव्यूह (गणित) के निर्धारक को ढूंढता है जो एक पहचान आव्यूह के दौरान बनता है, <math>I_n</math>, कुछ स्थिर से गुणा {{mvar|λ}}, उपरोक्त गुणांक आव्यूह से घटाया जाता है ताकि इसकी [[ विशेषता बहुपद ]] प्राप्त हो सके,
एक बार जब दो चर के गुणांक ऊपर प्रदर्शित [[ मैट्रिक्स (गणित) |आव्यूह]] रूप A में लिखे गए हों, तो आइगेनमान का मूल्यांकन किया जा सकता है। उस अंत तक, एक आव्यूह के निर्धारक को प्राप्त करता है जो तब बनता है जब एक तत्समक आव्यूह, <math>I_n</math>, कुछ स्थिरांक λ से गुणा किया जाता है, इसे उपरोक्त गुणांक आव्यूह से घटाया जाता है ताकि इसकी [[ विशेषता बहुपद |अभिलाक्षणिक बहुपद]] प्राप्त हो सके,
: <math>\det\left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} - \lambda\begin{bmatrix} 1 & 0\\0 & 1 \end{bmatrix}\right)~,</math>
: <math>\det\left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} - \lambda\begin{bmatrix} 1 & 0\\0 & 1 \end{bmatrix}\right)~,</math>
और इसके शून्य के लिए हल करें।
और इसके शून्य के लिए हल करें।


[[ मैट्रिक्स जोड़ | आव्यूह जोड़]] पैदावार के और सरलीकरण और बुनियादी नियमों को लागू करना
आगे सरलीकरण और [[ मैट्रिक्स जोड़ |आव्यूह योग]] प्राप्ति के मूल नियमों को लागू करना
: <math>\det\begin{bmatrix} 3-\lambda & -4\\4 & -7-\lambda \end{bmatrix}~.</math>
: <math>\det\begin{bmatrix} 3-\lambda & -4\\4 & -7-\lambda \end{bmatrix}~.</math>
एकल 2×2 आव्यूह के निर्धारक को खोजने के नियमों को लागू करने से निम्नलिखित प्राथमिक [[ द्विघात समीकरण ]] प्राप्त होता है,
एकल 2×2 आव्यूह के निर्धारक को ज्ञात करने के नियमों को लागू करने से निम्नलिखित प्रारंभिक [[ द्विघात समीकरण |द्विघात समीकरण]] उत्पन्न होता है,
: <math>\det\begin{bmatrix} 3-\lambda & -4\\4 & -7-\lambda \end{bmatrix} = 0</math>
: <math>\det\begin{bmatrix} 3-\lambda & -4\\4 & -7-\lambda \end{bmatrix} = 0</math>
: <math>-21 - 3\lambda + 7\lambda + \lambda^2 + 16 = 0  \,\!</math>
: <math>-21 - 3\lambda + 7\lambda + \lambda^2 + 16 = 0  \,\!</math>
उपरोक्त का एक सरल संस्करण प्राप्त करने के लिए इसे और कम किया जा सकता है,
जिसे उपरोक्त का सरल संस्करण प्राप्त करने के लिए और सरलीकृत किया जा सकता है,
: <math>\lambda^2 + 4\lambda - 5 = 0  ~.</math>
: <math>\lambda^2 + 4\lambda - 5 = 0  ~.</math>
अब दो मूल ज्ञात करना, <math>\lambda_1</math> और <math>\lambda_2</math> दिए गए द्विघात समीकरण का [[ गुणन ]]खंडन विधि लागू करने से प्राप्त होता है
अब [[ गुणन |गुणनखंड]] विधि को लागू करके दिए गए द्विघात समीकरण के दो मूल, <math>\lambda_1</math> और <math>\lambda_2</math> को ज्ञात करने पर प्राप्त होता है
: <math>\lambda^2 + 5\lambda - \lambda - 5 = 0</math>
: <math>\lambda^2 + 5\lambda - \lambda - 5 = 0</math>
: <math>\lambda (\lambda + 5) - 1 (\lambda + 5) = 0</math>
: <math>\lambda (\lambda + 5) - 1 (\lambda + 5) = 0</math>
: <math>(\lambda - 1)(\lambda + 5) = 0</math>
: <math>(\lambda - 1)(\lambda + 5) = 0</math>
: <math>\lambda = 1, -5 ~.</math>
: <math>\lambda = 1, -5 ~.</math>
मूल्य <math>\lambda_1 = 1</math> और <math>\lambda_2 = -5</math>, ऊपर परिकलित A के आवश्यक eigenvalues ​​​​हैं।
ऊपर परिकलित किए गए मान <math>\lambda_1 = 1</math> और <math>\lambda_2 = -5</math>, A के आवश्यक आइगेनमान हैं। कुछ स्थितियों में, माना अन्य आव्यूह ओडीई, आइगेनमान समिश्र हो सकते हैं, इस स्थिति में हल प्रक्रिया के अगले चरण के साथ-साथ अंतिम रूप और हल प्रभावशाली रूप से परिवर्तित कर सकते हैं।
कुछ मामलों में, अन्य आव्यूह ओडीई का कहना है, ईगेनवेल्यूज जटिल संख्याएं हो सकती हैं, इस मामले में हल करने की प्रक्रिया के अगले चरण के साथ-साथ अंतिम रूप और समाधान नाटकीय रूप से बदल सकते हैं।


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


गणना किए गए प्रत्येक eigenvalues ​​​​के लिए, हमारे पास एक अलग eigenvector है। पहले eigenvalue के लिए, जो है <math>\lambda_1 = 1</math>, अपने पास
गणना किए गए प्रत्येक आइगेनमान के लिए, हमारे पास एक विशिष्ट आइगेनसदिश है। पहले आइगेनमान के लिए, जो <math>\lambda_1 = 1</math> है, हमें निम्न प्राप्त है
: <math>\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix}\begin{bmatrix} \alpha\\\beta \end{bmatrix} = 1\begin{bmatrix} \alpha\\\beta \end{bmatrix}. </math>
: <math>\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix}\begin{bmatrix} \alpha\\\beta \end{bmatrix} = 1\begin{bmatrix} \alpha\\\beta \end{bmatrix}. </math>
मूल आव्यूह गुणन नियम लागू करके उपरोक्त व्यंजक को सरल करने पर प्राप्त होता है
मूल आव्यूह गुणन नियमों को लागू करके उपरोक्त व्यंजक को सरल बनाने पर प्राप्त होता है
: <math>3\alpha - 4\beta = \alpha</math>
: <math>3\alpha - 4\beta = \alpha</math>
: <math>\alpha = 2\beta~.</math>
: <math>\alpha = 2\beta~.</math>
ये सभी गणनाएँ केवल अंतिम अभिव्यक्ति प्राप्त करने के लिए की गई हैं, जो हमारे मामले में है {{math|1=''α'' = 2''β''}}. अब कुछ मनमाना मूल्य लेना, संभवतः, एक छोटा महत्वहीन मूल्य, जिसके साथ काम करना बहुत आसान है, दोनों के लिए {{mvar|α}} या {{mvar|β}} (ज्यादातर मामलों में, यह वास्तव में कोई फर्क नहीं पड़ता), हम इसे प्रतिस्थापित करते हैं {{math|1=''α'' = 2''β''}}. ऐसा करने से एक साधारण सदिश उत्पन्न होता है, जो कि इस विशेष eigenvalue के लिए आवश्यक eigenvector है। हमारे मामले में, हम चुनते हैं {{math|1=''α'' = 2}}, जो बदले में यह निर्धारित करता है {{math|1=''β'' = 1}} और, मानक सदिश संकेतन का उपयोग करते हुए, हमारा सदिश जैसा दिखता है
ये सभी गणना केवल अंतिम व्यंजक प्राप्त करने के लिए की गई है, जो हमारी स्थिति में {{math|1=''α'' = 2''β''}} है। अब कुछ यादृच्छिक मान लेते हुए, संभवतः, एक छोटा महत्वहीन मान, जिसके साथ कार्य करना बहुत सरल होता है, {{mvar|α}} या {{mvar|β}} के लिए (अधिकतर स्थितियों में, यह वास्तव में कोई मायने नहीं रखता), हम इसे {{math|1=''α'' = 2''β''}} में प्रतिस्थापित करते हैं। ऐसा करने से एक साधारण सदिश उत्पन्न होता है, जो इस विशेष आइगेनमान के लिए आवश्यक आइगेनसदिश होता है। हमारी स्थिति में, हम {{math|1=''α'' = 2}} चयनित करते हैं, जो बदले में यह निर्धारित करता है कि {{math|1=''β'' = 1}} और, मानक सदिश संकेतन का उपयोग करके, हमारा सदिश कुछ इस प्रकार का होता है
: <math>\mathbf{\hat{v}}_1 = \begin{bmatrix} 2\\1 \end{bmatrix}. </math>
: <math>\mathbf{\hat{v}}_1 = \begin{bmatrix} 2\\1 \end{bmatrix}. </math>
हमारे द्वारा गणना की गई दूसरी आइगेनवेल्यू का उपयोग करके उसी ऑपरेशन को करना, जो है <math>\lambda = -5</math>, हम अपना दूसरा ईजेनवेक्टर प्राप्त करते हैं। इस यूक्लिडियन सदिश को निकालने की प्रक्रिया नहीं दिखाई गई है, लेकिन अंतिम परिणाम है
हमारे द्वारा परिकलित किए गए दूसरे आइगेनमान, जो कि <math>\lambda = -5</math> है, का उपयोग करके उसी संक्रिया को निष्पादित करते हुए, हम अपना दूसरा आइगेनमान प्राप्त करते हैं। इस सदिश को प्राप्त करने की प्रक्रिया नहीं दिखाई गई है, लेकिन अंतिम परिणाम निम्न है
: <math>\mathbf{\hat{v}}_2 = \begin{bmatrix} 1\\2 \end{bmatrix}. </math>
: <math>\mathbf{\hat{v}}_2 = \begin{bmatrix} 1\\2 \end{bmatrix}. </math>
===तृतीय चरण===
यह अंतिम चरण उन आवश्यक फलनों को ज्ञात करता है जो मूल रूप से हमें दिए गए अवकलज के पीछे 'छिपे हुए' हैं। दो फलन हैं, क्योंकि हमारे अवकल समीकरण दो चर के साथ काम करते हैं।


 
वह समीकरण जिसमें वह सभी जानकारी सम्मिलित है जो हमने पहले प्राप्त की थी,  इसका का निम्न रूप है:
===तीसरा चरण===
यह अंतिम चरण उन आवश्यक कार्यों को ढूंढता है जो मूल रूप से हमें दिए गए डेरिवेटिव के पीछे 'छिपे हुए' हैं। दो कार्य हैं, क्योंकि हमारे अवकल समीकरण दो चरों से संबंधित हैं।
 
वह समीकरण जिसमें वह सभी जानकारी शामिल है जो हमें पहले मिली थी, उसका निम्न रूप है:
: <math>\begin{bmatrix} x\\y \end{bmatrix} = Ae^{\lambda_1t}\mathbf{\hat{v}}_1 + Be^{\lambda_2t}\mathbf{\hat{v}}_2. </math>
: <math>\begin{bmatrix} x\\y \end{bmatrix} = Ae^{\lambda_1t}\mathbf{\hat{v}}_1 + Be^{\lambda_2t}\mathbf{\hat{v}}_2. </math>
eigenvalues ​​​​और eigenvectors पैदावार के मूल्यों को प्रतिस्थापित करना
आइगेनमानों ​​और आइगेनसदिशों के मूल्यों को प्रतिस्थापित करने से प्राप्ति होती है
: <math>\begin{bmatrix} x\\y \end{bmatrix} = Ae^{t}\begin{bmatrix} 2\\1 \end{bmatrix} + Be^{-5t}\begin{bmatrix} 1\\2 \end{bmatrix}. </math>
: <math>\begin{bmatrix} x\\y \end{bmatrix} = Ae^{t}\begin{bmatrix} 2\\1 \end{bmatrix} + Be^{-5t}\begin{bmatrix} 1\\2 \end{bmatrix}. </math>
आगे सरलीकरण लागू करना,
और अधिक सरलीकरण लागू करना,
: <math>\begin{bmatrix} x\\y \end{bmatrix} = \begin{bmatrix} 2 & 1\\1 & 2 \end{bmatrix}\begin{bmatrix} Ae^{t}\\Be^{-5t} \end{bmatrix}. </math>
: <math>\begin{bmatrix} x\\y \end{bmatrix} = \begin{bmatrix} 2 & 1\\1 & 2 \end{bmatrix}\begin{bmatrix} Ae^{t}\\Be^{-5t} \end{bmatrix}. </math>
आगे सरल करना और कार्यों के लिए समीकरण लिखना {{mvar|x}} और {{mvar|y}} अलग से,
और अधिक सरल करना और फलन {{mvar|x}} और {{mvar|y}} के लिए अलग-अलग समीकरण लिखना,
: <math>x = 2Ae^{t} + Be^{-5t} </math>
: <math>x = 2Ae^{t} + Be^{-5t} </math>
: <math>y = Ae^{t} + 2Be^{-5t}.</math>
: <math>y = Ae^{t} + 2Be^{-5t}.</math>
उपरोक्त समीकरण, वास्तव में, मांगे गए सामान्य कार्य हैं, लेकिन वे अपने सामान्य रूप में हैं (अनिर्दिष्ट मूल्यों के साथ) {{mvar|A}} और {{mvar|B}}), जबकि हम वास्तव में उनके सटीक रूप और समाधान खोजना चाहते हैं। तो अब हम समस्या की दी गई प्रारंभिक स्थितियों पर विचार करते हैं (दिए गए प्रारंभिक शर्तों सहित समस्या तथाकथित [[ प्रारंभिक मूल्य समस्या ]] है)। मान लीजिए हमें दिया गया है <math>x(0) = y(0) = 1</math>, जो हमारे साधारण अवकल समीकरण के लिए शुरुआती बिंदु की भूमिका निभाता है; इन शर्तों का अनुप्रयोग स्थिरांक निर्दिष्ट करता है, {{mvar|A}} और {{mvar|B}}. जैसा कि हम से देखते हैं <math>x(0) = y(0) = 1</math> शर्तें, जब {{math|1=''t'' = 0}}, उपरोक्त समीकरणों के बाएँ पक्ष 1 के बराबर हैं। इस प्रकार हम रैखिक समीकरणों की निम्नलिखित प्रणाली का निर्माण कर सकते हैं,
उपर्युक्त समीकरण, वास्तव में, वांछित सामान्य फलन हैं, लेकिन वे अपने सामान्य रूप में होते हैं ({{mvar|A}} और {{mvar|B}} के अनिर्दिष्ट मूल्यों के साथ), जबकि हम वास्तव में उनके सटीक रूप और हल प्राप्त करना चाहते हैं। तो अब हम समस्या की दी गई प्रारंभिक शर्तों पर विचार करते हैं (दिए गए प्रारंभिक शर्तों सहित समस्या तथाकथित [[ प्रारंभिक मूल्य समस्या |प्रारंभिक मूल्य समस्या]] है)। मान लीजिए कि हमें <math>x(0) = y(0) = 1</math> दिया गया है, जो हमारे साधारण अवकल समीकरण के लिए शुरुआती बिंदु की भूमिका निभाता है; इन शर्तों का अनुप्रयोग स्थिरांक, {{mvar|A}} और {{mvar|B}} को निर्दिष्ट करता है। जैसा कि हम <math>x(0) = y(0) = 1</math> स्थितियों से देखते हैं, जब {{math|1=''t'' = 0}}, उपरोक्त समीकरणों के बाएँ पक्ष 1 के बराबर होते हैं। इस प्रकार हम रैखिक समीकरणों की निम्नलिखित प्रणाली का निर्माण कर सकते हैं,
: <math>1 = 2A + B </math>
: <math>1 = 2A + B </math>
: <math>1 = A + 2B~.</math>
: <math>1 = A + 2B~.</math>
इन समीकरणों को हल करने पर, हम पाते हैं कि दोनों स्थिरांक {{mvar|A}} और {{mvar|B}} बराबर 1/3। इसलिए इन मूल्यों को इन दो कार्यों के सामान्य रूप में प्रतिस्थापित करना
इन समीकरणों को हल करने पर, हम प्राप्त करते हैं कि दोनों स्थिरांक {{mvar|A}} और {{mvar|B}} 1/3 के बराबर हैं। इसलिए इन मानों को इन दो कार्यों के सामान्य रूप में प्रतिस्थापित करने से उनके सटीक रूप निर्दिष्ट होते हैं,<math display="block">x = \tfrac{2}{3}e^{t} + \tfrac{1}{3}e^{-5t} </math><math display="block">y = \tfrac{1}{3}e^{t} + \tfrac{2}{3}e^{-5t}~,</math>
उनके सटीक रूपों को निर्दिष्ट करता है,
 
<math display="block">x = \tfrac{2}{3}e^{t} + \tfrac{1}{3}e^{-5t} </math>
 
<math display="block">y = \tfrac{1}{3}e^{t} + \tfrac{2}{3}e^{-5t}~,</math>
वांछित दो फलन।
दो कार्यों की मांग की।


=== आव्यूह घातांक का उपयोग करना ===
=== आव्यूह घातांक का उपयोग करना ===
उपरोक्त समस्या को आव्यूह एक्सपोनेंशियल के सीधे आवेदन के साथ हल किया जा सकता था। यानी हम ऐसा कह सकते हैं
उपरोक्त समस्या को आव्यूह घातांक के अनुप्रयोग के साथ हल किया जा सकता था। अर्थात् हम ऐसा कह सकते हैं


  <math>\begin{bmatrix} x(t)\\y(t) \end{bmatrix} = \exp \left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} t\right) \begin{bmatrix} x_0(t)\\y_0(t) \end{bmatrix} </math>
  <math>\begin{bmatrix} x(t)\\y(t) \end{bmatrix} = \exp \left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} t\right) \begin{bmatrix} x_0(t)\\y_0(t) \end{bmatrix} </math>
यह देखते हुए (जिसे [[ MATLAB ]]'s जैसे किसी भी उपयुक्त उपकरण का उपयोग करके गणना की जा सकती है <code>expm</code> उपकरण, या [[ मैट्रिक्स विकर्णकरण | आव्यूह विकर्णकरण]] करके और संपत्ति का लाभ उठाकर कि एक विकर्ण आव्यूह का आव्यूह घातांक इसके तत्वों के तत्व-वार घातांक के समान है)
यह देखते हुए कि (जिसे किसी भी उपयुक्त उपकरण का उपयोग करके गणना की जा सकती है, जैसे कि [[ MATLAB |एमएटीएलएबी]] का <code>expm</code> टूल, या [[ मैट्रिक्स विकर्णकरण |आव्यूह विकर्णकरण]] करके और संपत्ति का लाभ उठाकर कि विकर्ण आव्यूह का आव्यूह घातांक इसके तत्वों के तत्व-वार घातांक के समान है)


  <math>\exp \left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} t\right) = \begin{bmatrix} 4 e^t/3 - e^{-5t}/3 & 2e^{-5t}/3 - 2e^t/3\\2e^t/3 - 2e^{-5t}/3 & 4e^{-5t}/3 - e^t/3 \end{bmatrix} </math>
  <math>\exp \left(\begin{bmatrix} 3 & -4\\4 & -7 \end{bmatrix} t\right) = \begin{bmatrix} 4 e^t/3 - e^{-5t}/3 & 2e^{-5t}/3 - 2e^t/3\\2e^t/3 - 2e^{-5t}/3 & 4e^{-5t}/3 - e^t/3 \end{bmatrix} </math>
Line 151: Line 147:


<math>\begin{bmatrix} x(t)\\y(t) \end{bmatrix} = \begin{bmatrix} e^{-5t}/3 + 2e^t/3\\ e^t/3 + 2e^{-5t}/3 \end{bmatrix} </math>
<math>\begin{bmatrix} x(t)\\y(t) \end{bmatrix} = \begin{bmatrix} e^{-5t}/3 + 2e^t/3\\ e^t/3 + 2e^{-5t}/3 \end{bmatrix} </math>
यह पहले दिखाए गए ईजेनवेक्टर दृष्टिकोण के समान है।
 
यह पहले दिखाए गए आइगेनसदिश दृष्टिकोण जैसा ही है।


== यह भी देखें ==
== यह भी देखें ==
* साधारण अवकल समीकरण#गैर-सजातीय समीकरण
* असमघातीय समीकरण
* [[ आव्यूह [[ अंतर समीकरण | अवकल समीकरण]]  ]]
* [[ अंतर समीकरण |आव्यूह अवकल समीकरण]]
* न्यूटन का शीतलन का नियम
* न्यूटन का शीतलन का नियम
* [[ फाइबोनैचि संख्या ]]
* [[ फाइबोनैचि संख्या |फिबोनाची अनुक्रम]]
* अवकल समीकरण
* अवकल समीकरण
* [[ तरंग समीकरण ]]
* [[ तरंग समीकरण |तरंग समीकरण]]
* [[ स्वायत्त प्रणाली (गणित) ]]
* [[ स्वायत्त प्रणाली (गणित) |स्वायत्त प्रणाली (गणित)]]
 


*
*

Revision as of 16:09, 11 January 2023

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

उदाहरण के लिए, एक प्रथम-कोटि आव्यूह साधारण अवकल समीकरण है

जहां अंतर्निहित चर के फलनों का सदिश होता है, इन फलनों के प्रथम अवकलज का सदिश है, और गुणांक का आव्यूह है।

ऐसी स्थिति में जहाँ नियतांक है और इसमें n एकघाततः स्वतंत्र आइगेनसदिश होता हैं, इस अवकल समीकरण का निम्नलिखित सामान्य हल है,

जहाँ λ1, λ2, …, λn A के आइगेनमान हैं; u1, u2, …, un A के संबंधित आइगेनवेक्टर हैं; और c1, c2, …, cn स्थिरांक हैं।

अत्यधिक सामान्य रूप से, यदि अपने समाकल के साथ विनिमय करता है तो मैग्नस विस्तार अग्रणी क्रम में कम हो जाता है, और अवकल समीकरण का सामान्य हल होता है

जहाँ एक एक नियत सदिश है।

कैले-हैमिल्टन प्रमेय और वैंडरमोंड-प्रकार आव्यूहों के उपयोग से, यह औपचारिक आव्यूह घातीय हल एक साधारण रूप में सरलीकृत किया जा सकता है।[1] नीचे, यह हल पुट्ज़ेर के एल्गोरिथम के संदर्भ में प्रदर्शित किया गया है।[2]

आव्यूह प्रणाली की स्थिरता और स्थिर अवस्था

आव्यूह समीकरण

n×1 पैरामीटर के साथ सतत सदिश b स्थिर है यदि और केवल तभी जब स्थिर आव्यूह A के सभी आइगेनमान ​​में ऋणात्मक वास्तविक भाग हो।

स्थिर स्थिति x* जिसमें स्थिर होने पर यह कन्वर्ज करता है, सेटिंग द्वारा पाया जाता है

इस प्रकार प्राप्त होता है

A को व्युत्क्रमणीय मानते हुए।

इस प्रकार, मूल समीकरण को सजातीय अवस्था से विचलन के संदर्भ में सजातीय रूप में लिखा जा सकता है,

इसे व्यक्त करने का एक समतुल्य तरीका यह है कि x* असमघात समीकरण का एक विशेष हल है, जबकि सभी हल इस रूप में हैं

के साथ सजातीय समीकरण (b=0) का हल।

दो-अवस्था-चर स्थिति में स्थिरता

n = 2 स्थिति में (दो अवस्था चर के साथ), स्थिरता की स्थिति है कि परिवर्ती आव्यूह A के दो आइगेनमान ​​प्रत्येक में एक ऋणात्मक वास्तविक भाग होता है, इस स्थिति के समकक्ष होता है कि A का ट्रेस ऋणात्मक होना चाहिए और इसकी सारणिक धनात्मक होनी चाहिए।

आव्यूह रूप में हल

के औपचारिक हल में आव्यूह घातीय रूप निम्नलिखित है

विभिन्न तकनीकों में से किसी का उपयोग करके मूल्यांकन किया गया।

eAt की गणना के लिए पुट्जर एल्गोरिथम

आइगेनमानों वाला आव्यूह A दिया गया है,

जहाँ

के समीकरण साधारण प्रथम कोटि के विषम ओडीई हैं।

ध्यान दें कि एल्गोरिथम के लिए यह आवश्यक नहीं है कि आव्यूह A विकर्णीय हो और सामान्य रूप से उपयोग किए जाने वाले जॉर्डन कैनोनिकल रूपों की जटिलताओं को बाईपास कर दिया जाता है।

आव्यूह साधारण अवकल समीकरण का विखंडित उदाहरण

दो फलनों x(t) और y(t) में एक प्रथम-कोटि सघातीय आव्यूह साधारण अवकल समीकरण, जब आव्यूह रूप में लिया जाता है, तो इसका निम्नलिखित रूप प्राप्त होता है:

जहाँ , , , और कोई भी यादृच्छिक अदिश हो सकते हैं।

उच्च कोटि आव्यूह ओडीई का अधिक जटिल रूप हो सकता है।

विघटित आव्यूह साधारण अवकल समीकरणों को हल करना

उपरोक्त समीकरणों को हल करने और इस विशेष क्रम और रूप के आवश्यक फलनों को प्राप्त करने की प्रक्रिया में 3 मुख्य चरण होते हैं। इन चरणों में से प्रत्येक का संक्षिप्त विवरण नीचे सूचीबद्ध किया गया हैं:

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

आव्यूह ओडीई का हल किया गया उदाहरण

प्रक्रिया में सरल आव्यूहों का उपयोग करते हुए, ऊपर दिए गए तीन चरणों के अनुसार आव्यूह ओडीई को हल करने के लिए, मान लीजिए, एक फलन x और एक फलन y दोनों को एक स्वतंत्र चर t के संदर्भ में प्रथम कोटि के निम्नलिखित सघातीय रैखिक अवकल समीकरण में प्राप्त होते हैं।

इस विशेष सामान्य अवकल समीकरण प्रणाली को हल करने के लिए, हल प्रक्रिया में किसी बिंदु पर, हमें दो प्रारंभिक मानों के एक समुच्चय की आवश्यकता होगी (शुरुआती बिंदु पर दो अवस्था चर के अनुरूप)। इस स्थिति में, हम x(0) = y(0) = 1 चयनित करते हैं।

प्रथम चरण

प्रथम चरण, जिसका पहले ही ऊपर उल्लेख किया जा चुका है, में A का आइगेनमान ज्ञात करना है

उपरोक्त सदिशों में से किसी एक में देखे गए अवकलज अंकन x' आदि को लाग्रेंज के संकेतन के रूप में जाना जाता है (सबसे पहले जोसेफ लुइस लाग्रेंज द्वारा प्रस्तुत किया गया। यह अवकलज संकेतन dx/dt के बराबर है जिसका उपयोग पूर्व समीकरण में किया गया था, जिसे लेबनिज़ के संकेतन के रूप में जाना जाता है, गॉटफ्राइड लाइबनिज़ के नाम का सम्मान करते हुए।)

एक बार जब दो चर के गुणांक ऊपर प्रदर्शित आव्यूह रूप A में लिखे गए हों, तो आइगेनमान का मूल्यांकन किया जा सकता है। उस अंत तक, एक आव्यूह के निर्धारक को प्राप्त करता है जो तब बनता है जब एक तत्समक आव्यूह, , कुछ स्थिरांक λ से गुणा किया जाता है, इसे उपरोक्त गुणांक आव्यूह से घटाया जाता है ताकि इसकी अभिलाक्षणिक बहुपद प्राप्त हो सके,

और इसके शून्य के लिए हल करें।

आगे सरलीकरण और आव्यूह योग प्राप्ति के मूल नियमों को लागू करना

एकल 2×2 आव्यूह के निर्धारक को ज्ञात करने के नियमों को लागू करने से निम्नलिखित प्रारंभिक द्विघात समीकरण उत्पन्न होता है,

जिसे उपरोक्त का सरल संस्करण प्राप्त करने के लिए और सरलीकृत किया जा सकता है,

अब गुणनखंड विधि को लागू करके दिए गए द्विघात समीकरण के दो मूल, और को ज्ञात करने पर प्राप्त होता है

ऊपर परिकलित किए गए मान और , A के आवश्यक आइगेनमान हैं। कुछ स्थितियों में, माना अन्य आव्यूह ओडीई, आइगेनमान समिश्र हो सकते हैं, इस स्थिति में हल प्रक्रिया के अगले चरण के साथ-साथ अंतिम रूप और हल प्रभावशाली रूप से परिवर्तित कर सकते हैं।

द्वितीय चरण

जैसा कि ऊपर उल्लेख किया गया है, इस चरण में मूल रूप से प्रदान की गई जानकारी से A के आइजनसदिशों को ज्ञात करना सम्मिलित है।

गणना किए गए प्रत्येक आइगेनमान के लिए, हमारे पास एक विशिष्ट आइगेनसदिश है। पहले आइगेनमान के लिए, जो है, हमें निम्न प्राप्त है

मूल आव्यूह गुणन नियमों को लागू करके उपरोक्त व्यंजक को सरल बनाने पर प्राप्त होता है

ये सभी गणना केवल अंतिम व्यंजक प्राप्त करने के लिए की गई है, जो हमारी स्थिति में α = 2β है। अब कुछ यादृच्छिक मान लेते हुए, संभवतः, एक छोटा महत्वहीन मान, जिसके साथ कार्य करना बहुत सरल होता है, α या β के लिए (अधिकतर स्थितियों में, यह वास्तव में कोई मायने नहीं रखता), हम इसे α = 2β में प्रतिस्थापित करते हैं। ऐसा करने से एक साधारण सदिश उत्पन्न होता है, जो इस विशेष आइगेनमान के लिए आवश्यक आइगेनसदिश होता है। हमारी स्थिति में, हम α = 2 चयनित करते हैं, जो बदले में यह निर्धारित करता है कि β = 1 और, मानक सदिश संकेतन का उपयोग करके, हमारा सदिश कुछ इस प्रकार का होता है

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

तृतीय चरण

यह अंतिम चरण उन आवश्यक फलनों को ज्ञात करता है जो मूल रूप से हमें दिए गए अवकलज के पीछे 'छिपे हुए' हैं। दो फलन हैं, क्योंकि हमारे अवकल समीकरण दो चर के साथ काम करते हैं।

वह समीकरण जिसमें वह सभी जानकारी सम्मिलित है जो हमने पहले प्राप्त की थी,  इसका का निम्न रूप है:

आइगेनमानों ​​और आइगेनसदिशों के मूल्यों को प्रतिस्थापित करने से प्राप्ति होती है

और अधिक सरलीकरण लागू करना,

और अधिक सरल करना और फलन x और y के लिए अलग-अलग समीकरण लिखना,

उपर्युक्त समीकरण, वास्तव में, वांछित सामान्य फलन हैं, लेकिन वे अपने सामान्य रूप में होते हैं (A और B के अनिर्दिष्ट मूल्यों के साथ), जबकि हम वास्तव में उनके सटीक रूप और हल प्राप्त करना चाहते हैं। तो अब हम समस्या की दी गई प्रारंभिक शर्तों पर विचार करते हैं (दिए गए प्रारंभिक शर्तों सहित समस्या तथाकथित प्रारंभिक मूल्य समस्या है)। मान लीजिए कि हमें दिया गया है, जो हमारे साधारण अवकल समीकरण के लिए शुरुआती बिंदु की भूमिका निभाता है; इन शर्तों का अनुप्रयोग स्थिरांक, A और B को निर्दिष्ट करता है। जैसा कि हम स्थितियों से देखते हैं, जब t = 0, उपरोक्त समीकरणों के बाएँ पक्ष 1 के बराबर होते हैं। इस प्रकार हम रैखिक समीकरणों की निम्नलिखित प्रणाली का निर्माण कर सकते हैं,

इन समीकरणों को हल करने पर, हम प्राप्त करते हैं कि दोनों स्थिरांक A और B 1/3 के बराबर हैं। इसलिए इन मानों को इन दो कार्यों के सामान्य रूप में प्रतिस्थापित करने से उनके सटीक रूप निर्दिष्ट होते हैं,


वांछित दो फलन।

आव्यूह घातांक का उपयोग करना

उपरोक्त समस्या को आव्यूह घातांक के अनुप्रयोग के साथ हल किया जा सकता था। अर्थात् हम ऐसा कह सकते हैं


यह देखते हुए कि (जिसे किसी भी उपयुक्त उपकरण का उपयोग करके गणना की जा सकती है, जैसे कि एमएटीएलएबी का expm टूल, या आव्यूह विकर्णकरण करके और संपत्ति का लाभ उठाकर कि विकर्ण आव्यूह का आव्यूह घातांक इसके तत्वों के तत्व-वार घातांक के समान है)


अंतिम परिणाम है

यह पहले दिखाए गए आइगेनसदिश दृष्टिकोण जैसा ही है।

यह भी देखें

संदर्भ

  1. Moya-Cessa, H.; Soto-Eguibar, F. (2011). विभेदक समीकरण: एक परिचालन दृष्टिकोण. New Jersey: Rinton Press. ISBN 978-1-58949-060-4.
  2. Putzer, E. J. (1966). "निरंतर गुणांक वाले रैखिक प्रणालियों की चर्चा में जॉर्डन कैननिकल फॉर्म से बचना". The American Mathematical Monthly. 73 (1): 2–7. doi:10.1080/00029890.1966.11970714. JSTOR 2313914.

श्रेणी:साधारण अवकल समीकरण