数学建模【GM(1, 1)灰色预测】

一、GM(1, 1)灰色预测简介

乍一看,这个名字好奇怪,其实是有含义的

  • G:Grey(灰色)
  • M:Model(模型)
  • (1, 1):只含有一个变量的一阶微分方程模型

提到灰色,就得先说说白色和黑色

白色系统

  • 典型例子:电阻器件
  • 特点:内部特征是完全已知的,电压、电流和电阻之间的关系(欧姆定律)是已知的
  • 知道电阻大小后,输入电压值,就能算出电流值
  • 白色系统课“计算”出想要的结果

黑色系统

  • 典型例子:一辆车
  • 特点:内部特征是完全未知的,开车并不需要懂发动机设计和工作原理等
  • 非专业人士虽然可以控制汽车,然而内部出故障时并不会修
  • 黑色系统具有“不可知性”

灰色系统

  • 典型例子:GDP就是灰色系统
  • 特点:介于黑色和白色之间,部分已知,部分未知,具有小样本数据的不确定系统
  • 我们有往年的数据和一定的理论基础(白色)
  • 但无法精确计算得出下一年的值(黑色)
  • 灰色无法“计算”,但并不是完全“不可知”,可以进行“预测”

二、适用赛题

数列预测

  • 特点:定时求量,已知xx年到xx年的数据,请预测下一年的数值
  • 常见GDP、人口数量、耕地面积、粮食产量等问题
  • 针对的问题往往短期波动小、可预测,但长期可能变化大、难以准确预测

灾变预测

  • 特点:定量求时,已知xx年到xx年的数据和某灾变的阈值,预测下一次灾变发生的时间
  • 常见洪涝灾害、虫灾等问题
  • 模型中需要把超出阈值的数据(异常数据)对应的时间组成新序列

拓扑预测

  • 特点:对数据波形进行预测,求的是多个模型构成的模型群,等于求解多个灾变预测
  • 与灾变预测类似,不过有较详细的分级,例如虫灾“轻微”“中度”“重度”

注意事项

  • 需要的数据量少,而且数据量太多了没意义,例如用近100年去预测下一年毫无意义
  • 只能短期预测,究竟多短没有严格限制

三、模型流程

四、流程分析

这里以一个例子贯穿流程分析

问题是预测下一年噪声数据

这个题目特点:数据少、看不出明显规律,适合用灰色预测

1.级比检验和平移变换

为了不让我们一顿操作猛如虎,一看预测不可靠,我们可以在开始之前对数据进行检验,确定原始数据GM(1, 1)模型的可行性。

这里

是原始数据的意思

如果在区间外,可尝试平移变换。也就是给每个数据都加上任意常数c后看是否在区间内,求解后再减去c。如果尝试多次平移变换后始终无法在区间内,说明题目不适合灰色预测。

2.累加数据

因为原始数据看不出什么规律,所以我们要进行操作制造规律,一般是累加数据

弱化其随机性,显现其规律性

3.构造方程

生成的新序列,看起来像一个指数曲线。因此可用一个指数曲线的表达式来逼近这个新序列,相应可构建一个一阶常微分方程来求解拟合指数曲线的函数表达式。

要预测下一年数值,就需要知道新序列的表达式,那就要解出微分方程。要解微分方程,就要先知道参数a和u

4.求参数

5.求解方程和预测值

拟合值与预测值

  • 因为第1到7年的噪声数据是已知的,那么通过微分方程的函数解求出的值就是拟合值
  • 拟合值与实际值的偏差,代表了模型的优劣
  • 偏差越小,模型越好。偏差过大说明有问题,所以下一步需要拟合值进行模型检验
6.检验

至此,预测完成。

注:笔者在CSDN上打不出那些变量符号,所以一些难打的符号都是以截图形式呈现

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

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

相关文章

大数据开发项目--音乐排行榜

环境:windows10,centos7.9,hadoop3.2、hbase2.5.3和zookeeper3.8完全分布式; 环境搭建具体操作请参考以下文章: CentOS7 Hadoop3.X完全分布式环境搭建 Hadoop3.x完全分布式环境搭建Zookeeper和Hbase 1. 集成MapReduce…

消息中间件之RocketMQ源码分析(十八)

Broker CommitLog索引机制中的构建过程 1.创建ConsumeQueue和IndexFile。 ConsumeQueue和IndexFile两个索引都是由ReputMessageService类创建的 RequestMessageService类图 ReputMessageService服务启动后的执行过程。 doReput()方法用于创建索引的入口,通常通过…

运用工具Postman快速导出python接口测试脚本

Postman的脚本可以导出多种语言的脚本,方便二次维护开发。 Python的requests库,支持python2和python3,用于发送http/https请求 使用unittest进行接口自动化测试 一、环境准备 1、安装python(使用python2或3都可以)…

【数据结构与算法】常用算法 前缀和

🎉🎉欢迎光临🎉🎉 🏅我是苏泽,一位对技术充满热情的探索者和分享者。🚀🚀 🌟特别推荐给大家我的最新专栏《数据结构与算法:初学者入门指南》📘&am…

离散数学——树思维导图

离散数学——树思维导图 文章目录 前言内容大纲参考 前言 这是当初学习离散数学时整理的笔记大纲,其中包含了自己对于一些知识点的体悟。现将其放在这里作为备份,也希望能够对你有所帮助。 当初记录这些笔记只是为了在复习时更快地找到对应的知识点。…

