Mathc matrices/c05n

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 c05n.c
'
/* ------------------------------------ */
/*  Save as :   c05n.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r,int rn)
{
double **A      = r_mR(i_mR(r+rn,r),9.);
double **Q      =      i_mR(r+rn,r);
double **R      =      i_mR(r,r);

  clrscrn();
  printf(" A :");
  p_mR(A,S3,P3,C6);
  
  QR_uv_mR(A,Q,R);
    
  printf(" Q :");
  p_mR(Q,S3,P5,C6);
  
  printf(" R :");
  p_mR(R,S10,P5,C6);
  stop();
  
  clrscrn();
  printf(" A :");
  p_mR(A,S3,P3,C6);
  
  printf(" A = Q * R  :");  
  mul_mR(Q,R, A);
  p_mR(A,S3,P3,C6);  
  
  f_mR(A);
  f_mR(Q);
  f_mR(R);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

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

} while(stop_w());

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


Le travail précédent nous conduit à la factorisation QR : QRuv_mR(A,Q,R);



Exemple de sortie écran :
 A :
+4.000 -5.000 +8.000 
+8.000 -3.000 -9.000 
-1.000 -3.000 +6.000 
+4.000 +2.000 +2.000 

 Q :
+0.40614 -0.60845 +0.49132 
+0.81228 -0.04654 -0.55125 
-0.10153 -0.55846 +0.16779 
+0.40614 +0.56191 +0.65312 

 R :
  +9.84886   -3.35064   -3.85832 
  +0.00000   +5.98107   -6.67570 
  +0.00000   -0.00000  +11.20484 

 Press return to continue.