Εκθετικό ολοκλήρωμα

Από testwiki
Μετάβαση στην πλοήγηση Πήδηση στην αναζήτηση
Plot of the exponential integral function E n(z) with n=2 in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
Γραφική παράσταση του εκθετικού ολοκληρώματος E n(z) με n=2 στο μιγαδικό επίπεδο από-2-2i έως 2+2i με χρώματα που δημιουργήθηκαν με τη συνάρτηση ComplexPlot3D του Mathematica 13.1 function ComplexPlot3D

Στα μαθηματικά, το εκθετικό ολοκλήρωμα Ei είναι μια ειδική συνάρτηση στο μιγαδικό επίπεδο.

Ορίζεται ως ένα συγκεκριμένο ορισμένο ολοκλήρωμα του λόγου μεταξύ μιας εκθετικής συνάρτησης και του ορίσματός της.

Ορισμός

Για πραγματικές μη μηδενικές τιμές του x, το εκθετικό ολοκλήρωμα Ei(x) ορίζεται ως εξής

Ei(x)=xettdt=xettdt.

Ο αλγόριθμος Ρις[1] δείχνει ότι η Ei δεν είναι στοιχειώδης συνάρτηση. Ο παραπάνω ορισμός μπορεί να χρησιμοποιηθεί για θετικές τιμές του x, αλλά το ολοκλήρωμα πρέπει να γίνει κατανοητό ως προς την κύρια τιμή Κωσύ λόγω της ιδιομορφίας του ολοκληρώματος στο μηδέν.

Για μιγαδικές τιμές του ορίσματος, ο καθορισμός γίνεται ασαφής λόγω των σημείων διακλάδωσης στο 0 και Πρότυπο:Nowrap [2] Αντί του Ei, χρησιμοποιείται ο ακόλουθος συμβολισμός,[3]

E1(z)=zettdt,|Arg(z)|<π
Plot of the exponential integral function Ei(z) in the complex plane from -2-2i to 2+2i with colors created with Mathematica 13.1 function ComplexPlot3D
Διάγραμμα του εκθετικού ολοκληρώματος Ei(z) στο μιγαδικό επίπεδο από -2-2i έως 2+2i με χρώματα που δημιουργήθηκε με το Mathematica 13.1 function ComplexPlot3D

Για θετικές τιμές του x, έχουμε Πρότυπο:Nowrap

Γενικά, η αποκοπή του κλάδου γίνεται στον αρνητικό πραγματικό άξονα και το E1 μπορεί να οριστεί με αναλυτική συνέχεια σε άλλο σημείο του μιγαδικού επιπέδου.

Για θετικές τιμές του πραγματικού μέρους του z, αυτό μπορεί να γραφεί[4]

E1(z)=1etztdt=01ez/uudu,(z)0.

Η συμπεριφορά του E1 κοντά στην αποκοπή της διακλάδωσης μπορεί να φανεί από την ακόλουθη σχέση:[5]

limδ0+E1(x±iδ)=Ei(x)iπ,x>0.

Ιδιότητες

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

Συγκλίνουσα σειρά

Γραφική παράσταση της συνάρτησης E1 (πάνω) και της συνάρτησης Ei (κάτω).

Για πραγματικά ή μιγαδικά ορίσματα εκτός του αρνητικού πραγματικού άξονα, το E1(z) μπορεί να εκφραστεί ως εξής [6]

E1(z)=γlnzk=1(z)kkk!(|Arg(z)|<π)

όπου γ είναι η σταθερά Όιλερ-Μαστσερόνι. Το άθροισμα συγκλίνει για όλα τα μιγαδικά z, και παίρνουμε τη συνήθη τιμή του μιγαδικού λογαρίθμου που έχει μια αποκοπή κλάδου κατά μήκος του αρνητικού πραγματικού άξονα.

