Discussion :

[membreComplexe12]

Bonjour tous,

j'ai un nuage de points et je souhaite l'interpoler avec un modèle de type "circuit RL" (en identifiant A et B):

Code :

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

y=A(1-exp(-x/B))

actuellement j'ai fais une interpolation sur matlab de type moindre carré avec une exponentielle :

Code :

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

1
2
3
4
5
Y = load('test.txt','-ascii');
Xi = linspace(Y(1,1),Y(end,1),100);
p = polyfit(Y(:,1),log(Y(:,2)),1);
m = p(1); n = exp(p(2));
Yi = n*exp(m*Xi);

mais je ne sais pas comment tranformer ce bout de code pour avoir une interpolation avec le modèle que je vous ai montré plus haut...

Pourriez vous m'aider svp ?

ps: j'aimerai bien aussi identifier ceci : A.exp(-B.X) mais en imposant une valeur à A ....

[Aleph69]

Bonjour,

je n'ai pas essayé mais en annulant le gradient de l'erreur quadratique entre y et A(1-exp(-x/B)), tu devrais normalement pouvoir identifier A et B.

[membreComplexe12]

oui mais en fait je voulais utiliser la methode polyfit de matlab mais ça n'a pas l'air possible tel quel

merci quand même