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

From Vigyanwiki
No edit summary
No edit summary
 
(6 intermediate revisions by 3 users not shown)
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> \partial \over \partial x </math> के लिए <math> ik </math> को प्रतिस्थापित करके एक संख्या में परिवर्तित किया जा सकता है। जहां <math> k</math> आवृत्ति (या अधिक ठीक से, तरंग संख्या है, क्योंकि हम एक स्थानिक चर के साथ काम कर रहे हैं और इस प्रकार स्थानिक आवृत्तियों के स्थान में परिवर्तित हो रहे हैं - अर्थात् तरंग संख्या) जो कुछ भी संचालित किया जा रहा है उसके फूरियर रूपांतरण से जुड़ा हुआ है। इस प्रकार, हम फूरियर रूपांतरण लेते हैं
:<math>e^{-idt\hat N/\hbar}\psi(x, t)</math>,
:<math>e^{-idt\hat N/\hbar}\psi(x, t)</math>,
संबंधित तरंग संख्या पुनर्प्राप्त करें, मात्रा की गणना करें
संबंधित तरंग संख्या पुनर्प्राप्त करें, मात्रा की गणना करें
:<math> e^{idtk^2}</math>,
:<math> e^{idtk^2}</math>,
और इसमें शामिल जटिल घातांकों का गुणनफल खोजने के लिए इसका उपयोग करें <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</math> के लिए <math>dt^3</math> क्रम की है। इस [[कलन विधि|एल्गोरिदम]] के फूरियर ट्रांसफॉर्म की गणना [[फास्ट फूरियर ट्रांसफॉर्म]] (एफएफटी) का उपयोग करके अपेक्षाकृत तेजी से की जा सकती है। इसलिए स्प्लिट-स्टेप फूरियर विधि विशिष्ट [[परिमित अंतर विधि|परिमित अंतर विधियों]] की तुलना में बहुत तेज़ हो सकती है।<ref name="Taha1984">{{cite journal
इस [[कलन विधि]] के फूरियर ट्रांसफॉर्म की गणना [[फास्ट फूरियर ट्रांसफॉर्म]] | फास्ट फूरियर ट्रांसफॉर्म (एफएफटी) का उपयोग करके अपेक्षाकृत तेजी से की जा सकती है। इसलिए स्प्लिट-स्टेप फूरियर विधि विशिष्ट [[परिमित अंतर विधि]]यों की तुलना में बहुत तेज़ हो सकती है।<ref name="Taha1984">{{cite journal
  | author = T. R. Taha and M. J. Ablowitz
  | author = T. R. Taha and M. J. Ablowitz
  | year = 1984
  | year = 1984
Line 69: Line 69:
श्रेणी:फाइबर ऑप्टिक्स
श्रेणी:फाइबर ऑप्टिक्स


 
[[Category:CS1 English-language sources (en)]]
[[Category: Machine Translated Page]]
[[Category:Collapse templates]]
[[Category:Created On 24/07/2023]]
[[Category:Created On 24/07/2023]]
[[Category:Machine Translated Page]]
[[Category:Navigational boxes| ]]
[[Category:Navigational boxes without horizontal lists]]
[[Category:Pages with script errors]]
[[Category:Sidebars with styles needing conversion]]
[[Category:Template documentation pages|Documentation/doc]]
[[Category:Templates Vigyan Ready]]
[[Category:Templates generating microformats]]
[[Category:Templates that are not mobile friendly]]
[[Category:Templates using TemplateData]]
[[Category:Wikipedia metatemplates]]

Latest revision as of 14:06, 14 August 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.


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

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