特異スペクトル変換による変化検知 - プログラム

違いを生んだ、データの生成部分だけを示す。

3パルス目以降の offsetの値だけが前回と異なる

pat_len = 3000

pat = np.empty(3600)

for i in range (5):

    if (i<2):

        offset = 720*i

    else :

        offset = 720*i - 360

    for p in range (720):

        pat[offset+p] = (math.exp(-p/60)+0.01) * (1-math.cos(math.pi * p / 60))