Mathc matrices/c05b2
Installer et compiler ces fichiers dans votre répertoire de travail.
|
b02.c |
|---|
/* ------------------------------------ */
/* Save as : b02.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r,int c)
{
double **A = r_mR(i_mR(r,c),9);
double **A_t = i_mR(c,r);
double **AA_t = i_mR(r,r);
double **A_tA = i_mR(c,c);
clrscrn();
printf(" A \n");
p_mR(A,3,0,C6);
printf(" A^t (transpose)\n");
p_mR(transpose_mR(A,A_t),S3,P0,C6);
stop();
clrscrn();
printf(" AA_t : ");
p_mR(mul_mR(A,A_t,AA_t),S8,P0,C6);
printf(" A_tA : ");
p_mR(mul_mR(A_t,A,A_tA),S8,P0,C6);
f_mR(A);
f_mR(A_t);
f_mR(AA_t);
f_mR(A_tA);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
fun(rp_I(R6),rp_I(C6));
while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Multiplions une matrice par sa transposé. Nous obtenons une matrice carré, avec des valeurs positives sur la diagonale. Exemple de sortie écran :
A
-8 +1 +5 -8 -2 -4
-3 -5 +6 -8 +3 -9
-9 -5 +6 -1 -1 +3
+5 +1 -8 +8 -8 -3
+5 +7 +1 -6 +5 +5
A^t (transpose)
-8 -3 -9 +5 +5
+1 -5 -5 +1 +7
+5 +6 +6 -8 +1
-8 -8 -1 +8 -6
-2 +3 -1 -8 +5
-4 -9 +3 -3 +5
Press return to continue.
------------------------------
AA_t :
+174 +143 +95 -115 -10
+143 +224 +66 -129 -26
+95 +66 +153 -107 -58
-115 -129 -107 +227 -79
-10 -26 -58 -79 +161
A_tA :
+204 +92 -147 +107 +1 +42
+92 +101 -56 +3 +15 +58
-147 -56 +162 -164 +71 -27
+107 +3 -164 +229 -101 +47
+1 +15 +71 -101 +103 +27
+42 +58 -27 +47 +27 +140
Press return to continue
Press X to stop