様々な移動平均・レポート-No.3
移動平均のレポートでは、表計算ソフトを用いて、移動平均の範囲のとり方などを変えながら、平均をとった結果に、どう影響するのかを考える。
Excel で様々な移動平均の式を入力
表計算ソフトに、移動平均の式を入力する際には、以下の資料を参考にせよ。
上図のB4〜E4にできた移動平均の式は、B5以下にコピーすればいい。
レポート内容
以下のような移動平均を、Excel にて計算し、その結果の違いについて考察せよ。
移動平均で用いる点の数は、自分の出席番号の3の余りによって、条件を与える。
| 計算方法 | 出席%3=0 | 出席%3=1 | 出席%3=2 |
|---|---|---|---|
| 単純移動平均(狭) | n=3 | n=4 | n=5 |
| 単純移動平均(広) | n=6 | n=8 | n=10 |
| 過去の値による単純移動平均(狭) | n=3 | n=4 | n=5 |
| 過去の値による単純移動平均(広) | n=6 | n=8 | n=10 |
| 加重移動平均(狭) | n=3 | n=4 | n=5 |
| 加重移動平均(広) | n=6 | n=8 | n=10 |
| 指数移動平均(基本) | α=1/2 | α=1/2 | α=1/2 |
| 指数移動平均(広) | α=1/3 | α=1/4 | α=1/5 |
平均に用いる値は、以下のデータとする。
-
- 2018-06-05-wave.csv 時刻tとx(t)のコンマ区切りファイル
- 時間の遅延がわかるような波形を用い、考察してあることが望ましい
1枚のグラフの中に、元波形+8波形=9本のグラフを記載すると、グラフの内容が分かりにくいので、複数のグラフ結果で図示すること。
プログラミングが得意な人は、上記をExcelで処理するのではなく、C言語にて移動平均を計算し、結果をExcelに取り込んでグラフとして表示することにチャレンジしてほしい。
提出レポートに、全データの計算結果は不要です。動作の確認の意味で、先頭10点ほどの計算結果をつけてください。
様々な移動平均
波形処理をハードウェアで行うかソフトウェアで行うか
組込み用の小型コンピュータが普及する以前は、このような波形に対する処理を行う場合には、電子回路的に波形のフィルタリングを行っていた。しかし電子回路的な方法では、回路の特性が変化してうまく処理ができなくなった場合に、回路の組み換えが発生するかもしれない。ただし回路の変更は基板の作り直しが必要であったりすることから、コストがかかる。
一方、A/D変換機能を内蔵した組込み用小型コンピュータも極めて安価になってきた。
こういったコンピュータの普及で、最近ではアナログ値をコンピュータに取り込んでデジタルの数値的な処理で余計な情報を取り除く。この方法であれば、プログラムを変更するだけなので、安価に変更が可能となる。ただし、アナログ値をA/D変換するのには時間がかかることから、周波数の高い信号には不向きである。
単純移動平均
前回説明を行った単純移動平均は、時刻ｔの平均を、その前後のデータで平均を求めた。この方式は、実際には与えられた波形のデータを全部記録した跡に、単純移動平均をとる場合に有効である。
しかし、時々刻々変化する測定値の平均をその都度使うことを考えると、上記の方法は、未来の測定値
を使っていることから、現実的ではない。
#define NS 3
int x[ SIZE ] ; // 入力値
int y[ SIZE ] ; // 出力値
for( int t = NS ; t < SIZE-NS ; t++ ) {
int s = 0 ;
for( int i = -NS ; i <= +NS ; i++ ) // 2*NS+1回の繰り返し
s += x[t+i] ;
y[t] = s / (2*NS + 1) ;
}
過去の値だけを使った移動平均
そこで、過去の値だけで移動平均をとることも考えられる。
この、単純移動平均と、過去の値だけを使う単純移動平均を、適当な測定値に対して適用した場合のグラフの変化を Excel によってシミュレーションした結果を以下に示す。
しかし、このグラフを見ると、波形後半の部分に注目するとよく分かるが、過去の値だけを使った移動平均では、測定値が立ち上がったのを追いかけて値が増えていく。これでは移動平均は時間的な遅れとなってしまう。
for( int t = NS ; t < SIZE ; t++ ) {
int s = 0 ;
for( int i = 0 ; i <= NS ; i++ ) // NS+1回の繰り返し
s += x[t-i] ;
y[t] = s / (NS+1) ;
}
加重移動平均
過去の値を使った移動平均では遅れが発生する。でも、平均を取る際に、「n回前の値」と「現在の値」を考えた時、「その瞬間の平均値」は「現在の値」の方が近い値のはず。であれば、平均を取る時に、「ｎ回前の値は少なめ」「現在の値は多め」に比重をかけて加算する方法がある。
for( int t = 3 ; t < SIZE ; t++ ) {
int s = x[t] * 3
+ x[t-1] * 2
+ x[t-2] * 1 ; // 数個の移動平均だし、
y[t] = s / (3+2+1) ; // ループを使わずに書いてみる。
}
この様に、過去に遡るにつれ、平均をとる比重を直線的に小さくしながら移動平均をとる方法は、加重移動平均と呼ばれる。以下にその変化をExcelでシミュレーションしたものを示す。
指数移動平均
ここまで説明してきた、単純移動平均や、加重移動平均は、平均をとる範囲の「過去の値」を記憶しておく必要がある。広い時間にわたる移動平均をとる場合は、それに応じてメモリも必要となる。これは、組み込み型の小型コンピュータであれば、メモリが足りず平均処理ができない場合もでてくる。
そこで、荷重移動平均の重みを、は、100%,
は50%,
は25%… というように、過去に遡るにつれ、半分にして平均をとる。
しかし、以降の項で、
を使うと以下のように書き換えることができる。
for( int t = 1 ; t < SIZE ; t++ ) {
y[t] = ( x[t] + y[t-1] ) / 2 ;
}
この方法であれば、直前の平均値を記録しておくだけで良い。このような移動平均を、指数移動平均と呼ぶ。
ここで示した指数移動平均は、過去を遡るにつれとなっているが、これをさらに一般化した指数移動平均は、以下の式で示される。前述の移動平均は、
とみなすことができる。
#define ALPHA 0.5
for( int t = 1 ; t < SIZE ; t++ ) {
y[t] = ALPHA * x[t] + (1.0 - ALPHA) * y[t-1] ;
}
以下のプログラムは、うまく動かない。理由を説明せよ。
#define RVA 4 for( int t = 1 ; t < SIZE ; t++ ) { // 以下はy[t]は全部ゼロになる。 y[t] = 1/RVA * x[t] + (1.0 - 1/RVA) * y[t-1] ; // 以下は、整数型演算だけで、正しく動くだろう。 // y[t] = ( x[t] + (RVA-1) * y[t-1] ) / RVA ; }