C++

Programme pour lire une matrice triangulaire inférieure à partir

Recherche :

Mot : Pseudo : Filtrer
Page : 1 2 Page Suivante Page Précédente Bas de page
Auteur	Sujet : Programme pour lire une matrice triangulaire inférieure à partir

iamora

Bonjour,

Je voudrais lire une matrice inférieure à partir d'un fichier .txt.
La forme de la matrice en question est la suivante :

Code :

4
1.23E-02
-1.24E-02 1.26E-02
1.27E-02 -1.23E-04 1.23E-03
1.27E-02 -1.23E-04 1.23E-03 1.23E-03

Le premier élément de cette matrice qui est 4 représente la dimension, le reste des éléments constituent la matrice.
N'étant pas expert en c++ je suis passé par deux codes différents mais aucun ne me donne un réusltat satisfaisant.

code 1 :

Code :

#include "essai.h"
#include <iostream>
#include <fstream>
#include <iomanip>
using namespace std;
using namespace std;
int main()
{
// Ouverture du fichier
ifstream fichierEntree("c:\\test1.txt", ios::in);
// Test d'ouverture du fichier
if (!fichierEntree) {
cerr << "Problème d'ouverture de fichier" << endl;
exit(1);
}
char* ligne;
int i = 0;
char * err = NULL;
//double z;
ligne = new char[10];
fichierEntree.getline(ligne, 4);
int conv1 = strtoul(ligne,0,10);
ArrayLengthSpecifier dimension(conv1);
cout << conv1 << endl;
SymmetricMatrix covariance(dimension);
int longueur;
do
{
i = i + 1;
longueur = 9 * i + 1;
ligne = new char [longueur];//format "0.00E+00"
//donc 8 caractères; 9i=8i+ (i-1) espaces + élt nul terminal
fichierEntree.getline(ligne,longueur);
cout << ligne << endl;
for (int j=1;j<=i;j++)
{
char temp[9];
for (int k=0;k<=7;k++)
{
temp[8]=0;
if (j==1)
temp[k] = ligne[k];
else
temp[k] = ligne[9* (j - 1) + k ];//attention dimension de 0 à n-1
}
const char* caractere=temp;
cout << caractere << endl;
covariance(i,j) = 2 * strtod(caractere, &err);
cout << covariance(i,j) << endl;
cout << covariance(j,i) << endl;
}
delete[] ligne;
}
while (!fichierEntree.eof());
fichierEntree.close();
}

code 2 :

[code]
#include "essai.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <string>

using namespace std;

int main()
{
ifstream FichierMatrice("c:\\test1.txt",ios::in);
if (!FichierMatrice)
{
cerr << "Ouverture du fichier impossible" << endl;
exit(1);
}

int val1;
int j = 0;
double* valf;
double val;

FichierMatrice >> val1;
ArrayLengthSpecifier dimension(val1);
SymmetricMatrix covariance(dimension);
cout << val1 << endl;

int dime = (val1*(val1+1))/2;
valf = new double [dime];

do
{
i = i + 1;
FichierMatrice >> valf[j];
cout << valf[j] <<endl;
//cout << valf[1]<<endl;
//cout << covariance (1,2) <<endl;
//cout << covariance (2,1) <<endl;
//cout << covariance (3,2) <<endl;
}
while (!FichierMatrice.eof());

return 0;
}
[code]
La classe SymmetricMatrix et la matrice covariance sont correctement définis à l'aide de la bibliothèque Newmat.
Je suis vraiment bloqué parceque je n'arrive pas à traiter le fichier de donnée. Pourriez-vous m'aider sur ce point en me conseillant laquelle des deux voies d'abord suivre puis comment écire ce code???

Merci d'avance

Publicité

franceso

tu aurais pu continuer sur le même topic plutôt que d'en ouvrir un autre

que donne un truc dans ce genre là :

Code :

ifstream fichierEntree("c:\\test1.txt", ios::in);
int dim;
fichierEntree >> dim;
ArrayLengthSpecifier dimension(dim);
SymmetricMatrix covariance(dimension);
/* lecture de la matrice */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=i ; ++j )
{
double val;
fichierEntree >> val;
covariance(i,j) = val;
}
}
fichierEntree.close();
/* affichage pour tester */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
std::cout << setw(10) << covariance(i,j);
}
std::cout << std::endl;
}

---------------
TriScale innov

iamora

Bonjour,
Merci pour ton message, je n'ai malheureusement pas pu le lire auparavant,
Le code que tu m'a laissé marche bieb je l'ai testé.
Par contre j'aurais une petite question concercant les lignes suivantes pour les tests car j'aimerai savoir comment ça marche pour pouvoir l'exploiter plus tard :

Code :

for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
std::cout << setw(10) << covariance(i,j);
}
std::cout << std::endl;
}

a quoi sert le setw(10)???
je voudrais utiliser ce qui sort en test dans un fichier c'est pour ça que je te pose la question, j'ai l'impression que c'est assez proche !

Message cité 1 fois

_darkalt3_

Proctopathe

iamora a écrit :

a quoi sert le setw(10)?

http://www.cplusplus.com/ref/iostr [...] /setw.html

---------------
Töp of the plöp

iamora

oui ok c'est l'espace entre les sorties je m'en suis aperçu quad j'ai testé je suis débutant en C++, alors excusez mon manque d'expérience.
Par contre j'ai encore deux petites questions. Je pense que la première est un peu évidente pour vous !

Pensez-vous que cette syntaxe est correcte pour l'écriture du résultat extrait dans un fichier (bien sûr ici il n'a subi aucun traitement, le traitement en question est le sujet de ma deuxième question, c'est juste pour voir si j'arrive à écrie sans problème).

Cette partie de code est à mettre juste après :

Code :

fichierEntree.close();

Code :

ofstream fichierSortie("C:\\sortie.txt",ios::in);
if (!fichierSortie)
{
cerr << "Création du fichier impossible" << endl;
exit(1);
}
for (int i=1;i<=dim;i++)
{
for (int j=1;j<=i;j++)
{
double vals;
vals << fichierSortie ;
}
}
fichierSortie.close();

Normalement ca devarit marché mais il y a un problème sur l'opérande << que je ne comprend pas!

Ma deuxième question porte plutôt sur le une fonction issu d'une bibliothèque mathématique la fonction est la suivante :

Code :

ReturnMatrix simind(const int& nbre_noms)
{
//créer une matrice colonne epsilon
ColumnVector epsilon;
//redimensionne la matrice espsilon avec la dimension correpondante
//au nombre de titre
epsilon.ReSize(nbre_noms);
//générateur de loie normale
for (int i=1; i<=nbre_noms; ++i)
{
epsilon(i)=(double)normal_random();
}
epsilon.Release();
return epsilon;
}

