Discussion :

[phyfab]

Bonjour à tous !

Voici mon problème

J'ai une fonction assez compliquée qui dépend de deux paramètres x et t.

Pour l'instant j'ai crée une liste x, j'ai fixé t et j'ai tracé ma fonction. J'aimerais pouvoir tracer ma fonction à différents t sans avoir à créer des t1,t2, ... et recopier ma fonction à chaque fois, j'ai essayé diverses méthodes mais ça ne marche pas

D'autre part je souhaiterais faire un spectrogramme avec t en abcisse,x en ordonnée, et les différentes valeurs de f correspondant à un code couleur. Le problème est que je ne peux pas inverser ma relation pour avoir x en fonction de f. Y a t'il moyen de faire autrement ?

[Jean Dumoncel]

Bonjour,

Pour l'instant j'ai crée une liste x, j'ai fixé t et j'ai tracé ma fonction. J'aimerais pouvoir tracer ma fonction à différents t sans avoir à créer des t1,t2, ... et recopier ma fonction à chaque fois, j'ai essayé diverses méthodes mais ça ne marche pas

Justment, tu peux créer une fonction (function) qui te permettra de ré-utiliser ta fonction sans ré-écrire tout le code.

D'autre part je souhaiterais faire un spectrogramme avec t en abcisse,x en ordonnée, et les différentes valeurs de f correspondant à un code couleur. Le problème est que je ne peux pas inverser ma relation pour avoir x en fonction de f. Y a t'il moyen de faire autrement ?

Pourquoi as-tu besoin de x en fonction de y?

[phyfab]

Ah merci c'est ça que je cherchais !

En fait ma fonction c'est la température en fonction de la position x et du temps T(x,t), je veux faire un spectrogramme de T(x,t), j'aimerais avoir x en ordonnée, t en abcisse, en rouge les zones de plus hautes températures, en bleu les zones de plus basses températures, quelque chose comme ça

un exemple : http://questions2physique.files.word...ectrogram3.gif

[Jean Dumoncel]

Donc il te suffit d'avoir les données x, t et T (que tu es en train de calculer si j'ai bien compris), puis utiliser une fonction comme imagesc, mais cela dépend du format des données. Quand tu auras ces données, pourrais-tu nous les fournir pour faire des tests?

[phyfab]

Ce sont des listes, j'ai:

Code :

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

1
2
x=0:0.001:e;
t=0:10:1000;

et

Code :

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

T=T1-4*((T1-T0)/(pi)).*sin(pi*x/e)*(exp(-D*t*(1*pi/e)^2))+4*((T1-T0)/(3*pi))*sin(3*pi*x/e)*exp(-D*t*(3*pi/e)^2)+4*((T1-T0)/(5*pi))*sin(5*pi*x/e)*exp(-D*t*(5*pi/e)^2)+4*((T1-T0)/(7*pi))*sin(7*pi*x/e)*exp(-D*t*(7*pi/e)^2);

avec :

Code :

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

1
2
3
4
T1=220;
T0=20;
D=0.0004;
e=1;

[Jean Dumoncel]

Du coup, si tu calcules T pour chaque x et t :

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
T1=220;
T0=20;
D=0.0004;
e=1;
x=0:0.001:e;
t=0:10:1000;
 
T = zeros(length(x),length(t));
for k=1:length(x)
    for p=1:length(t)
        T(k,p)=T1-4*((T1-T0)/(pi)).*sin(pi*x(k)/e)*(exp(-D*t(p)*(1*pi/e)^2))+4*((T1-T0)/(3*pi))*sin(3*pi*x(k)/e)*exp(-D*t(p)*(3*pi/e)^2)+4*((T1-T0)/(5*pi))*sin(5*pi*x(k)/e)*exp(-D*t(p)*(5*pi/e)^2)+4*((T1-T0)/(7*pi))*sin(7*pi*x(k)/e)*exp(-D*t(p)*(7*pi/e)^2);
    end
end
 
imagesc(x,t,T)

ou en vectorisant le code (attention à bien utiliser l'opérateur '.*' pour faire des multiplications terme à terme) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
T1=220;
T0=20;
D=0.0004;
e=1;
x=0:0.001:e;
t=0:10:1000;

T = zeros(length(x),length(t));
[tt xx] = meshgrid(t,x);
T=T1-4*((T1-T0)/(pi))*sin(pi*xx/e).*(exp(-D*tt*(1*pi/e)^2))+4*((T1-T0)/(3*pi))*sin(3*pi*xx/e).*exp(-D*tt*(3*pi/e)^2)+4*((T1-T0)/(5*pi))*sin(5*pi*xx/e).*exp(-D*tt*(5*pi/e)^2)+4*((T1-T0)/(7*pi))*sin(7*pi*xx/e).*exp(-D*tt*(7*pi/e)^2);

imagesc(x,t,T)

A adapter selon ce que tu souhaites afficher