[PDF] Méthodes numériques et langage Python Intégration numérique





Previous PDF Next PDF



Analyse numérique en Python Intégration et dérivation

Ce chapitre est le premier d'une série s'intéressant à la question du calcul numérique. Nous verrons dans ces chapitres comment Python peut être utilisé 



Méthodes numériques et langage Python Intégration numérique

8 oct. 2019 ? On appelle une formule de quadrature une expression linéaire fournissant une intégration approchée sur un intervalle. Raisons. ? f n'est ...



Intégration numérique déquations différentielles

On veut assurer que la solution numérique est stable dans le sens que l'erreur ne diverge pas. Appelons donc. 1. yi la solution exacte au point xi.



Analyse numérique avec Python

22 mai 2014 2 Intégration numérique. Dans cette section nous nous intéresserons aux algorithmes permettant le calcul numérique d'intégrales. Il s'agit bien ...



Informatique en PCSI et MPSI Champollion 2013-2014 Méthodes d

23 avr. 2014 Champollion 2013-2014. Méthodes d'Analyse Numériques. Implémentation et Application en Python. Intégration Numérique. A. HASSAN.



Quelques méthodes dintégration numérique

15 avr. 2018 En ordonnée l'erreur commise. On choisit d'intégrer l'exponentielle entre 0 et 1



Intégration numérique déquations différentielles

3 Principe d'intégration numérique. Méthode d'Euler. Intégrateur de type Leapfrog. Stabilité. Cl`audia Payrató A. Torcini et A. Honecker.



Chapitre 5 - Méthodes dintégration numérique

Méthodes d'intégration numérique. Le but. Le but de ce chapitre est d'aborder le calcul général de l'intégrale d'une fonction f(x) sur un domaine.



Analyse numérique en Python Résolution numérique déquations

1. Ou bien définie par morceaux avec de telles fonctions. 2.2 Schéma d'intégration d'Euler explicite. Un schéma d 



Python MP PC

TSI Oral



Donner son temps : les bénévoles dans la vie - INSEE

2 Idées de base de l’intégration numérique Nousconsidéronsl’intégrale Z b a f(x)dx (1) La plupart des méthodes numériques de calcul de cette intégrale divisent l’intégraled’origineenunesommedeplusieursintégraleschacunecouvrantune partiepluspetitedel’intervalled’intégrationd’origine[ab] Cetteréécriturede



Quelques méthodes d'intégration numérique

Le but de ce notebook est de comparer quelques méthodes d'intégration numérique Nous nous intéresserons à trois métodes : les méthodes des rectangles et des trapèzes vues en



Méthodes et simulation numériques - langage Python IPSA

1 Chapitre 1 : Intégration numérique Méthode rectangle Méthode du point milieu Méthode de trapèze Méthode de Simpson Mise en application en TP par Python 2 Chapitre 2 : Introduction à la résolution d’EDO Méthode d’Euler Méthode de Range-Kutta Mise en application en TP par Python 3 Chapitre 3 : Simulation de variables aléatoires



Searches related to python intégration numérique PDF

Intégration et dérivation Ce chapitre est le premier d’une série s’intéressant à la question du calcul numérique Nous verrons dans ces chapitres comment Python peut être utilisé comme un outil aidant à résoudre des problèmes de mathématiques de physique de sciences de l’ingénieur ou bien encore de chimie

Comment fonctionne l’intégration numérique?

D’un point de vue empirique, 20 points de support sont utilisés pour l’intégration numérique et les écarts-types sont obtenus à partir de l’estimation numérique de l’inverse de l’opposée de la matrice hessienne.

Quels sont les principes de dérivation et d’intégration numériques?

Introduction Dérivation et intégration numériques Déterminer avec précision : 1. La vitesse à chaque instant 2. L’accélération de la fusée 3. La consommation de carburant Évaluer les dérivées premières et secondes ainsi que l’intégrale de cette fonction. Ift2421 3 Chapitre 5 Principe général de dérivation et d’intégration numériques

Comment faire une résolution numérique d'équations differentielles ordinaires en Python ?

Le but de cette page est présenter quelques applications possibles en cours de physique de la résolution numérique d'équations differentielles ordinaires en python. Pour cela, nous allons utiliser la fonction odeint du module scipy.

Comment intégrer le numérique dans les pratiques pédagogiques?

Des thématiques de formation pour intégrer le numérique dans les pratiques pédagogiques Equiper chaque salle de classe, chaque enseignant et chaque élève avec du matériel neuf et performant. Des équipements performants pour favoriser l’usage du numérique en classe Mettre à disposition une sélection de ressources et d’outils numériques pédagogiques.

M´ethodes num´eriques et langage Python

Int´egration num´erique

R. Flamary

8 octobre 2019Int´egration num´erique

I=? b a f(x)dx=F(b)-F(a)abxy f(x) S

Objectif

Calcul num´erique de la surfaces`a partir d"un nombre fini d"appels `a la fonction. Cas particulier lorsque l"on a acc`es `a un ´echantillonnage r´egulier. On appelle une formule de quadrature une expression lin´eaire fournissant une int´egration approch´ee sur un intervalle.Raisons fn"est connue qu"en certains points. PrimitiveFconnue mais pas une fonction ´el´ementaire. (int´egrale deexp(-x2)). Primitive trop difficile `a calculer num´eriquement.2/16D´efinitions I=? b a f(x)dxHypoth`eses aetbdeux r´eels aveca < b.

On supposefint´egrable sur[a,b].

fne dispose pas de singularit´es sur[a,b].Mise en oeuvre

Approche classique :

D´ecomposition du domaine en morceaux (un intervalle en sous-intervalles contigus). Int´egration approch´ee de la fonction sur chaque morceau. Sommation des r´esultats num´eriques ainsi obtenus. Utilisation de polynˆomes pour approcher la fonction sur chaque morceau.

3/16abxy

f(x) S

Int´egration par quadrature simple (2)

Formule de Simpson

Popularis´ee pas Simpson mais utilis´ee par Kepler 100 ans plus tˆot.

Interpolation defpar un polynˆome d"ordre 2.

On calcule la valeur de la fonction ena,betm=a+b2

Interpolation polynomiale de Lagrange :

P(x) =f(a)(x-m)(x-b)(a-m)(a-b)+f(m)(x-a)(x-b)(m-a)(m-b)+f(b)(x-a)(x-m)(b-a)(b-m). Int´egrale approch´ee obtenue en int´egrant le polynome :

I=hf(a) + 4f?a+b2

?+f(b)6 +O?h590?25f(4)(μ)?

M´ethode d"ordre 2, exacte pour des fonctions

quadratiques et cubiques.Quadrature de Gauss

On approche encore une fois l"int´egrale par

I=? b a f(x)dx≈n? i=1w if(xi) Leswisont les coefficients de quadrature et lesxisont choisis comme les racine de polynˆomes orthogonaux. On n"utilise pas d"´echantillonnage r´egulier, possibilit´e d"avoir une meilleure approximation def. Int´egration exacte pour des polynˆomes de degr´e2n-1.

Exemple pour[a,b] = [-1,1]avec les polynˆomes de Legendre :Nb de pointsPoidswiPointsxiPoly. de Legendre

120x
21,1-

1⎷3

,1⎷3(3x2-1)/235 9 ,89 ,59- ⎷3⎷5 ,0,⎷3⎷5(5x3-3x)/28/16f ab

5/16Formules de Newton-Cotes

Int´egration num´erique sur[a,b]surn+ 1points r´eguli`erement ´echantillonn´es.

