应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab 粒子布朗运动逸出区域问题：
81/1返回列表
查看: 2448  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhouyuxing

新虫 (初入文坛)

[求助] 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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-05-04 11:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与，应助指数 +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

新虫 (初入文坛)
引用回帖:
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的回帖

牧马人man

新虫 (初入文坛)

【答案】应助回帖

读博真的很累很那毕业吗？本人研二正在就业与考博之间犹豫，听说读博很难毕业压力很大，现在有点胆却。正在读博的大虾们能给点指点吗？说说读博的真实生活与学习，谢谢啦
赞一下 回复此楼
4楼2013-07-05 20:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿萌宝贝1987

银虫 (正式写手)
引用回帖:
3楼: Originally posted by zhouyuxing at 2013-07-05 12:26:46
自己弄明白了，还是感谢你。...

楼主 你好 你是在做颗粒成长模型吗 我也是研究相关领域的 可以他探讨一下哈
赞一下 回复此楼
5楼2013-08-17 23:31:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿萌宝贝1987

银虫 (正式写手)
引用回帖:
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);
------------------------
这不是更新位置吗?判断一下上下界,如果是碰到边界,让速度方向改变就行了

你好 关于这个例子 我想请问下 判断上下界的程序语言应该怎么写呢？以什么标准判断粒子是否碰到边界呢？是粒子和边界之间的距离吗？多谢回复，我也再做这样的模型
赞一下(1人) 回复此楼
6楼2013-08-17 23:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

引用回帖:
6楼: Originally posted by 阿萌宝贝1987 at 2013-08-17 23:37:36
你好 关于这个例子 我想请问下 判断上下界的程序语言应该怎么写呢？以什么标准判断粒子是否碰到边界呢？是粒子和边界之间的距离吗？多谢回复，我也再做这样的模型...

axis([xmin xmax ymin ymax zmin zmax]);
这是三维的边界,
需要判断粒子坐标的x,y,z在这个区域内,否则按照碰撞处理,重新计算位置
赞一下 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
7楼2013-08-19 22:18:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zymasxcm

新虫 (小有名气)
引用回帖:
3楼: Originally posted by zhouyuxing at 2013-07-05 12:26:46
自己弄明白了，还是感谢你。...

您好，我想知道，你采用的是周期性边界条件吗
赞一下 回复此楼
8楼2014-10-14 10:52:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhouyuxing 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护，明天很可能放榜！ +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