Mathc complexes/03i

Quelques propriétés des valeurs propres. Matrices non symétrique conjuguée

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

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A            = i_mZ(r,r );
double **EigsValue    = i_mZ(r,C1);
double **ACT          = i_mZ(r,r );
double **EigsValueACT = i_mZ(r,C1);

   do{
   r_mZ(A,9);
   printf(".");
  }while(!det_Z(A).r||!det_Z(A).i);
  
  eigs_mZ(A,EigsValue);
  ctranspose_mZ(A,ACT);
  eigs_mZ(ACT,EigsValueACT);
 
  clrscrn();

  printf(" A :");
  p_mZ(A,S5,P0,S4,P0,C6);
  
  printf(" EigsValue of A :");
  p_mZ(EigsValue,S9,P3,S8,P3,C6);
      
  printf(" Conjugate transpose of A ;");
  p_mZ(ACT,S5,P0,S4,P0,C6);
  
  printf(" EigsValue of Conjugate transpose of A :");
  p_mZ(EigsValueACT,S9,P3,S8,P3,C6);   
      
  f_mZ(EigsValueACT); 
  f_mZ(EigsValue);  
  f_mZ(ACT);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(RC2)+RC2);

} while(stop_w());


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

Les valeurs propres de La matrice A et sa transposé conjugué sont conjugués  :

Exemple de sortie écran :

 A :
   +4  -6i    -8  +6i    -2  -8i    +2  +8i 
   +1  -4i    +1  -2i    +2  +2i    -9  +3i 
   -8  +9i    -8  +9i    +1  -7i    +3  -9i 
   -6  +2i    +9  -3i    -6  +5i    +9  -6i 

 EigsValue of A :
   -5.974 -19.168i 
  +16.874  -4.022i 
   +6.646  +4.018i 
   -2.546  -1.828i 

 Conjugate transpose of A ;
   +4  +6i    +1  +4i    -8  -9i    -6  -2i 
   -8  -6i    +1  +2i    -8  -9i    +9  +3i 
   -2  +8i    +2  -2i    +1  +7i    -6  -5i 
   +2  -8i    -9  -3i    +3  +9i    +9  +6i 

 EigsValue of Conjugate transpose of A :
   -5.974 +19.168i 
  +16.874  +4.022i 
   +6.646  -4.018i 
   -2.546  +1.828i 


 Press   return to continue
 Press X return to stop