Soitxi=a+iΔavecΔ =(b-a)n

=hn etfi=f(xi),?i.

La formule de degr´enest d´efinie par

I=? b a f(x) dx≈n? i=0w if(xi) o`u leswisont appel´es coefficients de quadrature. On d´eduit les poidswid"une interpolation de Lagrange de la fonction : f(x)≈L(x) =n? i=0f(xi)li(x),avecli(x) =n? j=0,j?=ix-xjx i-xj

L"int´egrale devient donc

I≈?

b a

L(x) dx=?

b an i=0f(xi)li(x) dx=n? i=0f(xi)? b a l i(x) dx w i.(1)

6/16Formule de Newton-Cotes (2)

I=? b a f(x) dx≈n? i=0w if(xi) Permet de retrouver les formules de quadrature simple :Degr´eNomFormuleErreur

1Trap`ezeh

2 (f0+f1)- h312 f(2)(μ)2Simpsonh 6 (f0+ 4f1+f2)- h52880 f(4)(μ)4Booleh 90
(7f0+ 32f1+ 12f2+ 32f3+ 7f4)- h71935360 f(6)(μ)pourμ?[a,b]eth= (b-a). Formule d´efinie pour n"importe quel degr´en. Probl`eme d"instabilit´e num´erique pour de grandn(Ph´enom`ene de Runge). En pratique, on pr´ef`ere d´ecouper la fonction en petits intervalles et utiliser des quadratures de degr´e faible sur ces intervalles.

7/16Quadrature de Gauss

On approche encore une fois l"int´egrale par

I=? b a f(x)dx≈n? i=1w if(xi) Leswisont les coefficients de quadrature et lesxisont choisis comme les racine de polynˆomes orthogonaux. On n"utilise pas d"´echantillonnage r´egulier, possibilit´e d"avoir une meilleure approximation def. Int´egration exacte pour des polynˆomes de degr´e2n-1.

Exemple pour[a,b] = [-1,1]avec les polynˆomes de Legendre :Nb de pointsPoidswiPointsxiPoly. de Legendre

120x
21,1-

1⎷3

,1⎷3(3x2-1)/235 9 ,89 ,59- ⎷3⎷5 ,0,⎷3⎷5(5x3-3x)/28/16

M´ethodes compositesf

a b x0x1x2x3m0m1m2 h? Les formules de Newton Cotes ont toutes une erreur sous la forme d"une puissance de(b-a). En pratique on d´ecoupe[a,b]ennsous intervalles et on utilise les formules de

Newton-Cotes sur les petits intervalles.

La longueur de l"intervalle d"int´egration devient donch=b-an Nous d´efinissons les points d"´echantillonnage r´egulier suivants : x k=a+kh, mk=a+ (k+ 1/2)h,?k?0,...,n 9/16 f a b x0x1x2x3m0m1m2 h f f

M´ethodes adaptatives

Principe

´Echantillonnage fin pas n´ecessaire sur tout l"intervalle[a,b]. On adapte l"´echantillonnage le long de l"intervalle de mani`ere r´ecursive. integrate(f,a,b,tau):

1.Calcul deI≈?

a,bf(x)dx

2.Estimation de l"erreur?≈ |I-?

a,bf(x)dx|

3.Si? > τ,

4.RetournerI

Approches classiques :

M´ethode de Romberg (Trap`eze+ extrapolation de Richardson).

?M´ethode de Simpson adaptative (Simpson+ extrapolation de Richardson).M´ethode de Simpson adaptative

Sur l"intervalle[a,b]avecm=a+b2

, la m´ethode de Simpson s"appelle avecS(a,b).

1.CalculerS(a,b),S(a,m)etS(m,b).

2.Si l"erreur|S(a,b)-S(a,m)-S(m,b)|/15> τon divise l"intervalle[a,b].

3.Sinon, on retourneS(a,m) +S(m,b) + (S(a,m) +S(m,b)-S(a,b))/15.13/16M´ethode de Simpson adaptative en Python

Impl´ementation simplifi´ee

1defsimpson(f,a,b): 2return(b-a)/6 *(f(a)+4*f((a+b)/2)+f(b))

3

4defint_adaptsimpson( f,a,b,tau):5m=(a+b)/26Sab=simpson(f,a,b)7Sam=simpson(f,a,m)8Smb=simpson(f,m,b)9ifabs (Sab-Sam-Smb)/15

11else:12returnint_adaptsimp son(f,a,m,tau)+int_adaptsimpson(f

,m,b,tau)Discussion

Fonction r´ecursive (complexit´e

d´epend defetτ).

Param`etre de pr´ecisionτdirectement

li´e `a l"erreur acceptable.

Impl´ementation non efficace,

comment faire mieux??

Minimiser les appels `af.

Stocker les valeurs de la fonction lors

des calculs de Simpson.

14/16M´ethodes de Monte Carlo

I=? b a f(x)dx≈b-an n i=1f(xi), xi≂U(a,b)?i Int´egration num´erique utilisant des r´ealisations de variables al´eatoires. La m´ethode converge vers la bonne valeur lorsquen→ ∞

L"erreur en esp´erance?=O?

1⎷n

d´ecroˆıt tr`es lentement. Tr`es utilis´ee pour les int´egrales multiples car erreur ind´ependante de la dimension.

Version plus efficaces bas´ees sur l"´echantillonnage pr´ef´erentiel (VEGAS, MISER).Impl´ementation

1defint_montecarlo( f,a,b,n):2res=03fori in range (n):4res+=f(a+(b-a)*np.random.rand())

5returnres/n Discussion

?nnombre de r´ealisations. np.random.rand()r´ealisation d"une variable al´eatoire uniforme sur[0,1]. a+np.random.rand() (b-a) r´ealisation d"une variable al´eatoire uniforme sur[a,b].15/16Comparaison num´erique I=? 1 0

4?1-x2dx=π100101102103104105106Nb d'appel de la fonction f10-1610-1410-1210-1010-810-610-410-2100ErreurVitesse de convergence

Rectangle

Pt. Milieu

Simpson

Gauss (deg 3)

Gauss (deg n)

Monte Carlo

Trapèze

Adapt. Simpson

Adapt. Gauss

16/16quotesdbs_dbs27.pdfusesText_33

[PDF] exercice python euler

[PDF] le médecin malgré lui acte 2 scène 4

[PDF] méthode dichotomie python

[PDF] le message andrée chedid résumé détaillé

[PDF] résolution équation différentielle matlab ode45

[PDF] le message andrée chedid genre

[PDF] algorithme méthode d'euler implicite matlab

[PDF] méthode de tir équation différentielle

[PDF] le message andrée chedid quiz

[PDF] le message andrée chedid extrait

[PDF] méthode euler implicite matlab

[PDF] le message andrée chedid texte intégral

[PDF] fonction ode45 matlab

[PDF] memoire de fin detude en telecommunication

[PDF] grille évaluation projet