Γκαουσιανό ολοκλήρωμα

Από testwiki
Μετάβαση στην πλοήγηση Πήδηση στην αναζήτηση
Η γραφική παράσταση της συνάρτησης f(x)=ex2 και το εμβαδόν μεταξύ αυτής και του άξονα x, (δηλ. ολόκληρη η πραγματική ευθεία) το οποίο είναι ίσο με π.

Το Γκαουσιανό ολοκλήρωμα, επίσης γνωστό ως ολοκλήρωμα Όιλερ-Πουασόν[1], είναι το ολοκλήρωμα της γκαουσιανής συνάρτησης f(x)=ex2 σε ολόκληρη την πραγματική γραμμή. Πήρε το όνομά του από τον Γερμανό μαθηματικό Καρλ Φρίντριχ Γκάους, το ολοκλήρωμα είναι

ex2dx=π.

Ο Αβραάμ ντε Μοίβρ ανακάλυψε αρχικά αυτόν τον τύπο ολοκληρώματος το 1733, ενώ ο Γκάους δημοσίευσε το ακριβές ολοκλήρωμα το 1809,[2] αποδίδοντας την ανακάλυψή του στον Λαπλάς. Το ολοκλήρωμα έχει ευρύ φάσμα εφαρμογών. Επί παραδείγματι, με μια μικρή αλλαγή των μεταβλητών χρησιμοποιείται για τον υπολογισμό της σταθεράς κανονικοποίησης της κανονικής κατανομής. Το ίδιο ολοκλήρωμα με πεπερασμένα όρια σχετίζεται στενά τόσο με τη συνάρτηση σφάλματος όσο και με τη συνάρτηση αθροιστικής κατανομής της κανονικής κατανομής. Στη φυσική αυτός ο τύπος ολοκληρώματος εμφανίζεται συχνά, όπως επί παραδείγματι στην κβαντομηχανική, για την εύρεση της πυκνότητας πιθανότητας της βασικής κατάστασης του αρμονικού ταλαντωτή. Το ολοκλήρωμα αυτό χρησιμοποιείται επίσης στη διατύπωση του ολοκληρώματος διαδρομής, για να βρεθεί ο πολλαπλασιαστής του αρμονικού ταλαντωτή, και στη στατιστική μηχανική, στη συνάρτηση κατανομής του.

Αν και δεν υπάρχει στοιχειώδης συνάρτηση για τη συνάρτηση σφάλματος, όπως μπορεί να αποδειχθεί με τον αλγόριθμο Ρις,[3] το γκαουσιανό ολοκλήρωμα μπορεί να επιλυθεί αναλυτικά μέσω των μεθόδων του πολυμεταβλητού λογισμού. Δηλαδή, δεν υπάρχει στοιχειώδες αόριστο ολοκλήρωμα για

ex2dx,

αλλά το ορισμένο ολοκλήρωμα

ex2dx

μπορεί να αξιολογηθεί. Το ορισμένο ολοκλήρωμα μιας αυθαίρετης γκαουσιανής συνάρτησης είναι

ea(x+b)2dx=πa.

Υπολογισμός

Με πολικές συντεταγμένες

Ένας συνηθισμένος τρόπος υπολογισμού του ολοκληρώματος Γκάους, η ιδέα του οποίου ανάγεται στον Πουασόν,[4] είναι να χρησιμοποιήσουμε την ιδιότητα ότι:

(ex2dx)2=ex2dxey2dy=e(x2+y2)dxdy.

Ας θεωρήσουμε τη συνάρτηση e(x2+y2)=er2 στο επίπεδο 2 και υπολογίζουμε το ολοκλήρωμά της με δύο τρόπους:

  1. αφενός, με διπλή ολοκλήρωση στο καρτεσιανό σύστημα συντεταγμένων, το ολοκλήρωμά του είναι τετράγωνο:

