Mathc complexes/a294

Quelques propriétés des valeurs propres. Matrices symétriques conjugués

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            = rcsymmetric_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,S9,P0,S5,P0,C6);

  printf(" A_%.0f :",p);
  p_mZ(A_p,S8,P0,S8,P0,C6);
       
  printf(" EigsValueA_%.0f :",p);
  p_mZ(EigsValueA_p,S9,P0,S5,P0,C6); 
  
  stop(); 
   
  clrscrn();
  
  printf(" EigsValue A_%.0f :",p);
  p_mZ(EigsValueA_p,S9,P0,S3,P0,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),S9,P0,S3,P0);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 :
 +83 +0i  -19+22i  +10+15i  -39+22i 
 -19-22i +106 +0i  -35+20i  +13+74i 
 +10-15i  -35-20i  +48 +0i  +12 -6i 
 -39-22i  +13-74i  +12 +6i  +85 +0i 

 EigsValue :
     +196   +0i 
     +100   -0i 
      +24   +0i 
       +2   -0i 

 A_2 :
  +10064      +0i    -2520   +6605i     +935    +845i    -8217   +2696i 
   -2520   -6605i   +19351      +0i    -5538   +3541i    +3408  +15024i 
    +935    -845i    -5538   -3541i    +4434      +0i    +2561   -2843i 
   -8217   -2696i    +3408  -15024i    +2561   +2843i   +15055      +0i 

 EigsValueA_2 :
   +38256   -0i 
   +10063   -0i 
     +581   +0i 
       +4   +0i 

 Press return to continue. 


 EigsValue A_2 :
   +38256 -0i 
   +10063 -0i 
     +581 +0i 
       +4 +0i 

  (EigsValues of A)**2 * :
   +38256 +0i 
   +10063 -0i 
     +581 +0i 
       +4 -0i 

 Press   return to continue
 Press X return to stop