【机器学习之---数学】随机游走

every blog every motto: You can do more than you think.
https://blog.csdn.net/weixin_39190382?type=blog

0. 前言

随机游走

1. 概念

1.1 例1

在你的饮食俱乐部度过了一个富有成效的晚上后,你在不太清醒的状态下离开了。因此,你会醉醺醺地在展望大道上闲逛。首先,你可能会随机向右走一步。然后,忘记你在哪里,你可能会向左转一步。完全迷失了方向,你又随机迈出了一步…例如,对于一些步骤,您可能最终会像这样走

image-5

让我们把它当作一个随机过程来建模。饮食俱乐部的位置是S_0(你在零步后的初始位置)。在第k步中,移动的量为 X k X_k Xk,可以是正的,也可以是负的(取决于向左还是向右移动)。n步后的位置 S n S_n Sn

S n = S 0 + x 1 + . . . + X n S_n = S_0+x_1+...+X_n Sn=S0+x1+...+Xn

当你喝醉的时候,你会完全忘记你每走一步都在哪里。因此,我们对步骤 X 1 X_1 X1, X 2 X_2 X2,…进行建模。为独立同分布(i.i.d)随机变量,与初始位置 S 0 S_0 S0无关。这给了我们第一个随机漫步的例子

1.2 例2

假设你和你的朋友玩下面的游戏。在每一轮中,每个玩家下注1美元,并独立地投掷一枚均匀硬币。如果硬币正面朝上,那么你赢了,也就是说,你拿回你的1美元,同时你也得到你朋友的1美元。如果硬币反面出现,那么你的朋友就赢了钱。你反复玩这个游戏,直到你们中的任何一个人决定是时候减少损失了

为了将你在游戏中的财富建模为一个随机过程,让我们假设你以50美元的初始财富开始游戏。请注意,你在每一轮的净收益是+1或1。让Xi成为你在第i轮的净收益。因为硬币是公平的,你输或赚1美元的概率是一样的

P ( X i = 1 ) = P ( X i = − 1 ) = 1 2 P(X_i =1) = P(X_i=-1) = \frac{1}{2} P(Xi=1)=P(Xi=1)=21

我们假设, S 0 S_0 S0, X 1 X_1 X1, X 2 X_2 X2,…是独立的。因此,第n轮之后的财富 S n S_n Sn由,

S n = S 0 + x 1 + . . . + X n S_n = S_0+x_1+...+X_n Sn=S0+x1+...+Xn

一个可能的结果如下:

image-6

在这两个例子中,我们从起始点S0开始,制作了一系列独立同分布(i.i.d)步骤 X k X_k Xk。这样的过程被称为随机漫步。在第二个例子中,随机漫步只产生大小为1的步的情况特别重要;这叫做简单随机漫步

1.3 定义

定义如下:

随机游走 S n , n ≥ 0 {S_n},n\geq0 Sn,n0为随机过程

S n = S 0 + x 1 + . . . + X n S_n = S_0+x_1+...+X_n Sn=S0+x1+...+Xn
其中, X 1 , X 2 , . . . X_1,X_2,... X1,X2,...是独立同分布(i.i.d)随机变量,与 S 0 S_0 S0无关。

2. 边界

2.1 会达到边界吗?

我们必须研究的第一个问题是,随机漫步是否真的会碰到边界 { a , b } \{a, b\} {a,b},也就是说,是否 T < ∞ T< \infin T< 。当然,到达边界的概率原则上可能取决于起始点S0。因此,让我们引入函数:

f ( i ) : = P ( T = ∞ , S 0 = i ) , a ≤ i ≤ b f(i) := P(T=\infin, S_0=i) , \qquad a\leq i \leq b f(i):=P(T=,S0=i),aib

也就是说, f ( i ) f(i) f(i) 是我们从未碰到边界a或b中的任何一个的概率 (换句话说,随机漫步永远严格地停留在这两个边界之间),假设我们从点 S 0 = i S0 = i S0=i 开始。我们想要计算这些概率。为此,我们将引入一个非常有用的工具来研究随机漫步的行为:第一步分析(the first-step analysis)

