 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 565 | 回復(fù): 1 | ||
[求助]
在Matlab上模擬DLA的吸附，請(qǐng)大神幫忙修改下代碼，謝謝
|
|
% clc;clear;close all; rand('state',0); set(gcf,'DoubleBuffer','on'); N=256; % 生成狀態(tài)矩陣大小的控制參數(shù) S=ones(N); % 生成狀態(tài)矩陣S % S(N/2,N/4:N*3/4)=0; % 設(shè)置狀態(tài)矩陣的初始值 S(N/2:N*3/4,N/2)=0; % 設(shè)置狀態(tài)矩陣的初始值 Ii=imshow(S); % 顯示狀態(tài)矩陣 ti=title(['time = 0',', N=0'],'Fontsize',14,'Fontname','Times New Roman'); % 顯示時(shí)間 T=0; % 記錄時(shí)間的參數(shù) for k=1:30000; % 循環(huán)計(jì)算 xt=4; % 產(chǎn)生粒子的位置 yt=N/2; % 產(chǎn)生粒子的位置 Ss=0; % 控制下面循環(huán)是否終止的參數(shù) while Ss<1; % 計(jì)算當(dāng)前粒子的吸附過(guò)程 [dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 計(jì)算粒子位移量 xt=xt+dx; % 計(jì)算下一時(shí)刻粒子的位置 yt=yt+dy; % 計(jì)算下一時(shí)刻粒子的位置 if xt<2.5|xt>N-1.5||yt<2.5|yt>N-1.5; % 判斷粒子是否逸出邊界 Ss=2; % 把Ss賦值為2,從而可以退出循環(huán) end if [S(xt+1,yt)+S(xt-1,yt)+S(xt,yt+1)+S(xt,yt-1)]<3.5; % 判斷近鄰是否有粒子 S(xt,yt)=0; % 固定當(dāng)前運(yùn)動(dòng)的粒子 Ss=2; % 把Ss賦值為2,從而可以退出循環(huán) end end set(Ii,'CData',S); % 顯示當(dāng)前狀態(tài)矩陣 T=T+1; % 累加時(shí)間參數(shù) set(ti,'String',['time = ',num2str(T),', N=',num2str(sum(sum(1-S))-[N/4+1])]); % 更新時(shí)間參數(shù) pause(0.01); % 暫停一下，顯示動(dòng)畫效果 end |
回復(fù)此樓» 猜你喜歡
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有80人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有120人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 在matlab上運(yùn)行，大神幫忙調(diào)試下，是搞分形模擬的DLA
已經(jīng)有10人回復(fù)
2樓2013-10-11 08:35:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 痞子ppt 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑
+6
|
18419759900 2026-03-25 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +5 | 曼111 2026-03-24 | 6/300 |
|
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 081700 調(diào)劑 267分 +10 | 迷人的哈哈 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +6 | NIFFFfff 2026-03-20 | 6/300 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研]
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|