Mathc matrices/05a
Apparence
Installer et compiler ces fichiers dans votre répertoire de travail.
c00a.c |
---|
/* ------------------------------------ */
/* Save as : c00a.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define RCA RC3
/* ------------------------------------ */
#define EV R1
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+0.390117647059, +0.169411764706, +0.457411764706,
+0.169411764706, +0.952941176471, -0.127058823529,
+0.457411764706, -0.127058823529, +0.656941176471
};
double evalue[RCA*C1]={
+1,
+1,
+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 = i_mR(RCA,RCA);
double **AmnssIde = i_mR(RCA,RCA);
double **Ab = i_Abr_Ac_bc_mR(RCA,RCA,C1);
double **b = i_mR(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(" A :");
p_mR(A,S10,P5,C10);
smul_mR(EValue[EV][C1],Ide,sIde);
sub_mR(A,sIde,AmnssIde);
c_A_b_Ab_mR(AmnssIde,b,Ab);
gj_PP_mR(Ab,NO);
// Ab : You can find the equation of the plan when the eigenvalues equal 1 (EV=R1)
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);
Normalize_mR(V);
printf(" V%d : Copy V%d into the last file",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;
}
/* ------------------------------------ */
/* ------------------------------------ */
Projection de l'espace sur un plan * Calculer un vecteur propre : * #define EV R1 pour calculer V1 et V2 * #define EV R3 pour calculer V3
Exemple de sortie écran :
A :
+0.39012 +0.16941 +0.45741
+0.16941 +0.95294 -0.12706
+0.45741 -0.12706 +0.65694
V1 : Copy V1 into the last file
+0.267643863786, +0.600000000000, -nan,
+0.963517909630, +0.000000000000, -nan,
+0.000000000000, +0.800000000000, -nan
Press return to continue.