云模型代码:
N=3000;
Ex=1;En=1/3;He=0. 1;
CloudDrp = zeros(1,N);
for i=1:N
E_n = normrnd(En,He,1,1);
CloudDrp(1,i) = normrnd(Ex,E_n,1,1);
CloudDrp(2,i) = exp(-(CloudDrp(1,i)-Ex)^2/(2*E_n^2));
end
plot(CloudDrp(1,:),CloudDrp(2,:),'.')
hold on
N=3000;
Ex=5. /2;En=1.5/6;He=0. 1;
CloudDrp = zeros(1,N);
for i=1:N
E_n = normrnd(En,He,1,1);
CloudDrp(1,i) = normrnd(Ex,E_n,1,1);
CloudDrp(2,i) = exp(-(CloudDrp(1,i)-Ex)^2/(2*E_n^2));
end
plot(CloudDrp(1,:),CloudDrp(2,:),'.')
hold on
N=3000;
Ex=8.5/2;En=2.5/6;He=0. 1;
CloudDrp = zeros(1,N);
for i=1:N
E_n = normrnd(En,He,1,1);
CloudDrp(1,i) = normrnd(Ex,E_n,1,1);
CloudDrp(2,i) = exp(-(CloudDrp(1,i)-Ex)^2/(2*E_n^2));
end
plot(CloudDrp(1,:),CloudDrp(2,:),'.')
hold on
N=3000;
Ex=11/2;En=1/6;He=0.1;
CloudDrp = zeros(1,N);
for i=1:N
E_n = normrnd(En,He,1,1);
CloudDrp(1,i) = normrnd(Ex,E_n,1,1);
CloudDrp(2,i) = exp(-(CloudDrp(1,i)-Ex)^2/(2*E_n^2));
end
plot(CloudDrp(1,:),CloudDrp(2,:),'.')
ylabel('确定度 ');
xlabel('岩石弹性变形能指数 Wet');
axis([0,9,0,1]);
legend('Ⅰ级', 'Ⅱ级', 'Ⅲ级','Ⅳ级');
grid on