📄 caizhong7.m
字号:
clc
clear all
format long
caizhong7 =[1 7 10 15 18 26 29
11 15 17 18 20 22 28
9 15 18 19 21 25 28
6 8 12 15 23 25 27
11 13 15 18 19 22 29
3 8 10 14 15 17 22
4 11 15 17 18 20 27
2 3 11 13 15 24 28
8 9 17 18 19 24 25
1 4 9 16 18 20 26
1 2 3 12 15 22 24
3 10 11 19 20 26 27
1 2 5 11 14 16 21
12 14 17 21 22 29 30
1 4 12 19 22 23 25
1 4 12 14 22 25 26
6 8 10 21 25 28 29
3 4 10 12 20 28 30
2 10 12 19 23 25 27
3 7 12 17 19 26 28
2 3 13 17 21 23 29
4 5 6 16 24 25 27
4 10 11 14 21 27 28
3 17 19 20 22 26 30
1 3 8 12 18 22 23
1 6 22 23 24 26 29
7 20 22 24 27 29 30
7 8 10 11 13 15 19
1 2 4 9 11 25 30
1 2 4 13 19 27 28
3 6 8 9 19 21 27
4 11 13 17 28 29 30
5 6 9 17 20 22 27
13 15 16 18 19 22 23
6 9 22 25 26 27 29
2 6 7 13 20 22 29
10 13 14 15 20 23 30
7 8 12 19 21 24 27
1 3 10 14 22 25 26
1 4 8 16 17 18 24
3 5 6 12 13 21 25
1 5 8 11 13 25 29
7 12 19 27 28 29 30
4 9 12 18 20 22 27
2 3 7 11 19 26 29
2 12 17 18 24 26 29
3 6 17 21 25 26 28
8 14 16 17 19 27 30
2 7 10 16 20 25 29
1 4 8 10 13 15 25
1 4 12 13 19 20 28
3 6 10 18 22 23 24
3 7 12 16 17 26 28
3 4 12 15 20 21 28
5 8 13 15 19 26 29
2 3 8 13 16 19 23
3 11 12 14 18 20 28
4 8 13 19 21 27 28
3 6 7 11 20 29 30
4 9 13 18 19 27 29
1 20 23 24 25 27 29
1 8 19 20 24 28 29
1 3 11 19 20 21 26
10 17 20 22 27 28 30
1 6 17 18 23 28 29
3 4 7 10 22 23 30
11 18 20 22 25 26 27
2 4 10 12 16 26 29
1 2 5 18 22 23 25
1 3 12 13 15 23 24
6 10 13 23 26 28 30
1 9 10 13 17 20 30
1 3 4 11 18 25 28
3 5 12 14 20 21 24
1 17 18 23 24 26 28
5 8 16 19 26 28 29
1 3 6 9 12 25 27
10 11 15 18 20 21 27
3 14 15 16 17 18 26
5 8 17 21 22 24 26
2 3 7 9 25 26 27
2 8 23 27 28 29 30
1 5 13 22 23 27 28
2 3 5 8 11 19 23
2 5 14 21 22 23 27
1 6 11 16 22 27 29
4 22 23 24 25 26 29
5 6 13 16 18 21 28
1 2 10 20 22 23 27
1 7 14 16 22 27 30
2 3 4 10 14 23 29
7 9 11 18 22 27 28
9 10 20 22 24 26 27
1 4 11 12 20 22 30
5 11 13 18 19 22 23
2 12 15 16 18 24 27
2 4 6 10 21 25 30
1 3 4 5 13 20 30
4 13 19 21 25 27 29
6 7 13 16 19 24 29
10 13 19 20 24 25 30
4 9 17 20 22 25 26
5 6 7 14 21 27 28
6 15 17 21 22 27 29
2 5 17 21 25 27 29
1 2 6 8 11 12 20
3 4 6 12 16 21 29
1 6 11 16 18 22 29
3 4 10 16 23 25 29
8 10 11 12 15 20 21
1 5 8 12 17 18 23
4 7 12 14 21 26 28
10 23 24 27 28 29 30
2 3 4 7 17 19 20
2 4 14 16 22 24 27
5 9 15 17 23 26 28
2 4 14 17 20 25 28
1 3 7 11 12 14 17
3 4 10 17 18 23 27
3 6 15 16 18 28 29
1 3 7 10 12 14 29
2 5 12 17 19 20 27
2 6 7 8 9 21 26
7 16 17 18 19 21 25
4 7 10 14 22 27 28
2 9 12 13 24 25 27
9 13 15 17 19 22 29
3 4 5 6 9 10 18
2 10 14 20 22 27 30
2 5 7 9 23 28 30
5 6 12 18 19 27 28
6 7 8 9 11 19 21
6 7 10 20 25 26 28
5 6 13 15 21 26 27
3 9 11 12 13 18 19
13 16 21 22 23 25 28
1 4 8 10 18 20 26
1 13 14 16 27 28 29
7 11 12 16 17 18 19
1 4 6 11 13 17 29
8 10 11 15 17 22 26
1 4 6 7 8 10 11
1 3 5 9 13 20 26
3 9 19 20 21 24 28
4 14 21 24 26 29 30
6 7 9 18 21 24 29
5 7 8 10 17 28 29
5 7 10 19 20 25 30
3 9 11 13 15 18 21
7 16 18 19 20 22 28
5 9 13 16 19 26 29
2 7 11 13 17 18 25
1 12 14 17 20 25 29
7 11 13 15 18 26 28
3 7 12 14 20 22 24
1 9 11 15 16 28 30
1 5 9 14 24 27 30
2 7 14 16 19 20 26
4 6 7 11 19 20 24
2 6 7 8 13 26 27
1 3 4 12 25 26 27
9 12 14 15 24 25 26
3 8 9 11 13 17 29
2 6 18 19 24 25 26
1 2 10 12 17 24 25
1 6 9 12 15 18 30
5 11 13 16 20 21 28
3 20 21 22 25 27 28
1 3 15 23 27 29 30
1 5 6 9 15 20 22
7 9 10 18 19 23 30
5 13 17 26 27 29 30
6 7 9 13 14 15 23
7 12 17 23 24 27 29
8 12 13 14 16 19 21
2 6 12 19 20 22 25
2 5 14 16 22 25 28
3 10 14 17 20 21 25
3 4 8 13 24 26 29
2 4 8 10 13 25 28
2 5 10 20 21 27 28
5 9 11 19 23 24 25
3 12 17 18 23 24 28
3 8 9 12 18 19 23
6 7 9 10 22 28 30
2 8 9 10 18 20 24
5 7 11 12 13 19 21
8 9 11 15 16 17 30
4 5 7 16 18 20 25
3 12 16 18 24 27 30
6 8 10 12 22 24 27
4 8 13 14 16 17 24
5 8 16 21 22 29 30
1 3 10 11 15 27 29
3 4 8 13 14 15 21
4 9 13 15 22 25 28
7 9 10 18 21 24 27
1 7 12 17 22 25 29
1 9 15 16 21 25 28
9 10 12 17 27 28 29
1 3 7 12 13 18 23
6 11 12 13 19 21 24
1 3 4 11 15 19 26
7 10 11 17 18 20 30
11 12 14 15 17 19 29
5 9 13 17 21 22 28
1 2 3 16 17 19 24
7 9 15 18 20 23 30
8 10 12 17 23 25 28
5 6 10 12 22 27 29
3 4 10 11 15 20 25
5 14 15 16 18 21 30
6 12 13 17 19 24 28
1 9 10 19 20 25 26
1 2 7 12 14 29 30
2 3 4 6 11 20 29
2 5 11 16 19 22 30
5 7 8 9 13 14 21
7 9 12 13 15 19 28
1 7 11 13 14 21 24
1 5 12 14 20 27 30
1 2 6 13 22 26 30
8 13 17 19 21 22 24
8 9 10 11 18 22 25
2 5 6 7 10 14 23
1 3 10 11 17 25 30
1 4 8 11 16 18 28
2 8 13 15 18 21 28
3 6 8 19 23 24 26
1 4 7 11 18 27 30
5 8 10 11 13 17 18
1 9 12 16 17 28 29
6 10 11 16 21 28 29
2 5 6 10 17 21 29
5 6 8 9 12 13 24
10 12 13 16 21 28 30
1 5 6 7 17 19 27
5 8 9 10 18 20 22
3 12 18 19 21 26 30
5 9 17 18 20 27 29
10 16 20 22 23 25 28
2 4 5 14 15 18 24
3 4 13 14 16 17 20
6 9 11 15 18 22 25
3 4 9 11 12 23 25
3 5 11 14 21 22 24
5 13 14 19 22 25 26
2 6 12 14 21 23 27
1 4 11 13 16 23 26
3 6 11 13 14 21 23
2 9 16 17 23 27 30
1 5 8 17 19 23 29
2 7 11 22 25 27 28
4 9 13 15 22 27 28
1 3 5 8 15 23 26
3 9 10 14 21 25 28
2 3 12 13 20 26 27
3 5 6 9 17 21 27
2 3 17 18 23 25 28
2 4 5 6 8 14 26
14 17 19 23 24 25 26
7 18 20 22 23 26 28
9 12 16 18 19 22 30
6 9 13 18 20 24 29
6 14 18 23 24 28 30
1 2 7 8 13 14 16
3 4 9 17 21 28 29
15 17 18 21 22 24 29
9 19 20 23 25 27 29
4 13 14 23 26 28 29
3 6 7 17 26 28 30
2 3 11 13 15 20 26
5 6 7 8 11 20 26
2 4 6 10 17 19 28
5 7 8 12 13 16 21
3 12 16 18 23 25 30
9 15 21 22 23 29 30];
nnii=zeros(100,7);
for ii=2:100
cai=cai11(1:ii,:);
[m,n]=size(cai);
h=zeros(30,1);
for i=1:m
for j=1:n
if cai(i,j)==0
else
h(cai(i,j))=h(cai(i,j))+1;
end
end
end
h1=h/(m*n);
h=zeros(30,1);
for i=1:30
x=[1 0 0 0 0 0 0];
x(1)=i*x(1);
dd=x;
cai1=([dd;cai]);
h(i)=caientropy(cai1);
end
d=max(h);
n11=find(h==d);
n1=n11(1);
h=zeros(30,1);
for i=1:30
if i~=n1
y=[1 0 0 0 0 0];
y(1)=i*y(1);
dd=y;
cai1=([n1 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n22=find(h==d);
n2=n22(1);
h=zeros(30,1);
for i=1:30
if i~=n1& i~=n2
y=[1 0 0 0 0 ];
y(1)=i*y(1);
dd=y;
cai1=([n1 n2 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n33=find(h==d);
n3=n33(1);
h=zeros(30,1);
for i=1:30
if i~=n1& i~=n2 & i~=n3
y=[1 0 0 0 ];
y(1)=i*y(1);
dd=y;
cai1=([n1 n2 n3 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n44=find(h==d);
n4=n44(1);
h=zeros(30,1);
for i=1:30
if i~=n1& i~=n2 & i~=n3 & i~=n4
y=[1 0 0];
y(1)=i*y(1);
dd=y;
cai1=([n1 n2 n3 n4 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n55=find(h==d);
n5=n55(1);
h=zeros(30,1);
for i=1:30
if i~=n1& i~=n2 & i~=n3 & i~=n4 & i~=n5
y=[1 0];
y(1)=i*y(1);
dd=y;
cai1=([n1 n2 n3 n4 n5 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n66=find(h==d);
n6=n66(1);
h=zeros(30,1);
for i=1:30
if i~=n1& i~=n2 & i~=n3 & i~=n4 & i~=n5 & i~=n6
y=ones(1,1);
y(1)=i*y(1);
dd=y;
cai1=([n1 n2 n3 n4 n5 n6 dd;cai]);
h(i)=caientropy(cai1);
end
end
d=max(h);
n7=find(h==d);
n7=n7(1);
% ii1=ii-9;
nnii(ii,:)=[n1 n2 n3 n4 n5 n6 n7];
end
nn12=sort(nnii,2);⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -