Mathc matrices/08w

Application

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

	c00b.c

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */     
#define RCA          RC3  
/* ------------------------------------ */       
/* ------------------------------------ */
int main(void)
{                          
double NVN[RCA*RCA] ={   
+0.784464540553, -0.408248290464, -0.557086014531, 
+0.588348405415, -0.408248290464, -0.742781352708, 
-0.196116135138, +0.816496580928, +0.371390676354  
};
 
double E[RCA*RCA] ={   
+1, 0, +0, 
+0, 1, +0, 
+0, 0, +0    
};
                       
double **Nvn    =  ca_A_mR(NVN,        i_mR(RCA,RCA));
double **invNvn = invgj_mR(Nvn,        i_mR(RCA,RCA));

double **EValue =  ca_A_mR(E,          i_mR(RCA,RCA));

double **T      =   mul_mR(Nvn,EValue, i_mR(RCA,RCA));
double **A      =   mul_mR(T,  invNvn, i_mR(RCA,RCA));

  clrscrn(); 
  printf(" Nv1,Nv2: The normalized plan vectors\n" 
         "     Nv3: The normalized projector vector\n\n");
  
  printf(" Nvn:  v1        v2        v3");
  p_mR(Nvn, S8,P6, C3);     

  printf(" EValue:     E1              E2              E3");
  p_mR(EValue, S9,P12, C4); 
 
  printf(" Copy this code into the files c00c.c and c00d.c\n\n"
         " A = Nvn EValue invNvn");
  P_mR(A, S8,P12, C3);
  stop();
  
  f_mR(A);
  f_mR(Nvn);
  f_mR(invNvn);
  f_mR(EValue);
  f_mR(T);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Vous calculez la matrice A qui projette le plan dans la direction du vecteur propre (Nv3) sur un plan déterminé par les vecteurs propres (Nv1,Nv2). Ensuite, vous devez copier cette matrice dans les fichiers c00c.c et c00d.c.

Exemple de sortie écran :

  
 Nv1,Nv2: The normalized plan vectors
     Nv3: The normalized projector vector

 Nvn:  v1        v2        v3
+0.784465 -0.408248 -0.557086 
+0.588348 -0.408248 -0.742781 
-0.196116 +0.816497 +0.371391 

 EValue:     E1              E2              E3
+1.000000000000 +0.000000000000 +0.000000000000 
+0.000000000000 +1.000000000000 +0.000000000000 
+0.000000000000 +0.000000000000 +0.000000000000 

 Copy this code into the files c00c.c and c00d.c

 A = Nvn EValue invNvn
+2.363636363638, -1.909090909092, -0.272727272727, 
+1.818181818184, -1.545454545456, -0.363636363636, 
-0.909090909092, +1.272727272728, +1.181818181818  

 Press return to continue.