बहुस्तरीय मोंटे कार्लो विधि: Difference between revisions

From Vigyanwiki
(Created page with "संख्यात्मक विश्लेषण में बहुस्तरीय मोंटे कार्लो (MLMC) विधियाँ स्ट...")
 
No edit summary
Line 1: Line 1:
[[संख्यात्मक विश्लेषण]] में बहुस्तरीय मोंटे कार्लो (MLMC) विधियाँ [[स्टोचैस्टिक सिमुलेशन]] में उत्पन्न होने वाले [[अपेक्षित मूल्य]]ों की गणना के लिए [[कलन विधि]] हैं। [[मोंटे कार्लो विधि]]यों की तरह, वे बार-बार सरल यादृच्छिक नमूने पर भरोसा करते हैं, लेकिन इन नमूनों को सटीकता के विभिन्न स्तरों पर लिया जाता है। MLMC विधियाँ कम सटीकता और इसी कम लागत के साथ अधिकांश नमूने लेकर मानक मोंटे कार्लो विधियों की कम्प्यूटेशनल लागत को बहुत कम कर सकती हैं, और केवल बहुत कम नमूने उच्च सटीकता और इसी उच्च लागत पर लिए जाते हैं।
[[संख्यात्मक विश्लेषण]] में बहुस्तरीय मोंटे कार्लो (एमएलएमसी) विधियाँ [[स्टोचैस्टिक सिमुलेशन|संयोजनात्मक सिमुलेशन]] में उत्पन्न होने वाले [[अपेक्षित मूल्य|अपेक्षित मूल्यों]] की गणना के लिए एक [[कलन विधि]] हैं। [[मोंटे कार्लो विधि]]यों की तरह, बहुस्तरीय मोंटे कार्लो विधियाँ भी दोहरे प्रक्रिया आधारित यादृच्छिक प्रतिरूप चयन पर आधारित होती हैं, परंतु इन प्रतिरूपो को विभिन्न सत्यता स्तरों पर लिया जाता है। एमएलएमसी विधियाँ मुख्य रूप से मानक मोंटे कार्लो विधियों की गणना के गणितीय लागत को अत्यधिक कम कर सकती हैं, क्योंकि इसमें अधिकांश प्रतिरूपो को कम सत्यता और उसके संबंधित कम लागत के साथ लिया जाता है, और मात्र बहुत कम संख्या में प्रतिरूपो को उच्च सत्यता और उसके संबंधित उच्च लागत के साथ लिया जाता है।


== लक्ष्य ==
== लक्ष्य ==
Line 6: Line 6:
{{align|center|<math> \operatorname{E}[G_{L}] = \operatorname{E}[G_{0}] + \sum_{\ell=1}^L \operatorname{E}[G_\ell - G_{\ell-1}],</math>}}
{{align|center|<math> \operatorname{E}[G_{L}] = \operatorname{E}[G_{0}] + \sum_{\ell=1}^L \operatorname{E}[G_\ell - G_{\ell-1}],</math>}}


अपेक्षा ऑपरेटर की रैखिकता के कारण यह तुच्छ रूप से संतुष्ट है। हर एक उम्मीद <math>\operatorname{E}[G_\ell - G_{\ell-1}]</math> इसके बाद मोंटे कार्लो विधि द्वारा अनुमान लगाया जाता है, जिसके परिणामस्वरूप बहुस्तरीय मोंटे कार्लो विधि होती है। ध्यान दें कि अंतर का एक नमूना लेना <math>G_\ell - G_{\ell-1}</math> स्तर पर <math>\ell</math> दोनों के अनुकरण की आवश्यकता है <math>G_{\ell}</math> और <math>G_{\ell-1}</math>.
अपेक्षा ऑपरेटर की रैखिकता के कारण यह तुच्छ रूप से संतुष्ट है। हर एक उम्मीद <math>\operatorname{E}[G_\ell - G_{\ell-1}]</math> इसके बाद मोंटे कार्लो विधि द्वारा अनुमान लगाया जाता है, जिसके परिणामस्वरूप बहुस्तरीय मोंटे कार्लो विधि होती है। ध्यान दें कि अंतर का एक प्रतिरूप  लेना <math>G_\ell - G_{\ell-1}</math> स्तर पर <math>\ell</math> दोनों के अनुकरण की आवश्यकता है <math>G_{\ell}</math> और <math>G_{\ell-1}</math>.


