एन्सेम्बल कलमैन फ़िल्टर: Difference between revisions
(Created page with "एन्सेम्बल कल्मन फ़िल्टर (EnKF) एक पुनरावर्ती फ़िल्टर है जो बड़ी संख...") |
m (7 revisions imported from alpha:एन्सेम्बल_कलमैन_फ़िल्टर) |
||
(6 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
एन्सेम्बल | '''एन्सेम्बल कलमैन फ़िल्टर''' (ईएनकेएफ) [[पुनरावर्ती फ़िल्टर|रिकर्सिव फ़िल्टर]] है जो बड़ी संख्या में वैरिएबल वाली समस्याओं के लिए उपयुक्त है, जैसे कि भूभौतिकीय मॉडल में आंशिक अंतर समीकरणों का विवेकीकरण ईएनकेएफ की उत्पत्ति बड़ी समस्याओं के लिए [[कलमन फ़िल्टर]] के संस्करण के रूप में हुई (अनिवार्य रूप से, सहसंयोजकता आव्यूह को [[नमूना सहप्रसरण मैट्रिक्स|प्रारूपिक सहसंयोजकता]] आव्यूह द्वारा प्रतिस्थापित किया जाता है), और यह अब संयोजन पूर्वानुमान का महत्वपूर्ण [[डेटा आत्मसात|डेटा एसीमिलेसन]] घटक है। ईएनकेएफ [[कण फिल्टर|पार्टिकल फिल्टर]] से संबंधित है (इस संदर्भ में, पार्टिकल समूह सदस्य के समान है) किन्तु ईएनकेएफ यह धारणा बनाता है कि इसमें सम्मिलित सभी संभाव्यता वितरण [[सामान्य वितरण|गाऊसी]] हैं; जब यह प्रयुक्त होता है, जिससे यह पार्टिकल फिल्टर की तुलना में बहुत अधिक उत्तम होता है। | ||
==परिचय== | ==परिचय== | ||
एन्सेम्बल | एन्सेम्बल कलमैन फ़िल्टर (ईएनकेएफ) [[बायेसियन अनुमान|बायेसियन अपडेट]] समस्या का [[मोंटे कार्लो विधि]] कार्यान्वयन है: मॉडल किए गए प्रणाली की स्थिति ([[पूर्व संभावना|प्रियोर]], जिसे अधिकांशतः भूविज्ञान में पूर्वानुमान कहा जाता है) और डेटा संभावना की संभाव्यता घनत्व फलन (पीडीएफ) दी गई है। बेयस प्रमेय का उपयोग डेटा संभावना को ध्यान में रखने के पश्चात पीडीएफ प्राप्त करने के लिए किया जाता है (पोस्टीरियर संभावना, जिसे अधिकांशतः विश्लेषण कहा जाता है)। इसे बायेसियन अपडेट कहा जाता है। बायेसियन अपडेट को समय-समय पर नए डेटा को सम्मिलित करते हुए मॉडल को आगे बढ़ाने के साथ जोड़ा जाता है। मूल कलमैन फ़िल्टर, 1960 में प्रस्तुत किया गया था,<ref name="Kalman-1960-NAL">{{cite journal |first=R. E. |last=Kalman |s2cid=1242324 |title=रेखीय छानने और भविष्यवाणी की समस्याओं के लिए एक नया दृष्टिकोण|journal= Journal of Basic Engineering|volume=82 |year=1960 |issue=1 |pages=35–45 |doi=10.1115/1.3662552 }}</ref> मानता है कि सभी पीडीएफ सामान्य वितरण (गॉसियन धारणा) हैं और बायेसियन अपडेट द्वारा माध्य और सहसंयोजकता आव्यूह के परिवर्तन के लिए बीजगणितीय सूत्र प्रदान करता है, साथ ही समय में माध्य और सहसंयोजकता को आगे बढ़ाने के लिए सूत्र प्रदान करता है, किन्तु प्रणाली रैखिक हो। चूंकि, उच्च-आयामी प्रणालियों के लिए सहसंयोजकता आव्यूह को बनाए रखना कम्प्यूटेशनल रूप से संभव नहीं है। इस कारण से, ईएनकेएफएस विकसित किए गए थे।<ref name="Evensen-1994-SDA">{{cite journal |first=G. |last=Evensen |title=त्रुटि आँकड़ों का पूर्वानुमान लगाने के लिए मोंटे कार्लो विधियों का उपयोग करके गैर-रेखीय अर्ध-जियोस्ट्रोफिक मॉडल के साथ अनुक्रमिक डेटा आत्मसात|journal=Journal of Geophysical Research |volume=99 |issue=C5 |year=1994 |pages=143–162 |doi=10.1029/94JC00572 |hdl=1956/3035 |bibcode=1994JGR....9910143E |hdl-access=free }}</ref><ref name="Houtekamer-1998-DAE">{{cite journal |first1=P. |last1=Houtekamer |first2=H. L. |last2=Mitchell |title=एन्सेम्बल कलमैन फ़िल्टर तकनीक का उपयोग करके डेटा आत्मसात करना|journal=[[Monthly Weather Review]] |volume=126 |year=1998 |issue= 3|pages=796–811 |doi=10.1175/1520-0493(1998)126<0796:DAUAEK>2.0.CO;2 |citeseerx=10.1.1.3.1706 |bibcode=1998MWRv..126..796H }}</ref> ईएनकेएफएस स्टेट वैक्टरों के संग्रह का उपयोग करके प्रणाली स्थिति के वितरण का प्रतिनिधित्व करते हैं, जिसे संख्यात्मक मौसम पूर्वानुमान या एनसेम्बल कहा जाता है, और सहसंयोजकता आव्यूह को समुच्चय से गणना किए गए [[नमूना सहप्रसरण|प्रारूपिक सहसंयोजकता]] द्वारा प्रतिस्थापित करते हैं। इस प्रकार समूह को ऐसे संचालित किया जाता है जैसे कि यह [[यादृच्छिक नमूना|यादृच्छिक प्रारूपिक]] हो, किन्तु समूह के सदस्य वास्तव में [[सांख्यिकीय स्वतंत्रता]] नहीं रखते हैं, क्योंकि वह सभी ईएनकेएफ संगठित करते हैं। ईएनकेएफएस का लाभ यह है कि पीडीएफ को समय पर आगे बढ़ाने का कार्य केवल समूह के प्रत्येक सदस्य को आगे बढ़ाना है।<ref name="Evensen-2007-DAE">For a survey of EnKF and related data assimilation techniques, see {{cite book |first=G. |last=Evensen |title=Data Assimilation : The Ensemble Kalman Filter |publisher=Springer |location=Berlin |year=2007 |isbn=978-3-540-38300-0 }}</ref> | ||
==व्युत्पत्ति== | ==व्युत्पत्ति== | ||
===कलमैन फ़िल्टर=== | ===कलमैन फ़िल्टर=== | ||
{{main| | {{main|कलमैन फ़िल्टर }} | ||
मान लीजिए कि <math>\mathbf{x}</math> एक मॉडल के <math>n</math>-आयामी स्टेट वेक्टर को दर्शाता है, और मान लेता है कि इसमें माध्य <math>\mathbf{\mu}</math> और सहसंयोजकता <math>Q</math> के साथ गॉसियन संभाव्यता वितरण है, अर्थात, इसका पीडीएफ है | |||
:<math> p(\mathbf{x})\propto\exp\left( -\frac{1}{2}(\mathbf{x}-\mathbf{\mu })^{\mathrm{T}}Q^{-1}(\mathbf{x}-\mathbf{\mu})\right) . </math> | :<math> p(\mathbf{x})\propto\exp\left( -\frac{1}{2}(\mathbf{x}-\mathbf{\mu })^{\mathrm{T}}Q^{-1}(\mathbf{x}-\mathbf{\mu})\right) . </math> | ||
यहां और नीचे, <math>\propto</math> का अर्थ आनुपातिक है; एक पीडीएफ को सदैव स्केल किया जाता है जिससे पूर्ण स्थान पर इसका अभिन्न अंग एक होता है। यह <math>p(\mathbf{x})</math>, जिसे पूर्व कहा जाता है, मॉडल को चलाकर समय पर विकसित किया गया था और अब इसे नए डेटा के लिए अपडेट किया जाना है। यह मान लेना स्वाभाविक है कि डेटा का त्रुटि वितरण ज्ञात है; डेटा को त्रुटि अनुमान के साथ आना होगा, अन्यथा वह अर्थहीन हैं। यहां, डेटा <math>\mathbf{d}</math> को सहसंयोजकता R और माध्य <math>H\mathbf{x}</math> के साथ गॉसियन पीडीएफ माना जाता है, जहां <math>H</math> तथाकथित अवलोकन आव्यूह है। सहसंयोजकता आव्यूह आर डेटा की त्रुटि के अनुमान का वर्णन करता है; यदि डेटा वेक्टर <math>\mathbf{d}</math> की प्रविष्टियों में यादृच्छिक त्रुटियां स्वतंत्र हैं, तो <math>R</math> विकर्ण है और इसकी विकर्ण प्रविष्टियां डेटा वेक्टर की संबंधित प्रविष्टियों की त्रुटि के मानक विचलन ("त्रुटि आकार") के वर्ग हैं <math>\mathbf{d}</math> मान <math>H\mathbf{x}</math> वह है जो डेटा त्रुटियों के अभाव में स्टेट <math>\mathbf{x}</math> के लिए डेटा का मान होगा। फिर प्रणाली स्थिति <math>\mathbf{x}</math> की नियमबद्ध डेटा <math>\mathbf{d}</math> की संभाव्यता घनत्व <math>p(\mathbf{d}|\mathbf{x})</math>, जिसे डेटा संभावना कहा जाता है, है | |||
:<math> p\left( \mathbf{d}|\mathbf{x}\right) \propto\exp\left( -\frac{1}{2}(\mathbf{d}-H\mathbf{x})^{\mathrm{T}}R^{-1}(\mathbf{d}-H\mathbf{x})\right) . </math> | :<math> p\left( \mathbf{d}|\mathbf{x}\right) \propto\exp\left( -\frac{1}{2}(\mathbf{d}-H\mathbf{x})^{\mathrm{T}}R^{-1}(\mathbf{d}-H\mathbf{x})\right) . </math> | ||
इस प्रकार स्थिति की पीडीएफ और डेटा संभावना को बेयस प्रमेय द्वारा डेटा <math>\mathbf{d}</math> के मूल्य पर नियमबद्ध प्रणाली स्थिति <math>\mathbf{x}</math> की नई संभाव्यता घनत्व देने के लिए संयोजित किया गया है, | |||
:<math> p\left( \mathbf{x}|\mathbf{d}\right) \propto p\left( \mathbf{d}|\mathbf{x}\right) p(\mathbf{x}). </math> | :<math> p\left( \mathbf{x}|\mathbf{d}\right) \propto p\left( \mathbf{d}|\mathbf{x}\right) p(\mathbf{x}). </math> | ||
आंकड़ा <math>\mathbf{d}</math> एक बार प्राप्त होने के | इस प्रकार आंकड़ा <math>\mathbf{d}</math> एक बार प्राप्त होने के पश्चात यह निश्चित हो जाता है, इसलिए पिछली स्थिति को इससे निरूपित करें <math>\mathbf{\hat{x}}</math> के अतिरिक्त <math>\mathbf{x}|\mathbf{d}</math> और पिछला पीडीएफ़ द्वारा <math>p\left( \mathbf{\hat{x}}\right) </math> इसे बीजगणितीय जोड़ द्वारा दिखाया जा सकता है <ref name="Anderson-1979-OF">{{cite book |first1=B. D. O. |last1=Anderson |first2=J. B. |last2=Moore |title=इष्टतम फ़िल्टरिंग|publisher=Prentice-Hall |location=Englewood Cliffs, NJ |year=1979 |isbn=978-0-13-638122-8 }}</ref> पिछला पीडीएफ भी गॉसियन है, | ||
:<math> p\left( \mathbf{\hat{x}}\right) \propto\exp\left( -\frac{1}{2}(\mathbf{\hat{x}}-\mathbf{\hat{\mu}})^{\mathrm{T}}\hat{Q}^{-1}(\mathbf{\hat{x}}-\mathbf{\hat{\mu}})\right) , </math> | :<math> p\left( \mathbf{\hat{x}}\right) \propto\exp\left( -\frac{1}{2}(\mathbf{\hat{x}}-\mathbf{\hat{\mu}})^{\mathrm{T}}\hat{Q}^{-1}(\mathbf{\hat{x}}-\mathbf{\hat{\mu}})\right) , </math> | ||
पश्च माध्य के साथ | इस प्रकार कलमैन अपडेट फ़ार्मुलों द्वारा दिए गए पश्च माध्य @ और सहसंयोजकता @ के साथ किया जाता है | ||
:<math> \mathbf{\hat{\mu}}=\mathbf{\mu}+K\left( \mathbf{d}-H\mathbf{\mu}\right) ,\quad\hat{Q}=\left( I-KH\right) Q, </math> | :<math> \mathbf{\hat{\mu}}=\mathbf{\mu}+K\left( \mathbf{d}-H\mathbf{\mu}\right) ,\quad\hat{Q}=\left( I-KH\right) Q, </math> | ||
जहाँ | |||
:<math> K=QH^{\mathrm{T}}\left( HQH^{\mathrm{T}}+R\right) ^{-1}</math> | :<math> K=QH^{\mathrm{T}}\left( HQH^{\mathrm{T}}+R\right) ^{-1}</math> | ||
तथाकथित कलमन फिल्टर | तथाकथित कलमन फिल्टर कलमैन लाभ व्युत्पत्ति आव्यूह है। | ||
=== | ===एन्सेम्बल कलमैन फ़िल्टर === | ||
ईएनकेएफ कलमैन फ़िल्टर का मोंटे कार्लो सन्निकटन है, जो स्टेट वेक्टर <math>\mathbf{x}</math> के पीडीएफ के सहसंयोजकता आव्यूह को विकसित करने से बचाता है . इसके अतिरिक्त, पीडीएफ को समूह द्वारा दर्शाया जाता है | |||
:<math> X=\left[ \mathbf{x}_{1},\ldots,\mathbf{x}_{N}\right] =\left[ \mathbf{x}_{i}\right]. </math> | :<math> X=\left[ \mathbf{x}_{1},\ldots,\mathbf{x}_{N}\right] =\left[ \mathbf{x}_{i}\right]. </math> | ||
डेटा | |||
इस प्रकार X एक <math>n\times N</math> आव्यूह है जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से एक प्रारूप बनाएंगे। चूँकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से स्वतंत्र नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ चरण उन्हें एक साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वह वास्तव में स्वतंत्र होंते है। | |||
डेटा <math>\mathbf{d}</math> को <math>m\times N</math> आव्यूह में दोहराएं जाते है | |||
:<math> D=\left[ \mathbf{d}_{1},\ldots,\mathbf{d}_{N}\right] =\left[ \mathbf{d}_{i}\right], \quad \mathbf{d}_{i}=\mathbf{d}+\mathbf{\epsilon_{i}}, \quad \mathbf{\epsilon_{i}} \sim N(0,R), </math> | :<math> D=\left[ \mathbf{d}_{1},\ldots,\mathbf{d}_{N}\right] =\left[ \mathbf{d}_{i}\right], \quad \mathbf{d}_{i}=\mathbf{d}+\mathbf{\epsilon_{i}}, \quad \mathbf{\epsilon_{i}} \sim N(0,R), </math> | ||
जिससे प्रत्येक कॉलम <math>\mathbf{d}_{i}</math> में डेटा वेक्टर <math>\mathbf{d}</math> प्लस एम-आयामी सामान्य वितरण <math>N(0,R)</math> से एक यादृच्छिक वेक्टर सम्मिलित होते है। यदि इसके अतिरिक्त, <math>X</math> के कॉलम पूर्व संभाव्यता वितरण से एक प्रारूप हैं, | |||
:<math> \hat{X}=X+K(D-HX) </math> | :<math> \hat{X}=X+K(D-HX) </math> | ||
पश्च संभाव्यता वितरण से एक | इस प्रकार पश्च संभाव्यता वितरण से एक प्रारूप बनाएं। इसे <math>H=1</math> के साथ अदिश स्थिति में देखने के लिए: मान लीजिए <math>x_i = \mu + \xi_i, \; \xi_i \sim N(0, \sigma_x^2)</math>, और <math>d_i = d + \epsilon_i, \; \epsilon_i \sim N(0, \sigma_d^2).</math> तब | ||
:<math>\hat{x}_i = \left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2} \mu + \frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2} d \right)+ \left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2} \xi_i + \frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2} \epsilon_i \right) </math>. | :<math>\hat{x}_i = \left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2} \mu + \frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2} d \right)+ \left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2} \xi_i + \frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2} \epsilon_i \right) </math>. | ||
पहला योग पश्च माध्य है, और दूसरे योग में, स्वतंत्रता को ध्यान में रखते हुए, | पहला योग पश्च माध्य है, और दूसरे योग में, स्वतंत्रता को ध्यान में रखते हुए, भिन्नता है | ||
:<math>\left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2}\right)^2 \sigma_x^2 + \left(\frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2}\right)^2 \sigma_d^2 = \frac{1}{1/\sigma_x^2 + 1/\sigma_d^2}</math>, | :<math>\left(\frac{1/\sigma_x^2}{1/\sigma_x^2 + 1/\sigma_d^2}\right)^2 \sigma_x^2 + \left(\frac{1/\sigma_d^2}{1/\sigma_x^2 + 1/\sigma_d^2}\right)^2 \sigma_d^2 = \frac{1}{1/\sigma_x^2 + 1/\sigma_d^2}</math>, | ||
जो पश्च | जो पश्च भिन्नता है। | ||
इस प्रकार ईएनकेएफ अब कलमैन गेन आव्यूह <math>K</math> में स्टेट सहसंयोजकता <math>Q</math> को संयोजन सदस्यों से गणना किए गए प्रारूप सहसंयोजकता <math>C</math> (जिसे समुच्चय सहसंयोजकता कहा जाता है) द्वारा प्रतिस्थापित करके प्राप्त किया जाता है, <ref name="Johns-2005-CEK">{{cite journal |first1=C. J. |last1=Johns |first2=J. |last2=Mandel |s2cid=14820232 |title=सहज डेटा सम्मिश्रण के लिए दो चरणों वाला कलमैन फ़िल्टर|journal=Environmental and Ecological Statistics |year=2008 |volume=15 |issue=1 |pages=101–110 |doi=10.1007/s10651-007-0033-0 |citeseerx=10.1.1.67.4916 }}</ref> जो कि <math>K=CH^{\mathrm{T}}\left( HCH^{\mathrm{T}}+R\right) ^{-1}</math> है | |||
==कार्यान्वयन == | |||
===मूल सूत्रीकरण=== | |||
यहां हम अनुसरण करते हैं।<ref name="Burgers-1998-ASE">{{cite journal |first1=G. |last1=Burgers |first2=P. J. |last2=van Leeuwen |first3=G. |last3=Evensen |title=एन्सेम्बल कलमन फ़िल्टर में विश्लेषण योजना|journal=Monthly Weather Review |volume=126 |issue= 6|year=1998 |pages=1719–1724 |doi=10.1175/1520-0493(1998)126<1719:ASITEK>2.0.CO;2 |citeseerx=10.1.1.41.5827 |bibcode=1998MWRv..126.1719B }}</ref><ref name="Evensen-2003-EKF">{{cite journal |first=G. |last=Evensen |s2cid=129233333 |title=The Ensemble Kalman Filter: Theoretical Formulation and Practical Implementation |journal=Ocean Dynamics |volume=53 |year=2003 |issue=4 |pages=343–367 |doi=10.1007/s10236-003-0036-9 |citeseerx=10.1.1.5.6990 |bibcode=2003OcDyn..53..343E }}</ref> मान लीजिए कि संयोजन आव्यूह <math>X</math> और डेटा आव्यूह <math>D</math> उपरोक्तानुसार हैं समुच्चय माध्य और सहसंयोजकता हैं | |||
यहां हम अनुसरण करते हैं।<ref name="Burgers-1998-ASE">{{cite journal |first1=G. |last1=Burgers |first2=P. J. |last2=van Leeuwen |first3=G. |last3=Evensen |title=एन्सेम्बल कलमन फ़िल्टर में विश्लेषण योजना|journal=Monthly Weather Review |volume=126 |issue= 6|year=1998 |pages=1719–1724 |doi=10.1175/1520-0493(1998)126<1719:ASITEK>2.0.CO;2 |citeseerx=10.1.1.41.5827 |bibcode=1998MWRv..126.1719B }}</ref><ref name="Evensen-2003-EKF">{{cite journal |first=G. |last=Evensen |s2cid=129233333 |title=The Ensemble Kalman Filter: Theoretical Formulation and Practical Implementation |journal=Ocean Dynamics |volume=53 |year=2003 |issue=4 |pages=343–367 |doi=10.1007/s10236-003-0036-9 |citeseerx=10.1.1.5.6990 |bibcode=2003OcDyn..53..343E }}</ref> मान लीजिए कि संयोजन | |||
:<math> E\left( X\right) =\frac{1}{N}\sum_{k=1}^{N}\mathbf{x}_{k},\quad C=\frac{AA^{T}}{N-1}, </math> | :<math> E\left( X\right) =\frac{1}{N}\sum_{k=1}^{N}\mathbf{x}_{k},\quad C=\frac{AA^{T}}{N-1}, </math> | ||
जहाँ | |||
:<math> A=X-E\left( X\right) \mathbf{e}_{1\times N} =X-\frac{1}{N}\left( X\mathbf{e}_{N\times1}\right) \mathbf{e}_{1\times N}, </math> | :<math> A=X-E\left( X\right) \mathbf{e}_{1\times N} =X-\frac{1}{N}\left( X\mathbf{e}_{N\times1}\right) \mathbf{e}_{1\times N}, </math> | ||
और <math>\mathbf{e}</math> संकेतित आकार के सभी के | और <math>\mathbf{e}</math> संकेतित आकार के सभी के आव्यूह को दर्शाता है। | ||
इसके पश्चात पिछला एन्सेम्बल <math>X^{p}</math> द्वारा दिया जाता है | |||
:<math> X^{p}=X+CH^{T}\left( HCH^{T}+R\right) ^{-1}(D-HX), </math> | :<math> X^{p}=X+CH^{T}\left( HCH^{T}+R\right) ^{-1}(D-HX), </math> | ||
जहां | जहां विकृत डेटा आव्यूह <math>D</math> ऊपर जैसा है। | ||
ध्यान दें कि तब से <math>R</math> | ध्यान दें कि तब से <math>R</math> सहसंयोजकता आव्यूह है, यह सदैव [[सकारात्मक अर्धनिश्चित मैट्रिक्स|धनात्मक अर्धनिश्चित]] आव्यूह होता है और सामान्यतः धनात्मक अर्धनिश्चित आव्यूह होता है, इसलिए उपरोक्त व्युत्क्रम उपस्थित है और सूत्र कोलेस्की अपघटन द्वारा कार्यान्वित किया जा सकता है।<ref name="Mandel-2006-EIE">{{cite journal |first=J. |last=Mandel |title=एन्सेम्बल कलमन फ़िल्टर का कुशल कार्यान्वयन|journal=Center for Computational Mathematics Reports |volume=231 |publisher=University of Colorado at Denver and Health Sciences Center |url=http://www.math.ucdenver.edu/ccm/reports/rep231.pdf |date=June 2006 }}</ref><ref name="Burgers-1998-ASE"/><ref name="Evensen-2003-EKF"/> <math>R</math> प्रारूपिक सहसंयोजकता <math>\tilde{D} \tilde{D}^{T}/\left( N-1\right) </math> द्वारा प्रतिस्थापित किया जाता है जहाँ <math>\tilde{D} = D - \frac{1}{N} d \, \mathbf{e}_{1\times N}</math> और व्युत्क्रम को छद्म व्युत्क्रम से परिवर्तित कर दिया जाता है, जिसकी गणना विलक्षण मान अपघटन (एसवीडी) का उपयोग करके की जाती है। | ||
चूँकि | चूँकि यह सूत्र प्रमुख ब्लास स्तर 3 ऑपरेशन वाले आव्यूह ऑपरेशन हैं,<ref name="Golub-1989-MAC">{{cite book |first1=G. H. |last1=Golub |author-link=Gene H. Golub |first2=C. F. V. |last2=Loan |title=मैट्रिक्स संगणना|publisher=Johns Hopkins Univ. Press |location=Baltimore |year=1989 |edition=Second |isbn=978-0-8018-3772-2 }}</ref> वह [[LAPACK|लैपैक]] (सीरियल और [[साझा स्मृति वास्तुकला|शेयर्ड मेमोरी]] कंप्यूटर पर) और [[ScaLAPACK|स्कालैपैक]] (वितरित मेमोरी कंप्यूटर पर) जैसे सॉफ़्टवेयर पैकेजों का उपयोग करके कुशल कार्यान्वयन के लिए उपयुक्त हैं।<ref name="Mandel-2006-EIE"/> किसी आव्यूह के व्युत्क्रम आव्यूह की गणना करने और उससे गुणा करने के अतिरिक्त [[उलटा मैट्रिक्स|विपरीत]] आव्यूह के चॉलेस्की अपघटन की गणना करना और व्युत्क्रम द्वारा गुणन को विभिन्न के साथ रैखिक प्रणाली के समाधान के रूप में मानना बहुत उत्तम (विभिन्न गुना सस्ता और अधिक स्पष्ट) है <ref name="Golub-1989-MAC"/> | ||
===अवलोकन आव्यूह-मुक्त कार्यान्वयन=== | |||
चूँकि हमने सहसंयोजकता आव्यूह को असेंबल सहसंयोजकता से परिवर्तित कर दिया है, इससे एक सरल सूत्र तैयार होता है जहाँ आव्यूह <math>H</math> को स्पष्ट रूप से निर्दिष्ट किए बिना सीधे सम्मिलन अवलोकनों का उपयोग किया जाता है। अधिक विशेष रूप से, फॉर्म के एक फलन <math>h(\mathbf{x})</math> को परिभाषित करें | |||
चूँकि हमने | |||
:<math> h(\mathbf{x})=H\mathbf{x}. </math> | :<math> h(\mathbf{x})=H\mathbf{x}. </math> | ||
इस प्रकार प्रोग्राम <math>h</math> इसे अवलोकन फलन कहा जाता है या, व्युत्क्रम समस्याओं के संदर्भ में, व्युत्क्रम समस्या व्युत्क्रम समस्याओं का संभाव्य सूत्रीकरण कहा जाता है। <math>h(\mathbf{x})</math> का मान है स्टेट के लिए डेटा <math>\mathbf{x}</math> का मूल्य क्या होगा यह मानते हुए कि माप स्पष्ट है। फिर पश्च संयोजन को पुनः लिखा जा सकता है | |||
:<math> X^{p}=X+\frac{1}{N-1}A\left( HA\right) ^{T}P^{-1}(D-HX) </math> | :<math> X^{p}=X+\frac{1}{N-1}A\left( HA\right) ^{T}P^{-1}(D-HX) </math> | ||
जहाँ | |||
:<math> HA=HX-\frac{1}{N}\left( \left( HX\right) \mathbf{e}_{N\times1}\right) \mathbf{e}_{1\times N}, </math> | :<math> HA=HX-\frac{1}{N}\left( \left( HX\right) \mathbf{e}_{N\times1}\right) \mathbf{e}_{1\times N}, </math> | ||
Line 101: | Line 95: | ||
:<math>\left[ HA\right] _{i} =H\mathbf{x}_{i}-H\frac{1}{N}\sum_{j=1}^{N}\mathbf{x}_{j}\ =h\left( \mathbf{x}_{i}\right) -\frac{1}{N}\sum_{j=1}^{N}h\left( \mathbf{x}_{j}\right) . </math> | :<math>\left[ HA\right] _{i} =H\mathbf{x}_{i}-H\frac{1}{N}\sum_{j=1}^{N}\mathbf{x}_{j}\ =h\left( \mathbf{x}_{i}\right) -\frac{1}{N}\sum_{j=1}^{N}h\left( \mathbf{x}_{j}\right) . </math> | ||
परिणाम स्वरुप, प्रत्येक संयोजन सदस्य पर एक बार अवलोकन फलन <math>H</math> का मूल्यांकन करके संयोजन अद्यतन की गणना की जा सकती है और आव्यूह <math>H</math> को स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह सूत्र एक निश्चित ऑफसेट <math>\mathbf{f}</math> के साथ एक अवलोकन फलन <math>h(\mathbf{x})=H\mathbf{x+f}</math> के लिए भी <ref name="Mandel-2006-EIE" /> प्रयुक्त करता है, जिसे जानने की भी आवश्यकता नहीं है स्पष्ट रूप से. उपरोक्त सूत्र का उपयोग सामान्यतः एक गैर-रेखीय अवलोकन फलन <math>H</math> के लिए किया गया है, जैसे कि तूफान भंवर की स्थिति <ref name="Chen-2006-AVP">{{cite journal |first1=Y. |last1=Chen |first2=C. |last2=Snyder |title=एन्सेम्बल कलमन फ़िल्टर के साथ भंवर स्थिति को आत्मसात करना|journal=Monthly Weather Review |volume=135 |year=2007 |issue=5 |pages=1828–1845 |doi=10.1175/MWR3351.1 |bibcode=2007MWRv..135.1828C }}</ref> उस स्थिति में, अवलोकन फलन को अनिवार्य रूप से समूह सदस्यों पर इसके मूल्यों से एक रैखिक फलन द्वारा अनुमानित किया जाता है। | |||
===बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन=== | ===बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन=== | ||
बड़ी संख्या के लिए | बड़ी संख्या में डेटा बिंदुओं के लिए, <math>P^{-1}</math> से गुणा एक बाधा बन जाता है। निम्न वैकल्पिक सूत्र तब लाभप्रद होता है जब डेटा बिंदुओं की संख्या <math>m</math> बड़ी होती है (जैसे कि ग्रिड या पिक्सेल डेटा को एसिमिलेसन करते समय) और डेटा त्रुटि सहसंयोजकता आव्यूह <math>R</math> विकर्ण होता है (जो कि तब होता है जब डेटा त्रुटियाँ असंबद्ध होती हैं), या सस्ता होता है विघटित (जैसे कि सीमित सहसंयोजकता दूरी के कारण बैंडेड) शर्मन-मॉरिसन-वुडबरी फॉर्मूला का उपयोग करता है <ref name="Hager-1989-UIM">{{cite journal |first=W. W. |last=Hager |title=मैट्रिक्स के व्युत्क्रम को अद्यतन करना|journal=[[SIAM Review]] |volume=31 |issue=2 |year=1989 |pages=221–239 |doi=10.1137/1031049 }}</ref> | ||
:<math> (R+UV^{T})^{-1}=R^{-1}-R^{-1}U(I+V^{T}R^{-1}U)^{-1}V^{T}R^{-1}, </math> | :<math> (R+UV^{T})^{-1}=R^{-1}-R^{-1}U(I+V^{T}R^{-1}U)^{-1}V^{T}R^{-1}, </math> | ||
साथ | साथ | ||
Line 114: | Line 110: | ||
:<math>\begin{align} P^{-1} & =\left( R+\frac{1}{N-1}HA\left( HA\right) ^{T}\right) ^{-1}\ = \\ | :<math>\begin{align} P^{-1} & =\left( R+\frac{1}{N-1}HA\left( HA\right) ^{T}\right) ^{-1}\ = \\ | ||
& =R^{-1}\left[ I-\frac{1}{N-1}\left( HA\right) \left( I+\left( HA\right) ^{T}R^{-1}\frac{1}{N-1}\left( HA\right) \right) ^{-1}\left( HA\right) ^{T}R^{-1}\right] , \end{align}</math> | & =R^{-1}\left[ I-\frac{1}{N-1}\left( HA\right) \left( I+\left( HA\right) ^{T}R^{-1}\frac{1}{N-1}\left( HA\right) \right) ^{-1}\left( HA\right) ^{T}R^{-1}\right] , \end{align}</math> | ||
= | जिसके लिए केवल आव्यूह <math>R</math> (सस्ता माना जाता है) और <math>m</math> दाहिनी ओर के आकार <math>N</math> के प्रणाली के समाधान की आवश्यकता होती है। ऑपरेशन गणना के लिए <ref name="Mandel-2006-EIE" /> देखें। | ||
==इसके अतिरिक्त विस्तार== | |||
इस प्रकार यहां वर्णित ईएनकेएफ संस्करण में डेटा का यादृच्छिककरण सम्मिलित है। डेटा के यादृच्छिकीकरण के बिना फ़िल्टर के लिए, देखें।<ref name="Anderson-2001-EAK">{{cite journal |first=J. L. |last=Anderson |title=डेटा सम्मिश्रण के लिए एक संयोजन समायोजन कलमन फ़िल्टर|journal=Monthly Weather Review |volume=129 |year=2001 |issue=12 |pages=2884–2903 |doi=10.1175/1520-0493(2001)129<2884:AEAKFF>2.0.CO;2 |citeseerx=10.1.1.5.9952 |bibcode=2001MWRv..129.2884A }}</ref><ref name="Evensen-2004-SSR">{{cite journal |first=G. |last=Evensen |s2cid=120171951 |title=EnKF के लिए नमूनाकरण रणनीतियाँ और वर्गमूल विश्लेषण योजनाएँ|journal=Ocean Dynamics |volume=54 |year=2004 |issue=6 |pages=539–560 |doi=10.1007/s10236-004-0099-2 |citeseerx=10.1.1.3.6213 |bibcode=2004OcDyn..54..539E }}</ref><ref name="Tippett-2003-ESR">{{cite journal |first1=M. K. |last1=Tippett |first2=J. L. |last2=Anderson |first3=C. H. |last3=Bishop |first4=T. M. |last4=Hamill |first5=J. S. |last5=Whitaker |title=वर्गमूल फिल्टरों को इकट्ठा करें|journal=Monthly Weather Review |volume=131 |year=2003 |issue=7 |pages=1485–1490 |doi=10.1175/1520-0493(2003)131<1485:ESRF>2.0.CO;2 |citeseerx=10.1.1.332.775 |bibcode=2003MWRv..131.1485T }}</ref> चूंकि समूह सहसंयोजकता में [[रैंक की कमी|रैंक डेफीसिएंट]] है (समूह सदस्यों की तुलना में विभिन्न अधिक स्टेट वैरिएबल हैं, सामान्यतः लाखों, सामान्यतः सौ से कम), इसमें उन बिंदुओं के जोड़े के लिए बड़े शब्द हैं जो स्थानिक रूप से दूर हैं। चूँकि वास्तव में दूर के स्थानों पर भौतिक क्षेत्रों के मान इतने [[सहसंबद्ध]] नहीं हैं, दूरी के आधार पर सहसंयोजकता आव्यूह को कृत्रिम रूप से कम कर दिया जाता है, इस प्रकार जो स्थानीयकृत कलमैन फ़िल्टर एल्गोरिदम को जन्म देता है।<ref name="Anderson-2003-LLS">{{cite journal |first=J. L. |last=Anderson |title=संयोजन फ़िल्टरिंग के लिए एक स्थानीय न्यूनतम वर्ग ढाँचा|journal=Monthly Weather Review |volume=131 |year=2003 |issue=4 |pages=634–642 |doi=10.1175/1520-0493(2003)131<0634:ALLSFF>2.0.CO;2 |citeseerx=10.1.1.10.6543 |bibcode=2003MWRv..131..634A }}</ref><ref name="Ott-2003-LEK">{{cite journal |author-link1=Edward Ott |first1=E. |last1=Ott |first2=B. R. |last2=Hunt |first3=I. |last3=Szunyogh |first4=A. V. |last4=Zimin |first5=E. J. |last5=Kostelich |first6=M. |last6=Corazza |author-link7=Eugenia Kalnay |first7=E. |last7=Kalnay |first8=D. |last8=Patil |author-link9=James A. Yorke |first9=J. A. |last9=Yorke |s2cid=218577557 |title=वायुमंडलीय डेटा सम्मिलन के लिए एक स्थानीय पहनावा कलमन फ़िल्टर|journal=[[Tellus A]] |volume=56 |year=2004 |issue=5 |pages=415–428 |doi=10.3402/tellusa.v56i5.14462 |arxiv=physics/0203058 |bibcode=2004TellA..56..415O }}</ref> ये विधियां गणनाओं में उपयोग किए जाने वाले सहसंयोजकता आव्यूह को संशोधित करती हैं और परिणामस्वरूप, पोस्टीरियर एन्सेम्बल अब केवल पोस्टीरियर संयोजन के रैखिक संयोजनों से नहीं बना है। | |||
इस प्रकार गैर-रेखीय समस्याओं के लिए, ईएनकेएफ गैर-भौतिक अवस्थाओं के साथ पश्च संयोजन बना सकता है। इसे [[तिखोनोव नियमितीकरण|नियमितीकरण]] द्वारा कम किया जा सकता है, जैसे कि बड़े स्थानिक [[ ग्रेडियेंट |ग्रेडियेंट]] वाले स्टेट की दंड विधि है।<ref name="Johns-2005-CEK"/> | |||
[[सुसंगत विशेषता|कोहेरेंट विशेषता]]ओं वाली समस्याओं, जैसे कि हरिकेन, थंडरस्टॉर्म, [[ आग रेखा |फायरलाइन]] , [[ तूफान रेखा |स्क्वॉल लाइन]] और रेन फ्रंट के लिए, अंतरिक्ष (इसके ग्रिड) में स्टेट को विकृत करके और साथ ही स्टेट के आयामों को सही करके संख्यात्मक मॉडल स्थिति को समायोजित करने की आवश्यकता है। इस प्रकार 2007 में, रवेला एट अल एनसेंबल का उपयोग करके संयुक्त स्थिति-आयाम समायोजन मॉडल प्रस्तुत करें, और व्यवस्थित रूप से अनुक्रमिक सन्निकटन प्राप्त करें जिसे एनकेएफ और अन्य फॉर्मूलेशन दोनों पर प्रयुक्त किया जा सकता है।<ref name ="Ravela-2007">{{cite journal |first1=S. |last1=Ravela |first2=K. |last2=Emanuel |author-link2=Kerry Emanuel |first3=D. |last3=McLaughlin |title=फ़ील्ड संरेखण द्वारा डेटा सम्मिलन|journal=[[Physica (journal)|Physica]] |series=D: Nonlinear Phenomena |volume=230 |issue=1–2 |year=2007 |pages=127–145 |doi=10.1016/j.physd.2006.09.035 |bibcode=2007PhyD..230..127R }}</ref> इस प्रकार उनकी पद्धति यह धारणा नहीं बनाती है कि आयाम और स्थिति त्रुटियां स्वतंत्र या संयुक्त रूप से गॉसियन हैं, जैसा कि अन्य करते हैं। एनकेएफ स्टेट के रैखिक संयोजनों के अतिरिक्त, इमेज रजिस्ट्रेसन और मॉर्फिंग से उधार ली गई तकनीकों द्वारा प्राप्त मध्यवर्ती स्टेट को नियोजित करता है।<ref name="Beezley-2007-MEK">{{cite journal |first1=J. D. |last1=Beezley |first2=J. |last2=Mandel |s2cid=1009227 |title=मॉर्फ़िंग पहनावा कलमन फ़िल्टर|journal=Tellus A |year=2008 |volume=60 |issue=1 |pages=131–140 |doi=10.1111/j.1600-0870.2007.00275.x |arxiv=0705.3693 |bibcode=2008TellA..60..131B }}</ref><ref name="Mandel-2006-PME">{{cite conference |first1=J. |last1=Mandel |first2=J. D. |last2=Beezley |title=उच्च आयामी नॉनलाइनर सिस्टम में विरल डेटा को आत्मसात करने के लिए प्रीडिक्टर-करेक्टर और मॉर्फिंग एसेंबल फिल्टर|work=CCM Report 239 |publisher=University of Colorado at Denver and Health Sciences Center |url=http://www.math.ucdenver.edu/ccm/reports/rep239.pdf |date=November 2006 |conference=11th Symposium on Integrated Observing and Assimilation Systems for the Atmosphere, Oceans, and Land Surface (IOAS-AOLS), CD-ROM, Paper 4.12, 87th American Meteorological Society Annual Meeting, San Antonio, TX, January 2007 }}</ref> | |||
औपचारिक रूप से, ईएनकेएफ गॉसियन धारणा पर विश्वास करते हैं। इस प्रकार व्यवहार में इनका उपयोग अरैखिक समस्याओं के लिए भी किया जा सकता है, जहां गॉसियन धारणा संतुष्ट नहीं हो सकती है। इस प्रकार एनकेएफ में गॉसियन धारणा को शिथिल करने का प्रयास करने वाले संबंधित फिल्टर, इसके लाभ को संरक्षित करते हुए, इसमें ऐसे फिल्टर सम्मिलित हैं जो विभिन्न गॉसियन कर्नेल के साथ स्टेट पीडीएफ को फिट करते हैं,<ref name="Anderson-1999-MCI">{{cite journal |first1=J. L. |last1=Anderson |first2=S. L. |last2=Anderson |title=सामूहिक सम्मिलन और पूर्वानुमान उत्पन्न करने के लिए नॉनलाइनियर फ़िल्टरिंग समस्या का मोंटे कार्लो कार्यान्वयन|journal=Monthly Weather Review |volume=127 |issue=12 |year=1999 |pages=2741–2758 |doi=10.1175/1520-0493(1999)127<2741:AMCIOT>2.0.CO;2 |bibcode=1999MWRv..127.2741A }}</ref> फ़िल्टर जो [[गाऊसी मिश्रण]] द्वारा स्टेट पीडीएफ का अनुमान लगाते हैं,<ref name="Bengtsson-2003-NFE">{{cite journal |first1=T. |last1=Bengtsson |first2=C. |last2=Snyder |first3=D. |last3=Nychka |title=उच्च आयामी प्रणालियों के लिए एक अरेखीय संयोजन फ़िल्टर की ओर|journal=Journal of Geophysical Research: Atmospheres |volume=108 |issue=D24 |year=2003 |pages=STS 2–1–10 |doi=10.1029/2002JD002900 |bibcode=2003JGRD..108.8775B |doi-access=free }}</ref> इस प्रकार [[घनत्व अनुमान]] द्वारा पार्टिकल भार की गणना के साथ पार्टिकल फिल्टर का प्रकार,<ref name="Mandel-2006-PME" /> और पार्टिकल फ़िल्टर अनुक्रमिक महत्व पुनः प्रारूपिककरण (एसआईआर) को कम करने के लिए [[कॉची वितरण]] डेटा पीडीएफ के साथ पार्टिकल फ़िल्टर का प्रकार है।<ref name="vanLeeuwen-2003-VMF">{{cite journal |first=P. |last=van Leeuwen |title=बड़े पैमाने के अनुप्रयोगों के लिए एक विचरण-न्यूनीकरण फ़िल्टर|journal=Monthly Weather Review |volume=131 |year=2003 |issue=9 |pages=2071–2084 |doi=10.1175/1520-0493(2003)131<2071:AVFFLA>2.0.CO;2 |citeseerx=10.1.1.7.3719 |bibcode=2003MWRv..131.2071V }}</ref> | |||
==यह भी देखें== | ==यह भी देखें== | ||
* डेटा | * डेटा एसीमिलेसन | ||
* | *न्यूमेरिकल वेदर प्रेडिक्शन एसेम्बल | ||
* | * पार्टिकल फिल्टर | ||
* [[पुनरावर्ती बायेसियन अनुमान]] | * [[पुनरावर्ती बायेसियन अनुमान|रिकर्सिव बायेसियन एसटीमेसन]] | ||
==संदर्भ== | ==संदर्भ== | ||
{{reflist|2}} | {{reflist|2}} | ||
==बाहरी संबंध== | ==बाहरी संबंध== | ||
* [http://enkf.nersc.no | * [http://enkf.nersc.no ईएनकेएफ webpage] | ||
* [http://topaz.nersc.no TOPAZ, real-time forecasting of the North Atlantic ocean and Arctic sea-ice with the | * [http://topaz.nersc.no TOPAZ, real-time forecasting of the North Atlantic ocean and Arctic sea-ice with the ईएनकेएफ] | ||
* [https://github.com/sakov/enkf-c | * [https://github.com/sakov/enkf-c ईएनकेएफ-C, a compact framework for data assimilation into large-स्काle layered geophysical models with the ईएनकेएफ] | ||
* [http://pdaf.awi.de PDAF] – [[Parallel Data Assimilation Framework]] – an open-source software for data assimilation providing different variants of the | * [http://pdaf.awi.de PDAF] – [[Parallel Data Assimilation Framework]] – an open-source software for data assimilation providing different variants of the ईएनकेएफ | ||
{{DEFAULTSORT:Ensemble Kalman Filter}}[[Category: रैखिक फिल्टर]] [[Category: अरेखीय फिल्टर]] [[Category: बायेसियन आँकड़े]] [[Category: संकेत अनुमान]] [[Category: मोंटे कार्लो विधियाँ]] | {{DEFAULTSORT:Ensemble Kalman Filter}}[[Category: रैखिक फिल्टर]] [[Category: अरेखीय फिल्टर]] [[Category: बायेसियन आँकड़े]] [[Category: संकेत अनुमान]] [[Category: मोंटे कार्लो विधियाँ]] | ||
Line 149: | Line 142: | ||
[[Category: Machine Translated Page]] | [[Category: Machine Translated Page]] | ||
[[Category:Created On 15/08/2023]] | [[Category:Created On 15/08/2023]] | ||
[[Category:Vigyan Ready]] |
Latest revision as of 09:33, 1 December 2023
एन्सेम्बल कलमैन फ़िल्टर (ईएनकेएफ) रिकर्सिव फ़िल्टर है जो बड़ी संख्या में वैरिएबल वाली समस्याओं के लिए उपयुक्त है, जैसे कि भूभौतिकीय मॉडल में आंशिक अंतर समीकरणों का विवेकीकरण ईएनकेएफ की उत्पत्ति बड़ी समस्याओं के लिए कलमन फ़िल्टर के संस्करण के रूप में हुई (अनिवार्य रूप से, सहसंयोजकता आव्यूह को प्रारूपिक सहसंयोजकता आव्यूह द्वारा प्रतिस्थापित किया जाता है), और यह अब संयोजन पूर्वानुमान का महत्वपूर्ण डेटा एसीमिलेसन घटक है। ईएनकेएफ पार्टिकल फिल्टर से संबंधित है (इस संदर्भ में, पार्टिकल समूह सदस्य के समान है) किन्तु ईएनकेएफ यह धारणा बनाता है कि इसमें सम्मिलित सभी संभाव्यता वितरण गाऊसी हैं; जब यह प्रयुक्त होता है, जिससे यह पार्टिकल फिल्टर की तुलना में बहुत अधिक उत्तम होता है।
परिचय
एन्सेम्बल कलमैन फ़िल्टर (ईएनकेएफ) बायेसियन अपडेट समस्या का मोंटे कार्लो विधि कार्यान्वयन है: मॉडल किए गए प्रणाली की स्थिति (प्रियोर, जिसे अधिकांशतः भूविज्ञान में पूर्वानुमान कहा जाता है) और डेटा संभावना की संभाव्यता घनत्व फलन (पीडीएफ) दी गई है। बेयस प्रमेय का उपयोग डेटा संभावना को ध्यान में रखने के पश्चात पीडीएफ प्राप्त करने के लिए किया जाता है (पोस्टीरियर संभावना, जिसे अधिकांशतः विश्लेषण कहा जाता है)। इसे बायेसियन अपडेट कहा जाता है। बायेसियन अपडेट को समय-समय पर नए डेटा को सम्मिलित करते हुए मॉडल को आगे बढ़ाने के साथ जोड़ा जाता है। मूल कलमैन फ़िल्टर, 1960 में प्रस्तुत किया गया था,[1] मानता है कि सभी पीडीएफ सामान्य वितरण (गॉसियन धारणा) हैं और बायेसियन अपडेट द्वारा माध्य और सहसंयोजकता आव्यूह के परिवर्तन के लिए बीजगणितीय सूत्र प्रदान करता है, साथ ही समय में माध्य और सहसंयोजकता को आगे बढ़ाने के लिए सूत्र प्रदान करता है, किन्तु प्रणाली रैखिक हो। चूंकि, उच्च-आयामी प्रणालियों के लिए सहसंयोजकता आव्यूह को बनाए रखना कम्प्यूटेशनल रूप से संभव नहीं है। इस कारण से, ईएनकेएफएस विकसित किए गए थे।[2][3] ईएनकेएफएस स्टेट वैक्टरों के संग्रह का उपयोग करके प्रणाली स्थिति के वितरण का प्रतिनिधित्व करते हैं, जिसे संख्यात्मक मौसम पूर्वानुमान या एनसेम्बल कहा जाता है, और सहसंयोजकता आव्यूह को समुच्चय से गणना किए गए प्रारूपिक सहसंयोजकता द्वारा प्रतिस्थापित करते हैं। इस प्रकार समूह को ऐसे संचालित किया जाता है जैसे कि यह यादृच्छिक प्रारूपिक हो, किन्तु समूह के सदस्य वास्तव में सांख्यिकीय स्वतंत्रता नहीं रखते हैं, क्योंकि वह सभी ईएनकेएफ संगठित करते हैं। ईएनकेएफएस का लाभ यह है कि पीडीएफ को समय पर आगे बढ़ाने का कार्य केवल समूह के प्रत्येक सदस्य को आगे बढ़ाना है।[4]
व्युत्पत्ति
कलमैन फ़िल्टर
मान लीजिए कि एक मॉडल के -आयामी स्टेट वेक्टर को दर्शाता है, और मान लेता है कि इसमें माध्य और सहसंयोजकता के साथ गॉसियन संभाव्यता वितरण है, अर्थात, इसका पीडीएफ है
यहां और नीचे, का अर्थ आनुपातिक है; एक पीडीएफ को सदैव स्केल किया जाता है जिससे पूर्ण स्थान पर इसका अभिन्न अंग एक होता है। यह , जिसे पूर्व कहा जाता है, मॉडल को चलाकर समय पर विकसित किया गया था और अब इसे नए डेटा के लिए अपडेट किया जाना है। यह मान लेना स्वाभाविक है कि डेटा का त्रुटि वितरण ज्ञात है; डेटा को त्रुटि अनुमान के साथ आना होगा, अन्यथा वह अर्थहीन हैं। यहां, डेटा को सहसंयोजकता R और माध्य के साथ गॉसियन पीडीएफ माना जाता है, जहां तथाकथित अवलोकन आव्यूह है। सहसंयोजकता आव्यूह आर डेटा की त्रुटि के अनुमान का वर्णन करता है; यदि डेटा वेक्टर की प्रविष्टियों में यादृच्छिक त्रुटियां स्वतंत्र हैं, तो विकर्ण है और इसकी विकर्ण प्रविष्टियां डेटा वेक्टर की संबंधित प्रविष्टियों की त्रुटि के मानक विचलन ("त्रुटि आकार") के वर्ग हैं मान वह है जो डेटा त्रुटियों के अभाव में स्टेट के लिए डेटा का मान होगा। फिर प्रणाली स्थिति की नियमबद्ध डेटा की संभाव्यता घनत्व , जिसे डेटा संभावना कहा जाता है, है
इस प्रकार स्थिति की पीडीएफ और डेटा संभावना को बेयस प्रमेय द्वारा डेटा के मूल्य पर नियमबद्ध प्रणाली स्थिति की नई संभाव्यता घनत्व देने के लिए संयोजित किया गया है,
इस प्रकार आंकड़ा एक बार प्राप्त होने के पश्चात यह निश्चित हो जाता है, इसलिए पिछली स्थिति को इससे निरूपित करें के अतिरिक्त और पिछला पीडीएफ़ द्वारा इसे बीजगणितीय जोड़ द्वारा दिखाया जा सकता है [5] पिछला पीडीएफ भी गॉसियन है,
इस प्रकार कलमैन अपडेट फ़ार्मुलों द्वारा दिए गए पश्च माध्य @ और सहसंयोजकता @ के साथ किया जाता है
जहाँ
तथाकथित कलमन फिल्टर कलमैन लाभ व्युत्पत्ति आव्यूह है।
एन्सेम्बल कलमैन फ़िल्टर
ईएनकेएफ कलमैन फ़िल्टर का मोंटे कार्लो सन्निकटन है, जो स्टेट वेक्टर के पीडीएफ के सहसंयोजकता आव्यूह को विकसित करने से बचाता है . इसके अतिरिक्त, पीडीएफ को समूह द्वारा दर्शाया जाता है
इस प्रकार X एक आव्यूह है जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से एक प्रारूप बनाएंगे। चूँकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से स्वतंत्र नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ चरण उन्हें एक साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वह वास्तव में स्वतंत्र होंते है।
डेटा को आव्यूह में दोहराएं जाते है
जिससे प्रत्येक कॉलम में डेटा वेक्टर प्लस एम-आयामी सामान्य वितरण से एक यादृच्छिक वेक्टर सम्मिलित होते है। यदि इसके अतिरिक्त, के कॉलम पूर्व संभाव्यता वितरण से एक प्रारूप हैं,
इस प्रकार पश्च संभाव्यता वितरण से एक प्रारूप बनाएं। इसे के साथ अदिश स्थिति में देखने के लिए: मान लीजिए , और तब
- .
पहला योग पश्च माध्य है, और दूसरे योग में, स्वतंत्रता को ध्यान में रखते हुए, भिन्नता है
- ,
जो पश्च भिन्नता है।
इस प्रकार ईएनकेएफ अब कलमैन गेन आव्यूह में स्टेट सहसंयोजकता को संयोजन सदस्यों से गणना किए गए प्रारूप सहसंयोजकता (जिसे समुच्चय सहसंयोजकता कहा जाता है) द्वारा प्रतिस्थापित करके प्राप्त किया जाता है, [6] जो कि है
कार्यान्वयन
मूल सूत्रीकरण
यहां हम अनुसरण करते हैं।[7][8] मान लीजिए कि संयोजन आव्यूह और डेटा आव्यूह उपरोक्तानुसार हैं समुच्चय माध्य और सहसंयोजकता हैं
जहाँ
और संकेतित आकार के सभी के आव्यूह को दर्शाता है।
इसके पश्चात पिछला एन्सेम्बल द्वारा दिया जाता है
जहां विकृत डेटा आव्यूह ऊपर जैसा है।
ध्यान दें कि तब से सहसंयोजकता आव्यूह है, यह सदैव धनात्मक अर्धनिश्चित आव्यूह होता है और सामान्यतः धनात्मक अर्धनिश्चित आव्यूह होता है, इसलिए उपरोक्त व्युत्क्रम उपस्थित है और सूत्र कोलेस्की अपघटन द्वारा कार्यान्वित किया जा सकता है।[9][7][8] प्रारूपिक सहसंयोजकता द्वारा प्रतिस्थापित किया जाता है जहाँ और व्युत्क्रम को छद्म व्युत्क्रम से परिवर्तित कर दिया जाता है, जिसकी गणना विलक्षण मान अपघटन (एसवीडी) का उपयोग करके की जाती है।
चूँकि यह सूत्र प्रमुख ब्लास स्तर 3 ऑपरेशन वाले आव्यूह ऑपरेशन हैं,[10] वह लैपैक (सीरियल और शेयर्ड मेमोरी कंप्यूटर पर) और स्कालैपैक (वितरित मेमोरी कंप्यूटर पर) जैसे सॉफ़्टवेयर पैकेजों का उपयोग करके कुशल कार्यान्वयन के लिए उपयुक्त हैं।[9] किसी आव्यूह के व्युत्क्रम आव्यूह की गणना करने और उससे गुणा करने के अतिरिक्त विपरीत आव्यूह के चॉलेस्की अपघटन की गणना करना और व्युत्क्रम द्वारा गुणन को विभिन्न के साथ रैखिक प्रणाली के समाधान के रूप में मानना बहुत उत्तम (विभिन्न गुना सस्ता और अधिक स्पष्ट) है [10]
अवलोकन आव्यूह-मुक्त कार्यान्वयन
चूँकि हमने सहसंयोजकता आव्यूह को असेंबल सहसंयोजकता से परिवर्तित कर दिया है, इससे एक सरल सूत्र तैयार होता है जहाँ आव्यूह को स्पष्ट रूप से निर्दिष्ट किए बिना सीधे सम्मिलन अवलोकनों का उपयोग किया जाता है। अधिक विशेष रूप से, फॉर्म के एक फलन को परिभाषित करें
इस प्रकार प्रोग्राम इसे अवलोकन फलन कहा जाता है या, व्युत्क्रम समस्याओं के संदर्भ में, व्युत्क्रम समस्या व्युत्क्रम समस्याओं का संभाव्य सूत्रीकरण कहा जाता है। का मान है स्टेट के लिए डेटा का मूल्य क्या होगा यह मानते हुए कि माप स्पष्ट है। फिर पश्च संयोजन को पुनः लिखा जा सकता है
जहाँ
और
साथ
परिणाम स्वरुप, प्रत्येक संयोजन सदस्य पर एक बार अवलोकन फलन का मूल्यांकन करके संयोजन अद्यतन की गणना की जा सकती है और आव्यूह को स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह सूत्र एक निश्चित ऑफसेट के साथ एक अवलोकन फलन के लिए भी [9] प्रयुक्त करता है, जिसे जानने की भी आवश्यकता नहीं है स्पष्ट रूप से. उपरोक्त सूत्र का उपयोग सामान्यतः एक गैर-रेखीय अवलोकन फलन के लिए किया गया है, जैसे कि तूफान भंवर की स्थिति [11] उस स्थिति में, अवलोकन फलन को अनिवार्य रूप से समूह सदस्यों पर इसके मूल्यों से एक रैखिक फलन द्वारा अनुमानित किया जाता है।
बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन
बड़ी संख्या में डेटा बिंदुओं के लिए, से गुणा एक बाधा बन जाता है। निम्न वैकल्पिक सूत्र तब लाभप्रद होता है जब डेटा बिंदुओं की संख्या बड़ी होती है (जैसे कि ग्रिड या पिक्सेल डेटा को एसिमिलेसन करते समय) और डेटा त्रुटि सहसंयोजकता आव्यूह विकर्ण होता है (जो कि तब होता है जब डेटा त्रुटियाँ असंबद्ध होती हैं), या सस्ता होता है विघटित (जैसे कि सीमित सहसंयोजकता दूरी के कारण बैंडेड) शर्मन-मॉरिसन-वुडबरी फॉर्मूला का उपयोग करता है [12]
साथ
देता है
जिसके लिए केवल आव्यूह (सस्ता माना जाता है) और दाहिनी ओर के आकार के प्रणाली के समाधान की आवश्यकता होती है। ऑपरेशन गणना के लिए [9] देखें।
इसके अतिरिक्त विस्तार
इस प्रकार यहां वर्णित ईएनकेएफ संस्करण में डेटा का यादृच्छिककरण सम्मिलित है। डेटा के यादृच्छिकीकरण के बिना फ़िल्टर के लिए, देखें।[13][14][15] चूंकि समूह सहसंयोजकता में रैंक डेफीसिएंट है (समूह सदस्यों की तुलना में विभिन्न अधिक स्टेट वैरिएबल हैं, सामान्यतः लाखों, सामान्यतः सौ से कम), इसमें उन बिंदुओं के जोड़े के लिए बड़े शब्द हैं जो स्थानिक रूप से दूर हैं। चूँकि वास्तव में दूर के स्थानों पर भौतिक क्षेत्रों के मान इतने सहसंबद्ध नहीं हैं, दूरी के आधार पर सहसंयोजकता आव्यूह को कृत्रिम रूप से कम कर दिया जाता है, इस प्रकार जो स्थानीयकृत कलमैन फ़िल्टर एल्गोरिदम को जन्म देता है।[16][17] ये विधियां गणनाओं में उपयोग किए जाने वाले सहसंयोजकता आव्यूह को संशोधित करती हैं और परिणामस्वरूप, पोस्टीरियर एन्सेम्बल अब केवल पोस्टीरियर संयोजन के रैखिक संयोजनों से नहीं बना है।
इस प्रकार गैर-रेखीय समस्याओं के लिए, ईएनकेएफ गैर-भौतिक अवस्थाओं के साथ पश्च संयोजन बना सकता है। इसे नियमितीकरण द्वारा कम किया जा सकता है, जैसे कि बड़े स्थानिक ग्रेडियेंट वाले स्टेट की दंड विधि है।[6]
कोहेरेंट विशेषताओं वाली समस्याओं, जैसे कि हरिकेन, थंडरस्टॉर्म, फायरलाइन , स्क्वॉल लाइन और रेन फ्रंट के लिए, अंतरिक्ष (इसके ग्रिड) में स्टेट को विकृत करके और साथ ही स्टेट के आयामों को सही करके संख्यात्मक मॉडल स्थिति को समायोजित करने की आवश्यकता है। इस प्रकार 2007 में, रवेला एट अल एनसेंबल का उपयोग करके संयुक्त स्थिति-आयाम समायोजन मॉडल प्रस्तुत करें, और व्यवस्थित रूप से अनुक्रमिक सन्निकटन प्राप्त करें जिसे एनकेएफ और अन्य फॉर्मूलेशन दोनों पर प्रयुक्त किया जा सकता है।[18] इस प्रकार उनकी पद्धति यह धारणा नहीं बनाती है कि आयाम और स्थिति त्रुटियां स्वतंत्र या संयुक्त रूप से गॉसियन हैं, जैसा कि अन्य करते हैं। एनकेएफ स्टेट के रैखिक संयोजनों के अतिरिक्त, इमेज रजिस्ट्रेसन और मॉर्फिंग से उधार ली गई तकनीकों द्वारा प्राप्त मध्यवर्ती स्टेट को नियोजित करता है।[19][20]
औपचारिक रूप से, ईएनकेएफ गॉसियन धारणा पर विश्वास करते हैं। इस प्रकार व्यवहार में इनका उपयोग अरैखिक समस्याओं के लिए भी किया जा सकता है, जहां गॉसियन धारणा संतुष्ट नहीं हो सकती है। इस प्रकार एनकेएफ में गॉसियन धारणा को शिथिल करने का प्रयास करने वाले संबंधित फिल्टर, इसके लाभ को संरक्षित करते हुए, इसमें ऐसे फिल्टर सम्मिलित हैं जो विभिन्न गॉसियन कर्नेल के साथ स्टेट पीडीएफ को फिट करते हैं,[21] फ़िल्टर जो गाऊसी मिश्रण द्वारा स्टेट पीडीएफ का अनुमान लगाते हैं,[22] इस प्रकार घनत्व अनुमान द्वारा पार्टिकल भार की गणना के साथ पार्टिकल फिल्टर का प्रकार,[20] और पार्टिकल फ़िल्टर अनुक्रमिक महत्व पुनः प्रारूपिककरण (एसआईआर) को कम करने के लिए कॉची वितरण डेटा पीडीएफ के साथ पार्टिकल फ़िल्टर का प्रकार है।[23]
यह भी देखें
- डेटा एसीमिलेसन
- न्यूमेरिकल वेदर प्रेडिक्शन एसेम्बल
- पार्टिकल फिल्टर
- रिकर्सिव बायेसियन एसटीमेसन
संदर्भ
- ↑ Kalman, R. E. (1960). "रेखीय छानने और भविष्यवाणी की समस्याओं के लिए एक नया दृष्टिकोण". Journal of Basic Engineering. 82 (1): 35–45. doi:10.1115/1.3662552. S2CID 1242324.
- ↑ Evensen, G. (1994). "त्रुटि आँकड़ों का पूर्वानुमान लगाने के लिए मोंटे कार्लो विधियों का उपयोग करके गैर-रेखीय अर्ध-जियोस्ट्रोफिक मॉडल के साथ अनुक्रमिक डेटा आत्मसात". Journal of Geophysical Research. 99 (C5): 143–162. Bibcode:1994JGR....9910143E. doi:10.1029/94JC00572. hdl:1956/3035.
- ↑ Houtekamer, P.; Mitchell, H. L. (1998). "एन्सेम्बल कलमैन फ़िल्टर तकनीक का उपयोग करके डेटा आत्मसात करना". Monthly Weather Review. 126 (3): 796–811. Bibcode:1998MWRv..126..796H. CiteSeerX 10.1.1.3.1706. doi:10.1175/1520-0493(1998)126<0796:DAUAEK>2.0.CO;2.
- ↑ For a survey of EnKF and related data assimilation techniques, see Evensen, G. (2007). Data Assimilation : The Ensemble Kalman Filter. Berlin: Springer. ISBN 978-3-540-38300-0.
- ↑ Anderson, B. D. O.; Moore, J. B. (1979). इष्टतम फ़िल्टरिंग. Englewood Cliffs, NJ: Prentice-Hall. ISBN 978-0-13-638122-8.
- ↑ 6.0 6.1 Johns, C. J.; Mandel, J. (2008). "सहज डेटा सम्मिश्रण के लिए दो चरणों वाला कलमैन फ़िल्टर". Environmental and Ecological Statistics. 15 (1): 101–110. CiteSeerX 10.1.1.67.4916. doi:10.1007/s10651-007-0033-0. S2CID 14820232.
- ↑ 7.0 7.1 Burgers, G.; van Leeuwen, P. J.; Evensen, G. (1998). "एन्सेम्बल कलमन फ़िल्टर में विश्लेषण योजना". Monthly Weather Review. 126 (6): 1719–1724. Bibcode:1998MWRv..126.1719B. CiteSeerX 10.1.1.41.5827. doi:10.1175/1520-0493(1998)126<1719:ASITEK>2.0.CO;2.
- ↑ 8.0 8.1 Evensen, G. (2003). "The Ensemble Kalman Filter: Theoretical Formulation and Practical Implementation". Ocean Dynamics. 53 (4): 343–367. Bibcode:2003OcDyn..53..343E. CiteSeerX 10.1.1.5.6990. doi:10.1007/s10236-003-0036-9. S2CID 129233333.
- ↑ 9.0 9.1 9.2 9.3 Mandel, J. (June 2006). "एन्सेम्बल कलमन फ़िल्टर का कुशल कार्यान्वयन" (PDF). Center for Computational Mathematics Reports. University of Colorado at Denver and Health Sciences Center. 231.
- ↑ 10.0 10.1 Golub, G. H.; Loan, C. F. V. (1989). मैट्रिक्स संगणना (Second ed.). Baltimore: Johns Hopkins Univ. Press. ISBN 978-0-8018-3772-2.
- ↑ Chen, Y.; Snyder, C. (2007). "एन्सेम्बल कलमन फ़िल्टर के साथ भंवर स्थिति को आत्मसात करना". Monthly Weather Review. 135 (5): 1828–1845. Bibcode:2007MWRv..135.1828C. doi:10.1175/MWR3351.1.
- ↑ Hager, W. W. (1989). "मैट्रिक्स के व्युत्क्रम को अद्यतन करना". SIAM Review. 31 (2): 221–239. doi:10.1137/1031049.
- ↑ Anderson, J. L. (2001). "डेटा सम्मिश्रण के लिए एक संयोजन समायोजन कलमन फ़िल्टर". Monthly Weather Review. 129 (12): 2884–2903. Bibcode:2001MWRv..129.2884A. CiteSeerX 10.1.1.5.9952. doi:10.1175/1520-0493(2001)129<2884:AEAKFF>2.0.CO;2.
- ↑ Evensen, G. (2004). "EnKF के लिए नमूनाकरण रणनीतियाँ और वर्गमूल विश्लेषण योजनाएँ". Ocean Dynamics. 54 (6): 539–560. Bibcode:2004OcDyn..54..539E. CiteSeerX 10.1.1.3.6213. doi:10.1007/s10236-004-0099-2. S2CID 120171951.
- ↑ Tippett, M. K.; Anderson, J. L.; Bishop, C. H.; Hamill, T. M.; Whitaker, J. S. (2003). "वर्गमूल फिल्टरों को इकट्ठा करें". Monthly Weather Review. 131 (7): 1485–1490. Bibcode:2003MWRv..131.1485T. CiteSeerX 10.1.1.332.775. doi:10.1175/1520-0493(2003)131<1485:ESRF>2.0.CO;2.
- ↑ Anderson, J. L. (2003). "संयोजन फ़िल्टरिंग के लिए एक स्थानीय न्यूनतम वर्ग ढाँचा". Monthly Weather Review. 131 (4): 634–642. Bibcode:2003MWRv..131..634A. CiteSeerX 10.1.1.10.6543. doi:10.1175/1520-0493(2003)131<0634:ALLSFF>2.0.CO;2.
- ↑ Ott, E.; Hunt, B. R.; Szunyogh, I.; Zimin, A. V.; Kostelich, E. J.; Corazza, M.; Kalnay, E.; Patil, D.; Yorke, J. A. (2004). "वायुमंडलीय डेटा सम्मिलन के लिए एक स्थानीय पहनावा कलमन फ़िल्टर". Tellus A. 56 (5): 415–428. arXiv:physics/0203058. Bibcode:2004TellA..56..415O. doi:10.3402/tellusa.v56i5.14462. S2CID 218577557.
- ↑ Ravela, S.; Emanuel, K.; McLaughlin, D. (2007). "फ़ील्ड संरेखण द्वारा डेटा सम्मिलन". Physica. D: Nonlinear Phenomena. 230 (1–2): 127–145. Bibcode:2007PhyD..230..127R. doi:10.1016/j.physd.2006.09.035.
- ↑ Beezley, J. D.; Mandel, J. (2008). "मॉर्फ़िंग पहनावा कलमन फ़िल्टर". Tellus A. 60 (1): 131–140. arXiv:0705.3693. Bibcode:2008TellA..60..131B. doi:10.1111/j.1600-0870.2007.00275.x. S2CID 1009227.
- ↑ 20.0 20.1 Mandel, J.; Beezley, J. D. (November 2006). उच्च आयामी नॉनलाइनर सिस्टम में विरल डेटा को आत्मसात करने के लिए प्रीडिक्टर-करेक्टर और मॉर्फिंग एसेंबल फिल्टर (PDF). 11th Symposium on Integrated Observing and Assimilation Systems for the Atmosphere, Oceans, and Land Surface (IOAS-AOLS), CD-ROM, Paper 4.12, 87th American Meteorological Society Annual Meeting, San Antonio, TX, January 2007. CCM Report 239. University of Colorado at Denver and Health Sciences Center.
- ↑ Anderson, J. L.; Anderson, S. L. (1999). "सामूहिक सम्मिलन और पूर्वानुमान उत्पन्न करने के लिए नॉनलाइनियर फ़िल्टरिंग समस्या का मोंटे कार्लो कार्यान्वयन". Monthly Weather Review. 127 (12): 2741–2758. Bibcode:1999MWRv..127.2741A. doi:10.1175/1520-0493(1999)127<2741:AMCIOT>2.0.CO;2.
- ↑ Bengtsson, T.; Snyder, C.; Nychka, D. (2003). "उच्च आयामी प्रणालियों के लिए एक अरेखीय संयोजन फ़िल्टर की ओर". Journal of Geophysical Research: Atmospheres. 108 (D24): STS 2–1–10. Bibcode:2003JGRD..108.8775B. doi:10.1029/2002JD002900.
- ↑ van Leeuwen, P. (2003). "बड़े पैमाने के अनुप्रयोगों के लिए एक विचरण-न्यूनीकरण फ़िल्टर". Monthly Weather Review. 131 (9): 2071–2084. Bibcode:2003MWRv..131.2071V. CiteSeerX 10.1.1.7.3719. doi:10.1175/1520-0493(2003)131<2071:AVFFLA>2.0.CO;2.
बाहरी संबंध
- ईएनकेएफ webpage
- TOPAZ, real-time forecasting of the North Atlantic ocean and Arctic sea-ice with the ईएनकेएफ
- ईएनकेएफ-C, a compact framework for data assimilation into large-स्काle layered geophysical models with the ईएनकेएफ
- PDAF – Parallel Data Assimilation Framework – an open-source software for data assimilation providing different variants of the ईएनकेएफ