Pythonでガウス分布を持つノイズの作り方と調整方法

  • このエントリーをはてなブックマークに追加

信号処理はノイズと切っても切れない縁があります。様々なモデル同定やテストに意思入れの無い信号としてノイズを用いることがあります。ノイズを用いる際は「どんな」ノイズか素性を知っておく必要があります。

こんにちは。wat(@watlablog)です。
ノイズは何も邪魔な存在だけではありません。
ここでは信号処理分野では代表的なガウシアンノイズの生成とその特徴を説明します

信号処理として必要とされるノイズ

ノイズとは?

信号処理分野で言うノイズとは、計測した信号以外の不要な雑音を意味しており、世間で頻繁に行われている「ノイズ除去」に代表されるようにできるだけ無い方が良いものとされています。

ノイズを作り出す理由とは?

本来ノイズとは邪魔な存在であり、除去する対象でもあるわけですが、本ページではあえてそのノイズを作り出そうとしています。

ノイズが必要が理由は様々ありますが、素性の確かなノイズは工学的にとても有用であることが、ノイズをあえて作り出す1つの理由です。

工学の世界ではノイズを用いて対象となるモデルの特性を得る実験や解析がよく行われています。

代表的な物に、構造物の共振調査があります。簡単に説明すると、あるモデル(実際の構造物、もしくは計算モデル)の任意の点を加振し、ある点の応答を振動として計測する実験です。

この計測したデータと入力した加振力の比をとることで、構造物の固有振動数やどんな形で動いているかといった振動特性を知ることができます。

ここで、加振力は調査したい周波数帯域までできるだけクセを持たずに、フラットな特性を持つことが望ましいとされています。

適当に手で揺らせば良いわけじゃないんだね。
じゃあある程度こちらで意図したノイズを用意すれば良さそうだ。

ガウシアンノイズを作ろう!

Advertisements

ガウシアンノイズとは?

このページでは数あるノイズの中で、一般的に扱いが容易なガウシアンノイズを作る方法を学びます。

ガウシアンノイズとは、信号のデータ点が正規分布(ガウス分布)を持つノイズのことを指します。

ガウシアンノイズのパラメータ

正規分布を持つノイズは一つでは無く、「平均値」と「標準偏差」というパラメータで調整可能です。

以下に平均=0, 標準偏差\(\sigma\)=1の標準正規分布の例を示します。

正規分布は以下のコードで簡単に書くことができます。

こちらは参考までにグラフ設定です。

ガウシアンノイズ生成のコード

ガウシアンノイズは以下のコードで生成することができます。平均値loc, 標準偏差scaleでパラメータ設定します。データのサイズsizeでガウシアンノイズを作る時間範囲を決めます。

ヒストグラムとFFTで確認!

生成したノイズ信号が意図したものになっているか、信号処理をして確かめてみましょう!

尚、FFTの方法に関しては以下の記事をご覧下さい。

PythonでFFT!SciPyのFFTまとめ

標準偏差\(\sigma\)を振ってみる

パラメータの一つである標準偏差を3段階に振った結果を以下の図に示します。

若干ヒストグラムがいびつな形をしていますが、FFT結果から周波数特性はほぼフラットになっていることがわかります。
※周波数波形で「完全に」フラットな波形にするためには、ディラックのデルタ関数を実現させるか、サイン波を周波数スイープさせるかしかないかも。

平均を振ってみる

続いて、平均を振ってみた結果を以下の図に示します。
今度はヒストグラムの山がずれていきましたね。

まとめ

本ページでは、本来邪魔な存在であるはずのノイズにも、工学的に利用価値があることを説明しました。

また、ガウシアンノイズの概要と、その生成方法を示し、パラメータによる波形変化を観察しました。

これでまた1つ信号処理としての武器を得たぞ!ノイズの世界は本当はもっと奥が深そうだけど、これでまずは道具として使えるようになったかな!

Twitterでも関連情報をつぶやいているので、wat(@watlablog)のフォローお待ちしています!

  • このエントリーをはてなブックマークに追加

SNSでもご購読できます。

コメント

コメントを残す

*