《随机过程》布朗运动理论中的两个反常问题

全世界只有3.14 % 的人关注了

爆炸吧知识

1827 年,英国植物学家布朗(Brown)用显微镜观察悬浮在液体中的花粉微粒时,发现花粉微粒总是在做无规则运动。后来人们发现,这是一种广泛存在于自然界、工程技术和社会经济等领域中的随机现象,如空气污染扩散、陀螺随机游走和股票价格波动等,因此这类随机运动现象被称为“布朗运动”。

1905年,爱因斯坦首先使用概率方法对一维布朗运动进行了定量研究,并从热分子运动扩散方程推导出了大量布朗粒子在任意时刻的位置服从正态分布的性质。

1908年,法国物理学家朗之万(Langevin)依据牛顿运动定律,建立了单个布朗粒子的动力学方程,并对大量布朗粒子的位移求平均,得出了与爱因斯坦布朗运动理论相同的统计规律,但是并未给出揭示单个布朗粒子运动规律的运动学方程。

1923年,维纳(Wiener)将单个布朗粒子的位移假设为随机变量,建立了描述单个布朗粒子位移与时间之间数量关系的运动学数学模型,得出了布朗粒子运动路径的一系列数学性质,从而为现代随机过程理论的建立奠定了概念、方法和理论基础,因此布朗运动在《随机过程》教科书中也被称为“维纳过程”。

一维布朗运动物理现象及统计规律

假设一个布朗粒子从x轴的原点出发,则布朗粒子在t时刻的位移x(t)是时间t的函数,x(t)的函数图像(质点位移轨迹)见图1(b)。

图1一维布朗粒子位移曲线

图2(a)和图2(b)分别为MEMS微机械陀螺的随机游走误差曲线和特斯拉(Tesla)股票价格收盘价曲线,均与图1(b)所示的布朗粒子位移曲线相似。

 

图2陀螺随机游走与股票价格曲线

美国爱达荷大学的哈蒙(Harmon)教授在《Phylogenetic ComparativeMethods》一书中描述了同时观察100个布朗粒子位移的仿真实验结果。哈蒙教授分别做了时间周期为10s、50s和100s时的三个实验,得到图3中A、B、C所示的三组实验结果曲线。

图3布朗粒子位移曲线(100个)

 

哈蒙教授对100个布朗粒子在t=10s、50s和100s时的位置进行了统计分析,对应的频率直方图分别列于实验结果的右边。

从频率直方图可以看出,直方图轮廓与正态分布的钟形曲线相似,呈现中间高,两端逐渐下降且完全对称的状态。

爱因斯坦在1905年根据扩散方程推导出了大量布朗粒子在t时刻空间位置的概率分布函数

式中D为扩散系数。

显然,所有布朗粒子在t时刻的位置服从数学期望为零、方差为2Dt的正态分布。图4给出了t=0.1秒、0.2秒、0.5秒和2秒时的正态分布曲线。

图4 布朗运动正态分布曲线

维纳过程数学性质

维纳将图1所示的单个布朗粒子作为研究对象,并假设布朗粒子在t时刻的位移X(t)为随机变量,得出了布朗粒子位移X(t)的两个重要性质:

(1)布朗粒子位移X(t)服从(0,σ2t)正态分布;

(2)布朗粒子路径处处不可导(瞬时速度无穷大或不存在)。

维纳过程反常问题

1、正态分布性质与经验事实不符

如果布朗粒子在t时刻的位移X(t)服从正态分布,则图1(b)所示的布朗粒子位移曲线应具有如下正态分布的两个特点:

(1)对称性。在每一时刻t,绝对值相等的正、负位移出现的次数大致相等。

(2)集中性。布朗粒子在0点附近出现的次数最多。

但是从图1(b)所示的布朗粒子位移曲线可以看出,布朗粒子随时间远离原点,其位移曲线x(t)既不符合正态分布的对称性,也不符合正态分布的集中性

事实上,服从(0,σ2t)正态分布的质点位移曲线如图5所示,是均值为零、方差与时间成正比的高斯噪声。高斯噪声与图1(b)所示的布朗粒子位移曲线有天壤之别,因此,实际的布朗粒子位移X(t)不服从正态分布。

图5 服从正态分布的质点位移曲线

维纳显然误解了爱因斯坦“布朗运动服从正态分布”的结论,爱因斯坦布朗运动理论描述的是大量布朗粒子在t时刻的位置服从正态分布(图2),而不是指一个布朗粒子的位移服从正态分布。

维纳将单个布朗粒子的位移假设为随机变量,因而只能用刻画大量布朗粒子集体行为的正态分布,来描述单个布朗粒子的个体行为,从而得出了与事实完全不符的错误结论。

2、处处不可导性质与物理学理论和实验不符

根据爱因斯坦“同一个布朗粒子在不同时间间隔中的运动相互独立”假设,可知单个布朗粒子的瞬时速度v(t)在不同时刻互不相关,因此单个布朗粒子瞬时速度v(t)的自相关函数可表示为

Rv(τ)=N0δ(τ)

式中τ为时间间隔,N0为正实常数,δ(τ)为单位冲击函数。

