Discussion :

[sebastien06]

Bonjour,

je me permets de vous écrire car j'ai une question sur la résolution d'équation. J'ai fais des recherches sur le help de matlab mais je ne trouve pas ce que je veux. Les forum et la faq non plus ne répondent pas à ma question.

Voilà ma problématique:
soit A une variable (vecteur de 45 lignes)
soit B une autre variable (même taille que A)

Je dois, en vue de déterminer une loi de corrélation entre A et B, résoudre l'équation:

B=a*B^n, où a et n seront à déterminer.

Cette équation répond à un modèle empirique que je teste sur mes données.

Je dois déterminer la meilleure solution pour a & n afin de coller au mieux à mes données. Enfin l'idée est de tracer le modèle créé par cette équation au dessus de mes données pour comparaison.

J'ai exploré les fonction relatives au polynomes (que j'utilise par ailleurs dans une seconde modélisation). A priori, de ce que j'ai lu, la seule fonction qui pourrait correspondre à mes besoins est 'solve'. Malheureusement elle n'est disponible que sur une toolbox payante.


Je reste à votre disposition si vous avez une solution.

En vous souhaitant une bonne soirée,

cordialement.

[phryte]

Bonjour.

Code :

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

B=a*B^n, où a et n seront à déterminer.

Où est A ?

[sebastien06]

Bonjour,

excusez moi il y a une erreur effectivement:

c'est B=a*A^n

merci

[phryte]

On peut utiliser les moindres carrés :

- non linéaire : Gauss-Newton
- Linéariser l'équation :

B=aA^n
log(B)=log(a) + n*log(A)
puis mettre sous forme matricielle :
|log(B(1)| |1 log(A(1)| |log(a)|
|log(B(2)| = |1 log(A(2)|*| n |
|--------| |-----------|
où :
Y=phi*X

et X = inv(phi'*phi)*phi'*Y
....

[sebastien06]

Bonjour,

pour info je viens de résoudre mon problème: mes recherches m'ont mené sur le site du mhs où un exemple de code est disponible pour différents types de fitting:

pour ceux que cela pourrait intéresser, voici le lien:

http://www.wcboe.k12.md.us/custom_pa...dex.html#curve

Cordialement

[sebastien06]

Bonjour Phryte,

merci pour ton idée. Je venais de le faire en fait avant de trouver l'exemple sur internet.

Cordialement