Αυτός ο τύπος μπορεί να χρησιμοποιηθεί για τον υπολογισμό του E1(x) με πράξεις κινητής υποδιαστολής για πραγματικό x μεταξύ 0 και 2,5. Για x>2,5, το αποτέλεσμα είναι ανακριβές λόγω ακύρωσης.

Μια σειρά με ταχύτερη σύγκλιση βρέθηκε από τον Ραμανουτζάν:[7]

Ei(x)=γ+lnx+exp(x/2)n=1(1)n1xnn!2n1k=0(n1)/212k+1

Ασυμπτωτικές (αποκλίνουσες) σειρές

Σχετικό σφάλμα της ασυμπτωτικής προσέγγισης για διαφορετικό αριθμό N όρων στο αποκομμένο άθροισμα

Δυστυχώς, η σύγκλιση της παραπάνω σειράς είναι αργή για ορίσματα μεγαλύτερου συντελεστή. Επί παραδείγματι, απαιτούνται περισσότεροι από 40 όροι για να πάρουμε μια σωστή απάντηση με τρία σημαντικά ψηφία για το E1(10).[8] Ωστόσο, για θετικές τιμές του x, υπάρχει μια προσέγγιση αποκλίνουσας σειράς που μπορεί να ληφθεί με ολοκλήρωση xexE1(x) κατά μέρη:[9]

E1(x)=exp(x)x(n=0N1n!(x)n+O(N!xN))

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

Ασυμπτωτική πέραν όλων των τάξεων

Χρησιμοποιώντας ολοκλήρωση κατά μέρη, μπορούμε να λάβουμε έναν ρητό τύπο a[10]

Ei(z)=ezz(k=0nk!zk+en(z)),en(z)(n+1)! zezzettn+2dt Για οποιοδήποτε σταθερό z,

η απόλυτη τιμή του όρου σφάλματος |en(z)| μειώνεται και στη συνέχεια αυξάνεται. Το ελάχιστο εμφανίζεται στο n|z|,σε αυτό το σημείο

|en(z)|2π|z|e|z|. Αυτό το όριο λέγεται ότι είναι «ασυμπτωτικό πέραν όλων των τάξεων».

Εκθετική και λογαριθμική συμπεριφορά: παρένθεση

Παρένθεση της E1 με στοιχειώδεις συναρτήσεις

Από τις δύο σειρές που προτάθηκαν στις προηγούμενες υποενότητες, προκύπτει ότι η E1 συμπεριφέρεται σαν αρνητικός εκθετικός για μεγάλες τιμές του ορίσματος και σαν λογάριθμος για μικρές τιμές. Για θετικές πραγματικές τιμές του ορίσματος, η E1 μπορεί να τεθεί σε παρένθεση με στοιχειώδεις συναρτήσεις ως εξής:[11]

12exln(1+2x)<E1(x)<exln(1+1x)x>0

Η αριστερή πλευρά αυτής της ανισότητας παρουσιάζεται στο γράφημα στα αριστερά με μπλε χρώμα- το κεντρικό τμήμα E1(x) παρουσιάζεται με μαύρο χρώμα και η δεξιά πλευρά παρουσιάζεται με κόκκινο χρώμα.

Ορισμός από τον Έιν

Τόσο η Ei όσο και η E1 μπορούν να γραφούν πιο απλά χρησιμοποιώντας την ακέραια συνάρτηση Ein [12] defined as

Ein(z)=0z(1et)dtt=k=1(1)k+1zkkk!

(να σημειωθεί ότι πρόκειται απλώς για την εναλλασσόμενη σειρά στον παραπάνω ορισμό του E1). Τότε έχουμε

E1(z)=γlnz+Ein(z)|Arg(z)|<π
Ei(x)=γ+lnxEin(x)x0

Η συνάρτηση Ein σχετίζεται με την εκθετική συνάρτηση παραγωγής των αρμονικών αριθμών:

Ein(z)=ezn=1znn!Hn

