 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 273 | 回復(fù): 0 | ||
[求助]
在Matlab上模擬DLA的吸附，請大神幫忙修改下代碼，謝謝
|
|
% 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'); % 顯示時間 T=0; % 記錄時間的參數(shù) for k=1:30000; % 循環(huán)計算 xt=4; % 產(chǎn)生粒子的位置 yt=N/2; % 產(chǎn)生粒子的位置 Ss=0; % 控制下面循環(huán)是否終止的參數(shù) while Ss<1; % 計算當(dāng)前粒子的吸附過程 [dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 計算粒子位移量 xt=xt+dx; % 計算下一時刻粒子的位置 yt=yt+dy; % 計算下一時刻粒子的位置 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)前運動的粒子 Ss=2; % 把Ss賦值為2,從而可以退出循環(huán) end end set(Ii,'CData',S); % 顯示當(dāng)前狀態(tài)矩陣 T=T+1; % 累加時間參數(shù) set(ti,'String',['time = ',num2str(T),', N=',num2str(sum(sum(1-S))-[N/4+1])]); % 更新時間參數(shù) pause(0.01); % 暫停一下，顯示動畫效果 end |
回復(fù)此樓» 猜你喜歡
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有281人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有79人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
點擊這里搜索更多相關(guān)資源
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 痞子ppt 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
|
[考研] 271求調(diào)劑 +4 | 生如夏花… 2026-03-22 | 4/200 |
|
|
|
[考研] 上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +6 | Stella_Yao 2026-03-20 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 | ￣^￣゜汗 2026-03-19 | 9/450 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|