根据维纳-辛钦定理,平稳随机过程的功率谱密度是其自相关函数的傅立叶变换,可得单个布朗粒子瞬时速度v(t)的功率谱密度

Sv(f)=N0

即v(t)的功率谱密度在整个频率轴上均匀分布,表明v(t)为平均功率为N0的白噪声。

2010年,美国得克萨斯大学的李统藏成功地用激光光镊技术首次实验测量到了悬浮布朗粒子的瞬时速度(图6),实验结果证明了布朗粒子的瞬时速度波形为白噪声,表明布朗运动的导数(瞬时速度)不仅存在,而且可观测。

图6布朗粒子瞬时速度测量实验及结果

 

李统藏的研究论文“Measurement of the Instantaneous Velocity of a Brownian Particle”在《科学(Science)》杂志上发表后,在全球引起了极大的轰动。《科学》杂志专门为李统藏的论文配发了录音采访,《自然》杂志也迅速报道了该实验。

李统藏的布朗粒子瞬时速度测量实验被《科学》杂志推荐为大学及高中教学内容,美国明尼苏达大学等学校的相关课程已经将该实验作为教学内容。

维纳过程对数理金融学的影响

维纳将单个布朗粒子位移抽象为随机变量的研究方法及结论,给数理金融学及华尔街金融市场带来了灾难性的影响。

股票价格(对数)随时间的变化过程与一个布朗粒子的运动过程在数学形式上完全相同,因此数理金融学也将股票价格假设为随机变量,并根据“布朗粒子位移X(t)服从(0,σ2t)正态分布”的性质,认为股票价格也服从正态分布,并用未来T时刻正态分布的标准差来描述未来T时刻的股票价格波动程度,因此建立的股票价格模型必然与事实不符,无法正确描述并预测股票价格的波动趋势。

数理金融学基于维纳过程建立的股票价格模型和BS期权定价公式在大规模应用于金融市场时,竟成为直接导致 1987、1997 和 2007 年三次重大金融危机的罪魁祸首。麦肯齐(Mackenzie)在《无言的宇宙(隐藏在24个数学公式背后的故事)》书中总结到:1987、1998和2007年三次重大金融危机的一般共性是股票价格布朗运动模型未能预期市场波动。

被誉为“中国金融数学开创者”、获得2020未来科学大奖“数学与计算机科学奖”的彭实戈院士,在《中国基础研究发展报告》第二章中国数学前沿进展中也明确指出:BS 期权定价理论是造成以前历次重大金融危机的关键性原因。

畅销书《黑天鹅》作者塔勒布(Taleb)在《金融时报》上发表了题为“破坏市场的伪科学”专栏文章,对数理金融学进行了严厉的批判。塔勒布在文章中指出:人们从一次又一次的金融危机中得出了“数理金融学的有效性与占星术一样不靠谱”和“数理金融学通过创造风险来危害金融系统”的结论。塔勒布痛斥数理金融学是破坏市场的伪科学,数理金融学理论获得诺贝尔奖不仅是对科学的侮辱,数理金融学一直使金融体系面临崩溃的风险。

结论

“与实际结合,问题驱动”是随机过程等应用数学学科发展的不竭动力和重要特征。随机过程布朗运动理论出现上述反常现象,表示现有《随机过程》教科书内容与经验事实不符,无法用来描述并解决实际问题,意味着《随机过程》教科书内容将面临重大范式变革,新的布朗运动概念、方法及理论将会推翻并替代现有教科书内容,从而把人类对布朗运动的认识提高到一个崭新的水平。

写在最后

大家有什么想法,欢迎发表到留言区,或者一键转发分享。

同时,也希望大家关注一下超模君精心推出的年度知识周边同理可得系列文化衫!!

穿在身上的

不仅仅是信仰!

原价158元

4月1日至31日下单
+1元换购杯垫一套(3枚)

可特惠价128元!!!

猛戳下单《同理可得·系列文化衫》

作者简介

高宏,毕业于清华大学精密仪器系,分别获工学学士、硕士和博士学位,留校任教从事测试信号分析与处理的教学与科研工作,现任紫光股份有限公司CTO,北京市科协委员。

点这里????关注我,记得标星~

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

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

相关文章

linux之用2张图片描述vim常见命令

对了,使得光标跳转到最后一行是这个命令 G

读《好好学习:个人知识管理精进指南》

关于学习的文章之前写过两篇:《掌握好的学习方法,让你在职场更有竞争力》《程序员是终身学习的职业,应该怎么学习?》我们都是终身学习者,我深知学习的重要性,所以每隔一段时间,有些新的心得和想…

跨域解决方案大全

什么是跨域 注:本文完整示例地址先来说一个概念就是同源,同源指的是协议,端口,域名全部相同。 同源策略(Same origin policy)是一种约定,它是浏览器最核心也最基本的安全功能,如果缺…

2013年下半年信息系统项目管理师考试试卷(回忆版)

2013年下半年信息系统项目管理师上午试卷(网友回忆版)信息系统的生命周期可以分为立项、开发、运维和消亡4个阶段,应在信息系统建设的(1)考虑系统消亡的条件和时机。(1)A.初期B&…

