- 締切済み
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)) %グラフを作成 本当に困っているので、解答お願いします。
- hikaru1003
- お礼率0% (0/2)
- その他（プログラミング・開発）
- 回答数1
- ありがとう数6
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- rabbit_cat
- ベストアンサー率40% (829/2062)
何が思ったとおりにいかないんでしょうか。 とりあえず、そのプログラムで何かは表示されるんでは、と思いますが。 横軸があってない、ってことですか？ ＞現在40kHzの正弦波 とのことですが、そのプログラムを見る限り、サンプリング周波数40kHzで取り込んだデータとして扱っているようですが、それで正しいのですか？サンプリング周波数はいくつなのでしょう？
関連するQ&A
- MATLABによるFFTと位相スペクトルについて
純粋な正弦波信号 y=sin2πft ただし周波数f=500MHz をサンプリング周波数2.5GHz,サンプリング数128,256,512,1024の四パターンMATLABでFFTし、パワースペクトルと位相スペクトルを求めました。 結果、パワースペクトルは予想通り500MHzでピーク値をとったのですが、位相スペクトルは各サンプリング数ごとに全く違う値を取りました。 というわけで、 （１）理論上、正弦波の位相スペクトルはどのような値になるのか （２）なぜサンプリング数ごとに位相スペクトルの結果が大きく変わったのか この二点をご教授ください。よろしくおねがいします。
- 締切済み
- 大学・短大
- MATLABによるFFT
MATLABのfft関数を使用して，フーリエ変換してPSDを算出すると，スペクトルの周波数領域での積分値が元データの二乗平均値と一致しないのですが，原因は何が考えられるでしょうか？FFTの出力がフーリエ係数であれば一致するはずなのですが．
- ベストアンサー
- 数学・算数
- パワースペクトルの0dBの値について
matlabを用いて時系列データのパワースペクトルを計算、対数表示でグラフにプロットしています。 matlabのfft関数のヘルプに記載されていた例を参考に、以下のような自作の関数を組んで計算しているのですが、2つ質問があります。 1. この場合、パワースペクトルの0dBは、時系列の値ではいくつになるのでしょうか？ 0dB = 1でしょうか？ 2. 関数の中において、power_fftをN_fft（高速フーリエ変換時のデータ長さ）で割っている事の理由がわからなかったのですが、なぜ行っているのでしょうか？ ------------------------------------------------ 以下、自作のパワースペクトル表示関数を示します function Power_db(Data_xx) fs = 1000; % サンプリング周波数 [Hz] N_fft = 2^12; % 4096個 xx_fft = fft(Data_xx,N_fft); power_xx = xx_fft.* conj(xx_fft) / N_fft; frequency = fs*(0:2047)/N_fft; % (0:2047)で500Hzまで plot(frequency,20*log10(power_xx(1:2048))) xlabel('Frequency [Hz]'); ylabel('Power Spectrum [dB]'); grid on; --------------------------------------------------------- 以上です。 何か足らない情報などがありましたらお申し付けください。 よろしくお願い致します。
- ベストアンサー
- 物理学
- FFT フーリエ変換 について
FFTについて質問があります。 あるシミュレーションのために ・サンプリング周波数48kHz ・量子化24ビット ・1kHzのサイン波 のデータを作って、それをFFTにかけてみると 裾野が広がり、フロアの高いスペクトラム となってしまいました。（データは4096ポイント、Hanning窓） ポイント数が少ないことによる劣化だと考え、データポイントを 65536ポイントで取ってみると綺麗なスペクトラムとなっており、 データの作り方としては間違っていないのではないかと考えています。 一般的にFFTはデータポイントが多い方が、FFTをとったときに、 裾野広がらず、シャープなスペクトルを得ることはわかっています。 が、そこを何とか少ないポイント数で綺麗なFFTを得られるような 入力データは作れないものでしょうか？ 出所は不明ですが、私の持っている別のデータを入力として用いると 少ないポイントでも綺麗なFFTを得ることが出来ているので、 何かデータを作るコツのようなものがあると思って投稿させて いただきました。 どなたか詳しい方がいらっしゃいましたら、ご教授お願いします。
- ベストアンサー
- 数学・算数
- matlab計算について
matlab独学で進めていますが、時間がないのと頭が足りないので 質問させていただきます。 p1=0.1; p2=2.663*10^6; D= 1:10:1000; A=(log(D)+12.7)/1.1; f(D)=(1/power((2*pi),0.5)*1.1*D)*exp(-0.5*power(A,2)); Y=(p/pash*4/3*pi*power((D/2),3))*f(D) plot(D,Y) Ｘ軸がＤ、Ｙ軸がf(D)のグラフを作成したいです。 内部行列の次元は一致しなければなりません。 とエラーが出ます。このエラーの原因を教えてください。 間違いだらけかもしれませんが...ご指摘よろしくお願いします。
- ベストアンサー
- その他（プログラミング・開発）
- FFTを使って信号から特定の周波数を抜き去る方法
FFTを用いてある時系列データから特定の周波数の成分を抜き去ろうと思っていますが、 うまくいかないのでお伺いします。 調べてみると、スペクトルを出すとき、例えば同じ振幅で周波数の異なるの正弦波を １０個入れても全て同じ振幅にならず、振幅の大きさにばらつきが出ていて周波数成 分の成分の推定がうまくいっていないようです。FFTの2乗のデータを小さな周波数区 間で積分してスペクトル密度にしても、先ほどの１０個の正弦波の振幅がそろわない のと、周波数推定が悪くなりうまく、引き去れなくなっています。 よくある状況だと思いますが、どういう工夫が必要なのでしょうか。 ちなみに、FFTのルーチンはnumerical recipeとFFTWの二つで試しました。
- ベストアンサー
- C・C++・C#
- Matlab 2行の行列をプロットしたいのですが
こんにちは、Matlabを使い始めのものです。 二行の行列（2x10とか2x40など）にについて、一行目をx座標、二行目をy座標としてグラフ化したいのですが、うまくいきません。 自動的にx座標は列の番号となり、二つのグラフが現れてしまいます：y座標に一行目、二行目、両方がプロットされてしまいます。どうにか、一行目をx、二行目をy座標としてグラフ化できないでしょうか。おそらく、とても単純なことなのかもしれないのですが、どうしてもわからず、とても歯がゆい思いをしています・・・。どうかご教示ください。 何をしたいかといいますと、音声データをMatlabで取り込んだのですが、振幅データがずらっと並んだ、1x500000とかの行列が得られます。 これを、プロットすると、x軸はサンプル番号が並ぶわけですが、本当はx軸は時間であって欲しいため、音声データの行列にひとつ行を足して、時間軸を作ってあげました。つまり、今手元にあるのは、2x500000の行列で、一行目が時間、二行目が振幅となっています。これをx軸が時間軸、でy軸が振幅となるようにグラフ化してあげたいと考えています。 文章がわかり辛くて申し訳ないのですが、何卒よろしくお願いします。
- ベストアンサー
- その他(ソフトウェア)
- matlabのｙ軸を2つ利用したグラフについて
matlabのグラフ作成について質問です。 y軸を左右に表示し、グラフを作成したいのですが、データが6個あります。 2個であればplotyyを使えばいいと思いますが、今回は、データのうち4つを左の軸、2つを右の軸にしたい場合はどのようにグラフを作成したらいいでしょうか？ わかる方、よろしくお願いいたします。
- 締切済み
- その他（プログラミング・開発）
- FFTデータ数が少ないと直流成分が大きく変動する？
FFTの振幅スペクトルから実効値算出のために単一正弦波をFFTしてみると、(1)のように計算データ点数が少ない場合、窓関数の位置によって「0Hz直流成分」のFFT結果値が大きく変動してしまいます。（(2)(3)は窓関数がどの位置でも問題ありません。） この(1)（1サイクル64点＝FFTデータ数64個）のような計算の場合、どのように実効値算出すればよいでしょうか？ また、正常な場合でもこれらのサンプル間隔(0、30Hz、60Hz等)の"すき間の部分”は考慮しなくても、全て計算結果に出力されているということでよいでしょうか？ 単一正弦波(60Hz) （(最大値5388.9v=実効値3810.5v)、サンプリング3.84KHz=1サイクル64点） （実効値換算式＝前半の√(実数^2＋虚数^2)の合計／(FFTデータ数／2)／√2） (1)FFTデータ数 64個⇒0Hz(0.00～1905.2)???、60Hz(1905.2),120Hz(952.6)、180Hz(0.0)、 ⇒合計( ??? ) (2)FFTデータ数128個⇒0Hz(0.0), 30Hz(952.6),60Hz(1905.2),90Hz(952.6), 120Hz(0.0)、 ⇒合計(3810.4) ≒実効値 (3)FFTデータ数256個⇒0Hz(0.0),15Hz(0.0),30Hz(0.0)、 45Hz(952.6),60Hz(1905.2),75Hz(952.6), 90Hz(0.0), ⇒合計(3810.4) ≒実効値
- 締切済み
- 電気・電子工学
