Discussion :

[SO-2016]

Bonjour,esque quand on donne un poids zero au points dans un intervalle de temps bien précis équivaut a ne pas prendre en considération ces points ( ils n'existent pas dans cet intervalle ) lors d'un fit linéaire donc la droite tracée dans cet intervalle n'est pas un fit linéaire ??Merci

[le fab]

Salut

non, mais il suffit de les enlever de l'échantillon pour le fitting

Fabien

[SO-2016]

Merci Fabien, mais j'aurai quoi dans cet intervalle de temps si je les enlèvent ? normalement y'aura des points avec un poids zero et non des zéro donc normalement il y aura rien c'est ça?

[le fab]

heu oui je crois (si j'ai bien compris ce que tu écris ;-))

[SO-2016]

Merci , mais je veux des explications convaincantes

[phryte]

Bonjour,

Comme je te l'ai indiqué dans l'autre discussion, des données de poids nuls et des données de valeurs nulles n'est pas la même chose.

Dans le cas d'exclusions (poids nuls) les données sont ignorées (courbe rouge) dans l'autre cas les données sont prises (courbe bleue)

Un exemple simple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear
xdata=[1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 1900 1910 1920,...
1930 1940 1950 1960 1970 1980 1990]';
ydata=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21]';
% Zones de données de poids nuls
raz1=[1:5];
raz2=[11:15];
Indice_nul=[raz1 raz2];
outliers = excludedata(xdata, ydata, 'indices',Indice_nul);
f=fit(xdata,ydata,'poly1','Exclude',outliers,'Robust','on')
% Données avec des valeurs nulles aux mêmes rangs
ydata=[0 0 0 0 0 6 7 8 9 10 0 0 0 0 0 16 17 18 19 20 21]';
f1=fit(xdata,ydata,'poly1')
h1=plot(f,xdata,ydata)
set(h1,'Color','r','LineWidth',2)
hold on
h2=plot(f1,xdata,ydata)
set(h2,'Color','b','LineWidth',2)
grid

[SO-2016]

Je vous remercie infiniment phryte pour votre réponse et explications ,mais quand on ignorent les données (poids nul) on fit quoi alors ?poids nul=données ignorées= non existence des données dans ce cas on fait un fit entre ces deux points seulements [0 et6] au départ et [10 et 16] a la fin (je suis confuse)
Une dernière question es que on peut dire tout simplement que la droite rouge est un fit linéaire pour des données avec un poids nuls qui sont ignorées ?
Merci beaucoup phyrte pour votre aide

[phryte]

Bonjour,

on peut dire tout simplement que la droite rouge est un fit linéaire pour des données avec un poids nuls qui sont ignorées ?

OUI si tu utilises l'option 'Exclude'

[SO-2016]

Merci phryte pour votre aide mais si je veux éliminer cette fois ci un seul intervalle [1:5] et fiter aprés , comment je dois procéder ?normalement

Code :

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

1
2
3
4
5
6
7
xdata=[1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 1900 1910 1920,...
1930 1940 1950 1960 1970 1980 1990]';
ydata=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21]';
 
Indice_nul=[1 :5];
outliers = excludedata(xdata, ydata, 'indices', Indice_nul);
f=fit(xdata,ydata,'poly8','Exclude',outliers,'Robust','on');

mais ça ne marche pas ,je ne comprends pas pourquoi?? help me please

[phryte]

Bonjour,

poly8 n'existe pas.

'poly1' ou 'poly2'

Fais : help fit

Si tu veux l'ordre 8 il faut faire :

Code :

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

1
2
3
P = polyfit(xdata,ydata,8);
 V=polyval(P,xdata);
 plot(xdata,V)

[SO-2016]

Bonsoir phryte mais je voulais donné un poids nul a un seul intervalle cette fois ci (pas deux) et fiter avec un polynôme de degré 8
Comment utiliser outliers avec un fit de degré 8

Code :

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

1
2
3
4
5
Indice_nul=[1 :5];
outliers = excludedata(xdata, ydata, 'indices', Indice_nul);
P = polyfit(xdata,ydata,8,outliers);
 V=polyval(P,xdata);
 plot(xdata,V)

mais ça ne marche pas?
Sincèrement je ne saurai comment vous remercier phryte

[phryte]

Bonjour,

polyfit n'accepte pas excludedata (hel polyfit)

Avec fit tu as ces cas :

FITTYPE DESCRIPTION
'poly1' Linear polynomial curve
'poly11' Linear polynomial surface
'poly2' Quadratic polynomial curve
'linearinterp' Piecewise linear interpolation
'cubicinterp' Piecewise cubic interpolation
'smoothingspline' Smoothing spline (curve)
'lowess' Local linear regression (surface)




