findpeaksの利用方法

findpeaksで局所最大点を求め、最大値の1/2の値までを得たいときのプログラムができません。

y = [1 9 5 7 5 6 1 3 2 5 1];

x = 0:numel(y)-1;

plot(x,y)

[py,px] = findpeaks(y)

pyは局所最大点のyの値

pxは局所最大点のxの値

この例ですと

py = [ 9 7 6 3 5]

px = [ 2 4 6 8 10]

となります。

そして求めたい結果がmax9の半分(4.5)以上の値のyとx

py2 = [ 9 7 6 5]

px2 = [2 4 6 10]

回答のほどよろしくお願いします。

Answer 1

[Shunichi Kusano]

findpeaks関数では様々なプロパティを指定して、いろんな条件付けの元でfindpeakしてくれます。

今回の場合ですとMinPeakHeightプロパティを指定することで解決できそうです。

[py,px] = findpeaks(y,'MinPeakHeight',max(y)/2)

[貴大 石川]

解決できました。ありがとうございます。

ちゃんと見るとfindpeaksの下のほうに掲載されていましたね・・・私の確認ミスでお手数おかけしました。

この度は回答ありがとうございました。

[Shunichi Kusano]

いえいえ、お役に立てて良かったです。