वर्णक्रमीय विधि: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
 
(8 intermediate revisions by 3 users not shown)
Line 13: Line 13:
===एक ठोस, रैखिक उदाहरण===
===एक ठोस, रैखिक उदाहरण===


यहां हम बुनियादी बहुभिन्नरूपी [[ गणना ]] और फूरियर श्रृंखला की समझ का अनुमान लगाते हैं। अगर <math>g(x,y)</math> दो वास्तविक चरों का एक ज्ञात, जटिल-मूल्यवान फ़ंक्शन है, और g, x और y में आवधिक है (अर्थात्, <math>g(x,y)=g(x+2\pi,y)=g(x,y+2\pi)</math>) तो हम एक फ़ंक्शन f(x,y) खोजने में रुचि रखते हैं ताकि
यहां हम आधारभूत बहुभिन्नरूपी [[ गणना | कलन(कैल्कुलस)]] और फूरियर श्रृंखला की समझ का अनुमान लगाते हैं। यदि <math>g(x,y)</math> दो वास्तविक चरों का एक ज्ञात, जटिल-मान फलन है, और g, x और y में आवधिक है (अर्थात्, <math>g(x,y)=g(x+2\pi,y)=g(x,y+2\pi)</math>) तो हम एक फलन f(x,y) खोजने में रुचि रखते हैं जिससे


:<math>\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right)f(x,y)=g(x,y)\quad \text{for all } x,y</math>
:<math>\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right)f(x,y)=g(x,y)\quad \text{for all } x,y</math>
<!--math>f_{xx}(x,y)+f_{yy}(x,y)=g(x,y)\quad \text{for all} x,y</math-->
जहां बाईं ओर की अभिव्यक्ति क्रमशः x और y में f के दूसरे आंशिक व्युत्पन्न को दर्शाती है। यह [[पॉइसन समीकरण]] है, और इसे भौतिक रूप से किसी प्रकार की ऊष्मा चालन समस्या, या अन्य संभावनाओं के बीच संभावित सिद्धांत में एक समस्या के रूप में व्याख्यायित किया जा सकता है।
जहां बाईं ओर की अभिव्यक्ति क्रमशः x और y में f के दूसरे आंशिक व्युत्पन्न को दर्शाती है। यह [[पॉइसन समीकरण]] है, और इसे भौतिक रूप से किसी प्रकार की ऊष्मा चालन समस्या, या अन्य संभावनाओं के बीच संभावित सिद्धांत में एक समस्या के रूप में व्याख्या किया जा सकता है।


यदि हम फूरियर श्रृंखला में f और g लिखते हैं:
यदि हम फूरियर श्रृंखला में f और g लिखते हैं:
Line 30: Line 29:
{{NumBlk|:| <math>a_{j,k}=-\frac{b_{j,k}}{j^2+k^2}</math>|{{EquationRef|<nowiki>*</nowiki>}}}}
{{NumBlk|:| <math>a_{j,k}=-\frac{b_{j,k}}{j^2+k^2}</math>|{{EquationRef|<nowiki>*</nowiki>}}}}


जो फूरियर गुणांक के लिए एक स्पष्ट सूत्र है<sub>''j'',''k''</sub>.
जो फूरियर गुणांक <sub>''j'',''k''</sub>. के लिए एक स्पष्ट सूत्र है


आवधिक सीमा स्थितियों के साथ, पॉइसन समीकरण का कोई समाधान केवल तभी होता है जब बी<sub>0,0</sub> = 0. इसलिए, हम स्वतंत्र रूप से a चुन सकते हैं<sub>0,0</sub> जो संकल्प के माध्य के बराबर होगा। यह एकीकरण स्थिरांक को चुनने के अनुरूप है।
आवधिक सीमा स्थितियों के साथ, पॉइसन समीकरण का कोई समाधान केवल तभी होता है जब b<sub>0,0</sub> = 0 होता है। इसलिए, हम स्वतंत्र रूप से a<sub>0,0</sub> चुन सकते हैं जो विश्लेषण के माध्य के बराबर होगा। यह एकीकरण स्थिरांक को चुनने के अनुरूप है।


