Discussion :

[auverland]

Bonjour à tous

Je me permet de vous solliciter pour un calcul qui à mon avis avis vous prendra que quelques minutes
Je débute completement en matlab je sais juste charger mes variables et faire un graph de base pour tout dire.

Voila j'ai une variable don't j'aimerais connaitre le nombre d'oscilation positive don't le deltat d'oscilation est supérieur à 2000.



exepmle dans mon image resultat attendu
6 amplitudes
et après si possible les classer par tranche

1 à 3300
2 à 3250
3 à 3200

Merci beaucoup si vous pouvez m'aider

[soft001]

Tu peux utiliser la fonction findpeaks + une condition sur l'amplitude.

et après si possible les classer par tranche

1 à 3300
2 à 3250
3 à 3200

J'ai pas compris !!

[Jean Dumoncel]

Bonjour,

cela semble être un travail pour la fonction findpeaks (nécessite la toolbox signal processing):

Code :

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

>> doc findpeaks

[auverland]

Tu peux utiliser la fonction findpeaks + une condition sur l'amplitude.


J'ai pas compris !!

En recherchant j'ai trouvé : findpeaks(select,Fs,'MinPeakProminence',2000)
Je met ou cette commande et je recupère mon résultat ou

Je voulais aussi si c'est possible le nombre de pics dans chaque plages (entre 2000 et 3000 et 3000 et 4000)
Merci d'avance

[soft001]

Tu peux essayer ce code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
load sunspot.dat
 
year = sunspot(:,1);
avSpots = sunspot(:,2)*20; % Ton signal, je pense que c'est toto dans ton cas
plot(year,avSpots);
[pks,locs]=findpeaks(avSpots,'MinPeakProminence',2000)
hold on
plot(year(locs),pks,'ro')
 
IndexPeakparplage(1)=numel(find(locs>0& locs<100)); % nombre de pics entre 0 et 100
IndexPeakparplage(2)=numel(find(locs>100& locs<1000)); % nombre de pics entre 0 et 100

[auverland]

Tu peux essayer ce code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
load sunspot.dat
 
year = sunspot(:,1);
avSpots = sunspot(:,2)*20; % Ton signal, je pense que c'est toto dans ton cas
plot(year,avSpots);
[pks,locs]=findpeaks(avSpots,'MinPeakProminence',2000)
hold on
plot(year(locs),pks,'ro')
 
IndexPeakparplage(1)=numel(find(locs>0& locs<100)); % nombre de pics entre 0 et 100
IndexPeakparplage(2)=numel(find(locs>100& locs<1000)); % nombre de pics entre 0 et 100

Bonjour,

Lorsque j'essais le code j'ai le message d'erreur (matlab R2013b) :

Error using uddpvparse (line 103)
The name 'MinPeakProminence' is not an accessible property for an
instance of class 'findpeaks'.

Error in findpeaks>parse_inputs (line 84)
hopts = uddpvparse('dspopts.findpeaks',varargin{:});

Error in findpeaks (line 59)
[X,Ph,Pd,Th,Np,Str,infIdx] = parse_inputs(X,varargin{:});

Y a t il une option à activer ?

Si je comprend le code, [pks,locs]=findpeaks(avSpots,'MinPeakProminence',2000)
celui-ci me donnera tout les pick supérieur à 2000 (en liste dans la fenetre de commande)
est-il possible de récupérer la liste sur une variable ?
Autre point, je n'ai qu'une variable 'toto' qui à comme valeur 1664617x1 double

Merci d'avance du temps que vous m'accordez

[Jean Dumoncel]

Bonjour,

les noms des options dans la R2013b sont légèrement différents, il suffit de lire la documentation pour corriger le nom.

Les variables pks et locs contiendront les valeurs des pics et leurs indices (encore une fois, tout cela est expliqué dans la documentation de la fonction).