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

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

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

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

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

相关文章

数据结构(五)单链表专题

在开始之前&#xff0c;我先来给大家讲一下顺序表与链表的区别&#xff1a; 它们在堆上存储的差异&#xff1a; 我们可以很容易的知道&#xff0c;循序表是连续的有序的&#xff0c;但链表是杂乱的&#xff0c;它们通过地址彼此联系起来。 1. 链表的概念及结构 概念&#xff1…

智慧交通(代码实现案例)

1.项目简介 目标: 了解智慧交通项目的架构知道智慧交通项目中的模块能够完成智慧交通项目的环境搭建 该项目是智慧交通项目&#xff0c;通过该项目掌握计算机视觉的方法在交通领域的相关应用&#xff0c;包括车道线检测的方法&#xff0c;多目标车辆追踪及流量统计方法&#…

Linux Tomcat的服务器如何查看接口请求方式?

问题描述 最近在和安卓开发对接接口&#xff0c;遇到一个接口总是报405错误&#xff0c;有对接经验的开发应该都知道是请求方式不对&#xff0c;假如接口定义为POST请求的&#xff0c;但是客户端却用GET请求&#xff0c;这时候就会报这个错误。Android客户端那边使用xUtils框架…

【小白入门篇3】还是GPT4更香

上一节文章《【小白入门篇2】总有一款AI工具适合你》介绍了很多ai产品给大家&#xff0c;有同学私信我&#xff0c;国内工具还是比较差&#xff0c;还是想用gpt4模型。这个章节介绍一些gpt4工具给大家, 其中大部分都只有一些免费的次数, 而且都需要kx上网才能访问。 OpenAI ch…

浙大版《C语言程序设计(第4版)》题目集-练习4-7 求e的近似值

自然常数 e 可以用级数 1 1 / 1 ! 1 / 2 ! ⋯ 1 / n ! ⋯ 11/1!1/2!⋯1/n!⋯ 11/1!1/2!⋯1/n!⋯来近似计算。本题要求对给定的非负整数 n&#xff0c;求该级数的前 n1 项和。 输入格式: 输入第一行中给出非负整数 n&#xff08;≤1000&#xff09;。 输出格式: 在一行…

记录一下安装ubuntu子系统的pycharm遇到的问题

sudo su #切换为root用户获取管理员权限用于新建用户 adduser username #新建用户&#xff08;例如用户名为username&#xff09; adduser username sudo #将用户添加到 sudo 组同样遇到这个问题&#xff0c;解决方法是&#xff1a;先新建一个用户名&#xff0c;然后再切换到这…

Android 性能优化实例分享-内存优化 兼顾效率与性能

背景 项目上线一段时间后,回顾重要页面 保证更好用户体验及生产效率&#xff0c;做了内存优化和下载导出优化&#xff0c;具体效果如最后的一节的表格所示。 下面针对拍摄流程的两个页面 预览页 导出页优化实例进行介绍&#xff1a; 一.拍摄前预览页面优化 预览效果问题 存在…

Quartus II仿真出现错误

ModelSim executable not found in D:/intelFPGA/18.0/quartus/bin64/modelsim_ase/win32aloem/ Error. 找不到modelsim地址&#xff0c;原来是我下载了.exe,但没有双击启动安装ase文件夹呀&#xff01;&#xff01;&#xff01;&#xff01;晕&#xff0c;服了我自己

Python7:接口自动化学习1 RPC

API&#xff08;Application Programmming Interface&#xff09; 应用编程接口&#xff0c;简称“接口” 接口&#xff1a;程序之间约定的通信方法 特点&#xff1a;约定了调用方法&#xff0c;以及预期的行为&#xff0c;但是不透露具体细节 意义&#xff1a;程序能解耦&…

【No.20】蓝桥杯简单数论下|寻找整数|素数的判断|笨小猴|最大最小公倍数|素数筛|埃氏筛|欧氏线性筛|质数|分解质因子(C++)

寻找整数 【题目描述】 有一个不超过 1 0 1 7 10^17 1017的正整数n&#xff0c;知道这个数除以2至49后的余数如下表所示&#xff0c;求这个正整数最小是多少 解法一&#xff1a;模拟 暴力法&#xff1a;一个个检验 1 … 1 0 17 1\dots 10^{17} 1…1017的每个数 由于这个数n…

