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.