इसे एक एल्गोरिदम में बदलने के लिए, केवल सीमित आवृत्तियों को हल किया जाता है। यह एक त्रुटि प्रस्तुत करता है जिसे आनुपातिक दिखाया जा सकता है <math>h^n</math>, कहाँ <math>h := 1/n</math> और <math>n</math> उपचारित उच्चतम आवृत्ति है।
इसे एक विधिकलन में परिवर्तित करने के लिए, केवल परिमित आवृत्तियों को हल किया जाता है। यह एक त्रुटि प्रस्तुत करता है जिसे <math>h^n</math> के आनुपातिक दिखाया जा सकता है , जहाँ <math>h := 1/n</math> और <math>n</math> उपचारित उच्चतम आवृत्ति है।


==== एल्गोरिथम ====
==== विधिकलन ====


# फूरियर रूपांतरण की गणना करें (बी<sub>j,k</sub>) जी का.
# ''g के'' फूरियर रूपांतरण (''b<sub>j,k</sub>'') की गणना करें।
# फूरियर रूपांतरण की गणना करें (a<sub>j,k</sub>) का f सूत्र के माध्यम से ({{EquationNote|*}}).
# ('''*''') सूत्र के माध्यम से f के फूरियर रूपांतरण (a<sub>j,k</sub>) की गणना करें। .
# (ए) का व्युत्क्रम फूरियर रूपांतरण लेकर एफ की गणना करें<sub>j,k</sub>).
# (''a<sub>j,k</sub>'') का व्युत्क्रम फूरियर रूपांतरण लेकर f की गणना करें


चूँकि हम केवल आवृत्तियों की एक सीमित विंडो (जैसे आकार n,) में रुचि रखते हैं, यह एक [[फास्ट फूरियर ट्रांसफॉर्म]] एल्गोरिदम का उपयोग करके किया जा सकता है। इसलिए, विश्व स्तर पर एल्गोरिथ्म चलता है {{nowrap|time ''O''(''n'' log ''n'').}}
चूँकि हम केवल आवृत्तियों की एक सीमित क्षेत्र (जैसे आकार n,) में रुचि रखते हैं, यह एक [[फास्ट फूरियर ट्रांसफॉर्म]] विधिकलन का उपयोग करके किया जा सकता है। इसलिए, विश्व स्तर पर विधिकलन {{nowrap|time ''O''(''n'' log ''n'').}} समय में चलता है


===अरेखीय उदाहरण===
===अरेखीय उदाहरण===


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


दिया गया <math>u(x,0)</math> आवधिक डोमेन पर
पुर्तस्य क्षेत्र <math>x\in\left[0,2\pi\right)</math> पर <math>t>0</math> के लिए दिए गए <math>u(x,0)</math> के साथ ऐसे <math>u \in \mathcal{U}</math> का पता लगाएं जो निम्नलिखित समीकरण को पूरा करते हैं:
<math>x\in\left[0,2\pi\right)</math>, पाना <math>u \in \mathcal{U}</math> ऐसा है कि
 
:<math>\partial_{t} u + u \partial_{x} u = \rho \partial_{xx} u + f \quad \forall x\in\left[0,2\pi\right), \forall t>0</math>
<math>\partial_{t} u + u \partial_{x} u = \rho \partial_{xx} u + f \quad \forall x\in\left[0,2\pi\right), \forall t>0</math>
जहाँ ρ [[श्यानता]] गुणांक है। कमजोर रूढ़िवादी रूप में यह बन जाता है
यहां <math>\rho</math> विशेषता (विस्कोसिटी) संकेतक है, और <math>f</math> प्रेरण (फोर्सिंग) संकेतक है।
 
