गतिशील मोड अपघटन: Difference between revisions
No edit summary |
No edit summary |
||
Line 1: | Line 1: | ||
{{Short description|Dimensionality reduction algorithm}}गतिशील मोड अपघटन (डीएमडी) आयामिकता घटाने का एल्गोरिदम है जिसका विकास 2008 में पीटर जे. श्मिड और जोर्न सेस्टरहेन ने किया था। डेटा की समय श्रृंखला के आधार पर, डीएमडी सेट के मोड की गणना करता है, जिनमें से प्रत्येक निश्चित दोलन आवृत्ति और क्षय/विकास दर से जुड़ा संबंधित होती है। विशेष रूप से रैखिक प्रणालियों के लिए, ये मोड और आवृत्तियाँ | {{Short description|Dimensionality reduction algorithm}}गतिशील मोड अपघटन (डीएमडी) आयामिकता घटाने का एल्गोरिदम है जिसका विकास 2008 में पीटर जे. श्मिड और जोर्न सेस्टरहेन ने किया था। डेटा की समय श्रृंखला के आधार पर, डीएमडी सेट के मोड की गणना करता है, जिनमें से प्रत्येक निश्चित दोलन आवृत्ति और क्षय/विकास दर से जुड़ा संबंधित होती है। विशेष रूप से रैखिक प्रणालियों के लिए, ये मोड और आवृत्तियाँ प्रणाली के [[सामान्य मोड]] के समानांतर होते हैं, किन्तु अधिक सासामान्य रूप से, वे [[रचना संचालक]] (जिसे कोपमैन ऑपरेटर भी कहा जाता है) के मोड और आइगेनवैल्यू के अनुमान होते हैं। प्रत्येक मोड से जुड़े आंतरिक अस्थायी व्यवहारों के कारण, डीएमडी प्रमुख घटक विश्लेषण जैसे आयामी कमी के तरीकों से भिन्न होता है, जो प्राधिकृती प्रतिस्थापन विधियों (प्रिंसिपल कॉम्पोनेंट एनालिसिस) के अतिरिक्त होते हैं। क्योंकि इसके मोड ऑर्थोगोनल नहीं हैं, डीएमडी-आधारित अभ्यावेदन पीसीए द्वारा उत्पन्न किए गए अभ्यावेदन की तुलना में कम सरल हो सकते हैं। चूंकि, वे अधिक शारीरिक रूप से सार्थक भी हो सकते हैं क्योंकि प्रत्येक मोड समय में नम (या संचालित) साइनसॉइडल व्यवहार से जुड़ा होता है। | ||
==अवलोकन== | ==अवलोकन== | ||
डायनेमिक मोड अपघटन को पहली बार श्मिड द्वारा फ्लो डेटा से गतिजंक्य विशेषताएं निकालने के लिए संख्यात्मक प्रक्रिया के रूप में | डायनेमिक मोड अपघटन को पहली बार श्मिड द्वारा फ्लो डेटा से गतिजंक्य विशेषताएं निकालने के लिए संख्यात्मक प्रक्रिया के रूप में प्रस्तुत किया गया था।<ref name=schmid2010>P.J. Schmid. "Dynamic mode decomposition of numerical and experimental data." Journal of Fluid Mechanics 656.1 (2010): 5–28.</ref> | ||
डेटा स्नैपशॉट अनुक्रम के रूप में प्राप्त होता है, जहां वी1 से एन तक के लिए: | डेटा स्नैपशॉट अनुक्रम के रूप में प्राप्त होता है, जहां वी1 से एन तक के लिए: | ||
: <math> V_1^N = \{v_1, v_2, \dots, v_N\}, </math> | : <math> V_1^N = \{v_1, v_2, \dots, v_N\}, </math> | ||
यहाँ <math> v_i\in \mathbb{R}^M</math> व्युत्पन्न फ्लो फ़ील्ड की <math>i</math>-प्रवाह स्नैपशॉट है, और <math>V_1^N\in\mathbb{R}^{M\times N}</math> डेटा मैट्रिक्स है जिसके कॉलम व्यक्तिगत स्नैपशॉट होते हैं। इन स्नैपशॉट कोरैखिक गणना के माध्यम से संबंधित होते हैं जो रैखिक गतिजंकीय प्रणाली को परिभाषित करता है: | यहाँ <math> v_i\in \mathbb{R}^M</math> व्युत्पन्न फ्लो फ़ील्ड की <math>i</math>-प्रवाह स्नैपशॉट है, और <math>V_1^N\in\mathbb{R}^{M\times N}</math> डेटा मैट्रिक्स है जिसके कॉलम व्यक्तिगत स्नैपशॉट होते हैं। इन स्नैपशॉट कोरैखिक गणना के माध्यम से संबंधित होते हैं जो रैखिक गतिजंकीय प्रणाली को परिभाषित करता है: | ||
: <math> v_{i+1} = A v_i, </math>जो समय के सैंपलिंग अवधि के | : <math> v_{i+1} = A v_i, </math>जो समय के सैंपलिंग अवधि के समय लगभग समान रहता है। मैट्रिक्स रूप में लिखें, इससे यह सिद्ध होता है कि: | ||
:<math> V_{2}^N = A V_{1}^{N-1} + re_{N-1}^T,</math> | :<math> V_{2}^N = A V_{1}^{N-1} + re_{N-1}^T,</math> | ||
यहाँ <math>r</math> अवशेषों का वेक्टर है जो ऐसे व्यवहारों का कारण बनता है जिनका पूरी तरह से <math>A</math>वर्णन नहीं किया जा सकता, <math>e_{N-1}=\{0,0,\ldots,1\}\in\mathbb{R}^{N-1}</math>है, <math>V_1^{N-1}=\{v_1, v_2, \dots, v_{N-1}\}</math>, और <math>V_2^{N}=\{v_2, v_3, \dots, v_{N}\}</math> चाहे किसी भी दृष्टिकोण से, डीएमडी का आउटपुट <math>A</math>, के इगेनवैल्यूज और इगेनवेक्टर होते हैं, जिन्हें डीएमडी इगेनवैल्यूज और डीएमडी मोड के रूप में संदर्भित किया जाता है। | यहाँ <math>r</math> अवशेषों का वेक्टर है जो ऐसे व्यवहारों का कारण बनता है जिनका पूरी तरह से <math>A</math>वर्णन नहीं किया जा सकता, <math>e_{N-1}=\{0,0,\ldots,1\}\in\mathbb{R}^{N-1}</math>है, <math>V_1^{N-1}=\{v_1, v_2, \dots, v_{N-1}\}</math>, और <math>V_2^{N}=\{v_2, v_3, \dots, v_{N}\}</math> चाहे किसी भी दृष्टिकोण से, डीएमडी का आउटपुट <math>A</math>, के इगेनवैल्यूज और इगेनवेक्टर होते हैं, जिन्हें डीएमडी इगेनवैल्यूज और डीएमडी मोड के रूप में संदर्भित किया जाता है। | ||
Line 12: | Line 12: | ||
== एल्गोरिथम == | == एल्गोरिथम == | ||
इन इगेनवैल्यूज और मोड्स को प्राप्त करने की दो विधियाँ हैं। पहला विधि अर्नोल्डी जैसा है, जो सिद्धांतिक विश्लेषण के लिए उपयुक्त है क्योंकि इसका क्रायलोव विधियों से संबंध होता है। दूसरा सिंगुलर वैल्यू विचलन (एसवीडी) पर आधारित विधि है जो डेटा में | इन इगेनवैल्यूज और मोड्स को प्राप्त करने की दो विधियाँ हैं। पहला विधि अर्नोल्डी जैसा है, जो सिद्धांतिक विश्लेषण के लिए उपयुक्त है क्योंकि इसका क्रायलोव विधियों से संबंध होता है। दूसरा सिंगुलर वैल्यू विचलन (एसवीडी) पर आधारित विधि है जो डेटा में ध्वनि से अधिक ट्रुटियों और संख्यात्मक त्रुटियों के प्रति शक्तिशाली है। | ||
=== अर्नोल्डी दृष्टिकोण === | === अर्नोल्डी दृष्टिकोण === | ||
तरल पदार्थ अनुप्रयोगों में, स्नैपशॉट का आकार, <math>M</math>, स्नैपशॉट की संख्या <math>N</math> से बहुत अधिक बड़ा माना जाता है, इसलिए समान रूप से कई वैध विकल्प <math>A</math> | तरल पदार्थ अनुप्रयोगों में, स्नैपशॉट का आकार, <math>M</math>, स्नैपशॉट की संख्या <math>N</math> से बहुत अधिक बड़ा माना जाता है, इसलिए समान रूप से कई वैध विकल्प <math>A</math> उपस्थित हैं। मूल डीएमडी एल्गोरिदम <math>A</math> चुनता है जिससे प्रत्येक स्नैपशॉट में <math>V_2^N</math> में <math>V_1^{N-1}</math>के स्नैपशॉटों के रूप में लिनियर संयोजन के रूप में व्यक्त किया जा सकता है। क्योंकि अधिकांश स्नैपशॉट दोनों डेटा सेट में दिखाई देते हैं, इस प्रतिनिधित्व को स्नैपशॉट <math>v_N</math>, को छोड़कर सभी स्नैपशॉटों के लिए त्रुटिरहित बनाता है, जो इस तरह से लिखा जाता है: | ||
:<math> v_N = a_1 v_1 + a_2 v_2 + \dots + a_{N-1}v_{N-1} + r = V_1^{N-1}a + r,</math> | :<math> v_N = a_1 v_1 + a_2 v_2 + \dots + a_{N-1}v_{N-1} + r = V_1^{N-1}a + r,</math> | ||
यहाँ <math>a={a_1, a_2, \dots, a_{N-1}}</math> गुणांकों का सेट है जिसे डीएमडी को पहचानना चाहिए और <math>r</math> अवशेष है। | यहाँ <math>a={a_1, a_2, \dots, a_{N-1}}</math> गुणांकों का सेट है जिसे डीएमडी को पहचानना चाहिए और <math>r</math> अवशेष है। | ||
Line 32: | Line 32: | ||
सदिश <math>a</math> न्यूनतम वर्ग समस्या को हल करके गणना की जा सकती है, जो समग्र अवशिष्ट को न्यूनतम करती है। विशेष रूप से यदि हम <math>V_1^{N-1} = QR</math> का QR विचलन लें, तो <math>a = R^{-1}Q^Tv_N</math> होता है। | सदिश <math>a</math> न्यूनतम वर्ग समस्या को हल करके गणना की जा सकती है, जो समग्र अवशिष्ट को न्यूनतम करती है। विशेष रूप से यदि हम <math>V_1^{N-1} = QR</math> का QR विचलन लें, तो <math>a = R^{-1}Q^Tv_N</math> होता है। | ||
इस रूप में, डीएमडी प्रकार का अर्नोल्डी पुनरावृत्ति है, और इसलिए <math>S</math> के इगेनवैल्यूज के <math>A</math> अनुमान होते हैं। इसके | इस रूप में, डीएमडी प्रकार का अर्नोल्डी पुनरावृत्ति है, और इसलिए <math>S</math> के इगेनवैल्यूज के <math>A</math> अनुमान होते हैं। इसके अतिरिक्त, यदि <math>y</math>, <math>S</math> का इगेनवेक्टर है, तब <math>V_1^{N-1}y</math>, <math>A</math> का अउपस्थित इगेनवेक्टर होता है। यहां इगेनवैल्यूज विचलन <math>A</math> के बजाए <math>S</math> पर किया जाता है क्योंकि <math>S</math>, <math>A</math> से अधिक छोटा होता है, इसलिए डीएमडी की संगणना निवेश स्नैपशॉट के आकार के अतिरिक्त स्नैपशॉट की संख्या से निर्धारित होती है। | ||
=== एसवीडी-आधारित दृष्टिकोण === | === एसवीडी-आधारित दृष्टिकोण === | ||
साथी मैट्रिक्स की गणना करने के | साथी मैट्रिक्स की गणना करने के अतिरिक्त <math> S </math>, एसवीडी-आधारित दृष्टिकोण मैट्रिक्स उत्पन्न करता है <math>\tilde S</math> वह संबंधित है <math>A</math> समानता परिवर्तन के माध्यम से. ऐसा करने के लिए, मान लें कि हमारे पास एसवीडी है <math>V_1^{N-1} = U\Sigma W^T</math>. तब | ||
:<math> V_{2}^N = A V_1^{N-1} + re_{N-1}^T = AU\Sigma W^T + re_{N-1}^T.</math> | :<math> V_{2}^N = A V_1^{N-1} + re_{N-1}^T = AU\Sigma W^T + re_{N-1}^T.</math> | ||
अर्नोल्डी-आधारित दृष्टिकोण द्वारा बनाई गई धारणा के बराबर, हम <math>A</math> चुनते हैं जैसे कि <math>V_2^N</math> स्नैपशॉट में <math>U</math> में स्तंभों के रैखिक सुपरपोजिशन के रूप में लिखा जा सकता है , जो इस आवश्यकता के बराबर है कि उन्हें प्रमुख घटक विश्लेषण के सुपरपोजिशन के रूप में लिखा जा सकता है। इस प्रतिबंध के साथ, अवशिष्ट को कम करने के लिए आवश्यक है कि यह POD आधार पर ऑर्थोगोनल हो (अर्थात्, <math>U^Tr = 0</math>)। फिर समीकरण को <math>U^T</math> से दोनों तरफ़ से गुणा करने से इसे मिलाया जा सकता है: <math> U^TV_2^N = U^T A U\Sigma W^T </math>, जिसे स्पष्ट किया जा सकता है: | अर्नोल्डी-आधारित दृष्टिकोण द्वारा बनाई गई धारणा के बराबर, हम <math>A</math> चुनते हैं जैसे कि <math>V_2^N</math> स्नैपशॉट में <math>U</math> में स्तंभों के रैखिक सुपरपोजिशन के रूप में लिखा जा सकता है , जो इस आवश्यकता के बराबर है कि उन्हें प्रमुख घटक विश्लेषण के सुपरपोजिशन के रूप में लिखा जा सकता है। इस प्रतिबंध के साथ, अवशिष्ट को कम करने के लिए आवश्यक है कि यह POD आधार पर ऑर्थोगोनल हो (अर्थात्, <math>U^Tr = 0</math>)। फिर समीकरण को <math>U^T</math> से दोनों तरफ़ से गुणा करने से इसे मिलाया जा सकता है: <math> U^TV_2^N = U^T A U\Sigma W^T </math>, जिसे स्पष्ट किया जा सकता है: | ||
:<math> U^T A U = U^TV_2^N W \Sigma^{-1} \equiv \tilde S. </math> | :<math> U^T A U = U^TV_2^N W \Sigma^{-1} \equiv \tilde S. </math> | ||
क्योंकि <math>A</math> और <math>\tilde S</math> समानता संवेदी रूप से जुड़े हुए हैं, इसलिए <math>S</math> के इगेनवैल्यूज <math>A</math> के इगेनवैल्यूज होते हैं, और | क्योंकि <math>A</math> और <math>\tilde S</math> समानता संवेदी रूप से जुड़े हुए हैं, इसलिए <math>S</math> के इगेनवैल्यूज <math>A</math> के इगेनवैल्यूज होते हैं, और यदि <math>y</math>, <math>\tilde S</math> का इगेनवेक्टर है , तब <math>Uy</math>, <math>A</math> का इगेनवेक्टर होता है। . | ||
संक्षेप में, एसवीडी-आधारित | संक्षेप में, एसवीडी-आधारित विधि का अनुसरण निम्नलिखित तरीके से किया जाता है: | ||
# डेटा की समय श्रृंखला को विभाजित करें <math>V_1^N</math> दो आव्यूहों में <math>V_1^{N-1}</math> और <math>V_2^N</math>. | # डेटा की समय श्रृंखला को विभाजित करें <math>V_1^N</math> दो आव्यूहों में <math>V_1^{N-1}</math> और <math>V_2^N</math>. | ||
# के एसवीडी की गणना करें <math>V_1^{N-1} = U\Sigma W^T</math>. | # के एसवीडी की गणना करें <math>V_1^{N-1} = U\Sigma W^T</math>. | ||
Line 48: | Line 48: | ||
# <math>i</math>th>-th डीएमडी इगेनवैल्यूज है <math>\lambda_i</math> और <math>i</math>-वें डीएमडी मोड है <math>Uy_i</math>. | # <math>i</math>th>-th डीएमडी इगेनवैल्यूज है <math>\lambda_i</math> और <math>i</math>-वें डीएमडी मोड है <math>Uy_i</math>. | ||
अर्नोल्डी-जैसे दृष्टिकोण पर एसवीडी-आधारित दृष्टिकोण का लाभ यह है कि डेटा में | अर्नोल्डी-जैसे दृष्टिकोण पर एसवीडी-आधारित दृष्टिकोण का लाभ यह है कि डेटा में ध्वनि और संख्यात्मक ट्रंकेशन मुद्दों की भरपाई एसवीडी को काटकर की जा सकती है। <math>V_1^{N-1}</math>. जैसा कि उल्लेख किया गया है <ref name=schmid2010 />इस ट्रंकेशन चरण के बिना प्रयोगात्मक डेटा सेट पर पहले जोड़े मोड और आइगेनवैल्यू से अधिक की स्पष्ट गणना करना कठिनाई हो सकता है। | ||
== सैद्धांतिक और एल्गोरिथम प्रगति == | == सैद्धांतिक और एल्गोरिथम प्रगति == | ||
इसकी प्रारंभिकता से 2010 में, डीएमडी को समझने और सुधारने पर | इसकी प्रारंभिकता से 2010 में, डीएमडी को समझने और सुधारने पर अधिक काम किया गया है।<ref>C.W. Rowley, I Mezic, S. Bagheri, P. Schlatter, and D.S. Henningson, "Spectral analysis of nonlinear flows." Journal of Fluid Mechanics 641 (2009): 85-113</ref> Rowley एवं उनके सहयोगियों ने डीएमडी का पहला विश्लेषण किया जिसमें डीएमडी और Koopman ऑपरेटर के बीच संबंध स्थापित किया गया और यह समझने में सहायता की जब डीएमडी को गैर-लीनियर प्रणालियों पर प्रयुक्त किया गया। इसके बाद से, कई संशोधन विकसित किए गए हैं जो इस संबंध को और मज़बूत करते हैं या इस दृष्टिकोन को अधिक उत्कृष्ट और व्यापक बनाने में सहायता करते हैं। | ||
*अनुकूलित डीएमडी: अनुकूलित डीएमडी मूल डीएमडी एल्गोरिदम का संशोधन है जिसे उस दृष्टिकोण की दो सीमाओं का समाधान करने के लिए डिज़ाइन किया गया है: (i) डीएमडी मोड चयन की कठिनाई, और (ii) डीएमडी की संवेदनशीलता | *अनुकूलित डीएमडी: अनुकूलित डीएमडी मूल डीएमडी एल्गोरिदम का संशोधन है जिसे उस दृष्टिकोण की दो सीमाओं का समाधान करने के लिए डिज़ाइन किया गया है: (i) डीएमडी मोड चयन की कठिनाई, और (ii) डीएमडी की संवेदनशीलता ध्वनि या समय शृंखला के अंतिम स्नैपशॉट में अन्य त्रुटियों के प्रति।<ref>K.K. Chen, J.H. Tu, and C.W. Rowley, "Variants of dynamic mode decomposition: boundary condition, Koopman, and Fourier analyses." Journal of Nonlinear Science 22 (2012): 887-915.</ref> ऑप्टिमाइज़्ड डीएमडी डीएमडी प्रक्रिया को अनुकूलन समस्या के रूप में पुनर्स्थापित करता है जहाँ पहचानी गई लीनियर ऑपरेटर की निश्चित श्रेणी होती है। इसके अतिरिक्त, जिस तरह डीएमडी ने पिछले सभी स्नैपशॉट को पूरी तरह से पुनर्जन्मित किया है, ऑप्टिमाइज़्ड डीएमडी द्वारा रीकंस्ट्रक्शन त्रुटियों को डेटा सेट में वितरित किया जा सकता है, जिससे यह विधि व्यावसायिक रूप से अधिक शक्तिशाली दिखता है। | ||
*इष्टतम मोड अपघटन: इष्टतम मोड अपघटन (ओएमडी) डीएमडी प्रक्रिया को अनुकूलन समस्या के रूप में पुनर्गठित करता है और उपयोगकर्ता को सीधे पहचाने गए | *इष्टतम मोड अपघटन: इष्टतम मोड अपघटन (ओएमडी) डीएमडी प्रक्रिया को अनुकूलन समस्या के रूप में पुनर्गठित करता है और उपयोगकर्ता को सीधे पहचाने गए प्रणाली की रैंक प्रयुक्त करने की अनुमति देता है।<ref>A. Wynn, D. S. Pearson, B. Ganapathisubramani and P. J. Goulart, "Optimal mode decomposition for unsteady flows." Journal of Fluid Mechanics 733 (2013): 473-503</ref> यदि यह श्रेणी योग्य ढंग से चुनी जाए, तो OMD सांदर्भिक और प्रायोगिक डेटा सेटों पर छोटे शेष त्रुटि और अधिक स्पष्ट इजेनवैल्यूज़ के साथ लीनियर मॉडल उत्पन्न कर सकता है। | ||
* | *स्पष्ट डीएमडी: स्पष्ट डीएमडी एल्गोरिदम मूल डीएमडी एल्गोरिदम को दो तरीकों से सामान्यीकृत करता है। सबसे पहले, मूल डीएमडी एल्गोरिदम में डेटा स्नैपशॉट की समय श्रृंखला होनी चाहिए, किन्तु स्पष्ट डीएमडी स्नैपशॉट जोड़े के डेटा सेट को स्वीकार करता है।<ref name=tu2014>{{Cite journal|url = http://www.aimsciences.org/journals/displayArticlesnew.jsp?paperID=10631|title = On Dynamic Mode Decomposition: Theory and Applications|last = Tu, Rowley, Luchtenburg, Brunton, and Kutz|date = December 2014|journal = American Institute of Mathematical Sciences| volume=1 | issue=2 | pages=391–421 |doi = 10.3934/jcd.2014.1.391|arxiv = 1312.0041| s2cid=46419148 }}</ref> जोड़ी में स्नैपशॉट को निश्चित <math>\Delta t</math> द्वारा अलग किया जाना चाहिए , किन्तु ही समय श्रृंखला से निकालने की आवश्यकता नहीं है। विशेष रूप से, स्पष्ट डीएमडी कई प्रयोगों के डेटा को ही डेटा सेट में एकत्रित करने की अनुमति दे सकता है। दूसरा, मूल डीएमडी एल्गोरिदम पीओडी मोड के सेट पर प्रोजेक्ट करके डेटा को प्रभावी ढंग से प्री-प्रोसेस करता है। स्पष्ट डीएमडी एल्गोरिदम इस प्री-प्रोसेसिंग चरण को हटा देता है, और डीएमडी मोड का उत्पादन कर सकता है जिसे पीओडी मोड के सुपरपोजिशन के रूप में नहीं लिखा जा सकता है। | ||
*स्पार्सिटी प्रमोशन डीएमडी: स्पार्सिटी प्रोमोटिंग डीएमडी, डीएमडी मोड और इजेनवैल्यू चयन के लिए पोस्ट प्रोसेसिंग प्रक्रिया है।<ref>M.R. Jovanovic, P.J. Schmid, and J.W. Nichols, "Sparsity-promoting dynamic mode decomposition." Physics of Fluids 26 (2014)</ref> स्पार्सिटी प्रोमोटिंग डीएमडी <math>\ell_1</math> पेनल्टी का उपयोग करता है जो महत्वपूर्ण डीएमडी मोड्स का छोटा सेट पहचानने के लिए उपयुक्त होता है, और यह डीएमडी मोड चयन समस्या के वैकल्पिक दृष्टिकोन है जिसे [[संवर्धित लैग्रेंजियन विधि]] का उपयोग करके कुशलतापूर्वक हल किया जा सकता है। | *स्पार्सिटी प्रमोशन डीएमडी: स्पार्सिटी प्रोमोटिंग डीएमडी, डीएमडी मोड और इजेनवैल्यू चयन के लिए पोस्ट प्रोसेसिंग प्रक्रिया है।<ref>M.R. Jovanovic, P.J. Schmid, and J.W. Nichols, "Sparsity-promoting dynamic mode decomposition." Physics of Fluids 26 (2014)</ref> स्पार्सिटी प्रोमोटिंग डीएमडी <math>\ell_1</math> पेनल्टी का उपयोग करता है जो महत्वपूर्ण डीएमडी मोड्स का छोटा सेट पहचानने के लिए उपयुक्त होता है, और यह डीएमडी मोड चयन समस्या के वैकल्पिक दृष्टिकोन है जिसे [[संवर्धित लैग्रेंजियन विधि]] का उपयोग करके कुशलतापूर्वक हल किया जा सकता है। | ||
*मल्टी-रिज़ॉल्यूशन डीएमडी: मल्टी-रिज़ॉल्यूशन डीएमडी (एमआरडीएमडी) | *मल्टी-रिज़ॉल्यूशन डीएमडी: मल्टी-रिज़ॉल्यूशन डीएमडी (एमआरडीएमडी) स्पष्ट डीएमडी के साथ मल्टी-रिज़ॉल्यूशन विश्लेषण में उपयोग की जाने वाली विधिों का संयोजन है, जिसे कई टाइमस्केल वाले डेटा सेट से डीएमडी मोड और आइजेनवैल्यू को शक्तिशाली करने के लिए डिज़ाइन किया गया है।<ref>J.N. Kutz, X. Fu, and S.L. Brunton, "Multi-resolution dynamic mode decomposition." arXiv preprint [https://arxiv.org/abs/1506.00564 arXiv:1506.00564] (2015).</ref> एमआरडीएमडी दृष्टिकोण वैश्विक सतह तापमान डेटा पर प्रयुक्त किया गया था, और एल नीनो वर्षों के समय दिखाई देने वाले डीएमडी मोड की पहचान करता है। | ||
*विस्तारित डीएमडी: विस्तारित डीएमडी | *विस्तारित डीएमडी: विस्तारित डीएमडी स्पष्ट डीएमडी का संशोधन है जो डीएमडी और कूपमैन ऑपरेटर के बीच संबंध को शक्तिशाली करता है।<ref>M.O. Williams , I.G. Kevrekidis, C.W. Rowley, "A Data–Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition." Journal of Nonlinear Science 25 (2015): 1307-1346.</ref> जैसा कि नाम से पता चलता है, विस्तारित डीएमडी डीएमडी का विस्तार है जो कोपमैन ऑपरेटर के अधिक स्पष्ट अनुमान उत्पन्न करने के लिए अवलोकन योग्य कार्यों के समृद्ध सेट का उपयोग करता है। इस विस्तारित सेट को प्राथमिकता से चुना जा सकता है या डेटा से सीखा जा सकता है।<ref>{{Cite journal|last1=Li|first1=Qianxiao|first2=Felix|last2=Dietrich|first3=Erik M.|last3=Bollt|first4=Ioannis G|last4=Kevrekidis|date=2018|title=Extended dynamic mode decomposition with dictionary learning: A data-driven adaptive spectral decomposition of the Koopman operator|journal=Chaos: An Interdisciplinary Journal of Nonlinear Science|volume=27|number=10|pages=103111|doi=10.1063/1.4993854|pmid=29092410 |arxiv=1707.00225|s2cid=41957686 }}</ref><ref>{{Cite journal|last1=Gulina|first1=Marvyn|first2=Alexandre|last2=Mauroy|date=2021|title=जलाशय कंप्यूटर के साथ कूपमैन ऑपरेटर का अनुमान लगाने की दो विधियाँ|journal=Chaos: An Interdisciplinary Journal of Nonlinear Science|volume=31|number=2|pages=023116|doi=10.1063/5.0026380|pmid=33653036 |arxiv=2008.10263|bibcode=2021Chaos..31b3116G |s2cid=221266743 }}</ref> इसने यह भी प्रदर्शित किया कि डीएमडी और संबंधित विधियाँ अधिक सामान्यतः उपयोग किए जाने वाले आइगेनवैल्यू और मोड के अतिरिक्त कोपमैन आइजेनफंक्शन के सन्निकटन का उत्पादन करती हैं। | ||
*नियंत्रण के साथ डीएमडी: नियंत्रण के साथ डायनामिक मोड अपघटन (डीएमडीसी) <ref>J.L. Proctor, S.L. Brunton, and J.N. Kutz, "Dynamic mode decomposition with control." arXiv preprint arXiv:1409.6358 (2014).</ref> इनपुट आउटपुट | *नियंत्रण के साथ डीएमडी: नियंत्रण के साथ डायनामिक मोड अपघटन (डीएमडीसी) <ref>J.L. Proctor, S.L. Brunton, and J.N. Kutz, "Dynamic mode decomposition with control." arXiv preprint arXiv:1409.6358 (2014).</ref> इनपुट आउटपुट प्रणाली से प्राप्त डेटा के लिए डिज़ाइन की गई डीएमडी प्रक्रिया का संशोधन है। डीएमडीसी की अनूठी विशेषता ओपन लूप डायनेमिक्स से प्रणाली एक्चुएशन के प्रभावों को स्पष्ट करने की क्षमता है, जो तब उपयोगी होती है जब एक्चुएशन की उपस्थिति में डेटा प्राप्त किया जाता है। | ||
*कुल न्यूनतम वर्ग डीएमडी: कुल न्यूनतम वर्ग डीएमडी हाल ही की बदलाव है जो एक्सैक्ट डीएमडी की संशोधन है, जिसका उद्देश्य डेटा में मापन | *कुल न्यूनतम वर्ग डीएमडी: कुल न्यूनतम वर्ग डीएमडी हाल ही की बदलाव है जो एक्सैक्ट डीएमडी की संशोधन है, जिसका उद्देश्य डेटा में मापन ध्वनि के प्रति स्थायित्व से संबंधित मुद्दों का सामना करना है।<ref>M.S. Hemati, C.W. Rowley, E.A. Deem, and L.N. Cattafesta, "De-Biasing the Dynamic Mode Decomposition for Applied Koopman Spectral Analysis of Noisy Datasets." arXiv preprint arXiv:1502.03854 (2015).</ref> में, लेखक एक्सैक्ट डीएमडी को प्रेक्षण समस्या के रूप में व्याख्या करते हैं जो साधारण लीस्ट स्क्वेयर्स (ओएलएस) का उपयोग करके हल किया जाता है, जो इसमें अनुमानित करता है कि प्रेक्षक ध्वनियुक्त मुक्त हैं। यह धारणा डीएमडी आइगेनवैल्यू में पूर्वाग्रह उत्पन्न करती है जब इसे प्रयोगात्मक डेटा सेट पर प्रयुक्त किया जाता है जहां सभी अवलोकन ध्वनि होते हैं। [[कुल न्यूनतम वर्ग]] डीएमडी ओएलएस समस्या को कुल न्यूनतम वर्ग से प्रतिस्थापित करता है, जो इस पूर्वाग्रह को समाप्त करता है। | ||
*गतिशील वितरण अपघटन: डीडीडी निरंतर समय में आगे की समस्या, | *गतिशील वितरण अपघटन: डीडीडी निरंतर समय में आगे की समस्या, अर्थात [[ स्थानांतरण ऑपरेटर |स्थानांतरण ऑपरेटर]] पर ध्यान केंद्रित करता है। चूंकि विकसित विधि का उपयोग निरंतर समय में डीएमडी समस्याओं को ठीक करने के लिए भी किया जा सकता है।<ref>{{Cite journal|last1=Taylor-King|first1=Jake P.|last2=Riseth|first2=Asbjørn N.|last3=Macnair|first3=Will|last4=Claassen|first4=Manfred|date=2020-01-10|title=एकल-सेल स्नैपशॉट समय श्रृंखला के लिए गतिशील वितरण अपघटन आईपीएससी रीप्रोग्रामिंग के दौरान उप-जनसंख्या और प्रक्षेपवक्र की पहचान करता है|journal=PLOS Computational Biology|language=en|volume=16|issue=1|pages=e1007491|doi=10.1371/journal.pcbi.1007491|issn=1553-7358|pmc=6953770|pmid=31923173|bibcode=2020PLSCB..16E7491T }}</ref> | ||
यहां सूचीबद्ध एल्गोरिदम के | यहां सूचीबद्ध एल्गोरिदम के अतिरिक्त, समान एप्लिकेशन-विशिष्ट विधिें विकसित की गई हैं। उदाहरण के लिए, डीएमडी की तरह, प्रोनी की विधि [[नम साइनसॉइड]] के सुपरपोजिशन के रूप में सिग्नल का प्रतिनिधित्व करती है। जलवायु विज्ञान में, रैखिक व्युत्क्रम मॉडलिंग भी डीएमडी के साथ दृढ़ता से जुड़ा हुआ है।<ref>{{Cite journal|url=http://journals.ametsoc.org/doi/abs/10.1175/1520-0442(1993)006%3C1067:PONSST%3E2.0.CO%3B2| title=Prediction of Niño 3 Sea Surface Temperatures Using Linear Inverse Modeling|last=Penland, Magorian|first=Cecile, Theresa|date=1993|journal=J. Climate|volume=6| issue=6| page=1067| doi=10.1175/1520-0442(1993)006<1067:PONSST>2.0.CO;2| bibcode=1993JCli....6.1067P}}</ref> अधिक विस्तृत सूची के लिए, तू एट अल देखें।<ref name=tu2014 /> | ||
== उदाहरण == | == उदाहरण == | ||
=== प्रोफ़ाइल का पिछला किनारा === | === प्रोफ़ाइल का पिछला किनारा === | ||
[[File:Joukowsky Karman Vortex Street s.png|thumb|300px|चित्र 1 ट्रेलिंग एज भंवर (एंट्रॉपी)]]प्रवाह में किसी बाधा के पीछे कार्मन वर्टेक्स स्ट्रीट विकसित हो सकता है। चित्र 1 में प्रोफ़ाइल के पीछे वर्टेक्स के बिछले समाप्ति का दृश्य दिखाया गया है। डीएमडी विश्लेषण को 90 क्रमिक Entropy फ़ील्ड्स (animated gif (1.9MB)) पर | [[File:Joukowsky Karman Vortex Street s.png|thumb|300px|चित्र 1 ट्रेलिंग एज भंवर (एंट्रॉपी)]]प्रवाह में किसी बाधा के पीछे कार्मन वर्टेक्स स्ट्रीट विकसित हो सकता है। चित्र 1 में प्रोफ़ाइल के पीछे वर्टेक्स के बिछले समाप्ति का दृश्य दिखाया गया है। डीएमडी विश्लेषण को 90 क्रमिक Entropy फ़ील्ड्स (animated gif (1.9MB)) पर प्रयुक्त किया गया और इससे अनुमानित इज़नवैल्यू-स्पेक्ट्रम प्राप्त हुआ है जैसा कि नीचे चित्रित है। इस विश्लेषण को संख्यात्मक परिणामों पर प्रयुक्त किया गया था, सरकारी समीक्षा से किसी संबंधित इकाइयों का संदर्भ नहीं किया गया था। प्रोफ़ाइल सफेद रंग में दिखाई देती है। सफेद वृत्ताकार हड्डियाँ प्रोसेसर सीमाएँ हैं क्योंकि गणना पैरलेल कंप्यूटर पर विभिन्न गणनात्मक ब्लॉक्स का उपयोग करके किया गया था। | ||
[[File:Joukowsky Karman Vortex Street Spectrum.png|center|300px|चित्र 2 डीएमडी-स्पेक्ट्रम]]स्पेक्ट्रम का लगभग तिहाई भाग अत्यधिक नम (बड़ा, नकारात्मक<math>\lambda_r</math> ) और उसे दिखाया नहीं गया है।प्रमुख वर्टेक्स शेडिंग मोड निम्नलिखित तस्वीरों में दिखाया गया है। बाईं तस्वीर वास्तविक भाग है, दाईं तस्वीर भावित भाग है। | [[File:Joukowsky Karman Vortex Street Spectrum.png|center|300px|चित्र 2 डीएमडी-स्पेक्ट्रम]]स्पेक्ट्रम का लगभग तिहाई भाग अत्यधिक नम (बड़ा, नकारात्मक<math>\lambda_r</math> ) और उसे दिखाया नहीं गया है।प्रमुख वर्टेक्स शेडिंग मोड निम्नलिखित तस्वीरों में दिखाया गया है। बाईं तस्वीर वास्तविक भाग है, दाईं तस्वीर भावित भाग है। | ||
Line 82: | Line 82: | ||
</math> | </math> | ||
यहाँ <math> x_1 </math> समस्या के किसी भी स्वतंत्र चर हो सकता है, | यहाँ <math> x_1 </math> समस्या के किसी भी स्वतंत्र चर हो सकता है, किन्तु इसे पहले से चुना जाना होता है। उदाहरण के रूप में इस पैटर्न को ले: | ||
:<math> | :<math> | ||
q(x,y,t)=e^{-i \omega t} \hat q (x,t) e^{-(y/b)^2} \Re \left\{ e^{i (k x - \omega t)} \right\} + \text{random noise} | q(x,y,t)=e^{-i \omega t} \hat q (x,t) e^{-(y/b)^2} \Re \left\{ e^{i (k x - \omega t)} \right\} + \text{random noise} | ||
Line 88: | Line 88: | ||
यहाँ समय को पूर्व-चयनित घनांकीय कारक माना जाता है। | यहाँ समय को पूर्व-चयनित घनांकीय कारक माना जाता है। | ||
निम्नलिखित चित्र में नमूना दिया गया है जहाँ <math> \omega = 2\pi /0.1 </math>, <math> b=0.02 </math> और <math> k = 2\pi/ b </math>हैं। बाईं तस्वीर बिना पैटर्न के दिखाती है, जबकि दाएं चित्र में | निम्नलिखित चित्र में नमूना दिया गया है जहाँ <math> \omega = 2\pi /0.1 </math>, <math> b=0.02 </math> और <math> k = 2\pi/ b </math>हैं। बाईं तस्वीर बिना पैटर्न के दिखाती है, जबकि दाएं चित्र में ध्वनि जोड़ा गया है। रैंडम ध्वनि का अधिकार के अम्पलीट्यूड नमूने के अम्पलीट्यूड के बराबर है। | ||
[[File:Q periodic.png|400px]][[File:Q periodic noise.png|400px]] समय डीएमडी विश्लेषण 21 सिंथेटिक रूप से उत्पन्न क्षेत्रों के साथ किया जाता है, जिसमें समय अंतराल<math> \Delta t =1/90\text{ s}</math>है, और विश्लेषण को सीमित किया जाता है <math> f =45\text{ Hz}</math>। | [[File:Q periodic.png|400px]][[File:Q periodic noise.png|400px]] समय डीएमडी विश्लेषण 21 सिंथेटिक रूप से उत्पन्न क्षेत्रों के साथ किया जाता है, जिसमें समय अंतराल<math> \Delta t =1/90\text{ s}</math>है, और विश्लेषण को सीमित किया जाता है <math> f =45\text{ Hz}</math>। | ||
[[File:Synth Spec.png|center|600px]]स्पेक्ट्रम सममित है और तीन लगभग अवमंदित मोड (छोटा नकारात्मक वास्तविक भाग) दिखाता है,जबकि अन्य मोड गहरे रूप से अवकाशित होते हैं। उनके अंकीय मान निम्नलिखित हैं: | [[File:Synth Spec.png|center|600px]]स्पेक्ट्रम सममित है और तीन लगभग अवमंदित मोड (छोटा नकारात्मक वास्तविक भाग) दिखाता है,जबकि अन्य मोड गहरे रूप से अवकाशित होते हैं। उनके अंकीय मान निम्नलिखित हैं: | ||
<math> \omega_1=-0.201, \omega_{2/3}=-0.223 \pm i 62.768</math> यहाँ, वास्तविक मोड में माध्यम है जो क्षेत्र का औसत है, जबकि <math> \omega_{2/3}</math> निर्दिष्ट पैटर्न से मेल खाता है जिसमें <math> f = 10\text{ Hz} </math> है। इससे अधिकतम त्रुटि का मात्रांशीय त्रुटि है -1/1000। | <math> \omega_1=-0.201, \omega_{2/3}=-0.223 \pm i 62.768</math> यहाँ, वास्तविक मोड में माध्यम है जो क्षेत्र का औसत है, जबकि <math> \omega_{2/3}</math> निर्दिष्ट पैटर्न से मेल खाता है जिसमें <math> f = 10\text{ Hz} </math> है। इससे अधिकतम त्रुटि का मात्रांशीय त्रुटि है -1/1000। ध्वनि को दस गुना बढ़ाकर भी लगभग ही त्रुटि होती है। निम्नलिखित चित्र में दोनों में से नीचे दिए गए मोड के वास्तविक और काल्पनिक भाग दिखाए गए हैं। | ||
[[File:Detected Eigenmode.png|center|600px]] | [[File:Detected Eigenmode.png|center|600px]] | ||
== यह भी देखें == | == यह भी देखें == | ||
प्रायोगिक डेटा के कई अन्य अपघटन | प्रायोगिक डेटा के कई अन्य अपघटन उपस्थित हैं। यदि शासी समीकरण उपलब्ध हैं, तो आइगेनवैल्यू अपघटन संभव हो सकता है। | ||
* आइगेनवैल्यू अपघटन | * आइगेनवैल्यू अपघटन | ||
* [[अनुभवजन्य मोड अपघटन]] | * [[अनुभवजन्य मोड अपघटन]] |
Revision as of 07:44, 3 August 2023
गतिशील मोड अपघटन (डीएमडी) आयामिकता घटाने का एल्गोरिदम है जिसका विकास 2008 में पीटर जे. श्मिड और जोर्न सेस्टरहेन ने किया था। डेटा की समय श्रृंखला के आधार पर, डीएमडी सेट के मोड की गणना करता है, जिनमें से प्रत्येक निश्चित दोलन आवृत्ति और क्षय/विकास दर से जुड़ा संबंधित होती है। विशेष रूप से रैखिक प्रणालियों के लिए, ये मोड और आवृत्तियाँ प्रणाली के सामान्य मोड के समानांतर होते हैं, किन्तु अधिक सासामान्य रूप से, वे रचना संचालक (जिसे कोपमैन ऑपरेटर भी कहा जाता है) के मोड और आइगेनवैल्यू के अनुमान होते हैं। प्रत्येक मोड से जुड़े आंतरिक अस्थायी व्यवहारों के कारण, डीएमडी प्रमुख घटक विश्लेषण जैसे आयामी कमी के तरीकों से भिन्न होता है, जो प्राधिकृती प्रतिस्थापन विधियों (प्रिंसिपल कॉम्पोनेंट एनालिसिस) के अतिरिक्त होते हैं। क्योंकि इसके मोड ऑर्थोगोनल नहीं हैं, डीएमडी-आधारित अभ्यावेदन पीसीए द्वारा उत्पन्न किए गए अभ्यावेदन की तुलना में कम सरल हो सकते हैं। चूंकि, वे अधिक शारीरिक रूप से सार्थक भी हो सकते हैं क्योंकि प्रत्येक मोड समय में नम (या संचालित) साइनसॉइडल व्यवहार से जुड़ा होता है।
अवलोकन
डायनेमिक मोड अपघटन को पहली बार श्मिड द्वारा फ्लो डेटा से गतिजंक्य विशेषताएं निकालने के लिए संख्यात्मक प्रक्रिया के रूप में प्रस्तुत किया गया था।[1]
डेटा स्नैपशॉट अनुक्रम के रूप में प्राप्त होता है, जहां वी1 से एन तक के लिए:
यहाँ व्युत्पन्न फ्लो फ़ील्ड की -प्रवाह स्नैपशॉट है, और डेटा मैट्रिक्स है जिसके कॉलम व्यक्तिगत स्नैपशॉट होते हैं। इन स्नैपशॉट कोरैखिक गणना के माध्यम से संबंधित होते हैं जो रैखिक गतिजंकीय प्रणाली को परिभाषित करता है:
- जो समय के सैंपलिंग अवधि के समय लगभग समान रहता है। मैट्रिक्स रूप में लिखें, इससे यह सिद्ध होता है कि:
यहाँ अवशेषों का वेक्टर है जो ऐसे व्यवहारों का कारण बनता है जिनका पूरी तरह से वर्णन नहीं किया जा सकता, है, , और चाहे किसी भी दृष्टिकोण से, डीएमडी का आउटपुट , के इगेनवैल्यूज और इगेनवेक्टर होते हैं, जिन्हें डीएमडी इगेनवैल्यूज और डीएमडी मोड के रूप में संदर्भित किया जाता है।
एल्गोरिथम
इन इगेनवैल्यूज और मोड्स को प्राप्त करने की दो विधियाँ हैं। पहला विधि अर्नोल्डी जैसा है, जो सिद्धांतिक विश्लेषण के लिए उपयुक्त है क्योंकि इसका क्रायलोव विधियों से संबंध होता है। दूसरा सिंगुलर वैल्यू विचलन (एसवीडी) पर आधारित विधि है जो डेटा में ध्वनि से अधिक ट्रुटियों और संख्यात्मक त्रुटियों के प्रति शक्तिशाली है।
अर्नोल्डी दृष्टिकोण
तरल पदार्थ अनुप्रयोगों में, स्नैपशॉट का आकार, , स्नैपशॉट की संख्या से बहुत अधिक बड़ा माना जाता है, इसलिए समान रूप से कई वैध विकल्प उपस्थित हैं। मूल डीएमडी एल्गोरिदम चुनता है जिससे प्रत्येक स्नैपशॉट में में के स्नैपशॉटों के रूप में लिनियर संयोजन के रूप में व्यक्त किया जा सकता है। क्योंकि अधिकांश स्नैपशॉट दोनों डेटा सेट में दिखाई देते हैं, इस प्रतिनिधित्व को स्नैपशॉट , को छोड़कर सभी स्नैपशॉटों के लिए त्रुटिरहित बनाता है, जो इस तरह से लिखा जाता है:
यहाँ गुणांकों का सेट है जिसे डीएमडी को पहचानना चाहिए और अवशेष है।
पूरे मिलाकर,
यहाँ साथी मैट्रिक्स है
सदिश न्यूनतम वर्ग समस्या को हल करके गणना की जा सकती है, जो समग्र अवशिष्ट को न्यूनतम करती है। विशेष रूप से यदि हम का QR विचलन लें, तो होता है।
इस रूप में, डीएमडी प्रकार का अर्नोल्डी पुनरावृत्ति है, और इसलिए के इगेनवैल्यूज के अनुमान होते हैं। इसके अतिरिक्त, यदि , का इगेनवेक्टर है, तब , का अउपस्थित इगेनवेक्टर होता है। यहां इगेनवैल्यूज विचलन के बजाए पर किया जाता है क्योंकि , से अधिक छोटा होता है, इसलिए डीएमडी की संगणना निवेश स्नैपशॉट के आकार के अतिरिक्त स्नैपशॉट की संख्या से निर्धारित होती है।
एसवीडी-आधारित दृष्टिकोण
साथी मैट्रिक्स की गणना करने के अतिरिक्त , एसवीडी-आधारित दृष्टिकोण मैट्रिक्स उत्पन्न करता है वह संबंधित है समानता परिवर्तन के माध्यम से. ऐसा करने के लिए, मान लें कि हमारे पास एसवीडी है . तब
अर्नोल्डी-आधारित दृष्टिकोण द्वारा बनाई गई धारणा के बराबर, हम चुनते हैं जैसे कि स्नैपशॉट में में स्तंभों के रैखिक सुपरपोजिशन के रूप में लिखा जा सकता है , जो इस आवश्यकता के बराबर है कि उन्हें प्रमुख घटक विश्लेषण के सुपरपोजिशन के रूप में लिखा जा सकता है। इस प्रतिबंध के साथ, अवशिष्ट को कम करने के लिए आवश्यक है कि यह POD आधार पर ऑर्थोगोनल हो (अर्थात्, )। फिर समीकरण को से दोनों तरफ़ से गुणा करने से इसे मिलाया जा सकता है: , जिसे स्पष्ट किया जा सकता है:
क्योंकि और समानता संवेदी रूप से जुड़े हुए हैं, इसलिए के इगेनवैल्यूज के इगेनवैल्यूज होते हैं, और यदि , का इगेनवेक्टर है , तब , का इगेनवेक्टर होता है। .
संक्षेप में, एसवीडी-आधारित विधि का अनुसरण निम्नलिखित तरीके से किया जाता है:
- डेटा की समय श्रृंखला को विभाजित करें दो आव्यूहों में और .
- के एसवीडी की गणना करें .
- मैट्रिक्स बनाएं , और इसके इगेनवैल्यूज की गणना करें और इगेनवेक्टर्स .
- th>-th डीएमडी इगेनवैल्यूज है और -वें डीएमडी मोड है .
अर्नोल्डी-जैसे दृष्टिकोण पर एसवीडी-आधारित दृष्टिकोण का लाभ यह है कि डेटा में ध्वनि और संख्यात्मक ट्रंकेशन मुद्दों की भरपाई एसवीडी को काटकर की जा सकती है। . जैसा कि उल्लेख किया गया है [1]इस ट्रंकेशन चरण के बिना प्रयोगात्मक डेटा सेट पर पहले जोड़े मोड और आइगेनवैल्यू से अधिक की स्पष्ट गणना करना कठिनाई हो सकता है।
सैद्धांतिक और एल्गोरिथम प्रगति
इसकी प्रारंभिकता से 2010 में, डीएमडी को समझने और सुधारने पर अधिक काम किया गया है।[2] Rowley एवं उनके सहयोगियों ने डीएमडी का पहला विश्लेषण किया जिसमें डीएमडी और Koopman ऑपरेटर के बीच संबंध स्थापित किया गया और यह समझने में सहायता की जब डीएमडी को गैर-लीनियर प्रणालियों पर प्रयुक्त किया गया। इसके बाद से, कई संशोधन विकसित किए गए हैं जो इस संबंध को और मज़बूत करते हैं या इस दृष्टिकोन को अधिक उत्कृष्ट और व्यापक बनाने में सहायता करते हैं।
- अनुकूलित डीएमडी: अनुकूलित डीएमडी मूल डीएमडी एल्गोरिदम का संशोधन है जिसे उस दृष्टिकोण की दो सीमाओं का समाधान करने के लिए डिज़ाइन किया गया है: (i) डीएमडी मोड चयन की कठिनाई, और (ii) डीएमडी की संवेदनशीलता ध्वनि या समय शृंखला के अंतिम स्नैपशॉट में अन्य त्रुटियों के प्रति।[3] ऑप्टिमाइज़्ड डीएमडी डीएमडी प्रक्रिया को अनुकूलन समस्या के रूप में पुनर्स्थापित करता है जहाँ पहचानी गई लीनियर ऑपरेटर की निश्चित श्रेणी होती है। इसके अतिरिक्त, जिस तरह डीएमडी ने पिछले सभी स्नैपशॉट को पूरी तरह से पुनर्जन्मित किया है, ऑप्टिमाइज़्ड डीएमडी द्वारा रीकंस्ट्रक्शन त्रुटियों को डेटा सेट में वितरित किया जा सकता है, जिससे यह विधि व्यावसायिक रूप से अधिक शक्तिशाली दिखता है।
- इष्टतम मोड अपघटन: इष्टतम मोड अपघटन (ओएमडी) डीएमडी प्रक्रिया को अनुकूलन समस्या के रूप में पुनर्गठित करता है और उपयोगकर्ता को सीधे पहचाने गए प्रणाली की रैंक प्रयुक्त करने की अनुमति देता है।[4] यदि यह श्रेणी योग्य ढंग से चुनी जाए, तो OMD सांदर्भिक और प्रायोगिक डेटा सेटों पर छोटे शेष त्रुटि और अधिक स्पष्ट इजेनवैल्यूज़ के साथ लीनियर मॉडल उत्पन्न कर सकता है।
- स्पष्ट डीएमडी: स्पष्ट डीएमडी एल्गोरिदम मूल डीएमडी एल्गोरिदम को दो तरीकों से सामान्यीकृत करता है। सबसे पहले, मूल डीएमडी एल्गोरिदम में डेटा स्नैपशॉट की समय श्रृंखला होनी चाहिए, किन्तु स्पष्ट डीएमडी स्नैपशॉट जोड़े के डेटा सेट को स्वीकार करता है।[5] जोड़ी में स्नैपशॉट को निश्चित द्वारा अलग किया जाना चाहिए , किन्तु ही समय श्रृंखला से निकालने की आवश्यकता नहीं है। विशेष रूप से, स्पष्ट डीएमडी कई प्रयोगों के डेटा को ही डेटा सेट में एकत्रित करने की अनुमति दे सकता है। दूसरा, मूल डीएमडी एल्गोरिदम पीओडी मोड के सेट पर प्रोजेक्ट करके डेटा को प्रभावी ढंग से प्री-प्रोसेस करता है। स्पष्ट डीएमडी एल्गोरिदम इस प्री-प्रोसेसिंग चरण को हटा देता है, और डीएमडी मोड का उत्पादन कर सकता है जिसे पीओडी मोड के सुपरपोजिशन के रूप में नहीं लिखा जा सकता है।
- स्पार्सिटी प्रमोशन डीएमडी: स्पार्सिटी प्रोमोटिंग डीएमडी, डीएमडी मोड और इजेनवैल्यू चयन के लिए पोस्ट प्रोसेसिंग प्रक्रिया है।[6] स्पार्सिटी प्रोमोटिंग डीएमडी पेनल्टी का उपयोग करता है जो महत्वपूर्ण डीएमडी मोड्स का छोटा सेट पहचानने के लिए उपयुक्त होता है, और यह डीएमडी मोड चयन समस्या के वैकल्पिक दृष्टिकोन है जिसे संवर्धित लैग्रेंजियन विधि का उपयोग करके कुशलतापूर्वक हल किया जा सकता है।
- मल्टी-रिज़ॉल्यूशन डीएमडी: मल्टी-रिज़ॉल्यूशन डीएमडी (एमआरडीएमडी) स्पष्ट डीएमडी के साथ मल्टी-रिज़ॉल्यूशन विश्लेषण में उपयोग की जाने वाली विधिों का संयोजन है, जिसे कई टाइमस्केल वाले डेटा सेट से डीएमडी मोड और आइजेनवैल्यू को शक्तिशाली करने के लिए डिज़ाइन किया गया है।[7] एमआरडीएमडी दृष्टिकोण वैश्विक सतह तापमान डेटा पर प्रयुक्त किया गया था, और एल नीनो वर्षों के समय दिखाई देने वाले डीएमडी मोड की पहचान करता है।
- विस्तारित डीएमडी: विस्तारित डीएमडी स्पष्ट डीएमडी का संशोधन है जो डीएमडी और कूपमैन ऑपरेटर के बीच संबंध को शक्तिशाली करता है।[8] जैसा कि नाम से पता चलता है, विस्तारित डीएमडी डीएमडी का विस्तार है जो कोपमैन ऑपरेटर के अधिक स्पष्ट अनुमान उत्पन्न करने के लिए अवलोकन योग्य कार्यों के समृद्ध सेट का उपयोग करता है। इस विस्तारित सेट को प्राथमिकता से चुना जा सकता है या डेटा से सीखा जा सकता है।[9][10] इसने यह भी प्रदर्शित किया कि डीएमडी और संबंधित विधियाँ अधिक सामान्यतः उपयोग किए जाने वाले आइगेनवैल्यू और मोड के अतिरिक्त कोपमैन आइजेनफंक्शन के सन्निकटन का उत्पादन करती हैं।
- नियंत्रण के साथ डीएमडी: नियंत्रण के साथ डायनामिक मोड अपघटन (डीएमडीसी) [11] इनपुट आउटपुट प्रणाली से प्राप्त डेटा के लिए डिज़ाइन की गई डीएमडी प्रक्रिया का संशोधन है। डीएमडीसी की अनूठी विशेषता ओपन लूप डायनेमिक्स से प्रणाली एक्चुएशन के प्रभावों को स्पष्ट करने की क्षमता है, जो तब उपयोगी होती है जब एक्चुएशन की उपस्थिति में डेटा प्राप्त किया जाता है।
- कुल न्यूनतम वर्ग डीएमडी: कुल न्यूनतम वर्ग डीएमडी हाल ही की बदलाव है जो एक्सैक्ट डीएमडी की संशोधन है, जिसका उद्देश्य डेटा में मापन ध्वनि के प्रति स्थायित्व से संबंधित मुद्दों का सामना करना है।[12] में, लेखक एक्सैक्ट डीएमडी को प्रेक्षण समस्या के रूप में व्याख्या करते हैं जो साधारण लीस्ट स्क्वेयर्स (ओएलएस) का उपयोग करके हल किया जाता है, जो इसमें अनुमानित करता है कि प्रेक्षक ध्वनियुक्त मुक्त हैं। यह धारणा डीएमडी आइगेनवैल्यू में पूर्वाग्रह उत्पन्न करती है जब इसे प्रयोगात्मक डेटा सेट पर प्रयुक्त किया जाता है जहां सभी अवलोकन ध्वनि होते हैं। कुल न्यूनतम वर्ग डीएमडी ओएलएस समस्या को कुल न्यूनतम वर्ग से प्रतिस्थापित करता है, जो इस पूर्वाग्रह को समाप्त करता है।
- गतिशील वितरण अपघटन: डीडीडी निरंतर समय में आगे की समस्या, अर्थात स्थानांतरण ऑपरेटर पर ध्यान केंद्रित करता है। चूंकि विकसित विधि का उपयोग निरंतर समय में डीएमडी समस्याओं को ठीक करने के लिए भी किया जा सकता है।[13]
यहां सूचीबद्ध एल्गोरिदम के अतिरिक्त, समान एप्लिकेशन-विशिष्ट विधिें विकसित की गई हैं। उदाहरण के लिए, डीएमडी की तरह, प्रोनी की विधि नम साइनसॉइड के सुपरपोजिशन के रूप में सिग्नल का प्रतिनिधित्व करती है। जलवायु विज्ञान में, रैखिक व्युत्क्रम मॉडलिंग भी डीएमडी के साथ दृढ़ता से जुड़ा हुआ है।[14] अधिक विस्तृत सूची के लिए, तू एट अल देखें।[5]
उदाहरण
प्रोफ़ाइल का पिछला किनारा
प्रवाह में किसी बाधा के पीछे कार्मन वर्टेक्स स्ट्रीट विकसित हो सकता है। चित्र 1 में प्रोफ़ाइल के पीछे वर्टेक्स के बिछले समाप्ति का दृश्य दिखाया गया है। डीएमडी विश्लेषण को 90 क्रमिक Entropy फ़ील्ड्स (animated gif (1.9MB)) पर प्रयुक्त किया गया और इससे अनुमानित इज़नवैल्यू-स्पेक्ट्रम प्राप्त हुआ है जैसा कि नीचे चित्रित है। इस विश्लेषण को संख्यात्मक परिणामों पर प्रयुक्त किया गया था, सरकारी समीक्षा से किसी संबंधित इकाइयों का संदर्भ नहीं किया गया था। प्रोफ़ाइल सफेद रंग में दिखाई देती है। सफेद वृत्ताकार हड्डियाँ प्रोसेसर सीमाएँ हैं क्योंकि गणना पैरलेल कंप्यूटर पर विभिन्न गणनात्मक ब्लॉक्स का उपयोग करके किया गया था।
स्पेक्ट्रम का लगभग तिहाई भाग अत्यधिक नम (बड़ा, नकारात्मक ) और उसे दिखाया नहीं गया है।प्रमुख वर्टेक्स शेडिंग मोड निम्नलिखित तस्वीरों में दिखाया गया है। बाईं तस्वीर वास्तविक भाग है, दाईं तस्वीर भावित भाग है।
File:Joukowsky Karman Vortex Street EV real.png | File:Joukowsky Karman Vortex Street EV imag.png |
फिर से, इस चित्र में एन्ट्रॉपी-ईजेनवेक्टर दिखाया गया है। इसी मोड की ध्वनिक सामग्री अगले प्लॉट के नीचे के आधार में दिखाई देती है। शीर्ष भाग ऊपर के भाग के समान एंट्रोपी मोड से संबंधित है।
यात्रा पैटर्न का सिंथेटिक उदाहरण
डीएमडी विश्लेषण मानता है कि ऐसे पैटर्न की संरचना होती है जिसका रूप है
यहाँ समस्या के किसी भी स्वतंत्र चर हो सकता है, किन्तु इसे पहले से चुना जाना होता है। उदाहरण के रूप में इस पैटर्न को ले:
यहाँ समय को पूर्व-चयनित घनांकीय कारक माना जाता है।
निम्नलिखित चित्र में नमूना दिया गया है जहाँ , और हैं। बाईं तस्वीर बिना पैटर्न के दिखाती है, जबकि दाएं चित्र में ध्वनि जोड़ा गया है। रैंडम ध्वनि का अधिकार के अम्पलीट्यूड नमूने के अम्पलीट्यूड के बराबर है।
समय डीएमडी विश्लेषण 21 सिंथेटिक रूप से उत्पन्न क्षेत्रों के साथ किया जाता है, जिसमें समय अंतरालहै, और विश्लेषण को सीमित किया जाता है ।
स्पेक्ट्रम सममित है और तीन लगभग अवमंदित मोड (छोटा नकारात्मक वास्तविक भाग) दिखाता है,जबकि अन्य मोड गहरे रूप से अवकाशित होते हैं। उनके अंकीय मान निम्नलिखित हैं:
यहाँ, वास्तविक मोड में माध्यम है जो क्षेत्र का औसत है, जबकि निर्दिष्ट पैटर्न से मेल खाता है जिसमें है। इससे अधिकतम त्रुटि का मात्रांशीय त्रुटि है -1/1000। ध्वनि को दस गुना बढ़ाकर भी लगभग ही त्रुटि होती है। निम्नलिखित चित्र में दोनों में से नीचे दिए गए मोड के वास्तविक और काल्पनिक भाग दिखाए गए हैं।
यह भी देखें
प्रायोगिक डेटा के कई अन्य अपघटन उपस्थित हैं। यदि शासी समीकरण उपलब्ध हैं, तो आइगेनवैल्यू अपघटन संभव हो सकता है।
- आइगेनवैल्यू अपघटन
- अनुभवजन्य मोड अपघटन
- वैश्विक मोड
- सामान्य मोड
- उचित ओर्थोगोनल अपघटन
- विलक्षण मान अपघटन
संदर्भ
- ↑ 1.0 1.1 P.J. Schmid. "Dynamic mode decomposition of numerical and experimental data." Journal of Fluid Mechanics 656.1 (2010): 5–28.
- ↑ C.W. Rowley, I Mezic, S. Bagheri, P. Schlatter, and D.S. Henningson, "Spectral analysis of nonlinear flows." Journal of Fluid Mechanics 641 (2009): 85-113
- ↑ K.K. Chen, J.H. Tu, and C.W. Rowley, "Variants of dynamic mode decomposition: boundary condition, Koopman, and Fourier analyses." Journal of Nonlinear Science 22 (2012): 887-915.
- ↑ A. Wynn, D. S. Pearson, B. Ganapathisubramani and P. J. Goulart, "Optimal mode decomposition for unsteady flows." Journal of Fluid Mechanics 733 (2013): 473-503
- ↑ 5.0 5.1 Tu, Rowley, Luchtenburg, Brunton, and Kutz (December 2014). "On Dynamic Mode Decomposition: Theory and Applications". American Institute of Mathematical Sciences. 1 (2): 391–421. arXiv:1312.0041. doi:10.3934/jcd.2014.1.391. S2CID 46419148.
{{cite journal}}
: CS1 maint: multiple names: authors list (link) - ↑ M.R. Jovanovic, P.J. Schmid, and J.W. Nichols, "Sparsity-promoting dynamic mode decomposition." Physics of Fluids 26 (2014)
- ↑ J.N. Kutz, X. Fu, and S.L. Brunton, "Multi-resolution dynamic mode decomposition." arXiv preprint arXiv:1506.00564 (2015).
- ↑ M.O. Williams , I.G. Kevrekidis, C.W. Rowley, "A Data–Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition." Journal of Nonlinear Science 25 (2015): 1307-1346.
- ↑ Li, Qianxiao; Dietrich, Felix; Bollt, Erik M.; Kevrekidis, Ioannis G (2018). "Extended dynamic mode decomposition with dictionary learning: A data-driven adaptive spectral decomposition of the Koopman operator". Chaos: An Interdisciplinary Journal of Nonlinear Science. 27 (10): 103111. arXiv:1707.00225. doi:10.1063/1.4993854. PMID 29092410. S2CID 41957686.
- ↑ Gulina, Marvyn; Mauroy, Alexandre (2021). "जलाशय कंप्यूटर के साथ कूपमैन ऑपरेटर का अनुमान लगाने की दो विधियाँ". Chaos: An Interdisciplinary Journal of Nonlinear Science. 31 (2): 023116. arXiv:2008.10263. Bibcode:2021Chaos..31b3116G. doi:10.1063/5.0026380. PMID 33653036. S2CID 221266743.
- ↑ J.L. Proctor, S.L. Brunton, and J.N. Kutz, "Dynamic mode decomposition with control." arXiv preprint arXiv:1409.6358 (2014).
- ↑ M.S. Hemati, C.W. Rowley, E.A. Deem, and L.N. Cattafesta, "De-Biasing the Dynamic Mode Decomposition for Applied Koopman Spectral Analysis of Noisy Datasets." arXiv preprint arXiv:1502.03854 (2015).
- ↑ Taylor-King, Jake P.; Riseth, Asbjørn N.; Macnair, Will; Claassen, Manfred (2020-01-10). "एकल-सेल स्नैपशॉट समय श्रृंखला के लिए गतिशील वितरण अपघटन आईपीएससी रीप्रोग्रामिंग के दौरान उप-जनसंख्या और प्रक्षेपवक्र की पहचान करता है". PLOS Computational Biology (in English). 16 (1): e1007491. Bibcode:2020PLSCB..16E7491T. doi:10.1371/journal.pcbi.1007491. ISSN 1553-7358. PMC 6953770. PMID 31923173.
- ↑ Penland, Magorian, Cecile, Theresa (1993). "Prediction of Niño 3 Sea Surface Temperatures Using Linear Inverse Modeling". J. Climate. 6 (6): 1067. Bibcode:1993JCli....6.1067P. doi:10.1175/1520-0442(1993)006<1067:PONSST>2.0.CO;2.
{{cite journal}}
: CS1 maint: multiple names: authors list (link)
- Schmid, P. J. & Sesterhenn, J. L. 2008 Dynamic mode decomposition of numerical and experimental data. In Bull. Amer. Phys. Soc., 61st APS meeting, p. 208. San Antonio.
- Hasselmann, K., 1988. POPs and PIPs. The reduction of complex dynamical systems using principal oscillation and interaction patterns. J. Geophys. Res., 93(D9): 10975–10988.