基于蒙特卡罗法的二维随机裂隙模拟Matlab程序

clear clc
N=input（‘裂隙组数=‘）;
chang=input（‘生成域长=‘）;
kuan=input（‘生成域宽=‘）;
mianji=chang*kuan;                                                         %面积
LXDD=[];
for k=1:N
    midu=input（‘密度=‘）;
    zouxjz=input（‘走向均值=‘）;
    zouxbzc=input（‘走向标准差=‘）;
    jicjz=input（‘迹长均值=‘）;
    jicbzc=input（‘迹长标准差=‘）;
    xikanjz=input（‘隙宽均值=‘）;
    xikanbzc=input（‘隙宽标准差=‘）;
    n=mianji*midu;                                                         %条数
    zhongxidian=unifrnd（-chang/2chang/2n2）;                             %中心点
    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）;                                         %迹长
    xikuang=normrnd（xikan