Différences
Ci-dessous, les différences entre deux révisions de la page.
Prochaine révision | Révision précédenteDernière révisionLes deux révisions suivantes | ||
teaching:progappchim:polynomes-11 [2012/11/22 11:22] – créée villersd | teaching:progappchim:polynomes-11 [2016/02/23 13:29] – villersd | ||
---|---|---|---|
Ligne 14: | Ligne 14: | ||
cf. http:// | cf. http:// | ||
""" | """ | ||
- | 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, | for i in range(n, | ||
- | p=p*x+a[i] | + | p = p*x+a[i] |
return p | return p | ||
def polyscal(s, | def polyscal(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): | ||
""" | """ | ||
- | b=[0]+a | + | b = [0]+a # cela revient à " |
return b | return b | ||
Ligne 35: | Ligne 35: | ||
""" | """ | ||
""" | """ | ||
- | 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): | ||
| | ||
Ligne 55: | Ligne 55: | ||
cf. http:// | cf. http:// | ||
""" | """ | ||
- | rep=[[1.], | + | rep = [[1.], |
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., | + | x = arange(-1., |
- | chebs=polycheby(10) | + | chebs = polycheby(10) |
print chebs | print chebs | ||
Ligne 84: | Ligne 84: | ||
À ce stade, il est utile de s' | À ce stade, il est utile de s' | ||
+ | De plus, des modules de calcul scientifique utilisant les familles classiques de polynômes orthogonaux existent dans [[http:// | ||
- | [[polynomes-12|Réponse | + | [[polynomes-12|Suite à la page suivante !]] |