Elle retourne un vecteur colonnes dont les composantes sont des variables normales générés aléatoirment. Pourtant quand j'applique cette fonction à mon code je trouve qu'elle génère toujours les mêmes variables, ce que je ne comprend pas :
voici mon code intégrale :

Code :

//fichier main principal
#include "essai.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <string>
using namespace std;
int main()
{
ifstream fichierEntree("c:\\test1.txt", ios::in);
int dim;
fichierEntree >> dim;
ArrayLengthSpecifier dimension(dim);
SymmetricMatrix covariance(dimension);
/* lecture de la matrice */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=i ; ++j )
{
double val;
fichierEntree >> val;
covariance(i,j) = val;
}
}
fichierEntree.close();
//décompose la matrice de covariance suivant cholesky
Matrix C;
C = chol(covariance);
//génère un vecteur colonne de variable normale
ColumnVector E;
E = simind (dim);
ofstream fichierSortie("C:\\sortie.txt",ios::in);
if (!fichierSortie)
{
cerr << "Création du fichier impossible" << endl;
exit(1);
}
for (int i=1;i<=dim;i++)
{
for (int j=1;j<=i;j++)
{
double vals;
vals << fichierSortie ;
}
}
fichierSortie.close();

Merci d'avance pour votre aide

Message cité 2 fois

franceso

Pour la question 1 :

iamora a écrit :

Code :

ofstream fichierSortie("C:\\sortie.txt",ios::in); // si tu veux écrire dans le fichier, c'est ios::out qu'il faut mettre
if (!fichierSortie)
{
cerr << "Création du fichier impossible" << endl;
exit(1);
}
for (int i=1;i<=dim;i++)
{
for (int j=1;j<=i;j++)
{
double vals; // si tu n'initialise pas vals, tu vas écrire n'importe quoi dans ton fichier
vals << fichierSortie ; // c'est dans l'autre sens qu'on écrit : monFlux << maValeurAEcrireDansLeFlux
// comme dans : std::cout << "hello world\n"
}
}
fichierSortie.close();

Pour la question 2:
que fait la fonction normal_random() ? D'où sort elle (une bibliothèque, ou bien c'est toi qui l'as programmée) ? Quand tu dis qu'elle génère toujours la même chose, a veut dire que tu obtiens toujours le même vecteur, ou bien que toutes les composantes de ton vecteur sont égales ? Es-tu sûr de l'utiliser correctement ? habituellement, les fonctions qui génèrent des nombres pseudo-aléatoires doivent être initialisées avec une graine (sinon elles renvoient toujours la même suite de nombres)

Message édité par franceso le 09-08-2006 à 15:43:10

---------------
TriScale innov

_darkalt3_

Proctopathe

iamora a écrit :

oui ok c'est l'espace entre les sorties je m'en suis aperçu quad j'ai testé je suis débutant en C++, alors excusez mon manque d'expérience.

Tout les débutants sont excusés à partir du moment où qu'il mettent les mains dans le cambouis, la documentation et google :jap:

---------------
Töp of the plöp

iamora

Tout d'abord merci pour la présicion de tes réponses.

Pour la question1 : Il me manque encore la notion d'initialiser une variable, que cela veut-il dire ?, et pourquoi ne l'a t-on pas fait lors de la lecture du fichier

Pour la question 2:

que fait la fonction normal_random() ? D'où sort elle (une bibliothèque, ou bien c'est toi qui l'as programmée) ?

La fonction normal_random() génére aléatoirement une valeur d'une variable aléatoire. Elle est définie ci-dessus et provient de la bibliothèque newmat

Code :

double normal_random()
{
return cndev(genrand_real1());
}

genrand_real1 est défini quant à lui ci-dessus

Code :

double genrand_real1(void)
{
return genrand_int32()*(1.0/4294967295.0);
/* divided by 2^32-1 */
}

genrand_int32 est défini quant à lui ci-dessus:

Code :

unsigned long genrand_int32(void)
{
unsigned long y;
static unsigned long mag01[2]={0x0UL, MATRIX_A};
/* mag01[x] = x * MATRIX_A for x=0,1 */
if (mti >= N) { /* generate N words at one time */
int kk;
if (mti == N+1) /* if init_genrand() has not been called, */
init_genrand(5489UL); /* a default initial seed is used */
for (kk=0;kk<N-M;kk++) {
y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1UL];
}
for (;kk<N-1;kk++) {
y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[y & 0x1UL];
}
y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1UL];
mti = 0;
}
y = mt[mti++];
/* Tempering */
y ^= (y >> 11);
y ^= (y << 7) & 0x9d2c5680UL;
y ^= (y << 15) & 0xefc60000UL;
y ^= (y >> 18);
return y;
}

Je crois que ces fonctions partent d'un générateur de nombre aléatoires entre 0 et 1 puis grâce à la fonction de répartition réciproque de la loi normale il permette de générer une variable aléatoire.

Quand tu dis qu'elle génère toujours la même chose, a veut dire que tu obtiens toujours le même vecteur, ou bien que toutes les composantes de ton vecteur sont égales ? Es-tu sûr de l'utiliser correctement ? habituellement, les fonctions qui génèrent des nombres pseudo-aléatoires doivent être initialisées avec une graine (sinon elles renvoient toujours la même suite de nombres)[/quotemsg]

Je veux dire par là que la composante 1 de mon vecteur est toujours la même quand je compile un nouvelle fois ce qui ne devrait pas être le cas normalement. Tout à fait d'accord au sujet de la graine j'ai pensé pareil il faut qu'elle soit intégrer dans le fonction.
Je sais qu'en VB c'est Randomize la graine mais en C++ je n'ai aucune idée???
J'espère que mon problème est un peu plus clair???

Message cité 1 fois
Message édité par iamora le 09-08-2006 à 16:14:34

franceso

iamora a écrit :

Pour la question1 : Il me manque encore la notion d'initialiser une variable, que cela veut-il dire ?, et pourquoi ne l'a t-on pas fait lors de la lecture du fichier

lorsque tu écris

Code :

double vals

tu es en train de dire à ton compilateur : "je vais utiliser à partir de maintenant une variable ne nom 'vals', de type double". Rien de plus, rien de moins. En particulier, la valeur de la variable vals est indéfinie. Lorsque à la ligne suivante, tu écris :

Code :

fichierSortie << vals;

tu es en train d'écrire dans le fichier la valeur de vals, à savoir n'importe quoi.

Il faut que tu écrives dans la variable avant de l'utiliser. Par exemple:

Code :

for( int i ...)
{
for( int j ... )
{
double val;
val = covariance(i,j);
fichierSortie << val;
}
}

