ss_random_thicknesses.m

来自「地震、测井方面matlab代码,解释的比较详细」· M 代码 · 共 19 行

function rnr=ss_random_thicknesses(n,m,minthick,maxthick)
% Create "n" by "m" matrix of random thicknesses between 
% "minthick" and "maxthick". P(x) ~ 1/x

lmax=length(maxthick);
rnr=rand(n,m);

if lmax == 1
   rnr=minthick*((maxthick/minthick).^rnr);

elseif lmax == m
   maxthick=maxthick./minthick;
   rnr=minthick*(maxthick(ones(n,1),:).^rnr);
   
else
   error([' "maxthick" must be a constant or a row vector of length ',num2str(m)])

end