(ex2dx)2;

  1. αφετέρου, με ολοκλήρωση κελύφους (περίπτωση διπλής ολοκλήρωσης σε πολικές συντεταγμένες), το ολοκλήρωμά του υπολογίζεται ως εξής π

Συγκρίνοντας αυτούς τους δύο υπολογισμούς προκύπτει το ολοκλήρωμα, αν και θα πρέπει να προσέξουμε τα καταχρησικά ολοκληρώματα που εμπλέκονται.

2e(x2+y2)dxdy=02π0er2rdrdθ=2π0rer2dr=2π012esdss=r2=π0esds=limxπ(e0ex)=π,

όπου ο παράγοντας του Πρότυπο:Mvar είναι η Ιακωβιανή ορίζουσα που εμφανίζεται λόγω του μετασχηματισμού σε πολικές συντεταγμένες (Πρότυπο:Math είναι το τυπικό μέτρο στο επίπεδο, εκφρασμένο σε πολικές συντεταγμένες Γενίκευσης), και η αντικατάσταση περιλαμβάνει τη λήψη Πρότυπο:Math, οπότε Πρότυπο:Math.

Συνδυάζοντας αυτά τα αποτελέσματα (ex2dx)2=π,

οπότε

ex2dx=π.

Πλήρης απόδειξη

Για να δικαιολογήσουμε τα καταχρηστικά διπλά ολοκληρώματα και την εξίσωση των δύο εκφράσεων, ξεκινάμε με μια προσεγγιστική συνάρτηση:

I(a)=aaex2dx.

Εάν το ολοκλήρωμα

ex2dx

ήταν απολύτως συγκλίνουσα, θα είχαμε ότι η κύρια τιμή Κωσύ, δηλαδή το όριο

limaI(a)

συμπίπτει με

ex2dx.

Για να διαπιστώσετε ότι αυτό ισχύει, ας θεωρήσουμε ότι

|ex2|dx<1xex2dx+11ex2dx+1xex2dx<.

Οπότε μπορούμε να υπολογίσουμε

ex2dx

λαμβάνοντας απλώς το όριο

limaI(a).

Παίρνοντας το τετράγωνο του I(a) προκύπτει

I(a)2=(aaex2dx)(aaey2dy)=aa(aaey2dy)ex2dx=aaaae(x2+y2)dydx.

Χρησιμοποιώντας το θεώρημα του Φουμπίνι, το παραπάνω διπλό ολοκλήρωμα μπορεί να θεωρηθεί ως ολοκλήρωμα εμβαδού

[a,a]×[a,a]e(x2+y2)d(x,y),

που λαμβάνεται σε ένα τετράγωνο με κορυφές Πρότυπο:Math στο xy-επίπεδο.

Το Γκαουσιανό ολοκλήρωμα, επίσης γνωστό ως ολοκλήρωμα Όιλερ-Πουασόν, είναι το ολοκλήρωμα της γκαουσιανής συνάρτησης f(x)=ex2 σε ολόκληρη την πραγματική γραμμή. Το ολοκλήρωμα πήρε το όνομά του από τον Γερμανό μαθηματικό Καρλ Φρίντριχ Γκάους και είναι

x=rcosθy=rsinθ 𝐉(r,θ)=[xrxθyryθ]=[cosθrsinθsinθrcosθ] d(x,y)=|J(r,θ)|d(r,θ)=rd(r,θ). 02π0arer2drdθ<I2(a)<02π0a2rer2drdθ.

( Βλ. σε πολικές συντεταγμένες από καρτεσιανές συντεταγμένες για βοήθεια σχετικά με τον πολικό μετασχηματισμό.)

Ολοκλήρωση,

π(1ea2)<I2(a)<π(1e2a2).

Σύμφωνα με το Κριτήριο παρεμβολής, αυτό δίνει το γκαουσιανό ολοκλήρωμα

ex2dx=π.

Με καρτεσιανές συντεταγμένες

