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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab 粒子布朗运动逸出区域问题：
51/1返回列表
查看: 2523  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

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的回帖

牧马人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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