एन्सेम्बल कलमैन फ़िल्टर: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
Line 1: Line 1:
एन्सेम्बल कल्मन फ़िल्टर (EnKF) [[पुनरावर्ती फ़िल्टर]] है जो बड़ी संख्या में चर वाली समस्याओं के लिए उपयुक्त है, जैसे कि भूभौतिकीय मॉडल में आंशिक अंतर समीकरणों का विवेकीकरण। EnKF की उत्पत्ति बड़ी समस्याओं के लिए [[कलमन फ़िल्टर]] के संस्करण के रूप में हुई (अनिवार्य रूप से, सहप्रसरण मैट्रिक्स को [[नमूना सहप्रसरण मैट्रिक्स]] द्वारा प्रतिस्थापित किया जाता है), और यह अब संयोजन पूर्वानुमान का महत्वपूर्ण [[डेटा आत्मसात]] घटक है। EnKF [[कण फिल्टर]] से संबंधित है (इस संदर्भ में, कण समूह सदस्य के समान है) लेकिन EnKF यह धारणा बनाता है कि इसमें शामिल सभी संभाव्यता वितरण [[सामान्य वितरण]] हैं; जब यह लागू होता है, तो यह कण फिल्टर की तुलना में बहुत अधिक कुशल होता है।
'''एन्सेम्बल कल्मन फ़िल्टर''' (ईएनकेएफ) [[पुनरावर्ती फ़िल्टर|रिकर्सिव फ़िल्टर]] है जो बड़ी संख्या में वैरिएबल वाली समस्याओं के लिए उपयुक्त है, जैसे कि भूभौतिकीय मॉडल में आंशिक अंतर समीकरणों का विवेकीकरण ईएनकेएफ की उत्पत्ति बड़ी समस्याओं के लिए [[कलमन फ़िल्टर]] के संस्करण के रूप में हुई (अनिवार्य रूप से, सहसंयोजकता आव्यूह को [[नमूना सहप्रसरण मैट्रिक्स|प्रारूपिक सहसंयोजकता]] आव्यूह द्वारा प्रतिस्थापित किया जाता है), और यह अब संयोजन पूर्वानुमान का महत्वपूर्ण [[डेटा आत्मसात|डेटा एसीमिलेसन]] घटक है। ईएनकेएफ [[कण फिल्टर|पार्टिकल फिल्टर]] से संबंधित है (इस संदर्भ में, पार्टिकल समूह सदस्य के समान है) किन्तु ईएनकेएफ यह धारणा बनाता है कि इसमें सम्मिलित सभी संभाव्यता वितरण [[सामान्य वितरण|गाऊसी]] हैं; जब यह प्रयुक्त होता है, जिससे यह पार्टिकल फिल्टर की तुलना में बहुत अधिक उत्तम होता है।




Line 5: Line 5:
==परिचय==
==परिचय==