Exemple avec l'ordre 3 (l'ordre 8 n'est pas forcément adapté au fit (==> message d'erreur) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
clear
xdata=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21]';
ydata=10+2*randn(21,1);
plot(xdata,ydata,'.r')
hold on
%Indice_nul=[1:5];
%fit avec les échantillons de 1 à 5 en moins
[P,S] = polyfit(xdata(6:end),ydata(6:end),3);
V=polyval(P,xdata);
plot(xdata,V,'k')
grid
axis([0 21 0 20])
%fit avec tous les échantillons
[P1,S1] = polyfit(xdata ,ydata ,3);
V1=polyval(P1,xdata);
plot(xdata,V1,'g')

[SO-2016]

Merci phryte mais ce que je veux vraiment c'est de donner un poids nul a un intervalle bien précis (entre 10:12 mn) puis fiter avec un fit polynomial de degrés 7 ou 8

Code :

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

1
2
3
4
5
6
7
% Data from data set "Moyenne vs. t with w":
%     x= t:
%     y = Moyenne:
%     Weights = w:
w=ones(8401,1);
t=linspace(0,35,length(y));%le temps en minutes 
y = Moyenne:% moyenne est une matrice (8401x1)

donc
1-je veux donner un poids zero aux données dans l'intervalle de temps (10:12) puis fiter le résultats avec un fit polynomial de degrés 7ou 8
2- je donne cette fois ci un poids zero au 2 intervalle (10:12) et (24:29) puis fiter le résultat toujours avec un fit polynomial de degrés 7ou 8
(sachant qu'il faut donner un poids =1 pour toutes les données sauf dans les intervalles (10:12 ou 24:29)ou il faut donner un poids zero)

[le fab]

avec un vecteur temps t en minute (de la même taille de xdata et ydata), il suffit d'écrire qqch comme ça :

Code :

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

[P,S] = polyfit(xdata(t>10 & t<=12),ydata(t>10 & t<=12),3);

du coup les données comprises entre 10 et 12 mn ne sont pas transférées à la fonction polyfit
ce qui reviens à un poids nul

[SO-2016]

Merci le fab , mais ça ne marche pas avec un polynôme de degré 7 ou 8

[phryte]

Bonjour,
Même solution que le fab .
Dis-nous plutôt ce que tu cherches à faire. Il y a d'autres solutions que fi ou polyfit :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
clear
Nx=100;
w=ones(Nx,1);
t=linspace(0,35,Nx);%le temps en minutes 
%y = Moyenne:% moyenne est une matrice (8401x1)
 y= 10*rand(1,Nx)+5;
 
plot(t,y)
axis([0 35 0 20])
grid
hold on
[P,S] = polyfit(t,y,8);
yf=polyval(P,t);
plot(t,yf,'r')
% Suppression des valeurs de 10 à 12
 ym=[y(1:9) y(13:end)];
 tm=[t(1:9) t(13:end)];
 [Pm,Sm] = polyfit(tm,ym,8);
yfm=polyval(Pm,tm);
plot(tm,yfm,'.k','MarkerSize',13)
%
% Solution avec fit et smoothingspline
Indice_nul=[10:12];
outliers = excludedata(t,y,'indices',Indice_nul);
f=fit(t',y','smoothingspline','Exclude',outliers);
hm=plot(f,t,y)
set(hm,'Color','g','LineWidth',2)

[SO-2016]

Bonjour phryte, j'ai un vecteur de 8401 lignes en fonction du temps (35mn) et je veux :

1- donner un poids zero aux données de l' intervalle de temps (10:12) puis fiter le résultats avec un fit polynomial de degrés 7ou 8 (c'est important que le degré soit grand 5,6,7,8,9) et en même temps il faut que l'écart type soit minimum

2- donner un poids zero aux 2 intervalles (10:12) et (24:29) puis fiter le résultat toujours avec un fit polynomial de degrés 7ou 8

3- tracer le différence entre le signal d'origine et son fit pour l'utiliser après

sachant qu'au départ il faut donner un poids =1 pour toutes les données( sauf pour les intervalles (10:12 ou 24:29)ou il faut donner un poids zero).

Remarque importante :on veut pas supprimer les données entre (10:12) mais leurs donner un poids zéro pour ne pas les prendre en considération lorsqu'on fit

Je tiens vraiment à vous remercier pour votre aide.Ci joint mes deux vecteurs: Moyenne.rar et t.rar

[phryte]

Bonjour,
On a compris ce que tu veux faire, mais pas le pourquoi.
Le lissage d'ordre élevé est délicat à utiliser et entraîne des erreurs.
Ce que l'on t'a donné te suffit ?

[SO-2016]

Mais phryte si je prends votre dernier code la taille du fit "yfm "est différente de la taille de" y "du coup je ne peux avoir leur différence(y-yfm) puisque vous avez supprimer toutes les données de l'intervalle (10:12) pour yfm( il faut juste leur donner un poids nul )

[phryte]

Si tu veux rephaser les signaux il faut faire (en ignorant la zone 10,11,12 :

Code :

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

1
2
yfmm=[yfm(1:9) zeros(1,3) yfm(10:end)];
tmm=[tm(1:9) zeros(1,3) tm(10:end)];

et utiliser yfmm.

Tout cela n'est peut-être pas très rigoureux !