स्प्लिट-स्टेप विधि: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
Line 1: Line 1:
[[संख्यात्मक विश्लेषण]] में, स्प्लिट-स्टेप (फूरियर) विधि एक स्यूडो-वर्णक्रमीय संख्यात्मक विधि है जिसका उपयोग नॉनलाइनियर श्रोडिंगर समीकरण जैसे नॉनलाइनियर आंशिक अंतर समीकरणों का समाधान करने के लिए किया जाता है। यह नाम दो कारणों से उत्पन्न हुआ है। सबसे पहले, विधि छोटे चरणों में समाधान की गणना करने और रैखिक और गैर-रेखीय चरणों को अलग-अलग करने पर निर्भर करती है (नीचे देखें)। दूसरा, फूरियर को आगे और पीछे बदलना आवश्यक है क्योंकि रैखिक चरण [[आवृत्ति डोमेन]] में बनाया जाता है जबकि गैर-रेखीय चरण समय डोमेन में बनाया जाता है।
[[संख्यात्मक विश्लेषण]] में, स्प्लिट-स्टेप (फूरियर) विधि एक स्यूडो-वर्णक्रमीय संख्यात्मक विधि है जिसका उपयोग नॉनलाइनियर श्रोडिंगर समीकरण जैसे नॉनलाइनियर आंशिक अंतर समीकरणों का समाधान करने के लिए किया जाता है। यह नाम दो कारणों से उत्पन्न हुआ है। सबसे पहले, विधि छोटे चरणों में समाधान की गणना करने और रैखिक और गैर-रेखीय चरणों को अलग-अलग करने पर निर्भर करती है (नीचे देखें)। दूसरा, फूरियर को आगे और पीछे बदलना आवश्यक है क्योंकि रैखिक चरण [[आवृत्ति डोमेन]] में बनाया जाता है जबकि गैर-रेखीय चरण समय डोमेन में बनाया जाता है।


