二维离散裂隙matlab

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/2,chang/2,n,2); %中心点,n表示2维，生成在区域-chang/2,chang/2中的随机数 zhongxidian zouxiang=normrnd(zouxjz,zouxbzc,n,1); %走向，生成正态分布的随机数 m=jicjz;%迹长均值 v=jicbzc^2;%迹长标准差的平方 mu = log((m^2)/sqrt(v+m^2));%迹长,对数正态分布函数的MU sigma = sqrt(log(v/(m^2)+1));%迹长,对数正态分布的sigma jichang=lognrnd(mu,sigma,n,1); %迹长,产生一维的对数正态分布随机数 xikuang=normrnd(xikanjz,xikanbzc,n,1); %隙宽，生成正态分布的随机数 x1=zhongxidian(:,1)+(jichang./2).*cos(zouxiang*pi/180);%(:,1)表示zhongxidian坐标的第一列为x1轴数据 y1=zhongxidian(:,2)+(jichang./2).*sin(zouxiang*pi/180);%(:,1)表示zhongxidian坐标的第二列为y1轴数据 x2=zhongxidian(:,1)+(jichang./2).*cos((zouxiang+180)*pi/180);%(:,1)表示zhongxidian坐标的第一列为x2轴数据 y2=zhongxidian(:,2)+(jichang./2).*sin((zouxiang+180)*pi/180);%(:,1)表示zhongxidian坐标的第二列为y2轴数据 for i=1:n if x1(i)>x2(i) w=x2(i); x2(i)=x1(i); q=y2(i); y2(i)=y1(i); x1(i)=w; y1(i)=q; end end LXDD1=[x1,y1,x2,y2,zouxiang,xikuang]; if isempty(LXDD)%判断矩阵是否为空 LXDD=LXDD1;%若为空矩阵则直接将LXDD1的值赋给LXDD else LXDD=[LXDD;LXDD1];%若不为空矩阵则将LXDD1的值添加到矩阵LXDD中 end for i=1:n subplot(1,3,1) plot([x1(i),x2(i)],[y1(i),y2(i)]);%定义两点不同的坐标，并绘制直线 xlim([-chang/2,chang/2]);%X轴的坐标范围 ylim([-kuan/2,kuan/2]);%Y轴的坐标范围 hold on end end LXDD b=zeros(2*size(LXDD, 1),2);%产生一个2n*2的0矩阵 c_q=LXDD(:,1);%取出X1的元素 d_q=LXDD(:,2);%取出Y1的元素 e_q=LXDD(:,3);%取出X2的元素 f_q=LXDD(:,4);%取出Y2的元素 for i=1:size(LXDD, 1) b(2*i-1,1)=c_q(i); b(2*i,1)=e_q(i); b(2*i-1,2)=d_q(i); b(2*i,2)=f_q(i); end b