Funció de Bessel

De testwiki
La revisió el 04:34, 1 març 2025 per imported>EVA3.0 (bot) (Tipografia)
(dif.) ← Versió més antiga | Versió actual (dif.) | Versió més nova → (dif.)
Salta a la navegació Salta a la cerca
La part radial dels modes de vibració d'un tambor circular segueixen la funció de Bessel.

Les funcions de Bessel són les solucions canòniques y(x) de l'equació diferencial de Bessel:[1]

x2d2ydx2+xdydx+(x2ν2)y=0

que tenen com a punt singular regular x=0 i una singularitat essencial a x=. El paràmetre ν és un nombre donat que es pot considerar positiu sense cap pèrdua de generalitat.

Aplicacions de la funció de Bessel

L'equació de Bessel sorgeix quan es busquen solucions separables a l'equació de Laplace i l'equació de Helmholtz en coordenades cilíndriques o esfèriques. Les funcions de Bessel són, per tant, especialment importants per molts problemes de propagació d'ones i potencials estàtics. En la resolució de problemes en sistemes de coordenades cilíndriques, s'obtenen funcions de Bessel d'ordre enter (Plantilla:Mvar = n); en problemes esfèrics, se n'obtenen d'ordre semi-enter (Plantilla:Mvar = n+1/2). Per exemple:

  • Ones electromagnètiques en una guia d'ones cilíndriques
  • Amplituds de pressió dels fluxos rotatacionals no viscosos
  • La conducció tèrmica en un objecte cilíndric
  • Modes de vibració d'una membrana acústica prima i circular (o annular) (com un tambor o un altre membranòfon)
  • Problemes de difusió en enreixats
  • Solucions de l'equació de Schrödinger radial (en coordenades esfèriques i cilíndriques) per una partícula lliure
  • Solució de patrons de radiació acústica
  • Fricció depenent de la freqüència en canoncades circulars
  • Dinàmica de cossos flotants
  • Resolució angular

Les funcions de Bessel també apareixen en altres problemes, com en el processament de senyal (per exemple, vegeu finestra de Kaiser o filtre de Bessel).

Definicions

Com que es tracta d'una equació diferencial de segon ordre, hi ha d'haver dues solucions linealment independents. Depenent de les circumstàncies, tanmateix, són convenients formulacions diverses d'aquestes solucions. En la següent taula, es resumeixen les diferents variacions, i es descriuen en les seccions següents.

Tipus Primer tipus Segon tipus
Funcions de Bessel Jα Yα
Funcions de Bessel modificades Iα Kα
Funcions de Hankel Hα(1) = Jα + iYα Hα(2) = Jα - iYα
Funcions de Bessel esfèriques jn yn
Funcions de Hankel esfèriques hn(1) = jn + iyn hn(2) = jn - iyn

Les funcions de Bessel de segon tipus i les funcions de Bessel esfèriques de segon tipus sovint s'anoten com a Nn i nn, respectivament, més que no tant Yn i yn.[2][3]

Funcions de Bessel de primer tipus: Plantilla:MvarPlantilla:Mvar

Les funcions de Bessel de primer tipus, denotades com a Jα(x), són solucions de l'equació diferencial de Bessel que són finites a l'origen (x = 0) per valors d'α enters o positius, i divergeixen a mesura que x tendeix a 0 per valors d' α no enters. Es pot definir la funció a través de la seva expansió en sèrie al voltant de x = 0, que es pot trobar aplicant el mètode de Frobenius a l'equació de Bessel:Plantilla:Sfn

Jα(x)=m=0(1)mm!Γ(m+α+1)(x2)2m+α

on Γ(z) és la funció gamma, una generalització canviada de la funció factorial per valors no enters. La funció de Bessel de primer tipus és una funció entera si Plantilla:Mvar és un enter, altrament és una funció multivaluada amb singularitat en el zero. La representació gràfica de les funcions de Bessel s'assemblen força funcions sinusoidals o cosinusoidals que decauen proporcionalment a 1/√x (vegin-se també les seves formes asimptòtiques més endavant), tot i que les seves arrels no són generalment periòdiques, excepte asimptòticament amb valors grans de x. (Les sèries mostren que −J1(x) és la derivada de J0(x), tal com −sin(x) és la derivada de cos(x); més generalment, la derivada de Jn(x) es pot expressar en termes de Jn±1(x) per les identitats que es mostren més endavant.)

