matlab 指定区域随机游走_了解随机游走模型和移动平均过程(Python)
在本文中,我将介绍两个能够对时间序列进行建模的模型:随机游走和移动平均过程。
随机游走模型
随机游走模型由以下公式表示:
换句话说,当前时刻t的位置是前一时刻(t-1)的位置与噪声(用z表示)之和。这里我们假设噪声是正态分布的(均值为0,方差为1)。
我们从0开始随机游走,也就是说,任何时间点都是该时间之前所有噪声的和。数学上表示为:
让我们在Python中模拟随机游走。
首先,我们导入所需的Python库:
from statsmodels.graphics.tsaplots import plot_acffrom statsmodels.tsa.arima_process import ArmaProcessfrom statsmodels.tsa.stattools import acfimport matplotlib.pyplot as pltimport numpy as np%matplotlib inline
然后,我们生成一个包含1000个数据点的数据集。起点是0,我们将随机噪声添加到上一个点以生成下一个点:
steps = np.random.standard_normal(1000)steps[0]=0random_walk = np.cumsum(steps)
绘制数据集的Python代码如下:
plt.figure(figsize=[10, 7.5]); # Set dimensions for figureplt.plot(random_walk)plt.title("Simulated Random Walk")plt.show()
模拟随机游走
你的随机行走可能与上面的图不同,因为噪声是随机的。
现在,让我们看看我们的随机游走的自相关图(或相关图):
random_walk_acf_coef = acf(random_walk)plot_acf(random_walk, lags=20);
随机游走的相关图
不管你的随机游走看起来像什么,你都应该得到一个非常相似的相关图。
现在,一切都指向数据集中的趋势。我们可以改变这种趋势吗?答案是肯定的。
让我们在Python中进行验证。
random_walk_diff = np.diff(random_walk, n=1)
然后我们绘制结果:
plt.figure(figsize=[10, 7.5]); # Set dimensions for figureplt.plot(random_walk_diff)plt.title('Noise')plt.show()
如您所见,上面的图没有趋势,也没有季节性,是一个完全随机的过程。
查看相关图的python代码如下:
plot_acf(random_walk_diff,lags = 20);
我们看到这是一个纯随机过程的相关图,其中自相关系数在滞后1处下降。
移动平均过程
我们先来直观了解一下什么是移动平均过程。
假设你把一块石头扔进一个池塘里,你要追踪水面上一滴水的位置。当石头撞击水面时,会形成波纹,所以我们要跟踪的水滴会上下移动。让我们假设波纹只持续两秒钟,在这两秒钟之后,水面就会完全变平。
我们的水滴位置可以表示为:
上面的方程表示,X在t时刻的位置取决于t时刻的噪声,加上t-1时刻的噪声(有一定的权重THETA),加上t-2时刻的噪声(有一定的权值)。
这被称为二阶移动平均过程,可以表示为MA(2)。
通用表示法是MA(q)。在上面的示例中,q = 2。
让我们在Python中模拟此过程。具体来说,我们将模拟以下过程:
这是一个二阶移动平均过程,我们指定了权重。您可以随意更改权重,并对参数进行试验。
我们从指定滞后开始,我们用的是2。
ar2 = np.array([2])
然后,我们指定权重,权重为[1、0.9、0.3]。
ma2 = np.array([1, 0.9, 0.3])
最后,我们模拟该过程并生成1000个数据点:
MA2_process = ArmaProcess(ar2, ma2).generate_sample(nsample=1000)
现在,让我们可视化该过程及其相关图:
plt.figure(figsize=[10, 7.5]); # Set dimensions for figureplt.plot(MA2_process)plt.title('Moving Average Process of Order 2')plt.show()plot_acf(MA2_process, lags=20);
由于噪声是随机产生的,因此您的图形可能与我的不同。但相关图应与下图类似:
正如您所注意到的,相关性在滞后2之前是显著的。这很有意义,因为我们指定了滞后为2。
这意味着您可以使用相关图来推断时间序列的滞后。如果您看到滞后q之后相关性并不显著,那么您可以将时间序列建模为MA(q)过程。
最后
通过本文,您了解了随机游走是什么以及如何对其进行模拟。此外,您还学习了移动平均过程,并了解了如何对其建模。
matlab 指定区域随机游走_了解随机游走模型和移动平均过程(Python)相关推荐
- matlab 指定区域随机游走_统计学原理与matlab——(5)正态分布
统计学原理中的matlab应用--方差,平均数,二项分布,排列组合,几何与超几何分布.正态分布.... 这一篇文章主要聊一聊matlab与正态分布那些事儿. 说正态分布之前,先区分一组概念:离散型随机 ...
- python随机生成字符串_如何随机生成大写字母和数字组成的字符串
需求 随机生成定长的大写字母和数字组合. 实现#!/usr/bin/env python # -*- coding:utf-8 -*- import random def getRandomSet(b ...
- python 梦幻西游手游脚本_梦幻西游手游抢怪有什么技巧 梦幻西游抢怪技巧
今天小编为大家带来了梦幻西游手游抢怪有什么技巧 梦幻西游抢怪技巧,感兴趣的朋友们可以跟着小编去下文了解一下哦! 游戏中的资源就在那里不多不少,你抢或不抢都随你,但是你能得多少就不能随你啦!在游戏中抢怪 ...
- 随机网络构建_构建随机报价机
随机网络构建 by Ayo Isaiah 通过Ayo Isaiah 构建随机报价机 (Building a Random Quote Machine) I really wasn't entirely ...
- 开启事物_《原神》全新角色甘雨登场 「浮生孰来」活动祈愿开启_新闻资讯_最新手游时评_原创手游资讯...
璃月地区的方方面面,都由「璃月七星」治理.由他们决议得出的策略牵动着市场的走向,因此每一条都要慎重再三.而会议之后,是谁摘录出要点,编入长长的条例中的,人们便不得而知了. 璃月的千头万绪,就是这样被七 ...
- js 随机1-10随机数_寻找随机的错误-一个真实的故事
js 随机1-10随机数 几周前,我完成了RapidFTR开源项目的错误查找 ,这花了我三个晚上. 我认为可能值得分享狩猎的故事. 本文将介绍我的工作. 我将概述我的旅程,以便真正找到正在发生的事情的 ...
- mysql外建名是随机的吗_创建角色随机名字(mysql抽取随机记录)和mysql游标的使用_MySQL...
最近在开发中遇到了一些问题,在此记录一下解决的方法,已作备忘. 1.现在创建游戏角色的时候,基本上都是支持角色名字随机的,以前此功能在客户端用代码实现,然后向服务器请求并验证,后来发现有时候连续几次都 ...
- wx-jssdk中的config走完没走_跃叔故事(走西口):冻死鬼
我告别老赵坐车出戈壁滩的时候遇到了一场沙尘暴,不是想象中那么恐怖,同行的人说每年都会这样的日子,这算是最常见的那种,不用担心.糟糕的天气会影响人的心情,昏黄的天空让赶路的心思也变得无聊了,最后我决定在 ...
- 神雕侠侣手游服务器维修到多久,神雕侠侣手游常见问题_神雕侠侣手游问答_疑难解答_九游手机游戏...
[神雕侠侣]节日礼包开放领取啦,九游为大家提供了大量[神雕侠侣]节日礼包,欢迎大家点击领取[神雕侠侣]节日礼包,大家不要错过哦!~ 使用期限:202...[详情] <神雕侠侣>手游将在8月 ...
最新文章
- php自动加载类与路由,PHP实现路由与类自动加载步骤详解
- vue 背景透明度_一款媒体小白喜爱的视频编辑软件,vue视频编辑APP,想学就来...
- CTFshow php特性 web115
- tensorflow常见函数
- 如何通过js处理相同时间的信息整合到一起的问题
- python中convert函数用法_Python Pandas DataFrame.tz_convert用法及代码示例
- 下面哪个进制能表述 13*16=244是正确的?)[中国台湾某计算机硬件公司V2010年5月面试题]...
- python 三维曲线拟合_基于三维数据和参数的Scipy曲线拟合
- 好吧,你说简单就简单,但简单的事,不要变成本能,要常思常变
- C++认识容器的迭代器
- Qt 插件学习(一)
- [WPF]c#调用默认浏览器打开网址
- 关于MAC安装yarn
- twitter、facebook、pinterest、linkedin 分享代码
- 第三方可视化数据分析图表Pyecharts(下载保存图片(生成的html图片)、zip函数(将数据转换为列表加元组的格式)、南丁格尔玫瑰图、双y轴可视化、饼形图和环形图)
- 电脑怎么搜索文本内容? 电脑快速搜索内容的技巧
- 凡泰极客成为W3C成员并加入MiniApps工作组
- 基金股市理财经验分享(个人体会)
- [Git]如何撤销上次commit
- 【征集+投票】微信投票小程序,公众号投票,免费创建投票打分平台,免费微信投票平台、微信投票小程序、微信群投票