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;
---------------------------------------------------------
以上です。
何か足らない情報などがありましたらお申し付けください。
よろしくお願い致します。
No.1ベストアンサー
- 回答日時:
dB値は相対的なものですから、基準となる信号の時系列データを作って、それを同じプログラムで処理したときのパワースペクトルの電力のd
B値を差引いてやれば、基準となる信号を基準としたパワースペクトルのdB値が得られるかと思いますが如何でしょうか?info22様
回答して頂きありがとうございます。
基準信号として片側振幅1のcos波を作成して上のプログラムを実行したところ、約60dBの値になりました。
これを基準値として考えていきたいと思います。
どうもありがとうございました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 フーリエ級数展開をExcelのFFTでシミュレートする 5 2023/07/03 22:02
- 数学 Pythonのspecgramに関する質問です。 1 2022/12/04 19:18
- 物理学 風力発電での音 1 2023/04/16 08:55
- 数学 「FFTの基本は、DFTはサンプル数Nが偶数なら 2つのDFTに分解できるということ。 分解するとD 3 2022/03/31 21:01
- Excel(エクセル) VBAで組み合わせ算出やCOUNTIFSの処理を高速化したいです。 4 2022/04/07 02:38
- レトロゲーム Fftというゲームでnpcとかあったようななんでしたっけよろしくお願いしますm(_ _)m 2 2023/04/16 15:50
- その他(Microsoft Office) パワークエリの複数ファイルのデータ統合について 3 2022/07/14 17:06
- Excel(エクセル) Power Query でのデータの一括修正について 2 2022/05/10 02:00
- 数学 5組のデータ (xx,yy) = (1.0 , 1.1), (2.0 , 1.0), (3.0 , 2 2023/01/09 16:52
- Java Java 年数計算 3 2023/01/28 10:52
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
三次元のフーリエ変換
-
1/3オクターブバンドについて質...
-
FFTのDC成分って、なんで大きく...
-
どなたか、S45Cの電気伝導率(S...
-
Igorの使い方について
-
近似直線について
-
分解能の値が大きいと小さいも...
-
隣接平均と移動平均
-
静的・動的の意味
-
MOSFETのgm-Vgs特性
-
1/f ゆらぎ スペクトル解析の...
-
半値幅の測り方
-
スペクトル解析と周波数解析の違い
-
フーリエ変換、FFT演算について
-
時系列に不等間隔なデータのス...
-
指数関数のカーブフィッティング
-
有効数字…?
-
MTFを求める際に、平滑化や隣接...
-
減衰振動のグラフが書けない…
-
savitzky-golayのフィルタについて
おすすめ情報