Discussion :

[maritec]

bonjour
je veux determiner la fct de transfert par la methode de strejc sous matlab j'ai commencé par ce programme

Code :

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

1
2
3
4
5
t=[2 5 8 10 12 15 17 20 25 30 40];
y=[0.02 0.05 0.15 0.25 0.35 0.52 0.65 0.85 1.12 1.35 1.7];
plot(t,y,'b');
hold on
grid on

comment je trace la tangente pour determiner T1 et T2
merci

[phryte]

Bonjour.
Une idée (à la main. Avec plus de points il est peut-être possible d'automatiser).

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
26
27
28
t=[2 5 8 10 12 15 17 20 25 30 40];
y=[0.02 0.05 0.15 0.25 0.35 0.52 0.65 0.85 1.12 1.35 1.7];
p=diff(y)/3;
p2=diff(diff(y))/3;
tt=7;
ytg=p(tt)*t+y(tt)-(p(tt)*t(tt));
figure(1)
plot(t,y,'b');
hold on
plot(t,ytg,'r')
grid on
 
ym=polyfit(t(1:10),p,3);
poly=polyval(ym,t(1:10));
figure(2)
plot(t(1:10),poly,'r')
hold on
plot(t(1:10),p)
grid
title('Dérivée première')
ym=polyfit(t(1:9),p2,3);
poly=polyval(ym,t(1:9));
figure(3)
plot(t(1:9),poly,'r')
hold on
plot(t(1:9),p2)
grid
title('Dérivée seconde')

[maritec]

salut
pour ces instructions :

Code :

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

1
2
3
p=diff(y)/3;
p2=diff(diff(y))/3;
tt=7;

pourquoi vous avez divisé par 3 et c quoi tt=7

[phryte]

Bonjour.
1) Je croyais que 3 était la période. En fait elle n'est pas régulière. Il faut donc faire (pour les dérivées) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
t=[2 5 8 10 12 15 17 20 25 30 40];
y=[0.02 0.05 0.15 0.25 0.35 0.52 0.65 0.85 1.12 1.35 1.7];
p=diff(y)./diff(t);
p2=diff(p)./diff(t(1:10));
tt=7;
ytg=p(tt)*t+y(tt)-(p(tt)*t(tt));
figure(1)
plot(t,y,'b');
hold on
plot(t,ytg,'r')
grid on

2) tt=7 correspond au rand temporel du point d'inflexion (entre 15 et 20).

[maritec]

salut
pour le point d'inflexion es ce que il y a une instruction ou une formule qui permet de retirer ses coordonnés

[phryte]

Bonjour.

pour le point d'inflexion es ce que il y a une instruction ou une formule qui permet de retirer ses coordonnés

Oui, c'est le point où la dérivée seconde s'annule.
C'est pour cela que j'ai calculé son équation avec p2. Mais tu n'as pas assez de points pour faire ce calcul correctement.

[maritec]

Bonjour
es ce que on peux calculer le point d'intersection de la tangente avec la courbe pour retirer le point d'inflexion

[phryte]

Bonjour.
Le point d'inflexion est connu par la valeur zéro de la dérivée seconde.
Il te reste à calculer les intersection de la tangente avec y0 et yinf (mais là il manque des points)
Une version automatique :

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
clear
t=[2 5 8 10 12 15 17 20 25 30 40];
y=[0.02 0.05 0.15 0.25 0.35 0.52 0.65 0.85 1.12 1.35 1.7];
py=polyfit(t,y,4);
polyy=polyval(py,t);
dspy1=diff(polyy)./diff(t);%Dérivée première
dspy2=diff(dspy1)./diff(t(1:10));%Dérivée seconde
rang_zero=find(dspy2 < 0);%Calcul du rang du zero de la dérivée seconde
p=diff(y)./diff(t);
p2=diff(p)./diff(t(1:10));
tt=rang_zero(1);
ytg=p(tt)*t+y(tt)-(p(tt)*t(tt));
figure(1)
plot(t,y,'b');%Tracé de la fonction brute
hold on
plot(t,ytg,'r')%Tracé de la tangente au point d'inflexion
plot(t,polyy,'g')%Tracé de la fonction lissée (vérification du degré = 4 choisi)
grid on