首記についてご質問させて頂きます。
あるデジタル信号の周波数解析を行う前の最適な窓関数を設計しております。
その際、メインローブとサイドローブを見る必要があると思います。
添付画像(Wikipediaより抜粋)のようなグラフです。
https://ja.wikipedia.org/wiki/%E7%AA%93%E9%96%A2 …
お恥ずかしいお話ですがこれの求め方がよく分かっておりません。
C#のシステム内に組み込む予定なので、求め方はプログラミングを考えています。
データ長Nの実数配列をFFTし、複素数配列に変換するメソッドはすでに実装してあるので、早速、データ長N窓関数配列をFFTした(以下、W[N])のですが、この後どのようにするかが分からず困っております。
W[N]について、パワースペクトル(実部虚部の平方和)、振幅スペクトル(複素数の絶対値)等試しに求めてみましたが、上記Wikipediaのようなグラフにはなってくれませんでした。
ソースコードを下さいとは言いませんので、回答を見て手で実装できる程度の粒度で構いませんので、詳しいお方がいらっしゃいましたら、何卒ご教授の程宜しくお願い申し上げます。
No.1
- 回答日時:
まずWikipediaの図は矩形窓のパワー・スペクトルと書いてあるので、この例は矩形窓ですね。
ついでにパワースペクトルはデシベル(dB)単位なので、図はパワースペクトルを最大値で割って対数を取って10倍した値ですね。
求めたパワースペクトルをp[n],n=0,...,N-1としたとき、
y[n]=10*log10(p[n]/max(p)) n=0,...,N-1
としてyのグラフを描いてください。
No.2ベストアンサー
- 回答日時:
矩形窓の幅を短めに取ってみるとか。
窓の外の点もないと窓にならないから。手元に良いプログラム環境がないのでExcelのFFTで試してみたけど、0,...,0,1,...,1,0,..,0と全512データでやっても少しはそれっぽくなったし。
ただFFT結果そのままだと負のDFT binsは出てこないので、周期性を利用して結果の後半を負のDFT binsとみなすようにしないといけませんね。
ご回答ありがとうございます。
窓の外の点を付け加えましたところ、無事に作成することができました。
大変助かりました。
ありがとうございました!
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Ruby 初心者プログラミング 3 2022/10/12 11:31
- 数学 フーリエ変換後の負の周波数成分の扱いについて 4 2022/09/03 10:18
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- Visual Basic(VBA) 複数ファイルのデータの統合について 12 2022/05/14 12:03
- 物理学 風力発電での音 1 2023/04/16 08:55
- C言語・C++・C# C言語初心者 ポインタについて、お助けください、、 2 2023/03/15 23:50
- C言語・C++・C# このプログラミングの問題を教えてほしいです。 キーボードからデータ数nとn個のデータを入力し、平均値 3 2022/12/19 22:51
- 統計学 Rを用いた「繰り返しがある直交表実験計画法」の分析方法 8 2023/08/01 17:58
- 統計学 t検定について教えてください 2 2023/02/23 16:35
- その他(プログラミング・Web制作) パイソンのプログラミングについての質問です 2 2023/05/22 12:39
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
サンプル数の異なる2群間にお...
-
応答で層別した場合の検定について
-
EXCELにてローパスフィルタを作...
-
心理機能診断をしたのですが、...
-
下の対数表示のグラフから低域...
-
エクセルのグラフから半値幅を...
-
脳波(EEG)の周波数は0.5~60Hzで...
-
検量線の決定係数について
-
極値をもつ時と持たない時、単...
-
ポアソン回帰でのカウントデー...
-
検定統計量の値がマイナス
-
統計学 カイ二乗検定とt検定の...
-
数3の複素数です 12分の25πが12...
-
t検定の使いどころ
-
片対数グラフで…
-
両対数グラフでの直線式の求め方
-
高校 数学 aを実数の定数とする...
-
最小二乗法を反比例の式を元に...
-
t検定の記述の仕方を教えてくだ...
-
パイロットサンプルって何ですか?
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
EXCELにてローパスフィルタを作...
-
サンプル数の異なる2群間にお...
-
心理機能診断をしたのですが、...
-
検量線の決定係数について
-
下の対数表示のグラフから低域...
-
エクセルのグラフから半値幅を...
-
x^2+y^2はどのような分布をする?
-
検定統計量の値がマイナス
-
統計について
-
【統計】有意に「高い」?「低...
-
統計学の基本的なことについて...
-
最小二乗法を反比例の式を元に...
-
心理学の統計について
-
ポアソン回帰でのカウントデー...
-
[Excel] リストからの無作為抽出
-
対数グラフの対数とは・・・?
-
極値をもつ時と持たない時、単...
-
対数目盛の読み方を教えてください
-
溶解度の問題 理科
-
正規分布について
おすすめ情報
ご回答ありがとうございます。
私もそう思い試してみたのですが、うまくいきませんでした。
矩形窓はw[n]のnによらず値が1ですので直流波形と言えると思います。
これに対しFFTを実行しパワースペクトルを求めると、直流成分にのみスペクトルが立ち、それ以外は0になってしまうのです。
そもそもFFTを実行するのが間違っているのでしょうか…。