Mathc complexes/a310

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

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

	c07a.c

/* ------------------------------------ */
/*  Save as :  c07a.c                   */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define  TAB  R4 
#define  RCA  R3  
/* ------------------------------------ */
void fun(void)
{
double **A[TAB];
double **P[TAB];
double **InvP[TAB];

double **T         = i_mZ(RCA,RCA);
double **EigsValue = i_mZ(RCA,C1);

int c;
int Yes = 0;

  for(c=C0; c<TAB; c++)
     {   
         A[c] = i_mZ(RCA,RCA); 
         P[c] = r_mZ(i_mZ(RCA,RCA),999); 
      InvP[c] = inv_mZ(P[c],i_mZ(RCA,RCA));       
	   } 

  r_mZ(A[0],99);
  
  for(c=C0; c<(TAB-C1); c++)
     {   
        mul_mZ(InvP[c],A[c],T); 
        mul_mZ(T,P[c],A[c+C1]);
	   } 

  clrscrn();
  for(c=C0; c<TAB; c++)
     {
	  printf(" A[%d] : ",c); 
      p_mZ(A[c],S10,P2,S10,P2,C6);
      
	  printf(" EigsValue  of A[%d] : ",c); 
	  eigs_mZ(A[c],EigsValue); 
      p_mZ(EigsValue,S12,P2,S12,P2,C6); 
      if(Yes){
			    stop();
			  clrscrn();
			   Yes=0; }
      else{       
		       Yes=1;}    
     } 

  for(c=C0; c<TAB; c++)
     { 
	   f_mZ(A[c]);
	   f_mZ(P[c]);		 
	   f_mZ(InvP[c]); 
	   }
	   
  f_mZ(T);       
  f_mZ(EigsValue); 
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
  do{
        fun();
        
  }while(stop_w());

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

Quatre Matrices semblables. Elles ont toutes les mêmes valeurs propres.

Exemple de sortie écran :

 A[0] : 
    +90.00    +61.00i     +69.00    -91.00i     +71.00    +41.00i 
    -10.00    -94.00i     +25.00    +35.00i     -18.00    +74.00i 
     +2.00    +78.00i      +6.00    -56.00i     +72.00    -91.00i 

 EigsValue  of A[0] : 
      +76.90     +129.96i 
      +75.37     -106.91i 
      +34.73      -18.05i 

 A[1] : 
   +109.37    +50.83i    +358.80     -9.78i    +235.12    +82.33i 
    -53.02    -45.07i     -54.41    -49.20i     -67.56     -9.97i 
    +42.72    +23.19i    +150.19    -11.97i    +132.04     +3.37i 

 EigsValue  of A[1] : 
      +76.90     +129.96i 
      +75.37     -106.91i 
      +34.73      -18.05i 

 Press return to continue. 


 A[2] : 
    -48.60   -257.43i    +348.63   -328.38i    +137.33    +75.65i 
   +172.32    -31.09i    +203.95   +297.08i     -80.06    +72.80i 
    -31.98     -5.24i     -11.70     -8.98i     +31.64    -34.64i 

 EigsValue  of A[2] : 
      +76.90     +129.96i 
      +75.37     -106.91i 
      +34.73      -18.05i 

 A[3] : 
   +143.74   +160.63i     +74.80   +183.78i     +65.30   -211.37i 
    +24.47   +216.90i     +33.76   +162.80i    +104.10   -148.07i 
   +224.07   +160.07i    +187.53   +189.28i      +9.50   -318.42i 

 EigsValue  of A[3] : 
      +76.90     +129.96i 
      +75.37     -106.91i 
      +34.73      -18.05i 

 Press return to continue.