相似文件
换一批
经管之家送您一份
应届毕业生专属福利!
求职就业群
感谢您参与论坛问题回答
经管之家送您两个论坛币!
+2 论坛币
- shuju=xlsread('modelerbookP165.xlsx','a2:c15');
- %计算分组前
- p1=sum(shuju(:,3))/14;
- p2=1-p1;
- EntS=-(p1*log2(p1)+p2*log2(p2));
- %数组准备
- for i=1:14
- h=shuju(i,2);
- d1=zeros(14,1);
- d2=zeros(14,1);
- for j=1:14
- if h-shuju(j,2)>=0;
- d1(j,1)=j;
- else
- d2(j,1)=j;
- end
- d1(d1==0)=[];
- d2(d2==0)=[];
- end
- s1=length(d1);
- s2=length(d2);
- s=14;
- p11=sum(shuju(d1,3))/s1;
- p22=sum(shuju(d2,3))/s2;
- EntATS(i,1)=(s1/s)*((-p11*log2(p11)-(1-p11)*log2(1-p11)))+(s2/s)*((-p22*log2(p22)-(1-p22)*log2(1-p22)));
- k=log2(7)-(2*EntS-2*((-p11*log2(p11)-(1-p11)*log2(1-p11)))-2*((-p22*log2(p22)-(1-p22)*log2(1-p22))));
- cost(i,1)=(log2(s-1))/s+k/s;
- end
- gain=zeros(14,1);
- aaa=zeros(14,1);
- for i=1:14
- gain(i,1)=EntS-EntATS(i,1);
- end
- for i=1:14
- aaa(i,1)=gain(i,1)-cost(i,1);
- end
复制代码
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
|