版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

愛(ài)笑的眼睛na

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1118.1
帖子: 17
在線: 80.7小時(shí)
蟲(chóng)號(hào): 2027185
注冊(cè): 2012-09-24
性別: MM
專(zhuān)業(yè): 煤炭地下開(kāi)采

[求助] 隨機(jī)生長(zhǎng)法模擬多孔介質(zhì)的matlab二維代碼改寫(xiě)三維，結(jié)果出來(lái)還是二維已有2人參與

本人根據(jù)http://www.gaoyang168.com/html/201111/3749191.html隨機(jī)生長(zhǎng)法模擬多孔介質(zhì)的二維代碼編寫(xiě)了一個(gè)三維代碼，但結(jié)果出來(lái)還是二維的，請(qǐng)求各位大神查一下錯(cuò)誤在哪兒，并給出解決的方法，謝謝啦！
代碼如下：
%======設(shè)定初值=============
clc;clear;
rand('state',0);
set(gcf,'DoubleBuffer','on');%消除振動(dòng)
Cd=0.01;%生長(zhǎng)核分布概率數(shù)
Vr=0.01;%生長(zhǎng)相的體積分?jǐn)?shù)
global kn
kn=500;
p=[1/8 1/8 1/8 1/8 1/4 1/4 1/4 1/4];%方向增長(zhǎng)概率數(shù)
A=rand(kn,kn);%初始網(wǎng)格kn×kn
A(find(A<=Cd))=0;%黑色，生長(zhǎng)核
A(find(A>Cd))=1;%白色
m=0;%迭代步驟
N=kn*kn*kn;%格點(diǎn)數(shù)
[x,y,z]=find(A==0);%生長(zhǎng)核腳標(biāo)
[y1,z] = ind2sub([kn,kn],y);
Ixy1z = sub2ind([kn,kn,kn],x,y1,z);
gn=length(Ixy1z);%生長(zhǎng)核計(jì)數(shù)
Vs=gn/N %體積分?jǐn)?shù)
%===========================
while Vs<0.001
   xW=x-1;xE=x+1;
xW(xW==0)=1;% 對(duì)近鄰進(jìn)行邊界條件處理
xE(xE==kn+1)=kn;% 對(duì)近鄰進(jìn)行邊界條件處理
y1N=y1-1;y1S=y1+1;
y1N(y1N==0)=1;% 對(duì)近鄰進(jìn)行邊界條件處理
y1S(y1S==kn+1)=kn;% 對(duì)近鄰進(jìn)行邊界條件處理
zB=z-1;zT=z+1;
zB(zB==0)=1;% 對(duì)近鄰進(jìn)行邊界條件處理
zT(zT==kn+1)=kn;% 對(duì)近鄰進(jìn)行邊界條件處理

InE=sub2ind([kn,kn],xE,z); % 右方轉(zhuǎn)化腳標(biāo)為索引
InB=sub2ind([kn,kn],x,zB); % 下方轉(zhuǎn)化腳標(biāo)為索引
InW=sub2ind([kn,kn],xW,z); % 左方轉(zhuǎn)化腳標(biāo)為索引
InT=sub2ind([kn,kn],x,zT); % 上方轉(zhuǎn)化腳標(biāo)為索引
   InS=sub2ind([kn,kn],x,y1N); % 下方轉(zhuǎn)化腳標(biāo)為索引
InN=sub2ind([kn,kn],x,y1S); % 上方轉(zhuǎn)化腳標(biāo)為索引

   InET=sub2ind([kn,kn],xE,zT); % 右上方轉(zhuǎn)化腳標(biāo)為索引
InEB=sub2ind([kn,kn],xE,zB); % 右下方轉(zhuǎn)化腳標(biāo)為索引
InWB=sub2ind([kn,kn],xW,zB); % 左下方轉(zhuǎn)化腳標(biāo)為索引
InWT=sub2ind([kn,kn],xW,zT); % 左上方轉(zhuǎn)化腳標(biāo)為索引
  InEN=sub2ind([kn,kn],xE,y1N); % 右上方轉(zhuǎn)化腳標(biāo)為索引
