logo资料库

正向云发生器云图代码.docx

第1页 / 共2页
第2页 / 共2页
资料共2页,全文预览结束
云模型代码: 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
分享到:
收藏