clear;
max_j?=?200;?max_i=200;?d14=?0.1;?d58=?d14/4;?n=?0.6;
cdd=0.01;?numtotal_need=?（1-?n）?*?max_i*?max_j;
numsoild=0;?arrgrid=0;
for?i=1:?max_i?%第1?步?遍歷所有網(wǎng)格?生成固相內(nèi)核
???for?j?=?1:max_j
?????if?rand（?）????????numsoild=?numsoild+?1;arrgrid（ij）?=?1;
???????soild（numsoild1）?=?i;soild（numsoild?2）?=?j;
???????
?????end?%?end?if
???end?%?end?j
end?%?end?i
Tnumsoild=?numsoild;
%?第2?步?從固相內(nèi)核向周圍8?個(gè)方向生長(zhǎng)
while?Tnumsoild??for?index_soild=?1:?Tnumsoild
????index_i=?soild（index_soild?1）;
????index_j=?soild（index_soild?2）;
?????%?固相內(nèi)核向右邊生長(zhǎng)?圖1?中右1?方向
?????if?index_j???????i=?index_i;?j?=?index_j?+?1;?%?right?1
???????if?（arrgrid（?i?j?）==?0?&?rand（?）???????????numsoild=?numsoild+?1;?arrgrid（?i?j?）?=?1;
??????????soild（?numsoild?1）?=?i;?soild（numsoild?2）?=?j;
???????end?
?????end?%index_j
%?固相內(nèi)核向上邊生長(zhǎng)?圖1?中右2?方向
?????if?index_i>?1
???????i=?index_i-1;?j?=?index_j;?%?2
???????if?（arrgrid（?i?j?）==?0?&?rand（?）???????????numsoild=?numsoild+?1;?arrgrid（?i?j?）?=?1;
??????????soild（?numsoild?1）?=?i;?soild（numsoild?2）?=?j;
???????end?
?????end?%index_j
%?固相內(nèi)核向左邊生長(zhǎng)?圖1?中右3?方向
?????if?index_j>?1
???????i=?index_i;?j?=?index_j-1;?%?3
???????if?（arrgrid（?i?j?）==?0?&?rand（?）???????????numsoild=?numsoild+?1;?arrgrid（?i?j?）?=?1;
??????????soild（?numsoild?1）?=?i;?soild（numsoild?2）?=?j;
???????end?
?????end?%index_j
%?固相內(nèi)核向下邊生長(zhǎng)?圖1?中右4?方向
?????if?index_i???????i=?index_i+1;?j?=?index_j;?%?4
?????