[csaps] Pour chaque x, extraire les y d'une fonction

**mikael.jaffre** · 18/10/2011, 14h31

Bonjour à tous,

en partant d'une vecteur y(120,1) j'effectue un lissage par 'Cubic spine smoothing' que je représente graphiquement ci-dessous.
En fait je cherche à extraire pour chaque x (1:120) une valeur de y issue de mon lissage... et je n'y arrive pas.

Merci d'avance pour votre aide !

Mikaël

Code :

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

1
2
3
4
5
6
7
8
9
10
11
clear all
y=[0.0670241286863271, 0, 0, 0, 0, 0, 0, 0.109051254089422, 0, 0.0794912559618442, 0.0731528895391368, 0, 0, 0, 0, 0, 0.0723065798987708, 0.0915750915750916, 0,0.0769230769230769, 0.287356321839080, 0, 0, 0.452830188679245, 0.874403815580286, 0.0769230769230769, 0.311769290724864, 0.575184880854560, 2.41404535479151, 3.62264150943396, 5.47636909227307, 7.84313725490196, 1.31578947368421, 15.9245283018868, 33.5094339622642, 65.0769230769231, 77.1569433032046, 64.7051298452092, 65.4199011997177, 65.7358490566038, 83.5360908353609, 113.688888888889, 131.644583008574, 46.7125382262997, 17.5094339622642, 59.6153846153846, 149.551020408163, 277.811320754717, 63.3671083398285, 18.6544342507645, 6.52173913043478, 1.85967878275571, 7.48792270531401, 0.541516245487365, 18.1531176006314, 11.2307692307692, 2.55957634598411, 8.62745098039216, 5.31511009870919, 4.53846153846154, 0, 0, 4.32653061224490, 6.96000000000000, 1.89423835832676, 0.868192580899763, 0.784313725490196, 0, 0, 0.251677852348993, 0.434782608695652, 0.338123415046492, 0.340425531914894, 0.0821692686935086, 0, 0.166666666666667, 0.0856898029134533, 0.0869565217391304, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.255319148936170, 0, 0,  0, 0, 0, 0, 0,  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0];
x=1:120;
 
fnplt(csaps(x,y,0.2),2)
hold on
fnplt(csaps(x,y,0.5),2,'k--')
fnplt(csaps(x,y,0.8),2,'r:')
plot(x,y,'+g')
legend('p = 0.2','p = 0.5','p = 0.8','données brutes'),
hold off