Ou plus simplement, puisque ici tu n'as pas besoin d'une variable intermédiaire vals :

Code :

for( int i ...)
{
for( int j ... )
{
fichierSortie << covariance(i,j);
}
}

iamora a écrit :

Pour la question 2:
Je veux dire par là que la composante 1 de mon vecteur est toujours la même quand je compile un nouvelle fois ce qui ne devrait pas être le cas normalement. Tout à fait d'accord au sujet de la graine j'ai pensé pareil il faut qu'elle soit intégrer dans le fonction.
Je sais qu'en VB c'est Randomize la graine mais en C++ je n'ai aucune idée???

D'après les commentaires de la fonction genrand_int32, l'initialisation de la graine se fait avec la fonction init_genrand.
Il n'est pas interdit de regarder la documentation de ta bibliothèque non plus. C'est là que tu trouveras toutes les informations utiles pour utiliser les fonctions qu'elle définit

---------------
TriScale innov

iamora

Merci pour tes réponses efficaces!

Publicité

iamora

En fait ce que je t'ai envoyé sur la fonction normal_random n'est pas tiré sur une bibliothèque Newmat, je viens de m'en apercevoir, je ne connais pas tellement sa proivenance. En tout cas il a l'aire de foinctionner correctement, je pense que c'est l'un de ces générateurs de nombres utilisés pour les simulations de Monté Carlo!
En tout cas voilà le code

Code :

#include "stdafx.h"
#ifndef PI
#define PI 3.141592653589793238462643
#endif
double n(const double& z)
{
return (1.0/sqrt(2.0*PI))*exp(-0.5*z*z);
}
double N(const double& z) {
if (z > 6.0)
return 1.0; // this guards against overflow
if (z < -6.0)
return 0.0;
double b1 = 0.31938153;
double b2 = -0.356563782;
double b3 = 1.781477937;
double b4 = -1.821255978;
double b5 = 1.330274429;
double p = 0.2316419;
double c2 = 0.3989423;
double a=fabs(z);
double t = 1.0/(1.0+a*p);
double b = c2*exp((-z)*(z/2.0));
double n = ((((b5*t+b4)*t+b3)*t+b2)*t+b1)*t;
n = 1.0-b*n;
if ( z < 0.0 ) n = 1.0 - n;
return n;
}
double random_0_1(void)
{
return double(rand())/double(RAND_MAX);
}
double box_muller(void)
{
double U1, U2, V1, V2;
double S=2;
while (S>=1)
{
U1 = random_0_1();
U2 = random_0_1();
V1 = 2.0*U1-1.0;
V2 = 2.0*U2-1.0;
S = (V1*V1)+(V2*V2);
}
double X1=V1*sqrt((-2.0*log(S))/S);
return X1;
}
double normal_random()
{
return cndev(genrand_real1());
}
double lognormal_random(const double& S,const double& r,const double& sigma,
const double& time)
{
double R = (r-0.5 * (sigma*sigma) )*time;
double SD = sigma * sqrt(time);
return S * exp(R + SD * normal_random());
}
#define N 624
#define M 397
#define MATRIX_A 0x9908b0dfUL
#define UPPER_MASK 0x80000000UL
#define LOWER_MASK 0x7fffffffUL
static unsigned long mt[N];
static int mti=N+1;
void init_genrand(unsigned long s)
{
mt[0]= s & 0xffffffffUL;
for (mti=1; mti<N; mti++) {
mt[mti] =
(1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti);
mt[mti] &= 0xffffffffUL;
}
}
void init_by_array(unsigned long init_key[], int key_length)
{
int i, j, k;
init_genrand(19650218UL);
i=1; j=0;
k = (N>key_length ? N : key_length);
for (; k; k--) {
mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525UL))
+ init_key[j] + j; /* non linear */
mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
i++; j++;
if (i>=N) { mt[0] = mt[N-1]; i=1; }
if (j>=key_length) j=0;
}
for (k=N-1; k; k--) {
mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941UL))
- i; /* non linear */
mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
i++;
if (i>=N) { mt[0] = mt[N-1]; i=1; }
}
mt[0] = 0x80000000UL; /* MSB is 1; assuring non-zero initial array */
}
/* generates a random number on [0,0xffffffff]-interval */
unsigned long genrand_int32(void)
{
unsigned long y;
static unsigned long mag01[2]={0x0UL, MATRIX_A};
/* mag01[x] = x * MATRIX_A for x=0,1 */
if (mti >= N) { /* generate N words at one time */
int kk;
if (mti == N+1) /* if init_genrand() has not been called, */
init_genrand(5489UL); /* a default initial seed is used */
for (kk=0;kk<N-M;kk++) {
y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1UL];
}
for (;kk<N-1;kk++) {
y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[y & 0x1UL];
}
y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1UL];
mti = 0;
}
y = mt[mti++];
/* Tempering */
y ^= (y >> 11);
y ^= (y << 7) & 0x9d2c5680UL;
y ^= (y << 15) & 0xefc60000UL;
y ^= (y >> 18);
return y;
}
/* generates a random number on [0,1]-real-interval */
double genrand_real1(void)
{
return genrand_int32()*(1.0/4294967295.0);
/* divided by 2^32-1 */
}
double cndev(double u)
{
static double a[4]={2.50662823884,-18.61500062529,41.39119773534,-25.44106049637};
static double b[4]={-8.47351093090,23.08336743743,-21.06224101826,3.13082909833};
static double c[9]={0.3374754822726147,0.9761690190917186,0.1607979714918209,0.0276438810333863,
0.0038405729373609,0.0003951896511919,0.0000321767881768,0.0000002888167364,0.0000003960315187};
double x,r;
x=u-0.5;
if (fabs(x)<0.42)
{
r=x*x;
r=x*(((a[3]*r+a[2])*r+a[1])*r+a[0])/((((b[3]*r+b[2])*r+b[1])*r+b[0])*r+1.0);
return(r);
}
r=u;
if(x>0.0)
r=1.0-u;
r=log(-log(r));
r=c[0]+r*(c[1]+r*(c[2]+r*(c[3]+r*(c[4]+r*(c[5]+r*(c[6]+r*(c[7]+r*c[8])))))));
if(x<0.0)
r=-r;
return(r);
}

Normalement à chaque fois que je compile je dois tomber sur une nouvelle valeur, mais là c'a à l'aire de garder la même valeur. J'ai remarqué par contre que si je mettai pour autant un :
init_genrand(n) avec n un nombre que je chose à chaque fois avant que je compile, cela change le résultat!

Code :

