当前位置：首页 > news >正文

基于GRU、CNN-GRU、SSA-CNN-GRU和ICEEMDAN-SSA-CNN-GRU模...

news 2026/7/5 7:23:45

GRU，CNN-GRU，SSA-CNN-GRU，ICEEMDAN-SSA-CNN-GRU做时间序列预测，多输入单输出预测。精度以ICEEMDAN-SSA-CNN-GRU。评价指标： RMSE = 0.08024 MSE = 0.0064385 MAE = 0.071505 MAPE = 0.05383

时间序列预测这活儿就像在炒菜——火候不够夹生，火大了容易糊。最近折腾了几个组合模型，从基础GRU到魔改版ICEEMDAN-SSA-CNN-GRU，发现这玩意儿真是越堆料效果越香。今天咱们就掰开揉碎了聊聊这些模型的实操细节，顺便看看怎么用代码把它们串起来。

先说说基础款的GRU网络。这兄弟比LSTM轻量，门控机制简化后对短期记忆处理很友好。下面这段PyTorch实现看着清爽：

class GRUModel(nn.Module): def __init__(self, input_dim, hidden_dim, output_dim): super().__init__() self.gru = nn.GRU(input_dim, hidden_dim, batch_first=True) self.fc = nn.Linear(hidden_dim, output_dim) def forward(self, x): out, _ = self.gru(x) out = self.fc(out[:, -1, :]) # 取最后一个时间步 return out

重点在forward里那个切片操作out[:, -1, :]，相当于只抓序列末尾的特征。但实测发现当输入序列较长时，这种简单粗暴的取法会漏掉局部特征，这时候就该CNN上场了。

CNN-GRU的杂交品种在特征提取层动了手脚：

class CNN_GRU(nn.Module): def __init__(self, cnn_dim, gru_dim): super().__init__() self.conv1 = nn.Conv1d(in_channels=cnn_dim, out_channels=64, kernel_size=3) self.gru = nn.GRU(64, gru_dim) self.fc = nn.Linear(gru_dim, 1) def forward(self, x): x = x.permute(0, 2, 1) # 调换维度适应卷积 x = F.relu(self.conv1(x)) x = x.permute(0, 2, 1) _, hn = self.gru(x) return self.fc(hn.squeeze(0))

这里有个维度变换的细节值得注意：原始输入是(batchsize, seqlen, features)，卷积层需要把features放在中间维度。卷积核沿着时间轴滑动，抓局部时间模式比单纯用GRU更带劲。

但这样还不够——参数初始化要是没弄好，模型可能困在局部最优里出不来。这时候SSA（麻雀搜索算法）就来救场了。这个元启发式算法调参比网格搜索高效得多，特别是处理超参数如卷积核数量、GRU隐藏层维度时：

def ssa_optimize(): param_ranges = { 'cnn_filters': (32, 128), 'gru_units': (64, 256) } # 麻雀种群初始化 population = initialize_sparrows() for epoch in range(100): # 评估当前参数组合的验证集损失 fitness = evaluate(population) # 更新麻雀位置（参数组合） population = update_positions(population, fitness)

这个优化过程就像让一群麻雀在参数空间里找食，发现好区域就呼朋引伴。实测下来比随机搜索快3倍左右，关键是不容易早熟。

最后的大杀器是ICEEMDAN分解。传统EMD容易产生模态混叠，改进版的自适应噪声完备分解把原始序列拆得更干净：

from PyEMD import ICEEMDAN def decompose(signal): iceemdan = ICEEMDAN() imfs = iceemdan(signal) return imfs # 对每个IMF分别训练模型 imfs = decompose(raw_sequence) models = [CNN_GRU() for _ in range(len(imfs))] # 各分量预测结果叠加 final_output = sum([model.predict(imf) for imf, model in zip(imfs, models)])

这个分治策略妙在把非平稳信号拆成相对平稳的子信号，每个CNN-GRU模块只需专注自己的频段。实测中MAPE从0.062直降到0.053，效果拔群。

看最终指标，ICEEMDAN-SSA-CNN-GRU在多个维度碾压基础模型：