Discussion :

[trantsyx]

Bonjour,

je réalise en ce moment une modélisation qui requiert l'utilisation de splines. Et je bloque sur un problème qui sera peut-être simpliste pour vous.

J'ai un ensemble de 1200 valeurs nommé DifMed le long de l'axe des x=1:1200, oscillant entre -10 et +20.
Mon objectif est de rechercher les valeurs de x pour lesquelles DifMed égale 0.

Afin de rendre DifMed continue, je dois utiliser des splines. J'arrive à faire passer ma spline par les points, mais je n'arrive pas à trouver la fonction ou la commande qui me permettrait de demander à Matlab de me donner les valeurs de x lorsque y=0.

J'ai bien essayé avec ppval, qui me donne des ordonnées pour une abscisse, soit l'inverse de ce que je veux, et toute mes tentatives avec fsolve se sont soldées par des échecs...
Si par hasard vous auriez une idée... Merci d'avance!

[FR119492]

Salut!

J'ai un ensemble de 1200 valeurs nommé DifMed le long de l'axe des x=1:1200

Tu as une fonction échantillonnée, ce qui signifie que tu n'as aucune information sur ce qui se passe entre deux points consécutifs. Le spline te donne l'illusion de le savoir. C'est pourquoi je te suggère une méthode beaucoup plus simple: tu parcours ta liste de 1200 valeurs et tu localises

• les valeurs nulles;
• les changements de signe.

Jean-Marc Blanc

[trantsyx]

Bonjour,

merci de votre réponse.
Je n'ai pas été très loquace sur mon projet, je m'en excuse.C'est un processus graduel, on peut donc croire sans trop se tromper que deux points consécutifs sont reliés par quelque chose de plus ou moins droit.

J'ai essayé le code x1=find(spline1==0) dans mes 1200 valeurs, mais aucune n'est strictement égale à 0, d'où ma volonté de chercher le changement de signe et les valeurs nulles le long de la spline. Mais c'est là mon problème, j'ignore comment le demander à matlbab.
Merci!

[Jerome Briot]

Dans un premier temps, on peut utiliser les fonctions SIGN et STRFIND pour déterminer les indices des valeurs consécutives qui changent de signe :

Code :

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

1
2
3
4
s = sign(X);
 
idx1 = strfind(s,[1 -1]);
idx2 = strfind(s,[-1 1]);

[trantsyx]

Merci beaucoup!
Je devrais pouvoir en faire quelque chose

Mais vraiment, n'existe-t-il aucun moyen simple de trouver les 0 d,une spline, un peu comme la commande fzero pour trouver les x pour que f(x)=0 ?
Je désespère

Mais merci pour votre aide, à bientôt!

[floopy]

Bonsoir,

As tu regardé ce post :
http://www.mathworks.com/matlabcentr..._thread/280148 ?
Je n'ai pas testé la fonction rootspline, mais apparement, elle est censée faire ce que tu veux, non ?

EDIT: je viens de tester l'exemple fournit avec le code, elle a l'air de marcher