Representació gràfica de la funció de Bessel de primer tipus, J Plantilla:Mvar(x), per ordres d'enters Plantilla:Mvar = 0, 1, 2

Per valors de Plantilla:Mvar no enters, les funcions J Plantilla:Mvar(x) i JPlantilla:Mvar(x) són linealment dependents, i són per tant les dues solucions de l'equació diferencial. D'altra banda, per Plantilla:Mvar d'ordre enter, la següent relació és vàlida (noti's que la funció gamma té pols simples per cadascun dels enters no positius):Plantilla:Sfn

Jn(x)=(1)nJn(x).

Això significa que les dues solucions ja no són linealment independents. En aquest cas, la segona de les solucions linealment independent a la primera es troba com la funció de Bessel de segon tipus, com es comentarà més endavant.

Integrals de Bessel

Una altra definició de la funció de Bessel, per valors enters de n, és possible mitjançant l'ús de la representació integral:[4]

Jn(x)=1π0πcos(nτxsin(τ))dτ.

Una altra representació integral és:[4]

Jn(x)=12πππei(nτxsin(τ))dτ.

Aquesta és l'aproximació que va fer Bessel, i a partir d'aquesta definció va derivar diverses propietats de la funció. La definició es pot estendre a ordres no enteres per una de les integrals de Schläfli, per (x)>0:[4]

Jα(x)=1π0πcos(ατxsinτ)dτsin(απ)π0exsinh(t)αtdt. Plantilla:Sfn[5][6][7]

Relació amb les sèries hipergeomètriques

Les funcions de Bessel es poden expressa en termes de les sèries hipergeomètriques generalitzades com:Plantilla:Sfn

Jα(x)=(x2)αΓ(α+1)0F1(α+1;x24).

Aquesta expressió està relacionada al desenvolupament de les funcions de Bessel en termes de la funció de Bessel-Clifford.

Relació amb els polinomis de Laguerre

En termes dels polinomis de Laguerre Lk i el paràmetre elegit arbitràriament t, la funció de Bessel es pot expressar com:[8]

Jα(x)(x2)α=etΓ(α+1)k=0Lk(α)(x24t)(k+αk)tkk!.

Funcions de Bessel de segon tipus: Plantilla:MvarPlantilla:Mvar

Les funcions de Bessel de segon tipus, anotades com a Yα(x), i ocasionalment anotades com a Nα(x), són solucions de l'equació diferencial de Bessel que té una singularitat a l'origen (x = 0) i són multivaluades. De vegades se les anomena funcions de Weber, ja que van ser introduïdes per H. M. Weber (1873), o també funcions de Neumann en honor de Carl Neumann.[9]

Representació gràfica de la funció de Bessel de segon tipus, Yα(x), per ordres enters Plantilla:Mvar = 0, 1, 2.

Per valors d'Plantilla:Mvar no enters, YPlantilla:Mvar(x) està relacionada amb JPlantilla:Mvar(x) com:

Yα(x)=Jα(x)cos(απ)Jα(x)sin(απ).

En el cas d'ordre enter Plantilla:Mvar, la funció es defineix prenent el límit quan Plantilla:Mvar no enter tendeix a Plantilla:Mvar,

Yn(x)=limαnYα(x).

Hi ha també una fórmula integral corresponent (per Re(Plantilla:Mvar) > 0),Plantilla:Sfn

Yn(x)=1π0πsin(xsinθnθ)dθ1π0[ent+(1)nent]exsinhtdt.

Yα(x) és necessàriament una segona solució linealment indepdendent de l'equació de Bessel quan α és un enter. Però Yα(x) pot adoptar altres significats més enllà d'aquest. Es pot considerar un company 'natural' de Jα(x). Vegeu també el subapartat de les funcions de Hankel més endavant.

Quan Plantilla:Mvar és un enter, a més, com ho era similarment en el cas de funcions del primer tipus, la següent relació és vàlida:

Yn(x)=(1)nYn(x).

Tant J Plantilla:Mvar(x) com Y Plantilla:Mvar(x) són funcions holomorfes de Plantilla:Mvar en el pla complex tallat al llarg de l'eix real negatiu. Quan Plantilla:Mvar és un enter, les funcions de Bessel Plantilla:Mvar són funcions enteres d'Plantilla:Mvar. Si Plantilla:Mvar es manté constant en un valor no zero, llavors les funcions de Bessel són funcions enteres d' Plantilla:Mvar.

