Mathc matrices/00t
Installer et compiler ces fichiers dans votre répertoire de travail.
|
c00a.c |
|---|
/* ------------------------------------ */
/* Save as : c00a.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **U = r_mR( i_mR(R1, r),9);
double **V = r_mR( i_mR(r ,C1),9);
double **A = rToeplitz_mR(U,V,i_mR(r , r) );
double **AInv = i_mR(r,r);
double **Q = i_mR(r,r);
double **Q_T = i_mR(r,r);
double **R = i_mR(r,r);
double **RInv = i_mR(r,r);
clrscrn();
printf(" A :");
P_mR(A,S5,P3,C7);
QR_mR(A,Q,R);
printf(" Q :");
p_mR(Q,S8,P4,C7);
printf(" R :");
p_mR(R,S8,P4,C7);
stop();
clrscrn();
printf(" Q_T = QInv :");
p_mR(transpose_mR(Q,Q_T),S8,P4,C7);
printf(" Rinv :");
p_mR(inv_mR(R,RInv),S8,P4,C7);
stop();
clrscrn();
printf(" Inverse of A :");
p_mR(inv_mR(A,AInv),S5,P4,C7);
printf(" RInv QInv = RInv Q_T :");
p_mR(mul_mR(RInv,Q_T,AInv),S5,P4,C7);
f_mR(U);
f_mR(V);
f_mR(A);
f_mR(AInv);
f_mR(Q);
f_mR(R);
f_mR(Q_T);
f_mR(RInv);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(rp_I(RC3)+C2);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Matrice de Toeplitz.
Exemple de sortie écran :
A :
-9.000, -2.000, +5.000, -3.000, -5.000,
-9.000, -9.000, -2.000, +5.000, -3.000,
+5.000, -9.000, -9.000, -2.000, +5.000,
-2.000, +5.000, -9.000, -9.000, -2.000,
+3.000, -2.000, +5.000, -9.000, -9.000
Q :
-0.6364 -0.0212 +0.2344 -0.3912 +0.6218
-0.6364 -0.5320 -0.1521 +0.1550 -0.5146
+0.3536 -0.7261 -0.4131 -0.1750 +0.3828
-0.1414 +0.3926 -0.7492 -0.4897 -0.1570
+0.2121 -0.1875 +0.4358 -0.7433 -0.4212
R :
+14.1421 +2.6870 -2.7577 -2.6163 +5.2326
+0.0000 +13.7033 +3.0219 -2.9898 -1.0260
+0.0000 +0.0000 +14.1161 +2.1833 -5.2047
-0.0000 +0.0000 -0.0000 +13.3959 +8.2851
-0.0000 -0.0000 +0.0000 +0.0000 +4.4538
Press return to continue.
Q_T = QInv :
-0.6364 -0.6364 +0.3536 -0.1414 +0.2121
-0.0212 -0.5320 -0.7261 +0.3926 -0.1875
+0.2344 -0.1521 -0.4131 -0.7492 +0.4358
-0.3912 +0.1550 -0.1750 -0.4897 -0.7433
+0.6218 -0.5146 +0.3828 -0.1570 -0.4212
Rinv :
+0.0707 -0.0139 +0.0168 +0.0080 -0.0815
-0.0000 +0.0730 -0.0156 +0.0188 -0.0365
-0.0000 -0.0000 +0.0708 -0.0115 +0.1043
+0.0000 +0.0000 +0.0000 +0.0746 -0.1389
-0.0000 -0.0000 -0.0000 -0.0000 +0.2245
Press return to continue.
Inverse of A :
-0.0946 +0.0030 -0.0045 -0.0191 +0.0533
-0.0353 -0.0148 -0.0638 +0.0369 -0.0191
+0.0859 -0.0662 +0.0127 -0.0638 -0.0045
-0.1155 +0.0830 -0.0662 -0.0148 +0.0030
+0.1396 -0.1155 +0.0859 -0.0353 -0.0946
RInv QInv = RInv Q_T :
-0.0946 +0.0030 -0.0045 -0.0191 +0.0533
-0.0353 -0.0148 -0.0638 +0.0369 -0.0191
+0.0859 -0.0662 +0.0127 -0.0638 -0.0045
-0.1155 +0.0830 -0.0662 -0.0148 +0.0030
+0.1396 -0.1155 +0.0859 -0.0353 -0.0946
Press return to continue
Press X return to stop