ReturnMatrix simind(const int& nbre_noms)
{
//créer une matrice colonne epsilon
ColumnVector epsilon;
//redimensionne la matrice espsilon avec la dimension correpondante
//au nombre de titre
epsilon.ReSize(nbre_noms);
//générateur de loie normale
for (int i=1; i<=nbre_noms; ++i)
{
init_genrand(n)
epsilon(i)= normal_random();
}

Décidement, c'est bien compliqué le C++

Amonchakai

Les méthodes aléatoires c'est pas si compliqué : ce sont en fait des suite récurentes du type

Un+1 = (Un^6) % un grand chiffre premier;

enfin mon exemple est peut être (sûrement ) bidon Mais l'idée c'est que c'est une suite. Donc tu doit d'initialiser a partir d'une certaine valeur et ensuite tu obtient des nombres pseudo aléatoire. C'est pour cela que tu as toujours les même nombres : tu part du même premier terme et les autres sont calculé par récurrence a partir du premier...

L'idée si tu veut avoir des résultats différents a chaque exécution serait de l'initialiser a partir d'un chiffre obtenu a partir de l'horloge de l'ordi. Comme ça a chaque fois que tu exécute ton prog tu a un point de départ différent pour ta suite et donc des résultats différents.

Message cité 1 fois

iamora

Amonchakai a écrit :

Merci pour ces explications,C'est plus clair maintenant.
Jstement cette fonction qui me permetterait de récuperer l'horloge de l'ordi, je pensais que c'était la fonction init_genrand(n) ais comme elle prend en argument un entier c pas possible.
Une telle fonction doit être sans argument apparamment, non?
Sauras tu quelle fonction utiliser?
Peut être qu'elle figure déjà dans le code que j'ai intégré dans mon dernier message?

franceso

c'est justement le n de init_genrand(n) qu'on choisit généralement sur la base de l'horloge de l'ordi.
En général, j'utilise la fonction C time().

man 2 time

Message cité 1 fois

---------------
TriScale innov

iamora

franceso a écrit :

c'est justement le n de init_genrand(n) qu'on choisit généralement sur la base de l'horloge de l'ordi.
En général, j'utilise la fonction C time().

Merci franseco pour ton aide, ça marche maintenant j'utiliser pour ça la fonction time() que tu m'as montré

Code :

init_genrand(time(NULL));

En fait c'est quoi le site que tu m'as montré ???, il est un peu bizarre à quoi sert la photo et c'est quoi man??

man 2 time

Message cité 1 fois

_darkalt3_

Proctopathe

"man" c'est pour avoir la doc de quelque chose dans les *nix

franceso

iamora a écrit :

En fait c'est quoi le site que tu m'as montré ???, il est un peu bizarre à quoi sert la photo et c'est quoi man??

Bah la photo, je suppose que c'est pour faire joli :spamafote:

"man", c'est l'abbréviation de "manuel". Il s'agit au départ d'une commande présente dans les unixoides pour donner de la documentation sur les fonctions C et les commandes que tu peux taper dans un terminal. Les pages de man sont toutes écrites de la même façon, et constituent une source centrale de documentation de référence. La plupart des pages de man sont aussi disponible sur internet sur le genre de site que je t'ai montré (donc la prochaine fois que tu te demandes comment utiliser une fonction C, tu peux essayer de taper 'man printf' (par exemple) dans google pour voir ce qui sort).

---------------
TriScale innov

skelter

'man 3 time' plutot ?

Message cité 1 fois

franceso

skelter a écrit :

'man 3 time' plutot ?

> man 3 time
No entry for time in section 3 of the manual

---------------
TriScale innov

iamora

OK!
Merci pour ces infos les gars, ça me sera surement utile par la suite.
L'opérateur << est certes puissant dans le code que je vous ai transmis et que tu as corrigé franseco, je voudrais faire une petit changement d'affichage, mais je me rends compte que c'est impossibel :

Code :

/* affichage pour tester */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
std::cout << setw(10) << covariance(i,j);
}
std::cout << std::endl;
}

C'était pour afficher ma matraice de covariance, seulement moi je ne voudrais pas d'espace au début de chaque ligne, alors je me suis dis en faisant ça :

Code :

std::cout << covariance(i,j) << setw(10);

Mais non ça ne marche pas, il le fait bien pour la première ligne mais pas pour les autres lignes qui suivent puisq'ils laissent un espace au début de chaque nouvelle ligne??

Peut être qu'avec l'opérateur << de sortie c'est pas possible de réaliser ce que je veux non?

Message cité 1 fois

franceso

iamora a écrit :

Peut être qu'avec l'opérateur << de sortie c'est pas possible de réaliser ce que je veux non?

Si si, c'est possible !
Essaie de comprendre le fonctionnement des manipulateurs. En particulier, essaie de comprendre pourquoi tu as un espace en début de chaque ligne, et pourquoi ta solution ne marche pas. Si tu n'y arrives pas, on te proposera une solution, mais il est important que tu arrives à trouver les informations par toi même.

EDIT: Cf le post d'Elmoricq sur les histoires de poissons et de canne à pêche

Message cité 1 fois
Message édité par franceso le 10-08-2006 à 11:16:41

---------------
TriScale innov

iamora

franceso a écrit :

Si si, c'est possible !
Essaie de comprendre le fonctionnement des manipulateurs. En particulier, essaie de comprendre pourquoi tu as un espace en début de chaque ligne, et pourquoi ta solution ne marche pas. Si tu n'y arrives pas, on te proposera une solution, mais il est important que tu arrives à trouver les informations par toi même.

Pas mal l'histoire,
Mais je dois dire que tu as bien raison, j'essaierai de trouver seul !

EDIT: Cf le post d'Elmoricq sur les histoires de poissons et de canne à pêche

iamora

Code :

/* affichage pour tester la matrice de cholesky */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
if (j==1)
std::cout << C(i,j);
else
{
std::cout << setw(12) << C(i,j);
}
}
std::cout << std::endl;
}

C'est mieux comme ça non?

franceso

C'est effectivement l'une des solutions. Bravo !

Message cité 1 fois

---------------
TriScale innov

iamora

franceso a écrit :

C'est effectivement l'une des solutions. Bravo !

Je pense que l'énigme suivante est plus diffcile et malgrès que je me suis planché la-dessus assez longtemps je bloque, je pense que j'ai besoin d'aide cette fois.
Voilà le but de mon programme est le suivant :

1) récupérer une matrice de variance-covariance depuis un fichier .txt (étape faite et testée)

Le fichier .txt a la forme suivante:

Code :

6
3.06E-01
1.83E-02 4.42E-02
6.33E-03 2.98E-03 3.11E-02
1.88E-04 1.88E-04 1.88E-04 6.25E-04
1.29E-02 1.36E-02 -5.03E-04 1.88E-04 9.51E-02
4.49E-03 2.91E-02 6.50E-03 1.88E-04 3.05E-02 2.27E-01