Les funcions de Bessel del segon tipus quan Plantilla:Mvar és un enter és un exemple del segon tipus de solució en el teorema de Fuchs.

Funcions de Hankel:Hα(1), Hα(2)

Una altra formulació important de les dues solucions linealment independents a l'equació de Bessel són les funció Hankel de primer i segon tipus, Plantilla:Math i Plantilla:Math, definides com:Plantilla:Sfn[10]

Hα(1)(x)=Jα(x)+iYα(x),Hα(2)(x)=Jα(x)iYα(x),

on Plantilla:Mvar és la unitat imaginària. Aquestes combinacions lineals també són conegudes com a funcions de Bessel de tercer tipus; són dues solucions linealment independents de l'equació diferencial de Bessel. Duen el nom de Hermann Hankel.

La importància de les funcions de Hankel de primer i segon tipus rau més en llur desenvolupament teòric més que no en la seva aplicació. Aquestes formes de combinació lineal satisfan diverses propietats d'aspecte simple, com les fórmules asimptòtiques o les representacions integrals. Aquí, "simple" significa una aparença del factor de forma Plantilla:Math. Es pot entendre, doncs, que la funció de Bessel de segon tipus apareix de manera natural com la part imaginària de les funcions de Hankel.

Les funcions de Hankel de primer i segon tipus són usades per representar les solucions d'ones entrants i sortints d'una equació d'ones en simetries cilíndriques respectivament (o viceversa depenent de la convecció de signe de la freqüència).

Usant les relacions prèvies, es poden expressar com:

Hα(1)(x)=Jα(x)eαπiJα(x)isinαπ,Hα(2)(x)=Jα(x)eαπiJα(x)isinαπ.

Si Plantilla:Mvar és un enter, s'ha de calcular el límit the limit. Les següents relacions són vàlides, tant si Plantilla:Mvar és enter com si no ho és:Plantilla:Sfn

Hα(1)(x)=eαπiHα(1)(x),Hα(2)(x)=eαπiHα(2)(x).

En particular, si Plantilla:Math amb Plantilla:Mvar enter no negatiu, les relacions superiors impliquen directament que:

J(m+12)(x)=(1)m+1Ym+12(x),Y(m+12)(x)=(1)mJm+12(x).

Això és útil en el desenvolupament de les funcions esfèriques de Bessel (veure més a baix).

Les funcions de Hankel admeten les següents representacions integrals per Plantilla:Math:Plantilla:Sfn

Hα(1)(x)=1πi++iπexsinhtαtdt,Hα(2)(x)=1πi+iπexsinhtαtdt,

on els límits de les integrals indiquen que la integració és al llarg d'un contorn que es pot triar segons: de Plantilla:Math a 0 al llarg de l'eix real negatiu, de 0 a Plantilla:Math al llarg de l'eix imaginari, i de Plantilla:Math a Plantilla:Math al llarg d'un contorn paral·lel a l'eix real.Plantilla:Sfn

Les funcions de Bessel són vàlides fins i tot amb arguments complexos Plantilla:Mvar, i un cas especialment particular és quan l'argument és imaginari pur, és a dir, quan la part real de l'argument és zero. En aquest cas, les solucions de l'equació de Bessel reben el nom funcions de Bessel modificades (o també funcions hiperbòliques de Bessel) de primer i segon tipus i estan definides com:Plantilla:Sfn

Iα(x)=iαJα(ix)=m=01m!Γ(m+α+1)(x2)2m+α,Kα(x)=π2Iα(x)Iα(x)sinαπ,

quan Plantilla:Mvar no és un enter; quan Plantilla:Mvar és un enter, llavors s'utilitza el límit. Es tria que tinguin valors reals per arguments reals i positius Plantilla:Mvar. L'expansió en sèrie per Plantilla:Math és doncs similar a la de Plantilla:Math, però sense el factor alternant Plantilla:Math.

Si Plantilla:Math, Plantilla:Math poden ser expressades com a funcions de Hankel de primer tipus:

Kα(x)=π2iα+1Hα(1)(ix),

i si Plantilla:Math, es pot expressar com una funció de Hankel de segon tipus:

Kα(x)=π2(i)α+1Hα(2)(ix).

Es poden expressar les funcions primera i segona de Bessel en termes de les funcions modificades de Bessel (són vàlides si Plantilla:Math):Plantilla:Sfn

Jα(iz)=eαiπ2Iα(z),Yα(iz)=e(α+1)iπ2Iα(z)2πeαiπ2Kα(z).

Plantilla:Math i Plantilla:Math són les dues solucions linealment independents de lPlantilla:'equació modificada de Bessel:Plantilla:Sfn

x2d2ydx2+xdydx(x2+α2)y=0.

A diferència del cas de les funcions ordinàries de Bessel, que oscil·len com a funció d'un argument real, Plantilla:Mvar i Plantilla:Mvar són funcions exponencialment creixents i decreixents respectivament. Com en la funció de Bessel ordinària Plantilla:Mvar, la funció Plantilla:Mvar val zero a Plantilla:Math per Plantilla:Math i és finita a Plantilla:Math per Plantilla:Math. De forma anàloga, Plantilla:Mvar divergeix a Plantilla:Math amb una singularitat de tipus logarítmica.[11]

Funcions de Bessel modificades del primer tipus, Plantilla:Math, per Plantilla:Math
Funcions de Bessel modificades de segon tipus, Plantilla:Math, per Plantilla:Math


Dues fórmules integrals per les funcions de Bessel modificades són (per Plantilla:Math):Plantilla:Sfn

Iα(x)=1π0πexcosθcosαθdθsinαππ0excoshtαtdt,Kα(x)=0excoshtcoshαtdt.

En alguns càlculs de física, pot ser útil saber el valor de la següent relació:

2K0(ω)=eiωtt2+1dt.

Es pot demostrar mostrant la igualtat amb la definició integral de més amunt per Plantilla:Math. Això es fa integrant una corba tancada en el primer quadrant del pla complex.

Es poden representar les funcions de Bessel modificades Plantilla:Math i Plantilla:Math en termes de les integrals ràpidament convergents:[12]

K1/3(ξ)=30exp(ξ(1+4x23)1+x23)dx,K2/3(ξ)=1303+2x21+x23exp(ξ(1+4x23)1+x23)dx.

Les funcions de Bessel modificades de segon tipus han rebut també els següents noms (actualment en desús):

  • Funció de Basset per Alfred Barnard Basset
  • Funció de Bessel modificada de tercer tipus
  • Funció de Hankel modificada[13]
  • Funció de Macdonald per Hector Munro Macdonald
Funcions de Bessel esfèriques de primer tipus, Plantilla:Math, per Plantilla:Math
Funcions de Bessel esfèriques de segon tipus, Plantilla:Math, per Plantilla:Math

Quan se soluciona l'equació de Helmholtz en coordenades esfèriques amb separació de variables, l'equació radial té la forma:

x2d2ydx2+2xdydx+(x2n(n+1))y=0.

Les dues solucions linealment independents a aquest equació s'anomenen funcions de Bessel esfèriques Plantilla:Mvar i Plantilla:Mvar, i estan relacionades amb les funcions de Bessel ordinàries Plantilla:Mvar i Plantilla:Mvar segons:Plantilla:Sfn

jn(x)=π2xJn+12(x),yn(x)=π2xYn+12(x)=(1)n+1π2xJn12(x).

Plantilla:Mvar també es denota Plantilla:Mvar o Plantilla:Mvar; alguns autors anomenen aquestes funcions the funcions de Neumann esfèriques.

Les funcions de Bessel esfèriques també es poden escriure com a (fórmules de Rayleigh)Plantilla:Sfn

jn(x)=(x)n(1xddx)nsinxx,yn(x)=(x)n(1xddx)ncosxx.

La primera funció de Bessel esfèrica Plantilla:Math també es coneix com la funció sinc (no normalitzada). Les primeres funcions de Bessel esfèriques són:Plantilla:Sfn

j0(x)=sinxx.j1(x)=sinxx2cosxx,j2(x)=(3x21)sinxx3cosxx2,j3(x)=(15x36x)sinxx(15x21)cosxx

i:Plantilla:Sfn

y0(x)=j1(x)=cosxx,y1(x)=j2(x)=cosxx2sinxx,y2(x)=j3(x)=(3x2+1)cosxx3sinxx2,y3(x)=j4(x)=(15x3+6x)cosxx(15x21)sinxx.

