1. 求一个MATLAB程序,通过蒙特卡洛方法模拟二维布朗运动的,谢谢
function [x,y,m,n]=br2(x0,xf,y0,yf,h)
x=x0:h:xf;
y=y0:h:yf;
a=randn(size(x));
b=randn(size(y));
m(1)=0;
n(1)=0;
for k=1:length(x)-1;
m(k+1)=m(k)+a(k);
n(k+1)=n(k)+b(k);
end;
再在命令窗口键入
x0=0;
xf=10;
h=0.01;
y0=0;
yf=10;
[x,y,m,n]=br2(x0,xf,y0,yf,h);
plot(m,n)
xlabel('m');
ylabel('n')
2. 如何用matlab模拟布朗运动
布朗运动的示例
[plain] view plain copy
%Brownian motion
clf;
n=20;
s=0.02;
x = rand(n,1)-0.5;
y = rand(n,1)-0.5;
h = plot(x,y,'.');
axis([-1 1 -1 1])
axis square
grid off
set(h,'EraseMode','xor','MarkerSize',18)
grid on;
title('Press Ctl-C to stop');
while 1
drawnow
x = x + s*randn(n,1);
y = y + s*randn(n,1);
set(h,'XData',x,'YData',y)
end
3. 关于布朗运动的一些问题 有大神会吗?求帮助啊啊!!!!
布朗运动是独立增量过程, 所以协方差, cov(Bs,Bt)=min(s,t), 可假设s>t证之.
Bt服从N(0,t), 积分即得原点反射的期望方差.
Bt^2依然是独立增量过程, 并且E(Bt^2) = t. 对于此的计算, 一者可以利用独立增量与正态分布随机变量四阶矩的相关计算完成, 二者可利用(对正态分布可用)
E(X1X2X3X4)=E(X1X2)E(X3X4)+E(X1X3)E(X2X4)+E(X1X4)E(X2X3)
依然是独立增量过程. 利用此性质可简化积分
4. 布朗运动怎么在matlab上模拟
clc;
clear all;
t=1;
x=100;
y=100;
vx=17;vy=0;
for k=1:40
p=2*pi*rand(1,1) ;
vx=17*cos(p);
vy=17*sin(p);
x=x+vx*t;
y=y+vy*t;
line('ydata',y,'xdata',x,'Color' ,[1 0 0], 'Marker' ,'.' , 'MarkerSize' ,12, 'EraseMode' , 'non');
plot(x,y);
axis([0 200 0 200])
if(x>200||y>200||x<0||y<0)
break
end
hold on
pause(0.4);
end
望采纳
5. matlab多粒子布朗运动记录坐标问题。
你编程生成的图太恶心了,不过符合蒙卡的思想,你可以从布朗运动的本质入手,定义粒子后让他的运动方向是随机的,如一维中定义前进后退(-1)^n让n成为那个随机数,然后拓展到空间,三维方向上去。。。很简单的
6. 怎么用matlab用标准布朗运动做利率模型
可以模拟多个质点的布朗运动动态模型,可以更好地观察布朗运动!
7. 如何用matlab编写蒙特卡洛模拟计算VaR的程序
这种资源建议去官网论坛或者一些其他比较知名的大论坛里咨询
如果曾经有人在网上发布,或有网站下载,一般都会被搜索引擎收录,如果搜不到,你可以找一下相关的论坛
最好是那种人气比较高的论坛,注册会员,发帖求助,会有高手帮你的。
8. matlab 模拟质点运动的问题
自由运动,速度不是恒定的么?
ps,万一正撞了,怎么样“旋转90度分开”呢?这个描述不物理啊。