MATLAB：データ反復からピークを見つける

Question

私は信号からfft関数の大きさをプロットする関数を持っています。

繰り返しごとに、2000未満の2つのピークのx値を求めたいと思います。これは、関数findpeaksを使用して比較的単純だと思っていましたが、正しい出力を得られませんでした。

私は出力に含まをプロットするつもりはないが、ちょうど説明のためにここにプロットです。ここに私のコードのビットが

：私だけ

（ピークの最初のセット）2000、以下のデータのピークを知るために一回の反復の例をしたいです。 Bは、解析が必要なデータのすべてのセグメントの開始インデックスを含むベクトルです。

function [number] = fourir_(data,sampling_rate) 

%Finds the approximate starting index of every peak segment 
%B is a vector containing the indeces 
[A,B] = findpeaks(double(abs(data) > 0.6), 'MinPeakDistance', 2500); 

Fs = sampling_rate; 
t = 0:1/Fs:0.25; 

C = zeros(size(B),2) 
    for i = 1:numel(B) 
    new_data = data(B(i):(B(i)+200)) 
    y = double(new_data)/max(abs(new_data)); 
    n = length(y); 
    p = abs(fft(y)); 
    f = (0:n-1)*(Fs/n); 

    end 

例データ：https://www.dropbox.com/s/zxypn3axoqwo2g0/signal%20%281%29.mat?dl=0は

Answer 1

は、あなたの答えは、これは@Edスミスは彼の最初のコメントで提案まさにです。メジャーピークを区別するためにスレッショルドを追加することができます。

%Finds the approximate starting index of every peak segment 
%B is a vector containing the indeces 
[A,B] = findpeaks(double(abs(data) > 0.6), 'MinPeakDistance', 2500); 

Fs = sampling_rate; 
t = 0:1/Fs:0.25; 

C = zeros(size(B),2) 
    for i = 1:numel(B) 
    new_data = data(B(i):(B(i)+200)) 
    y = double(new_data)/max(abs(new_data)); 
    n = length(y); 
    p = abs(fft(y)); 
    f = (0:n-1)*(Fs/n); 
    p1 = p(1:round(length(p)/2)); 
    p1(p1<10) = 0; %add a threshold 
    [~,ind] = findpeaks(p1); %index of where are the peaks 
    C(i,:) = f(ind); 
    hold on 
    plot(f,p,'b',C(i,:),p(ind),'ro') 
    end 

Answer 2

以下が役立つかもしれない、あなたの信号データの1つのFFTからのピークを取得しているようだ、

clear all 
close all 

%load sample data from https://www.dropbox.com/s/zxypn3axoqwo2g0/signal%20%281%29.mat?dl=0 
load('./signal (1).mat') 

%get an FFT and take half 
p = abs(fft(signal)); 
p = p(1:length(p)/2); 

%find peaks and plot 
[pk, loc] = findpeaks(p,'MINPEAKHEIGHT',100,'MINPEAKDISTANCE',100); 

plot(p,'k-') 
hold all 
plot(loc, pk, 'rx') 

のように見え、

ピークのいくつかは単離されているated ...ここ