灰色预测GM(1,1)

news2024/11/26 3:53:22

目录

 

一、灰色预测模型简介

 二、GM(1,1)灰色模型

 1、GM(1,1)模型预测方法

(1)原始数据(参考列)

(2)累加生成序列(Acumulated Generating Operator,1-AGO)

(3)生成均值数列

(4)建立灰微分方程

(5)利用最小二乘估计【2】中参数a,b

(5.1)

(5.2)构造

(6)求解方程【3】,得到预测方程

2、GM(1,1)模型的预测步骤

(1)数据的检验与处理

(2)建立模型

(3)检验预测值 

(3.1)残差检验

(3.2)级差比偏差检验

(4)预测预报公式​编辑【5】


 

一、灰色预测模型简介

灰色预测的主要特点就是模型使用的不是原始数据列,而是生成的数据序列。其核心是灰色模型(Grey Model,GM),即对原始数据作累加生成(累加或差分,都能弱化噪声)能得到近似指数规律再进行建模的方法。

优点:

  1. 所需数据少,一般只需要4个以上数据,能解决历史数据少,序列完整性和可靠性低的问题
  2. 利用微分方程来充分挖掘系统的本质,精度高
  3. 能将无规律的原始数据进行生成得到规律性较强的生成序列
  4. 运算简便,易于检验,不考虑分布规律,不考虑变化趋势

缺点:

只适合短期预测,只适合指数增长的预测。

 二、GM(1,1)灰色模型

GM(1,1),表示模型是一阶微分方程,且只含一个变量的灰色模型。

 1、GM(1,1)模型预测方法

(1)原始数据(参考列)

395d613b27424b67bcecd4ac29185480.png

(2)累加生成序列(Acumulated Generating Operator,1-AGO)

b224fb93716f454ab9f5215f20dfd749.png

其中,8275db3b63c54755998c49a8d8f234b9.png显然,数列x(1)比数列x(0)光滑,弱化噪声。

(3)生成均值数列

0fc3d37a11674302a3b0314822781d13.png(进一步弱化噪声)

其中ad4b83d6a0704f9aa4f3fd4770d435c1.png

(4)建立灰微分方程

d2b1c854619e4f17a6d6216c07393957.png【1】

此方程白化为微分方程06c9234594ff493b8a60ebdf350efac7.png【2】

【1】的来源:x(1)是x(0)的累加,反过来,x(0)就是x(1)的差分,连续化就是导数,【1】中用z(1)而不直接用x(1),就是使数据进一步更光滑化,避免奇异值。可以利用曲线拟合来求a和b得值(就是利用最小二乘法)。用这个模型之前,最好绘制x(0),z(1)的散点图,看看大概情况。

(5)利用最小二乘估计【2】中参数a,b

5269eaf73c024b01856a9726388cb8e3.png【2】

c95a1424599c4352bc3b495798d3ed5a.png的观测值序列为x(0)(k),k=2,3,…,n;

【2】中x(1)(t)的观测值为z(1)(k),k=2,3,…,n;

(5.1)

把x(0),z(1)作为【2】的因变量y和自变量x,作回归拟合,得到a,b的估计值;

(5.2)构造

53ba49749d6f499791f4379a6b9b7b19.pngb0ce4926d4374618ac8c087f86431b74.pngfea755f5105a498faa87af4bc958cdb1.png

根据最小二乘法,使得326f4526577a4b37984655c0936e4431.png达到最小的u的估计值为6695bb2447c84dfb98c4e612e5f5fd51.png求出a,b的估计值为4ae9885f60464e30a8e866dab6a0b248.png方程【2】变为bd23ef9cf092483fb11f345e9e774c15.png

(6)求解方程【3】,得到预测方程

为了求【3】的特解,引入初值de8036a508764bf7b39542cae5f5f8a3.png确定特解。

【3】为一阶线性常系数微分方程,代入通解公式,代入定解条件,得2388b02b68054856b1822ad628d05c2b.png

写成差分或者预测公式,为97cd65dac636439c9338e14604f99a0c.png【4】k=0,1,2,…,n-1,…

2、GM(1,1)模型的预测步骤

(1)数据的检验与处理

设原始数列为31298b1c75594fa0a3c8a2b68d6f7f4e.png计算数列的级比

85d2e35b4e8542de9e703f5ac8539488.png

b8d5f61fd7284ca5964bf152f455f7db.png(一个包含1的区间)

则认为原始数列适合作GM(1,1)灰色预测,否则需要对x(0)作平移:6a424b0901e24e5cb453ac6de5ac2b18.png使得数列y(0)满足级比检验。