InES=sub2ind([kn,kn],xE,y1S); % 右下方轉(zhuǎn)化腳標(biāo)為索引
InWS=sub2ind([kn,kn],xW,y1S); % 左下方轉(zhuǎn)化腳標(biāo)為索引
InWN=sub2ind([kn,kn],xW,y1N); % 左上方轉(zhuǎn)化腳標(biāo)為索引
InTS=sub2ind([kn,kn],y1S,zT); % 右上方轉(zhuǎn)化腳標(biāo)為索引
InBS=sub2ind([kn,kn],y1S,zB); % 右下方轉(zhuǎn)化腳標(biāo)為索引
InBN=sub2ind([kn,kn],y1N,zB); % 左下方轉(zhuǎn)化腳標(biāo)為索引
InTN=sub2ind([kn,kn],y1N,zT); % 左上方轉(zhuǎn)化腳標(biāo)為索引

%==========做圖========
CInE=setdiff(InE,Ixy1z); En=rand(1,length(CInE));[i,j,k]=find(En<=p(5));pE=sub2ind([1,length(CInE)],i,j,k);
A(CInE(pE))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInT= setdiff(InT,Ixy1z);Tn=rand(1,length(CInT));[i,j,k]=find(Tn<=p(6));pT=sub2ind([1,length(CInT)],i,j,k);
A(CInT(pT))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInW=setdiff(InW,Ixy1z);Wn=rand(1,length(CInW));[i,j,k]=find(Wn<=p(7));pW=sub2ind([1,length(CInW)],i,j,k);
   A(CInW(pW))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInB=setdiff(InB,Ixy1z);Bn=rand(1,length(CInB));[i,j,k]=find(Bn<=p(8));pB=sub2ind([1,length(CInB)],i,j,k);
   A(CInB(pB))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInEB=setdiff(InEB,Ixy1z);EBn=rand(1,length(CInEB));[i,j,k]=find(EBn<=p(8));pEB=sub2ind([1,length(CInEB)],i,j,k);
   A(CInEB(pEB))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInET=setdiff(InET,Ixy1z);ETn=rand(1,length(CInET));[i,j,k]=find(ETn<=p(8));pET=sub2ind([1,length(CInET)],i,j,k);
   A(CInET(pET))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInWT=setdiff(InWT,Ixy1z);WTn=rand(1,length(CInWT));[i,j,k]=find(WTn<=p(8));pWT=sub2ind([1,length(CInWT)],i,j,k);
   A(CInWT(pWT))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInWB=setdiff(InWB,Ixy1z);WBn=rand(1,length(CInWB));[i,j,k]=find(WBn<=p(8));pWB=sub2ind([1,length(CInWB)],i,j,k);
   A(CInWB(pWB))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInS= setdiff(InS,Ixy1z);Sn=rand(1,length(CInS));[i,j,k]=find(Sn<=p(5));pS=sub2ind([1,length(CInS)],i,j,k);