Σχέση με άλλες συναρτήσεις

Εξίσωση του Κούμερ

zd2wdz2+(bz)dwdzaw=0

επιλύεται συνήθως με τις συγκλίνουσες υπεργεωμετρικές συναρτήσεις M(a,b,z) και U(a,b,z). Αλλά όταν a=0 και b=1, δηλαδή,

zd2wdz2+(1z)dwdz=0

έχουμε

M(0,1,z)=U(0,1,z)=1

για όλα τα z. Μια δεύτερη λύση δίνεται τότε από την E1(−z). Πράγματι,

E1(z)=γiπ+[U(a,1,z)M(a,1,z)]a,0<Arg(z)<2π

με την παράγωγο να εκτιμάται στο a=0. Μια άλλη σύνδεση με τις συρροές των υπεργεωμετρικών συναρτήσεων είναι ότι η E1 είναι εκθετική επί τη συνάρτηση U(1,1,z):

E1(z)=ezU(1,1,z)

Το εκθετικό ολοκλήρωμα συνδέεται στενά με τη λογαριθμική ολοκληρωτική συνάρτηση li(x) με τον τύπο

li(ex)=Ei(x)

για μη μηδενικές πραγματικές τιμές του x.

Γενίκευση

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

En(x)=1exttndt,

η οποία μπορεί να γραφεί ως ειδική περίπτωση της ανώτερης ατελούς συνάρτησης γάμμα:[13]

En(x)=xn1Γ(1n,x).

Η γενικευμένη μορφή ονομάζεται μερικές φορές συνάρτηση Μίσρα [14] φm(x), που ορίζεται ως

φm(x)=Em(x).

Πολλές ιδιότητες αυτής της γενικευμένης μορφής μπορούν να βρεθούν στη Ψηφιακή Βιβλιοθήκη Μαθηματικών Συναρτήσεων NIST

Συμπεριλαμβάνοντας έναν λογάριθμο ορίζεται η γενικευμένη ολοκληρο-εκθετική συνάρτηση[15]

Esj(z)=1Γ(j+1)1(logt)jezttsdt.

Παράγωγα

Οι παράγωγοι των γενικευμένων συναρτήσεων En μπορούν να υπολογιστούν μέσω του τύπου [16]

En(z)=En1(z)(n=1,2,3,)

Ας σημειωθεί ότι η συνάρτηση E0 είναι εύκολο να εκτιμηθεί (καθιστώντας αυτή την αναδρομή χρήσιμη), αφού είναι απλώς ez/z. [17]

Εκθετικό ολοκλήρωμα φανταστικού επιχειρήματος

E1(ix) κατά x- πραγματικό μέρος μαύρο, φανταστικό μέρος κόκκινο.

Αν το z είναι φανταστικό, έχει μη αρνητικό πραγματικό μέρος, οπότε μπορούμε να χρησιμοποιήσουμε τον τύπο

E1(z)=1etztdt

για να προκύψει μια σχέση με τα τριγωνομετρικά ολοκληρώματα

Si and Ci:

E1(ix)=i[12π+Si(x)]Ci(x)(x>0)

Το πραγματικό και το φανταστικό μέρος του E1(ix) απεικονίζονται στο σχήμα στα δεξιά με μαύρες και κόκκινες καμπύλες.

Προσεγγίσεις

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

  • Η προσέγγιση των Σουάμι και Οχίγια.[18] E1(x)=(A7.7+B)0.13,

όπου

A=ln[(0.56146x+0.65)(1+x)]B=x4e7.7x(2+x)3.7

  • Η προσέγγιση Άλεν και Χέιστινγκς [18][19]

