Discussion :

[prince cta]

Salut, chers tous
svp j'ai une préoccupation:

J'ai utilisé la fonction préprogrammé ode15s ([t,Y] = ode15s(@f,[500 10^5],Y0);) pour la résolution d'équations différentielles (de type raide) de la forme dY/dt=f(x,t);
J'ai donc obtenu une matrice Y de dimension 450x450, et j'aimerai juste extraire quelques colonnes significatives (par exemple Y(:,5); Y(:,17), Y(:,39), Y(:,73); Y(:,279) et Y(:,426)) pour inscrire dans un nouveau fichier avec la variable temps t.

Mon problème: j'aimerai écrire dans un même fichier les couples (t, Y(:,i)), i=5, 17, ...,426 de sorte d'obtenir des vecteurs colonnes facilement exploitables.
j'ai utilisé la syntaxe

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fprintf(nomdufichierdesortie, '%f\n', t,  Y(:,5), t, Y(:,17), t, Y(:,39), t, Y(:,73), t, Y(:,279), t, Y(:,426));

mais j'obtiens plutôt un fichier dans lequel les valeurs se suivent et ne sont pas rangées en colonnes comme je le souhaite!

j'ai essayé d'autres syntaxes en changeant même l'extension du fichier de sortie, mais toujours rien. Merci de m'aider!

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
fidCTA4 = fopen('TestCsvCta.txt', 'w');
fidCTA44 = fopen('TestCsvCta_2.csv', 'w');
fidCTA5 = fopen('TestDatCta.csv', 'w');
%fidCTA6 = fopen('TestXlsxCta.xlsx', 'w');
fidCTA6 = fopen('TestXlsxCta.xls', 'w');
fidCTA7 = fopen('TestTxtCta.txt', 'w');
 
fprintf(fidCTA5, '%i\n', t', Y(:,5*n+3));
%fprintf(fidCTA4, '%i\n', t);
%fprintf(fidCTA4, '%f\v', Y(:,5*n+2));
fprintf(fidCTA4, '%f\n', t);
fprintf(fidCTA44, '%f\n %f\n', [t Y(:,5*n+7)]);
fprintf(fidCTA6, '%f\n', t,  Y(:,5), t, Y(:,17), t, Y(:,39), t, Y(:,73), t, Y(:,279), t, Y(:,426));
fprintf(fidCTA7, '%f\v', Y(:,5*n+2), Y(:,5*n+5), Y(:,5*n+7));
fprintf(fidCTA7, '%f\v', t);
fprintf(fidCTA7, '%f\v', Y(:,5*n+2));
fprintf(fidCTA7, '%f\v', t);
fprintf(fidCTA7, '%f\v', Y(:,5*n+7));
 
fclose(fidCTA4);
fclose(fidCTA44);
fclose(fidCTA5);
fclose(fidCTA6);
fclose(fidCTA7);

[Jerome Briot]

Comme ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fprintf(nomdufichierdesortie, '%f %f %f %f %f %f %f\n', t, Y(:,5), Y(:,17), Y(:,39), Y(:,73), Y(:,279), Y(:,426));

[prince cta]

Merci Jerome, ça marche!

Je viens aussi de m'en rendre compte que c'est bloc-notes (mon soft de lecture par défaut) qui n'affiche que des blocs qui se suivent toujours. J'ai donc essayé d'ouvrir avec WordPad et Notepad++, et l'alignement est respecté.

SVP j'ai deux petites préoccupations supplémentaires:

-Lorsque j'utilise l'extension ".csv" pour pouvoir ouvrir le fichier de sortie avec Excel, tous les résultats s'inscrivent dans une même colonne du tableur excel, en respectant toujours la séparation. Je ne peux donc tracer une variable en fonction de l'autre dans excel!

Ma question 1 : Est-ce possible, si oui comment faire pour que les différents blocs s'inscrivent dans différentes colonnes (avec excel)?

-Question 2: Comment insérer au début de chaque bloc le nom de la variable mesurée?

Merci d'avance...

[le fab]

salut

en utilisant le séparateur ; dans ton csv, ca marchera nickel en l'ouvrant dans Excel ou avec xlsread

Fabien

[prince cta]

ok Fabien,
mais en utilisant le ; comme séparateur dans le code matlab, ça marquera "??? Error: File:", Unbalanced or unexpected parenthesis or bracket.

et s'il faut séparer ainsi dans le fichier résultat, ce sera un travail fastidieux vu la taille de la matrice!

[Jerome Briot]

Comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fprintf(nomdufichierdesortie, '%f;%f;%f;%f;%f;%f;%f\n', t, Y(:,5), Y(:,17), Y(:,39), Y(:,73), Y(:,279), Y(:,426));

[prince cta]

dans la format '%f %f %f %f\n', j'ai remplacé \n par \r\n ==> '%f %f %f %f\r\n', et ça semble résoudre mou problème avec Bloc-notes.

J'ai même pensé à écrire un autre petit code qui prend en entrée les résultats du premier code après itérations, pour mieux les manipuler.

Et pour l'insertion des noms des variables, t'as une idée de la procédure?

Armand

[prince cta]

Jerome, pour le code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fprintf(nomdufichierdesortie, '%f;%f;%f;%f;%f;%f;%f\n', t, Y(:,5), Y(:,17), Y(:,39), Y(:,73), Y(:,279), Y(:,426));

ça marche au niveau de l'écriture, j'obtiens bel et bien 7 colonnes, mais à regarder les résultats de plus prêt, je m'en rends compte que toutes les "n" premières lignes des 7 colonnes (sous matrice nx7) dans le fichier résultat ne concernent que la variable t. les n linges suivantes (sous matrice (n+1:2*n)x7) correspondent à la variable Y(:,5), ainsi de suite. Dans chaque colonne du fichier résultat se trouvent donc les valeurs de toutes les 7 variables qui m'intéressent!

or je voudrais plutôt écrire, par exemple dans toute la première colonne du fichier de sortie (fichier résultat), uniquement les valeurs de à la première variable "t", ce qui correspond à la sous matrice nx1 dans le fichier résultat, et aussi écrire, dans toute la deuxième colonne, uniquement la 2ème variable "Y(:,5)" ==> la sous matrice nx2, ..., à la 7ème colonne, la 7ème variable "Y(:,426)" ==> la sous matrice nx7;

Armand

[Jerome Briot]

J'ai effectivement répondu un peu n'importe quoi
Je n'avais pas lu correctement ton premier message

Essaie plutôt ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fprintf(nomdufichierdesortie, '%f;%f;%f;%f;%f;%f;%f\n', [t Y(:,5) Y(:,17) Y(:,39) Y(:,73) Y(:,279) Y(:,426)].');

[prince cta]

Merci, Jerome.

j'ai oublié de le mentionner la dernière fois, mais le nouveau code a très bien marché. Cela me permettra d'évoluer plus vite. Merci.

Armand