Scribd Logo
Importer

Bienvenue sur Scribd !

  • 19 vues

    Détermination Des Trois Intégrales de Bessel - v1

    Transféré par

    VELONTSOA

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme DOCX, PDF, TXT ou lisez en ligne sur Scribd

    Détermination Des Trois Intégrales de Bessel - v1

    Transféré par

    VELONTSOA
    19 vues4 pages

    Titre original

    Détermination Des Trois Intégrales de Bessel_v1

    Copyright

    © © All Rights Reserved

    Formats disponibles

    DOCX, PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme DOCX, PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format docx, pdf ou txt
    19 vues4 pages

    Détermination Des Trois Intégrales de Bessel - v1

    Transféré par

    VELONTSOA

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme DOCX, PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format docx, pdf ou txt
    Vous êtes sur la page 1sur 4

    1

    Tentative de chercher une plus de précision aux quelques intégrales irrationnelles


    trigonométriques

    1. Position du problème

    On se propose de donner une plus de précision aux quelques intégrales numériques du type Bessel parmi
    lesquelles :

    (1)

    (2)

    (3)

    2. Intégrale connue

    On pense de pouvoir partir de la relation suivante

    (4)

    Avec :
    2

    3. Equations montées

    Actuellement, on a pu monter les huit équations suivantes :

    (5)

    (6)

    (7)

    (8)

    (9)

    (10)
    3

    (11)

    (12)

    La programmation et Les équations :

    clear all, close all


    clc
    %[X,Y] = meshgrid(-0.05:0.01:0.05,-0.05:0.01:0.05);
    X=.1;
    Y=.1;
    Z=.02;
    a1=.02;
    b1=.02;
    a=-2*b1*(X-a1);
    b=-2*Y.*b1;
    k=sqrt((X-a1).^2+Y.^2+Z.^2+b1^2);
    teta1=-pi/3;
    teta2=pi/3;
    p=1./(2*k.^2);
    %-------------les diverses données----------------
    c1prime=
    -2*(1./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^0.5-1./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^0
    .5);
    c11= -((a.*cos(teta2)+b.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^1.5-
    (a.*cos(teta1)+b.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^1.5);
    c12=1./k.^2*((a.*cos(teta2)+b.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^0.5-
    (a.*cos(teta1)+b.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^0.5);
    c1=c11+c12;
    c21=-((-a.*cos(teta2)+b.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^1.5-(-
    a.*cos(teta1)+b.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^1.5);
    c22=1./k.^2*((-a.*cos(teta2)+b.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^0.5-(-
    a.*cos(teta1)+b.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^0.5);
    c2=c21+c22;
    c31=-((b.*cos(teta2)+a.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^1.5-
    (b.*cos(teta1)+a.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^1.5);
    c32=1./k.^2*((b.*cos(teta2)+a.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^0.5-
    (b.*cos(teta1)+a.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^0.5);
    c3=c31+c32;
    c41=-((-b.*cos(teta2)+a.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^1.5-(-
    b.*cos(teta1)+a.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^1.5);
    c42=1./k.^2*((-b.*cos(teta2)+a.*sin(teta2))./(k.^2+a.*cos(teta2)+b.*sin(teta2)).^0.5-(-
    b.*cos(teta1)+a.*sin(teta1))./(k.^2+a.*cos(teta1)+b.*sin(teta1)).^0.5);
    c4=c41+c42;
    %-------------integrale de numérateur sinus carré---------------
    gsinuscarre=@(t)(sin(t)).^2./(k.^2+a.*cos(t)+b.*sin(t)).^1.5;
    V=quadv(gsinuscarre,teta1,teta2);
    gcosinuscarre=@(t)(cos(t)).^2./(k.^2+a.*cos(t)+b.*sin(t)).^1.5;
    U=quadv(gcosinuscarre,teta1,teta2);
    gsinuscos=@(t)(sin(t)).*cos(t)./(k.^2+a.*cos(t)+b.*sin(t)).^1.5;
    W=quadv(gsinuscos,teta1,teta2);
    %-----------------------------------------------------------------
    hsinuscarre=@(t)(sin(t)).^2./(k.^2+a.*cos(t)+b.*sin(t)).^2.5;
    Z=quadv(hsinuscarre,teta1,teta2);
    hcosinuscarre=@(t)(cos(t)).^2./(k.^2+a.*cos(t)+b.*sin(t)).^2.5;
    Y=quadv(hcosinuscarre,teta1,teta2);
    4
    hsinuscos=@(t)(sin(t)).*cos(t)./(k.^2+a.*cos(t)+b.*sin(t)).^2.5;
    X=quadv(hsinuscos,teta1,teta2);
    %-----------------les equations-------------------------------
    eq1=3/2*(-a.^2+b.^2)*X+3/2*a.*b.*Y-3/2*a.*b.*Z+p*(-a.^2+b.^2).*W+p*a.*b.*U-p*a.*b.*V
    c1=c11+c12
    eq2=3/2*(a.^2+b.^2)*X-3/2*a.*b.*Y-3/2*a.*b.*Z+p*(a.^2+b.^2).*W+3*p*a.*b.*U+3*p*a.*b.*V
    c2=c21+c22
    eq3=-3/2*a.^2*Z+3/2*b.^2*Y+p*(2*a.^2-b.^2).*U+p*(a.^2-2*b.^2).*V
    c3=c31+c32
    eq4=3*a.*b.*X-3/2*b.^2*Y-3/2*a.^2*Z+2*p*a.*b.*W+p*(2*a.^2+b.^2).*U+p*(a.^2+2*b.^2).*V
    c4=c41+c42
    eq5=3/2*(-a.^2+b.^2)*X+3/2*a.*b.*Y-3/2*a.*b.*Z
    c5=c1prime+c11
    eq6=p*(-a.^2+b.^2).*W+p*a.*b.*U-p*a.*b.*V
    c6=-c1prime+c12
    eq1 =1.0385, c1 =1.0385, eq2 =0.2020, c2 = 0.2020, eq3 = 1.0190,
    c3=1.0190, eq4 =-0.0266, c4 =-0.0266, eq5 =0.8263, c5 =0.8263,
    eq6 = 0.2122, c6 = 0.2122
    avec

    , , , ,

    , , , ,

    , , , ,

    , ,

    (13)

    (14)

    Vous aimerez peut-être aussi