相関.

Presentation on theme: "相関."— Presentation transcript:

1 相関

2 似ている

3 フーリエ係数は似ている度合 フーリエ係数が大きい＝似ている

4 一般的な 二つの関数の似ている度合を で判定する

5 ところが、、、、 であった。 この二つは、“まったく似ていない”　という判定になってしまう

6 と は同じもの

7 片方の関数をずらして、 判定する必要がある

8 相互相関係数 相互相関係数は、定数ではなく、ずらす量τの関数となる。 実際は、無限大まで考えるのは難しいので
として、Tをなるべく大きくとる。

9 Excelで計算する場合 fiはf(t)のi番目のデータ giはg(t＋τ）のi番目のデータ Rを関数ｆとgの相互相関係数という

10 Excelでの計算 と の相互相関係数

11 相互相関係数は、平均で出せる

12 τをいろいろ変えて、相互相関係数を見てみる

13 相互相関係数を 時間の関数としてみてみる

14 実際の信号の解析 信号が式で表せない、ノイズ入り信号の場合は、難しい これの相互相関係数を求めたい

15 表をずらして計算 これは かなり大変

16 自己相関係数 自分のコピーを作り、それをずらして相関係数を調べる

17 マクロの利用

18 VBA Visual　Basic　というプログラミング言語でプログラムを書く

19 VBA Dim f(1000) Dim g(1000) Dim t(1000) 'Nはでデータ数
Sub sougosoukan() Dim f(1000) Dim g(1000) Dim t(1000) 'Nはでデータ数 N = Cells(2, 1).End(xlDown).Row - 1 Cells(1, 5) = "データ数" Cells(2, 5) = N 'データ読み込み For i = 1 To N t(i) = Cells(i + 1, 1) f(i) = Cells(i + 1, 2) g(i) = Cells(i + 1, 3) Next i Cells(1, 6) = "τ" Cells(1, 7) = "相互相関係数" For j = 0 To N - 1 Cells(2 + j, 6) = t(j + 1) - t(1) s = 0 For k = 1 To N - j y = f(k) * g(k + j) s = s + y Next k soukan = s / (N - j) Cells(2 + j, 7) = soukan Next j End Sub Sub jikosoukan() Dim f(1000) Dim t(1000) 'Nはでデータ数 N = Cells(2, 1).End(xlDown).Row - 1 Cells(1, 4) = "データ数" Cells(2, 4) = N 'データ読み込み For i = 1 To N t(i) = Cells(i + 1, 1) f(i) = Cells(i + 1, 2) Next i Cells(1, 5) = "τ" Cells(1, 6) = "自己相関係数" For j = 0 To N - 1 Cells(2 + j, 5) = t(j + 1) - t(1) s = 0 For k = 1 To N - j y = f(k) * f(k + j) s = s + y Next k soukan = s / (N - j) Cells(2 + j, 6) = soukan Next j End Sub

20 使い方（自己相関係数） =SIN(2*PI()/3*A3)+(4*RAND()-2) 時間は少し長めに(20秒くらい）とる

21 マクロの実行

22 結果 3秒、6秒、9秒 ずらすとよく重なる →周期3秒

23 使い方（相互相関係数）

24 マクロの実行

25 結果 2.25秒ずらすと よく重なる