एन्सेम्बल कल्मन फ़िल्टर (EnKF) [[बायेसियन अनुमान]] समस्या का [[मोंटे कार्लो विधि]] कार्यान्वयन है: मॉडल किए गए सिस्टम की स्थिति ([[पूर्व संभावना]], जिसे अक्सर भूविज्ञान में पूर्वानुमान कहा जाता है) और डेटा संभावना की संभाव्यता घनत्व फ़ंक्शन (पीडीएफ) दी गई है। , बेयस प्रमेय का उपयोग डेटा संभावना को ध्यान में रखने के बाद पीडीएफ प्राप्त करने के लिए किया जाता है (पोस्टीरियर संभावना, जिसे अक्सर विश्लेषण कहा जाता है)। इसे बायेसियन अपडेट कहा जाता है। बायेसियन अपडेट को समय-समय पर नए डेटा को शामिल करते हुए मॉडल को आगे बढ़ाने के साथ जोड़ा जाता है। मूल कलमैन फ़िल्टर, 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> मानता है कि सभी पीडीएफ सामान्य वितरण (गॉसियन धारणा) हैं और बायेसियन अपडेट द्वारा माध्य और सहप्रसरण मैट्रिक्स के परिवर्तन के लिए बीजगणितीय सूत्र प्रदान करता है, साथ ही समय में माध्य और सहप्रसरण को आगे बढ़ाने के लिए सूत्र प्रदान करता है, बशर्ते सिस्टम रैखिक हो। हालाँकि, उच्च-आयामी प्रणालियों के लिए सहप्रसरण मैट्रिक्स को बनाए रखना कम्प्यूटेशनल रूप से संभव नहीं है। इस कारण से, EnKFs विकसित किए गए थे।<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> EnKFs राज्य वैक्टरों के संग्रह का उपयोग करके सिस्टम स्थिति के वितरण का प्रतिनिधित्व करते हैं, जिसे संख्यात्मक मौसम पूर्वानुमान#एनसेम्बल कहा जाता है, और सहप्रसरण मैट्रिक्स को समुच्चय से गणना किए गए [[नमूना सहप्रसरण]] द्वारा प्रतिस्थापित करते हैं। समूह को ऐसे संचालित किया जाता है जैसे कि यह [[यादृच्छिक नमूना]] हो, लेकिन समूह के सदस्य वास्तव में [[सांख्यिकीय स्वतंत्रता]] नहीं रखते हैं, क्योंकि वे सभी EnKF साझा करते हैं। EnKFs का फायदा यह है कि पीडीएफ को समय पर आगे बढ़ाने का काम केवल समूह के प्रत्येक सदस्य को आगे बढ़ाना है।<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>
एन्सेम्बल कल्मन फ़िल्टर (ईएनकेएफ) [[बायेसियन अनुमान|बायेसियन अपडेट]] समस्या का [[मोंटे कार्लो विधि]] कार्यान्वयन है: मॉडल किए गए प्रणाली की स्थिति ([[पूर्व संभावना|प्रियोर]], जिसे अधिकांशतः भूविज्ञान में पूर्वानुमान कहा जाता है) और डेटा संभावना की संभाव्यता घनत्व फलन (पीडीएफ) दी गई है। बेयस प्रमेय का उपयोग डेटा संभावना को ध्यान में रखने के पश्चात पीडीएफ प्राप्त करने के लिए किया जाता है (पोस्टीरियर संभावना, जिसे अधिकांशतः विश्लेषण कहा जाता है)। इसे बायेसियन अपडेट कहा जाता है। बायेसियन अपडेट को समय-समय पर नए डेटा को सम्मिलित करते हुए मॉडल को आगे बढ़ाने के साथ जोड़ा जाता है। मूल कलमैन फ़िल्टर, 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>




Line 11: Line 11:


===कलमैन फ़िल्टर===
===कलमैन फ़िल्टर===
{{main|Kalman filter}}
{{main|कलमन फ़िल्टर}}
होने देना <math>\mathbf{x}</math> निरूपित करें <math>n</math> मॉडल का -आयामी राज्य स्थान प्रतिनिधित्व, और मान लें कि इसका माध्य के साथ सामान्य वितरण है <math>\mathbf{\mu}</math> और सहप्रसरण <math>Q</math>यानी कि इसका पीडीएफ है
 
 
मान लीजिए कि <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> ऐसा माना जाता है कि इसमें सहप्रसरण के साथ गॉसियन पीडीएफ है <math>R</math> और मतलब <math>H\mathbf{x}</math>, कहाँ <math>H</math> तथाकथित [[टोपी मैट्रिक्स]] है। सहप्रसरण मैट्रिक्स <math>R</math> डेटा की त्रुटि के अनुमान का वर्णन करता है; यदि डेटा वेक्टर की प्रविष्टियों में यादृच्छिक त्रुटियाँ हैं <math>\mathbf{d}</math> स्वतंत्र हैं, <math>R</math> विकर्ण है और इसकी विकर्ण प्रविष्टियाँ डेटा वेक्टर की संबंधित प्रविष्टियों की त्रुटि के [[मानक विचलन]] ("त्रुटि आकार") के वर्ग हैं <math>\mathbf{d}</math>. मूल्य <math>H\mathbf{x}</math> राज्य के लिए डेटा का मूल्य क्या होगा <math>\mathbf{x}</math> डेटा त्रुटियों के अभाव में. फिर संभाव्यता घनत्व <math>p(\mathbf{d}|\mathbf{x})</math> डेटा का <math>\mathbf{d}</math> सिस्टम स्थिति की सशर्तता <math>\mathbf{x}</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{x}</math> डेटा के मूल्य पर सशर्त <math>\mathbf{d}</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{\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>\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}}</math> और सहप्रसरण <math>\hat{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> \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>
तथाकथित कलमन फिल्टर#कलमैन लाभ व्युत्पत्ति मैट्रिक्स है।
तथाकथित कलमन फिल्टर कलमैन लाभ व्युत्पत्ति आव्यूह है।


===कलामन फ़िल्टर इकट्ठा करें===
===एन्सेम्बल कलमैन फ़िल्टर ===


EnKF कलमैन फ़िल्टर का मोंटे कार्लो सन्निकटन है, जो राज्य वेक्टर के पीडीएफ के सहप्रसरण मैट्रिक्स को विकसित करने से बचाता है <math>\mathbf{x}</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>


<math>X</math> <math>n\times N</math> मैट्रिक्स जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से यादृच्छिक नमूना बनाएंगे। हालाँकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से सांख्यिकीय स्वतंत्रता में नहीं हैं, क्योंकि प्रत्येक EnKF कदम उन्हें साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वे वास्तव में स्वतंत्र हों।
<math>X</math> <math>n\times N</math> आव्यूह जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से यादृच्छिक प्रारूपिक बनाएंगे। चूंकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से सांख्यिकीय स्वतंत्रता में नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ कदम उन्हें साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वे वास्तव में स्वतंत्र हों।
 
X एक <math>n\times N</math> आव्यूह है जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से एक प्रारूप बनाएंगे। चूँकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से स्वतंत्र नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ चरण उन्हें एक साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वह वास्तव में स्वतंत्र होंते है।


डेटा को दोहराएँ <math>\mathbf{d}</math> में <math>m\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>m</math>-आयामी सामान्य वितरण <math>N(0,R)</math>. यदि, इसके अतिरिक्त, के कॉलम <math>X</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>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>.
Line 54: Line 64:
:<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>,


जो पश्च विचरण है।
जो पश्च भिन्नता है।
 
EnKF अब केवल राज्य सहप्रसरण को प्रतिस्थापित करके प्राप्त किया जाता है <math>Q</math> कलमन गेन मैट्रिक्स में <math>K</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>
 


ईएनकेएफ अब कल्मन गेन आव्यूह <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>X</math> और डेटा आव्यूह <math>D</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> 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}</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>D</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>और व्युत्क्रम को छद्म व्युत्क्रम से बदल दिया जाता है, जिसकी गणना एकवचन-मूल्य अपघटन (एसवीडी) का उपयोग करके की जाती है।
ध्यान दें कि तब से <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> और व्युत्क्रम को छद्म व्युत्क्रम से परिवर्तित कर दिया जाता है, जिसकी गणना विलक्षण मान अपघटन (एसवीडी) का उपयोग करके की जाती है।


चूँकि ये सूत्र प्रमुख BLAS#स्तर 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"/>
चूँकि ये सूत्र प्रमुख ब्लास स्तर 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</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>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 98: Line 106:


:<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> स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह फार्मूला भी लागू है<ref name="Mandel-2006-EIE"/> अवलोकन समारोह के लिए <math>h(\mathbf{x})=H\mathbf{x+f}</math> निश्चित ऑफसेट के साथ <math>\mathbf{f}</math>, जिसे स्पष्ट रूप से जानने की भी आवश्यकता नहीं है। उपरोक्त सूत्र का उपयोग आमतौर पर अरेखीय अवलोकन फ़ंक्शन के लिए किया गया है <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>h</math> प्रत्येक समूह सदस्य पर बार और आव्यूह <math>H</math> स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह फार्मूला भी प्रयुक्त है<ref name="Mandel-2006-EIE"/> अवलोकन समारोह के लिए <math>h(\mathbf{x})=H\mathbf{x+f}</math> निश्चित ऑफसेट के साथ <math>\mathbf{f}</math>, जिसे स्पष्ट रूप से जानने की भी आवश्यकता नहीं है। उपरोक्त सूत्र का उपयोग सामान्यतः अरेखीय अवलोकन फलन के लिए किया गया है <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>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" /> उस स्थिति में, अवलोकन फलन को अनिवार्य रूप से समूह सदस्यों पर इसके मूल्यों से एक रैखिक फलन द्वारा अनुमानित किया जाता है।


===बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन===
===बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन===


बड़ी संख्या के लिए <math>m</math> डेटा बिंदुओं का, से गुणा <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>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 111: Line 121:
:<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>N</math> साथ <math>m</math> दाहिनी ओर. देखना<ref name="Mandel-2006-EIE"/>ऑपरेशन गणना के लिए.
 
 
जिसके लिए केवल आव्यूह <math>R</math> (सस्ता माना जाता है) और <math>m</math> दाहिनी ओर के आकार <math>N</math> के प्रणाली के समाधान की आवश्यकता होती है। ऑपरेशन गणना के लिए <ref name="Mandel-2006-EIE" /> देखें।


==आगे का विस्तार==
==आगे का विस्तार==


यहां वर्णित EnKF संस्करण में डेटा का यादृच्छिककरण शामिल है। डेटा के यादृच्छिकीकरण के बिना फ़िल्टर के लिए, देखें।<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-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="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>


गैर-रेखीय समस्याओं के लिए, EnKF गैर-भौतिक अवस्थाओं के साथ पश्च संयोजन बना सकता है। इसे [[तिखोनोव नियमितीकरण]] द्वारा कम किया जा सकता है, जैसे कि बड़े स्थानिक [[ ग्रेडियेंट |ग्रेडियेंट]] वाले राज्यों की दंड विधि।<ref name="Johns-2005-CEK"/>
औपचारिक रूप से, ईएनकेएफ गॉसियन धारणा पर विश्वास करते हैं। व्यवहार में इनका उपयोग अरैखिक समस्याओं के लिए भी किया जा सकता है, जहां गॉसियन धारणा संतुष्ट नहीं हो सकती है। एनकेएफ में गॉसियन धारणा को शिथिल करने का प्रयास करने वाले संबंधित फिल्टर, इसके लाभ को संरक्षित करते हुए, इसमें ऐसे फिल्टर सम्मिलित हैं जो विभिन्न गॉसियन कर्नेल के साथ स्टेट पीडीएफ को फिट करते हैं,<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>


[[सुसंगत विशेषता]]ओं वाली समस्याओं, जैसे कि तूफान, तूफान, [[ आग रेखा |आग रेखा]] , [[ तूफान रेखा |तूफान रेखा]] और बारिश के मोर्चों के लिए, अंतरिक्ष (इसके ग्रिड) में राज्य को विकृत करके और साथ ही राज्य के आयामों को सही करके संख्यात्मक मॉडल स्थिति को समायोजित करने की आवश्यकता है। . 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>
औपचारिक रूप से, EnKF गॉसियन धारणा पर भरोसा करते हैं। व्यवहार में इनका उपयोग अरैखिक समस्याओं के लिए भी किया जा सकता है, जहां गॉसियन धारणा संतुष्ट नहीं हो सकती है। एनकेएफ में गॉसियन धारणा को शिथिल करने का प्रयास करने वाले संबंधित फिल्टर, इसके फायदों को संरक्षित करते हुए, इसमें ऐसे फिल्टर शामिल हैं जो कई गॉसियन कर्नेल के साथ राज्य पीडीएफ को फिट करते हैं,<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>




==यह भी देखें==
==यह भी देखें==
* डेटा आत्मसात
* डेटा एसीमिलेसन
*संख्यात्मक मौसम भविष्यवाणी#समूह
*न्यूमेरिकल वेदर प्रेडिक्शन एसेम्बल
* कण फिल्टर
* पार्टिकल फिल्टर
* [[पुनरावर्ती बायेसियन अनुमान]]
* [[पुनरावर्ती बायेसियन अनुमान|रिकर्सिव बायेसियन एसटीमेसन]]


