Discussion :

[Julien44]

Bonjour,

J'ai un signal temporel e=f(t) où t est une variable symbolique:

Code :

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

1
2
syms t positive
e = 2*t

Je voudrais évaluer cette fonction pour différentes valeur de t afin de la tracer dans un "plot" (je ne peux pas utiliser "ezplot" car je trace d'autres résultats sur la même figure)

Pour évaluer la fonction j'ai trouvé la fonction "subs" qui remplace le "t" par les valeurs voulues:

Code :

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

1
2
3
4
 
t_num=0:1e-2:10;
e_num=subs(e, t_num);
plot(t_num, e_num);

Jusqu'ici tout va bien mais j'ai un problème si la fonction "e" est constante (ne dépend pas du temps)

Code :

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

e=sym(2)

Dans ce cas la fonction "subs" me renvoie la valeur 2, alors que je voudrais un tableau de la même dimension que t_num où chaque élément serait égal à 2.
En fait j'aimerais pouvoir utiliser un même code que la fonction dépende du temps ou non, pour l'instant je doit les traiter différemment...
Y'a t'il une fonction que je n'aurais pas trouvé qui permette de faire cela ?

Merci

[phryte]

Bonjour.

je voudrais un tableau de la même dimension que t_num où chaque élément serait égal à 2.

Code :

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

e=2*ones(1,length(t_num))

[Julien44]

Bonjour,

Code :

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

e=2*ones(1,length(t_num))

Justement ça c'est ce que je fait actuellement si "e" est une constante, je cherche quelque chose qui marche dans tous les cas, que "e" soit une constante ou dépende de "tau"

Merci

[houcemc]

Bonjour,
Voilà un code qui peut t'aider.

Code :

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

1
2
3
4
5
6
7
8
 
clear
f=inline('2.*t','t');
for i=1:10,
    T(i)=i;
    F(i)=f(i);
end
plot(T,F,'r')

Vous pouvez aussi augmenter les valeurs de t:

Code :

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

1
2
3
4
5
6
7
8
9
 
clear
f=inline('2.*t','t');
e=0.1;
for i=1:100,
    T(i)=i*e;
    F(i)=f(i*e);
end
plot(T,F,'r')