当前位置: 首页 > news >正文

伊藤积分(Ito Integral):随机世界中的积分魔法

伊藤积分(Ito Integral):随机世界中的积分魔法

在研究随机微分方程(SDE)和布朗运动时,伊藤积分(Ito Integral)是一个绕不开的关键概念。它是处理布朗运动随机项 ( d W ( t ) dW(t) dW(t) ) 的数学工具,与传统的黎曼积分截然不同。对于深度学习研究者来说,伊藤积分不仅是 SDE 的核心,也是理解扩散模型(如 DDPM)中逆过程的基础。本篇博客将以直观的语言,面向具有一定数学和深度学习背景的读者,介绍伊藤积分的定义、特点及其在随机建模中的意义。

为什么需要伊藤积分?

普通的常微分方程(ODE)可以用黎曼积分求解,例如:
x ( t ) = x ( 0 ) + ∫ 0 t f ( s , x ( s ) ) d s x(t) = x(0) + \int_0^t f(s, x(s)) ds x(t)=x(0)+0tf(s,x(s))ds
这里的积分是对确定性函数 ( f ( s , x ( s ) ) f(s, x(s)) f(s,x(s)) ) 在时间 ( [ 0 , t ] [0, t] [0,t] ) 上的累积。但在 SDE 中,如:
d x ( t ) = f ( t , x ) d t + g ( t , x ) d W ( t ) dx(t) = f(t, x) dt + g(t, x) dW(t) dx(t)=f(t,x)dt+g(t,x)dW(t)
第二个项 ( g ( t , x ) d W ( t ) g(t, x) dW(t) g(t,x)dW(t) ) 涉及布朗运动 ( W ( t ) W(t) W(t) ) 的微分 ( d W ( t ) dW(t) dW(t) )。由于布朗运动路径连续但无处可微,传统的黎曼积分无法直接处理这种“随机抖动”。伊藤积分应运而生,为随机项提供了一种特殊的积分定义。

伊藤积分的直观定义

伊藤积分的目标是对形式 ( ∫ 0 t g ( s ) d W ( s ) \int_0^t g(s) dW(s) 0tg(s)dW(s)) 的积分赋值,其中 ( g ( s ) g(s) g(s) ) 是一个随机过程(通常依赖 ( W ( s ) W(s) W(s) )),( d W ( s ) dW(s) dW(s) ) 是布朗运动的微分增量。直观上,它类似于黎曼积分的分段求和,但有关键区别。

黎曼积分的类比

在黎曼积分中,( ∫ 0 t f ( s ) d s \int_0^t f(s) ds 0tf(s)ds) 被近似为:
∑ i = 0 n − 1 f ( s i ) ( s i + 1 − s i ) \sum_{i=0}^{n-1} f(s_i) (s_{i+1} - s_i) i=0n1f(si)(si+1si)
其中 ( s i s_i si ) 是时间区间 ( [ 0 , t ] [0, t] [0,t] ) 的划分点,( s i + 1 − s i = Δ t s_{i+1} - s_i = \Delta t si+1si=Δt )。

伊藤积分的构造

对于 ( ∫ 0 t g ( s ) d W ( s ) \int_0^t g(s) dW(s) 0tg(s)dW(s)),我们类似地构造分段和:
∫ 0 t g ( s ) d W ( s ) ≈ ∑ i = 0 n − 1 g ( s i ) [ W ( s i + 1 ) − W ( s i ) ] \int_0^t g(s) dW(s) \approx \sum_{i=0}^{n-1} g(s_i) [W(s_{i+1}) - W(s_i)] 0tg(s)dW(s)i=0n1g(si)[W(si+1)W(si)]

  • ( s i s_i si ):时间划分点,如 ( s i = i ⋅ t n s_i = i \cdot \frac{t}{n} si=int )。
  • ( W ( s i + 1 ) − W ( s i ) W(s_{i+1}) - W(s_i) W(si+1)W(si) ):布朗运动在 ( [ s i , s i + 1 ] [s_i, s_{i+1}] [si,si+1] ) 的增量,记作 ( Δ W i ∼ N ( 0 , s i + 1 − s i ) \Delta W_i \sim \mathcal{N}(0, s_{i+1} - s_i) ΔWiN(0,si+1si))。