इस विधि के उपयोग का एक उदाहरण ऑप्टिकल फाइबर में प्रकाश पल्स प्रसार के क्षेत्र में है, जहां रैखिक और गैर-रेखीय तंत्र की बातचीत से सामान्य विश्लेषणात्मक समाधान ढूंढना मुश्किल हो जाता है। चूँकि, स्प्लिट-स्टेप विधि समस्या का संख्यात्मक समाधान प्रदान करती है। स्प्लिट-स्टेप विधि का एक और अनुप्रयोग जो 2010 के बाद से बहुत अधिक लोकप्रियता प्राप्त कर रहा है वह [[ऑप्टिकल माइक्रोरेसोनेटर]] में [[केर आवृत्ति कंघी]] गतिशीलता का अनुकरण है।<ref>{{Cite journal|last1=Erkintalo|first1=Miro|author-link1=Miro Erkintalo |last2=Sylvestre|first2=Thibaut|last3=Randle|first3=Hamish G.|last4=Coen|first4=Stéphane|date=2013-01-01|title=Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato–Lefever model|journal=Optics Letters|language=EN|volume=38|issue=1|pages=37–39|doi=10.1364/OL.38.000037|pmid=23282830|issn=1539-4794|arxiv=1211.1697|bibcode=2013OptL...38...37C|s2cid=7248349 }}</ref><ref>{{Cite journal|last1=Maleki|first1=L.|last2=Seidel|first2=D.|last3=Ilchenko|first3=V. S.|last4=Liang|first4=W.|last5=Savchenkov|first5=A. A.|last6=Matsko|first6=A. B.|date=2011-08-01|title=मोड-लॉक केर फ्रीक्वेंसी कॉम्ब्स|journal=Optics Letters|language=EN|volume=36|issue=15|pages=2845–2847|doi=10.1364/OL.36.002845|pmid=21808332|issn=1539-4794|bibcode=2011OptL...36.2845M}}</ref><ref>{{Cite journal|last1=Hansson|first1=Tobias|last2=Wabnitz|first2=Stefan|date=2016|title=Dynamics of microresonator frequency comb generation: models and stability|journal=Nanophotonics|volume=5|issue=2|pages=231–243|doi=10.1515/nanoph-2016-0012|issn=2192-8606|bibcode=2016Nanop...5...12H|url=https://iris.unibs.it/bitstream/11379/477683/1/nanoph-2016-0012.pdf|doi-access=free}}</ref> उचित संख्यात्मक लागत के साथ लुगियाटो-लेफ़ेवर समीकरण के कार्यान्वयन की सापेक्ष आसानी, प्रयोगात्मक स्पेक्ट्रा को पुन: प्रस्तुत करने में इसकी सफलता के साथ-साथ इन माइक्रोरेसोनेटर में [[सॉलिटन|सॉलिटॉन]] व्यवहार की भविष्यवाणी ने विधि को बहुत लोकप्रिय बना दिया है।
इस विधि के उपयोग का एक उदाहरण ऑप्टिकल फाइबर में प्रकाश पल्स प्रसार के क्षेत्र में है, जहां रैखिक और गैर-रेखीय तंत्र की बातचीत से सामान्य विश्लेषणात्मक समाधान ढूंढना मुश्किल हो जाता है। चूँकि, स्प्लिट-स्टेप विधि समस्या का संख्यात्मक समाधान प्रदान करती है। स्प्लिट-स्टेप विधि का एक और अनुप्रयोग जो 2010 के बाद से बहुत अधिक लोकप्रियता प्राप्त कर रहा है वह [[ऑप्टिकल माइक्रोरेसोनेटर]] में [[केर आवृत्ति कंघी|केर आवृत्ति काम्ब]] गतिशीलता का अनुकरण है।<ref>{{Cite journal|last1=Erkintalo|first1=Miro|author-link1=Miro Erkintalo |last2=Sylvestre|first2=Thibaut|last3=Randle|first3=Hamish G.|last4=Coen|first4=Stéphane|date=2013-01-01|title=Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato–Lefever model|journal=Optics Letters|language=EN|volume=38|issue=1|pages=37–39|doi=10.1364/OL.38.000037|pmid=23282830|issn=1539-4794|arxiv=1211.1697|bibcode=2013OptL...38...37C|s2cid=7248349 }}</ref><ref>{{Cite journal|last1=Maleki|first1=L.|last2=Seidel|first2=D.|last3=Ilchenko|first3=V. S.|last4=Liang|first4=W.|last5=Savchenkov|first5=A. A.|last6=Matsko|first6=A. B.|date=2011-08-01|title=मोड-लॉक केर फ्रीक्वेंसी कॉम्ब्स|journal=Optics Letters|language=EN|volume=36|issue=15|pages=2845–2847|doi=10.1364/OL.36.002845|pmid=21808332|issn=1539-4794|bibcode=2011OptL...36.2845M}}</ref><ref>{{Cite journal|last1=Hansson|first1=Tobias|last2=Wabnitz|first2=Stefan|date=2016|title=Dynamics of microresonator frequency comb generation: models and stability|journal=Nanophotonics|volume=5|issue=2|pages=231–243|doi=10.1515/nanoph-2016-0012|issn=2192-8606|bibcode=2016Nanop...5...12H|url=https://iris.unibs.it/bitstream/11379/477683/1/nanoph-2016-0012.pdf|doi-access=free}}</ref> उचित संख्यात्मक लागत के साथ लुगियाटो-लेफ़ेवर समीकरण के कार्यान्वयन की सापेक्ष आसानी, प्रयोगात्मक स्पेक्ट्रा को पुन: प्रस्तुत करने में इसकी सफलता के साथ-साथ इन माइक्रोरेसोनेटर में [[सॉलिटन|सॉलिटॉन]] व्यवहार की भविष्यवाणी ने विधि को बहुत लोकप्रिय बना दिया है।


==विधि का विवरण==
==विधि का विवरण==
उदाहरण के लिए, अरेखीय श्रोडिंगर समीकरण पर विचार करें<ref name="NonlinearAgrawal">{{cite book |last=Agrawal |first=Govind P. |title=नॉनलीनियर फाइबर ऑप्टिक्स|edition=3rd|year=2001 |publisher=Academic Press |location=San Diego, CA, USA|isbn=0-12-045143-3}}</ref>
उदाहरण के लिए, अरेखीय श्रोडिंगर समीकरण पर विचार करें<ref name="NonlinearAgrawal">{{cite book |last=Agrawal |first=Govind P. |title=नॉनलीनियर फाइबर ऑप्टिक्स|edition=3rd|year=2001 |publisher=Academic Press |location=San Diego, CA, USA|isbn=0-12-045143-3}}</ref>
:<math>{\partial A \over \partial z} = - {i\beta_2 \over 2} {\partial^2 A \over \partial t^2} + i \gamma | A |^2 A = [\hat D + \hat N]A, </math>
:<math>{\partial A \over \partial z} = - {i\beta_2 \over 2} {\partial^2 A \over \partial t^2} + i \gamma | A |^2 A = [\hat D + \hat N]A, </math>
कहाँ <math>A(t,z)</math> समय में नाड़ी आवरण का वर्णन करता है <math>t</math> स्थानिक स्थिति पर <math>z</math>. समीकरण को रैखिक भाग में विभाजित किया जा सकता है,
जहां <math>A(t,z)</math> स्थानिक स्थिति <math>z</math> पर समय <math>t</math> में पल्स लिफाफे का वर्णन करता है। समीकरण को एक रैखिक भाग में विभाजित किया जा सकता है,
:<math>{\partial A_D \over \partial z} = - {i\beta_2 \over 2} {\partial^2 A \over \partial t^2} = \hat D A, </math>
:<math>{\partial A_D \over \partial z} = - {i\beta_2 \over 2} {\partial^2 A \over \partial t^2} = \hat D A, </math>
और अरैखिक भाग,
और अरैखिक भाग,
:<math>{\partial A_N \over \partial z} = i \gamma | A |^2 A = \hat N A. </math>
:<math>{\partial A_N \over \partial z} = i \gamma | A |^2 A = \hat N A. </math>
रैखिक और अरेखीय दोनों भागों में विश्लेषणात्मक समाधान होते हैं, लेकिन दोनों भागों वाले अरेखीय श्रोडिंगर समीकरण में कोई सामान्य विश्लेषणात्मक समाधान नहीं होता है।
रैखिक और अरेखीय दोनों भागों में विश्लेषणात्मक समाधान होते हैं, किन्तु दोनों भागों वाले अरेखीय श्रोडिंगर समीकरण में कोई सामान्य विश्लेषणात्मक समाधान नहीं होता है।


चूँकि, यदि केवल 'छोटा' कदम है <math>h</math> साथ ले जाया जाता है <math>z</math>, तो केवल 'छोटी' संख्यात्मक त्रुटि के साथ दोनों भागों को अलग-अलग माना जा सकता है। इसलिए सबसे पहले कोई छोटा सा अरैखिक कदम उठा सकता है,
चूँकि, यदि <math>z</math> के साथ केवल एक 'छोटा' कदम <math>h</math> उठाया जाता है, तो केवल 'छोटी' संख्यात्मक त्रुटि के साथ दोनों भागों को अलग-अलग माना जा सकता है। इसलिए कोई भी पहले विश्लेषणात्मक समाधान का उपयोग करके एक छोटा गैर-रैखिक कदम


:<math>A_N(t, z+h) = \exp\left[i \gamma |A(t, z)|^2 h \right] A(t, z), </math>
:<math>A_N(t, z+h) = \exp\left[i \gamma |A(t, z)|^2 h \right] A(t, z), </math>
विश्लेषणात्मक समाधान का उपयोग करना. ध्यान दें कि यह ansatz लगाता है <math>|A(z)|^2=const.</math> और इसके परिणामस्वरूप <math>\gamma \in \mathbb{R}</math>.
ले सकता है। ध्यान दें कि यह अंसत्ज़ <math>|A(z)|^2=const</math> लगाता है और इसके परिणामस्वरूप <math>\gamma \in \mathbb{R}</math> लगाता है।


फैलाव चरण में आवृत्ति डोमेन में विश्लेषणात्मक समाधान होता है, इसलिए फूरियर रूपांतरण के लिए यह सबसे पहले आवश्यक है <math>A_N</math> का उपयोग करते हुए
प्रसार चरण में आवृत्ति डोमेन में विश्लेषणात्मक समाधान होता है, इसलिए फूरियर को <math>A_N</math> का उपयोग करके रूपांतरित करना सबसे पहले आवश्यक है
:<math>\tilde A_N(\omega, z) = \int_{-\infty}^\infty A_N(t,z) \exp[i(\omega-\omega_0)t] dt </math>,
:<math>\tilde A_N(\omega, z) = \int_{-\infty}^\infty A_N(t,z) \exp[i(\omega-\omega_0)t] dt </math>,
कहाँ <math>\omega_0</math> नाड़ी की केंद्र आवृत्ति है.
जहाँ <math>\omega_0</math> नाड़ी की केंद्र आवृत्ति है।
 
यह दिखाया जा सकता है कि फूरियर ट्रांसफॉर्म की उपरोक्त परिभाषा का उपयोग करके, रैखिक चरण का विश्लेषणात्मक समाधान, गैर-रेखीय चरण के लिए आवृत्ति डोमेन समाधान के साथ परिवर्तित किया जाता है।
यह दिखाया जा सकता है कि फूरियर ट्रांसफॉर्म की उपरोक्त परिभाषा का उपयोग करके, रैखिक चरण का विश्लेषणात्मक समाधान, गैर-रेखीय चरण के लिए आवृत्ति डोमेन समाधान के साथ परिवर्तित किया जाता है।


:<math>\tilde{A}(\omega, z+h) = \exp\left[{i \beta_2 \over 2} (\omega-\omega_0)^2 h \right] \tilde{A}_N(\omega, z).</math>
:<math>\tilde{A}(\omega, z+h) = \exp\left[{i \beta_2 \over 2} (\omega-\omega_0)^2 h \right] \tilde{A}_N(\omega, z).</math>
का व्युत्क्रम फूरियर रूपांतरण लेकर <math>\tilde{A}(\omega, z+h)</math> प्राप्त होता है <math>A\left(t, z+h\right)</math>; इस प्रकार नाड़ी को छोटे कदम से प्रचारित किया गया है <math>h</math>. उपरोक्त को दोहराते हुए <math>N</math> कई बार, नाड़ी को लंबाई तक प्रसारित किया जा सकता है <math>N h</math>.
<math>\tilde{A}(\omega, z+h)</math> का व्युत्क्रम फूरियर रूपांतरण लेने से <math>A\left(t, z+h\right)</math> प्राप्त होता है; इस प्रकार पल्स को एक छोटे चरण <math>h</math> में प्रसारित किया गया है। उपरोक्त <math>N</math> बार दोहराकर, पल्स को <math>N h</math> की लंबाई में प्रसारित किया जा सकता है।


ऊपर दिखाया गया है कि किसी समाधान को अंतरिक्ष में आगे बढ़ाने के लिए विधि का उपयोग कैसे किया जाए; चूँकि, कई भौतिकी अनुप्रयोगों, जैसे कि कण का वर्णन करने वाले तरंग पैकेट के विकास का अध्ययन, के लिए अंतरिक्ष के बजाय समय में समाधान को आगे बढ़ाने की आवश्यकता होती है। गैर-रैखिक श्रोडिंगर समीकरण, जब तरंग फ़ंक्शन के समय विकास को नियंत्रित करने के लिए उपयोग किया जाता है, तो रूप लेता है
ऊपर दिखाया गया है कि किसी समाधान को अंतरिक्ष में आगे बढ़ाने के लिए विधि का उपयोग कैसे किया जाए; चूँकि, कई भौतिकी अनुप्रयोगों, जैसे कि कण का वर्णन करने वाले तरंग पैकेट के विकास का अध्ययन, के लिए अंतरिक्ष के अतिरिक्त समय में समाधान को आगे बढ़ाने की आवश्यकता होती है। गैर-रैखिक श्रोडिंगर समीकरण, जब तरंग फ़ंक्शन के समय विकास को नियंत्रित करने के लिए उपयोग किया जाता है, तो रूप लेता है


:<math>i \hbar {\partial \psi \over \partial t} = - {{\hbar}^2 \over {2m}} {\partial^2 \psi \over \partial x^2} +  \gamma | \psi|^2 \psi = [\hat D + \hat N]\psi, </math>
:<math>i \hbar {\partial \psi \over \partial t} = - {{\hbar}^2 \over {2m}} {\partial^2 \psi \over \partial x^2} +  \gamma | \psi|^2 \psi = [\hat D + \hat N]\psi, </math>
कहाँ <math>\psi(x, t)</math> स्थिति में तरंग फ़ंक्शन का वर्णन करता है <math>x</math> और समय <math>t</math>. ध्यान दें कि
जहां <math>\psi(x, t)</math> स्थिति <math>x</math> और समय <math>t</math> पर तरंग फ़ंक्शन का वर्णन करता है। ध्यान दें कि
:<math>\hat D=- {{\hbar}^2 \over {2m}} {\partial^2 \over \partial x^2}</math> और <math> \hat N =\gamma | \psi|^2 </math>, ओर वो <math> m </math> कण का द्रव्यमान है और <math> \hbar </math> प्लैंक का स्थिरांक है <math>2\pi</math>.
:<math>\hat D=- {{\hbar}^2 \over {2m}} {\partial^2 \over \partial x^2}</math> और <math> \hat N =\gamma | \psi|^2 </math>, और वह <math> m </math> कण का द्रव्यमान है और <math> \hbar </math> <math>2\pi</math> से अधिक प्लैंक स्थिरांक है।
इस समीकरण का औपचारिक समाधान जटिल घातांक है, इसलिए हमारे पास वह है
इस समीकरण का औपचारिक समाधान जटिल घातांक है, इसलिए हमारे पास वह है
:<math> \psi(x, t)=e^{-it(\hat D+\hat N)/\hbar}\psi(x, 0)</math>.
:<math> \psi(x, t)=e^{-it(\hat D+\hat N)/\hbar}\psi(x, 0)</math>.
तब से <math>\hat{D}</math> और <math>\hat{N}</math> ऑपरेटर हैं, वे सामान्य रूप से आवागमन नहीं करते हैं। चूँकि, बेकर-हॉसडॉर्फ फॉर्मूला यह दिखाने के लिए लागू किया जा सकता है कि उनके साथ ऐसा व्यवहार करने से त्रुटि व्यवस्थित होगी <math>dt^2</math> यदि हम छोटा लेकिन सीमित समय वाला कदम उठा रहे हैं <math>dt</math>. इसलिए हम लिख सकते हैं
तब से <math>\hat{D}</math> और <math>\hat{N}</math> ऑपरेटर हैं, वे सामान्य रूप से आवागमन नहीं करते हैं। चूँकि, बेकर-हॉसडॉर्फ फॉर्मूला यह दिखाने के लिए लागू किया जा सकता है कि यदि हम एक छोटा लेकिन सीमित समय कदम <math>dt</math> ले रहे हैं, तो उन्हें इस तरह मानने से त्रुटि <math>dt^2</math> क्रम की होगी। इसलिए हम लिख सकते हैं
:<math> \psi(x, t+dt) \approx e^{-idt\hat D/\hbar}e^{-idt\hat N/\hbar}\psi(x, t)</math>.
:<math> \psi(x, t+dt) \approx e^{-idt\hat D/\hbar}e^{-idt\hat N/\hbar}\psi(x, t)</math>.
इस समीकरण का वह भाग जिसमें शामिल है <math> \hat N </math> समय पर तरंग फ़ंक्शन का उपयोग करके सीधे गणना की जा सकती है <math> t </math>, लेकिन शामिल घातांक की गणना करने के लिए <math> \hat D </math> हम इस तथ्य का उपयोग करते हैं कि आवृत्ति स्थान में, आंशिक व्युत्पन्न ऑपरेटर को प्रतिस्थापित करके संख्या में परिवर्तित किया जा सकता है <math> ik </math> के लिए <math> \partial \over \partial x </math>, कहाँ <math> k</math> आवृत्ति (या अधिक ठीक से, तरंग संख्या, जैसा कि हम स्थानिक चर के साथ काम कर रहे हैं और इस प्रकार स्थानिक आवृत्तियों के स्थान में परिवर्तित हो रहे हैं - यानी तरंग संख्या) जो कुछ भी संचालित किया जा रहा है उसके फूरियर रूपांतरण से जुड़ा हुआ है। इस प्रकार, हम फूरियर रूपांतरण लेते हैं
इस समीकरण का वह भाग जिसमें शामिल है <math> \hat N </math> समय पर तरंग फ़ंक्शन का उपयोग करके सीधे गणना की जा सकती है <math> t </math>, किन्तु शामिल घातांक की गणना करने के लिए <math> \hat D </math> हम इस तथ्य का उपयोग करते हैं कि आवृत्ति स्थान में, आंशिक व्युत्पन्न ऑपरेटर को प्रतिस्थापित करके संख्या में परिवर्तित किया जा सकता है <math> ik </math> के लिए <math> \partial \over \partial x </math>, जहाँ <math> k</math> आवृत्ति (या अधिक ठीक से, तरंग संख्या, जैसा कि हम स्थानिक चर के साथ काम कर रहे हैं और इस प्रकार स्थानिक आवृत्तियों के स्थान में परिवर्तित हो रहे हैं - यानी तरंग संख्या) जो कुछ भी संचालित किया जा रहा है उसके फूरियर रूपांतरण से जुड़ा हुआ है। इस प्रकार, हम फूरियर रूपांतरण लेते हैं
:<math>e^{-idt\hat N/\hbar}\psi(x, t)</math>,
:<math>e^{-idt\hat N/\hbar}\psi(x, t)</math>,
संबंधित तरंग संख्या पुनर्प्राप्त करें, मात्रा की गणना करें
संबंधित तरंग संख्या पुनर्प्राप्त करें, मात्रा की गणना करें
Line 40: Line 41:
और इसमें शामिल जटिल घातांकों का गुणनफल खोजने के लिए इसका उपयोग करें <math> \hat N</math> और <math> \hat D </math> आवृत्ति स्थान में निम्नानुसार:
और इसमें शामिल जटिल घातांकों का गुणनफल खोजने के लिए इसका उपयोग करें <math> \hat N</math> और <math> \hat D </math> आवृत्ति स्थान में निम्नानुसार:
:<math> e^{idtk^2}F[e^{-idt\hat N}\psi(x, t)]</math>,
:<math> e^{idtk^2}F[e^{-idt\hat N}\psi(x, t)]</math>,
कहाँ <math> F</math> फूरियर रूपांतरण को दर्शाता है। फिर हम व्युत्क्रम फूरियर इस अभिव्यक्ति को भौतिक स्थान में अंतिम परिणाम खोजने के लिए रूपांतरित करते हैं, जिससे अंतिम अभिव्यक्ति प्राप्त होती है
जहाँ <math> F</math> फूरियर रूपांतरण को दर्शाता है। फिर हम व्युत्क्रम फूरियर इस अभिव्यक्ति को भौतिक स्थान में अंतिम परिणाम खोजने के लिए रूपांतरित करते हैं, जिससे अंतिम अभिव्यक्ति प्राप्त होती है
:<math>\psi(x, t+dt)=F^{-1}[e^{idtk^2}F[e^{-idt\hat N}\psi(x, t)]]</math>.
:<math>\psi(x, t+dt)=F^{-1}[e^{idtk^2}F[e^{-idt\hat N}\psi(x, t)]]</math>.
इस पद्धति का रूपांतर सममितीय विभाजन-चरण फूरियर विधि है, जो ऑपरेटर का उपयोग करके आधा समय कदम उठाती है, फिर केवल दूसरे के साथ पूर्णकालिक कदम उठाती है, और फिर केवल पहले के साथ फिर से दूसरा आधा समय कदम उठाती है। यह विधि सामान्य स्प्लिट-स्टेप फूरियर विधि का सुधार है क्योंकि इसकी त्रुटि क्रमानुसार है <math>dt^3</math> समय के कदम के लिए <math>dt</math>.
इस पद्धति का रूपांतर सममितीय विभाजन-चरण फूरियर विधि है, जो ऑपरेटर का उपयोग करके आधा समय कदम उठाती है, फिर केवल दूसरे के साथ पूर्णकालिक कदम उठाती है, और फिर केवल पहले के साथ फिर से दूसरा आधा समय कदम उठाती है। यह विधि सामान्य स्प्लिट-स्टेप फूरियर विधि का सुधार है क्योंकि इसकी त्रुटि क्रमानुसार है <math>dt^3</math> समय के कदम के लिए <math>dt</math>.

Revision as of 07:50, 30 July 2023

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

इस विधि के उपयोग का एक उदाहरण ऑप्टिकल फाइबर में प्रकाश पल्स प्रसार के क्षेत्र में है, जहां रैखिक और गैर-रेखीय तंत्र की बातचीत से सामान्य विश्लेषणात्मक समाधान ढूंढना मुश्किल हो जाता है। चूँकि, स्प्लिट-स्टेप विधि समस्या का संख्यात्मक समाधान प्रदान करती है। स्प्लिट-स्टेप विधि का एक और अनुप्रयोग जो 2010 के बाद से बहुत अधिक लोकप्रियता प्राप्त कर रहा है वह ऑप्टिकल माइक्रोरेसोनेटर में केर आवृत्ति काम्ब गतिशीलता का अनुकरण है।[1][2][3] उचित संख्यात्मक लागत के साथ लुगियाटो-लेफ़ेवर समीकरण के कार्यान्वयन की सापेक्ष आसानी, प्रयोगात्मक स्पेक्ट्रा को पुन: प्रस्तुत करने में इसकी सफलता के साथ-साथ इन माइक्रोरेसोनेटर में सॉलिटॉन व्यवहार की भविष्यवाणी ने विधि को बहुत लोकप्रिय बना दिया है।

विधि का विवरण

उदाहरण के लिए, अरेखीय श्रोडिंगर समीकरण पर विचार करें[4]

जहां स्थानिक स्थिति पर समय में पल्स लिफाफे का वर्णन करता है। समीकरण को एक रैखिक भाग में विभाजित किया जा सकता है,

और अरैखिक भाग,

रैखिक और अरेखीय दोनों भागों में विश्लेषणात्मक समाधान होते हैं, किन्तु दोनों भागों वाले अरेखीय श्रोडिंगर समीकरण में कोई सामान्य विश्लेषणात्मक समाधान नहीं होता है।

चूँकि, यदि के साथ केवल एक 'छोटा' कदम उठाया जाता है, तो केवल 'छोटी' संख्यात्मक त्रुटि के साथ दोनों भागों को अलग-अलग माना जा सकता है। इसलिए कोई भी पहले विश्लेषणात्मक समाधान का उपयोग करके एक छोटा गैर-रैखिक कदम

ले सकता है। ध्यान दें कि यह अंसत्ज़ लगाता है और इसके परिणामस्वरूप लगाता है।

प्रसार चरण में आवृत्ति डोमेन में विश्लेषणात्मक समाधान होता है, इसलिए फूरियर को का उपयोग करके रूपांतरित करना सबसे पहले आवश्यक है

,

जहाँ नाड़ी की केंद्र आवृत्ति है।

यह दिखाया जा सकता है कि फूरियर ट्रांसफॉर्म की उपरोक्त परिभाषा का उपयोग करके, रैखिक चरण का विश्लेषणात्मक समाधान, गैर-रेखीय चरण के लिए आवृत्ति डोमेन समाधान के साथ परिवर्तित किया जाता है।

का व्युत्क्रम फूरियर रूपांतरण लेने से प्राप्त होता है; इस प्रकार पल्स को एक छोटे चरण में प्रसारित किया गया है। उपरोक्त बार दोहराकर, पल्स को की लंबाई में प्रसारित किया जा सकता है।

ऊपर दिखाया गया है कि किसी समाधान को अंतरिक्ष में आगे बढ़ाने के लिए विधि का उपयोग कैसे किया जाए; चूँकि, कई भौतिकी अनुप्रयोगों, जैसे कि कण का वर्णन करने वाले तरंग पैकेट के विकास का अध्ययन, के लिए अंतरिक्ष के अतिरिक्त समय में समाधान को आगे बढ़ाने की आवश्यकता होती है। गैर-रैखिक श्रोडिंगर समीकरण, जब तरंग फ़ंक्शन के समय विकास को नियंत्रित करने के लिए उपयोग किया जाता है, तो रूप लेता है

जहां स्थिति और समय पर तरंग फ़ंक्शन का वर्णन करता है। ध्यान दें कि

और , और वह कण का द्रव्यमान है और से अधिक प्लैंक स्थिरांक है।

इस समीकरण का औपचारिक समाधान जटिल घातांक है, इसलिए हमारे पास वह है

.

तब से और ऑपरेटर हैं, वे सामान्य रूप से आवागमन नहीं करते हैं। चूँकि, बेकर-हॉसडॉर्फ फॉर्मूला यह दिखाने के लिए लागू किया जा सकता है कि यदि हम एक छोटा लेकिन सीमित समय कदम ले रहे हैं, तो उन्हें इस तरह मानने से त्रुटि क्रम की होगी। इसलिए हम लिख सकते हैं

.

इस समीकरण का वह भाग जिसमें शामिल है समय पर तरंग फ़ंक्शन का उपयोग करके सीधे गणना की जा सकती है , किन्तु शामिल घातांक की गणना करने के लिए हम इस तथ्य का उपयोग करते हैं कि आवृत्ति स्थान में, आंशिक व्युत्पन्न ऑपरेटर को प्रतिस्थापित करके संख्या में परिवर्तित किया जा सकता है के लिए , जहाँ आवृत्ति (या अधिक ठीक से, तरंग संख्या, जैसा कि हम स्थानिक चर के साथ काम कर रहे हैं और इस प्रकार स्थानिक आवृत्तियों के स्थान में परिवर्तित हो रहे हैं - यानी तरंग संख्या) जो कुछ भी संचालित किया जा रहा है उसके फूरियर रूपांतरण से जुड़ा हुआ है। इस प्रकार, हम फूरियर रूपांतरण लेते हैं

,

संबंधित तरंग संख्या पुनर्प्राप्त करें, मात्रा की गणना करें

,

और इसमें शामिल जटिल घातांकों का गुणनफल खोजने के लिए इसका उपयोग करें और आवृत्ति स्थान में निम्नानुसार:

,

जहाँ फूरियर रूपांतरण को दर्शाता है। फिर हम व्युत्क्रम फूरियर इस अभिव्यक्ति को भौतिक स्थान में अंतिम परिणाम खोजने के लिए रूपांतरित करते हैं, जिससे अंतिम अभिव्यक्ति प्राप्त होती है

.

इस पद्धति का रूपांतर सममितीय विभाजन-चरण फूरियर विधि है, जो ऑपरेटर का उपयोग करके आधा समय कदम उठाती है, फिर केवल दूसरे के साथ पूर्णकालिक कदम उठाती है, और फिर केवल पहले के साथ फिर से दूसरा आधा समय कदम उठाती है। यह विधि सामान्य स्प्लिट-स्टेप फूरियर विधि का सुधार है क्योंकि इसकी त्रुटि क्रमानुसार है समय के कदम के लिए . इस कलन विधि के फूरियर ट्रांसफॉर्म की गणना फास्ट फूरियर ट्रांसफॉर्म | फास्ट फूरियर ट्रांसफॉर्म (एफएफटी) का उपयोग करके अपेक्षाकृत तेजी से की जा सकती है। इसलिए स्प्लिट-स्टेप फूरियर विधि विशिष्ट परिमित अंतर विधियों की तुलना में बहुत तेज़ हो सकती है।[5]


संदर्भ

  1. Erkintalo, Miro; Sylvestre, Thibaut; Randle, Hamish G.; Coen, Stéphane (2013-01-01). "Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato–Lefever model". Optics Letters (in English). 38 (1): 37–39. arXiv:1211.1697. Bibcode:2013OptL...38...37C. doi:10.1364/OL.38.000037. ISSN 1539-4794. PMID 23282830. S2CID 7248349.
  2. Maleki, L.; Seidel, D.; Ilchenko, V. S.; Liang, W.; Savchenkov, A. A.; Matsko, A. B. (2011-08-01). "मोड-लॉक केर फ्रीक्वेंसी कॉम्ब्स". Optics Letters (in English). 36 (15): 2845–2847. Bibcode:2011OptL...36.2845M. doi:10.1364/OL.36.002845. ISSN 1539-4794. PMID 21808332.
  3. Hansson, Tobias; Wabnitz, Stefan (2016). "Dynamics of microresonator frequency comb generation: models and stability" (PDF). Nanophotonics. 5 (2): 231–243. Bibcode:2016Nanop...5...12H. doi:10.1515/nanoph-2016-0012. ISSN 2192-8606.
  4. Agrawal, Govind P. (2001). नॉनलीनियर फाइबर ऑप्टिक्स (3rd ed.). San Diego, CA, USA: Academic Press. ISBN 0-12-045143-3.
  5. T. R. Taha and M. J. Ablowitz (1984). "Analytical and numerical aspects of certain nonlinear evolution equations. II. Numerical, nonlinear Schrödinger equation". J. Comput. Phys. 55 (2): 203–230. Bibcode:1984JCoPh..55..203T. doi:10.1016/0021-9991(84)90003-2.


बाहरी सन्दर्भ

श्रेणी:संख्यात्मक अंतर समीकरण श्रेणी:फाइबर ऑप्टिक्स