Mathc complexes/00q
Installer et compiler ce fichier dans votre répertoire de travail.
|
c00a.c |
|---|
/* ------------------------------------ */
/* Save as : c00a.c */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int rc)
{
double **A = r_mZ( i_mZ(rc,rc),9);
double **Q = i_mZ(rc,rc);
double **R = i_mZ(rc,rc);
double **A_T = transpose_mZ(A, i_mZ(rc,rc));
double **Q_T = i_mZ(rc,rc);
double **R_T = i_mZ(rc,rc);
clrscrn();
printf(" A : ");
p_mZ(A, S8,P4, S8,P4, C6);
QR_mZ(A,Q,R);
printf(" Q :");
p_mZ(Q, S8,P4, S8,P4, C6);
printf(" R :");
p_mZ(R, S8,P4, S8,P4, C6);
stop();
clrscrn();
printf(" A_T : Transpose");
p_mZ(A_T, S8,P4, S8,P4, C6);
QR_mZ(A_T,Q_T,R_T);
printf(" Q_T :");
p_mZ(Q_T, S8,P4, S8,P4, C6);
printf(" R_T :");
p_mZ(R_T, S8,P4, S8,P4, C6);
stop();
clrscrn();
printf(" Q :");
p_mZ(Q, S8,P4, S8,P4, C6);
printf(" R :");
p_mZ(R, S8,P4, S8,P4, C6);
printf(" Q_T :");
p_mZ(Q_T, S8,P4, S8,P4, C6);
printf(" R_T :");
p_mZ(R_T, S8,P4, S8,P4, C6);
stop();
clrscrn();
printf(" A :");
p_mZ(A, S8,P4, S8,P4, C6);
printf(" A = Q * R :");
mul_mZ(Q,R, A);
p_mZ(A, S8,P4, S8,P4, C6);
printf(" A_T :");
p_mZ(A_T, S8,P4, S8,P4, C6);
printf(" A_T = Q_T * R_T :");
mul_mZ(Q_T,R_T, A_T);
p_mZ(A_T, S8,P4, S8,P4, C6);
f_mZ(A);
f_mZ(Q);
f_mZ(R);
f_mZ(A_T);
f_mZ(Q_T);
f_mZ(R_T);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
fun(rp_I(RC3)+RC1);
while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
La QR décomposition d'une matrice et de sa transposée non aucun lien.
Exemple de sortie écran :
A :
+7.0000 -6.0000i -4.0000 +2.0000i -4.0000 +5.0000i
-7.0000 -2.0000i -4.0000 +2.0000i +5.0000 +6.0000i
-4.0000 -3.0000i +8.0000 +9.0000i -2.0000 -8.0000i
Q :
+0.5483 -0.4700i +0.0714 +0.0959i -0.1945 +0.6530i
-0.5483 -0.1567i -0.5675 -0.0688i -0.5285 +0.2624i
-0.3133 -0.2350i +0.5920 +0.5555i -0.3981 -0.1709i
R :
+12.7671 -0.0000i -5.8745 -3.4463i -5.7178 +0.3916i
-0.0000 -0.0000i +11.7734 +0.0000i -8.6839 -5.9451i
-0.0000 -0.0000i -0.0000 +0.0000i +5.1380 -0.0000i
Press return to continue.
A_T : Transpose
+7.0000 -6.0000i -7.0000 -2.0000i -4.0000 -3.0000i
-4.0000 +2.0000i -4.0000 +2.0000i +8.0000 +9.0000i
-4.0000 +5.0000i +5.0000 +6.0000i -2.0000 -8.0000i
Q_T :
+0.5793 -0.4966i -0.3081 +0.3602i -0.4394 -0.0073i
-0.3310 +0.1655i -0.7383 -0.1024i -0.1983 +0.5178i
-0.3310 +0.4138i +0.1590 +0.4410i -0.6494 -0.2787i
R_T :
+12.0830 -0.0000i -0.5793 -8.6899i -4.6346 -4.5518i
+0.0000 -0.0000i +7.6257 +0.0000i -10.5223 -3.8506i
+0.0000 +0.0000i +0.0000 +0.0000i +8.3818 -0.0000i
Press return to continue.
Q :
+0.5483 -0.4700i +0.0714 +0.0959i -0.1945 +0.6530i
-0.5483 -0.1567i -0.5675 -0.0688i -0.5285 +0.2624i
-0.3133 -0.2350i +0.5920 +0.5555i -0.3981 -0.1709i
R :
+12.7671 -0.0000i -5.8745 -3.4463i -5.7178 +0.3916i
-0.0000 -0.0000i +11.7734 +0.0000i -8.6839 -5.9451i
-0.0000 -0.0000i -0.0000 +0.0000i +5.1380 -0.0000i
Q_T :
+0.5793 -0.4966i -0.3081 +0.3602i -0.4394 -0.0073i
-0.3310 +0.1655i -0.7383 -0.1024i -0.1983 +0.5178i
-0.3310 +0.4138i +0.1590 +0.4410i -0.6494 -0.2787i
R_T :
+12.0830 -0.0000i -0.5793 -8.6899i -4.6346 -4.5518i
+0.0000 -0.0000i +7.6257 +0.0000i -10.5223 -3.8506i
+0.0000 +0.0000i +0.0000 +0.0000i +8.3818 -0.0000i
Press return to continue.
A :
+7.0000 -6.0000i -4.0000 +2.0000i -4.0000 +5.0000i
-7.0000 -2.0000i -4.0000 +2.0000i +5.0000 +6.0000i
-4.0000 -3.0000i +8.0000 +9.0000i -2.0000 -8.0000i
A = Q * R :
+7.0000 -6.0000i -4.0000 +2.0000i -4.0000 +5.0000i
-7.0000 -2.0000i -4.0000 +2.0000i +5.0000 +6.0000i
-4.0000 -3.0000i +8.0000 +9.0000i -2.0000 -8.0000i
A_T :
+7.0000 -6.0000i -7.0000 -2.0000i -4.0000 -3.0000i
-4.0000 +2.0000i -4.0000 +2.0000i +8.0000 +9.0000i
-4.0000 +5.0000i +5.0000 +6.0000i -2.0000 -8.0000i
A_T = Q_T * R_T :
+7.0000 -6.0000i -7.0000 -2.0000i -4.0000 -3.0000i
-4.0000 +2.0000i -4.0000 +2.0000i +8.0000 +9.0000i
-4.0000 +5.0000i +5.0000 +6.0000i -2.0000 -8.0000i
Press return to continue
Press X return to stop