No.8ベストアンサー
- 回答日時:
#2のものです。
「スペクトルが一様になりません。」ということですが、念のため私もやってみますと、スペクトルは一様になります。
思うに「一様」とか「白色」とかいう意味を誤解していらっしゃるのでは?
全く平坦なスペクトルを得たいのなら入力データ列の一点だけが有値で他は全てゼロの値とすればいいです。いわゆるデルタ関数。そのフーリエ変換は全くフラットです。こういうスペクトルをイメージしておられるのではないでしょうか。
しかし!!!! 今はノイズつまりランダムな数列を扱っているのですよね。そのフーリエ変換はやはりランダムな凹凸のスペクトル分布です。どこかの周波数が高いとか低いとか起きて当たり前、偶発事象ですから。ところがランダムな数列を新たに発生させて同じことをやるとその周波数は今度は必ずしもピークでもボトムでもなく、ピークやボトムの位置はまた別の周波数になっているでしょう。これを何度もやってスペクトルの平均をとってごらんなさい(やるのはめんどくさいからイメージするだけでも)。回数を重ねるごとにだんだん平坦になります。これが「ホワイトノイズ」です。
念のため、元のデータ列に移動平均を施してそのフーリエ変換をみてください。高周波部分は振幅減少しているでしょう。ホワイトノイズではなくなったのです。
ご回答ありがとうございます。
imoriimoriさんのおっしゃる通り、ホワイトノイズのパワースペクトルはフラットになるのだと思っていました。違うんですね…。
確かに、rand()で作成した波のパワースペクトルはランダムな凹凸のスペクトル分布になりました。これが一様ということなんですね。
これがホワイトノイズのスペクトル分布だったんですね。
No.7
- 回答日時:
皆さんの回答に重複するところもありますが、
1)"white" になる為には、時系列の個々の値が独立ランダムに決まっていれば十分です。振幅の分布に関しては、ANo.2さんの一様乱数でもかまいませんし、ANo.3さんのポアソン分布でもかまいません。任意です。時系列の個々の値達に関連がなければ、"white"で、すべての周波数において同じエネルギになります。
2)ガウシアン分布(正規乱数)を得る方法としては、ANo.4さん紹介の、and()+rand()+rand()+・・・ の他に、ANo.5さんの方法も便利です。
NORMSINV(RAND()) とセルに書込むだけのようです。
その他、一様乱数から正規乱数へ変換する有名な手法としては、ボックス・ミューラー法があります。
http://mathworld.wolfram.com/Box-MullerTransform …
[0,1)の一様乱数の一組の、x,y、に対して sqrt(-2*ln(x))*cos(2*π*y) は平均零、分散1の正規乱数を与えます。セルに書込むのは、
SQRT(-2*LN(RAND()))*COS(2*PI()*RAND()) です。
ご回答ありがとうございます。
rand()+rand()+rand()+・・・の方法
NORMSINV(RAND())の方法
で波を作り、パワースペクトル解析したのですが、
一様の分布になりません。
何か間違っているのでしょうか…。
No.6
- 回答日時:
#2のお礼を見ると、スペクトルとサンプリング周波数の関係が微妙なんだろうと思います。
とりあえず、ホワイトノイズというのは、全ての周波数でパワースペクトルが一様ということですが、これは、ウィナー・ヒンチンの定理より、自己相関が時間0以外では全て0ということと同値です。
つまり、ホワイトノイズをw(t)と書いたとき、t0≠t1のとき、w(t0)とw(t1)が独立であればいいことになるます。w(t0)の確率分布が正規分布であるか一様分布であるかによらないことは#4さんの仰るとおりです。もっと言えば、w(t0)とw(t1)で従う確率分布自体が異なってもかまいません。
有限のサンプリング周波数でサンプリングする場合は、ナイキスト周波数までは、スペクトルが一定になります。
シミューレーションとかで、あるシステムにホワイトノイズを与えたときの挙動を調べたいといった場合には、対象システムの最も早い周波数のさらに100倍程度のサンプリング周波数でノイズを与えないと、正確な結果が得られません。
ご回答ありがとうございます。
rabbit catさんがおっしゃるようにシュミレーションであるシステムにホワイトノイズを与えたときの挙動を調べたいんです。
これまでホワイトノイズなどに関係するような分野ではなかったので、ホワイトノイズをシュミレーションするときにどうやって作ればいいのかがわからなかったのです。
No.5
- 回答日時:
使えるかどうか分かりませんが、エクセルなら分析ツールで一様、ガウシアン、ポアソンなどの乱数を発生させることができます。
(Excel2000しか知りませんが、多分、新しいバージョンでも・・・・)一度試してみては?
No.4
- 回答日時:
私の理解するところでは。
。。まず第一に
確率密度分布がどうであるか(ポワソン or ガウシアン or 矩形一様 or その他)と、スペクトルが白色か否かとは別の話です。
エクセルのrandが発生する乱数は矩形一様ですが、矩形一様であっても、そのスペクトル(フーリエ変換)は特にどのへんの周波数で強大とか微弱とかいうことはなく、白色となります。
矩形一様は自然界にあまり見あたらないので、それっぽい乱数が欲しいならガウシアンが手頃かと思います。最も自然界にあるっぽいものといえばポワソンがいいかもしれませんが、めんどくさいです。それに多くの場合本質的現象はポワソンであっても大数近似で結局ガウシアンで扱うことが多いです。。
矩形一様乱数の加算でほぼガウシアンを作れます(中央極限定理)。例えば#2様のB列をrand()+rand()+rand()+rand()+rand()+rand()-3 とかにすればほぼガウシアン分布になります。
第2に
次にご質問でお求めの「白色」の「波」ということですが、ゆるやかな波(不規則だけど波として見えるような)をイメージしたらそれはもはや厳密な意味で白色ではありません。高周波域がカットされている雑音ですから。厳密に白色といえば、周波数帯域∞です。ところがこのような白色雑音は自然界はおろか人工現象にも存在しません。周波数帯域∞ということはエネルギー∞、あり得ないのです。
最後に
そういう厳密な白色雑音は要らないのだ、波っぽいものが欲しいのだ、ということでしょう。ならば#1様や#2様の言われるようにエクセルで乱数を並べ、その乱数列を軽く平滑化(running smoothing)してグラフ化してごらんなさい。「帯域制限された白色雑音」(変な言い方ですが)の波が見えるはずです。こういう波が欲しいのではないのでしょうか。
このやりかたで作った「白色雑音の波」はフーリエ変換すると高域がちょっと不自然に見えるかも知れません。もっと自然なものが欲しいというのであれば、どのような現象を模擬するものかによりけりで、そのモデルに従い乱数発生や平滑化とか特定周波数域の強調とかをしなければいけないでしょう(それは白色からほど遠い雑音を作ることになりますが)。というよりその現象の実データ例を探すことが最善でしょう。(そもそも地震波は白色なんでしょうか。私には白色からほど遠い波形に見えますが。)
ご回答ありがとうございます。
imoriimoriさんのおっしゃるようにrand()+rand()+rand()+rand()+rand()+rand()-3とB列に入力し、それをパワースペクトル解析をしたのですが、スペクトルが一様になりません。
なんででしょう…。
また地震波に関しては有名な地震波はインターネットなどで手に入れられるので、ホワイトノイズに関してもそのようなことがあるのかと思って、地震波を例に出しました。地震波が白色というわけではないんです…。
No.3
- 回答日時:
ホワイトノイズの時系列信号を作るには、Poisson分布に従う乱数を
発生させる必要があります。
ExcelのRAND()関数を単純に使っただけではできませn。
ExcelのRAND()関数は「一様乱数」なので、この乱数系列だけで
作った時系列信号の周波数特性は「ホワイト」にはなりません。
この検証は、先の方のような方法で作成したデータ系列を、同じく
Excelに備わっている「フーリエ解析」にかけると一発で分かって
しまいます。
Poisson乱数の発生方法は、そう簡単ではありませんが、ウェブ検索
をすると下記のような、Poisson乱数発生アルゴリズムについて書かれた
サイトが見つかりますので、」参考にされると良いでしょう。
http://hata.cc/docs/sim/poisson.html
参考URL:http://hata.cc/docs/sim/poisson.html
No.2
- 回答日時:
1さんがおっしゃるように
Excelの関数rand()を使えばできます。rand()は0から1までの乱数を返す関数です。
例えばA列を時間軸、B列を振幅軸にしたいなら
A列には
0
=A1+1
=A2+1
=A3+1
・
・
・
といったセルをコピー・貼りつけで作り、
B列には
=rand()-0.5
=rand()-0.5
=rand()-0.5
・
・
・
といったセルをこれもコピー・貼りつけで作ればいいわけです。
この場合最大周波数0.5 Hz、最大振幅0.5 のホワイトノイズが作れます(時間軸の単位が秒の時)。
ご回答ありがとうございます。
paddlerさんがおっしゃることを参考にするとrand()を使ってもできないように思えるのですが…。
すべての周波数において同じエネルギーになるのでしょうか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 統計学 お世話になっています. x軸は時間(期間)y軸はある値に対する2つのグラフ比較をしますが、私個人の考 2 2023/03/30 11:42
- 地球科学 津波のデータ(東日本大震災) 2 2023/06/02 04:23
- Wi-Fi・無線LAN モバイルWi-Fiとデータ通信の電波について 電波の強さはデータ通信に比べてモバイルWi-Fiの方が 9 2022/05/06 13:01
- 英語 強調構文だと思うのですが。。。 5 2022/10/22 13:31
- Excel(エクセル) エクセルで沢山のレコードの最後に追記するには? 7 2023/04/10 13:27
- SoftBank(ソフトバンク) 機種変更のSIMについて疑問です。 ソフトバンクでiPhoneXRからiPhone14へ機種変更をし 3 2022/10/01 21:34
- Android(アンドロイド) IIJmio 使用しています 職場は電波が入りにくいですが ドコモは一番電波が良いですか? 3 2022/09/27 13:36
- iPhone(アイフォーン) iPhoneのデータ通信、SIMカードについて質問です。 海外で使うためにSIMカードを変えたのです 2 2022/05/14 00:10
- 英語 "In a considerable effect to those..." の In は正しいか? 3 2023/06/22 07:58
- Wi-Fi・無線LAN 電波について 電波って音声データも飛び交ってるんでしょうか?電波あってたら建物の中でも他の国のワンセ 8 2023/08/23 20:02
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
スペクトル強度について
-
水銀灯のスペクトルは、なぜ4...
-
線スペクトルと連続スペクトル
-
"spectrum"と「スペクトラム」...
-
Arのスペクトル
-
ラマン測定におけるネオンラン...
-
偏光レンズの仕組み
-
レーザーの線幅
-
写真の数式の意味が分からない...
-
蛍光の輝度を測定する方法
-
フォノンレプリカとは何ですか
-
波長の求め方が分かりません
-
white noiseの波のデータ
-
黒体放射、蛍光、シンクロトロ...
-
メチルオレンジについて
-
o-ニトロアニリンとp-ニトロア...
-
ブランク値って・・・
-
ジベンジリデンアセトンの再結...
-
ポリエチレンが溶解可能な溶媒...
-
わかる方、この解説おねがいい...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報