yun_云模型隶属度_云模型_云模型隶属度

clear; N=1000; Ex1=0.666; En1=0.136; He1=0.062; CloudDrp = zeros(2,N); for i=1:N E_n1 = normrnd(En1,He1,1,1); CloudDrp(1,i) = normrnd(Ex1,E_n1,1,1); CloudDrp(2,i) = exp(-(CloudDrp(1,i)-Ex1)^2/(2*E_n1^2)); end plot(CloudDrp(1,:),CloudDrp(2,:),'.') 标准云 1.标准云的生成 clear all clc a1=[0,1]; a2=[1 2]; a3=[2,3]; a4=[3 4]; a5=[4,5]; Ex=(max(a1)+min(a1))/2 En=(max(a1)-min(a1))/6 He=0.06 dd=1000; % 云滴数量 for i=1:dd %E1=normrnd(En,He,1); %X=normrnd(Ex,E1,1); %Enn=randn(1)*He+En; %x(i)=randn(1)*Enn+Ex; Enn=normrnd(En,He,1); x1(i)=normrnd(Ex,Enn,1); y1(i)=exp(-(x1(i)-Ex)^2/(2*Enn^2)); end plot(x1,y1,'*g') hold on Ex=(max(a2)+min(a2))/2 En=(max(a2)-min(a2))/6 He=0.06 hold on for i=1:dd Enn=normrnd(En,He,1); x2(i)=normrnd(Ex,Enn,1); y2(i)=exp(-(x2(i)-Ex)^2/(2*Enn^2)); end plot(x2,y2,'Oc') Ex=(max(a3)+min(a3))/2 En=(max(a3)-min(a3))/6 He=0.06 for i=1:dd %Enn=randn(1)*He+En; %x(i)=randn(1)*Enn+Ex; Enn=normrnd(En,He,1); x3(i)=normrnd(Ex,Enn,1); y3(i)=exp(-(x3(i)-Ex)^2/(2*Enn^2)); end plot(x3,y3,'^y') Ex=(max(a4)+min(a4))/2 En=(max(a4)-min(a4))/6 He=0.06 for i=1:dd %Enn=randn(1)*He+En; %x(i)=randn(1)*Enn+Ex; Enn=normrnd(En,He,1); x4(i)=normrnd(Ex,Enn,1); y4(i)=exp(-(x4(i)-Ex)^2/(2*Enn^2)); end plot(x4,y4,'+m') Ex=(max(a5)+min(a5))/2 En=(max(a5)-min(a5))/6 He=0.06 for i=1:dd %Enn=randn(1)*He+En; %x(i)=randn(1)*Enn+Ex; Enn=normrnd(En,He,1); x5(i)=normrnd(Ex,Enn,1); y5(i)=exp(-(x5(i)-Ex)^2/(2*Enn^2)); end plot(x5,y5,'pr') %title('XXXXX') ylabel('隶属度'); xlabel('山岭公路隧道施工安全风险数域'); axis([0,5,0,1]) %grid on legend('低度风险','较低风险','中度风险','较高风险','高度风险') 熵权法求权值 %熵权法，数据一定要经过同一化 clc;clear; load('hh.mat') x=load('data.txt');%（nxm）数据矩阵，m个指标，n个评价对象; x1=x(:,1:3); x2=x(:,4:9); x3=x(:,10:13); x4=x(:,14:17); x5=x(:,18:20); %% 权重1 hang1=size(x1,1);%n个评价对象 lie1=size(x1,2);%m个评价指标 for j=1:lie1 %计算第j个指标下第i个项目的指标值的比重 for i=1:hang1 p1(i,j)=x1(i,j)/sum(x1(:,j)); end end k1=1/(log(hang1)); for j=1:lie1 % 计算第j个指标的熵值 e1(j)=-k1*sum(p1(:,j).*log(p1(:,j))); end for j=1:lie1 %计算第j个指标的熵权 w1(j)=(1-e1(j))/sum(1-e1); end w1 %最后计算出来的权值 %% 权重2 hang2=size(x2,1);%n个评价对象 lie2=size(x2,2);%m个评价指标 for j=1:lie2 %计算第j个指标下第i个项目的指标值的比重 for i=1:hang2 p2(i,j)=x2(i,j)/sum(x2(:,j)); end end k2=1/(log(hang2)); for j=1:lie2 % 计算第j个指标的熵值 e2(j)=-k2*sum(p2(:,j).*log(p2(:,j))); end for j=1:lie2 %计算第j个指标的熵权 w2(j)=(1-e2(j))/sum(1-e2); end w2 %最后计算出来的权值 %% 权重3 hang3=size(x3,1);%n个评价对象 lie3=size(x3,2);%m个评价指标 for j=1:lie3 %计算第j个指标下第i个项目的指标值的比重 for i=1:hang3 p3(i,j)=x3(i,j)/sum(x3(:,j)); end end k3=1/(log(hang3)); for j=1:lie3 % 计算第j个指标的熵值 e3(j)=-k3*sum(p3(:,j).*log(p3(:,j))); end for j=1:lie3 %计算第j个指标的熵权 w3(j)=(1-e3(j))/sum(1-e3); end w3 %最后计算出来的权值 %% 权重4 hang4=size(x4,1);%n个评价对象 lie4=size(x4,2);%m个评价指标 for j=1:lie4 %计算第j个指标下第i个项目的指标值的比重 for i=1:hang4 p4(i,j)=x4(i,j)/sum(x4(:,j)); end end k4=1/(log(hang4)); for j=1:lie4 % 计算第j个指标的熵值 e4(j)=-k4*sum(p4(:,j).*log(p4(:,j))); end for j=1:lie4 %计算第j个指标的熵权 w4(j)=(1-e4(j))/sum(1-e4); end w4 %最后计算出来的权值 %% 权重5 hang5=size(x5,1);%n个评价对象 lie5=size(x5,2);%m个评价指标 for j=1:lie5 %计算第j个指标下第i个项目的指标值的比重 for i=1:hang5 p5(i,j)=x5(i,j)/sum(x5(:,j)); end end k5=1/(log(hang5)); for j=1:lie5 % 计算第j个指标的熵值 e5(j)=-k5*sum(p5(:,j).*log(p5(:,j))); end for j=1:lie5 %计算第j个指标的熵权 w5(j)=(1-e5(j))/sum(1-e5); end w5 %最后计算出来的权值 %% 第二层权重 x=load('data1.txt');%（nxm）数据矩阵，m个指标，n个评价对象; % 第二层权重 hang=size(x,1);%n个评价对象 lie=size(x,2);%m个评价指标 for j=1:lie %计算第j个指标下第i个项目的指标值的比重 for i=1:hang p(i,j)=x(i,j)/sum(x(:,j)); end end k1=1/(log(hang)); for j=1:lie % 计算第j个指标的熵值 e(j)=-k1*sum(p(:,j).*log(p(:,j))); end for j=1:lie %计算第j个指标的熵权 w0(j)=(1-e(j))/sum(1-e); end w0 %最后计算出来的权值 Ex1=sum(hh(1,1:3).*hh(2,1:3).*w1)/sum(hh(2,1:3).*w1); Ex2=sum(hh(1,4:9).*hh(2,4:9).*w2)/sum(hh(2,4:9).*w2); Ex3=sum(hh(1,10:13).*hh(2,10:13).*w3)/sum(hh(2,10:13).*w3); Ex4=sum(hh(1,14:17).*hh(2,14:17).*w4)/sum(hh(2,14:17).*w4); Ex5=sum(hh(1,18:20).*hh(2,18:20).*w5)/sum(hh(2,18:20).*w5); En1=sum(hh(2,1:3).*w1); En2=sum(hh(2,4:9).*w2); En3=sum(hh(2,10:13).*w3); En4=sum(hh(2,14:17).*w4); En5=sum(hh(2,18:20).*w5); He1=sum(hh(3,1:3).*hh(2,1:3).*w1)/sum(hh(2,1:3).*w1); He2=sum(hh(3,4:9).*hh(2,4:9).*w2)/sum(hh(2,4:9).*w2); He3=sum(hh(3,10:13).*hh(2,10:13).*w3)/sum(hh(2,10:13).*w3); He4=sum(hh(3,14:17).*hh(2,14:17).*w4)/sum(hh(2,14:17).*w4); He5=sum(hh(3,18:20).*hh(2,18:20).*w5)/sum(hh(2,18:20).*w5); hh' c1_5=[Ex1 Ex2 Ex3 Ex4 Ex5;En1 En2 En3 En4 En5;He1 He2 He3 He4 He5]; hhh=c1_5 hhh' Ex0=sum(hhh(1,1:end).*hhh(2,1:end).*w0)/sum(hhh(2,1:end).*w0); En0=sum(hhh(2,1:end).*w0); He0=sum(hhh(3,1:end).*hhh(2,1:end).*w0)/sum(hhh(2,1:end).*w0); c0=[Ex0 En0 He0] save c1_5 c1_5 save c0 c0