गॉसियन अभिन्न

आँकड़ों और भौतिकी से इस समाकलन को गौसियन चतुर्भुज, संख्यात्मक समाकलन की एक विधि के साथ भ्रमित नहीं होना है।

गॉसियन समाकलन, जिसे यूलर-पॉइसन समाकलन के रूप में भी जाना जाता है, गौसियन फलन $f(x)=e^{-x^{2}}$ का समाकलन है जो पूरी वास्तविक रेखा पर है। जर्मन गणितज्ञ कार्ल फ्रेडरिक गॉस के नाम पर, समाकलन है

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}.$$

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

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

$$\int e^{-x^{2}}\,dx,$$

लेकिन निश्चित समाकलन

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx$$

मूल्यांकन किया जा सकता है। एकपक्षीय गॉसियन फलन का निश्चित समाकलन है

$$\int _{-\infty }^{\infty }e^{-a(x+b)^{2}}\,dx={\sqrt {\frac {\pi }{a}}}.$$

संगणना

ध्रुवीय निर्देशांक द्वारा

गॉसियन समाकलन की गणना करने का एक मानक तरीका, जिसका विचार पोइसन तक जाता है, गुण का उपयोग करना है कि:

$$\left(\int _{-\infty }^{\infty }e^{-x^{2}}\,dx\right)^{2}=\int _{-\infty }^{\infty }e^{-x^{2}}\,dx\int _{-\infty }^{\infty }e^{-y^{2}}\,dy=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }e^{-\left(x^{2}+y^{2}\right)}\,dx\,dy.$$

फलन $e^{-\left(x^{2}+y^{2}\right)}=e^{-r^{2}}$तल $\mathbb {R} ^{2}$ पर विचार करें, और इसके समाकलन दो तरीकों की गणना करें:

1. एक ओर, कार्तीय समन्वय प्रणाली में दोहरे समाकलन द्वारा, इसका समाकलन वर्ग है:
   $$\left(\int e^{-x^{2}}\,dx\right)^{2};$$
2. दूसरी ओर, शेल समाकलन (ध्रुवीय निर्देशांक में दोहरे समाकलन की स्थिति) द्वारा, इसके समाकलन की गणना $\pi$ के रूप में की जाती है

इन दो संगणनाओं की तुलना करने से समाकलन प्राप्त होती है, हालांकि इसमें सम्मिलित अनुपयुक्त समाकलनो के बारे में ध्यान रखना चाहिए।

$${\begin{aligned}\iint _{\mathbb {R} ^{2}}e^{-\left(x^{2}+y^{2}\right)}dx\,dy&=\int _{0}^{2\pi }\int _{0}^{\infty }e^{-r^{2}}r\,dr\,d\theta \\[6pt]&=2\pi \int _{0}^{\infty }re^{-r^{2}}\,dr\\[6pt]&=2\pi \int _{-\infty }^{0}{\tfrac {1}{2}}e^{s}\,ds&&s=-r^{2}\\[6pt]&=\pi \int _{-\infty }^{0}e^{s}\,ds\\[6pt]&=\pi \left(e^{0}-e^{-\infty }\right)\\[6pt]&=\pi ,\end{aligned}}$$

जहां r का कारक जैकबियन निर्धारक है जो ध्रुवीय निर्देशांक में परिवर्तन के कारण प्रकट होता है (r dr dθ समतल पर मानक माप है, जिसे ध्रुवीय निर्देशांकों विकीबुक्स: गणना/ध्रुवीय समाकलन#सामान्यीकरण सामान्यीकरण में व्यक्त किया गया है, और प्रतिस्थापन में s = −r² इसलिए ds = −2r dr लेना सम्मिलित है।

इससे उत्पन्न का संयोजन

$$\left(\int _{-\infty }^{\infty }e^{-x^{2}}\,dx\right)^{2}=\pi ,$$

इसलिए

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}.$$

पूरा प्रमाण

अनुपयुक्त डबल समाकलन्स को सही ठहराने के लिए और दो भावों को समान करने के लिए, हम एक अनुमानित कार्य से शुरू करते हैं:

$$I(a)=\int _{-a}^{a}e^{-x^{2}}dx.$$

यदि समाकलन

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx$$

पूर्ण रूप से अभिसारी होते तो हमें उसका कौशी मूल मूल्य, यानी सीमा होती

$$\lim _{a\to \infty }I(a)$$

के साथ मेल खाएगा

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx.$$

यह देखने के लिए कि यह स्थिति है, उस पर विचार करें