Funcions generadores

Les funcions esfèriques de Bessel tenen les funcions generadores:Plantilla:Sfn

1zcos(z22zt)=n=0tnn!jn1(z),1zsin(z2+2zt)=n=0(t)nn!yn1(z).

Relacions diferencials

En la següent expressió, Plantilla:Mvar és qualsevol de Plantilla:Mvar, Plantilla:Mvar, Plantilla:Math, Plantilla:Math per Plantilla:MathPlantilla:Sfn

(1zddz)m(zn+1fn(z))=znm+1fnm(z),(1zddz)m(znfn(z))=(1)mznmfn+m(z).

Funcions de Hankel esfèriques:hn(1), hn(2)

Les funcions esfèriques de Hankel es defineixen de forma anàloga a les no esfèriques:

hn(1)(x)=jn(x)+iyn(x)
hn(2)(x)=jn(x)iyn(x).

De fet, això implica que existeixen expressions tancades de les funcions de Bessel d'ordre semienter en terme de funcions trigonomètriques i, per tant, també de les funcions esfèriques de Bessel. D'això es dedueix que, per n enter no negatiu es té:

hn(1)(x)=(i)n+1eixxm=0nimm!(2x)m(n+m)!(nm)!

i hn(2) és la funció complexa conjugada d'aquesta (per a x real). D'aquesta fórmula es poden deduir les formes tancades de les funcions esfèriques de Bessel ordinàries, per exemple, j0(x)=sin(x)/x i y0(x)=cos(x)/x, i en general per a qualsevol argument n.

Funcions de Riccati-Bessel: Sn, Cn, ξn, ζn

Les funcions de Riccati-Bessel són una petita modificació de les funció de Bessel esfèriques:

Sn(x)=xjn(x)=πx/2Jn+1/2(x)
Cn(x)=xyn(x)=πx/2Yn+1/2(x)
ξn(x)=xhn(1)(x)=πx/2Hn+1/2(1)(x)=Sn(x)iCn(x)
ζn(x)=xhn(2)(x)=πx/2Hn+1/2(2)(x)=Sn(x)+iCn(x)

Aquestes funcions satisfan la següent equació diferencial:

x2d2ydx2+[x2n(n+1)]y=0

Aquesta equació difeencial i les seves solucions, les equacions de Riccati-Bessel, s'usen per resoldre el problema de la dispersió d'ones electromagnètiques per una esfera, un problema conegut com a difusió de Mie, per la publicació per primer cop d'aquests resultat per Mie, l'any 1908. Vegeu, per exemple, Du (2004).[14]

Segons Debye (1909) s'utilitza de vegades la notació ψn,χn en lloc de Sn,Cn.

Funció de Bessel per a ν enter

El cas especial més habitual i important és aquell en què ν és un nombre enter; llavors ν es coneix com l'ordre de la funció de Bessel. Les funcions de Bessel de primer ordre són les solucions de l'equació diferencial de Bessel que són finites a l'origen per a enters ν no negatius i divergeixen quan x s'aproxima a zero per a ν negatius no enters. La solució general és:

y(x)=c1Jν(x)+c2Jν(x)

on J es defineix mitjançant la seva sèrie de Taylor al voltant de zero:

Jν(x)=n=0(1)nn!Γ(ν+n+1)(x2)ν+2n

i Γ és la funció gamma.

Referències

Plantilla:Commonscat Plantilla:Referències

Bibliografia

Plantilla:Autoritat

  1. Plantilla:Ref-llibre
  2. Plantilla:MathWorld
  3. Plantilla:MathWorld
  4. 4,0 4,1 4,2 Plantilla:Ref-llibre
  5. Plantilla:Ref-web
  6. Plantilla:Ref-web
  7. Arfken & Weber, exercise 11.1.17.
  8. Szegö, G. Orthogonal Polynomials, 4th ed. Providence, RI: Amer. Math. Soc., 1975.
  9. Plantilla:Ref-web
  10. Plantilla:Ref-llibre
  11. Plantilla:Ref-llibre
  12. Plantilla:Ref-publicació. Derived from formulas sourced to I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Fizmatgiz, Moscow, 1963; Academic Press, New York, 1980).
  13. Plantilla:Ref-publicació
  14. Plantilla:Ref-publicació