Mathc matrices/c12fn2

Application

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

	c00d.c

/* ------------------------------------ */
/*  Save as :   c00d.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define FACTOR_E        +1.E-5         
/* ------------------------------------ */
void fun(int r,int c)
{
double **B       =       r_mR(        i_mR(c,r),99999.);
double **Pinv    = Pinv_Rn_mR(B,      i_mR(r,c),FACTOR_E);           
double **PinvB   =     mul_mR(Pinv,B, i_mR(r,r)); 

  clrscrn();
  printf(" Copy/Paste into the octave windows \n\n\n");
  p_Octave_mR(B,"B",P2);
  printf("pinv(B)\n\n\n");
  stop();

  clrscrn();  
  printf(" Pinv = V invS_T U_T "); 
  pE_mR(Pinv,S12,P4,C10);    
  
  printf(" Ide = Pinv B ");   
  p_mR(PinvB,S10,P6,C10);  
    
  f_mR(B);  
  f_mR(Pinv);
  f_mR(PinvB); 
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

int i;

do
{
  i = rp_I(R2)+R1;
  
  fun(i,i+R2);

} while(stop_w());

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

Exemple de sortie écran :

 Copy/Paste into the octave windows 


 B=[
-68190.00,+54828.00,+46667.00;
+16277.00,+84306.00,+28568.00;
+67372.00,+82585.00,+45071.00;
+54863.00,-41159.00,+71427.00;
-71477.00,+81278.00,-74086.00]

pinv(B)


 Press return to continue. 


 Pinv = V invS_T U_T 
 -7.9358e-06  +1.3469e-06  +5.0941e-06  +3.5079e-07  -1.0422e-06 
 +7.2578e-07  +3.6124e-06  +4.2543e-06  -1.5167e-06  +2.9761e-06 
 +7.8459e-06  +1.0701e-06  -1.4956e-07  +4.3334e-06  -4.0562e-06 

 Ide = Pinv B 
 +1.000000  +0.000000  +0.000000 
 +0.000000  +1.000000  +0.000000 
 +0.000000  -0.000000  +1.000000 


 Press   return to continue
 Press X return to stop