让我们首先解释一下第一步分析背后的思想。假设我们从点 a < i < b a < i < b a<i<b严格地在边界之间开始。那么我们能到达边界的最早时间是在随机漫步的一步之后。此外,由于这是一个简单的对称随机漫步,在一步之后,我们可能到达的位置只有两个,每个位置出现的概率都是相等的:

image

因为随机漫步的每一步都是独立同分布的,所以随机漫步在每一步之后都会在当前位置重新开始。举个例子,如果在第一步中我们从 i → i + 1 i \rightarrow i+1 ii+1开始,那么从这一点开始,随机漫步就像从i+1开始的随机漫步一样。特别地,第一步是 i → i + 1 i \rightarrow i + 1 ii+1,那么之后不到达边界的概率是f (i + 1),对于其他可能的步骤也是如此

image-1

如图,我们可以得到:

f ( i ) = P ( T = ∞ ∣ S 1 = i + 1 , S 0 = i ) P ( S 1 = i + 1 ∣ S 0 = i ) + P ( T = ∞ ∣ S 1 = i − 1 , S 0 = i ) P ( S 1 = i − 1 ∣ S 0 = i ) f(i)= P(T=\infin|S_1=i+1, S_0=i)P(S_1=i+1|S_0=i)+P(T=\infin|S_1=i-1, S_0=i)P(S_1=i-1|S_0=i) f(i)=P(T=∞∣S1=i+1,S0=i)P(S1=i+1∣S0=i)+P(T=∞∣S1=i1,S0=i)P(S1=i1∣S0=i)

其中,
P ( S 1 = i + 1 ∣ S 0 = i ) = 1 2 P ( S 1 = i − 1 ∣ S 0 = i ) = 1 2 P(S_1=i+1|S_0=i)= \frac{1}{2} \qquad \qquad P(S_1=i-1|S_0=i)= \frac{1}{2} P(S1=i+1∣S0=i)=21P(S1=i1∣S0=i)=21

同时,有
P ( T = ∞ ∣ S 1 = i + 1 , S 0 = i ) = f ( i + 1 ) P ( T = ∞ ∣ S 1 = i − 1 , S 0 = i ) = f ( i − 1 ) \begin{align} P(T=\infin|S_1=i+1, S_0=i)=f(i+1) \\ P(T=\infin|S_1=i-1, S_0=i)=f(i-1) \end{align} P(T=∞∣S1=i+1,S0=i)=f(i+1)P(T=∞∣S1=i1,S0=i)=f(i1)

因此,

f ( i ) = 1 2 f ( i + 1 ) + 1 2 f ( i − 1 ) , a < i < b f(i)= \frac{1}{2}f(i+1)+\frac{1}{2}f(i-1),\qquad a < i < b f(i)=21f(i+1)+21f(i1),a<i<b

另一方面,假设我们从S0 = a或S0 = b开始,那么我们一开始就已经从边界开始了,所以我们从未到达边界的概率为零。特别地,这意味着f (a) = f (b) = 0。因此,我们的第一步分析使我们得到以下等式

