Discussion :

[Seikyo67]

Bonjour,
J'ai un enregistrement sous forme de matrice que je souhaite analyser
Pour cela, je les stocke dans une matrice A
Je retire la moyenne temporelle ( j'ai utilisé la fonction mean mais je ne suis pas sur que ça soit la bonne) et stocke ces nouvelles valeurs dans une matrice AM
J'annule les parties négatives dans une matrice AMM

Il me faut maintenant détecter les périodes
Pour cela, j'ai utilisé le code suivant

Code :

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

1
2
3
4
5
6
7
8
9
for i=2:size(A)-1  
    if AMM(i)>0
        if AMM(i-1)==0
            L=[L i*10e-3];
        elseif AMM(i+1)==0
            L=[L i*10e-3];
        end
    end
end

Je dois également trouver les valeurs maximales de l'enregistrement pour chacune de ces périodes, mais à ce stade là, je sèche un peu. J'aurai voulu savoir si quelqu'un avait une idée sachant que j'ai un grand nombre de valeurs à analyser
Merci

[awawawa]

Je ne suis pas sur d'avoir bien compris ton problème. Tu as un signal composé de plusieurs périodes et tu veux connaitre la "valeur" de chacune des périodes? Je te conseille de passer par un Transformée de Fourrier...

[Paul-Ceria]

Si tu as qu'un seul pic dans chaque période :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
valeur_max=[];
for i=2:size(A)-1
    if AMM(i)>0
        if AMM(i-1)==0
            L=[L i*10e-3];
        elseif AMM(i+1)==0
            L=[L i*10e-3];
        end
 
        if AMM(i-1)<AMM(i) && AMM(i+1)<AMM(i)
            valeur_max(size(valeur_max,2)+1)=AMM(i);
        end
    end
end

Si chaque période est constituée de plusieurs pics :

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
temp=[];
valeur_max=[];
for i=2:size(A)-1
    if AMM(i)>0
 
        if AMM(i-1)<AMM(i) && AMM(i+1)<AMM(i)
            temp(size(temp,2)+1)=AMM(i);
        end
 
        if AMM(i-1)==0
            L=[L i*10e-3];
        elseif AMM(i+1)==0
            L=[L i*10e-3];
 
            valeur_max(size(valeur_max,2)+1)=max(temp);
            temp=[];
        end
 
    end
end