Cubase学习:Cubase高效技巧,快速复制,快速写音符

1小节后再写 写歌时,建议在第3小节开始写。 快速复制 可以选中后CtrlD,进行快速的复制粘贴。 快速写音符 先是切换到选择状态,按Alt时写音符,松开时回到选择状态。 写一个音时,先确定是哪个区域的,再…

Python爬虫之极验滑动验证码的识别

极验滑动验证码的识别 上节我们了解了可以直接利用 tesserocr 来识别简单的图形验证码。近几年出现了一些新型验证码,其中比较有代表性的就是极验验证码,它需要拖动拼合滑块才可以完成验证,相对图形验证码来说识别难度上升了几个等级。本节将…

MinIO搭建指南:构建高性能的私有云存储服务

MinIO搭建指南:构建高性能的私有云存储服务 摘要:本文将指导您一步步搭建MinIO,一个高性能的开源对象存储服务,用于构建私有云存储解决方案。我们将介绍MinIO的基本概念、安装配置、以及如何与应用程序集成。 一、MinIO简介 Mi…

如何在Linux部署Portainer并结合内网穿透远程管理本地Docker容器

文章目录 前言1. 部署Portainer2. 本地访问Portainer3. Linux 安装cpolar4. 配置Portainer 公网访问地址5. 公网远程访问Portainer6. 固定Portainer公网地址 前言 Portainer 是一个轻量级的容器管理工具,可以通过 Web 界面对 Docker 容器进行管理和监控。它提供了可…

数据分析---Python与sql

目录 Python的pandas,如何实现SQL中的leftjoinPython的pandas,如何实现SQL中的unionPython的pandas,如何实现类似SQL中的where进行限制Python的pandas,如和实现SQL中的group byPython的pandas,如何删除某一列Python的pandas,如何实现SQL中的leftjoin 在Python的pandas库中…

[C++]std::function用法

从函数指针到std::function 函数指针的作用大家应该十分熟悉,它使得我们可以把函数当成参数传递。在C语言中,这种方法几乎非常完美,可以实现基本上所有的传递函数的操作。但是在C中,由于仿函数的出现,C语言中的函数指针…

python装饰器@runWithTry的解析

出发点: 在函数的开始位置有时会看到runWithTry的字样,这样的好处是什么呢? 解析: Python 装饰器,用于包装一个函数,可能为了提供异常处理、日志记录、性能监测等额外功能。在 Python 中,装饰…

如何选择使用哪种剪枝方法?

确定使用哪种剪枝方法通常需要考虑多个因素,包括你的模型类型、任务类型、数据集特性以及训练资源。以下是一些指导原则: 模型类型: 决策树: 对于决策树,可以考虑使用预剪枝或后剪枝。预剪枝通常在构建树的过程中根据节点属性进行判断,而后剪枝则在树已构建完成后根据一些…

用c# 自己封装的Modbus工具类库源码

前言 Modbus通讯协议在工控行业的应用是很多的,并且也是上位机开发的基本技能之一。相关的类库也很多也很好用。以前只负责用,对其并没有深入学习和了解。前段时间有点空就在这块挖了挖。想做到知其然还要知其所以然。所以就有了自己封装的Modbus工具类库…

Mysql 常用数据类型

数值型(整数)的基本使用 如何定义一个无符号的整数 数值型(bit)的使用 数值型(小数)的基本使用 字符串的基本使用 字符串使用细节 日期类型的基本使用

前缀和算法

一、简析前缀和 有一系列元素 A [ a 0 , a 1 , . . . , a n , . . . ] A[a_0,~a_1,~...,~a_n,~...] A[a0​, a1​, ..., an​, ...],前缀和 p r e _ s u m [ n ] A [ 0 ] A [ 1 ] ⋅ ⋅ ⋅ A [ n ] pre\_sum[n]A[0]A[1]A[n] pre_sum[n]A[0]A[1]⋅⋅⋅A[n]。 …

Orange3数据预处理(列选择组件)数据角色及类型描述

在Orange3的文件组件中,datetime、categorical、numeric以及text代表不同种类的数据类型,具体如下: datetime:代表日期和时间类型的数据。通常用于时间序列分析、生存分析和其他需要考虑时间因素的机器学习任务中。例如&#xff0…

图像读取裁剪与人脸识别

图像读取 Image read ⇒ \Rightarrow ⇒ torchvision.datasets from torchvision import datasets dataset datasets.ImageFolder(data_dir, transformtransforms.Resize((512, 512)))Return value illustration dataset[0][0]是PIL.Image objects,这利用IPyth…

小红书关键词爬虫

标题 1 统计要收集的关键词,制作一个文件夹2 爬取每一页的内容3 爬取标题和内容4 如果内容可以被查看,爬取评论内容5 将结果进行汇总,并且每个帖子保存为一个json文件,具体内容6 总结1 统计要收集的关键词,制作一个文件夹 例如,我要收集旅游相关的,就收集: 旅游、旅行…

模型训练基本结构

project_name/ │ ├── data/ │ ├── raw/ # 存放原始数据 │ ├── processed/ # 存放预处理后的数据 │ └── splits/ # 存放数据集划分(训练集、验证集、测试集等) │ ├── noteboo…