Mathc complexes/a311
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