当划分 ( n → ∞ n \to \infty n ) 时,这个和的极限定义了伊藤积分:
∫ 0 t g ( s ) d W ( s ) = lim ⁡ n → ∞ ∑ i = 0 n − 1 g ( s i ) Δ W i \int_0^t g(s) dW(s) = \lim_{n \to \infty} \sum_{i=0}^{n-1} g(s_i) \Delta W_i 0tg(s)dW(s)=nlimi=0n1g(si)ΔWi

关键区别:前点选择
  • 黎曼积分:( f ( s ) f(s) f(s) ) 的取值可以在区间 ( [ s i , s i + 1 ] [s_i, s_{i+1}] [si,si+1] ) 内任意点(如中点),极限结果一致。
  • 伊藤积分:( g ( s ) g(s) g(s) ) 固定取左端点 ( s i s_i si )(即前一时刻的值),这被称为“非预期性”(non-anticipating),因为 ( g ( s i ) g(s_i) g(si) ) 不能依赖未来的 ( W ( s i + 1 ) W(s_{i+1}) W(si+1) )。

这种前点选择让伊藤积分与布朗运动的马尔可夫性质一致,确保积分是可预测过程。

伊藤积分的特性
  1. 随机性
    由于 ( Δ W i \Delta W_i ΔWi) 是随机变量,伊藤积分的结果是一个随机变量,而不是固定值。例如:
    ∫ 0 t d W ( s ) = W ( t ) − W ( 0 ) = W ( t ) \int_0^t dW(s) = W(t) - W(0) = W(t) 0tdW(s)=W(t)W(0)=W(t)
    它本身就是布朗运动。

  2. 非平滑性
    伊藤积分的路径继承了布朗运动的连续但不可导特性,抖动剧烈。

  3. 二次变差
    对于普通积分,( ∑ ( s i + 1 − s i ) 2 → 0 \sum (s_{i+1} - s_i)^2 \to 0 (si+1si)20)(当 ( n → ∞ n \to \infty n ))。但对于伊藤积分:
    ∑ i = 0 n − 1 ( Δ W i ) 2 → t \sum_{i=0}^{n-1} (\Delta W_i)^2 \to t i=0n1(ΔWi)2t
    极限是时间长度 ( t t t ),而不是 0。这是因为布朗运动的增量方差累积随时间线性增长。

  4. 伊藤引理
    伊藤积分需要特殊的微分规则。例如,对于函数 ( F ( x ) F(x) F(x) ):
    d F ( x ( t ) ) = F ′ ( x ( t ) ) d x ( t ) + 1 2 F ′ ′ ( x ( t ) ) g ( t , x ) 2 d t dF(x(t)) = F'(x(t)) dx(t) + \frac{1}{2} F''(x(t)) g(t, x)^2 dt dF(x(t))=F(x(t))dx(t)+21F′′(x(t))g(t,x)2dt
    相比普通链式法则,多了一项二阶修正,反映了随机项的影响。

伊藤积分与 SDE

SDE 的解通常写成积分形式:
x ( t ) = x ( 0 ) + ∫ 0 t f ( s , x ( s ) ) d s + ∫ 0 t g ( s , x ( s ) ) d W ( s ) x(t) = x(0) + \int_0^t f(s, x(s)) ds + \int_0^t g(s, x(s)) dW(s) x(t)=x(0)+0tf(s,x(s))ds+0tg(s,x(s))dW(s)

  • 第一项:普通黎曼积分,计算确定性漂移。
  • 第二项:伊藤积分,处理随机扩散。

例如,几何布朗运动:
d x ( t ) = μ x d t + σ x d W ( t ) dx(t) = \mu x dt + \sigma x dW(t) dx(t)=μxdt+σxdW(t)
解为:
x ( t ) = x ( 0 ) exp ⁡ ( ( μ − σ 2 2 ) t + σ W ( t ) ) x(t) = x(0) \exp\left( (\mu - \frac{\sigma^2}{2}) t + \sigma W(t) \right) x(t)=x(0)exp((μ2σ2)t+σW(t))
其中伊藤积分 ( ∫ 0 t σ x ( s ) d W ( s ) \int_0^t \sigma x(s) dW(s) 0tσx(s)dW(s)) 通过伊藤引理推导。

在深度学习中的应用

