Matlab FFT 周波数の最大値、最小値を求める方法

Question

音声ファイルを読み込み、Matlabで最大値、最小値を求めるプログラムが知りたいです。
音声ファイルのデータをfftしてスペクトルにすることは出来ました。
[sig,Fs] = audioread('.wav');
​
T=1/Fs;​
L=length(sig);​
t=(0:L-1)T;​
​
Y=fft(sig);​
P2=abs(Y/L);​
P1=P2(1:L/2+1)​
P1(2:end-1)=2P1(2:end-1)​
f=Fs*(0:(L/2))/L​
​
plot(f,P1)​
​
title('task1')​
xlabel('f(kHz)')​
ylabel('|P1(f)|')​
​
ここまでは何とか出来たのですが、
Ma=max(f)
Mi=min(f)
とすると音声ファイルを変更しても毎回同じ数字(11025、0)と出力されます。
どなたかこの原因・解決法ご存知の方居られましたら宜しくお願いします。

Answer 1

いろいろな実装上の都合やらなんやらで、現実的にはそうなってしまいますね。
まずは、あなたのいう周波数の最大値、最小値をきっちり定義しましょう。
なにをするための最大／最小値なのか、どういうときの値をそれとするのか、をはっきりしないと、どうやったところで最小値は０、最大値は無限大ということになります。

なにも考えずにFFTかけて結果は出ますが、意味のある結果を出そうとするなら難しくなります。
とりあえず、ナイキスト定理、窓関数、あたりでぐぐって、出てきたものを一通り読みましょう。
道のりはまだまだ遠いです

Answer 2

恐らくこういうことがやりたいのでしょう。

[~,I]=max(P1);
Ma=f(I);
[~,I]=min(P1);
Mi=f(I);

Comments

[Nogis]

回答ありがとうございます。
求めていた値が出力されました！
理解が浅く申し訳ないのですが
Ma=f(I )とは何を意味していますか？

[WathMorks]

Matlabのドキュメントをみてください。