(2)建立模型

5f6eac8e6cdb4eceaf4887915f519c84.png【4】

 当初始数列或平移后的数列满足级比检验后,可用结论【4】进行预测,1357431af445421997f69643a717f517.png

(3)检验预测值 

(3.1)残差检验

计算残差92f9dc4d117d4c7fa5f067defaace224.png

这里,f49e3622734d47829ebcf8a3f8d4b1a3.png

70cf86e9ad5449c3b671d5789d9e5c8d.png则认为预测达到一般要求;

4f4596dbaccb42f09deb298a35260ce0.png则认为预测达到较高要求;

(3.2)级差比偏差检验

首先用原始数据x(0)(k),x(0)(k-1)计算级比λ(k),再用发展系数a求出相应的级比偏差95c28d289f5e4dc7943415ba595f2cbc.png

若ρ(k)<0.2,认为预测达到一般要求;若ρ(k)<0.1,则认为预测达到较高要求。这里的a就是最小二乘估计的第一个参数,称为发展系数

(4)预测预报公式8e8866264dec4722bae09ce55c5a1d62.png【5】

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1042818.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

一台电脑安装多个不同版本Python

1、前提 当前Windows电脑下已经安装了一个python3.11&#xff0c;现在需要安装一个python3.9。下载地址&#xff1a;Python Releases for Windows | Python.org 2、步骤 找到对应的版本&#xff0c;并下载安装包。下载后&#xff0c;打开安装包。按图片勾选&#xff0c;点击n…

【网络安全自学宝典】从零开始自学网络安全,按照这个路线就可以了

前景 很多零基础朋友开始将网络安全作为发展的大方向&#xff0c;的确&#xff0c;现如今网络安全已经成为了一个新的就业风口&#xff0c;不仅大学里开设相关学科&#xff0c;连市场上也开始大量招人。那么网络安全到底前景如何&#xff1f;大致从市场规模、政策扶持、就业方…

共赴华为全联接大会2023,新点软件携手加速数字政府建设

数字技术使能&#xff0c;加速行业智能化。9月20—22日&#xff0c;华为全联接大会2023在上海拉开序幕&#xff0c;来自全球ICT产业的业界思想领袖、商业精英、技术专家、先锋企业、生态合作伙伴、应用服务商以及开发者共襄盛会&#xff0c;共同探讨行业智能化的未来路径。 作…

【c语言】详解动态内存管理

目录 关于动态内存分配malloc和calloc函数介绍动态内存回收----freerealloc函数介绍常见的动态内存错误 关于动态内存分配 回想一下我们之前学过的内存开辟方式&#xff1a; int val 20;//在栈空间上开辟四个字节 char arr[10] {0};//在栈空间上开辟10个字节的连续空间在学…

Windows系统ping命令的c++实现

// ping.cpp : 定义控制台应用程序的入口点。 //#include <winsock2.h> #include <stdio.h> #include <stdlib.h> #define ICMP_ECHO 8 #define ICMP_ECHOREPLY 0 /* The IP header */ typedef struct iphdr { unsigned char h_len:4; // length of the…

什么是内容运营?

关于内容运营&#xff0c;在不同种类的公司&#xff0c;侧重点也不一样。 电商平台的内容运营岗更偏内容营销&#xff1b;产品功能比较简单的公司&#xff0c;内容运营和新媒体运营的岗位职责差不多&#xff1b;而内容平台的内容运营更多的是做内容的管理和资源整合。

使用ElementUI结合Vue完善主页的导航菜单和书籍管理的后台数据分页查询

目录 一、动态树 ( 1 ) 数据表 ( 2 ) 后端 ( 2 ) 前端 二、书籍管理 数据表 后端 前端 ElementUI的背景 是一套基于Vue.js的桌面端组件库&#xff0c;由饿了么前端团队开发维护。它提供了丰富的UI组件和交互效果&#xff0c;可以帮助开发者快速构建出美观、易用的We…

第77步 时间序列建模实战:多因素预测 vol-2(以ARIMA为例)

基于WIN10的64位系统演示 一、写在前面 上一期&#xff0c;我们构建了多变量的ARIMA时间序列预测模型&#xff0c;其实人家有名字的&#xff0c;叫做ARIMAX模型&#xff08;X就代表解释变量&#xff09;。 这一期&#xff0c;我们介绍其他机器学习回归模型如何建立多变量的时…

『C语言进阶』qsort函数及模拟实现