$$\int _{-\infty }^{\infty }\left|e^{-x^{2}}\right|dx<\int _{-\infty }^{-1}-xe^{-x^{2}}\,dx+\int _{-1}^{1}e^{-x^{2}}\,dx+\int _{1}^{\infty }xe^{-x^{2}}\,dx<\infty .$$

तो हम गणना कर सकते हैं

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx$$

केवल सीमा लेकर

$$\lim _{a\to \infty }I(a).$$

का वर्ग ले रहा है $I(a)$ पैदावार

$${\begin{aligned}I(a)^{2}&=\left(\int _{-a}^{a}e^{-x^{2}}\,dx\right)\left(\int _{-a}^{a}e^{-y^{2}}\,dy\right)\\[6pt]&=\int _{-a}^{a}\left(\int _{-a}^{a}e^{-y^{2}}\,dy\right)\,e^{-x^{2}}\,dx\\[6pt]&=\int _{-a}^{a}\int _{-a}^{a}e^{-\left(x^{2}+y^{2}\right)}\,dy\,dx.\end{aligned}}$$

फ़ुबिनी के प्रमेय का उपयोग करते हुए, उपरोक्त दोहरे समाकलन को एक क्षेत्र समाकलन के रूप में देखा जा सकता है

$$\iint _{[-a,a]\times [-a,a]}e^{-\left(x^{2}+y^{2}\right)}\,d(x,y),$$

शीर्षों के साथ एक वर्ग पर कब्जा कर लिया {(−a, a), (a, a), (a, −a), (−a, −a)} कार्टेशियन विमान पर।

चूँकि सभी वास्तविक संख्याओं के लिए घातीय फलन 0 से अधिक है, तो यह इस प्रकार है कि वर्ग के अंतर्वृत्त पर लिया गया समाकलन इससे कम होना चाहिए $I(a)^{2}$, और इसी प्रकार वर्ग के परिवृत्त पर लिया गया समाकलन इससे बड़ा होना चाहिए $I(a)^{2}$. कार्टेसियन निर्देशांक से विहित समन्वय परिवर्तनों की सूची में स्विच करके दो डिस्क पर समाकलन आसानी से गणना की जा सकती है:

$${\begin{aligned}x&=r\cos \theta \\y&=r\sin \theta \end{aligned}}$$

$$\mathbf {J} (r,\theta )={\begin{bmatrix}{\dfrac {\partial x}{\partial r}}&{\dfrac {\partial x}{\partial \theta }}\\[1em]{\dfrac {\partial y}{\partial r}}&{\dfrac {\partial y}{\partial \theta }}\end{bmatrix}}={\begin{bmatrix}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \end{bmatrix}}$$

$$d(x,y)=|J(r,\theta )|d(r,\theta )=r\,d(r,\theta ).$$

$$\int _{0}^{2\pi }\int _{0}^{a}re^{-r^{2}}\,dr\,d\theta <I^{2}(a)<\int _{0}^{2\pi }\int _{0}^{a{\sqrt {2}}}re^{-r^{2}}\,dr\,d\theta .$$

(ध्रुवीय परिवर्तन में सहायता के लिए विहित समन्वय परिवर्तनों की सूची देखें।)

एकीकृत,

$$\pi \left(1-e^{-a^{2}}\right)<I^{2}(a)<\pi \left(1-e^{-2a^{2}}\right).$$

निचोड़ प्रमेय द्वारा, यह गॉसियन समाकलन देता है

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}.$$

कार्तीय निर्देशांक द्वारा

एक अलग तकनीक, जो लाप्लास (1812) तक जाती है,निम्नलखित में से कोई। होने देना

$${\begin{aligned}y&=xs\\dy&=x\,ds.\end{aligned}}$$

चूंकि सीमाएँ हैं s जैसा y → ±∞ के चिह्न पर निर्भर करता है x, यह इस तथ्य का उपयोग करने के लिए गणना को सरल करता है कि e^−x2 एक सम फलन है, और, इसलिए, सभी वास्तविक संख्याओं का समाकलन शून्य से अनंत तक समाकलन का दुगुना है। वह है,

$$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx=2\int _{0}^{\infty }e^{-x^{2}}\,dx.$$

इस प्रकार, समाकलन की सीमा से अधिक, x ≥ 0, और चर y और s की समान सीमाएँ हैं। यह प्रदान करता है:

$${\begin{aligned}I^{2}&=4\int _{0}^{\infty }\int _{0}^{\infty }e^{-\left(x^{2}+y^{2}\right)}dy\,dx\\[6pt]&=4\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-\left(x^{2}+y^{2}\right)}\,dy\right)\,dx\\[6pt]&=4\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}\left(1+s^{2}\right)}x\,ds\right)\,dx\\[6pt]\end{aligned}}$$

