 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 570 | 回復(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ù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有203人回復(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ù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，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ā)表 | |
|---|---|---|---|---|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過(guò)A區(qū)國(guó)家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
|
[考研] 277跪求調(diào)劑 +3 | 1915668 2026-03-27 | 4/200 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+5
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 349求調(diào)劑 +5 | 杰斯塔里斯 2026-03-21 | 5/250 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +3 | 哇呼哼呼哼 2026-03-26 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +3 | Li李魚 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|