A(CInS(pS))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInN=setdiff(InN,Ixy1z);Nn=rand(1,length(CInN));[i,j,k]=find(Nn<=p(8));pN=sub2ind([1,length(CInN)],i,j,k);
   A(CInN(pN))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInEN=setdiff(InEN,Ixy1z);ENn=rand(1,length(CInEN));[i,j,k]=find(ENn<=p(8));pEN=sub2ind([1,length(CInEN)],i,j,k);
   A(CInEN(pEN))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInES=setdiff(InES,Ixy1z);ESn=rand(1,length(CInES));[i,j,k]=find(ESn<=p(8));pES=sub2ind([1,length(CInES)],i,j,k);
   A(CInES(pES))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInWS=setdiff(InWS,Ixy1z);WSn=rand(1,length(CInWS));[i,j,k]=find(WSn<=p(8));pWS=sub2ind([1,length(CInWS)],i,j,k);
   A(CInWS(pWS))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInWN=setdiff(InWN,Ixy1z);WNn=rand(1,length(CInWN));[i,j,k]=find(WNn<=p(8));pWN=sub2ind([1,length(CInWN)],i,j,k);
   A(CInWN(pWN))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInTS=setdiff(InTS,Ixy1z); TSn=rand(1,length(CInTS));[i,j,k]=find(TSn<=p(8));pTS=sub2ind([1,length(CInTS)],i,j,k);
   A(CInTS (pTS))=0;
   [x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInBS=setdiff(InBS,Ixy1z);BSn=rand(1,length(CInBS));[i,j,k]=find(BSn<=p(8));pBS=sub2ind([1,length(CInBS)],i,j,k);
   A(CInBS(pBS))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInBN=setdiff(InBN,Ixy1z); BNn=rand(1,length(CInBN));[i,j,k]=find(BNn<=p(8));pBN=sub2ind([1,length(CInBN)],i, j,k);
   A(CInBN (pBN))=0;
[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

CInTN=setdiff(InTN,Ixy1z); TNn=rand(1,length(CInTN));[i,j,k]=find(TNn<=p(8));pTN=sub2ind([1,length(CInTN)],i, j,k);
   A(CInTN (pTN))=0;
   [x,y1,z]=find(A==0);
   Ixy1z=sub2ind([kn,kn,kn],x,y1,z);

[x,y1,z]=find(A==0);
Ixy1z=sub2ind([kn,kn,kn],x,y1,z);
gn=length(Ixy1z);
Vs=gn/N %體積分?jǐn)?shù)

m=m+1
pause(0.1)
imshow(A,'InitialMagnification','fit')%更新圖像
title('time=','Fontsize',14,'Fontname','Times new roman'); % 更新時(shí)間參數(shù)

end

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專(zhuān)碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專(zhuān)碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

1樓 2015-01-22 11:12:02

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

mengyuqiya

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.4
帖子: 9
在線: 12.4小時(shí)
蟲(chóng)號(hào): 4292046
注冊(cè): 2015-12-16
專(zhuān)業(yè): 建筑物理

樓主程序?qū)懗鰜?lái)了么我最近也在弄這個(gè) 不知道怎么做

贊一下

回復(fù)此樓

2樓2016-06-22 19:51:30

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

mengyuqiya

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.4
帖子: 9
在線: 12.4小時(shí)
蟲(chóng)號(hào): 4292046
注冊(cè): 2015-12-16
專(zhuān)業(yè): 建筑物理

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

你最早的A生成的就是個(gè)二維數(shù)組所以后面全是二維計(jì)算的
inshow好像也只能畫(huà)二維的圖像

贊一下

回復(fù)此樓

3樓2016-06-22 19:53:48

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

abc5852362

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 2
在線: 4.9小時(shí)
蟲(chóng)號(hào): 9702222
注冊(cè): 2018-07-31

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

Z 變量沒(méi)有應(yīng)用上，基本不起作用

贊一下

回復(fù)此樓

4樓2018-08-06 17:03:37

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主愛(ài)笑的眼睛na 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	12/600	2026-03-18 15:27 by 咪嗚喵嗚
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +7	rare12345 2026-03-18	7/350	2026-03-18 14:31 by laoshidan
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +6	慕寒mio 2026-03-16	6/300	2026-03-18 14:26 by 007_lilei
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂(lè)
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料專(zhuān)碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門(mén)版塊排行榜

愛(ài)笑的眼睛na

[求助] 隨機(jī)生長(zhǎng)法模擬多孔介質(zhì)的matlab二維代碼改寫(xiě)三維，結(jié)果出來(lái)還是二維 已有2人參與

» 猜你喜歡

mengyuqiya

mengyuqiya

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

abc5852362

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

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 隨機(jī)生長(zhǎng)法模擬多孔介質(zhì)的matlab二維代碼改寫(xiě)三維，結(jié)果出來(lái)還是二維已有2人參與