Μια διαφορετική τεχνική, η οποία ανάγεται στον Λαπλάς (1812)[4] , είναι η ακόλουθη. Έστω

y=xsdy=xds.

Δεδομένου ότι τα όρια του Πρότυπο:Mvar ως Πρότυπο:Math εξαρτώνται από το πρόσημο του Πρότυπο:Mvar, απλοποιείται ο υπολογισμός με τη χρήση του γεγονότος ότι το Πρότυπο:Math είναι μια άρτια συνάρτηση και, επομένως, το ολοκλήρωμα σε όλους τους πραγματικούς αριθμούς είναι απλώς το διπλάσιο του ολοκληρώματος από το μηδέν στο άπειρο. Δηλαδή,

ex2dx=20ex2dx.

Έτσι, στο εύρος ολοκλήρωσης, Πρότυπο:Math, και οι μεταβλητές Πρότυπο:Mvar και Πρότυπο:Mvar έχουν τα ίδια όρια. Από αυτό προκύπτει:

I2=400e(x2+y2)dydx=40(0e(x2+y2)dy)dx=40(0ex2(1+s2)xds)dx

Στη συνέχεια, χρησιμοποιώντας το θεώρημα του Φουμπίνι για να αλλάξουμε τη σειρά ολοκλήρωσης:

I2=40(0ex2(1+s2)xdx)ds=40[ex2(1+s2)2(1+s2)]x=0x=ds=4(120ds1+s2)=2arctan(s)|0=π.

Επομένως, I=π, όπως αναμενόταν.

Με τη μέθοδο Λαπλάς

Στην προσέγγιση Λαπλάς, ασχολούμαστε μόνο με όρους μέχρι δεύτερης τάξης στο ανάπτυγμα Τέιλορ, οπότε θεωρούμε ex21x2(1+x2)1.

Στην πραγματικότητα, αφού

(1+t)et1 για όλα τα t, έχουμε τα ακριβή όρια:1x2ex2(1+x2)1

Μπορούμε τότε να φέρουμε το όριο στο όριο της προσέγγισης Λαπλάς:

[1,1](1x2)ndx[1,1]enx2dx[1,1](1+x2)ndx

Δηλαδή,

2n[0,1](1x2)ndx[n,n]ex2dx2n[0,1](1+x2)ndx

Με τριγωνομετρική αντικατάσταση, υπολογίζουμε ακριβώς αυτά τα δύο όρια: 2n(2n)!!/(2n+1)!! και 2n(π/2)(2n3)!!/(2n2)!!

Λαμβάνοντας την τετραγωνική ρίζα του τύπου του Γουόλις,

π2=n=1(2n)2(2n1)(2n+1)

έχουμε π=2limnn(2n)!!(2n+1)!!, το επιθυμητό κατώτερο όριο. Ομοίως μπορούμε να πάρουμε το επιθυμητό άνω όριο. Αντίστροφα, αν υπολογίσουμε πρώτα το ολοκλήρωμα με μία από τις άλλες μεθόδους που προαναφέρθηκαν, θα λάβουμε μια απόδειξη του τύπου Γουόλις.

Σχέση με τη συνάρτηση γάμμα

Το ολοκλήρωμα είναι μια άρτια συνάρτηση,

ex2dx=20ex2dx

Έτσι, μετά την αλλαγή της μεταβλητής x=t, αυτό μετατρέπεται στο ολοκλήρωμα Όιλερ

20ex2dx=2012 et t12dt=Γ(12)=π

όπου Γ(z)=0tz1etdt είναι η συνάρτηση γάμμα. Αυτό δείχνει γιατί το παραγοντικό ενός μισού ακέραιου αριθμού είναι ένα ρητό πολλαπλάσιο του π.

0xneaxbdx=Γ((n+1)/b)ba(n+1)/b,

η οποία μπορεί να προκύψει αντικαθιστώντας t=axb στο ολοκλήρωμα της συνάρτησης γάμμα για να πάρουμε Γ(z)=azb0xbz1eaxbdx.