कमजोर रूढ़िवादी रूप में यह बन जाता है
:<math>\left\langle \partial_{t} u , v \right\rangle = \left\langle  \partial_x \left(-\frac{1}{2} u^2 + \rho \partial_{x} u\right) , v \right\rangle + \left\langle f, v \right\rangle \quad \forall v\in \mathcal{V}, \forall t>0</math>
:<math>\left\langle \partial_{t} u , v \right\rangle = \left\langle  \partial_x \left(-\frac{1}{2} u^2 + \rho \partial_{x} u\right) , v \right\rangle + \left\langle f, v \right\rangle \quad \forall v\in \mathcal{V}, \forall t>0</math>
जहां [[आंतरिक उत्पाद स्थान]] संकेतन निम्नलिखित है। [[भागों द्वारा एकीकरण]] और आवधिकता अनुदान का उपयोग करना
जहां [[आंतरिक उत्पाद स्थान|आंतरिक गुणन]] संकेतन निम्नलिखित है। [[भागों द्वारा एकीकरण]] और आवधिकता अनुदान का उपयोग करने पर
:<math>\langle \partial_{t} u , v \rangle = \left\langle  \frac{1}{2} u^2 - \rho \partial_{x} u  ,  \partial_x v\right\rangle+\left\langle f, v \right\rangle \quad \forall v\in \mathcal{V}, \forall t>0.</math>
:<math>\langle \partial_{t} u , v \rangle = \left\langle  \frac{1}{2} u^2 - \rho \partial_{x} u  ,  \partial_x v\right\rangle+\left\langle f, v \right\rangle \quad \forall v\in \mathcal{V}, \forall t>0.</math>
फूरियर-गैलेरकिन विधि लागू करने के लिए, दोनों को चुनें
फूरियर-गैलेरकिन विधि लागू करने के लिए, दोनों को चुना जाता है
:<math>\mathcal{U}^N := \left\{ u : u(x,t)=\sum_{k=-N/2}^{N/2-1} \hat{u}_{k}(t) e^{i k x}\right\}</math>
:<math>\mathcal{U}^N := \left\{ u : u(x,t)=\sum_{k=-N/2}^{N/2-1} \hat{u}_{k}(t) e^{i k x}\right\}</math>
और
और
:<math>\mathcal{V}^N :=\operatorname{span}\left\{ e^{i k x} : k\in -N/2,\dots,N/2-1\right\}</math>
:<math>\mathcal{V}^N :=\operatorname{span}\left\{ e^{i k x} : k\in -N/2,\dots,N/2-1\right\}</math>
कहाँ <math>\hat{u}_k(t):=\frac{1}{2\pi}\langle u(x,t), e^{i k x} \rangle</math>. इससे खोजने में समस्या कम हो जाती है <math>u\in\mathcal{U}^N</math> ऐसा है कि
जहाँ <math>\hat{u}_k(t):=\frac{1}{2\pi}\langle u(x,t), e^{i k x} \rangle</math>. इससे खोजने में समस्या कम हो जाती है।
 