एमएलएमसी विधि काम करती है अगर भिन्नताएं <math>\operatorname{V}[G_\ell - G_{\ell-1}]\rightarrow0</math> जैसा <math>\ell\rightarrow\infty</math>, जो कि दोनों के मामले में होगा <math>G_{\ell}</math> और <math>G_{\ell-1}</math> लगभग एक ही यादृच्छिक चर <math>G</math>. [[केंद्रीय सीमा प्रमेय]] द्वारा, इसका तात्पर्य है कि अंतर की अपेक्षा को सटीक रूप से अनुमानित करने के लिए किसी को कम और कम नमूनों की आवश्यकता होती है <math>G_\ell - G_{\ell-1}</math> जैसा <math>\ell\rightarrow\infty</math>. इसलिए ज्यादातर सैंपल लेवल पर ही लिए जाएंगे <math>0</math>, जहां नमूने सस्ते हैं, और बेहतरीन स्तर पर बहुत कम नमूनों की आवश्यकता होगी <math>L</math>. इस अर्थ में, MLMC को एक पुनरावर्ती नियंत्रण भिन्न रणनीति के रूप में माना जा सकता है।
एमएलएमसी विधि काम करती है अगर भिन्नताएं <math>\operatorname{V}[G_\ell - G_{\ell-1}]\rightarrow0</math> जैसा <math>\ell\rightarrow\infty</math>, जो कि दोनों के मामले में होगा <math>G_{\ell}</math> और <math>G_{\ell-1}</math> लगभग एक ही यादृच्छिक चर <math>G</math>. [[केंद्रीय सीमा प्रमेय]] द्वारा, इसका तात्पर्य है कि अंतर की अपेक्षा को सटीक रूप से अनुमानित करने के लिए किसी को कम और कम प्रतिरूपो  की आवश्यकता होती है <math>G_\ell - G_{\ell-1}</math> जैसा <math>\ell\rightarrow\infty</math>. इसलिए ज्यादातर सैंपल लेवल पर ही लिए जाएंगे <math>0</math>, जहां नमूने सस्ते हैं, और बेहतरीन स्तर पर बहुत कम प्रतिरूपो  की आवश्यकता होगी <math>L</math>. इस अर्थ में, एमएलएमसी को एक पुनरावर्ती नियंत्रण भिन्न रणनीति के रूप में माना जा सकता है।