{ f ( i ) = 1 2 f ( i + 1 ) + 1 2 f ( i − 1 ) , a < i < b f ( a ) = f ( b ) = 0 \begin{cases} f(i)= \frac{1}{2}f(i+1)+\frac{1}{2}f(i-1),\qquad a < i < b \\ f(a)=f(b)=0 \end{cases} {f(i)=21f(i+1)+21f(i1),a<i<bf(a)=f(b)=0

通过变换方程我们知道,f(i)位于f(i-1)和f(i+1)中间位置。同时,在一条直线上。

image-2

所以,f(i)是直线。又,f(a)=f(b)=0,可得f(i)=0:

image-3

P ( T = ∞ ∣ S 0 = i ) = 0 , f o r a l l a ≤ i ≤ b P(T=\infin|S_0=i)=0, \qquad for all a \leq i \leq b P(T=∞∣S0=i)=0,forallaib

所以,随机漫步最终总是到达a或b,而不管它的起点是什么。

2.2 会达到哪一个边界

既然我们知道我们最终总会碰到边界a或b中的一个,我们可能想知道我们先碰到哪一个?例如,在赌博问题中, T < ∞ T <\infin T<表示你最终要么赢得100美元要么破产。在实践中,您可能对这些结果中哪一种更有可能发生非常感兴趣。

为了深入了解这个问题,让我们用数学术语来表达它。在我们到达a或b的第一时间T,值 S T S_T ST显然必须是a或b。我们感兴趣的是计算 S T = b S_T = b ST=b的概率,比如(在这种情况下,随机行走在到达a之前到达b)。由于这个概率可能再次取决于我们开始随机行走的位置,让我们定义这个函数

r ( i ) : = P ( S T = b ∣ S 0 = i ) r(i) := P(S_T=b|S_0 = i) r(i):=P(ST=bS0=i)

为了计算r(i),我们可以再次使用第一步分析。如果a < i < b,那么一步之后我们可以最早到达其中一个边界。一步之后,我们在i + 1或i -1处的概率相等,在a之前到达b的概率分别由r(i + 1)和r(i - 1)给出。正如上文所述,我们因此发现

r ( i ) = 1 2 r ( i + 1 ) + 1 2 r ( i − 1 ) , a < i < b r(i) = \frac{1}{2}r(i+1)+\frac{1}{2}r(i-1), \qquad a < i < b r(i)=21r(i+1)+21r(i1),a<i<b

这里不同的是边界条件。在这种情况下,如果我们从S0 = b开始,那么我们显然会在a之前到达b(因为我们已经从b开始了!),所以r(b) = 1。另一方面,如果S0 = a,那么显然a在b之前,所以r(b) = 0。因此,我们得到以下线性方程

{ r ( i ) = 1 2 r ( i + 1 ) + 1 2 r ( i − 1 ) , a < i < b f ( a ) = 0 , f ( b ) = 1 \begin{cases} r(i)= \frac{1}{2}r(i+1)+\frac{1}{2}r(i-1),\qquad a < i < b \\ f(a)=0, \quad f(b)=1 \end{cases} {r(i)=21r(i+1)+21r(i1),a<i<bf(a)=0,f(b)=1

同样,我们可以通过变换方程将r(i)表示为一条直线。同时r(a)=0,r(b)=1,如下:

image-4

最终公式如下:

P ( S T = b ∣ S 0 = i ) = r ( i ) = i − a b − a , a ≤ i ≤ b P(S_T=b|S_0=i)= r(i) = \frac{i-a}{b-a}, \qquad a \leq i \leq b P(ST=bS0=i)=r(i)=baia,aib

我们越接近b,我们越有可能在a之前碰到b,而我们越接近a,我们越有可能在b之前碰到a,这和直觉上一致。

上面的赌博游戏,直到你破产或达到100美元的财富。让它成为游戏结束的时间。我们已经证明了 T < ∞ T <\infin T<(游戏最终结束),你获得100美元的概率是:

P ( S T = 100 ∣ S 0 = i ) = i 100 , 0 ≤ i ≤ 100 P(S_T=100|S_0=i) = \frac{i}{100}, \quad 0 \leq i \leq 100 P(ST=100∣S0=i)=100i,0i100

即,如果你开始有i美元,你最终有100美元的概率是i/100。

让我们计算一下期望,

E ( S T ∣ S 0 = i ) = 0 ∗ P ( S T = 0 ∣ S 0 = i ) + 100 ∗ P ( S T = 100 ∣ S 0 = i ) = i E(S_T|S_0=i) = 0*P(S_T=0|S_0=i)+100*P(S_T=100|S_0=i) = i E(STS0=i)=0P(ST=0∣S0=i)+100P(ST=100∣S0=i)=i

最终,你的拥有的钱和初始状态一样。

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

【opencv】实时位姿估计(real_time_pose_estimation)—3D模型注册

相机成像原理图 物体网格、关键点&#xff08;局内点、局外点&#xff09;图像 box.ply resized_IMG_3875.JPG 主程序main_registration.cpp 主要实现了利用OpenCV库进行3D模型的注册。主要步骤包括加载3D网格模型、使用鼠标事件选择对应的3D点进行2D到3D的注册、利用solvePnP算…

在django中使用kindeditor出现转圈问题

在django中使用kindeditor出现转圈问题 【一】基础检查 【1】前端检查 确保修改了uploadJson的默认地址 该地址需要在路由层有映射关系 确认有加载官方文件 kindeditor-all-min.js确保有传递csrfmiddlewaretoken 或者后端关闭了csrf验证 <textarea name"content&qu…

无人驾驶矿卡整体解决方案(5g物联网通信方案)

​无人驾驶矿卡是智能矿山的重要组成部分,通过远程操控替代人工驾驶,可以显著提高采矿效率和作业安全性。但要实现无人驾驶矿卡,需要依赖于可靠高效的通信网络,来传输现场视频、控制指令和运行数据。以下是某大型煤矿在部署无人驾驶矿卡时,所采用的星创易联物联网整体解决方案。…

如何区分模型文件是稳定扩散模型和LORA模型

区分模型文件是否为稳定扩散模型&#xff08;Stable Diffusion Models&#xff09;或LORA模型&#xff08;LowRank Adaptation&#xff09;通常需要对模型的结构和内容有一定的了解。以下是一些方法来区分这两种模型文件&#xff1a; 1. 文件格式和结构 稳定扩散模型&#xff1…

词根词缀基础

一&#xff0e;词根词缀方法&#xff1a; 1. 类似中文的偏旁部首&#xff08;比如“休”单人旁木→一个人靠木头上休息&#xff09; 2. 把单词拆分后&#xff0c;每一个部分都有它自己的意思&#xff0c;拼凑在一起就构成了这个单词的意思 3. 一个规律&#xff0c;适用大部分…

基于nodejs+vue多媒体素材管理系统python-flask-django-php

该系统采用了nodejs技术、express 框架&#xff0c;连接MySQL数据库&#xff0c;具有较高的信息传输速率与较强的数据处理能力。包含管理员、教师和用户三个层级的用户角色&#xff0c;系统管理员可以对个人中心、用户管理、教师管理、资源类型管理、资源信息管理、素材类型管理…

论文阅读-《Lite Pose: Efficient Architecture Design for 2D Human Pose Estimation》

摘要 这篇论文主要研究了2D人体姿态估计的高效架构设计。姿态估计在以人为中心的视觉应用中发挥着关键作用&#xff0c;但由于基于HRNet的先进姿态估计模型计算成本高昂&#xff08;每帧超过150 GMACs&#xff09;&#xff0c;难以在资源受限的边缘设备上部署。因此&#xff0…

(三)Ribbon负载均衡

1.1.负载均衡原理 SpringCloud底层其实是利用了一个名为Ribbon的组件&#xff0c;来实现负载均衡功能的。 1.2.源码跟踪 为什么我们只输入了service名称就可以访问了呢&#xff1f;之前还要获取ip和端口。 显然有人帮我们根据service名称&#xff0c;获取到了服务实例的ip和…

GitLab更新失败(Ubuntu)

在Ubuntu下使用apt更新gitlab报错如下&#xff1a; An error occurred during the signature verification.The repository is not updated and the previous index files will be used.GPG error: ... Failed to fetch https://packages.gitlab.com/gitlab/gitlab-ee/ubuntu/d…

Leetcode 3.26

Leetcode Hot 100 一级目录1.每日温度 堆1.数组中的第K个最大元素知识点&#xff1a;排序复杂度知识点&#xff1a;堆的实现 2.前 K 个高频元素知识点&#xff1a;优先队列 一级目录 1.每日温度 每日温度 思路是维护一个递减栈&#xff0c;存储的是当前元素的位置。 遍历整个…

web学习笔记(四十五)Node.js

目录 1. Node.js 1.1 什么是Node.js 1.2 为什么要学node.js 1.3 node.js的使用场景 1.4 Node.js 环境的安装 1.5 如何查看自己安装的node.js的版本 1.6 常用终端命令 2. fs 文件系统模块 2.1引入fs核心模块 2.2 读取指定文件的内容 2.3 向文件写入指定内容 2.4 创…

app自动化-Appium学习笔记

使用Appium&#xff0c;优点&#xff1a; 1、支持语言比较多&#xff0c;例如&#xff1a;Java、Python、Javascript、PHP、C#等语言 2、支持跨应用&#xff08;windows、mac、linux&#xff09; 3、适用平台Android、iOS 4、支持Native App(原生app)、Web App、Hybird App…

canvas画图写文字,有0.5像素左右的位置偏差,无解决办法,希望有知道问题的大神告知一下

提示&#xff1a;canvas画图写文字 文章目录 前言一、写文字总结 前言 一、写文字 test.html <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><meta name"viewport" content"widthdevice-widt…

Fragment 与 ViewPager的联合应用(2)

5.创建底部布局bottom_layout <LinearLayout xmlns:android"http://schemas.android.com/apk/res/android"android:orientation"horizontal"android:layout_width"match_parent"android:layout_height"55dp"android:background&qu…

【算法】求最大公约数和最小公倍数

题目 输入两个数&#xff08;空格隔开&#xff09;分2行输出他们的最大公因数和最小公倍数 原理 辗转相除法计算最大公约数 将两个数中较大的数除以较小的数&#xff0c;并将较小的数作为除数&#xff0c;较大的数作为被除数。计算余数。若余数为零&#xff0c;则较小的数即…

深入探索MySQL高阶查询语句的艺术与实践

目录 引言 一、条件查询 &#xff08;一&#xff09;比较运算符查询 1.使用匹配符号查询 2.范围查找 &#xff08;二&#xff09;逻辑运算符 二、关键字排序 三、分组与聚合函数 四、限制查询 五、别名 &#xff08;一&#xff09;设置列别名 &#xff08;二&#x…

Dockerfile和Docker-compose

一、概述 Dockerfile和Docker Compose是用于构建和管理 Docker 容器的两个工具&#xff0c;但它们的作用和使用方式不同。 Dockerfile Dockerfile 是一个文本文件&#xff0c;用于定义 Docker 镜像的构建规则。它包含一系列指令&#xff0c;如 FROM&#xff08;指定基础镜像…

python(django)之单一接口管理功能后台开发

1、创建数据模型 在apitest/models.py下加入以下代码 class Apis(models.Model):Product models.ForeignKey(product.Product, on_deletemodels.CASCADE, nullTrue)# 关联产品IDapiname models.CharField(接口名称, max_length100)apiurl models.CharField(接口地址, max_…

uniapp微信小程序_computed_计算BMI

一、computed的用法还有它是什么&#xff1f; 首先它叫计算属性&#xff0c;顾名思义他是用来计算属性&#xff0c;计算你在data模板上定义的属性&#xff08;其实在插值表达式也能直接计算但是首先太长了在{{}}里面写那么多不好看&#xff0c;还有其他特点我在下面一起说&…

jupyter notebook导出含中文的pdf(LaTex安装和Pandoc、MiKTex安装)

用jupyter notebook导出pdf时&#xff0c;因为报错信息&#xff0c;需要用到Tex nbconvert failed: xelatex not found on PATH, if you have not installed xelatex you may need to do so. Find further instructions at https://nbconvert.readthedocs.io/en/latest/install…