IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Calcul scientifique Python Discussion :

Evaluation d'une intégrale


Sujet :

Calcul scientifique Python

  1. #1
    Nouveau Candidat au Club
    Homme Profil pro
    Enseignant
    Inscrit en
    Novembre 2023
    Messages
    1
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Nord (Nord Pas de Calais)

    Informations professionnelles :
    Activité : Enseignant

    Informations forums :
    Inscription : Novembre 2023
    Messages : 1
    Points : 1
    Points
    1
    Par défaut Evaluation d'une intégrale
    Nom : Capture d'écran 2023-11-24 104019.png
Affichages : 113
Taille : 81,9 Ko

    Voici ce que j'ai fait mais je n'obtiens pas un résultat pou mon intégrale.
    Pourriez-vous m'aider, s'il vous plaît ?

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
     
    import pywt
    import numpy as np
    import matplotlib.pyplot as plt
    import scipy.integrate as integrate
     
    k=0.1
     
    # Paramètres
    wavelet_name = 'dmey'  # Ondelette de Meyer discrète
    scale = 10  # Échelle de l'ondelette de Meyer
     
    # Générer la fonction d'échelle
    wavelet = pywt.Wavelet(wavelet_name)
    phi, _, x = wavelet.wavefun(level=scale)
     
    # Calculer la transformée de Fourier de la fonction d'échelle
    fourier_transform_phi = np.fft.fftshift(phi)
    modsquare = np.abs(fourier_transform_phi) ** 2
    freq = np.fft.fftfreq(len(x), d=x[1]-x[0])
     
    def integrand(t,n):
        intermediaire = np.exp(1j * t * n)
        moninter = np.abs(np.sinc(t/ 2)) ** (2 * k) * modsquare
        return intermediaire * moninter
     
     
    n = 2
    # Intégration numérique de la fonction intégrande
    res = integrate.quad(integrand, -4 * np.pi / 3, 4 * np.pi / 3 , args=(n,))
     
     
    print(res)

  2. #2
    Rédacteur/Modérateur

    Avatar de User
    Homme Profil pro
    Développeur informatique
    Inscrit en
    Août 2004
    Messages
    8 199
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 54
    Localisation : France, Ain (Rhône Alpes)

    Informations professionnelles :
    Activité : Développeur informatique

    Informations forums :
    Inscription : Août 2004
    Messages : 8 199
    Points : 19 297
    Points
    19 297
    Billets dans le blog
    59
    Par défaut
    Bonjour,

    ça peut venir entre autre de la fonction np.fft.fftshift qui retourne un objet tableau (nbarray) :

    Code Python : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    ...
    fourier_transform_phi = np.fft.fftshift(phi)
    ...

    https://numpy.org/doc/stable/referen....fftshift.html

    Du coup la fonction à intégrer va aussi retourner un array alors qu'il attend une valeur unique.

    Après il faut voir aussi s'il sait gérer les nombres complexes..

    Cdlt
    Vous trouverez dans la FAQ, les sources ou les tutoriels, de l'information accessible au plus grand nombre, plein de bonnes choses à consulter sans modération

    Des tutoriels pour apprendre à créer des formulaires de planning dans vos applications Access :
    Gestion sur un planning des présences et des absences des employés
    Gestion des rendez-vous sur un calendrier mensuel


    Importer un fichier JSON dans une base de données Access :
    Import Fichier JSON

Discussions similaires

  1. Probleme d'evaluation sur une liste numérique
    Par cryptorchild dans le forum Langage
    Réponses: 2
    Dernier message: 16/01/2006, 20h39
  2. calcul approché d'une intégrale
    Par safaoua dans le forum Algorithmes et structures de données
    Réponses: 9
    Dernier message: 10/01/2006, 00h38
  3. evaluation d'une chaine de caractere (2+3)*4/5%8
    Par ricardvince dans le forum Algorithmes et structures de données
    Réponses: 24
    Dernier message: 23/11/2005, 18h02
  4. [EXP] Evaluation dans une expression régulière
    Par SergentHeinz dans le forum Langage
    Réponses: 7
    Dernier message: 10/11/2005, 19h17
  5. Evaluation d'une variable string ou char* en C++
    Par Coelacanthe dans le forum SL & STL
    Réponses: 2
    Dernier message: 08/11/2005, 10h27

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo