Mathc matrices/c05k

Un livre de Wikilivres.
Aller à la navigation Aller à la recherche


Sommaire


Installer et compiler ces fichiers dans votre répertoire de travail.

Crystal Clear mimetype source c.png c05k.c
'
/* ------------------------------------ */
/*  Save as :   c05k.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **u       = r_mR(i_mR(r,C1),9.);
double **v       = r_mR(i_mR(r,C1),9.);
double **projuv  =      i_mR(r,C1);
double **umnsprojuv  =  i_mR(r,C1);

  clrscrn();
  
  printf(" u :");
  p_mR(u,S3,P0,C6);
  printf(" v :");
  p_mR(v,S3,P0,C6);
  stop();
  
  clrscrn();
  
  proj_uv_mR(u,v,projuv);
  
  printf(" vector component of u along v \n\n"); 
   
  printf(" projuv :           [ (<u,v> / ||v||^2) * v  ]");
  p_mR(projuv,S3,P3,C6);
  stop();
  
  clrscrn();
  printf(" vector component of (u-projuv) is orthogonal to v \n\n");
  printf(" u - projuv :");
  sub_mR(u,projuv,umnsprojuv);
  p_mR(umnsprojuv,S3,P3,C6);
  printf(" <v,u-projuv> = %.6f \n\n", dot_uv_R(v,umnsprojuv));

  f_mR(u);
  f_mR(v);
  f_mR(projuv);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(R3)+R2);


} while(stop_w());

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


Nous vérifions une des propriétés du projeté de u sur v : <v,u-projuv> = -0.000000 



Exemple de sortie écran :
 u :
 -7 
 -5 
 -9 
 +6 
 +8 

 v :
 +4 
 +2 
 +4 
 +2 
 -5 

 vector component of u along v 

 projuv :
-6.277 
-3.138 
-6.277 
-3.138 
+7.846 

 

 vector component of u orthogonal to v 

 u - projuv :
-0.723 
-1.862 
-2.723 
+9.138 
+0.154 

 <v,u - projuv> = -0.000000 


 Press return to continue
 Press X      to stop