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 :

Calcul intégrale avec scipy.integrate.quad


Sujet :

Calcul scientifique Python

  1. #1
    Débutant  
    Profil pro
    Inscrit en
    Novembre 2003
    Messages
    1 122
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Novembre 2003
    Messages : 1 122
    Points : 189
    Points
    189
    Par défaut Calcul intégrale avec scipy.integrate.quad
    Bonsoir

    j'utilise arrange de numpy pour trouver des valeurs de 0 a 32 par pas de 1.

    Dans l'exercice il est indiqué que l'on calcul l'intégrale de 0 à t, t contient les valeurs de np.arrange de 0 à 32.

    Voici ci-dessous mon petit programme.

    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
     
     
    from scipy import interpolate, integrate
    import numpy as np
     
    tx = [0, 0.7, 3.1, 5.7, 10.5, 13, 15.5, 18, 20.3, 21, 21.5, 22, 22.5, 22.9, 25.4, 27.8, 30.3, 32.8]
     
    vz = [0.0, 22.2, 39.9, 47.7, 52.2, 53.4, 53.8, 53.8, 53.8, 42.0, 33.6, 16.7, 10.0, 7.9, 7.4, 27.8, 7.4, 7.4]
     
    def fonct():
     
        f = interpolate.interp1d(tx, vz)
     
        t = np.arange(0, 32, 1)
     
        z = 1250 - integrate.quad(f, t,vz)
     
        return z
     
     
    def afficher():
     
     
        print('intégrale')
     
        print(fonct())
     
    afficher()
    Voici l'erreur que me renvoi python.

    Traceback (most recent call last):
    File "C:\Users\EDMOND\Documents\programme python\exam_ex1part2.py", line 36, in <module>
    afficher()
    File "C:\Users\EDMOND\Documents\programme python\exam_ex1part2.py", line 34, in afficher
    print(fonct())
    File "C:\Users\EDMOND\Documents\programme python\exam_ex1part2.py", line 24, in fonct
    z = 1250 - integrate.quad(f, t,vz)
    File "C:\Users\EDMOND\AppData\Local\Programs\Python\Python38\lib\site-packages\scipy\integrate\_quadpack_py.py", line 408, in quad
    flip, a, b = b < a, min(a, b), max(a, b)
    ValueError: operands could not be broadcast together with shapes (32,) (18,)
    Moi, je ne vois pas pourquoi l'exercice me de demande d'utiliser np.arrange, l'intégrale fait déjà le calcul sur un intervalle.

    A bientôt
    je suis un développeur debutant qui cherche à comprendre.

    Certain livre sont pas facile à comprendre.

  2. #2
    Expert éminent sénior
    Homme Profil pro
    Architecte technique retraité
    Inscrit en
    Juin 2008
    Messages
    21 287
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Manche (Basse Normandie)

    Informations professionnelles :
    Activité : Architecte technique retraité
    Secteur : Industrie

    Informations forums :
    Inscription : Juin 2008
    Messages : 21 287
    Points : 36 776
    Points
    36 776
    Par défaut
    Citation Envoyé par argon Voir le message
    Moi, je ne vois pas pourquoi l'exercice me de demande d'utiliser np.arrange, l'intégrale fait déjà le calcul sur un intervalle.
    C'est une question à poser à l'enseignant ou au site qui propose cet exercice.
    Pour l'erreur: quand l'appel à une fonction plante, on ouvre la documentation pour vérifier que les paramètres passés sont bien de la forme (et du type) de ceux qui sont attendus.

    - W
    Architectures post-modernes.
    Python sur DVP c'est aussi des FAQs, des cours et tutoriels

Discussions similaires

  1. [Python 3.X] intégrale simpson avec SCIPY
    Par argon dans le forum Calcul scientifique
    Réponses: 3
    Dernier message: 03/03/2023, 16h33
  2. Vitesse de calcul de scipy.integrate
    Par Kuchikhi dans le forum Calcul scientifique
    Réponses: 2
    Dernier message: 05/04/2014, 00h59
  3. Calcul d'intégrale avec bornes inconnues
    Par clem3900 dans le forum MATLAB
    Réponses: 1
    Dernier message: 19/07/2012, 12h15
  4. calcul d'une intégrale avec indicatrice
    Par Stalinou dans le forum MATLAB
    Réponses: 4
    Dernier message: 11/04/2011, 21h07
  5. Exécutable Java avec JRE intégré
    Par clawhammer dans le forum JBuilder
    Réponses: 2
    Dernier message: 06/10/2003, 16h26

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