<math>u\in\mathcal{U}^N</math> इस प्रकार है कि
:<math>\langle \partial_{t} u , e^{i k x} \rangle = \left\langle  \frac{1}{2} u^2 - \rho \partial_{x} u  ,  \partial_x e^{i k x}  \right\rangle + \left\langle f, e^{i k x} \right\rangle \quad \forall k\in \left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.</math>
:<math>\langle \partial_{t} u , e^{i k x} \rangle = \left\langle  \frac{1}{2} u^2 - \rho \partial_{x} u  ,  \partial_x e^{i k x}  \right\rangle + \left\langle f, e^{i k x} \right\rangle \quad \forall k\in \left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.</math>
[[ओर्थोगोनालिटी]] संबंध का उपयोग करना <math>\langle e^{i l x}, e^{i k x} \rangle = 2 \pi \delta_{lk}</math> कहाँ <math>\delta_{lk}</math> [[क्रोनकर डेल्टा]] है, हम प्रत्येक के लिए उपरोक्त तीन शब्दों को सरल बनाते हैं <math>k</math> देखने के लिए
[[ओर्थोगोनालिटी|लंबकोणीयता]] संबंध <math>\langle e^{i l x}, e^{i k x} \rangle = 2 \pi \delta_{lk}</math> का उपयोग जहाँ <math>\delta_{lk}</math> [[क्रोनकर डेल्टा]] है, हम प्रत्येक के लिए उपरोक्त तीन शब्दों को सरल बनाते हैं। <math>k</math> देखने के लिए
:<math>
:<math>
\begin{align}
\begin{align}
Line 115: Line 118:
\end{align}
\end{align}
</math>
</math>
प्रत्येक के लिए तीन पद एकत्रित करें <math>k</math> प्राप्त करने के लिए
प्रत्येक के लिए तीन पद एकत्रित करने प <math>k</math> प्राप्त करने के लिए
:<math>
:<math>
2 \pi \partial_t \hat{u}_k
2 \pi \partial_t \hat{u}_k
Line 124: Line 127:
\quad k\in\left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.
\quad k\in\left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.
</math>
</math>
द्वारा विभाजित करना <math>2\pi</math>, हम अंततः पहुँच गए
<math>2\pi</math> द्वारा विभाजित करने पर, हम अंततः निम्नलिखित निष्कर्ष पर पहुँचते हैं
:<math>
:<math>
\partial_t \hat{u}_k
\partial_t \hat{u}_k
Line 133: Line 136:
\quad k\in\left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.
\quad k\in\left\{ -N/2,\dots,N/2-1 \right\}, \forall t>0.
</math>
</math>
फूरियर के साथ प्रारंभिक स्थितियाँ बदल गईं <math>\hat{u}_{k}(0)</math> और जबरदस्ती <math>\hat{f}_{k}(t)</math>, सामान्य अभिविभाज्य समीकरणों की इस युग्मित प्रणाली को समाधान खोजने के लिए समय में एकीकृत किया जा सकता है (उदाहरण के लिए, रनगे कुट्टा तकनीक का उपयोग करके)। अरेखीय शब्द एक संलयन है, और इसे कुशलतापूर्वक मूल्यांकन करने के लिए कई परिवर्तन-आधारित तकनीकें हैं। बॉयड और कैनुटो एट अल के संदर्भ देखें। अधिक जानकारी के लिए।
फुरियर रूपांतरण के साथ प्रारंभिक उपबंध <math>\hat{u}{k}(0)</math> और प्रेरण <math>\hat{f}{k}(t)</math> के साथ, यह संयुक्त सामान्य अभिविभाज्य समीकरणों की प्रणाली को, उदाहरण के लिए, एक [[रुनगे कुट्टा]] तकनीक का उपयोग करके कोई समाधान ढूंढने हेतु, समय के साथ समन्वयित किया जा सकता है। अरेखीय शब्द एक संलयन है, और इसे कुशलतापूर्वक मूल्यांकन करने के लिए कई रूपांतरण-आधारित तकनीकें हैं। अधिक जानकारी के लिए बॉयड और कैनुटो एट अल के संदर्भ देखें।
 
 
 
 
 
 
 
 
 
 
 


== स्पेक्ट्रल तत्व विधि के साथ एक संबंध ==
== स्पेक्ट्रल तत्व विधि के साथ संबंध ==


कोई ऐसा दिखा सकता है अगर <math>g</math> असीम रूप से भिन्न है, तो फास्ट फूरियर ट्रांसफॉर्म का उपयोग करने वाला संख्यात्मक एल्गोरिदम ग्रिड आकार एच में किसी भी बहुपद की तुलना में तेजी से परिवर्तित हो जाएगा। अर्थात्, किसी भी n>0 के लिए, एक है <math>C_n<\infty</math> जिससे त्रुटि कम हो <math>C_nh^n</math> के सभी पर्याप्त छोटे मानों के लिए <math>h</math>. हम कहते हैं कि स्पेक्ट्रल विधि क्रमबद्ध है <math>n</math>, प्रत्येक n>0 के लिए।
यदि <math>g</math> अनंत बार अविभाज्य है, तो फास्ट फ़ौरीयर रूपांतरण का उपयोग करके संख्यात्मक विधिकलन को सिद्ध किया जा सकता है कि यह ग्रिड आकार h के किसी भी बहुपद से शीघ्रता से अभिसरित होगा। अर्थात्, किसी भी n>0 के लिए, एक ऐसी संख्या <math>C_n<\infty</math> है जिसके लिए त्रुटि, <math>h</math> के सभी पर्याप्त छोटे मानों के लिए <math>C_nh^n</math> से कम हो। तब हम कह सकते हैं कि स्पेक्ट्रल विधि प्रत्येक n>0 के लिए <math>n</math> के क्रम में है।


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


== यह भी देखें ==
== यह भी देखें ==
Line 167: Line 181:
{{SpectralTheory}}
{{SpectralTheory}}


{{DEFAULTSORT:Spectral Method}}[[Category: संख्यात्मक विश्लेषण]] [[Category: संख्यात्मक अंतर समीकरण]]
{{DEFAULTSORT:Spectral Method}}
 
 


[[Category: Machine Translated Page]]
[[Category:Collapse templates|Spectral Method]]
[[Category:Created On 24/07/2023]]
[[Category:Created On 24/07/2023|Spectral Method]]
[[Category:Machine Translated Page|Spectral Method]]
[[Category:Navigational boxes| ]]
[[Category:Navigational boxes without horizontal lists|Spectral Method]]
[[Category:Pages with script errors|Spectral Method]]
[[Category:Sidebars with styles needing conversion|Spectral Method]]
[[Category:Template documentation pages|Documentation/doc]]
[[Category:Templates Vigyan Ready|Spectral Method]]
[[Category:Templates generating microformats|Spectral Method]]
[[Category:Templates that are not mobile friendly|Spectral Method]]
[[Category:Templates using TemplateData|Spectral Method]]
[[Category:Wikipedia metatemplates|Spectral Method]]
[[Category:संख्यात्मक अंतर समीकरण|Spectral Method]]
[[Category:संख्यात्मक विश्लेषण|Spectral Method]]

Latest revision as of 14:43, 11 August 2023

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

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

स्पेक्ट्रल विधियों का उपयोग अभिविभाज्य समीकरणों (पीडीई, ओडीई, आइजेनवैल्यू, आदि) और अनुकूलन समस्याओं को हल करने के लिए किया जा सकता है। समय-निर्भर पीडीई के लिए स्पेक्ट्रल विधियों को लागू करते समय, समाधान सामान्यतः समय-निर्भर गुणांक के साथ आधार फलनों के योग के रूप में लिखा जाता है; इसे पीडीई में प्रतिस्थापित करने से गुणांकों में ओडीई की एक प्रणाली प्राप्त होती है जिसे साधारण अभिविभाज्य समीकरणों के लिए किसी भी संख्यात्मक विधियों का उपयोग करके हल किया जा सकता है। ओडीई के लिए आइजेनवैल्यू समस्याओं को इसी तरह आव्यूह आइजेनवैल्यू समस्याओं में परिवर्तित किया जाता है।

1969 में स्टीवन ओर्सज़ैग द्वारा पत्रों की एक लंबी श्रृंखला में स्पेक्ट्रल विधियां विकसित की गईं, जिनमें आवधिक ज्यामिति समस्याओं के लिए फूरियर श्रृंखला विधियां, परिमित और असीमित ज्यामिति समस्याओं के लिए बहुपद स्पेक्ट्रल विधियां, अत्यधिक गैर-रेखीय समस्याओं के लिए छद्मस्पेक्ट्रल विधियां, और स्थिर-अवस्था समस्याओं के तेज़ समाधान के लिए स्पेक्ट्रल पुनरावृत्ति विधियां शामिल हैं, लेकिन इन्हीं तक सीमित नहीं हैं। स्पेक्ट्रल विधि का कार्यान्वयन सामान्यतः या तो सहसंयोजन विधि या गैलेरकिन विधि या ताऊ विधि दृष्टिकोण के साथ पूरा किया जाता है। बहुत छोटी समस्याओं के लिए, स्पेक्ट्रल विधि इस मायने में अद्वितीय है कि समाधानों को प्रतीकात्मक रूप से लिखा जा सकता है, जिससे अभिविभाज्य समीकरणों के लिए श्रृंखला समाधानों का व्यावहारिक विकल्प मिलता है।

