Mr123456

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 73
帖子: 18
在線: 16.6小時
蟲號: 1650518
注冊: 2012-02-28
專業(yè): 光學(xué)

[求助] 二維隨機行走運動

一粒子從原點隨機選擇（上，下，左，右）方向開始移動，經(jīng)過N步移動后（每步移動的距離為1），求距離原點的平均距離。（N=10000）（要完整的程序，可以在matlab上模擬）

回復(fù)此樓

» 猜你喜歡

327求調(diào)劑已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有3人回復(fù)
309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab生成隨機行走網(wǎng)格模擬材料生長已經(jīng)有26人回復(fù)

1樓 2012-02-28 18:24:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖置頂 ( 共有1個 )

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Mr123456(金幣+1): ★★★★★最佳答案很不錯，很強大。 2012-03-01 20:48:57
臭水溝(金幣+5, 仿真EPI+1): 專家辛苦了~~ 2012-03-01 20:50:01

CODE:

function randwalk

ra=0;

for i=1:10

    ra=ra+walk(i);

end

mr=ra/10%10次運行平均值

function rr=walk(x)

clc;close all;

z=0; % 設(shè)置初始數(shù)值，這里利用一個復(fù)數(shù)來表示位置

rand('state',x); % 初始化隨機數(shù)的狀態(tài)

ph=plot(z);hold on; % 繪圖

T=1; % 計算時間的參數(shù)

z1=1; % 位置的中間變量

ti=title('time = 0','Fontsize',14,'Fontname','Times new roman'); % 在圖題處實時顯示時間

while T<1001; % 進行循環(huán)1000次進行計算

    [dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 隨機生成一個位移量

    z1=z+[dx+i*dy]; % 更新位置數(shù)值

    plot(real([z,z1]),imag([z,z1])); % 連接t時刻和t+1時刻兩點繪制線段

    z=z1; % 更新z的數(shù)值

    set(ti,'String',['time = ',num2str(T)]); % 更新圖題處的時間數(shù)值

    T=T+1; % 時間值更新

    pause(0.01); % 暫停一下，顯示動畫效果

end

rr=abs(z)

function [dx,dy]=choose_D(p1,p2,p3,p4);

% p1=0.2; % 向左運動的幾率

% p2=0.3; % 向上運動的幾率

% p3=0.1; % 向右運動的幾率

% p4=0.4; % 向下運動的幾率

r=rand; % 產(chǎn)生一個隨機數(shù)

if r
    dx=-1; % x方向的隨機移動量

    dy=0; % y方向的隨機移動量

elseif r
    dx=0; % x方向的隨機移動量

    dy=-1; % y方向的隨機移動量

elseif r
    dx=1; % x方向的隨機移動量

    dy=0; % y方向的隨機移動量

else

    dx=0; % x方向的隨機移動量

    dy=1; % y方向的隨機移動量

end

贊一下(3人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2012-03-01 19:38:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

運行如圖所示

贊一下(1人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

3樓2012-03-01 19:58:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

Mr123456

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 73
帖子: 18
在線: 16.6小時
蟲號: 1650518
注冊: 2012-02-28
專業(yè): 光學(xué)

引用回帖:

樓: Originally posted by dbb627 at 2012-03-01 19:58:18:
運行如圖所示
75/2a/149791_1330603095_490.gif

p1,p2,p3,p4 前面的%去掉有什么影響？

贊一下

回復(fù)此樓

4樓2012-03-01 21:34:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

★ ★
xiegangmai(金幣+2, 專家考核): 謝謝參與！ 2012-03-01 23:12:27

[dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 隨機生成一個位移量
這句表明四個方向概率相等p1=p2=p3=p4=0.25
function [dx,dy]=choose_D(p1,p2,p3,p4);
下面的可以去掉
% p1=0.2; % 向左運動的幾率
% p2=0.3; % 向上運動的幾率
% p3=0.1; % 向右運動的幾率
% p4=0.4; % 向下運動的幾率

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5樓2012-03-01 22:05:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Mr123456

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 73
帖子: 18
在線: 16.6小時
蟲號: 1650518
注冊: 2012-02-28
專業(yè): 光學(xué)

引用回帖:

樓: Originally posted by dbb627 at 2012-03-01 22:05:48:
=choose_D(0.25,0.25,0.25,0.25); % 隨機生成一個位移量
這句表明四個方向概率相等p1=p2=p3=p4=0.25
function =choose_D(p1,p2,p3,p4);
下面的可以去掉
% p1=0.2; % 向左運動的幾率
% p2=0.3; % 向上運動的幾 ...

哦，明白了。很感謝。

回復(fù)此樓

6樓2012-03-02 13:50:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Mr123456

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 73
帖子: 18
在線: 16.6小時
蟲號: 1650518
注冊: 2012-02-28
專業(yè): 光學(xué)

引用回帖:

5樓: Originally posted by dbb627 at 2012-03-01 22:05:48:
=choose_D(0.25,0.25,0.25,0.25); % 隨機生成一個位移量
這句表明四個方向概率相等p1=p2=p3=p4=0.25
function =choose_D(p1,p2,p3,p4);
下面的可以去掉
% p1=0.2; % 向左運動的幾率
% p2=0.3; % 向上運動的幾 ...

想再畫條步數(shù)和平均距離的曲線，要怎么寫呢？

贊一下

回復(fù)此樓

7樓2012-03-23 16:38:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mengzqy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.5
帖子: 1
在線: 1.1小時
蟲號: 2634120
注冊: 2013-09-04

為什么這個程序在我的電腦沒法運行出來，而是沿著X軸運動

贊一下

回復(fù)此樓

8樓2013-09-04 21:53:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

胡北

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 2
在線: 1.6小時
蟲號: 3111630
注冊: 2014-04-03

引用回帖:

8樓: Originally posted by mengzqy at 2013-09-04 21:53:46
為什么這個程序在我的電腦沒法運行出來，而是沿著X軸運動

我也是

回復(fù)此樓

9樓2014-12-18 14:42:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

2361640267

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 1
在線: 29分鐘
蟲號: 4528003
注冊: 2016-03-21

引用回帖:

9樓: Originally posted by 胡北 at 2014-12-18 14:42:26
我也是...

我也是，是不是代碼哪里錯了。樓主能否再看下

贊一下

回復(fù)此樓

10樓2016-03-21 20:41:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +4	小小機靈蟲 2026-03-29	10/500	2026-03-29 18:20 by 無際的草原
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 295求調(diào)劑 +4	wei-5 2026-03-26	4/200	2026-03-28 23:20 by 小木蟲tim
[考研] 085701求調(diào)劑初試286分 +4	secret0328 2026-03-28	4/200	2026-03-28 21:09 by 15366876211
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護地球你我做?/a>
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[有機交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 085602化學(xué)工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 材料考研求調(diào)劑 +3	Dendel 2026-03-23	6/300	2026-03-26 17:51 by fmesaito
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 0854電子信息求調(diào)劑 324 +4	Promise-jyl 2026-03-23	4/200	2026-03-25 11:36 by Sugarlight
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯條 2026-03-23	5/250	2026-03-24 16:52 by 星空星月