Maximum local dans Matlab en présence d'un BRUIT

Version imprimable

21/05/2007, 16h55
shinobida

bonjour,

j'ai un graphique où il y a 3 dos d'ânes

Le 1er et le 3ème sont dus au signal (donc utile) , le 2ème (celui au milieu) est du à l'interférence co-existante entre eux. Donc ce que je veux c'est obtenir la localisation de ces 2 dos d'ânes (obtenir la localisation des 3 et supprimer celui au milieu)

Comment faire ?
Merci

Salut,

Essaie cet exemple:

Code:

1
2
3
4
5
6
7
8
 
z=peaks;
z=z(30,:);
d=diff(sign(diff(z)));
j=find(d==-2)+1;
plot(z)
hold on
plot(j,z(j),'r*')

22/05/2007, 11h10
shinobida

1 pièce(s) jointe(s)

Merci boss, Mais j'ai vu le coté théorique et que j'ai essayer d'implementer mon graphe est le suivant, et la ya deux GRANDs pics aux extrimité (je veux pouvoir avoir cela), et le pb c'est que j'ai 64 graphe comme ça, et il faut que je reli les max (les deux) entre eux pour que j'obtienne deux droite sur le graphe, et donc à priori j'ai aucune infi sur le singal.

Merci de m'aider

MMMMMMMMMMMerci rostumus tu m'a bcp aidé, donc voici le code final ( pour ceux qui sont intéréssé)
merci encore rostumus:king:

Code:

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
for i=10:N-1
    z=Wr(:,i);
    d=diff(sign(diff(z)));
    j=find(d==-2)+1;
    %[val,posit]=max(z(j));
%    plot(z)
%hold on
%plot(j,z(j),'r*')
%    for o=1:length(z(j))
 %       M(o)=z(j(o));
  %  end
  %m=z(j)
    for u=1:3
 
        [val,posit]=max(z(j))
        posipeak(u)=j(posit);
        z(j(posit))=0;
    end
 
    finalposipeak(1,i)=min(posipeak);
    finalposipeak(2,i)=max(posipeak);
 
end
 
%axes(handles.PlotIF);
 
 
%finalposipeak=finalposipeak+ones(2,63);
for i=1:N-1;
    %i
       %%plot(absc,abs(Wr(i,:))+20*(i/164));
       %%celle la
 axes(handles.PlotIF);
 
       plot(freqc,(abs(Wr(:,i))+6*(i/taille_temps)));
hold on
if  finalposipeak(1,i)~=0;
plot(freqc(finalposipeak(1,i)),(abs(Wr(finalposipeak(1,i),i))+6*(i/taille_temps)),'r*');
end
 
if  finalposipeak(2,i)~=0;
%plot(freqc(finalposipeak(1,i)),(abs(Wr(finalposipeak(1,i),i))+6*(i/taille_temps)),'r*');
hold on
plot(freqc(finalposipeak(2,i)),(abs(Wr(finalposipeak(1,i),i))+6*(i/taille_temps)),'g*');
end
hold on
 
end