Mathc matrices/c22u

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

	`c00d.c`

/* ------------------------------------ */
/*  Save as :   c00d.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
#define   RA R3
#define   CA C2
#define   Cb C1 
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double ab[RA*(CA+Cb)]={
 1, -1, 4,
 3,  2, 1,
-2,  4, 3
};

double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(RA,CA,Cb));
double **A  = c_Ab_A_mR(Ab,i_mR(RA,CA));
double **b  = c_Ab_b_mR(Ab,i_mR(RA,Cb));

double **A_T        = i_mR(CA,RA);
double **A_TA       = i_mR(CA,CA); //         A_T*A
double **invA_TA    = i_mR(CA,CA); //     inv(A_T*A)
double **invA_TAA_T = i_mR(CA,RA); //     inv(A_T*A)*A_T

double **x          = i_mR(CA,Cb); // x = inv(A_T*A)*A_T*b

  clrscrn();
  printf("Unique Least Squares Solution :\n\n");
  printf(" A :");
  p_mR(A,S5,P1,C7);
  printf(" b :");
  p_mR(b,S5,P1,C7);
  printf(" Ab :");
  p_mR(Ab,S5,P1,C7);
  stop();
  
  clrscrn();
  printf(" A_T :");
  p_mR(transpose_mR(A,A_T),S5,P1,C7);
  printf(" A_TA :");
  p_mR(mul_mR(A_T,A,A_TA),S5,P1,C7);
  printf(" inv(A_TA) :");
  p_mR(inv_mR(A_TA,invA_TA),S5,P4,C7);  
  printf(" inv(A_TA)*A_T :");
  p_mR(mul_mR(invA_TA,A_T,invA_TAA_T),S5,P4,C7);
  printf("\n x = inv(A_TA)*A_T*b :");
  p_mR(mul_mR(invA_TAA_T,b,x),S5,P4,C7);     
  stop();  
  
  f_mR(A);
  f_mR(b);
  f_mR(Ab);

  f_mR(A_T);
  f_mR(A_TA);       //         A_T*A
  f_mR(invA_TA);    //     inv(A_T*A)
  f_mR(invA_TAA_T); //     inv(A_T*A)*A_T
    
  f_mR(x); 
}
/* ------------------------------------ */
int main(void)
{
	
  fun();

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

Ici le nombre de lignes doit être supérieurs aux nombres de colonnes. La fonction gaussjordan ne fonctionne pas dans cette situation.



Exemple de sortie écran :

 -----------------------------------
Unique Least Squares Solution :

 A :
 +1.0  -1.0 
 +3.0  +2.0 
 -2.0  +4.0 

 b :
 +4.0 
 +1.0 
 +3.0 

 Ab :
 +1.0  -1.0  +4.0 
 +3.0  +2.0  +1.0 
 -2.0  +4.0  +3.0 

 Press return to continue. 

 -----------------------------------
 A_T :
 +1.0  +3.0  -2.0 
 -1.0  +2.0  +4.0 

 A_TA :
+14.0  -3.0 
 -3.0 +21.0 

 inv(A_TA) :
+0.0737 +0.0105 
+0.0105 +0.0491 

 inv(A_TA)*A_T :
+0.0632 +0.2421 -0.1053 
-0.0386 +0.1298 +0.1754 


 x = inv(A_TA)*A_T*b :
+0.1789 
+0.5018 

 Press return to continue.