Code :

//fichier main principal
#include "essai.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <string>
using namespace std;
int main()
{
ifstream fichierEntree("D:\\Documents and Settings\\AMOR\\Desktop\\VaRaction\\covariance.txt", ios::in);
int dim;
fichierEntree >> dim;
ArrayLengthSpecifier dimension(dim);
SymmetricMatrix covariance(dimension);
/* lecture de la matrice */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=i ; ++j )
{
double val;
fichierEntree >> val;
covariance(i,j) = val;
}
}
fichierEntree.close();

2) Traiter la matrice récupérée, cette étape se fait sous deux sous-étapes :

2-1) Appliquer une fonction nommée chol (cette fonction réalise la décomposition de Cholesky) Cette fonction c'est moi qui l'ai codée et elle est en annexe en à la fin de mon message (étape faite et testée)

Code :

//décompose la matrice de covariance suivant cholesky
int nbre_simul;
Matrix C;
C = chol(covariance);

2-2) étape ou je bloque:

2-2-1) Prendre la matrice transformée par la fonction Chol puis la multiplier par un vecteur de variable aléatoire normale généré par la fonction simind, dont mainteant tu connais la fonction:
simind : retourne à partir d'un générateur de nombre entre 0 et 1, la valeur correspondante à une variable normale

(étape ok)

2-2-2) Mon but dans cet étape que je n'arrive pas à réaliser construire dans une matrice de dimension (nombre_simul, dim_vecteur généré) ou bien entendu la ligne i représente le produit de la matrice issue de la fonction chol par un vecteur dont les composantes sont générés aléatoirement.
Bien sur chaque ligne i représente une nouvelle simulation .

voici ce que j'ai codé et qui me semble être un début mais ça bloque pour des raisons que je ne comprend pas !

Code :

//décompose la matrice de covariance suivant cholesky
int nbre_simul;
Matrix C;
C = chol(covariance);
//saisie du nombre de simulation
cout << "Veuillez saisir le nombre de simulations souhaitees : ";
cin >> nbre_simul;
// définie la matrice de sortie
Matrix simul_epsilon (nbre_simul,dim);
for( int i=1 ; i<=nbre_simul ; i++ )
{
for( int j=0 ; j=dim ; ++j )
{
//génère un vecteur colonne de variable normale
ColumnVector E;
E = simind (dim);
std::cout << E(1) << std::endl;
//multiplie la matrice de cholesky par le vecteur de varibles normales (simulation d'un vecteur normale)
ColumnVector epsilon_correl;
epsilon_correl = C * E;
std::cout << epsilon_correl(0) << std::endl;
//std::cout << epsilon_correl(2) << std::endl;
simul_epsilon(i,j) = epsilon_correl(j);
}
}
/* affichage pour tester la matrice de cholesky */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
if (j==1)
std::cout << C(i,j);
else
{
std::cout << " " << C(i,j);
}
}
std::cout << std::endl;
}
/* affichage pour tester la matrice simul_epsilon */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
if (j==1)
std::cout << simul_epsilon(i,j);
else
{
std::cout << " " << simul_epsilon(i,j);
}
}
std::cout << std::endl;
}
/* affichage pour tester la génération d'une variable normale*/
//std::cout << E(1) << std::endl;
/* affichage pour tester la matrice de epsilon correl */
//std::cout << epsilon_correl(1) << std::endl;
//écriture de la matrice de Cholesky transformée dans un fichier .txt
ofstream fichierSortie("C:\\Matrice de cholesky transformée.txt",ios::out);
if (!fichierSortie)
{
cerr << "Création du fichier impossible" << endl;
exit(1);
}
for (int i=1;i<=dim;i++)
{
for (int j=1;j<=i;j++)
{
if (j==1)
fichierSortie << C(i,j);
else
{
fichierSortie << " " << C(i,j);
}
}
fichierSortie << endl;
}
return 0;
}

code pour la fonction simind et chol

Code :

// Il s'agit du fichier DLL principal.
#include "stdafx.h"
#include "CorrelSimul.h"
ReturnMatrix chol(const SymmetricMatrix& correlation)
{
//renvoie la dimension de la matrice de corrélation
ArrayLengthSpecifier dimension(correlation.Nrows());
int dimension1 = correlation.Nrows();
//initialise le test de défini positivité à faux
bool test=false;
//diagonalise la matrice de corrélation
DiagonalMatrix D(dimension);
Matrix V(dimension1,dimension1);
LowerTriangularMatrix L(dimension);
Matrix resultat;
//renvoie les valeurs propres
EigenValues(correlation,D,V);
//
if (D(1)<=0) test=true;
if (test=false)
{
L = Cholesky(correlation);
resultat=L;
}
else
{
for (int i = 1; i <= dimension1; i++)
{
D(i)=sqrt(max(D(i),0)); /* on transforme D en D+ puis en (D+)^1/2
pour prendre la racine carrée */
}
resultat = V*D*V.t();
}
D.Release();
V.Release();
L.Release();
return resultat;
}
ReturnMatrix simind(const int& nbre_noms)
{
//créer une matrice colonne epsilon
ColumnVector epsilon;
//redimensionne la matrice espsilon avec la dimension correpondante
//au nombre de titre
epsilon.ReSize(nbre_noms);
//générateur de loie normale
for (int i=1; i<=nbre_noms; ++i)
{
init_genrand(time(NULL));
epsilon(i)= normal_random();
}
epsilon.Release();
return epsilon;
}

Message édité par iamora le 10-08-2006 à 15:47:48

franceso

plusieurs remarques :

1- il faut initialiser la graine de ton générateur alatoire une seule fois au début du programme (par exemple dans la fonction main()).

2- dans chol(), je ne vois pas trop l'intéret de la variable booléenne temporaire test

3- tu devrais faire une fonction qui affiche une matrice, plutôt que de faire du copier/coller de code

4- tel quel, ton code ne doit sûrement pas compiler : la variable E est locale à la boucle du haut, mais tu demandes de l'afficher en dehors de cette boucle

5- tu as à mon avis un gros problème dans ta première double boucle : tel que tu es en train de le faire, tu génères un nouveau vecteur aléatoire pour chaque couple (i,j), ce qui veut dire qu'au final tu auras généré dim x nbre_simul vecteurs alatoires. Ceci ne semble pas cohérent avec le fait que chaque vecteur aléatoire est censé reprsenter une simulation.
D'après ce que j'ai compris (arrête moi si je me trompe), C est la transformation de Choleski de la matrice correlation. Donc C est de dimension dim.
Une 'simulation' est modélisée par un produit entre la matrice C et un vecteur E de dim variables alatoires indpendantes et identiquement distribues selon une loi normale. Le résultat d'une simulation epsilon_correl est donc un vecteur de dimension dim.
Tu ranges tes nbre_simul résultats de simulation en lignes dans une matrice simul_epsilon, qui est donc de dimension (nbre_simul, dim)

