• 售前

  • 售后

热门帖子
入门百科

Matlab实现时间序列猜测分类实例代码

[复制链接]
123456865 显示全部楼层 发表于 2021-8-13 22:01:29 |阅读模式 打印 上一主题 下一主题
目录


  • 一、数据预备
  • 二、时间序列预测分类

    • 1、输入为xt,输出是yt
    • 2、有x值,有y值:NARX

      • (1)选择模型类型
      • (2)选择输出,只有y_t
      • (3)选择70%用来作为训练数据,15%用来作为验证利用,15%用来测试
      • (4)选择delay
      • (5)开始训练
      • (6)得到参数
      • (7)将神经网络导出代码

    • 3、无x,有y值:NAR

  • 三、总结
Matlab从2010b版本以后,神经网络工具箱已经升级为7.0,功能大大加强。在之前的版本做时间预测是比力贫苦操纵的,MathWorks公司对时间序列预测做了具体的办理,跑模型非常简便。
下面通过一个例子演示在Matlab实现时间序列预测。

一、数据预备

极客范儿在夏天吹电扇的体温变革
            时间            风速            温度                                    0            1            37.21405                            0.12457            1.01            37.26016                            0.24915            1.02            37.26324                            0.37373            1.03            37.31242                            0.4983            1.04            37.3155                            0.62258            1.05            37.36468                            0.74745            1.06            37.36776                            0.87203            1.07            37.41694                            0.99661            1.08            37.42002                            …            …            …        
  1. % 原始数据读入到Matlab中
  2. rawData=xlsread('time_series_data.xlsx','sheet1','A2:C52);
  3. % 第一列时间,第二列风速,第三列温度
  4. % yt第三列
  5. y_t=rawData(:,3);
  6. % xt第二列
  7. x_t=rawData(:,2);
复制代码
二、时间序列预测分类

时间序列预测分为三类:

1、输入为xt,输出是yt

即有过去的输入xt,也有过去的输出yt,同时当前的输出不但依赖于过去的输入,也同时依赖于过去的输出
过去时间段温度的变革,预测未来某个时间温度的变革,这种环境就是只有过去的输出
  1. %   x_t - 时间序列输入
  2. %   y_t - 反馈时间序列
  3. X = tonndata(x_t,false,false);
  4. T = tonndata(y_t,false,false);
  5. % 选择训练功能
  6. % 'trainlm'通常是最快
  7. % 'trainbr'耗时较长,但可能更适合解决挑战性的问题
  8. % 'trainscg'使用更少的内存。适用于低内存情况
  9. trainFcn = 'trainlm';  % Levenberg-Marquardt反向传播
  10. % 创建一个非线性自回归网络
  11. feedbackDelays = 1:6;
  12. hiddenLayerSize = 20;
  13. net = narnet(feedbackDelays,hiddenLayerSize,'open',trainFcn);
  14. % 为训练和模拟准备数据
  15. % PREPARETS函数为特定网络准备时间序列数据
  16. % 移动时间的最小量,以声明填充输入状态和层
  17. % 使用PREPARETS允许保留原始的时间序列数据不变,同时轻松定制它的网络与不同
  18. % 具有开环或闭环反馈模式的延迟数
  19. [x,xi,ai,t] = preparets(net,{},{},T);
  20. % 建立训练,验证,测试的数据
  21. net.divideParam.trainRatio = 70/100;
  22. net.divideParam.valRatio = 15/100;
  23. net.divideParam.testRatio = 15/100;
  24. % 训练静态神经网络
  25. [net,tr] = train(net,x,t,xi,ai);
  26. % 测试神经网络
  27. y = net(x,xi,ai);
  28. e = gsubtract(t,y);
  29. performance = perform(net,t,y)
  30. % 查看神经网络
  31. view(net)
  32. % Plots
  33. % Uncomment these lines to enable various plots.
  34. %figure, plotperform(tr)
  35. %figure, plottrainstate(tr)
  36. %figure, ploterrhist(e)
  37. %figure, plotregression(t,y)
  38. %figure, plotresponse(t,y)
  39. %figure, ploterrcorr(e)
  40. %figure, plotinerrcorr(x,e)
  41. % 提前预测网络
  42. % 利用该网络进行多步预测
  43. % CLOSELOOP函数将反馈输入替换为直接输入
  44. % 从外部层连接
  45. nets = removedelay(net);
  46. nets.name = [net.name ' - Predict One Step Ahead'];
  47. view(netc)
  48. [xs,xis,ais,ts] = preparets(nets,X,{},T);
  49. ys = nets(xs,xis,ais);
  50. stepAheadPerformance = perform(nets,ts,ys)
复制代码
2、有x值,有y值:NARX

只有过去的输出
假如给环境加一个风扇,这时间有了风速,过去时间风速在改变,同时也在影响温度的改变
Matlab如今提供时间序列预测工具箱,可以在图形界面上举行调参选择,利用命令
  1. ntstool
复制代码
打开时间序列预测工具箱
雷同股票的模型,只知道早上9:30开市到11:30的股票行情,预测11:30之后的股票行情,不思量任何的输入

(1)选择模型类型



(2)选择输出,只有y_t



(3)选择70%用来作为训练数据,15%用来作为验证利用,15%用来测试



(4)选择delay



(5)开始训练



(6)得到参数




(7)将神经网络导出代码



3、无x,有y值:NAR

没有线性的输入输出,很少遇到这种环境

三、总结

Matlab从2010b版本以后,利用图形界面训练网络调参,天生的代码与手敲的功能无异,Matlab时间序列预测工具箱实用而且好用。
到此这篇关于Matlab实现时间序列预测分类的文章就介绍到这了,更多相干Matlab时间序列预测内容请搜刮脚本之家从前的文章或继续欣赏下面的相干文章盼望各人以后多多支持脚本之家!

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x

帖子地址: 

回复

使用道具 举报

分享
推广
火星云矿 | 预约S19Pro,享500抵1000!
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

草根技术分享(草根吧)是全球知名中文IT技术交流平台,创建于2021年,包含原创博客、精品问答、职业培训、技术社区、资源下载等产品服务,提供原创、优质、完整内容的专业IT技术开发社区。
  • 官方手机版

  • 微信公众号

  • 商务合作