E1(x)={lnx+aTx5,x1exxbTx3cTx3,x1

όπου

a[0.57722,0.99999,0.24991,0.05519,0.00976,0.00108]Tb[0.26777,8.63476,18.05902,8.57333]Tc[3.95850,21.09965,25.63296,9.57332]Txk[x0,x1,,xk]T

  • Η συνεχής επέκταση του κλάσματος [19]

E1(x)=exx+11+1x+21+2x+3.

  • Η προσέγγιση του Μπάρι και λοιποί[20]

E1(x)=exG+(1G)ex1Gln[1+Gx1G(h+bx)2],

όπου

h=11+xx+hq1+qq=2047x3126h=(1G)(G26G+12)3G(2G)2bb=2(1G)G(2G)G=eγ

με γ τη σταθερά Όιλερ-Μαστσερόνι.

Αντίστροφη συνάρτηση του εκθετικού ολοκληρώματος

Η αντίστροφη συνάρτηση του εκθετικού ολοκληρώματος μπορεί να εκφραστεί σε μορφή δυναμοσειράς:[21]

|x|<μln(μ),Ei1(x)=n=0xnn!Pn(ln(μ))μn

όπου μ είναι η σταθερά Ραμανουτζάν-Σόλντνερ και (Pn) είναι πολυωνυμική ακολουθία που ορίζεται από την ακόλουθη σχέση αναδρομής:

P0(x)=x, Pn+1(x)=x(Pn(x)nPn(x)).

Για n>0, degPn=n και έχουμε τον εξής τύπο :

Pn(x)=(ddt)n1(texEi(t+x)Ei(x))n|t=0.

Εφαρμογές

  • Χρονοεξαρτώμενη μεταφορά θερμότητας
  • Ροή εδαφικών υδάτων χωρίς ισορροπία στη λύση Θέις (γνωστή ως συνάρτηση πηγής)
  • Μεταφορά ακτινοβολίας σε αστρικές και πλανητικές ατμόσφαιρες
  • Εξίσωση ακτινικής διαχυτότητας για ροή σε μεταβατική ή ασταθή κατάσταση με πηγές και καταβόθρες γραμμής
  • Λύσεις της εξίσωσης μεταφοράς νετρονίων σε απλουστευμένες μονοδιάστατες γεωμετρίες [22]
  • Λύσεις της μη γραμμικής διαφορικής εξίσωσης Τραχένκο-Ζακόνε για την τεντωμένη εκθετική συνάρτηση στη χαλάρωση των άμορφων στερεών και τη υαλώδη μετάβαση[23][24].

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

Δείτε επίσης

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

Παραπομπές

Πρότυπο:Reflist

Πηγές

Πρότυπο:Κατάλογοι ολοκληρωμάτων Πρότυπο:Portal bar Πρότυπο:Authority control

  1. Πρότυπο:Cite web
  2. Abramowitz and Stegun, p. 228
  3. Abramowitz and Stegun, p. 228, 5.1.1
  4. Abramowitz and Stegun, p. 228, 5.1.4 with n = 1
  5. Abramowitz and Stegun, p. 228, 5.1.7
  6. Abramowitz and Stegun, p. 229, 5.1.11
  7. Andrews and Berndt, p. 130, 24.16
  8. Bleistein and Handelsman, p. 2
  9. Bleistein and Handelsman, p. 3
  10. Πρότυπο:Citation
  11. Abramowitz and Stegun, p. 229, 5.1.20
  12. Abramowitz and Stegun, p. 228, see footnote 3.
  13. Abramowitz and Stegun, p. 230, 5.1.45
  14. After Misra (1940), p. 178
  15. Milgram (1985)
  16. Abramowitz and Stegun, p. 230, 5.1.26
  17. Abramowitz and Stegun, p. 229, 5.1.24
  18. 18,0 18,1 Πρότυπο:Cite journal
  19. 19,0 19,1 Πρότυπο:Cite journal
  20. Πρότυπο:Cite journal
  21. Πρότυπο:Cite web
  22. Πρότυπο:Cite book
  23. Πρότυπο:Cite journal
  24. Πρότυπο:Cite journal