Coursera MOOC : analyse numrique pour ingnieurs e e Sections PH, ME, IN, SC, printemps 2013 Enseignant : Marco
Picasso Chapitre 1 - Exercice 1 : nonc e e Le but de cet exercice est de complter un programme permettant dinterpoler e une fonction continue par un polynme (sect. 1.4 du livre). Soit f : [1, 1] R o une fonction continue donne et soit t0 < t1 < t2 < . . . < tn , (n + 1) points de e [1, 1]. On cherche un polynme pn de degr infrieur ou gal ` n tel que o e e e a pn (tj ) = f (tj ), 0 j n.
Question 1 (1 pt). Rappeler la formule du cours permettant dexpliciter pn en fonction de la base de Lagrange de Pn associe aux points t0 , t1 , t2 , . . . , tn . e Question 2 (1 pt). Le programme intlag.m est ` votre disposition ` la n de ce document. a a A partir de la donne de f et n, ce programme permet de calculer le e polynme pn qui interpole f aux points tj , 0 j n, quidistribus o e e uniformment sur lintervalle [1, 1]. Ce programme est incomplet, vous e devez en particulier programmer les formules issues de la question 1 cidessus. Pour lancer le programme avec n = 10, tapez [err,t,f,x,p]=intlag(10); dans la fentre de commande matlab ou octave. Ceci aura pour eet e de crer les vecteurs t,f,x,p. Les rsultats pourront tre visualiss en e e e e tapant ensuite plot(t,f,"o",x,p); Taper err dans la fentre de commande matlab ou octave pour obtenir une approxe imation de max |f (t) pn (t)|.
1t1
On consid`re le cas o` f (t) = sin(t). Les 3 lignes ` complter sont .... e u a e Question 3 (2 pts). On consid`re le cas o` f (t) = sin(t). Remplir le tableau suivant (1pt). e u
n 5 10 20 Que vaut lim
n 1t1
err
max |f (t) pn (t)|? Expliquer ce rsultat grce au thor`me e a e e
1.1 du livre (1pt). Question 4 (1 pt). On consid`re le cas o` f (t) = e u
1 1+25t2 .
Remplir le tableau suivant err
n 5 10 20 40 Que vaut lim
n 1t1
max |f (t) pn (t)|?
1 1+25t2
Question 5 (1 pt) On consid`re le cas o` f (t) = e u tj = cos
mais avec j = 1, ..., n + 1, max |f (t) pn (t)|?
(2j 1) , 2(n + 1)
qui sont les points dits de Chebyshev. Que vaut lim
n 1t1
Le chier octave intlag.m %Etant donne un entier n et une fonction continue f, %le programme interpole la fonction f par un polynome p %de degre n aux points d interpolation t(1),t(2),...,t(n). %parametres: entree : n % sortie : err : erreur max entre la fonction f et l interpolant p % t : (n+1) vecteur contenant les points d interpolation % f : (n+1) vecteur contenant les valeurs de la fonction f % aux points d interpolation % x : vecteur contenant 1001 points uniformement distribues % sur [-1,1] % p : vecteur contenant les valeurs du polynome p au point x(i) % function [err,t,f,x,p]=intlag(n) % %initialisation des vecteurs t et f % for i=1:n+1 t(i)=-1+(2.*(i-1))/n; f(i)=funct(t(i)); end % %calcul de la valeur du polynome d interpolation au point x(i) % m=1000; err=0; for i=1:m+1 x(i)=-1+(2.*(i-1))/m; p(i)=0; for j=1:n+1 p(i) = p(i) + f(j) * ?????; end err = max(err,abs(p(i)-funct(x(i)))); end %calcul de phi_j (la j_ieme fonction de la base de Lagrange) function basis = phi(j,n,xx,t) basis=1; for k=1:n+1 if k ~= j basis = basis * ?????; end end %fonction a interpoler function f = funct(xx) f = ?????;
