s = self.ps*self.skip + self.Ck-1 + self.h-1 X1 = np.zeros((self.n-s, self.w, self.m)) # short term dataset X2 = np.zeros((self.n-s, self.ps*self.Ck, self.m)) # long term dataset Y = np.zeros((self.n-s, self.m)) # output dataset for i in range(s, self.n): t = i-self.h+1 X1[i-s] = self.raw[t-self.w:t].copy() idx = [] for k in range(self.ps): # 执行4次，每次间隔skip长，取Ck=6， total = 4 * 6 = 24 idx = list(range(t-self.Ck-k*self.skip, t-k*self.skip)) + idx idx = np.array(idx, dtype=int) X2[i-s] = self.raw[idx].copy() Y[i-s] = self.raw[i].copy() return X1, X2, Y 把这个单步预测改为多步预测结构

为了将单步预测改为多步预测结构，我们需要对代码进行以下修改：

1. 首先，我们需要定义一个变量 `horizon`，表示我们要预测的时间步数。

2. 然后，在循环中，我们需要将 `Y` 的维度从 `(self.n-s, self.m)` 改为 `(self.n-s-horizon, horizon, self.m)`，表示我们要预测多个时间步。

3. 接下来，我们需要在循环中添加一个内嵌循环，用于预测每个时间步。我们可以使用和单步预测相同的方法进行预测。

4. 最后，我们需要将预测结果添加到 `Y` 中，以便在循环结束后返回。

下面是修改后的代码：

s = self.ps*self.skip + self.Ck-1 + self.h-1                  
X1 = np.zeros((self.n-s, self.w, self.m))  # short term dataset         
X2 = np.zeros((self.n-s, self.ps*self.Ck, self.m)) # long term dataset         
horizon = 3  # 预测的时间步数
Y  = np.zeros((self.n-s-horizon, horizon, self.m))  # output dataset                  

for i in range(s, self.n-horizon):  # 预测第 i+h 时间步
    t = i-self.h+1
    X1[i-s] = self.raw[t-self.w:t].copy()
    idx = []
    for k in range(self.ps): # 执行4次，每次间隔skip长，取Ck=6， total = 4 * 6 = 24
        idx = list(range(t-self.Ck-k*self.skip, t-k*self.skip)) + idx
    idx = np.array(idx, dtype=int)
    X2[i-s] = self.raw[idx].copy()
    
    for j in range(horizon):
        y_pred = self.model.predict([X1[i-s][None, ...], X2[i-s][None, ...]])[0]  # 预测下一个时间步
        Y[i-s, j] = y_pred  # 将预测结果添加到 Y 中
        X1[i-s, :-1] = X1[i-s, 1:]  # 更新 X1
        X1[i-s, -1] = y_pred
        X2[i-s, :-self.Ck] = X2[i-s, self.Ck:]  # 更新 X2
        X2[i-s, -self.Ck:] = X1[i-s, -self.Ck:]  
        
return X1, X2, Y