24小时热门版块排行榜

返回列表

zhouyuxing

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 2
在线: 1.9小时
虫号: 2368982
注册: 2013-03-22
专业: 环境工程

[求助] matlab 粒子布朗运动逸出区域问题：

以下是1000个例子在三维空间中的布朗运动，运动过程中有粒子溢出坐标轴区域，如何使运动到边界的粒子向内部继续做布朗运动？

n=1000;%指布朗运动的点数

s=0.02;%指温度或速率

%产生n个随机点(x,y,z)，处于-0.5~0.5之间

x=2*rand(n,1)-1;

y=2*rand(n,1)-1;

z=2*rand(n,1)-1;

h=plot3(x,y,z,'.');

xmin=-1;

xmax=1;

ymin=-1;

ymax=1;

zmin=-1;

zmax=1;

axis([xmin xmax ymin ymax zmin zmax]);

axis square

grid on

set(h, 'MarkerSize',5);%设置擦除模式

%循环5000次，产生动画效果

for i=linspace(1,10,5000)

drawnow

x=x+s*randn(n,1);

y=y+s*randn(n,1);

z=z+s*randn(n,1);

set(h,'XData',x,'YData',y,'ZData',z);

end

1.png

回复此楼

» 猜你喜欢

之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
博士申请都是内定的吗？已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有12人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【原创】关于MATLAB 描述分子布朗运动等行为的过程已经有59人回复

1楼2013-05-04 11:02:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

感谢参与，应助指数 +1

x=x+s*randn(n,1);
y=y+s*randn(n,1);
z=z+s*randn(n,1);
------------------------
这不是更新位置吗?判断一下上下界,如果是碰到边界,让速度方向改变就行了

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2013-05-06 17:41:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhouyuxing

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 2
在线: 1.9小时
虫号: 2368982
注册: 2013-03-22
专业: 环境工程

引用回帖:

2楼: Originally posted by libralibra at 2013-05-06 17:41:35
x=x+s*randn(n,1);
y=y+s*randn(n,1);
z=z+s*randn(n,1);
------------------------
这不是更新位置吗?判断一下上下界,如果是碰到边界,让速度方向改变就行了

自己弄明白了，还是感谢你。

赞一下

回复此楼

3楼2013-07-05 12:26:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