MATLABを用いたFFT

質問

現在40kHzの正弦波をMATLABを用いてFFTをかけましたがしかし思った結果になりません。作成したmファイルの内容は以下のとおりです。

data = xlsread('test001.xls','Sheet1');　　　%excelファイルの読み込み
data = data(:,2);　　 %２列目を選択
Y = fft(data);　　　 %fft
Pyy = Y.*conj(Y)/2500; %スペクトルを求める
f = 40000 * (0:2499)/2500;　　　　　　　　　 %範囲を指定
plot(f,Pyy(1:2500))　　　　　　　　　　　　　%グラフを作成

本当に困っているので、解答お願いします。

回答No.1

rabbit_cat

何が思ったとおりにいかないんでしょうか。
とりあえず、そのプログラムで何かは表示されるんでは、と思いますが。
横軸があってない、ってことですか？
＞現在40kHzの正弦波
とのことですが、そのプログラムを見る限り、サンプリング周波数40kHzで取り込んだデータとして扱っているようですが、それで正しいのですか？サンプリング周波数はいくつなのでしょう？