matlabにてFFTしたのち周波数平均, 周波数中央値を求めたい

Question

前提・実現したいこと

matlabにてあるデータをFFTしたのち周波数平均, 周波数中央値を求めたいです。
matlabにはそれぞれ求めるmeanfreq, medfreqがあるのですが, それを用いても欲しい値が出ません。
引数に問題があるのか, 使い方に問題があるのかご教授いただきたいです。よろしくお願いします。

ソースコード

function [MeanFreq, MedFreq] = FFT(X, num)
%X:FFTしたいデータ(●行1列の行列), num:Xのデータ数

Fs = 1000;            % Sampling frequency 
L = num;             % Length of signal

Y = fft(X);
p2 = abs(Y/L);
p1 = p2(1:L/2+1);
p1(2:end-1) = 2*p1(2:end-1);

f = Fs*(0:(L/2))/L;

plot(f,p1);

title('振幅スペクトル')
xlabel('f(Hz)')
ylabel('|p1(t)|')

MeanFreq = meanfreq(p1,Fs); % 平均周波数
MedFreq = medfreq(p1,Fs); % 中央周波数


end

Answer 1

周波数平均 - MATLAB meanfreq - MathWorks 日本

freq = meanfreq(x) は、時間領域信号 x のパワー スペクトルにおける正規化周波数の平均 freq を推定します。
freq = meanfreq(x,fs) は、周波数平均をサンプルレート fs で推定します。

というように、

時間領域信号 x

なのでFFTしたもの(周波数領域信号)を入れるのではありません

Comments

[alastashim]

回答ありがとうございます。
もう一つ質問があり恐縮なのですが、FFTした周波数の平均などを求める関数は存在しないということですか？

[ozwk]

ぶっちゃけmatlabとか使ったこと無いので知りません

リファレンス読んで、ないならないんじゃないですかね？

ただ、多分なんですけど、あなたのやりたいことは、Y=fft(X)とせずにそのままXをmeanfreqにつっこめばいいだけな気がします

[alastashim]

ありがとうございます