परिमित तत्व विधियों की तुलना में स्पेक्ट्रल विधियाँ कम्प्यूटेशनल रूप से कम महंगी और लागू करने में सरल हो सकती हैं; जब सहज समाधानों के साथ सरल क्षेत्र में उच्च सटीकता की मांग की जाती है तो वे सबसे अच्छे विकल्प के रूप में उभरते हैं। यद्यपि, उनकी वैश्विक प्रकृति के कारण, चरण गणना से जुड़े आव्यूह सघन हैं और स्वतंत्रता की कई डिग्री होने पर संगणनीय दक्षता शीघ्रता से प्रभावित होगी। बड़ी समस्याओं और गैर-सुचारू समाधानों के लिए, विरल आव्यूह और असंतुलन और तीव्र घूर्णन के बेहतर प्रारूपण के कारण परिमित तत्व सामान्यतः बेहतर कार्य करते है।

स्पेक्ट्रल विधियों के उदाहरण

एक ठोस, रैखिक उदाहरण

यहां हम आधारभूत बहुभिन्नरूपी कलन(कैल्कुलस) और फूरियर श्रृंखला की समझ का अनुमान लगाते हैं। यदि दो वास्तविक चरों का एक ज्ञात, जटिल-मान फलन है, और g, x और y में आवधिक है (अर्थात्, ) तो हम एक फलन f(x,y) खोजने में रुचि रखते हैं जिससे

जहां बाईं ओर की अभिव्यक्ति क्रमशः x और y में f के दूसरे आंशिक व्युत्पन्न को दर्शाती है। यह पॉइसन समीकरण है, और इसे भौतिक रूप से किसी प्रकार की ऊष्मा चालन समस्या, या अन्य संभावनाओं के बीच संभावित सिद्धांत में एक समस्या के रूप में व्याख्यायित किया जा सकता है।

यदि हम फूरियर श्रृंखला में f और g लिखते हैं:

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

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

 

 

 

 

(*)

जो फूरियर गुणांक j,k. के लिए एक स्पष्ट सूत्र है

आवधिक सीमा स्थितियों के साथ, पॉइसन समीकरण का कोई समाधान केवल तभी होता है जब b0,0 = 0 होता है। इसलिए, हम स्वतंत्र रूप से a0,0 चुन सकते हैं जो विश्लेषण के माध्य के बराबर होगा। यह एकीकरण स्थिरांक को चुनने के अनुरूप है।

इसे एक विधिकलन में परिवर्तित करने के लिए, केवल परिमित आवृत्तियों को हल किया जाता है। यह एक त्रुटि प्रस्तुत करता है जिसे के आनुपातिक दिखाया जा सकता है , जहाँ और उपचारित उच्चतम आवृत्ति है।

विधिकलन

  1. g के फूरियर रूपांतरण (bj,k) की गणना करें।
  2. (*) सूत्र के माध्यम से f के फूरियर रूपांतरण (aj,k) की गणना करें। .
  3. (aj,k) का व्युत्क्रम फूरियर रूपांतरण लेकर f की गणना करें

चूँकि हम केवल आवृत्तियों की एक सीमित क्षेत्र (जैसे आकार n,) में रुचि रखते हैं, यह एक फास्ट फूरियर ट्रांसफॉर्म विधिकलन का उपयोग करके किया जा सकता है। इसलिए, विश्व स्तर पर विधिकलन time O(n log n). समय में चलता है