फिर, समाकलन के क्रम (कलन) को बदलने के लिए फ़ुबिनी के प्रमेय का उपयोग करना:

$${\begin{aligned}I^{2}&=4\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}\left(1+s^{2}\right)}x\,dx\right)\,ds\\[6pt]&=4\int _{0}^{\infty }\left[{\frac {e^{-x^{2}\left(1+s^{2}\right)}}{-2\left(1+s^{2}\right)}}\right]_{x=0}^{x=\infty }\,ds\\[6pt]&=4\left({\frac {1}{2}}\int _{0}^{\infty }{\frac {ds}{1+s^{2}}}\right)\\[6pt]&=2\arctan(s){\Big |}_{0}^{\infty }\\[6pt]&=\pi .\end{aligned}}$$

इसलिए, $I={\sqrt {\pi }}$, आशा के अनुसार।

लाप्लास की विधि से

लाप्लास सन्निकटन में, हम टेलर विस्तार में केवल दूसरे क्रम की शर्तों तक ही व्यवहार करते हैं, इसलिए हम विचार करते हैं $e^{-x^{2}}\approx 1-x^{2}\approx (1+x^{2})^{-1}$.

वास्तव में, चूंकि $(1+t)e^{-t}\leq 1$ सभी के लिए $t$, हमारे पास परिशुद्ध रूप से सीमाएँ हैं:

$$1-x^{2}\leq e^{-x^{2}}\leq (1+x^{2})^{-1}$$

तब हम लाप्लास सन्निकटन सीमा पर बाध्य कर सकते हैं:

$$\int _{[-1,1]}(1-x^{2})^{n}dx\leq \int _{[-1,1]}e^{-nx^{2}}dx\leq \int _{[-1,1]}(1+x^{2})^{-n}dx$$

वह है,

$$2{\sqrt {n}}\int _{[0,1]}(1-x^{2})^{n}dx\leq \int _{[-{\sqrt {n}},{\sqrt {n}}]}e^{-x^{2}}dx\leq 2{\sqrt {n}}\int _{[0,1]}(1+x^{2})^{-n}dx$$

त्रिकोणमितीय प्रतिस्थापन द्वारा, हम वास्तव में दो सीमाओं की गणना करते हैं: $2{\sqrt {n}}(2n)!!/(2n+1)!!$, $2{\sqrt {n}}(\pi /2)(2n-3)!!/(2n-2)!!$ वालिस सूत्र द्वारा, दो सीमाओं का भागफल 1 में परिवर्तित होता है। प्रत्यक्ष गणना द्वारा, दो सीमाओं का उत्पाद अभिसरण होता है $\pi$.

$${\frac {\pi }{2}}=\prod _{n=1}{\frac {(2n)^{2}}{(2n-1)(2n+1)}}$$

इसके विपरीत, यदि हम पहले उपरोक्त अन्य विधियों में से एक के साथ समाकलन की गणना करते हैं, तो हमें वालिस सूत्र का एक प्रमाण प्राप्त होगा।

गामा फलन से संबंध

इंटीग्रैंड एक समान कार्य है,

$$\int _{-\infty }^{\infty }e^{-x^{2}}dx=2\int _{0}^{\infty }e^{-x^{2}}dx$$

इस प्रकार, चर के परिवर्तन के बाद ${\textstyle x={\sqrt {t}}}$, यह यूलर समाकलन में बदल जाता है

$$2\int _{0}^{\infty }e^{-x^{2}}dx=2\int _{0}^{\infty }{\frac {1}{2}}\ e^{-t}\ t^{-{\frac {1}{2}}}dt=\Gamma \left({\frac {1}{2}}\right)={\sqrt {\pi }}$$

कहाँ ${\textstyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}dt}$ गामा फलन है। इससे पता चलता है कि आधे पूर्णांक का क्रमगुणन का परिमेय गुणक क्यों होता है ${\textstyle {\sqrt {\pi }}}$. आम तौर पर अधिक,

$$\int _{0}^{\infty }x^{n}e^{-ax^{b}}dx={\frac {\Gamma \left((n+1)/b\right)}{ba^{(n+1)/b}}},$$

जिसे प्रतिस्थापित करके प्राप्त किया जा सकता है $t=ax^{b}$ प्राप्त करने के लिए गामा फलन के समाकलन में ${\textstyle \Gamma (z)=a^{z}b\int _{0}^{\infty }x^{bz-1}e^{-ax^{b}}dx}$.

सामान्यीकरण

गौसियन फलन का समाकलन

Main article: Integral of a Gaussian function

एक स्वेच्छ गौसियन फलन का समाकलन है

$$\int _{-\infty }^{\infty }e^{-a(x+b)^{2}}\,dx={\sqrt {\frac {\pi }{a}}}.$$

