Discussion :

[lea6969]

Bonjour,

j'ai un petit problème avec un programme de matlab, je résouds un système de 7 équations différentielles couplées et j'aimerais que la matrice des résultats s'affiche dans un fichier de sortie. Mais le problème est que je n'arrive pas à afficher la matrice comme il faut !!
Voici quelques lignes :

Code :

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

1
2
3
4
[t,y] = ode45(@hote,[To Tmax],[S0,E0,I0,R0,X0,Y0,Z0]); 
fid=fopen('fichiersortie.m','w');
fprintf(fid,'*******matrice obtenue : \n\n');
fprintf(fid,'%1.4f\t %1.4f\t %1.4f\t %1.4f\t %1.4f\t %1.4f\t %1.4f\t %1.4f\t\n',[t,y]);

mais avec ça, j'obtiens comme matrice :

t0 t1 t2 t3 t4 t5 t6 t7
t8 t9 t10......................
....................................
y1(t0) y1(t2) y1(t3)...........
......................................
Y2(t0) y2(t1).................
...................................

Alors que j'aimerais :

t0 y1(t0) y2(t0) y3(t0) y4(t0) y5(t0) y6(t0) y7(t0)
t1 y1(t1) y2(t1) y3(t1) y4(t1) y5(t1) y6(t1) y7(t1)
......................................................................

Merci pour votre aide.

[Jerome Briot]

Bonjour,

fid=fopen('fichiersortie.m','w');

Premièrement, on utilise toujours l'argument t de FOPEN lorsque l'on travaille sur des fichiers ascii (sauf sous Unix)

Code :

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

fid=fopen('fichiersortie.m','wt');

Ensuite, on peut lire dans la documentation de FPRINTF :

Code :

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

1
2
3
4
FPRINTF behaves like ANSI C with certain exceptions and extensions...
    4. FPRINTF is "vectorized" for the case when A is nonscalar. The
       format string is recycled through the elements of A (columnwise)
       until all the elements are used up.

Ce qui donne (par exemple) :

Code :

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

1
2
3
t=1:10
y=rand(10,3)
fprintf('%2d %.4f %.4f %.4f\n',[t(:) y]')