मोंटे कार्लो एकीकरण

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

मोंटे कार्लो एकीकरण करने के लिए अलग-अलग विधि हैं, जैसे समान वितरण (निरंतर), स्तरीकृत प्रतिचयन, महत्व प्रतिचयन , कण निस्यंदक (कण निस्यंदक के रूप में भी जाना जाता है), और माध्य-क्षेत्र कण विधियाँ।

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

समस्या मोंटे कार्लो एकीकरण पतों एक एकाधिक अभिन्न


 * $$I = \int_{\Omega}f(\overline{\mathbf{x}}) \, d\overline{\mathbf{x}}$$

की गणना है, जहाँ Ω, Rm का एक उपसमुच्चय है, आयतन


 * $$V = \int_{\Omega}d\overline{\mathbf{x}}$$ है

अनुभवहीन मोंटे कार्लो दृष्टिकोण Ω पर समान रूप से प्रतिदर्श बिंदुओं के लिए है: दिए गए N एकसमान प्रतिदर्श,


 * $$\overline{\mathbf{x}}_1, \cdots, \overline{\mathbf{x}}_N\in \Omega,$$

I को


 * $$ I \approx Q_N \equiv V \frac{1}{N} \sum_{i=1}^N f(\overline{\mathbf{x}}_i) = V \langle f\rangle$$ द्वारा अनुमानित किया जा सकता है।

ऐसा इसलिए है क्योंकि बड़ी संख्या का नियम यह सुनिश्चित करता है कि


 * $$ \lim_{N \to \infty} Q_N = I$$।

QN से I का अनुमान दिया गया है, QN की त्रुटि पट्टियों को भिन्नता


 * $$ \mathrm{Var}(f)\equiv\sigma_N^2 = \frac{1}{N-1} \sum_{i=1}^N \left (f(\overline{\mathbf{x}}_i) - \langle f \rangle \right )^2. $$

के निष्पक्ष अनुमान का उपयोग करके प्रतिदर्श भिन्नता से अनुमान लगाया जा सकता है। जो


 * $$ \mathrm{Var}(Q_N) = \frac{V^2}{N^2} \sum_{i=1}^N \mathrm{Var}(f) = V^2\frac{\mathrm{Var}(f)}{N} = V^2\frac{\sigma_N^2}{N}$$ की ओर जाता है।

जब तक क्रम


 * $$ \left \{ \sigma_1^2, \sigma_2^2, \sigma_3^2, \ldots \right \} $$

हुआ है, तब तक यह भिन्नता 1/N के रूप में शून्य से कम हो जाती है। QN की त्रुटि का अनुमान इस प्रकार


 * $$\delta Q_N\approx\sqrt{\mathrm{Var}(Q_N)}=V\frac{\sigma_N}{\sqrt{N}}$$

है, जो $$\tfrac{1}{\sqrt{N}}$$ के रूप में घटता है। यह $$V$$ से गुणा किए गए माध्य की मानक त्रुटि है। यह परिणाम अभिन्न के आयामों की संख्या पर निर्भर नहीं करता है, जो कि आयाम पर घातीय रूप से निर्भर करने वाले अधिकांश नियतात्मक विधियों के विरुद्ध मोंटे कार्लो एकीकरण का वचन दिया गया लाभ है। यह ध्यान रखना महत्वपूर्ण है कि, नियतात्मक विधियों के विपरीत, त्रुटि का अनुमान जटिल त्रुटि सीमा नहीं है; यादृच्छिक प्रतिचयन एकीकृत की सभी महत्वपूर्ण विशेषताओं को उजागर नहीं कर सकता है जिसके परिणामस्वरूप त्रुटि का अनुमान लगाया जा सकता है।

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

उदाहरण
thumb|right|स्केलिंग दिखाते हुए प्रतिदर्शों की संख्या के एक फलन के रूप में सापेक्ष त्रुटि $\tfrac{1}{\sqrt{N}}$ |link=|alt= {\displaystyle {\tfrac {1}{\sqrt {N}}}}

मोंटे कार्लो एकीकरण का एक आदर्श उदाहरण π का ​​अनुमान है। फलन


 * $$H\left(x,y\right)=\begin{cases}

1 & \text{if }x^{2}+y^{2}\leq1\\ 0 & \text{else} \end{cases}$$ और V = 4 के साथ समुच्चय Ω = [−1,1] × [−1,1] पर विचार करें। ध्यान दें कि


 * $$I_\pi = \int_\Omega H(x,y) dx dy = \pi.$$

