Mathc matrices/04i
Installer et compiler ces fichiers dans votre répertoire de travail.
|
c04b.c |
|---|
/* ------------------------------------ */
/* Save as : c04b.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define RCA RC4
/* ------------------------------------ */
#define EV R2
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+9.825339522090, +7.309094034726, -3.954100051573, -1.198212136840,
+7.309094034726, +5.437252879491, -2.941464672512, -0.891352931064,
-3.954100051573, -2.941464672512, +1.591284167096, +0.482207323363,
-1.198212136840, -0.891352931064, +0.482207323363, +0.146123431322
};
double evalue[RCA*C1]={
+17,
-0,
-0,
-0
};
double **A = ca_A_mR(a, i_mR(RCA,RCA));
double **EValue = ca_A_mR(evalue, i_mR(RCA,C1));
double **Ide = eye_mR( i_mR(RCA,RCA));
double **sIde = smul_mR(EValue[EV][C1],Ide, i_mR(RCA,RCA));
double **AmnssIde = sub_mR(A,sIde, i_mR(RCA,RCA));
double **b = m0_mR(i_mR(RCA,C1));
double **Ab = c_A_b_Ab_mR(AmnssIde,b, i_Abr_Ac_bc_mR(RCA,RCA,C1));
double **Ab_free = i_Abr_Ac_bc_mR(RCA,RCA,RCA);
double **b_free = i_mR(RCA,RCA);
double **V = i_mR(RCA,RCA);
clrscrn();
printf(" The eigenvectors associated with eigenvalues\n"
" with multiple multiplicity are not unique. \n"
" There exist several set of unit vectors that\n"
" correspond to the same eigenspace. \n\n\n"
" A:");
p_mR(A,S10,P5,C10);
gj_PP_mR(Ab,YES);
put_zeroR_mR(Ab,Ab_free);
put_freeV_mR(Ab_free);
gj_PP_mR(Ab_free,YES);
c_Ab_b_mR(Ab_free,b_free);
c_c_mR(b_free,C2,V,EV);
c_c_mR(b_free,C3,V,EV+C1);
c_c_mR(b_free,C4,V,EV+C2);
Normalize_mR(V);
printf("\n\n\n"
" Copy V%d into the last file\n\n"
" V%d:",EV,EV);
P_mR(V,S10,P12,C10);
stop();
f_mR(A);
f_mR(EValue);
f_mR(Ide);
f_mR(sIde);
f_mR(AmnssIde);
f_mR(Ab);
f_mR(b);
f_mR(Ab_free);
f_mR(b_free);
f_mR(V);
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Calculer les trois premiers vecteurs propres
Exemple de sortie écran :
The eigenvectors associated with eigenvalues
with multiple multiplicity are not unique.
There exist several set of unit vectors that
correspond to the same eigenspace.
A:
+9.82534 +7.30909 -3.95410 -1.19821
+7.30909 +5.43725 -2.94146 -0.89135
-3.95410 -2.94146 +1.59128 +0.48221
-1.19821 -0.89135 +0.48221 +0.14612
Copy V2 into the last file
V2:
-nan, -0.596864313004, +0.373340433764, +0.121054375311,
-nan, +0.802342191251, +0.000000000000, +0.000000000000,
-nan, +0.000000000000, +0.927694411171, +0.000000000000,
-nan, +0.000000000000, +0.000000000000, +0.992645877551
Press return to continue.