Discussion :

[angelmn]

bonjour,
je ne sais pas si mon explication va être clair mais je vais essayer de vous expliquer mon problème.
en fait j'ai un signal semi-périodique avec un motif qui se répète disons dans ce genre:

je cherche au final à rendre mon signal périodique, donc j'ai détecté les pics et j'ai essayé de ré-échantillonner la distance entre chaque deux pics de façon à avoir la même distance partout.
j'utilise la fonction resample de matlab
mais voici un exemple d'un échantillonnage d'une période:

Le problème c'est que ça me déforme mon signal , surtout le pic que j'ai le reste m'importe peu mais je veux garder le même aspect pour mon motif répétitif.
Auriez vous une idée de ce que je devrais faire??
Merci

[Jerome Briot]

Une idée plus 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
20
21
22
23
24
25
26
27
28
29
% Donnees initiales
x = ...;
y = ...;
 
% Indices des pics
idx = ...; 
 
xper = [];
yper = [];
 
for n = 1:numel(idx)-1
 
    subx = x(idx(n):idx(n+1));
    suby = y(idx(n):idx(n+1));
 
    xi = linspace(subx(1), subx(end), 30);
    yi = interp1(subx,suby,xi);
 
    nper = numel(xper);
    xper = [xper (1:numel(yi))+nper];
    yper = [yper yi]; 
 
end
 
figure
subplot(121)
plot(x,y,'b-',x(idx),y(idx),'r*')
subplot(122)
plot(xper,yper,'b-');

[le fab]

angelmn, peux tu fournir ton code et un jeux de données ?
je n'utilise pas (encore) les time series , je vais m'y mettre ... du coup ça m’intéresse de regarde de plus près

Dut : juste au cas ou angelmn soit intéressée par une interpolation de type Zoh (zero order hold), ta solution ne peut pas marcher
(c'est pour moi le principale intérêt de la fonction resample, sinon autant utiliser interp1)

[angelmn]

merci de m'avoir répondu. DUT je n'ai pas vraiment compris votre code je vais essayer de faire ca tout de suite. Sinon voici mon code:

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
%%%%%%%%%%%%detection des pics%%%%%%%%%%%
j=1;
for i=1:length(x)
    if x(i)>15
        T(j)=i;
        j=j+1;
    end
end
%%%%%%%%%%%réechantillonage%%%%%%%%%%%%%%%
p=0;
for i=1:length(T)-1
    s=x(T(i):T(i+1))% récupère le bout de signal entre 2 pics successifs
    sg(1+p:25+p)=resample(s,25,length(s))%reechatillone pour avoir une periode de 25
    p=p+25;
end
%%%%%%%%%%affichage%%%%%%%%%%%%%%%
subplot(2,1,1)
plot(x);
axis([0 80 0 25])
hold on 
plot(T,x(T),'r*');
subplot(2,1,2)
plot(sg);

et ça me donne ça:

ce qui me dérange le plus dans le résultat c'est que les maxima locaux ont changé et donc mon signal n'est pas devenu periodique. Moi je croyais pouvoir allonger la distance entre deux pics tout en les gardant intacts

[Jerome Briot]

Donne nous les valeurs de x dans un fichier .mat (lui-même dans une archive zip ou rar)

[angelmn]

le voici mon signal x ( ce n'est pas le signal final sur lequel je souhaite travailler c plus compliqué mais j'essai dejà de resoudre le problème avec ce signal)
Merci de votre aide