Discussion :

[mimich77]

Bonjour,

Alors je suis actuellement bloqué sur mon code et j'ai cherché sans trouver une solution à mon problème.
Je m'explique, j'ai tracé un graph tout bête avec des coordonnées prédéfinies auparavant. Mais j'ai besoin d'obtenir, la valeur de y à x=0.
J'arrive à le faire manuellement avec la toolbox, en utilisant "basic fiitting" puis "shape-preserving interpolant" et en utilisant le "Data cursor" j'arrive aux coordonnées voulu (cad x=0).
J'aimerais avoir ça automatiquement dans mon code, malheureusement lorsque je fit avec "pchip" celui-ci commence à ma première valeur de x qui est supérieur à x=0. Du coup quand j'utilise la fonction find ou interp1, cela ne fonctionne pas car je sors de ma plage de données de x.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
x=[handles.Ti1;handles.Ti2;handles.Ti3;handles.Ti4;handles.Ti5;handles.Ti6;handles.Ti7];
y=[Mt_1;Mt_2;Mt_3;Mt_4;Mt_5;Mt_6;Mt_7];
title('\bf \fontsize{15} Signal = f (Inversion Time)')
xlabel('\fontsize{12} \color{red} Ti')
ylabel('\fontsize{12} \color{red} M(t)')
hold on
scatter(x,y,100,'x');
hold off
f=fit(x,y,'pchip');
plot(f,x,y);
grid on

Si quelqu'un peut m'éclairer sur mon problème, merci d'avance.

[gaspar.koenig]

Bonjour, la valeur de y en x=0 tend vers la valeur -Inf, Matlab peut avoir quelques difficultés avec +/-Inf. Tu peux essayer de créer un point supplémentaire en insérant une valeur 0 au début du vecteur x et -Inf au début du vecteur y. Mais je ne sais pas ce que donnera l'insertion de -Inf dans ton programme. D'autre part tu seras peut-être de toute façon obligé de remplacer -Inf par une valeur autre valeur très basse. Dis nous si cela a marché.

[mimich77]

Bonjour,
Tout d'abord merci de votre réponse.
Malheureusement, j'avais déjà essayé cela aussi. En mettant -Inf en y, il me dit que c'est impossible. Et de toute façon, il faut que Matlab me détermine la valeur de y en x=0. Cela n'aurait pas de sens d'insérer une valeur en début du vecteur y pour correspondre à x=0.



Ce que je disais : manuellement en ajoutant un "basic fitting" (courbe verte sur mon graph) et en utilsant le cursor, il me recréer la suite de la courbe jusqu'à x=0, je peux donc en interpoler la valeur en y.
Je pensais que comme Matlab était capable de le faire manuellement disons, cela aurait été aussi possible de le coder automatiquement.
Je cherche beaucoup et je ne trouve nulle part lorsque qu'il s'agit de trouver un point qui sort de la plage de données. Peut-être cela n'est pas possible à coder ?

Merci d'avance si vous avez d'autres solutions à me proposer.