本文介绍: 1D和2D布朗运动matlab
布朗运动是一种随机现象,下面的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
原文地址:https://blog.csdn.net/y15520833229/article/details/134788003
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如若转载,请注明出处:http://www.7code.cn/show_50877.html
如若内容造成侵权/违法违规/事实不符,请联系代码007邮箱:suwngjj01@126.com进行投诉反馈,一经查实,立即删除!
声明:本站所有文章,如无特殊说明或标注,均为本站原创发布。任何个人或组织,在未征得本站同意时,禁止复制、盗用、采集、发布本站内容到任何网站、书籍等各类媒体平台。如若本站内容侵犯了原著者的合法权益,可联系我们进行处理。