Mathc complexes/a309

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

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

	c06a.c

/* ------------------------------------ */
/*  Save as :  c06a.c                   */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double p = rp_I(2)+1;

double **A            =    r_mZ(    i_mZ(r, r),5.);
double **EigsValueA   = eigs_mZ(A,  i_mZ(r,C1));
double **A_p          =  pow_mZ(p,A,i_mZ(r, r));
double **EigsValueA_p = eigs_mZ(A_p,i_mZ(r,C1));

int  i =         0;
nb_Z z = i_Z(0, 0);

  clrscrn();

  printf(" A :");
  p_mZ(A,S4,P0,S3,P0,C6);
       
  printf(" EigsValue :");
  p_mZ(EigsValueA,S10,P3,S10,P3,C6);

  printf(" A_%.0f :",p);
  p_mZ(A_p,S8,P0,S8,P0,C6);
       
  printf(" EigsValueA_%.0f :",p);
  p_mZ(EigsValueA_p,S10,P3,S10,P3,C6); 
  
  stop(); 
   
  clrscrn();
  
  printf(" EigsValue A_%.0f :",p);
  p_mZ(EigsValueA_p,S10,P3,S10,P3,C6);   

  printf("  (EigsValues of A)**%.0f * :\n",p);    
  for(i=R1;i<=r;i++)
     {
	  z = i_Z(EigsValueA[i][C1], EigsValueA[i][C2]);
	  	 
       p_Z(pow_Z(p,z),S10,P3,S10,P3);printf("\n");
      } 

  f_mZ(A);
  f_mZ(EigsValueA);
  f_mZ(A_p);
  f_mZ(EigsValueA_p);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
 fun(rp_I(R3)+R1);

} while(stop_w());

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

Les valeurs propres de A**P correspondent aux valeurs propres de A à la puissance P :

Exemple de sortie écran :

 A :
  -4 +2i   -2 +2i   -3 -5i 
  -5 +3i   -4 -3i   +1 +2i 
  +5 -3i   -4 -3i   +4 +3i 

 EigsValue :
    -1.426    +7.178i 
    -6.228    -0.218i 
    +3.654    -4.960i 

 A_3 :
     +76      -8i      -18     -26i      -86    +248i 
     -50    +200i     -217     -49i     -131     -57i 
    -374     -68i       -5     +55i     -103    -371i 

 EigsValueA_3 :
  +217.513  -326.015i 
  -240.649   -25.356i 
  -220.864   -76.630i 

 Press return to continue. 


 EigsValue A_3 :
  +217.513  -326.015i 
  -240.649   -25.356i 
  -220.864   -76.630i 

  (EigsValues of A)**3 * :
  +217.513  -326.015i 
  -240.648   -25.353i 
  -220.866   -76.628i 

 Press   return to continue
 Press X return to stop