1D和2D布朗运动matlab

最新推荐文章于 2024-11-11 19:42:18 发布

自我o救赎

最新推荐文章于 2024-11-11 19:42:18 发布

阅读量464

点赞数 1

分类专栏： Matlab语言专栏文章标签： matlab 布朗运动

本文链接：https://blog.csdn.net/y15520833229/article/details/134788003

版权

Matlab语言专栏同时被 2 个专栏收录

15 篇文章 0 订阅

订阅专栏

随机微积分专栏1

10 篇文章 3 订阅

订阅专栏

布朗运动是一种随机现象，下面的M函数brwnm2.m给出了二维Brown运动，其中[t0,tf]是时间区间，h是采样步长，w(t)，z(t)是布朗运动。

function [t,w,z]=brwnm2(t0,tf,h)
t=t0:h:tf;
x=randn(size(t))*sqrt(h);
y=randn(size(t))*sqrt(h);
w(1)=0;
z(1)=0;
  for k=1:length(t)-1
     w(k+1)=w(k)+x(k);
     z(k+1)=z(k)+y(k);
  end
  plot(w,z);

 

在command window下运行 brwnm2(0,10,0.01)可得到如下图像 作者：ocean_tide https://www.bilibili.com/read/cv27155811/ 出处：bilibili

在这里插入图片描述

下面的M函数brwnm3.m给出了三维Brown运动。
function [t,w,z,i]=brwnm3(t0,tf,h)
t=t0:h:tf;
x=randn(size(t))*sqrt(h);
y=randn(size(t))*sqrt(h);
s=randn(size(t))*sqrt(h);
w(1)=0;
z(1)=0;
i(1)=0;
  for k=1:length(t)-1
     w(k+1)=w(k)+x(k);
     z(k+1)=z(k)+y(k);
     i(k+1)=i(k)+s(k);
  end
  plot3(w,z,i);grid on ;box on;

在command window下运行 brwnm3(0,10,0.01)可得到如下图像 作者：ocean_tide https://www.bilibili.com/read/cv27155811/ 出处：bilibili