== अनुप्रयोग ==
== अनुप्रयोग ==
[[File:Multilevel monte carlo sample paths for an SDE.png|thumb|विभिन्न स्तरों पर एक SDE के नमूना पथ का अनुमान।|सही]]MLMC के पहले आवेदन का श्रेय माइक जाइल्स को दिया जाता है,<ref>{{cite journal |last=Giles |first=M. B. |date=2008 |title=बहुस्तरीय मोंटे कार्लो पथ सिमुलेशन|journal=Operations Research |volume=56 |issue=3 |pages=607–617|url=https://ora.ox.ac.uk/objects/uuid:d9d28973-94aa-4179-962a-28bcfa8d8f00/datastreams/ATTACHMENT01 |doi=10.1287/opre.1070.0496|citeseerx=10.1.1.121.713 |s2cid=3000492 }}</ref> [[मोंटे कार्लो विकल्प मॉडल]] के लिए [[स्टोचैस्टिक अंतर समीकरण]] (एसडीई) के संदर्भ में, हालांकि, पैरामीट्रिक एकीकरण के संदर्भ में हेनरिक के काम में पहले के निशान पाए जाते हैं।<ref>{{cite journal |last=Heinrich |first=S. |date=2001 |title=बहुस्तरीय मोंटे कार्लो तरीके|journal=Lecture Notes in Computer Science (Multigrid Methods) |series=Lecture Notes in Computer Science |publisher=Springer |volume=2179 |pages=58–67|doi=10.1007/3-540-45346-6_5 |isbn=978-3-540-43043-8 }}</ref> यहाँ, यादृच्छिक चर <math>G=f(X(T))</math> अदायगी समारोह, और सन्निकटन के अनुक्रम के रूप में जाना जाता है <math>G_\ell</math>, <math>\ell=0,\ldots,L</math> नमूना पथ के सन्निकटन का उपयोग करें <math>X(t)</math> समय कदम के साथ <math>h_\ell=2^{-\ell}T</math>.
[[File:Multilevel monte carlo sample paths for an SDE.png|thumb|विभिन्न स्तरों पर एक SDE के नमूना पथ का अनुमान।|सही]]एमएलएमसी के पहले आवेदन का श्रेय माइक जाइल्स को दिया जाता है,<ref>{{cite journal |last=Giles |first=M. B. |date=2008 |title=बहुस्तरीय मोंटे कार्लो पथ सिमुलेशन|journal=Operations Research |volume=56 |issue=3 |pages=607–617|url=https://ora.ox.ac.uk/objects/uuid:d9d28973-94aa-4179-962a-28bcfa8d8f00/datastreams/ATTACHMENT01 |doi=10.1287/opre.1070.0496|citeseerx=10.1.1.121.713 |s2cid=3000492 }}</ref> [[मोंटे कार्लो विकल्प मॉडल]] के लिए [[स्टोचैस्टिक अंतर समीकरण]] (एसडीई) के संदर्भ में, हालांकि, पैरामीट्रिक एकीकरण के संदर्भ में हेनरिक के काम में पहले के निशान पाए जाते हैं।<ref>{{cite journal |last=Heinrich |first=S. |date=2001 |title=बहुस्तरीय मोंटे कार्लो तरीके|journal=Lecture Notes in Computer Science (Multigrid Methods) |series=Lecture Notes in Computer Science |publisher=Springer |volume=2179 |pages=58–67|doi=10.1007/3-540-45346-6_5 |isbn=978-3-540-43043-8 }}</ref> यहाँ, यादृच्छिक चर <math>G=f(X(T))</math> अदायगी समारोह, और सन्निकटन के अनुक्रम के रूप में जाना जाता है <math>G_\ell</math>, <math>\ell=0,\ldots,L</math> प्रतिरूप  पथ के सन्निकटन का उपयोग करें <math>X(t)</math> समय कदम के साथ <math>h_\ell=2^{-\ell}T</math>.


अनिश्चितता परिमाणीकरण (यूक्यू) में समस्याओं के लिए एमएलएमसी का अनुप्रयोग अनुसंधान का एक सक्रिय क्षेत्र है।<ref>{{cite journal |last1=Cliffe |first1=A. |last2=Giles |first2=M. B. |last3=Scheichl |first3=R. |last4=Teckentrup |first4=A.|author4-link=Aretha Teckentrup |date=2011 |title=बहुस्तरीय मोंटे कार्लो के तरीके और रैंडम गुणांक वाले अण्डाकार पीडीई के अनुप्रयोग|journal=Computing and Visualization in Science |volume=14 |issue=1 |pages=3–15|url=http://people.maths.ox.ac.uk/~gilesm/files/cgst.pdf |doi=10.1007/s00791-011-0160-x|s2cid=1687254 }}</ref><ref>{{cite journal |last1=Pisaroni |first1=M. |last2=Nobile |first2=F. B. |last3=Leyland |first3=P. |date=2017 |title=कंप्रेसिबल इनविसिड एरोडायनामिक्स में अनिश्चितता मात्रा के लिए एक निरंतरता बहु स्तरीय मोंटे कार्लो विधि|journal=Computer Methods in Applied Mechanics and Engineering |volume=326 |issue=C |pages=20–50|url=https://pdfs.semanticscholar.org/6dbc/8dde601b1757c42a4c54fa9cfd69317c82c8.pdf |archive-url=https://web.archive.org/web/20180214142033/https://pdfs.semanticscholar.org/6dbc/8dde601b1757c42a4c54fa9cfd69317c82c8.pdf |url-status=dead |archive-date=2018-02-14 |doi=10.1016/j.cma.2017.07.030|s2cid=10379943 }}</ref> इन समस्याओं का एक महत्वपूर्ण प्रोटोटाइपिकल उदाहरण [[आंशिक अंतर समीकरण]] (पीडीई) हैं जो [[स्टोकेस्टिक आंशिक अंतर समीकरण]] के साथ हैं। इस संदर्भ में, यादृच्छिक चर <math>G</math> ब्याज की मात्रा के रूप में जाना जाता है, और सन्निकटन का क्रम अलग-अलग जाल आकारों के साथ पीडीई के [[विवेक]] से मेल खाता है।
अनिश्चितता परिमाणीकरण (यूक्यू) में समस्याओं के लिए एमएलएमसी का अनुप्रयोग अनुसंधान का एक सक्रिय क्षेत्र है।<ref>{{cite journal |last1=Cliffe |first1=A. |last2=Giles |first2=M. B. |last3=Scheichl |first3=R. |last4=Teckentrup |first4=A.|author4-link=Aretha Teckentrup |date=2011 |title=बहुस्तरीय मोंटे कार्लो के तरीके और रैंडम गुणांक वाले अण्डाकार पीडीई के अनुप्रयोग|journal=Computing and Visualization in Science |volume=14 |issue=1 |pages=3–15|url=http://people.maths.ox.ac.uk/~gilesm/files/cgst.pdf |doi=10.1007/s00791-011-0160-x|s2cid=1687254 }}</ref><ref>{{cite journal |last1=Pisaroni |first1=M. |last2=Nobile |first2=F. B. |last3=Leyland |first3=P. |date=2017 |title=कंप्रेसिबल इनविसिड एरोडायनामिक्स में अनिश्चितता मात्रा के लिए एक निरंतरता बहु स्तरीय मोंटे कार्लो विधि|journal=Computer Methods in Applied Mechanics and Engineering |volume=326 |issue=C |pages=20–50|url=https://pdfs.semanticscholar.org/6dbc/8dde601b1757c42a4c54fa9cfd69317c82c8.pdf |archive-url=https://web.archive.org/web/20180214142033/https://pdfs.semanticscholar.org/6dbc/8dde601b1757c42a4c54fa9cfd69317c82c8.pdf |url-status=dead |archive-date=2018-02-14 |doi=10.1016/j.cma.2017.07.030|s2cid=10379943 }}</ref> इन समस्याओं का एक महत्वपूर्ण प्रोटोटाइपिकल उदाहरण [[आंशिक अंतर समीकरण]] (पीडीई) हैं जो [[स्टोकेस्टिक आंशिक अंतर समीकरण]] के साथ हैं। इस संदर्भ में, यादृच्छिक चर <math>G</math> ब्याज की मात्रा के रूप में जाना जाता है, और सन्निकटन का क्रम अलग-अलग जाल आकारों के साथ पीडीई के [[विवेक]] से मेल खाता है।
Line 20: Line 20:
दोहराना
दोहराना
     स्तर पर वार्म-अप के नमूने लें <math>L</math>
     स्तर पर वार्म-अप के नमूने लें <math>L</math>
सभी स्तरों पर नमूना प्रसरण की गणना करें <math>\ell=0,\ldots,L</math>
सभी स्तरों पर प्रतिरूप  प्रसरण की गणना करें <math>\ell=0,\ldots,L</math>
नमूनों की इष्टतम संख्या को परिभाषित करें <math>N_\ell</math> सभी स्तरों पर <math>\ell=0,\ldots,L</math>
प्रतिरूपो  की इष्टतम संख्या को परिभाषित करें <math>N_\ell</math> सभी स्तरों पर <math>\ell=0,\ldots,L</math>
प्रत्येक स्तर पर अतिरिक्त नमूने लें <math>\ell</math> के अनुसार <math>N_\ell</math>
प्रत्येक स्तर पर अतिरिक्त नमूने लें <math>\ell</math> के अनुसार <math>N_\ell</math>
अगर <math>L\geq2</math> तब
अगर <math>L\geq2</math> तब

Revision as of 10:05, 1 July 2023

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

लक्ष्य

बहुस्तरीय मोंटे कार्लो पद्धति का लक्ष्य अपेक्षित मूल्य का अनुमान लगाना है यादृच्छिक चर का यह एक स्टोकेस्टिक सिमुलेशन का आउटपुट है। मान लीजिए कि यह यादृच्छिक चर बिल्कुल अनुकरण नहीं किया जा सकता है, लेकिन सन्निकटन का एक क्रम है बढ़ती सटीकता के साथ, लेकिन बढ़ती लागत के साथ, जो कि अभिसरण करता है जैसा . बहुस्तरीय पद्धति का आधार दूरबीन राशि पहचान है,[1]

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

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

अनुप्रयोग

सही

एमएलएमसी के पहले आवेदन का श्रेय माइक जाइल्स को दिया जाता है,[2] मोंटे कार्लो विकल्प मॉडल के लिए स्टोचैस्टिक अंतर समीकरण (एसडीई) के संदर्भ में, हालांकि, पैरामीट्रिक एकीकरण के संदर्भ में हेनरिक के काम में पहले के निशान पाए जाते हैं।[3] यहाँ, यादृच्छिक चर अदायगी समारोह, और सन्निकटन के अनुक्रम के रूप में जाना जाता है , प्रतिरूप पथ के सन्निकटन का उपयोग करें समय कदम के साथ .

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

एमएलएमसी अनुकरण के लिए एक एल्गोरिथ्म

एमएलएमसी सिमुलेशन के लिए एक सरल स्तर-अनुकूली एल्गोरिदम छद्म कोड में नीचे दिया गया है। दोहराना

    स्तर पर वार्म-अप के नमूने लें 

सभी स्तरों पर प्रतिरूप प्रसरण की गणना करें प्रतिरूपो की इष्टतम संख्या को परिभाषित करें सभी स्तरों पर प्रत्येक स्तर पर अतिरिक्त नमूने लें के अनुसार अगर तब

        अभिसरण के लिए परीक्षण
    अंत
    अगर नहीं मिला तो         

अंत

अभिसरण होने तक

एमएलएमसी का विस्तार

बहुस्तरीय मोंटे कार्लो पद्धति के हाल के विस्तार में मल्टी-इंडेक्स मोंटे कार्लो शामिल हैं,[6] जहां शोधन की एक से अधिक दिशाओं पर विचार किया जाता है, क्वासी-मोंटे कार्लो विधि पद्धति के साथ एमएलएमसी का संयोजन।[7][8]


यह भी देखें

संदर्भ

  1. Giles, M. B. (2015). "बहुस्तरीय मोंटे कार्लो तरीके". Acta Numerica. 24: 259–328. arXiv:1304.5472. doi:10.1017/s096249291500001x. S2CID 13805654.
  2. Giles, M. B. (2008). "बहुस्तरीय मोंटे कार्लो पथ सिमुलेशन". Operations Research. 56 (3): 607–617. CiteSeerX 10.1.1.121.713. doi:10.1287/opre.1070.0496. S2CID 3000492.
  3. Heinrich, S. (2001). "बहुस्तरीय मोंटे कार्लो तरीके". Lecture Notes in Computer Science (Multigrid Methods). Lecture Notes in Computer Science. Springer. 2179: 58–67. doi:10.1007/3-540-45346-6_5. ISBN 978-3-540-43043-8.
  4. Cliffe, A.; Giles, M. B.; Scheichl, R.; Teckentrup, A. (2011). "बहुस्तरीय मोंटे कार्लो के तरीके और रैंडम गुणांक वाले अण्डाकार पीडीई के अनुप्रयोग" (PDF). Computing and Visualization in Science. 14 (1): 3–15. doi:10.1007/s00791-011-0160-x. S2CID 1687254.
  5. Pisaroni, M.; Nobile, F. B.; Leyland, P. (2017). "कंप्रेसिबल इनविसिड एरोडायनामिक्स में अनिश्चितता मात्रा के लिए एक निरंतरता बहु स्तरीय मोंटे कार्लो विधि" (PDF). Computer Methods in Applied Mechanics and Engineering. 326 (C): 20–50. doi:10.1016/j.cma.2017.07.030. S2CID 10379943. Archived from the original (PDF) on 2018-02-14.
  6. Haji-Ali, A. L.; Nobile, F.; Tempone, R. (2016). "Multi-Index Monte Carlo: When Sparsity Meets Sampling". Numerische Mathematik. 132 (4): 767–806. arXiv:1405.3757. doi:10.1007/s00211-015-0734-5. S2CID 253742676.
  7. Giles, M. B.; Waterhouse, B. (2009). "बहुस्तरीय अर्ध-मोंटे कार्लो पथ अनुकरण" (PDF). Advanced Financial Modelling, Radon Series on Computational and Applied Mathematics. De Gruyter: 165–181.
  8. Robbe, P.; Nuyens, D.; Vandewalle, S. (2017). "लॉगनॉर्मल डिफ्यूजन प्रॉब्लम के लिए एक मल्टी-इंडेक्स क्वैसी-मोंटे कार्लो एल्गोरिथम". SIAM Journal on Scientific Computing. 39 (5): A1811–C392. arXiv:1608.03157. Bibcode:2017SJSC...39S.851R. doi:10.1137/16M1082561. S2CID 42818387.