伊藤积分在生成模型中至关重要:

  • 扩散模型(DDPM)
    • 前向过程添加噪声,逆过程用 SDE 表示:
      d x = f ( x , t ) d t + g ( t ) d W ( t ) dx = f(x, t) dt + g(t) dW(t) dx=f(x,t)dt+g(t)dW(t)
    • 逆向积分涉及伊藤积分,描述从噪声到数据的随机路径。
  • Langevin 动力学
    • NCSN 的采样公式 ( x t + 1 = x t + α s θ ( x t ) + α z t x_{t+1} = x_t + \alpha s_\theta(x_t) + \sqrt{\alpha} z_t xt+1=xt+αsθ(xt)+α zt ) 是 SDE 的离散近似,( z t z_t zt) 对应 ( d W ( t ) dW(t) dW(t) )。
总结

伊藤积分是处理布朗运动随机项 ( d W ( t ) dW(t) dW(t) ) 的数学魔法,与黎曼积分不同,它通过前点和的形式定义,适应了布朗运动的不可导性和随机性。形式上:
∫ 0 t g ( s ) d W ( s ) = lim ⁡ n → ∞ ∑ i = 0 n − 1 g ( s i ) [ W ( s i + 1 ) − W ( s i ) ] \int_0^t g(s) dW(s) = \lim_{n \to \infty} \sum_{i=0}^{n-1} g(s_i) [W(s_{i+1}) - W(s_i)] 0tg(s)dW(s)=nlimi=0n1g(si)[W(si+1)W(si)]
它不仅是 SDE 解的关键,还为扩散模型等深度学习方法提供了理论支持。对于研究者来说,理解伊藤积分就像掌握了随机世界的“积分钥匙”,打开了从噪声到数据的建模之门。


注:本文以直观解释为主,未深入严格证明,适合快速入门。

后记

2025年3月8日20点56分于上海,在Grok 3大模型辅助下完成。

相关文章:

伊藤积分(Ito Integral):随机世界中的积分魔法

