0

Henon写像のESN学習を行なっているのですが上記のエラーが起こってしまいます。
エラーの箇所は予測性能評価の部分で起こっていると思うのですが原因がわかりません。
よろしければ回答の程をお願いします。

import numpy as np
import matplotlib.pyplot as plt
from model import ESN, Tikhonov


np.random.seed(seed=0)

# Henon写像によるデータ生成
class Henon:
    # パラメータの指定
    def __init__(self, a, b):
        self.a = a
        self.b = b

    def f1(self, t, x, y):
        return 1 - self.a*x**2 + y

    def f2(self, t, x, y):
        return self.b*x


    def Henon(self, t, X):
        '''
        :param t: 時間
        :param X: 2次元ベクトル
        :return: 2次元ベクトル
        '''
        next_X = [self.f1(t, X[0], X[1]), 
                  self.f2(t, X[0], X[1])]
        return np.array(next_X)


    def Date_Henon(self, x0, T, dt):
        '''
        :param x0: 初期値
        :param T: 時間
        :param dt: ステップ幅
        :return: Henon写像
        '''
        X = x0
        t = 0
        data = []

        while t < T:
            next_X =  self.Henon(t,X)
            data.append(next_X)
            X = next_X
            t = t + dt

        return np.array(data)

if __name__ == '__main__':

    # Lorenz方程式から時系列データ生成
    T_train = 1000  # 学習データの長さ
    T_test = 100# テストデータの長さ
    dt = 1.0
    x0 = np.array([1, 1])  # 初期値
    
    dynamics = Henon(1.4, 0.3)
    data = dynamics.Date_Henon(x0, T_train + T_test, dt)

    # 訓練・検証用情報
    train_U = data[:T_train]
    train_D = data[1:T_train+1]

    test_U = data[T_train:T_train + T_test]
    test_D = data[1 + T_train:T_train + T_test + 1]
    
    N_x = 1000  # リザバーのノード数
    model = ESN(train_U.shape[1], train_D.shape[1], N_x, \
                density=0.1, input_scale=0.1, rho=0.95)

    # 学習(リッジ回帰)
    train_Y = model.train(train_U, train_D, 
                          Tikhonov(N_x, train_D.shape[1], 1e-4))

    # モデル出力(自律系のフリーラン)
    test_Y = model.predict(test_U)
    T0 = test_D[:, 0]
    T1 = test_Y[:, 0]
    plt.plot(T0,label= 'テストデータ')
    plt.plot(T1)


    # 予測性能評価(valid time)
    eps = 1  # 許容誤差
    for n in range(int(T_test/dt)):
        dif = np.sqrt(((test_D[n,:] - test_Y[n,:]) ** 2).mean())
        if dif > eps:
            valid_time = float(n)*dt
            break

    print('valid time = ', valid_time)
1
  • 3
    ja.stackoverflow.com/questions/53823/… こちらにおそらく同様のエラーの質問内容があります。おそらく配列の要素数として99個も無いよ、という意味なのでfor文部分などで参照しようとしている要素数が元々の配列の要素数を超えていないか確認できますか?(for n in range(int(T_test/dt)):のところのrange内部とかか、配列の要素数をカウントしてみてください)
    – MSAII
    Commented 2022年10月3日 3:51

1 件の回答 1

1

既に解決済みかもしれませんが,念のため回答いたします。
コメント欄のアドバイスに従い,68行目

test_D = data[1 + T_train:T_train + T_test + 1]

の直後に

print(test_D.shape, data.shape)

を挿入して調べてみると

(99, 2) (1100, 2)

が得られるので, test_D の( axis 0 の)サイズが 99 であり,これがエラーの原因と考えられます。
なお,68行目の記述を実際の数値( T_train=1000, T_test=100 )で置き換えてみると

test_D = data[1001:1101]

なので, test_D のサイズを 100 にするには,data のサイズが 1101 必要です。

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。