随机生成多孔介质matlab程序

%ÀîÈÊÃñclc;clearall;closeallmax_j=193;max_i=81;d11=0.005;d13=0.005;d14=0.001;d12=0.001;d15=0.0001;d16=0.0001;d17=0.0001;d18=0.0001;n=0.7;%¿×϶¶Ècdd=0.001;%Pcd¹ÌÏàÉú³¤ºË¸ÅÂÊnumtotal_need=(1-n)*max_i*max_j;numsoild=0;arrgrid=zeros(max_i,max_j);whilenumsoildnumtotal_need/40fori=1:max_iforj=1:max_jifabs(rand())cddnumsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endendendendTnumsoild=numsoild;%×ÜÉú³¤ºËÊý%figure(1)%imshow(~arrgrid);%1°×É«ÊÇ¿×϶£¬0ºÚÉ«ÊǹÌÌå%µÚ2²½£¬´Ó¹ÌÏàÄÚºËÏòÖÜΧ8¸ö·½ÏòÉú³¤whileTnumsoildnumtotal_needforindex_soild=1:Tnumsoildindex_i=soild(index_soild,1);index_j=soild(index_soild,2);%ºá×ŵÄÊÇy·½Ïò%Êú×ŵÄÊÇx·½Ïò%¹ÌÏàÄÚºËÏòÓÒÉú³¤1·½Ïòifindex_jmax_j-1i=index_i;j=index_j+1;%ÓÒ1ifarrgrid(i,j)==0&&rand()d11numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤2·½Ïòifindex_imax_i-1i=index_i+1;j=index_j;%ÓÒ1ifarrgrid(i,j)==0&&rand()d12numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤3·½Ïòifindex_j2+1i=index_i;j=index_j-1;ifarrgrid(i,j)==0&&rand()d13numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤4·½Ïòifindex_i2+1i=index_i-1;j=index_j;ifarrgrid(i,j)==0&&rand()d14numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤5·½Ïòifindex_jmax_j-1&&index_imax_i-1i=index_i+1;j=index_j+1;ifarrgrid(i,j)==0&&rand()d15numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤6·½Ïòifindex_imax_i-1&&index_j2+1i=index_i+1;j=index_j-1;%ÓÒ1ifarrgrid(i,j)==0&&rand()d16numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤7·½Ïòifindex_j2+1&&index_i2+1i=index_i-1;j=index_j-1;%ÓÒ1ifarrgrid(i,j)==0&&rand()d17numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endend%¹ÌÏàÄÚºËÏòÓÒÉú³¤8·½Ïòifindex_jmax_j-1&&index_i2+1i=index_i-1;j=index_j+1;%ÓÒ1ifarrgrid(i,j)==0&&rand()d18numsoild=numsoild+1;arrgrid(i,j)=1;soild(numsoild,1)=i;soild(numsoild,2)=j;endendend%index_soilTnumsoild=numsoild;end%endwhile%x=grid(:,1);%y=grid(:,2);%z=grid(:,3);%Íø¸ñÊôÐÔ0ÒºÏà1¹ÌÏà%[x,y]=meshgrid(1:max_i,1:max_j);figure(2)imshow(~arrgrid);%°×É«ÊÇ¿×϶£¬ºÚÉ«ÊǹÌÌå%colormap(gray());imagesc(arrgrid);%forj=1:max_j%scatter(arrgrid(:,j),'filled');%end%holdon%fori=1:max_i%scatter(arrgrid(i,:),'filled')%end%%%»æÖƹÌÌåÂÖÀªÏßcΪÂÖÀªÊý¾Ý%figure(2)%[c,h]=contour(x,y,arrgrid,[11]);