伊藤积分(Ito Integral):随机世界中的积分魔法 在研究随机微分方程(SDE)和布朗运动时,伊藤积分(Ito Integral)是一个绕不开的关键概念。它是处理布朗运动随机项 ( d W ( t ) dW(t)…...

【Deepseek应用】Zotero+Deepseek 阅读和分析文献(下)

【Deepseek应用】Deepseek R1 本地部署(OllamaDockerOpenWebUI) 【Deepseek应用】ZoteroDeepseek 阅读和分析文献(上) 【Deepseek应用】ZoteroDeepseek 阅读和分析文献(下) 使用邀请码 cXfb9wOT 注册 硅基流…...

人工智能与深度学习的应用案例:从技术原理到实践创新

第一章 引言 人工智能(AI)作为21世纪最具变革性的技术之一,正通过深度学习(Deep Learning)等核心技术推动各行业的智能化进程。从计算机视觉到自然语言处理,从医疗诊断到工业制造,深度学习通过模拟人脑神经网络的层次化学习机制,实现了对复杂数据的高效分析与决策。本…...

Docker和DockerCompose基础教程及安装教程

Docker的应用场景 Web 应用的自动化打包和发布。自动化测试和持续集成、发布。在服务型环境中部署和调整数据库或其他的后台应用。从头编译或者扩展现有的 OpenShift 或 Cloud Foundry 平台来搭建自己的 PaaS 环境。 CentOS Docker 安装 使用官方安装脚本自动安装 安装命令…...

ArcGIS操作:13 生成最小外接矩阵

应用情景:筛选出屋面是否能放下12*60m的长方形,作为起降场候选点(一个不规则的形状内,判断是否能放下指定长宽的长方形) 1、面积初步筛选 Area ≥ 720 ㎡ 面积计算见 2、打开 ArcToolbox → Data Management Tools …...

Qt:事件

目录 处理事件 鼠标事件 键盘事件 定时器事件 窗口事件 虽然 Qt 是跨平台的 C 开发框架,Qt 的很多能力其实是操作系统提供的 只不过 Qt 封装了系统的 API 事件 前面学习过信号槽: 用户进行的各种操作,就可能会产生出信号,可以…...

python 程序一次启动有两个进程的问题(flask)

0. 背景 写了一个使用 flask 作为服务框架的程序,发现每次启动程序的时候,使用 ps 都能观察到两个 python 进程。 此外,这个程序占用了 GPU 资源,我发现有两个 python 进程,分别占用了完全相同的 GPU 显存 1. 原因 …...

ethtool的资料

ethtoolethtool(8) — Linux manual pageethtool(8) - Linux man pageUsing ethtool in LinuxLooking at your Linux system’s network interface with ethtoolHow to Change Speed & Duplex of Ethernet Card in Linux with ethtool CommandNVIDIA EthtoolRed Hat Enterp…...

SpringBoot过滤器(Filter)的使用:Filter接口、FilterRegistrationBean类配置、@WebFilter注释

1、过滤器(Filter)的介绍 Spring Boot 的过滤器用于对数据进行过滤处理。通过 Spring Boot 的过滤器,程序开发人员不仅可以对用户通过 URL 地址发送的请求进行过滤处理(例如:过滤一些错误的请求或者请求中的敏感词等),而且可以对服务器返回的数据进行过滤处理(例如:压…...

“此电脑”中删除WPS云盘方法(百度网盘通用)

📣此方法适用于卸载WPS云盘后,WPS云盘图标依然在此电脑中显示的问题。 原理:通过注册来进行删除 步骤: WIN键R,打开运行窗口,输入regedit命令,来打开【注册表编辑器】; 从左侧,依…...

Manus AI:开启Agent元年的ChatGPT时刻(附赠资料)

1. Manus AI:全球首个通用Agent Manus AI 是全球首个通用人工智能代理,连接思想与行动,不仅思考,还能交付成果。Manus 擅长处理工作和生活中的各种任务,帮助用户完成一切。其核心理念是“less structure, more intell…...

ChromeDriver下载 最新版本 134.0.6998.35

平时为了下个驱动,到处找挺麻烦,收集了很多无偿分享给需要的人,仅供学习和交流。 ChromeDriver及浏览器134.0.6998.35 ChromeDriver及浏览器133.0.6943.141 ChromeDriver 102.0.5005.61 ChromeDriver 105.0.5195.102 ChromeDriver 108.0…...

Sass进阶之路:@forward 的可见性控制与变量覆盖

文章目录 前言1. 转发导入2. 添加前缀3. 控制可见性4. 转发时修改默认值总结 前言 在上一篇中,我们深入探讨了 use 的使用, 也介绍了 use 在使用深层模块中的变量时具有一定的缺点。所以在本文中,我们将深入解析 forward 的核心用法。 1. 转…...

MySQL作业一

一、创建数据库 #创建数据库 mysql> create database db_ck; Query OK, 1 row affected (0.01 sec)mysql> show databases like "db_%"; ----------------- | Database (db_%) | ----------------- | db_ck | | db_system | ----------------…...

虚拟机总结| 关于虚拟机的一些配置总结

前言 每次安装新的虚拟机都需要重新在网上搜索如何配置网络,我需要写一个自己的部署步骤,增加工作效率,不用每次配置的时候再去网上去翻找。 1.只需要联网功能记录(不固定IP) 1.1 修改ifcfg-ens33 vi etc/sysconfig/network-scripts/ifcfg…...

leetcode-sql数据库面试题冲刺(高频SQL五十题)

题目: 577.员工奖金 表:Employee -------------------- | Column Name | Type | -------------------- | empId | int | | name | varchar | | supervisor | int | | salary | int | -------------------- empId 是该表中具有唯一值的列。 该表的每一行…...

OpenManus:解锁测试工程师的效率密码——实践与应用指南

随着软件行业的快速发展,测试工程师面临的挑战也日益增多:如何在有限的时间内保证产品质量、如何高效生成测试数据、如何快速定位问题根源?这些问题直接影响到产品上线的节奏和用户体验。而在这一背景下,开源项目 OpenManus 的出现…...

Mybatis中的设计模式

1. 工厂模式(Factory Pattern) 概念:工厂模式是一种创建对象的设计模式,它将对象的创建和使用分离,通过一个工厂类来负责创建对象。MyBatis 中的应用:MyBatis 使用 SqlSessionFactory 来创建 SqlSession 对…...

[PWNME 2025] PWN 复现

这种比赛得0也不容易,前边暖声还是能作的。 GOT 指针前溢出,可以溢出到GOT表,然后把后门写上就行 Einstein 这个拿到WP也没复现成,最后自己改了一下。 int __cdecl handle() {int offset; // [rsp8h] [rbp-38h] BYREFunsigne…...

【笔记】记一次easyExcel中注解ExcelProperty映射字段赋值无效问题

背景 最近在开发一个导入excel文件解析数据并入库的功能,导入过程中出现一个BUG:后端映射实体字段pType没有被赋值,但是导入的excel源文件是有数据。 定位问题 经过debug源码发现最终创建对象并赋值使用的是net.sf.cglib.beans.BeanMap类中…...

使用VSCode开发Django指南

使用VSCode开发Django指南 一、概述 Django 是一个高级 Python 框架,专为快速、安全和可扩展的 Web 开发而设计。Django 包含对 URL 路由、页面模板和数据处理的丰富支持。 本文将创建一个简单的 Django 应用,其中包含三个使用通用基本模板的页面。在此…...

MongoDB学习和应用(高效的非关系型数据库)

一丶 MongoDB简介 对于社交类软件的功能,我们需要对它的功能特点进行分析: 数据量会随着用户数增大而增大读多写少价值较低非好友看不到其动态信息地理位置的查询… 针对以上特点进行分析各大存储工具: mysql:关系型数据库&am…...

python/java环境配置

环境变量放一起 python: 1.首先下载Python Python下载地址:Download Python | Python.org downloads ---windows -- 64 2.安装Python 下面两个,然后自定义,全选 可以把前4个选上 3.环境配置 1)搜高级系统设置 2…...

Go 语言接口详解

Go 语言接口详解 核心概念 接口定义 在 Go 语言中,接口是一种抽象类型,它定义了一组方法的集合: // 定义接口 type Shape interface {Area() float64Perimeter() float64 } 接口实现 Go 接口的实现是隐式的: // 矩形结构体…...

Qwen3-Embedding-0.6B深度解析:多语言语义检索的轻量级利器

第一章 引言:语义表示的新时代挑战与Qwen3的破局之路 1.1 文本嵌入的核心价值与技术演进 在人工智能领域,文本嵌入技术如同连接自然语言与机器理解的“神经突触”——它将人类语言转化为计算机可计算的语义向量,支撑着搜索引擎、推荐系统、…...

使用van-uploader 的UI组件,结合vue2如何实现图片上传组件的封装

以下是基于 vant-ui&#xff08;适配 Vue2 版本 &#xff09;实现截图中照片上传预览、删除功能&#xff0c;并封装成可复用组件的完整代码&#xff0c;包含样式和逻辑实现&#xff0c;可直接在 Vue2 项目中使用&#xff1a; 1. 封装的图片上传组件 ImageUploader.vue <te…...

数据链路层的主要功能是什么

数据链路层&#xff08;OSI模型第2层&#xff09;的核心功能是在相邻网络节点&#xff08;如交换机、主机&#xff09;间提供可靠的数据帧传输服务&#xff0c;主要职责包括&#xff1a; &#x1f511; 核心功能详解&#xff1a; 帧封装与解封装 封装&#xff1a; 将网络层下发…...

【决胜公务员考试】求职OMG——见面课测验1

2025最新版&#xff01;&#xff01;&#xff01;6.8截至答题&#xff0c;大家注意呀&#xff01; 博主码字不易点个关注吧,祝期末顺利~~ 1.单选题(2分) 下列说法错误的是:&#xff08; B &#xff09; A.选调生属于公务员系统 B.公务员属于事业编 C.选调生有基层锻炼的要求 D…...

Maven 概述、安装、配置、仓库、私服详解

目录 1、Maven 概述 1.1 Maven 的定义 1.2 Maven 解决的问题 1.3 Maven 的核心特性与优势 2、Maven 安装 2.1 下载 Maven 2.2 安装配置 Maven 2.3 测试安装 2.4 修改 Maven 本地仓库的默认路径 3、Maven 配置 3.1 配置本地仓库 3.2 配置 JDK 3.3 IDEA 配置本地 Ma…...

深度学习水论文:mamba+图像增强

&#x1f9c0;当前视觉领域对高效长序列建模需求激增&#xff0c;对Mamba图像增强这方向的研究自然也逐渐火热。原因在于其高效长程建模&#xff0c;以及动态计算优势&#xff0c;在图像质量提升和细节恢复方面有难以替代的作用。 &#x1f9c0;因此短时间内&#xff0c;就有不…...