==संदर्भ==
==संदर्भ==
Line 135: Line 148:


==बाहरी संबंध==
==बाहरी संबंध==
* [http://enkf.nersc.no EnKF webpage]
* [http://enkf.nersc.no ईएनकेएफ webpage]
* [http://topaz.nersc.no TOPAZ, real-time forecasting of the North Atlantic ocean and Arctic sea-ice with the EnKF]
* [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 EnKF-C, a compact framework for data assimilation into large-scale layered geophysical models with the EnKF]
* [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 EnKF
* [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: मोंटे कार्लो विधियाँ]]  

Revision as of 18:52, 23 November 2023

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


परिचय

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


व्युत्पत्ति

कलमैन फ़िल्टर


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


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


स्थिति की पीडीएफ और डेटा संभावना को बेयस प्रमेय द्वारा डेटा के मूल्य पर नियमबद्ध प्रणाली स्थिति की नई संभाव्यता घनत्व देने के लिए संयोजित किया गया है,

आंकड़ा एक बार प्राप्त होने के पश्चात यह निश्चित हो जाता है, इसलिए पिछली स्थिति को इससे निरूपित करें के अतिरिक्त और पिछला पीडीएफ़ द्वारा इसे बीजगणितीय जोड़ द्वारा दिखाया जा सकता है [5] पिछला पीडीएफ भी गॉसियन है,

कल्मन अपडेट फ़ार्मुलों द्वारा दिए गए पश्च माध्य @ और सहसंयोजकता @ के साथ किया जाता है

जहाँ

तथाकथित कलमन फिल्टर कलमैन लाभ व्युत्पत्ति आव्यूह है।

एन्सेम्बल कलमैन फ़िल्टर

ईएनकेएफ कलमैन फ़िल्टर का मोंटे कार्लो सन्निकटन है, जो स्टेट वेक्टर के पीडीएफ के सहसंयोजकता आव्यूह को विकसित करने से बचाता है . इसके अतिरिक्त, पीडीएफ को समूह द्वारा दर्शाया जाता है

आव्यूह जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से यादृच्छिक प्रारूपिक बनाएंगे। चूंकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से सांख्यिकीय स्वतंत्रता में नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ कदम उन्हें साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वे वास्तव में स्वतंत्र हों।

X एक आव्यूह है जिसके कॉलम समूह सदस्य हैं, और इसे पूर्व समूह कहा जाता है। आदर्श रूप से, समूह के सदस्य पूर्व वितरण से एक प्रारूप बनाएंगे। चूँकि, प्रारंभिक समूह को छोड़कर समूह के सदस्य सामान्य रूप से स्वतंत्र नहीं हैं, क्योंकि प्रत्येक ईएनकेएफ चरण उन्हें एक साथ जोड़ता है। उन्हें लगभग स्वतंत्र माना जाता है, और सभी गणनाएँ ऐसे आगे बढ़ती हैं मानो वह वास्तव में स्वतंत्र होंते है।

डेटा को आव्यूह में दोहराएं जाते है


जिससे प्रत्येक कॉलम में डेटा वेक्टर प्लस एम-आयामी सामान्य वितरण से एक यादृच्छिक वेक्टर सम्मिलित होते है। यदि इसके अतिरिक्त, के कॉलम पूर्व संभाव्यता वितरण से एक प्रारूप हैं,

पश्च संभाव्यता वितरण से एक प्रारूप बनाएं। इसे के साथ अदिश स्थिति में देखने के लिए: मान लीजिए , और तब

.

पहला योग पश्च माध्य है, और दूसरे योग में, स्वतंत्रता को ध्यान में रखते हुए, भिन्नता है

,

जो पश्च भिन्नता है।

ईएनकेएफ अब कल्मन गेन आव्यूह में स्टेट सहसंयोजकता को संयोजन सदस्यों से गणना किए गए प्रारूप सहसंयोजकता (जिसे समुच्चय सहसंयोजकता कहा जाता है) द्वारा प्रतिस्थापित करके प्राप्त किया जाता है, [6] जो कि है

कार्यान्वयन

मूल सूत्रीकरण

यहां हम अनुसरण करते हैं।[7][8] मान लीजिए कि संयोजन आव्यूह और डेटा आव्यूह उपरोक्तानुसार हैं समुच्चय माध्य और सहसंयोजकता हैं

जहाँ

और संकेतित आकार के सभी के आव्यूह को दर्शाता है।

इसके पश्चात पिछला एन्सेम्बल द्वारा दिया जाता है

जहां विकृत डेटा आव्यूह ऊपर जैसा है।

ध्यान दें कि तब से सहसंयोजकता आव्यूह है, यह सदैव धनात्मक अर्धनिश्चित आव्यूह होता है और सामान्यतः धनात्मक अर्धनिश्चित आव्यूह होता है, इसलिए उपरोक्त व्युत्क्रम उपस्थित है और सूत्र कोलेस्की अपघटन द्वारा कार्यान्वित किया जा सकता है।[9][7][8] प्रारूपिक सहसंयोजकता द्वारा प्रतिस्थापित किया जाता है जहाँ और व्युत्क्रम को छद्म व्युत्क्रम से परिवर्तित कर दिया जाता है, जिसकी गणना विलक्षण मान अपघटन (एसवीडी) का उपयोग करके की जाती है।

चूँकि ये सूत्र प्रमुख ब्लास स्तर 3 ऑपरेशन वाले आव्यूह ऑपरेशन हैं,[10] वह लैपैक (सीरियल और शेयर्ड मेमोरी कंप्यूटर पर) और स्कालैपैक (वितरित मेमोरी कंप्यूटर पर) जैसे सॉफ़्टवेयर पैकेजों का उपयोग करके कुशल कार्यान्वयन के लिए उपयुक्त हैं।[9] किसी आव्यूह के व्युत्क्रम आव्यूह की गणना करने और उससे गुणा करने के अतिरिक्त विपरीत आव्यूह के चॉलेस्की अपघटन की गणना करना और व्युत्क्रम द्वारा गुणन को विभिन्न के साथ रैखिक प्रणाली के समाधान के रूप में मानना ​​​​बहुत उत्तम (विभिन्न गुना सस्ता और अधिक स्पष्ट) है [10]


अवलोकन आव्यूह-मुक्त कार्यान्वयन

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

प्रोग्राम इसे अवलोकन फलन कहा जाता है या, व्युत्क्रम समस्याओं के संदर्भ में, व्युत्क्रम समस्या व्युत्क्रम समस्याओं का संभाव्य सूत्रीकरण कहा जाता है। का मान है स्टेट के लिए डेटा का मूल्य क्या होगा यह मानते हुए कि माप स्पष्ट है। फिर पश्च संयोजन को पुनः लिखा जा सकता है

जहाँ

और

साथ

परिणाम स्वरुप, अवलोकन फलन का मूल्यांकन करके संयोजन अपडेट की गणना की जा सकती है प्रत्येक समूह सदस्य पर बार और आव्यूह स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह फार्मूला भी प्रयुक्त है[9] अवलोकन समारोह के लिए निश्चित ऑफसेट के साथ , जिसे स्पष्ट रूप से जानने की भी आवश्यकता नहीं है। उपरोक्त सूत्र का उपयोग सामान्यतः अरेखीय अवलोकन फलन के लिए किया गया है , जैसे कि तूफान भंवर की स्थिति।[11] उस स्थिति में, अवलोकन फलन को अनिवार्य रूप से समूह सदस्यों पर इसके मूल्यों से रैखिक फलन द्वारा अनुमानित किया जाता है।

परिणाम स्वरुप, प्रत्येक संयोजन सदस्य पर एक बार अवलोकन फलन का मूल्यांकन करके संयोजन अद्यतन की गणना की जा सकती है और आव्यूह को स्पष्ट रूप से जानने की आवश्यकता नहीं है। यह सूत्र एक निश्चित ऑफसेट के साथ एक अवलोकन फलन के लिए भी [9] प्रयुक्त करता है, जिसे जानने की भी आवश्यकता नहीं है स्पष्ट रूप से. उपरोक्त सूत्र का उपयोग सामान्यतः एक गैर-रेखीय अवलोकन फलन के लिए किया गया है, जैसे कि तूफान भंवर की स्थिति [11] उस स्थिति में, अवलोकन फलन को अनिवार्य रूप से समूह सदस्यों पर इसके मूल्यों से एक रैखिक फलन द्वारा अनुमानित किया जाता है।

बड़ी संख्या में डेटा बिंदुओं के लिए कार्यान्वयन

बड़ी संख्या में डेटा बिंदुओं के लिए, से गुणा एक बाधा बन जाता है। निम्न वैकल्पिक सूत्र तब लाभप्रद होता है जब डेटा बिंदुओं की संख्या बड़ी होती है (जैसे कि ग्रिड या पिक्सेल डेटा को एसिमिलेसन करते समय) और डेटा त्रुटि सहसंयोजकता आव्यूह विकर्ण होता है (जो कि तब होता है जब डेटा त्रुटियाँ असंबद्ध होती हैं), या सस्ता होता है विघटित (जैसे कि सीमित सहसंयोजकता दूरी के कारण बैंडेड) शर्मन-मॉरिसन-वुडबरी फॉर्मूला का उपयोग करता है [12]

साथ

देता है


जिसके लिए केवल आव्यूह (सस्ता माना जाता है) और दाहिनी ओर के आकार के प्रणाली के समाधान की आवश्यकता होती है। ऑपरेशन गणना के लिए [9] देखें।

आगे का विस्तार

यहां वर्णित ईएनकेएफ संस्करण में डेटा का यादृच्छिककरण सम्मिलित है। डेटा के यादृच्छिकीकरण के बिना फ़िल्टर के लिए, देखें।[13][14][15] चूंकि समूह सहसंयोजकता में रैंक डेफीसिएंट है (समूह सदस्यों की तुलना में विभिन्न अधिक स्टेट वैरिएबल हैं, सामान्यतः लाखों, सामान्यतः सौ से कम), इसमें उन बिंदुओं के जोड़े के लिए बड़े शब्द हैं जो स्थानिक रूप से दूर हैं। चूँकि वास्तव में दूर के स्थानों पर भौतिक क्षेत्रों के मान इतने सहसंबद्ध नहीं हैं, दूरी के आधार पर सहसंयोजकता आव्यूह को कृत्रिम रूप से कम कर दिया जाता है, जो स्थानीयकृत कलमैन फ़िल्टर एल्गोरिदम को जन्म देता है।[16][17] ये विधियां गणनाओं में उपयोग किए जाने वाले सहसंयोजकता आव्यूह को संशोधित करती हैं और परिणामस्वरूप, पोस्टीरियर एन्सेम्बल अब केवल पोस्टीरियर संयोजन के रैखिक संयोजनों से नहीं बना है।

गैर-रेखीय समस्याओं के लिए, ईएनकेएफ गैर-भौतिक अवस्थाओं के साथ पश्च संयोजन बना सकता है। इसे नियमितीकरण द्वारा कम किया जा सकता है, जैसे कि बड़े स्थानिक ग्रेडियेंट वाले स्टेट की दंड विधि है।[6]

कोहेरेंट विशेषताओं वाली समस्याओं, जैसे कि हरिकेन, थंडरस्टॉर्म, फायरलाइन , स्क्वॉल लाइन और रेन फ्रंट के लिए, अंतरिक्ष (इसके ग्रिड) में स्टेट को विकृत करके और साथ ही स्टेट के आयामों को सही करके संख्यात्मक मॉडल स्थिति को समायोजित करने की आवश्यकता है। 2007 में, रवेला एट अल एनसेंबल का उपयोग करके संयुक्त स्थिति-आयाम समायोजन मॉडल प्रस्तुत करें, और व्यवस्थित रूप से अनुक्रमिक सन्निकटन प्राप्त करें जिसे एनकेएफ और अन्य फॉर्मूलेशन दोनों पर प्रयुक्त किया जा सकता है।[18] उनकी पद्धति यह धारणा नहीं बनाती है कि आयाम और स्थिति त्रुटियां स्वतंत्र या संयुक्त रूप से गॉसियन हैं, जैसा कि अन्य करते हैं। एनकेएफ स्टेट के रैखिक संयोजनों के अतिरिक्त, इमेज रजिस्ट्रेसन और मॉर्फिंग से उधार ली गई तकनीकों द्वारा प्राप्त मध्यवर्ती स्टेट को नियोजित करता है।[19][20]

औपचारिक रूप से, ईएनकेएफ गॉसियन धारणा पर विश्वास करते हैं। व्यवहार में इनका उपयोग अरैखिक समस्याओं के लिए भी किया जा सकता है, जहां गॉसियन धारणा संतुष्ट नहीं हो सकती है। एनकेएफ में गॉसियन धारणा को शिथिल करने का प्रयास करने वाले संबंधित फिल्टर, इसके लाभ को संरक्षित करते हुए, इसमें ऐसे फिल्टर सम्मिलित हैं जो विभिन्न गॉसियन कर्नेल के साथ स्टेट पीडीएफ को फिट करते हैं,[21] फ़िल्टर जो गाऊसी मिश्रण द्वारा स्टेट पीडीएफ का अनुमान लगाते हैं,[22] घनत्व अनुमान द्वारा पार्टिकल भार की गणना के साथ पार्टिकल फिल्टर का प्रकार,[20] और पार्टिकल फ़िल्टर अनुक्रमिक महत्व पुनः प्रारूपिककरण (एसआईआर) को कम करने के लिए कॉची वितरण डेटा पीडीएफ के साथ पार्टिकल फ़िल्टर का प्रकार है।[23]


यह भी देखें

संदर्भ

  1. Kalman, R. E. (1960). "रेखीय छानने और भविष्यवाणी की समस्याओं के लिए एक नया दृष्टिकोण". Journal of Basic Engineering. 82 (1): 35–45. doi:10.1115/1.3662552. S2CID 1242324.
  2. Evensen, G. (1994). "त्रुटि आँकड़ों का पूर्वानुमान लगाने के लिए मोंटे कार्लो विधियों का उपयोग करके गैर-रेखीय अर्ध-जियोस्ट्रोफिक मॉडल के साथ अनुक्रमिक डेटा आत्मसात". Journal of Geophysical Research. 99 (C5): 143–162. Bibcode:1994JGR....9910143E. doi:10.1029/94JC00572. hdl:1956/3035.
  3. 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.
  4. 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.
  5. Anderson, B. D. O.; Moore, J. B. (1979). इष्टतम फ़िल्टरिंग. Englewood Cliffs, NJ: Prentice-Hall. ISBN 978-0-13-638122-8.
  6. 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. 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. 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. 9.0 9.1 9.2 9.3 9.4 Mandel, J. (June 2006). "एन्सेम्बल कलमन फ़िल्टर का कुशल कार्यान्वयन" (PDF). Center for Computational Mathematics Reports. University of Colorado at Denver and Health Sciences Center. 231.
  10. 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.
  11. 11.0 11.1 Chen, Y.; Snyder, C. (2007). "एन्सेम्बल कलमन फ़िल्टर के साथ भंवर स्थिति को आत्मसात करना". Monthly Weather Review. 135 (5): 1828–1845. Bibcode:2007MWRv..135.1828C. doi:10.1175/MWR3351.1.
  12. Hager, W. W. (1989). "मैट्रिक्स के व्युत्क्रम को अद्यतन करना". SIAM Review. 31 (2): 221–239. doi:10.1137/1031049.
  13. 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.
  14. 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.
  15. 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.
  16. 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.
  17. 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.
  18. 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.
  19. 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. 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.
  21. 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.
  22. 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.
  23. 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.


बाहरी संबंध