酱油和gbt酱油哪个好_酱油可不是越贵越好?看清瓶身上的5个字,教你1分钟买到好酱油...

今天妈妈做饭说酱油用完了,让我去打酱油的地方打5毛钱的酱油,我拿起塑料壶去打酱油,闻起来香香的,一个推自行车的过来,我问他这酱油是勾兑酱油吗?推自行车的小哥,看了我一眼,问我是监…

让 WPF 的 RadioButton 支持再次点击取消选中的功能

让 WPF 的 RadioButton 支持再次点击取消选中的功能目录让 WPF 的 RadioButton 支持再次点击取消选中的功能零、前言一、方法一:后台直接处理二、方法二:提取为自定义控件(用户控件)三、方法三:附加行为法独立观察员 2…

java数组转换成string_java面试复习重点:类的管理及常用工具,教你抓住面试重点

java复习: 类的管理及常用工具类包写在程序文件的第一行一个Java 源文件中只能声明一个包,且声明语句只能作为源文件的第一条指令导入类能导入非public类,但是不能用因为在其他包缺省的权限用不了package Testp;import Testpackage.*;public …

ubuntu 新建的用户 table 无法补全命令 解决办法

为什么80%的码农都做不了架构师?>>> vi /etc/passwd 用adduser命令新增了用户之后,发现在该新建用户下的命令终端,使用方向键无法调出历史命令,同时tab键也无法补全输入命令。 找到 你新增的用户xxx 修改/bin/bash…

【Unity3D基础】让物体动起来②--UGUI鼠标点击逐帧移动

背景 上一篇通过鼠标移动的代码很简单,所以看的人也不多,但是还是要感谢“武装三藏”在博客园给出的评论和支持,希望他也能看到第二篇,其实可以很简单,而且是精灵自控制,关键是代码少是我喜欢的方式&#x…

一个有趣的Go项目,3D界面管理k8s集群,真好玩!

大家好,我是小碗汤,今天分享一个用Golang开发,很好玩的工具KubeCraftAdmin[1]:用Minecraft方式管理k8s的工具,感兴趣的兄弟不妨玩一玩。文末有视频,供您鉴赏~Minecraft:我的世界,是微…

java 数组拼接_打印Java数组最优雅的方式是什么?这波操作闪瞎我

在 Java 中,数组虽然是一个对象,但并未明确的定义这样一个类,因此也就没有覆盖 toString() 方法的机会。如果尝试直接打印数组的话,输出的结果并不是我们预期的结果。那有没有一些简单可行的方式呢?如果大家也被这个问…

AmazeUI基本样式

2019独角兽企业重金招聘Python工程师标准>>> AmazeUI是一个轻量级、Mobile first的前端框架,基于开源社区流行的前端框架编写。 Normalize AmazeUI使用了normalize.css,但做了些调整:html添加了-webkit-font-smoothing:antialiase…

Oracle基础中的基础视频讲座录像(西安)供免费下载

下载地址播放器也在上述目录中。记得那年园博会,培训中间有一天参观:转载于:https://blog.51cto.com/botang/1323099

她13岁自己造飞机,17岁进麻省理工,3篇黑洞论文被霍金引用......

全世界只有3.14 % 的人关注了爆炸吧知识13岁的时候,你在干嘛?我想,绝大多数人都没有萨布丽娜这么硬核——独自组装飞机。之后,16岁的萨布丽娜完成了生涯首飞,历史第一年轻。17岁时,她考上麻省理工&#xff…

Windows 10 2022 年更新来了!

面向 Release Preview 频道的 Windows 10 预览体验成员,微软现已发布 Windows 10 Build 19044.1499。主要修复1.微软修复了阻止某些环绕声音频在 Microsoft Edge 中播放的问题。2.微软修复了使用中文输入法时,一些应用停止工作的意外错误。3.微软修复了在…

当代成年人的生活状态......

1 小喵咪举起了它的狙击枪▼2 简直有毒...▼3 哈哈哈哈哈▼4 电焊既视感▼5 请问需要小猫咪吗?▼6 哈哈哈哈哈▼7 一物降一物▼7 数学能有多好玩?▼你点的每个赞,我都认真当成了喜欢

c++ why can't class template hide its implementation in cpp file?

类似的问题还有: why cant class template use Handle Class Pattern to hide its implementation? || why there are linker problems (undefined reference) to my class template? 我出现问题的源码(见main.cpp,Stack.h,Stack.cpp)&…

C# using static 声明

许多实际的扩展可以通过扩展方法来实现,并非所有实际的扩展都有可以扩展的类型。对于某些场景,简单的静态方法比较适合。为了更容易调用这些方法,可以使用 using static 声明除去类名。例如,如果打开了 System.Console using stat…

PHP性能追踪及分析工具xhprof的安装与使用

PHP性能追踪及分析工具xhprof的安装与使用 对于本地开发环境来说,进行性能分析xdebug是够用了,但如果是线上环境的话,xdebug消耗较大,配置也不够灵活,因此线上环境建议使用xhprof进行PHP性能追踪及分析。 我们今天就简…