tel/TPmatlab/integration/TP3g/demo_compar_methode_poly_gauss.m

   1 % demo_compar_methode_poly_gauss : ce script permet de comparer les deux méthodes\r
   2 %  de calcul des polynômes orthogonaux,\r
   3 % (par récurrence et directement) en terme de temps de calcul en faisant un graphe\r
   4 %\r
   5 %\r
   6 % ************ Fonctions auxiliaires utilisées ************\r
   7 %\r
   8 %    calcul_poly_gauss_symb, calcul_poly_gauss_symb_dir\r
   9 %\r
  10 % *********************************************************\r
  11 \r
  12 clear all;\r
  13 c=input('entrez c : ');\r
  14 n=input('entrez n : ');\r
  15 temps1=zeros(1,n+1);\r
  16 nombreop1=zeros(1,n+1);\r
  17 temps2=zeros(1,n+1);\r
  18 nombreop2=zeros(1,n+1);\r
  19 for k=0:n\r
  20    tic;\r
  21    d1=calcul_poly_gauss_symb_dir(c,k);\r
  22    temps1(k+1)=toc;\r
  23    tic;\r
  24    d2=calcul_poly_gauss_symb(c,k);\r
  25    temps2(k+1)=toc;\r
  26 end\r
  27 plot(1:n+1,temps1,1:n+1,temps2,'-.');\r
  28 legend('méthode directe','méthode par récurrence')\r
  29 switch c\r
  30    case 1\r
  31       title('polynômes de Gauss-Legendre');   \r
  32    case 2\r
  33       title('polynômes de Gauss-Tchebychev');   \r
  34    case 3\r
  35       title('polynômes de Gauss-Hermite');   \r
  36    case 4\r
  37       title('polynômes de Gauss-Laguerre');   \r
  38 end \r
  39 xlabel('degré');\r
  40 ylabel('temps');\r