एक वैकल्पिक रूप है

$$\int _{-\infty }^{\infty }e^{-ax^{2}+bx+c}\,dx={\sqrt {\frac {\pi }{a}}}\,e^{{\frac {b^{2}}{4a}}+c}.$$

यह फॉर्म सामान्य वितरण से संबंधित कुछ सतत संभाव्यता वितरण की अपेक्षाओं की गणना के लिए उपयोगी है, उदाहरण के लिए लॉग-सामान्य वितरण।

एन-आयामी और कार्यात्मक सामान्यीकरण

Main article: multivariate normal distribution

मान लीजिए A एक सममित सकारात्मक-निश्चित है (इसलिए उलटा) n × n परिशुद्ध मैट्रिक्स, जो परिशुद्ध रूप से मैट्रिक्स का मैट्रिक्स व्युत्क्रम है। तब,

$$\int _{-\infty }^{\infty }\exp {\left(-{\frac {1}{2}}\sum \limits _{i,j=1}^{n}A_{ij}x_{i}x_{j}\right)}\,d^{n}x=\int _{-\infty }^{\infty }\exp {\left(-{\frac {1}{2}}x^{\mathsf {T}}Ax\right)}\,d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det A}}}={\sqrt {\frac {1}{\det(A/2\pi )}}}={\sqrt {\det(2\pi A^{-1})}}$$

जहां समाकलन को R के ऊपर समझा जाता है^एन. यह तथ्य बहुभिन्नरूपी सामान्य वितरण के अध्ययन में लागू होता है।

भी,

$$\int x_{k_{1}}\cdots x_{k_{2N}}\,\exp {\left(-{\frac {1}{2}}\sum \limits _{i,j=1}^{n}A_{ij}x_{i}x_{j}\right)}\,d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det A}}}\,{\frac {1}{2^{N}N!}}\,\sum _{\sigma \in S_{2N}}(A^{-1})_{k_{\sigma (1)}k_{\sigma (2)}}\cdots (A^{-1})_{k_{\sigma (2N-1)}k_{\sigma (2N)}}$$

जहां σ का क्रमचय है {1, …, 2N} और दाहिनी ओर का अतिरिक्त गुणनखंड के सभी संयोजी युग्मों का योग है {{math|{1, …, 2N}}ए की एन प्रतियों की }^-1.

वैकल्पिक रूप से,

$$\int f({\vec {x}})\exp {\left(-{\frac {1}{2}}\sum _{i,j=1}^{n}A_{ij}x_{i}x_{j}\right)}d^{n}x={\sqrt {(2\pi )^{n} \over \det A}}\,\left.\exp {\left({1 \over 2}\sum _{i,j=1}^{n}\left(A^{-1}\right)_{ij}{\partial \over \partial x_{i}}{\partial \over \partial x_{j}}\right)}f({\vec {x}})\right|_{{\vec {x}}=0}$$

कुछ विश्लेषणात्मक कार्य f के लिए, बशर्ते कि यह इसके विकास और कुछ अन्य तकनीकी मानदंडों पर कुछ उपयुक्त सीमाओं को पूरा करे। (यह कुछ कार्यों के लिए काम करता है और दूसरों के लिए विफल रहता है। बहुपद ठीक हैं।) एक अंतर ऑपरेटर पर घातांक को एक शक्ति श्रृंखला के रूप में समझा जाता है।

जबकि कार्यात्मक समाकलन की कोई कठोर परिभाषा नहीं है (या यहां तक ​​​​कि ज्यादातर मामलों में एक गैर-कठोर कम्प्यूटेशनल), हम परिमित-आयामी मामले के अनुरूप एक गॉसियन कार्यात्मक समाकलन को परिभाषित कर सकते हैं। हालांकि, अभी भी समस्या है $(2\pi )^{\infty }$ अनंत है और साथ ही, कार्यात्मक निर्धारक भी सामान्य रूप से अनंत होगा। यदि हम केवल अनुपातों पर विचार करें तो इसका ध्यान रखा जा सकता है:

$$\int \exp \left(-{\frac {1}{2}}\sum _{i,j=1}^{n}A_{ij}x_{i}x_{j}+\sum _{i=1}^{n}B_{i}x_{i}\right)d^{n}x=\int e^{-{\frac {1}{2}}{\vec {x}}^{\mathsf {T}}\mathbf {A} {\vec {x}}+{\vec {B}}^{\mathsf {T}}{\vec {x}}}d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det {A}}}}e^{{\frac {1}{2}}{\vec {B}}^{\mathsf {T}}\mathbf {A} ^{-1}{\vec {B}}}.$$