【pytest、playwright】构建POM项目,以及解决登录问题,allure环境问题

目录 前言 1、文件目录 2、安装依赖 3、POM项目实战-案例&#xff1a;打开指定页面 目录结构&#xff1a; pages中的代码&#xff1a; cases中的代码&#xff1a; 4、解决登录问题 问题&#xff1a; 解决方案&#xff1a; 获取登录的用户信息&#xff08;cookie&a…

DasViewer电脑客户端打开文件夹时,一直显示崩溃,该怎么解决?

问题如图 如若用的是DasViewer V3.2.4Beta版本&#xff0c;可以换回3.2.1版本进行尝试。 DasViewer是由大势智慧自主研发的免费的实景三维模型浏览器,采用多细节层次模型逐步自适应加载技术,让用户在极低的电脑配置下,也能流畅的加载较大规模实景三维模型,提供方便快捷的数据浏…

发送请求- header配置

请求头里是客户端的要求&#xff0c;把你的诉求告诉服务端&#xff0c;服务端按照你的要求返回数据 &#xff0c; 请求header需要严格全配置&#xff0c;把请求header全部传入&#xff0c;不能频繁访问&#xff0c;让后端知道它是正常请求 一般只配置User-Agent和Content Typ…

云原生靶场kebernetesGoat、Metarget

靶场 文章目录 靶场kebernetesGoat靶场安装Docker in DockerSSRF漏洞容器逃逸到主系统Docker CIS 基线分析Kubernetes CIS 安全基线分析分析被部署挖矿软件的容器镜像获取环境信息Hidden in layersRBAC最低权限配置错误使用 Sysdig Falco 进行运行时安全监控和检测 Metarget ke…

iOS App 上架指南及关键

引言 上架App Store是将iOS应用提交申请并上线的过程&#xff0c;旨在让应用在App Store上展示&#xff0c;吸引用户并获取流量。本文将介绍iOS上架的整体流程&#xff0c;并提供一些建议和注意事项。 一、iOS上架的整体流程 1. 申请开发者账号 首先&#xff0c;需要申请苹…

DevSecOps平台架构系列-亚马逊云AWS DevSecOps平台架构

目录 一、概述 二、AWS DevSecOps实施原则 2.1 尽早采用安全测试&#xff0c;加速问题反馈 2.2 优先考虑预防性安全控制 2.3 部署检测性安全控制时&#xff0c;确保有与之互补的响应性安全控制 2.4 安全自动化 2.5 总结 三、AWS DevSecOps关键组件 3.1 关键组件 3.2 关…

新能源汽车充电桩主板各模块成本占比解析

汽车充电桩主板是汽车充电桩的重要组件&#xff0c;主要由微处理器模块、通信模块、控制模块、安全保护模块、传感器模块等多个模块构成。深入探究各模块在总成本中的比重&#xff0c;我们可以更好地优化成本结构、提高生产效率&#xff0c;并为未来的技术创新和市场需求变化做…

TypeScript【入门】

TypeSript ts是js的一个超集&#xff0c;解决js类型支持 js: 动态类型 ts:静态类型 代码执行的顺序是先编译在执行 发现错误更早 为什么要添加类型支持&#xff1f; 优势有什么&#xff1f; 全局安装 npm i -g typescript tsc -v 简化运行ts的步骤 typescript常用类型 类型…

C++类和对象、面向对象编程 (OOP)

文章目录 一、封装1.抽象、封装2.类和对象(0)学习视频(1)类的构成(2)三种访问权限(3)struct和class的区别(4)私有的成员变量、共有的成员函数(5)类内可以直接访问私有成员&#xff0c;不需要经过对象 二、继承三、多态1.概念2.多态的满足条件3.多态的使用条件4.多态原理剖析5.纯…

MySQL详细教程

文章目录 前言一、数据库管理1.查看已有的数据库2.创建数据库3.删除数据库4.进入数据库 二、 数据表管理1.查看当前数据库下的所有数据表2.创建数据表3.删除表4.查看表结构 三、常用数据类型1.整型tinyintintbigint 2.浮点型floatdoubledecimal 3.字符型char(m)varchar(m)textm…