teaching:progappchim:polynomes-11

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

Les deux révisions précédentes Révision précédente
Dernière révisionLes deux révisions suivantes
teaching:progappchim:polynomes-11 [2015/03/06 11:15] – villersdteaching:progappchim:polynomes-11 [2016/02/23 13:29] – villersd
Ligne 14: Ligne 14:
     cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner     cf. http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Ruffini-Horner
     """     """
-    n=len(a)-1 # n = ordre du polynome +    n = len(a)-1 # n = ordre du polynome 
-    p=0.+    p = 0.
     for i in range(n,-1,-1):     for i in range(n,-1,-1):
-        p=p*x+a[i]+        p = p*x+a[i]
     return p     return p
  
 def polyscal(s,a): def polyscal(s,a):
     """polynôme multiplié par un scalaire s """     """polynôme multiplié par un scalaire s """
-    b=[]+    b = []
     for coef in a:     for coef in a:
         b.append(coef*s)         b.append(coef*s)
Ligne 29: Ligne 29:
 def polyshift(a): def polyshift(a):
     """Multiplication du polynôme par la variable x"""     """Multiplication du polynôme par la variable x"""
-    b=[0]+a   # cela revient à "shifter" la liste des coefficients en insérant un 0 "à gauche"+    b = [0]+a   # cela revient à "shifter" la liste des coefficients en insérant un 0 "à gauche"
     return b     return b
  
Ligne 35: Ligne 35:
     """ Addition de deux polynômes de coefficients a et b     """ Addition de deux polynômes de coefficients a et b
     """     """
-    r=a[:] # on travaille sur une copie de a pour ne pas le modifier +    r = a[:] # on travaille sur une copie de a pour ne pas le modifier 
-    t=b[:] # idem pour b    +    t = b[:] # idem pour b    
-    g=[]   # polynôme somme +    g = []   # polynôme somme 
-    n1=len(r) # ordre du premier polynôme +    n1 = len(r) # ordre du premier polynôme 
-    n2=len(t) # ordre du second polynôme +    n2 = len(t) # ordre du second polynôme 
-    if n1>n2: # premier polynôme de plus haut degré que le second+    if n1 > n2: # premier polynôme de plus haut degré que le second
         for i in range (n1-n2):         for i in range (n1-n2):
             t.append(0)             t.append(0)
-    elif n1<n2: # second polynôme de plus haut degré que le premier+    elif n1 < n2: # second polynôme de plus haut degré que le premier
          for i in range (n2-n1):          for i in range (n2-n1):
              r.append(0)              r.append(0)
Ligne 55: Ligne 55:
     cf. http://fr.wikipedia.org/wiki/Polyn%C3%B4me_de_Tchebychev pour la formule de récurrence     cf. http://fr.wikipedia.org/wiki/Polyn%C3%B4me_de_Tchebychev pour la formule de récurrence
     """     """
-    rep=[[1.],[0.,1.]] # les deux premiers polynômes (degrés 0 et 1) pour l'application de la formule de récurrence+    rep = [[1.],[0.,1.]] # les deux premiers polynômes (degrés 0 et 1) pour l'application de la formule de récurrence
     if nmax < 1: # si nmax est inférieur au degré 1, on renvoie le polynôme de degré 0     if nmax < 1: # si nmax est inférieur au degré 1, on renvoie le polynôme de degré 0
         rep=[[1.]]         rep=[[1.]]
Ligne 64: Ligne 64:
  
 # utilisation des objets numpy # utilisation des objets numpy
-x=arange(-1.,1.00001,0.01) +x = arange(-1.,1.00001,0.01) 
-chebs=polycheby(10)   # quelques premiers polynômes de Tchebyshev+chebs = polycheby(10)   # quelques premiers polynômes de Tchebyshev
 print chebs print chebs
  
  • teaching/progappchim/polynomes-11.txt
  • Dernière modification : 2017/02/24 11:58
  • de villersd