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

% à?èê?? clc;clear all;close all max_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; % Pcd1ì?àéú3¤o????ê numtotal_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 end Tnumsoild=numsoild; % ×üéú3¤o?êy % figure(1) % imshow(~arrgrid) ; %