tel/TPmatlab/interpolation/TP2f/test_diff_div.m

   1 function R=test_diff_div(T,X,Y) \r
   2 \r
   3 % test_diff_div : fonction test pour le calcul du polynôme interpolateur support quelconque\r
   4 %\r
   5 %\r
   6 % *********************************************************\r
   7 %\r
   8 %       R=test_diff_div(T,X,Y)  \r
   9 %     calcul de la valeur du polynôme d'interpolation défini par le nuage de point X,Y\r
  10 %     sur le tableau T (utilisation du calcul des différences divisées et \r
  11 %     de l'algorithme d'évaluation de Horner) pour support quelconque\r
  12 %\r
  13 %       variables d'entrées : \r
  14 %     * T : contient le tableau rectangulaire des réels où est évalué p_n\r
  15 %     * X : contient les valeurs x_i, pour 0 <=i<=n (dans l'ordre croissant)\r
  16 %     * Y : contient les valeurs f^(r)(x_i), pour 0 <=i<=n, de telle sorte que :\r
  17 %             si x_k est multiple avec x_{k-1}~=x_k=x_{k+1}=...=x_{k+r}~=x_{k+r+1} alors\r
  18 %                y_k=f(x_k), y_{k+1}=f'(x_k), ..., y_{k+r}=f^(r)(x_k).\r
  19 %\r
  20 %       variables de sortie\r
  21 %     * R : contient le tableau rectangulaire des images de T par p_n,\r
  22 %           polynôme d'interpolation de f sur le support x_i, pour 0 <=i<=n\r
  23 %\r
  24 % \r
  25 % ************ Fonctions auxiliaires utilisées ************\r
  26 %\r
  27 %     diff_div, eval_horner\r
  28 %\r
  29 % *********************************************************\r
  30 \r
  31 \r
  32 D=diff_div(X,Y);\r
  33 R=eval_horner(T,X(1:end-1),D);\r
  34 \r