Différences
Ci-dessous, les différences entre deux révisions de la page.
Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédente Prochaine révisionLes deux révisions suivantes | ||
teaching:progappchim:numpy_simple [2018/10/31 09:17] – [Références] villersd | teaching:progappchim:numpy_simple [2019/03/05 02:11] – villersd | ||
---|---|---|---|
Ligne 25: | Ligne 25: | ||
Les fonctions arange et shape sont bien pratiques pour générer des nombres en séquences et réarranger des listes de nombres. La fonction linspace est utile parce qu' | Les fonctions arange et shape sont bien pratiques pour générer des nombres en séquences et réarranger des listes de nombres. La fonction linspace est utile parce qu' | ||
- | Vous pouvez consulter [[http://wiki.scipy.org/Numpy_Example_List|cette page]] pour consulter d' | + | Vous pouvez consulter [[https://docs.scipy.org/doc/numpy/|cette page]] pour consulter d' |
</ | </ | ||
<code python arrays_01.py> | <code python arrays_01.py> | ||
Ligne 203: | Ligne 203: | ||
**poly1d & polynomial ordonnent les coefficients en sens inverses !!!**</ | **poly1d & polynomial ordonnent les coefficients en sens inverses !!!**</ | ||
- | <sxh python; title : arrays_polynomes_06.py> | + | <code python arrays_polynomes_06.py> |
#! / | #! / | ||
# -*- coding: utf-8 -*- | # -*- coding: utf-8 -*- | ||
Ligne 210: | Ligne 210: | ||
""" | """ | ||
import numpy as np | import numpy as np | ||
+ | from numpy.polynomial import Polynomial as P | ||
- | # les coefficients du polynômes sont donnés par ordre décroissance des degrés dans poly1d | + | # les coefficients du polynômes sont donnés par ordre décroissance des dégrés |
- | a=np.poly1d([1.,2.,3.,4.]) # = x³ + 2x² + 3x +4 | + | a = P([4., 3., 2., 1.]) # = x³ + 2x² + 3x + 4 |
- | print " | + | print(" |
# les coefficients de a : | # les coefficients de a : | ||
- | print " | + | print(" |
# les racines de a : | # les racines de a : | ||
- | print " | + | print(" |
# l' | # l' | ||
- | print "ordre : ",a.order | + | print("ordre : ", a.degree()) |
# évaluations sur un vecteur | # évaluations sur un vecteur | ||
- | x=np.linspace(0, | + | x = np.linspace(0, |
- | print "x = ",x | + | print("x = ", x) |
- | print " | + | print(" |
# dérivation | # dérivation | ||
- | print " | + | print(" |
+ | print(" | ||
+ | print(" | ||
+ | print(" | ||
# intégration | # intégration | ||
- | print " | + | print(" |
# création d'un polynôme par ses racines | # création d'un polynôme par ses racines | ||
- | b=a.roots | + | b = a.roots() |
- | c=np.poly1d(b,True) | + | c = P.fromroots(b) |
- | print " | + | print(" |
+ | # | ||
# fitting polynomial | # fitting polynomial | ||
- | xd=np.array([0., | + | # |
- | yd1=np.array([0.05, | + | # utilisation de poly1d (ancienne librairie) |
- | pfit=np.poly1d(np.polyfit(xd, | + | # |
- | print "fit d'une parabole (polynôme d' | + | # numpy.polyfit (poly1d) : |
- | print xd | + | # https:// |
- | print yd1 | + | # https:// |
- | print " | + | # |
- | </sxh> | + | xd = np.array([0., |
+ | yd = np.array([0.05, | ||
+ | pfit = np.poly1d(np.polyfit(xd, | ||
+ | print("fit d'une parabole (polynôme d' | ||
+ | print(xd) | ||
+ | print(yd) | ||
+ | print(" | ||
+ | # | ||
+ | # " | ||
+ | # in the opposite order of that for np.polyfit and np.polyval" | ||
+ | # → https:// | ||
+ | # | ||
+ | ################################################## | ||
+ | # Ajouter les fits utilisant numpy.polynomial... # | ||
+ | ################################################## | ||
+ | # | ||
+ | # numpy.polynomial.polynomial.Polynomial.fit : | ||
+ | # https:// | ||
+ | # https:// | ||
+ | # https:// | ||
+ | # https:// | ||
+ | # | ||
+ | # numpy.polynomial.polynomial.polyfit : | ||
+ | # https:// | ||
+ | # | ||
+ | </code> | ||
Autres fonctions : voir [[http:// | Autres fonctions : voir [[http:// |