Discussion :

[Kichott]

Bonjour,

J'ai une fonction de transfert pour un système donné.
A l'aide de matlab je passe en temporelle grâce a ilplace.

Maintenant je bloque sur quelque chose qui parait pourtant simple.

>> y = ilaplace(w)
y = v/exp(t) + t*v + (2*3^(1/2)*v*sin((3^(1/2)*t)/2))/(3*exp(t/2))

Maintenant je veux calculer y pour différentes valeur de v et de t.

Comment puis-je faire ?

Mon code :

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
25
% Initialization of variable
ms = 1 ;
mus = 1 ;
Ks = 1 ;
Bs = 1 ;
Kt = 1 ;
Bt = 1 ;
 
% Calculation of coefficients
a4 = ms*mus;
a3 = ms*Bs+ms*Bt+mus*Bs;
a2 = ms*Ks+ms*Kt+Bs*Bt+mus*Ks+Bs*Bs-Bs*Bs*mus/ms;
a1 = Bs*Ks+Bs*Kt+Bt*Ks-Bs*Ks*mus/ms;
a0 = Ks*Kt;
b5 = ms*Bt;
b4 = ms*Kt+Bs*Bt;
b3 = Bt*Ks+Bs*Kt;
b2 = Ks*Kt;
 
%Function transfert
syms s v;
 
w = (a4*s.^4 + a3*s.^3 + a2*s.^2 + a1*s + a0 ) * v / ( b5*s.^5 + b4*s.^4 + b3*s.^3 + b2*s.^2 ) ;
 
y = ilaplace(w)

[chris05]

Bonjour,

personnellement je ferai une double boucle :

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
 
%initialisation
v=[];
t=[];
y_sol=[];
 
 
v=[1:100]; %vecteur contenant les valeur de v (a remplacer par tes valeurs)
t=[1:100]; %vecteur contenant les valeur de t (a remplacer par tes valeurs)
 
for i=1:length(v)
     for j=1:length(t)
          y_sol(i,j)= v(i)/exp(t(j)) + t(j)*v(i) + (2*3^(1/2)*v(i)*sin((3^(1/2)*t(j))/2))/(3*exp(t(j)/2));
 
     end
end

y_sol est une matrice.
sur une meme ligne tu trouveras les valeurs de y pour un v fixé et t variable.
sur une meme colonne tu trouveras les valeurs de y pour un t fixé, v variable.

Pour tracer y pour v fixé (3eme valeur de v) et t varible :

Code :

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

1
2
 
plot(t,ysol(3, : ))

Pour tracer y pour t fixé (50eme valeur de t) et v varible :

Code :

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

1
2
 
plot(v,ysol( : ,50))

[Kichott]

Je suis d'accord que l'on peut faire ça...

C'est ce que je pensais faire si je ne trouvrais pas d'autres solutions.

Mais je voudrais éviter de réécrire la fonction, plutôt la récupérer de la fonction d'avant (ilplace).

Donc je cherche à modifier les valeurs de t et v, mais du fait que je sois en symbolique je ne peux pas changer les valeurs. En faisant un simple t=... ou v=...

Donc pour être plus précis comment puis-je changer les valeurs si je les ai déclarées avant en symbolique ?

Merci

[Caro-Line]

Je ne suis pas sure n'ayant pas la Symbolic toolbox mais je dirais tout simplement :

Code :

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

1
2
3
t=1;
v=1;
res=double(y);

Pour l'utilisation de double dans ce cadre : DOUBLE.

[Kichott]

Je ne suis pas sure n'ayant pas la Symbolic toolbox mais je dirais tout simplement :

Code :

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

1
2
3
t=1;
v=1;
res=double(y);

Pour l'utilisation de double dans ce cadre : DOUBLE.

Non malheureusement ça ne fonctionne pas.

Je comprends pas je fais :

Code :

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

1
2
3
4
5
 
>> clear t v syms
>> t=1
>> v=1
>> y

Et il me renvoie tjs la fonction sans modifier le t ou le v

[Jerome Briot]

Code :

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

1
2
3
4
5
 
>> clear t v syms
>> t=1
>> v=1
>> y

Je dirais plutôt :

Code :

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

1
2
s = subs(y,{t,v},{1,1})
double(s)

non ?

[Kichott]

Ca ne marche pas non plus.

Il me dit :

Code :

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

  ??? Undefined function or variable 't'.

Le t vient de l'inverse de la transformee de laplace qui transforme le s en t...
La variable v est par contre est bien modifie.

[Jerome Briot]

Ne connaissant pas bien cette Toolbox, peux-tu nous dire ce que renvoie ceci ?

Code :

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

1
2
y = ilaplace(w);
whos

Si y est de classe Char, il est peut être possible de faire :

Code :

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

1
2
3
4
5
 
y = ilaplace(w);
yy = inline(y,'v','t');
 
yy(1,1)

Je ne peux même pas dire si ça fonctionne mais ça ne coûte rien d'essayer