इस प्रकार, मोंटे कार्लो एकीकरण के साथ π के मान की गणना करने की एक अपरिष्कृत विधि Ω पर N यादृच्छिक संख्या चुनना और


 * $$Q_N = 4 \frac{1}{N}\sum_{i=1}^N H(x_{i},y_{i})$$ की गणना करना है।

दाईं ओर की आकृति में, सापेक्ष त्रुटि $$\tfrac{Q_N-\pi}{\pi}$$ को N के एक फलन के रूप में मापा जाता है, जो $$\tfrac{1}{\sqrt{N}}$$ की पुष्टि करता है।

C उदाहरण
ध्यान रखें कि एक वास्तविक यादृच्छिक संख्या जनक का उपयोग किया जाना चाहिए।

पायथन उदाहरण
पायथन (प्रोग्रामिंग भाषा) में निर्मित।

वोल्फ्राम गणित उदाहरण
नीचे दिया गया कोड गणित में मोंटे-कार्लो पद्धति का उपयोग करके फलन
 * $$f(x) = \frac{1}{1+\sinh(2x)\log(x)^2}$$

को $$0.8<x<3$$ से एकीकृत करने की प्रक्रिया का वर्णन करता है:

पुनरावर्ती स्तरीकृत प्रतिचयन
thumb|right|पुनरावर्ती स्तरीकृत प्रतिचयन का एक उदाहरण। इस उदाहरण में, फलन: उपरोक्त उदाहरण से

$f(x,y) = \begin{cases}1 & x^2+y^2<1 \\0 & x^2+y^2 \ge 1 \end{cases}$ सुझाए गए एल्गोरिथम का उपयोग करके एक इकाई वर्ग के भीतर एकीकृत किया गया था। प्रतिचयित किए गए बिंदुओं को दर्ज किया गया और आलेखित किया गया। स्पष्ट रूप से स्तरीकृत प्रतिचयन एल्गोरिदम उन क्षेत्रों में बिंदुओं को केंद्रित करता है जहां फलन की भिन्नता सबसे बड़ी है।|link=|alt= {\displaystyle f(x,y)={\begin{cases}1&x^{2}+y^{2}<1\\0&x^{2}+y^{2}\geq 1\end{cases}}}

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

सामान्य 'दो से विभाजित' रणनीति बहु-आयामों के लिए काम नहीं करती है क्योंकि पद चिन्ह रखने के लिए उप-खंडों की संख्या बहुत तीव्रता से बढ़ती है। इसके बजाय एक अनुमान है कि किस आयाम के साथ एक उपखंड को सबसे अधिक लाभांश लाना चाहिए और मात्र इस आयाम के साथ मात्रा को उपविभाजित करता है।

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

लोकप्रिय मिसेर नेमी एक समान एल्गोरिथम लागू करता है।

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

स्तरीकृत प्रतिचयन का विचार अवलोकन के साथ प्रारम्भ होता है कि दो असंयुक्त समुच्चय क्षेत्रों के लिए a और b मोंटे कार्लो के साथ $$E_a(f)$$ और $$E_b(f)$$ और प्रसरण $$\sigma_a^2(f)$$ और $$\sigma_b^2(f)$$ का अनुमान है, संयुक्त अनुमान
 * $$E(f) = \tfrac{1}{2} \left (E_a(f) + E_b(f) \right )$$

का भिन्नता Var(f),


 * $$\mathrm{Var}(f) = \frac{\sigma_a^2(f)}{4 N_a} + \frac{\sigma_b^2(f)}{4 N_b}$$ द्वारा दिया गया है।

यह दिखाया जा सकता है कि इस भिन्नता को अंक वितरित करके कम किया जाता है जैसे कि,


 * $$\frac{N_a}{N_a + N_b} = \frac{\sigma_a}{\sigma_a + \sigma_b}$$

इसलिए प्रत्येक उप-क्षेत्र में फलन के मानक विचलन के अनुपात में प्रतिदर्श बिंदु आवंटित करके सबसे छोटी त्रुटि अनुमान प्राप्त किया जाता है।

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

महत्व प्रतिचयन
विभिन्न प्रकार के महत्वपूर्ण प्रतिचयन एल्गोरिदम हैं, जैसे

