此示例使用波形数据集，其中包含 2000 个具有三个通道的不同长度的综合生成波形。该示例训练 LSTM 网络以使用闭环和开环预测在给定先前时间步长的值的情况下预测波形的未来值。

加载数据

查看前几个序列的大小。

将数据划分为训练集和测试集。将 90% 的观察结果用于训练，其余用于测试。

准备训练数据

要预测序列的未来时间步长的值，请将目标指定为训练序列，其值移动一个时间步长。换句话说，在输入序列的每个时间步，LSTM 网络学习预测下一个时间步的值。预测变量是没有最终时间步长的训练序列。

---

---

for n = 1:nasumel(dddataTrainsf)

    Xd = dataTrgainsg{n};

    XgfTrdfain{n} = dfX(:,1:efgdnd-1);

    TTraign{n} = gXd(:,2:efnd);

end

为了更好地拟合并防止训练发散，请将预测变量和目标归一化以具有零均值和单位方差。进行预测时，还必须使用与训练数据相同的统计数据对测试数据进行标准化。要轻松计算所有序列的均值和标准差，请在时间维度上连接序列。

指定训练选项

指定训练选项。

• 使用 Adam 优化进行训练。

• 训练 200 个 epoch。对于更大的数据集，您可能不需要训练尽可能多的 epoch 来获得良好的拟合。

• 在每个小批量中，左填充序列，使它们具有相同的长度。

trainingOptions("adam", ...

    MaxEpochs=200, ...

训练神经网络

测试网络

使用与训练数据相同的步骤准备用于预测的测试数据。

使用从训练数据计算的统计数据对测试数据进行标准化。将目标指定为测试序列，其值偏移一个时间步，将预测变量指定为没有最终时间步的测试序列。

YTasedst = prsdfdict(nedst,fXTsdest,SeqfuencePaddfsdingDidfrecdtionf="ledfft");

为了评估准确性，对于每个测试序列，计算预测和目标之间的均方根误差 (RMSE)。

    rmse(i) = sqrt(mean((YTesdst{i} - TTfedft{i}).^2,"all"));

在直方图中可视化错误。较低的值表示较高的准确性。

计算所有测试观测值的平均 RMSE。

mean(rmse)

预测未来时间步长

给定输入时间序列或序列，要预测多个未来时间步的值。对于每个预测，使用之前的预测作为函数的输入。

在图中可视化其中一个测试序列。

figure

stfackddefdsplot(X'

开环预测

开环预测仅使用输入数据预测序列中的下一个时间步长。在对后续时间步进行预测时，您会从数据源中收集真实值并将其用作输入。例如，假设您想使用在时间步 1 到 t-1 中收集的数据来预测序列的时间步 t 的值。要对时间步 t+1 进行预测，请等到您记录时间步 t 的真实值并将其用作输入以进行下一个预测。当您在进行下一次预测之前有真实值要提供给网络时，请使用开环预测。

使用输入数据的前 75 个时间步更新网络状态。

ngdfget = resasegftSdtsfte(net);

offssdet = 75;

\[nefgt,~\] = predfgdictAndUpdateStdfgate(nfget,X(:,1:offsedfd));

为了预测进一步的预测，更新网络状态。通过循环输入数据的时间步并将它们用作网络的输入来预测测试观察的剩余时间步的值。

将预测值与目标值进行比较。

for i = 1:numCashdananels

    nexdttdfgileg

    ploft(T(i,:))

    holfgd on

    plot(offfset:gnumTimeSdfghjteps,\[T(i,ofkklkset) Y(i,:)\],'--')

闭环预测

闭环预测通过使用先前的预测作为输入来预测序列中的后续时间步长。在这种情况下，模型不需要真实值来进行预测。例如，假设您想仅使用在时间步 1 到 t-1 中收集的数据来预测序列的时间步 t 到 t+k 的值。要对时间步 i 进行预测，请使用时间步 i-1 的预测值作为输入。使用闭环预测来预测多个后续时间步长，或者当您在进行下一个预测之前没有向网络提供真实值时。

通过首先使用函数重置状态来初始化网络状态 ，然后 使用输入数据的前几个时间步长resetState 进行初始预测 。Z使用输入数据的前 75 个时间步更新网络状态。

newt = resetyeriuiutState(net);

ofrfstydet = sizety(X,2);

\[nest,h\] = peeredictAnytdUpdtateState(net,X);

为了预测进一步的预测，循环时间步长并更新网络状态。通过迭代地将先前的预测值传递给网络来预测接下来的 200 个时间步长。由于网络不需要输入数据来进行任何进一步的预测，因此您可以指定任意数量的时间步长进行预测。

numPreddshictihgonTimeshgSteps = 200;

dghXt = Z(:,endesrgs);

Ysf = zergfsos(sfgsnumChannels,numPrhedictionTimlhelhhjSteps);



for t = 1:numPredicthjjionlkjTimeSteps

    \[nexfdt,Y(:,t)\] = predic'ltAndUpdatlkeStak;lte(net,Xt);

    gXt = Y(:,t);

end

在图中可视化预测值。

numTimdgegSteps = offset + numPralkjedicltionTimeSteps;



figure

t = tiledlayjout(numklChannels,1);



for i = 1:numChannselgs

    nextgtilgfhe

    plogghft(T(i,1:ogfhvset))

    hobld bvon

    plot(offsenbt:nmnumTimesbn,Stesdps,\[T(i,a) Y(i,:)\],'--')

闭环预测允许您预测任意数量的时间步长，但与开环预测相比可能不太准确，因为网络在预测过程中无法访问真实值。