Dans ce cas, ta première double boucle devrait s'écrire :

Code :

ColumnVector E;
ColumnVector epsilon_correl;
for(int i=1 ; i<=nbre_simul ; ++i )
{
E = imind(dim);
epsilon_correl = C * E;
for(int j=1 ; j<=dim ; ++j )
{
simul_epsilon(i,j) = epsilon_correl(j);
}
}

6- tu as un problème dans les dimensions de tes matrices utilisées pour l'affichage

---------------
TriScale innov

iamora

[quotemsg=1423777,26,383694]plusieurs remarques :

Merci, encore une fois tes indications sont trés clairs, et j'ajouterai que tu connais trés bien le sujet.
Pour le point 5, de ton précédent message, je n'avais pas fait attention effectivement mais j'ai bien redimensionner ma matrice.
Après l'application de tes conseils voilà le code que j'obtiens et qui marche!

Code :

//fichier main principal
#include "essai.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <string>
using namespace std;
int main()
{
//initialise la graine du générateur de variable normale
init_genrand(time(NULL));
ifstream fichierEntree("D:\\Documents and Settings\\AMOR\\Desktop\\VaRaction\\covariance.txt", ios::in);
int dim;
fichierEntree >> dim;
ArrayLengthSpecifier dimension(dim);
SymmetricMatrix covariance(dimension);
/* lecture de la matrice */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=i ; ++j )
{
double val;
fichierEntree >> val;
covariance(i,j) = val;
}
}
fichierEntree.close();
//décompose la matrice de covariance suivant cholesky
Matrix C;
C = chol(covariance);
int nbre_simul;
ColumnVector E;
ColumnVector epsilon_correl;
//saisie du nombre de simulation
cout << "Veuillez saisir le nombre de simulations souhaitees : ";
cin >> nbre_simul;
// définie la matrice de sortie
Matrix simul_epsilon (nbre_simul,dim);
for(int i=1 ; i<=nbre_simul ; ++i )
{
E = simind(dim);
epsilon_correl = C * E;
for(int j=1 ; j<=dim ; ++j )
{
simul_epsilon(i,j) = epsilon_correl(j);
}
}
std::cout << simul_epsilon(4,4) << std::endl;
/* affichage pour tester la matrice de cholesky */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
if (j==1)
std::cout << C(i,j);
else
{
std::cout << " " << C(i,j);
}
}
std::cout << std::endl;
}
/* affichage pour tester la matrice simul_epsilon */
for( int i=1 ; i<=nbre_simul ; ++i )
{
for( int j=1 ; j<=dim ; ++j )
{
if (j==1)
std::cout << simul_epsilon(i,j);
else
{
std::cout << " " << simul_epsilon(i,j);
}
}
std::cout << std::endl;
}
//écriture de la matrice de Cholesky transformée dans un fichier .txt
ofstream fichierSortie("C:\\Matrice de cholesky transformée.txt",ios::out);
if (!fichierSortie)
{
cerr << "Création du fichier impossible" << endl;
exit(1);
}
for (int i=1;i<=dim;i++)
{
for (int j=1;j<=i;j++)
{
if (j==1)
fichierSortie << C(i,j);
else
{
fichierSortie << " " << C(i,j);
}
}
fichierSortie << endl;
}
return 0;
}

Cependant comme tu commence à me donner goût à la rigueur en C++ je suis tout à fait d'accord avec toi pour representer les choses dans une fonction comme réalisée pour chol et simind.
J'arrange cette fonction et si tu n'y voi pas d'incovénient je te la montre pour me dire si c'est ok???

Message édité par iamora le 10-08-2006 à 17:39:14

iamora

Voilà la fonction que j'ai composé !
Elle s'apelle simcor, je te fournis tout mon fichier ou sont définies mes fonctions

Code :

// Il s'agit du fichier DLL principal.
#include "stdafx.h"
#include "CorrelSimul.h"
ReturnMatrix chol(const SymmetricMatrix& correlation)
{
//renvoie la dimension de la matrice de corrélation
ArrayLengthSpecifier dimension(correlation.Nrows());
int dimension1 = correlation.Nrows();
//initialise le test de défini positivité à faux
bool test=false;
//diagonalise la matrice de corrélation
DiagonalMatrix D(dimension);
Matrix V(dimension1,dimension1);
LowerTriangularMatrix L(dimension);
Matrix resultat;
//renvoie les valeurs propres
EigenValues(correlation,D,V);
//
if (D(1)<=0) test=true;
if (test=false)
{
L = Cholesky(correlation);
resultat=L;
}
else
{
for (int i = 1; i <= dimension1; i++)
{
D(i)=sqrt(max(D(i),0)); /* on transforme D en D+ puis en (D+)^1/2
pour prendre la racine carrée */
}
resultat = V*D*V.t();
}
D.Release();
V.Release();
L.Release();
return resultat;
}
ReturnMatrix simind(const int& nbre_noms)
{
//créer une matrice colonne epsilon
ColumnVector epsilon;
//redimensionne la matrice espsilon avec la dimension correpondante
//au nombre de titre
epsilon.ReSize(nbre_noms);
//générateur de loie normale
for (int i=1; i<=nbre_noms; ++i)
{
epsilon(i)= normal_random();
}
epsilon.Release();
return epsilon;
}
ReturnMatrix simcor(const int& nbre_simul, const int& dim, const SymmetricMatrix& covariance)
{
//décompose la matrice de covariance suivant cholesky
Matrix C;
C = chol(covariance);
//définie le vecteur de variables aléatoires normales
ColumnVector E;
//définie le vecteur issue de la multiplication de E et C
ColumnVector epsilon_correl;
// définie la matrice de sortie
Matrix simul_epsilon (nbre_simul,dim);
for(int i=1 ; i<=nbre_simul ; ++i )
{
E = simind(dim);
epsilon_correl = C * E;
for(int j=1 ; j<=dim ; ++j )
{
simul_epsilon(i,j) = epsilon_correl(j);
}
}
C.Release();
E.Release();
epsilon_correl.Release();
simul_epsilon.Release();
return simul_epsilon;
}

Puis voilà mon fichier principal ou j'utilise la fonction

Code :

//fichier main principal
#include "essai.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <string>
using namespace std;
int main()
{
//initialise la graine du générateur de variable normale
init_genrand(time(NULL));
ifstream fichierEntree("D:\\Documents and Settings\\AMOR\\Desktop\\VaRaction\\covariance.txt", ios::in);
int dim;
fichierEntree >> dim;
ArrayLengthSpecifier dimension(dim);
SymmetricMatrix covariance(dimension);
/* lecture de la matrice */
for( int i=1 ; i<=dim ; ++i )
{
for( int j=1 ; j<=i ; ++j )
{
double val;
fichierEntree >> val;
covariance(i,j) = val;
}
}
fichierEntree.close();
//saisie du nombre de simulation
int x;
cout << "Veuillez saisir le nombre de simulations souhaitees : ";
cin >> x;
//définie la matrice de sortie
Matrix sortie;
sortie = simcor(x, dim, covariance);
std::cout << sortie(4,4) << std::endl;
return 0;
}

Normalement ça devrait marcher, mais j'ai un message d'erreur le compilateur me dit que

Code :

essai.cpp(43) : error C2660: 'simcor' : la fonction ne prend pas 3 arguments

Ce que je trouve bizarre puisque j'ai bien définie une fonction à trois arguments???

Message édité par iamora le 11-08-2006 à 10:37:39

franceso

Oui, ce que tu as écrit me semble pas mal. Tu devrais peut-être encapsuler le code d'input de la matrice dans une fonction : ça rendrait les choses plus claires et plus llisibles.

Concernant ton erreur de compilation, tu n'as pas de messages de warnings autour de ton erreur ? As-tu défini les prototypes de tes fonctions dans un fichier d'en-tête ? As-tu inclus le fichier d'en-tête en question dans ton programme principal ?

Message cité 1 fois

---------------
TriScale innov

iamora

franceso a écrit :

Oui, ce que tu as écrit me semble pas mal.
Bonjour,

Tu devrais peut-être encapsuler le code d'input de la matrice dans une fonction : ça rendrait les choses plus claires et plus llisibles.

Je ne sais pas ce que veux dire le terme encapsuler un code (je suis encore débutant, j'apprend sur le tas...)

Concernant ton erreur de compilation, tu n'as pas de messages de warnings autour de ton erreur ?
Non pas de warning, la seule erreure que j'ai est celle que je t'ai envoyé !

As-tu défini les prototypes de tes fonctions dans un fichier d'en-tête ?
Encore une fois et je suis désolén mais je ne comprend pas ce que sont les prototypes de fonction??, est ce les includes???

As-tu inclus le fichier d'en-tête en question dans ton programme principal ?

En fait voici mon fichier d'en-tête précompilé stadfx.h, je ne sais pas si il faut mettre la dedans comme mis en commentaire les fonctions que j'ai codé dans le fichier correl_simul.cpp.
Mais apparament non puisque ça marche sans ??

Code :

// stdafx.h : Fichier Include pour les fichiers Include système standard,
// ou les fichiers Include spécifiques aux projets qui sont utilisés fréquemment,
// et sont rarement modifiés
#define WINVER 0x400
#pragma once
#include <Afx.h>
#include <afxdisp.h>
#include <afxwin.h>
#include <algorithm>
#include <numeric>
#include <fstream>
#include <vector>
#include <time.h>
#include <math.h>
#include <iostream>
#include <tchar.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <vector>
#include "newmat\NEWMATAP.h"
#include "newmat\NEWMAT.h"
#include "newmat\myexcept.h"
using namespace std;
void init_genrand(unsigned long s);
void init_by_array(unsigned long init_key[], int key_length);
unsigned long genrand_int32(void);
double genrand_real1(void);
double random_0_1(void);
double cndev(double u);
double N(const double& z);
double normal_random();
double box_muller();
double n(const double& z);
double lognormal_random();
/*ReturnMatrix cholesky (const SymmetricMatrix& S);*/
//ReturnMatrix simind(const int& nbre_noms)
/*ReturnMatrix simcor(const int& nbre_simul, const int& dim, const SymmetricMatrix& covariance);*/

franceso

Je ne sais pas trop comment marchent les headers précompilés (je suppose que tu es sous VC++ ou quelque chose du genre ?). Essaie peut-être de te documenter dans ce sens la...

Sinon, que contiennent tes headers "essai.h" et "CorrelSimul.h" ?

Message cité 1 fois

---------------
TriScale innov

iamora

Ne fait pas attention au début de ce message,regarde la fin c'est peut être la raison de mon erreur,

franceso a écrit :

OK
Je suis sous Visualstudio.net 2003.
Je suis en mode release et dans les proporiétés de mon projet sous entêtes précompilés je compile avec l'option :
Générer automatiquement (/YX)
Normalement ça ne devrait pas posé mon problème .
Un ami me dit que c'est peut être du au fait que ma fonction est définie avec un ReturnMatrix et que j'initialise sortie à Matrix dans le fichier priniicpal, mais je doute que ça soit la raison puisqu'avant ça marchait pour la fonction chol définie avec ReturnMatrix et renvoyant un resultat dans le fichier principal dans une variable déclarée comme Matrix
Matrix C;
C = chol(covariance);

En tout cas voilà le fichier "essai.h"

Code :

#include "stdafx.h"
#include "CorrelSimul.h"

et le fichier correlsimul.h

Code :

// CorrelSimul.h
//Specifies that the file will be included (opened) only once by the compiler in a build.
//This can reduce build times as the compiler will not open and read the file
//after the first #include of the module.
//#pragma once
//using namespace System
//namespace CorrelSimul
// public __gc class Class1
// {
// TODO : ajoutez ici vos méthodes pour cette classe.
// };
ReturnMatrix chol(const SymmetricMatrix& correlation);
ReturnMatrix simind(const int& nbre_noms);
ReturnMatrix simcor(const int& nbre_noms, const double& rho);
ReturnMatrix spreads(const Matrix& lambda, const ColumnVector& T, const Matrix& r, const ColumnVector& sigma, const ColumnVector& defaut,
const ColumnVector& x, const double& correlation, const int& nbre_noms,const int& no_steps);

Ah oui je ne l'ai pas changé la fdéfinition de ma focntion dans le fichier correlsimul.h, erreur de débutant je la change et je tiens au courant

[code]

Message édité par iamora le 11-08-2006 à 11:26:57

iamora

C'était ça l'erreur, la preuve que je ne maitrise pas du tout encore le c++.
Merci pour tes réflexions ça m'a mis sur la voie

franceso

on est tous passés par là. C'est les erreurs les plus bêtes qui prennent parfois le plus de temps à corriger.

De rien

---------------
TriScale innov

iamora

Franseco, t'es de l'ensta?
Si c'est le cas je peux te demander ton avis sur une question purement technique en mathématiques?
SI ça risque de te gêner t'as qu'à me le dire!

Message cité 1 fois

franceso

iamora a écrit :

Franseco, t'es de l'ensta?
Oui

Si c'est le cas je peux te demander ton avis sur une question purement technique en mathématiques?
Tu peux toujours demander... je ne garantis rien sur la réponse !

SI ça risque de te gêner t'as qu'à me le dire!
Pas de prob

Message cité 1 fois

---------------
TriScale innov

iamora

franceso a écrit :

Donc voilà

Je voudrais comme tu le sais simuler en C++ un vecteur de variables aléatoires normales. Pour cela la méthode conseillée est de générer un vecteur colonne de variables normales, que je noterai [epsilon] puis de le multiplier par la matrice de "Cholesky" qui est une matrice issue de la décomposition de Cholesky de la matrice de variance covariance :

[var-covar]=[chol]*[chol]t

[var-covar]: est la matrice de variance-covariance
[chol]: est la matrice de Cholesky
[chol]t: est la transposée de la matrice de Cholesky

Le problème est que pour que cette décomposition soit juste, il faut que [var-covar] soit définie positive or ce n'est pas toujours mon cas.
Mon objectif est donc de trouver une méthode qui puisse forcer la matrice [var-covar] à devenir une matrice symétrique définie positive.
Pour cela j'ai pensé à diagonaliser [var-covar] puis à remplacer toutes les valeurs propres négatives par des 0 :

[var-covar] = [v]-1*[d]*[v]
ou [d] est une matrice diagonale dont les valeurs>0

ensuite je prends la racine des éléments de [d], donc diag[di]=di^(1/2).
Puis je dis que la matrice de Cholesky est la suivante
chol = [v]-1*[d]^(1/2)*[v].

Je voudrais savoir si mon raisonnemment est correct, car j'ai un doute sur la validité de ma méthode.

Amonchakai

Salut !
Bon j'avoue tout de suite je n'y connait rien en matrice de variance-covaiance et Cholesky. Aussi j'ai été voir Google et j'ai trouvé ça :
http://fr.wikipedia.org/wiki/Matri [...] covariance

je sais pas si ce qui est dit dans Wikipedia est juste mais visiblement y'a un soucis : Ils disent que la matrice de variance-covariance est Symétrique positive. Donc y'a comme un soucis si tu trouve des Valeurs propre négatives... :??:
Ensuite ces valeurs propre (négatives...) tu veut les remplacer par des 0. Déja là tu rend ta matrice positive et PAS définie positive (tes valeurs propre doivent être strictement positive...) donc faudrait les remplacer par un nombre positif != de zéro... Mais moi ce qui me turlupine c'est qu'est-ce qui te donne le doit de changer ta matrice comme ça... Si tu change ses valeurs propre elle n'aurra plus rien a voir avec celle initiale (surtout qu'après tu va la multiplier par les matrices diagonalisantes...) Je sais pas ce que ça va donner sûrement qqch qui n'a rien a voir avec tes résultats initiaux...

Et puis juste par curiosité tu fait comment pour trouver les valeurs propres... Moi je chercherai les racines du polynôme caractéristique... Mais c'est du calcul formel : ça se fait bien a la main mais pour le faire faire a un ordinateur... (avec Maple Ok... mais sinon...)

[edit] J'espère que tu prend pas mal ma remarque, je veut pas répondre méchament c'est juste que je comprend pas bien ce que tu veut faire Et si j'ai dit un bêtise excuse moi

Message cité 1 fois

franceso

Amonchakai a écrit :

je sais pas si ce qui est dit dans Wikipedia est juste mais visiblement y'a un soucis : Ils disent que la matrice de variance-covariance est Symétrique positive. Donc y'a comme un soucis si tu trouve des Valeurs propre négatives... :??:
Ensuite ces valeurs propre (négatives...) tu veut les remplacer par des 0. Déja là tu rend ta matrice positive et PAS définie positive (tes valeurs propre doivent être strictement positive...) donc faudrait les remplacer par un nombre positif != de zéro... Mais moi ce qui me turlupine c'est qu'est-ce qui te donne le doit de changer ta matrice comme ça... Si tu change ses valeurs propre elle n'aurra plus rien a voir avec celle initiale (surtout qu'après tu va la multiplier par les matrices diagonalisantes...) Je sais pas ce que ça va donner sûrement qqch qui n'a rien a voir avec tes résultats initiaux...

+1
c'est ce qui me paraît être la plus grande faille dans le raisonnement. Tu ne peux pas comme ça décider de changer ta matrice pour qu'elle vérifie les propriétés dont tu as besoin.

Il y a aussi un autre problème dans les dernières lignes : à mon avis tu ne peux pas te passer de calculer la décomposition de choleski de ta matrice "transformée". Ton histoire à base de "racine carrée" de la matrice diagonale ne marche pas à mon avis (il suffit de calculer le C x Ct avec l'expression que tu proposes pour voir qu'on ne retombe pas sur la matrice de départ)

Amonchakai a écrit :

Et puis juste par curiosité tu fait comment pour trouver les valeurs propres... Moi je chercherai les racines du polynôme caractéristique... Mais c'est du calcul formel : ça se fait bien a la main mais pour le faire faire a un ordinateur... (avec Maple Ok... mais sinon...)

Il existe plein de méthodes numériques pour déterminer des valeurs approchées des valeurs propres d'une matrice (par exemple l'algorithme de Jacobi pour les matrices symétriques)

Message édité par franceso le 11-08-2006 à 14:00:08

---------------
TriScale innov

Amonchakai

Citation :

Il existe plein de méthodes numériques pour déterminer des valeurs approchées des valeurs propres d'une matrice (par exemple l'algorithme de Jacobi pour les matrices symétriques)

ha tien je connais pas.... Merci pour l'info je vais voir ça pour ma culture....

Publicité

Page : 1 2

Page Suivante

Page Précédente

Haut de page

FORUM HardWare.fr

Programmation

C++

Programme pour lire une matrice triangulaire inférieure à partir

Sujets relatifs
[Résolu] Lire les 10 dernières lignes d'un fichier	Problème sur un flux de lecture à partir d'un fichier .txt
Transformer une matrice à une image	Lire un fichier de properties
recharger une page à partir d'une autre	[Résolu] [C#] Mettre à jour une List à partir d'une autre
Création d'un Word ou Excel à partir d'une page ASP	Inserser et lire fichier word dans base de donnees SQL
[Fortran] Renommer un fichier à partir d'une variable	[OpenGL && FMOD] Lire un mp3
Plus de sujets relatifs à : Programme pour lire une matrice triangulaire inférieure à partir

Page générée en 0.338 secondes