随机生成多孔介质matlab程序.doc

% clc;clear all;close allmax_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; % Pcdnumtotal_need=(1-n)*max_i*max_j;numsoild=0;arrgrid=zeros(max_i,max_j); while numsoildnumtotal_need/40 for i=1:max_i for j=1:max_j if abs(rand()cdd numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end end endTnumsoild=numsoild; % % figure(1)% imshow(arrgrid) ; % 10% 28while Tnumsoildnumtotal_need for index_soild=1:Tnumsoild index_i=soild(index_soild,1); index_j=soild(index_soild,2); % y % x % 1 if index_jmax_j-1 i=index_i;j=index_j+1; % 1 if arrgrid(i,j)=0 & rand()d11 numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end % 2 if index_imax_i-1 i=index_i+1;j=index_j; % 1 if arrgrid(i,j)=0 & rand()2+1 i=index_i;j=index_j-1; if arrgrid(i,j)=0 & rand()2+1 i=index_i-1;j=index_j; if arrgrid(i,j)=0 & rand()d14 numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end % 5 if index_jmax_j-1 & index_imax_i-1 i=index_i+1;j=index_j+1; if arrgrid(i,j)=0 & rand()d15 numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end % 6 if index_i2+1 i=index_i+1;j=index_j-1; % 1 if arrgrid(i,j)=0 & rand()2+1 & index_i2+1 i=index_i-1;j=index_j-1; % 1 if arrgrid(i,j)=0 & rand()d17 numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end % 8 if index_j2+1 i=index_i-1;j=index_j+1; % 1 if arrgrid(i,j)=0 & rand()d18 numsoild=numsoild+1;arrgrid(i,j)=1; soild(numsoild,1)=i;soild(numsoild,2)=j; end end end % index_soil Tnumsoild=numsoild;end % end while% 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);% for j=1:max_j% scatter(arrgrid(:,j),filled);% end% hold on