महत्व प्रतिदर्श एल्गोरिथ्म
महत्व प्रतिचयन मोंटे-कार्लो एकीकरण करने के लिए एक बहुत ही महत्वपूर्ण उपकरण प्रदान करता है। इस पद्धति के महत्व के प्रतिदर्श का मुख्य परिणाम यह है कि $$\overline{\mathbf{x}}$$ का एकसमान प्रतिचयन अधिक सामान्य विकल्प की एक विशेष स्थिति है, जिस पर किसी भी वितरण $$p(\overline{\mathbf{x}})$$ से प्रतिदर्श लिए जाते हैं। विचार यह है कि माप QN के विचरण को कम करने के लिए $$p(\overline{\mathbf{x}})$$ को चुना जा सकता है।

निम्नलिखित उदाहरण पर विचार करें जहां कोई गॉसियन फलन को संख्यात्मक रूप से एकीकृत करना चाहता है, जो 0 पर केंद्रित है, σ = 1 के साथ -1000 से 1000 तक। स्वाभाविक रूप से, यदि प्रतिदर्श अंतराल [-1000, 1000] पर समान रूप से खींचे जाते हैं, तो मात्र एक उनमें से बहुत छोटा भाग अभिन्न के लिए महत्वपूर्ण होगा। जहां से प्रतिदर्श चुने गए हैं वहां से एक अलग वितरण चुनकर इसमें सुधार किया जा सकता है, उदाहरण के लिए σ = 1 के साथ 0 पर केंद्रित गाऊसी वितरण के अनुसार प्रतिचयन करके। निश्चित रूप से उचित विकल्प एकीकृत पर दृढ़ता से निर्भर करता है।

विधिवत रूप से, वितरण


 * $$p(\overline{\mathbf{x}}) : \qquad \overline{\mathbf{x}}_1, \cdots, \overline{\mathbf{x}}_N \in V $$

से चुने गए प्रतिदर्शों का एक समुच्चय दिया गया है, I के लिए अनुमानक


 * $$ Q_N \equiv \frac{1}{N} \sum_{i=1}^N \frac{f(\overline{\mathbf{x}}_i)}{p(\overline{\mathbf{x}}_i)}$$ द्वारा दिया गया है

सहज रूप से, यह कहता है कि यदि हम किसी विशेष प्रतिदर्श को अन्य प्रतिदर्शों की तुलना में दोगुना चुनते हैं, तो हम इसे अन्य प्रतिदर्शों की तुलना में आधा वजन देते हैं। यह अनुमानक समान प्रतिचयन के लिए स्वाभाविक रूप से मान्य है, जहां मामला है $$p(\overline{\mathbf{x}})$$ स्थिर है।

मेट्रोपोलिस-हेस्टिंग्स एल्गोरिदम उत्पन्न करने के लिए सबसे अधिक उपयोग किए जाने वाले एल्गोरिदम में से एक है $$\overline{\mathbf{x}}$$ से $$p(\overline{\mathbf{x}})$$, इस प्रकार पूर्णांक कंप्यूटिंग का एक कुशल विधिप्रदान करता है।

वेगास मोंटे कार्लो
VEGAS एल्गोरिथ्म एकीकरण क्षेत्र पर कई पास बनाकर यथार्थ वितरण का अनुमान लगाता है जो फलन f का हिस्टोग्राम बनाता है। प्रत्येक हिस्टोग्राम का उपयोग अगले पास के लिए प्रतिचयन वितरण को परिभाषित करने के लिए किया जाता है। असम्बद्ध रूप से यह प्रक्रिया वांछित वितरण में परिवर्तित हो जाती है। K की तरह बढ़ने वाले हिस्टोग्राम डिब्बे की संख्या से बचने के लिएd, प्रायिकता बंटन एक वियोज्य फलन द्वारा अनुमानित है:


 * $$g(x_1, x_2, \ldots) = g_1(x_1) g_2(x_2) \ldots $$

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

यह भी देखें

 * क्वासी-मोंटे कार्लो विधि
 * सहायक क्षेत्र मोंटे कार्लो
 * सांख्यिकीय भौतिकी में मोंटे कार्लो विधि
 * मोंटे कार्लो विधि
 * भिन्नता में कमी

बाहरी संबंध

 * Café math : Monte Carlo Integration : A blog article describing Monte Carlo integration (principle, hypothesis, confidence interval)
 * Boost.Math : Naive Monte Carlo integration: Documentation for the C++ naive Monte-Carlo routines
 * Monte Carlo applet applied in statistical physics problems