Aller au contenu

Mathc complexes/a311

Un livre de Wikilivres.


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         = r_mZ(i_mZ(r,r),999.);
double **EigsValue =      i_mZ(r,C1);

  clrscrn();

  printf(" A :");
  p_mZ(A,S5,P0,S4,P0,C6);
    
  clrscrn();
  printf(" Copy/Past into the octave windows \n\n\n");
  p_Octave_mZ(A,"a",P0,P0);  
  printf(" EigenValues  = eigs (a,%d) \n\n\n",r);
         
  printf(" EigsValue :");
  p_mZ(eigs_mZ(A,EigsValue),S10,P3,S10,P3,C6);    
       
  f_mZ(A);
  f_mZ(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
 fun(rp_I(R4)+R2);

} while(stop_w());

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


Vérifions avec octave les valeurs propres de A d'une matrice non symétrique conjuguée :

Exemple de sortie écran :
 Copy/Past into the octave windows 


 a=[
+391+858*i,-762-207*i,+958-224*i;
+989-160*i,+216-417*i,+780+316*i;
+714-215*i,+658-496*i,+142-856*i]

 EigenValues  = eigs (a,3) 


 EigsValue :
 +1332.661  -958.169i 
  +204.862 +1209.408i 
  -788.523  -666.239i 


 Press   return to continue
 Press X return to stop