clear?clc
N=input（‘裂隙組數(shù)=‘）;
chang=input（‘生成域長(zhǎng)=‘）;
kuan=input（‘生成域?qū)?‘）;
mianji=chang*kuan;?????????????????????????????????????????????????????????%面積
LXDD=[];
for?k=1:N
????midu=input（‘密度=‘）;
????zouxjz=input（‘走向均值=‘）;
????zouxbzc=input（‘走向標(biāo)準(zhǔn)差=‘）;
????jicjz=input（‘跡長(zhǎng)均值=‘）;
????jicbzc=input（‘跡長(zhǎng)標(biāo)準(zhǔn)差=‘）;
????xikanjz=input（‘隙寬均值=‘）;
????xikanbzc=input（‘隙寬標(biāo)準(zhǔn)差=‘）;
????n=mianji*midu;?????????????????????????????????????????????????????????%條數(shù)
????zhongxidian=unifrnd（-chang/2chang/2n2）;?????????????????????????????%中心點(diǎn)
????zouxiang=normrnd（zouxjzzouxbzcn1）;??????????????????????????????????%走向
????m=jicjz;
????v=jicbzc^2;
????mu=log（（m^2）/sqrt（v+m^2））;
????sigma=sqrt（log（v/（m^2）+1））;
????jichang=lognrnd（musigman1）;?????????????????????????????????????????%跡長(zhǎng)
????xikuang=normrnd（xikan