Discussion :

[thomMonteillet]

Bonjour,

je me permets de vous déranger car j'ai deux "petites" questions.

1) Tout d'abord,j'aimerai savoir si l'un d'entre vous sait comment modifier le calcul du résidu de la fonction "lsqnonlin" ?

Actuellement, si j'ai bien compris, on fournit par l'intermédiaire d'une fonction un vecteur que je notera R à "lsqnonlin". Ce vecteur est la différence entre un "vecteur objectif" et le vecteur calculé à l'instant "t".

Ensuite, "lsqnonlin" faire cette opération afin d'obtenir la valeur du résidu :

Code :

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

residu=R'*R;

Mon soucis est que j'aimerai garder le même algorithme mais modifier juste la façon avec laquelle il calcul le résidu, par exemple en lui disant de faire comme ceci :

Code :

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

residu=R'*T'*R*T;

ou

Code :

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

residu=abs(R);

savez vous si cela est possible ?
(j'ai pensé à aller dans le code source afin de mettre mon opération mais ça me parait une mauvaise idée...)

2) ma seconde question est sur le nombre de point d'optimisation...

Ma deuxième question est un peu plus difficile à formuler :

J'ai une courbe en bleu sur ma figure qui est la courbe objectif et en noir la courbe simulée.

On voit qu'après optimisation "lsqnonlin" à privilégié lors de l'identification
la partie de la courbe qui a la plus grande amplitude (car il y a plus de points que dans la zone de plus faible amplitude). Du coup, je voudrais savoir si il est possible de faire quelque chose pour que l'identification soit indépendante du nombre de points d'une ou l'autre des parties de courbes... ?

En fait j'aurais préféré que "lsqnonlin" me fasse un signal triangulaire d'amplitude moyenne constante plutôt que quelque chose qui est bon à la fin mais mauvais au début.

En espérant être clair dans mes questions.
merci.

[le fab]

pour le première question :
je ne vois pas comment modifier le calcul du résidu, par contre dans les options de lsqnonlin, tu as OutputFcn qui te permet d'appeler une fonction faite par tes soins à chaque itération

pour ta 2ème question :
- la forme de ta solution dépends de la fonction à fitter que tu as écris
si tu veux une amplitude constante il suffit décrire une fonction qui a une amplitude constante
- concernant le nombre de point qui est plus important dans une partie de tes données, tu peux arranger les choses soit en faisant ton optimisation sur une partie de tes données (tu coupes à droite) soit en augmentant le nombre de points à droite par interpolation (interp1)
- dans ton cas l’idéal est même d’écrire une fonction ayant deux amplitudes constantes (if x<x0,...amplitude1...;else,...amplitude2...;end) en augmentant le nombre de paramètres à faire trouver à lsqnonlin

[thomMonteillet]

Merci beaucoup pour tes réponses très intéressantes !!!
je pense que ça va débloquer la situation.

A+

ps: ma fonction est fixe, je ne peux pas la modifier