Γενικεύσεις

Ολοκλήρωμα μιας γκαουσιανής συνάρτησης

Το ολοκλήρωμα μιας αυθαίρετης γκαουσιανής συνάρτησης είναι

ea(x+b)2dx=πa.

Μια εναλλακτική μορφή είναι

e(ax2+bx+c)dx=πaeb24ac.

Αυτή η μορφή είναι χρήσιμη για τον υπολογισμό των προσδοκιών ορισμένων συνεχών κατανομών πιθανότητας που σχετίζονται με την κανονική κατανομή, όπως επί παραδείγματι η λογαριθμοκανονική κατανομή.

Μιγαδική μορφή

e12it2dt=eiπ/42π

και γενικότερα,

Ne12ixTAxdx=det(A)12(eiπ/42π)N

για κάθε θετικά ορισμένο συμμετρικό πίνακα A.

n-διάστατη και συναρτησιακή γενίκευση

Κύριο άρθρο: πολυμεταβλητή κανονική κατανομή

Ας υποθέσουμε ότι ο Α είναι ένας συμμετρικός θετικά ορισμένος (άρα αντιστρέψιμος) πίνακας ακριβείας Πρότυπο:Math ο οποίος είναι ο αντίστροφος πίνακας του πίνακα συνδιακύμανσης. Τότε,

nexp(12i,j=1nAijxixj)dnx=nexp(12x𝖳Ax)dnx=(2π)ndetA=1det(A/2π)=det(2πA1)

Ολοκληρώνοντας το τετράγωνο, αυτό γενικεύεται σε

nexp(12x𝖳Ax+b𝖳x+c)dnx=det(2πA1)e12b𝖳A1b+c

Το γεγονός αυτό εφαρμόζεται στη μελέτη της πολυμεταβλητής κανονικής κατανομής.

Επίσης,

xk1xk2Nexp(12i,j=1nAijxixj)dnx=(2π)ndetA12NN!σS2N(A1)kσ(1)kσ(2)(A1)kσ(2N1)kσ(2N)

όπου σ είναι μια μετάθεση του Πρότυπο:Math και ο επιπλέον παράγοντας στη δεξιά πλευρά είναι το άθροισμα πάνω σε όλα τα συνδυαστικά ζεύγη Πρότυπο:Math των N αντιγράφων του A-1.

Εναλλακτικά,[5]

f(x)exp(12i,j=1nAijxixj)dnx=(2π)ndetAexp(12i,j=1n(A1)ijxixj)f(x)|x=0

για κάποια αναλυτική συνάρτηση f, υπό την προϋπόθεση ότι ικανοποιεί κάποια κατάλληλα όρια για την ανάπτυξή της και κάποια άλλα τεχνικά κριτήρια. (Λειτουργεί για ορισμένες συναρτήσεις και αποτυγχάνει για άλλες. Τα πολυώνυμα είναι μια χαρά.) Το εκθετικό επί ενός διαφορικού τελεστή νοείται ως σειρά δυνάμεων.

Ενώ τα συναρτησιακά ολοκληρώματα δεν έχουν αυστηρό ορισμό (ούτε καν έναν μη αυστηρό υπολογιστικό στις περισσότερες περιπτώσεις), μπορούμε να ορίσουμε ένα γκαουσιανό συναρτησιακό ολοκλήρωμα κατ' αναλογία με την περίπτωση των πεπερασμένων διαστάσεων. Ωστόσο, εξακολουθεί να υπάρχει το πρόβλημα ότι (2π) είναι άπειρο και επίσης, ο συναρτησιακός προσδιοριστής θα ήταν επίσης άπειρος γενικά. Αυτό μπορεί να τακτοποιηθεί αν εξετάσουμε μόνο τις αναλογίες:

f(x1)f(x2N)exp[12A(x2N+1,x2N+2)f(x2N+1)f(x2N+2)ddx2N+1ddx2N+2]𝒟fexp[12A(x2N+1,x2N+2)f(x2N+1)f(x2N+2)ddx2N+1ddx2N+2]𝒟f=12NN!σS2NA1(xσ(1),xσ(2))A1(xσ(2N1),xσ(2N)).

Στον συμβολισμό DeWitt, η εξίσωση μοιάζει πανομοιότυπη με την περίπτωση των πεπερασμένων διαστάσεων.

n'-διάσταση με γραμμικό όρο

Αν ο A είναι πάλι ένας συμμετρικός θετικά ορισμένος πίνακας, τότε (υποθέτοντας ότι όλα είναι διανύσματα στήλης)

exp(12i,j=1nAijxixj+i=1nBixi)dnx=e12x𝖳𝐀x+B𝖳xdnx=(2π)ndetAe12B𝖳𝐀1B.

Ομοειδή ολοκληρώματα

0x2nex2a2dx=πa2n+1(2n1)!!2n+1 0x2n+1ex2a2dx=n!2a2n+2 0x2nebx2dx=(2n1)!!bn2n+1πb 0x2n+1ebx2dx=n!2bn+1 0xnebx2dx=Γ(n+12)2bn+12

όπου n είναι ένας θετικός ακέραιος αριθμός

Ένας εύκολος τρόπος για να τα εξάγουμε αυτά είναι να τα διαφοροποιήσουμε κάτω από το πρόσημο του ολοκληρώματος

x2neαx2dx=(1)nnαneαx2dx=(1)nnαneαx2dx=π(1)nnαnα12=πα(2n1)!!(2α)n

Θα μπορούσε κανείς επίσης να ολοκληρώσει κατά μέρη και να βρει μια αναδρομική σχέση για την επίλυση αυτού του προβλήματος.

Πολυώνυμα ανώτερης τάξης

Η εφαρμογή μιας γραμμικής αλλαγής της βάσης δείχνει ότι το ολοκλήρωμα του εκθετικού ενός ομογενούς πολυωνύμου σε n μεταβλητές μπορεί να εξαρτάται μόνο από τις SL(n)-παραλλαγές του πολυωνύμου. Μια τέτοια αναλλοίωτη είναι η διακριτική ικανότητα, τα μηδενικά της οποίας σηματοδοτούν τις ιδιομορφίες του ολοκληρώματος. Ωστόσο, το ολοκλήρωμα μπορεί επίσης να εξαρτάται από άλλες αναλλοίωτες[6].

Τα εκθετικά άλλων ζυγών πολυωνύμων μπορούν να επιλυθούν αριθμητικά χρησιμοποιώντας σειρές. Αυτές μπορούν να ερμηνευθούν ως τυπικοί υπολογισμοί όταν δεν υπάρχει σύγκλιση. Επί παραδείγματι, η λύση του ολοκληρώματος του εκθετικού ενός τεταρτοβάθμιου πολυωνύμου είναι

eax4+bx3+cx2+dx+fdx=12efn,m,p=0n+p=0mod2bnn!cmm!dpp!Γ(3n+2m+p+14)(a)3n+2m+p+14.

Η απαίτηση Πρότυπο:Math mod 2 οφείλεται στο γεγονός ότι το ολοκλήρωμα από το -∞ έως το 0 συνεισφέρει έναν παράγοντα Πρότυπο:Math σε κάθε όρο, ενώ το ολοκλήρωμα από το 0 έως το +∞ συνεισφέρει έναν παράγοντα 1/2 σε κάθε όρο. Αυτά τα ολοκληρώματα εμφανίζονται σε θέματα όπως η κβαντική θεωρία πεδίου.

Εξωτερικοί σύνδεσμοι

Δείτε επίσης

Βιβλιογραφία

Παραπομπές

Πρότυπο:Reflist

Πηγές

Πρότυπο:Portal bar Πρότυπο:Authority control