Discussion :

[amand237]

Bonjour,

J'utilise la fonction [fittype] pour fitter mes données en sigmoïde et trouver la valeur de la pente au point d'inflexion de ma sigmoïde. Mais j'aurai besoin de récupérer les différents points fittés pour d'autres calculs (trouver la valeur de x à y=0.5 et calculer les aires sous courbes de chaque côté de cette valeur de x).

Informations complémentaires: grâce au programme que j'ai de fit, je peux précupérer les valeurs de a et b de l'équation type utilisée pour le fit (je pense que ça peut éventuellement aider pour récupérer l'équation de la courbe et l'utiliser pour faire mes calculs.

Désolé si ce n'est pas clair mais je suis débutante et le premier programme a été développé pour moi par le service informatique de mon bureau.

Si ça peut aider voilà une partie du code de mon programme. Je ne sais pas si ça donne suffisament d'informations.

Code :

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

1
2
3
4
5
 
%calcul du fit :
%courbe a fitter : sigmoide
f = fittype('b/(1+exp(-a*x))');
fit1 = fit(X,Y2,f,'Startpoint',[1 1])

Egalement voilà ce que j'obtiens en retour sur ma fenêtre de commande:

fit1 =

General model:
fit1(x) = b/(1+exp(-a*x))
Coefficients (with 95% confidence bounds):
a = 2.956 (-1.748, 7.659)
b = 0.9252 (0.761, 1.089)

valeur de la pente en 0 : 0.68372

J'espère que vous pourrez me répondre.

Bonne journée

[phryte]

Bonjour.
Exemple de calcul de la pente à l'inflexion (dérivée seconde = 0) :

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
clear
 X=0:0.1:2;
 Y2=0:0.1:2;
f = fittype('b/(1+exp(-a*x))');
fit1 = fit(X',Y2',f,'Startpoint',[1 1])
 a=fit1.a
 b=fit1.b
 %Tracés
x=-10:0.1:10;
fit1 = b./(1+exp(-a*x)); 
plot(fit1)
grid
title('Sigmoide')
figure(2)
dfit1=diff(fit1);
plot(dfit1)
grid
title('Dérivée première')
dfit2=diff(dfit1);
figure(3)
plot(dfit2)
grid
title('Dérivée seconde')
x0=find(dfit2==0)
x0=find(dfit2==0)
ordonnee_inflexion=fit1(x0)
pente_inflexion=dfit1(x0)