Répondre à la discussion
Affichage des résultats 1 à 5 sur 5

Approximation, développement asymptotique - Fonction Bessel



  1. #1
    Lévesque

    Approximation, développement asymptotique - Fonction Bessel


    ------

    Bonjour,

    je révise mes notes de cours en Théorie des Perturbations, et je n'arrive pas à ce qui est écrit dans mes notes... C'est tellement simple comme calcul, que je ne vois pas mon erreur. C'est un exemple pour montrer que le développement asymptotique peut donner de meilleurs résultats que la série infinie.

    Voilà.

    La fonction peut se définir par une somme infinie:


    Son développement asymptotique est donné par

    ,

    pous z assez grand, avec et des séries divergentes.

    Mon prof illustre ensuite le comportement de l'approximation de la série si on la coupe avant l'infini. Cela revient à évaluer le reste défini comme la série infinie mois les termes qui manquent (parce qu'on coupe à N<infini):


    D'après ses résultats (s'ils sont bons), il trouve que pour z=15, le reste augmente avec N, jusqu'à N=10 environ, puis diminue jusqu'à une valeur pratiquement nulle à N=30 environ. Il lui reste à montrer qu'on fait une aussi bonne approximation avec moins de termes par le développement asymptotique.

    Or, lorsque je programme le reste en fonction de N, je trouve pour N=3, et assez petit pour que Mathematica indique 0. à N=6. Le reste que JE calcul diminue beaucoup plus vite que celui de mon prof, et beaucoup plus vite que celui du développement asymptotique.
    C'est pas normal tout ça...

    Quelqu'un voit un problème? Quelqu'un peu m'aider?

    Merci beaucoup à l'avance.

    Simon

    PS: voici mon expression du reste(N,Z) en language Mathematica, au cas ou j'aurais fait une erreur que je ne vois vraiment pas:
    In(1): Reste[n_, z_] := Sum[((-1)^k (z)^(2k))/(2^(2k)(k!)^2), {k, 0, Infinity}] - Sum[((-1)^k (z)^(2k))/(2^(2k)(k!)^2), {k, 0, n - 1}]

    In(2):Table[Reste[n,1/15],{n,0,10}]//N

    Out(2): (0.998889197492763, -0.001110802507237052, 3.0860387406998294*^-7, -3.810129760776037*^-11, 2.7617454674163618*^-15, 0.`, 0.`, 0.`,1.8099691943947416*^-16, 0.`, -1.6199032837520655*^-16)

    -----

  2. Publicité
  3. #2
    martini_bird

    Re : Approximation, développement asymptotique - Fonction Bessel

    Salut,

    comme mathematica conna&#238;t les fonctions de Bessel tu peux peut-&#234;tre essayer de remplacer
    Code:
    Sum[((-1)^k (z)^(2k))/(2^(2k)(k!)^2), {k, 0, Infinity}]
    par
    Code:
    BesselJ[0,z]
    Je ne vois pas trop, mais mathematica se choucroute peut-&#234;tre avec la s&#233;rie infinie?

    Cordialement.

  4. #3
    Lévesque

    Re : Approximation, développement asymptotique - Fonction Bessel

    Malheureusement, ce n'est pas ça. Regarde:
    Citation Envoyé par Mathematica
    In[10]:=
    Sum[((-1)^k (z)^(2k))/(2^(2k)(k!)^2),{k,0,Infinity}]

    Out[10]=
    BesselJ[0,z]
    Il la reconnait bien... Tu vois, c'est tellement simple comme calcul que ça devrait marcher. J'ai essayé de remplacer la somme par BesselJ, et j'obtient exactement (disons, par rapport à la précision de ma machine) les mêmes valeurs numériques:
    Citation Envoyé par Mathematica
    In[11]:=
    Reste1[n_,z_]:=BesselJ[0,z]-Sum[((-1)^k (z)^(2k))/(2^(2k)(k!)^2),{k,0,n-1}]

    In[12]:=
    Table[Reste1[n,1/15],{n,0,10}]//N

    Out[12]={0.998889197492763`,-0.001110802507237052`, 3.0860387401077105`*^-7, -3.81012998929009`*^-11, 2.6645352591003757`*^-15, 0.`, 0.`, 0.`,1.1102230246251565`*^-16, -1.1102230246251565`*^-16, -1.1102230246251565`*^-16}
    J'aimerais vraiment beaucoup trouver mon problème... Grrr!

    Merci mb!

    Simon

  5. #4
    Meumeul

    Re : Approximation, développement asymptotique - Fonction Bessel

    j'ai l'impression que le z=15 de ton cours s'est transform&#233; en z=1/15 dans ton calcul

  6. #5
    Lévesque

    Re : Approximation, développement asymptotique - Fonction Bessel

    Citation Envoyé par Meumeul
    j'ai l'impression que le z=15 de ton cours s'est transformé en z=1/15 dans ton calcul
    Youpi!!!!!!!!!

    J'ai commencé par répondre à ta question en croyant que c'était pas ça le problème. J'ai fait du code pour te le donner ici et te montrer que ce n'était pas ça. Je calculais jusqu'à N=10, et je trouvais que le reste augmente toujours jusqu'à N=10 pour atteindre 35000!!! Je te répondais alors que ça divergeais dans ce cas. En relisant, j'ai vu N=10, et j'ai décidé de bien vérifier avant d'envoyer. J'ai changer pour N=30 et le reste diminue!!!

    Jme doutais bien que c'était rien!

    Merci infiniment!

    Simon

  7. A voir en vidéo sur Futura

Sur le même thème :

Discussions similaires

  1. fonction de Bessel
    Par Ademma dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 05/07/2007, 10h05
  2. Suite/Série, Développement Asymptotique
    Par Ayrawhsia Aathsir Tia dans le forum Mathématiques du supérieur
    Réponses: 6
    Dernier message: 22/05/2007, 16h02
  3. developpement asymptotique et serie harmonique
    Par Ksilver dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 20/09/2006, 15h51
  4. Développement asymptotique de cotan
    Par doudache dans le forum Mathématiques du supérieur
    Réponses: 3
    Dernier message: 19/06/2006, 10h11
  5. Réponses: 0
    Dernier message: 28/01/2006, 11h34