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

多维时序 | Matlab实现基于VMD-DBO-BiLSTM、VMD-BiLSTM、BiLSTM的多变量时间序列预测

多维时序 | Matlab实现基于VMD-DBO-BiLSTM、VMD-BiLSTM、BiLSTM的多变量时间序列预测

目录

    • 多维时序 | Matlab实现基于VMD-DBO-BiLSTM、VMD-BiLSTM、BiLSTM的多变量时间序列预测
      • 预测效果
      • 基本介绍
      • 程序设计
      • 参考资料

预测效果

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

基本介绍

Matlab实现基于VMD-DBO-BiLSTM、VMD-BiLSTM、BiLSTM的多变量时间序列预测(完整程序和数据)
1.先运行vmdtest,进行vmd分解;
2.再运行VMD-DBO-BiLSTM,三个模型对比;
3.运行环境Matlab2018及以上。

  • VMD-DBO-BiLSTM:变分模态分解结合蜣螂算法优化双向长短期记忆神经网络;
  • VMD-BiLSTM:变分模态分解结合双向长短期记忆神经网络;
  • BiLSTM:双向长短期记忆神经网络。

程序设计

  • 完整程序和数据获取方式:私信博主回复Matlab实现基于VMD-DBO-BiLSTM、VMD-BiLSTM、BiLSTM的多变量时间序列预测
clc;clear;close all;format compact
tic
clc
clear all
fs=1;%采样频率,即时间序列两个数据之间的时间间隔,这里间隔1h采样
Ts=1/fs;%采样周期X = xlsread('北半球光伏数据.xlsx','C2:E296');save origin_data XL=length(X);%采样点数,即有多少个数据
t=(0:L-1)*Ts;%时间序列
STA=0; %采样起始位置,这里第0h开始采样%--------- some sample parameters forVMD:对于VMD样品参数进行设置---------------
alpha = 2500;       % moderate bandwidth constraint:适度的带宽约束/惩罚因子
tau = 0;          % noise-tolerance (no strict fidelity enforcement):噪声容限(没有严格的保真度执行)
K = 5;              % modes:分解的模态数
DC = 0;             % no DC part imposed:无直流部分
init = 1;           % initialize omegas uniformly  :omegas的均匀初始化
tol = 1e-7         
%--------------- Run actual VMD code:数据进行vmd分解---------------------------
[u, u_hat, omega] = VMD(X(:,end), alpha, tau, K, DC, init, tol);
%  参数设置
options = trainingOptions('adam', ...                 % 优化算法Adam'MaxEpochs', 70, ...                              % 最大训练次数'GradientThreshold', 1, ...                       % 梯度阈值'InitialLearnRate', 0.01, ...                     % 初始学习率'LearnRateSchedule', 'piecewise', ...             % 学习率调整'LearnRateDropPeriod', 60, ...                    % 训练850次后开始调整学习率'LearnRateDropFactor',0.2, ...                    % 学习率调整因子'L2Regularization', 0.01, ...                     % 正则化参数'ExecutionEnvironment', 'cpu',...                 % 训练环境'Verbose', 0, ...                                 % 关闭优化过程'Plots', 'training-progress');                    % 画出曲线%  训练
net = trainNetwork(vp_train, vt_train, layers, options);
%analyzeNetwork(net);% 查看网络结构

参考资料

[1] https://blog.csdn.net/kjm13182345320/category_11799242.html?spm=1001.2014.3001.5482
[2] https://blog.csdn.net/kjm13182345320/article/details/124571691

http://www.lryc.cn/news/305226.html

相关文章:

  • 实战营第四节笔记
  • 【网络】HTTPS协议原理
  • C#常识篇(一)
  • Linux使用C语言获取进程信息
  • tomcat通过JAVA_OPTS注入自定义变量 —— 筑梦之路
  • Linux——简单的Shell程序
  • 计算机网络-网络互联与互联网(一)
  • Android 解决后台服务麦克风无法录音问题
  • 【计网】TCP的三次握手四次挥手
  • android pdf框架-4,分析vudroid源码2
  • [算法沉淀记录] 排序算法 —— 希尔排序
  • conda 进入虚拟环境命令报错
  • 域名 SSL 证书信息解析 API 数据接口
  • 学习JAVA的第二天(基础)
  • 代码随想录算法训练营第二十七天|93.复原IP地址、78.子集、90.子集II
  • 【蓝桥备赛】字串简写
  • nios ii开发随笔
  • SpringBoot项目嵌入RabbitMQ
  • 提升网络质量:UDPspeeder 实现网络优化与提速
  • 为什么前端开发变得越来越复杂了?这可能是我们的错
  • VR系统的开发流程
  • 前端输入框校验限制不能输入中文
  • 强大的文本绘图——PlantUML
  • ES相关问题
  • 基于Linux直接安装的Nginx版本升级方法
  • 探索设计模式的魅力:状态模式揭秘-如何优雅地处理复杂状态转换
  • 力扣hot100题解(python版10-12题)
  • 【算法】复杂度分析
  • 车载电子测试学习内容
  • STM32F103x 的时钟源