अरेखीय उदाहरण

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

पुर्तस्य क्षेत्र पर के लिए दिए गए के साथ ऐसे का पता लगाएं जो निम्नलिखित समीकरण को पूरा करते हैं:

यहां विशेषता (विस्कोसिटी) संकेतक है, और प्रेरण (फोर्सिंग) संकेतक है।

कमजोर रूढ़िवादी रूप में यह बन जाता है

जहां आंतरिक गुणन संकेतन निम्नलिखित है। भागों द्वारा एकीकरण और आवधिकता अनुदान का उपयोग करने पर

फूरियर-गैलेरकिन विधि लागू करने के लिए, दोनों को चुना जाता है

और

जहाँ . इससे खोजने में समस्या कम हो जाती है।

इस प्रकार है कि

लंबकोणीयता संबंध का उपयोग जहाँ क्रोनकर डेल्टा है, हम प्रत्येक के लिए उपरोक्त तीन शब्दों को सरल बनाते हैं। देखने के लिए

प्रत्येक के लिए तीन पद एकत्रित करने प प्राप्त करने के लिए

द्वारा विभाजित करने पर, हम अंततः निम्नलिखित निष्कर्ष पर पहुँचते हैं

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







स्पेक्ट्रल तत्व विधि के साथ संबंध

यदि अनंत बार अविभाज्य है, तो फास्ट फ़ौरीयर रूपांतरण का उपयोग करके संख्यात्मक विधिकलन को सिद्ध किया जा सकता है कि यह ग्रिड आकार h के किसी भी बहुपद से शीघ्रता से अभिसरित होगा। अर्थात्, किसी भी n>0 के लिए, एक ऐसी संख्या है जिसके लिए त्रुटि, के सभी पर्याप्त छोटे मानों के लिए से कम हो। तब हम कह सकते हैं कि स्पेक्ट्रल विधि प्रत्येक n>0 के लिए के क्रम में है।

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

यह भी देखें

संदर्भ

  1. pp 235, Spectral Methods: evolution to complex geometries and applications to fluid dynamics, By Canuto, Hussaini, Quarteroni and Zang, Springer, 2007.
  • Bengt Fornberg (1996) A Practical Guide to Pseudospectral Methods. Cambridge University Press, Cambridge, UK
  • Chebyshev and Fourier Spectral Methods by John P. Boyd.
  • Canuto C., Hussaini M. Y., Quarteroni A., and Zang T.A. (2006) Spectral Methods. Fundamentals in Single Domains. Springer-Verlag, Berlin Heidelberg
  • Javier de Frutos, Julia Novo (2000): A Spectral Element Method for the Navier–Stokes Equations with Improved Accuracy
  • Polynomial Approximation of Differential Equations, by Daniele Funaro, Lecture Notes in Physics, Volume 8, Springer-Verlag, Heidelberg 1992
  • D. Gottlieb and S. Orzag (1977) "Numerical Analysis of Spectral Methods : Theory and Applications", SIAM, Philadelphia, PA
  • J. Hesthaven, S. Gottlieb and D. Gottlieb (2007) "Spectral methods for time-dependent problems", Cambridge UP, Cambridge, UK
  • Steven A. Orszag (1969) Numerical Methods for the Simulation of Turbulence, Phys. Fluids Supp. II, 12, 250–257
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Section 20.7. Spectral Methods". Numerical Recipes: The Art of Scientific Computing (3rd ed.). New York: Cambridge University Press. ISBN 978-0-521-88068-8.
  • Jie Shen, Tao Tang and Li-Lian Wang (2011) "Spectral Methods: Algorithms, Analysis and Applications" (Springer Series in Computational Mathematics, V. 41, Springer), ISBN 354071040X
  • Lloyd N. Trefethen (2000) Spectral Methods in MATLAB. SIAM, Philadelphia, PA