Mathc matrices/c12an6
Installer et compiler ces fichiers dans votre répertoire de travail.
cgn6e.c |
---|
/* ------------------------------------ */
/* Save as : cgn6e.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
void fun(int r,int c)
{
double **B = r_mR(i_mR(c,r),9.);
double **B_T = transpose_mR(B,i_mR(r,c));
double **U_TB = i_mR(r,r); /* rc cr */
double **U = i_mR(c,r);
double **V = i_mR(r,r);
double **V_T = i_mR(r,r);
double **S = i_mR(r,r);
double **U_T = i_mR(r,c);
double **US = i_mR(c,r); /* cr rr U*S */
double **B_2 = i_mR(c,r); /* cr rr rr B = U*S*V_T */
clrscrn();
svd_U_Rn_mR(B_T,U);
svd_V_Rn_mR(B_T,V);
/* S = U_T * B * V */
transpose_mR(U,U_T);
mul_mR(U_T,B,U_TB);
mul_mR(U_TB,V,S);
/* B = U * S * V_T */
transpose_mR(V,V_T);
mul_mR(U,S,US);
mul_mR(US,V_T,B_2);
printf(" B :");
p_mR(B ,S10,P3,C10);
printf(" B = U * S * V_T ");
p_mR(B_2,S10,P3,C10);
f_mR(B);
f_mR(B_T);
f_mR(U);
f_mR(V);
f_mR(S);
f_mR(US);
f_mR(V_T);
f_mR(U_T);
f_mR(U_TB);
f_mR(B_2);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
int i;
do
{
i = rp_I(R3)+R1;
fun(i,i+R6);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de sortie écran :
------------------------------------
B :
+8.000 -1.000 +8.000
+8.000 -3.000 -1.000
-1.000 -7.000 +2.000
-3.000 +4.000 +6.000
+4.000 -1.000 +2.000
-1.000 +4.000 -1.000
-5.000 +6.000 -1.000
+4.000 -1.000 +6.000
-1.000 -3.000 +8.000
B = U * S * V_T
+8.000 -1.000 +8.000
+8.000 -3.000 -1.000
-1.000 -7.000 +2.000
-3.000 +4.000 +6.000
+4.000 -1.000 +2.000
-1.000 +4.000 -1.000
-5.000 +6.000 -1.000
+4.000 -1.000 +6.000
-1.000 -3.000 +8.000
Press return to continue
Press X to stop