&#x1f525;博客主页&#xff1a; 小羊失眠啦 &#x1f516;系列专栏&#xff1a; C语言 &#x1f325;️每日语录&#xff1a;没有退路&#xff0c;只能让自己变得强大 ❤️感谢大家点赞&#x1f44d;收藏⭐评论✍️ 前言 在上篇指针进阶中&#xff0c;我们对函数指针、函数…

深入理解 Swift 新并发模型中 Actor 的重入(Reentrancy)问题

问题现象 我们知道,Swift 5.5 引入的新并发模型极大简化了并行逻辑代码的开发,更重要的是:使用新并发模型中的 Actor 原语可以大大降低并发数据竞争的可能性。 不过,即便 Actor 有如此神奇之功效,它也不是“万能药”,仍不能防止误用带来的问题。比如:Actor 重入(Reen…

132.【MySQL_进阶篇】

MySQL_进阶 (一)、存储引擎1.MySQL体系结构(1).连接层(2).服务层(3).引擎层(4).存储层 2.存储引擎简介(1).查看某张表的数据引擎(2).展示此版本支持的所有存储引擎(3).创建表my_myisam,并指定MyIASM存储引擎(4).存储引擎示列 3.存储引擎 _ Innodb(1).Innodb 介绍(2).Innodb 特点…

通过 HelpLook ChatBot AI自动问答机器人降低客户服务成本

在当今竞争激烈的商业环境中&#xff0c;提供卓越的客户服务对于维持忠诚的客户群和推动业务增长至关重要。客户服务涵盖了公司与其客户之间的所有互动&#xff0c;包括解答问题、解决问题和提供支持。它在塑造客户对品牌的看法方面起着关键作用&#xff0c;并且可以显著影响他…

NSSSCTF做题(2)

1.[BJDCTF 2020]easy_md5 打开页面发现没什么东西&#xff0c;只有一个提交表单&#xff0c;然后url会显示你提交的信息 源代码里也看不到什么 &#xff0c;用dirsearch扫一下&#xff0c;这些是扫出来的数据 抓包看到了提示 hint: select * from admin where passwordmd5($pas…

【文件操作——详细讲解】

1. 为什么使用文件&#xff1f;&#x1f9d0; 如果没有⽂件&#xff0c;我们写的程序的数据是存储在电脑的内存中&#xff0c;如果程序退出&#xff0c;内存回收&#xff0c;数据就丢失了&#xff0c;等再次运⾏程序&#xff0c;是看不到上次程序的数据的&#xff0c;如果要将数…

cmd下scala退出命令

cmd下scala退出命令 &#xff1a;quit

checksec使用

checksec Relro&#xff1a;Full Relro&#xff08;重定位表只读&#xff09; Relocation Read Only&#xff0c; 重定位表只读。重定位表即.got 和 .plt个表。 Stack&#xff1a;No Canary found&#xff08;能栈溢出&#xff09; 栈保护。栈溢出保护是一种缓冲区溢出攻击缓解…

Docker 部署前端项目(非自动化)

部署前端项目 1. nginx配置文件1.1 nginxConf 2. 创建容器2.1 添加项目2.2 下载项目依赖2.3 打包前端项目2.4 创建容器2.5 查看容器 3. 视频演示4. 注意 1. nginx配置文件 1.1 nginxConf 首先你需要有nginx配置文件&#xff0c;你可以执行以下命令获取配置文件 # 安装镜像-生成…

JS进阶-闭包

概念&#xff1a;一个函数对周围状态的引用捆绑在一起&#xff0c;内层函数中访问到其外层函数的作用域 简单理解&#xff1a;闭包&#xff08;Closure&#xff09;内层函数外层函数的变量 function outer() {const a 1function f() {console.log(a)}f()} outer() 闭包作用…

最新商道融绿ESG数据(2015-2023)

数据简介&#xff1a;环境、社会和公司治理&#xff08;environmental, social, and corporate governance, 简称ESG&#xff09;方面的非财务信息已经普遍被用于评估企业在可持续发展方面的绩效。投资者通过考察上市公司的ESG绩效&#xff0c;便于识别企业面临的风险&#xff…

L1-033 出生年 c++解法

一、题目再现 以上是新浪微博中一奇葩贴&#xff1a;“我出生于1988年&#xff0c;直到25岁才遇到4个数字都不相同的年份。”也就是说&#xff0c;直到2013年才达到“4个数字都不相同”的要求。本题请你根据要求&#xff0c;自动填充“我出生于y年&#xff0c;直到x岁才遇到n个…