jinlin.m

在matlab环境下

clear all;
close all;

x1(1,1)=1.8796; x1(1,2)=1.8041;
x1(2,1)=2.6801; x1(2,2)=2.7526;
x1(3,1)=3.6284; x1(3,2)=1.3313;
x1(4,1)=4.7302; x1(4,2)=3.0267;
x1(5,1)=5.7865; x1(5,2)=0.3089;
x1(6,1)=7.1831; x1(6,2)=2.9453;
x1(7,1)=7.2395; x1(7,2)=3.6268;
x1(8,1)=8.0763; x1(8,2)=4.9714;
x1(9,1)=9.9172; x1(9,2)=3.9551;
x1(10,1)=11.7397; x1(10,2)=3.9500;
x1(11,1)=12.8685; x1(11,2)=2.4619;
x1(12,1)=12.5289; x1(12,2)=3.5313;
x1(13,1)=13.3206; x1(13,2)=4.4376;
x1(14,1)=15.7457; x1(14,2)=0.9094;
x1(15,1)=15.4758; x1(15,2)=5.2898;
x1(16,1)=17.2917; x1(16,2)=5.2197;
x1(17,1)=18.9338; x1(17,2)=3.7324;
x1(18,1)=19.3299; x1(18,2)=2.0778;
x1(19,1)=20.7408; x1(19,2)=5.2698;
x1(20,1)=20.0199; x1(20,2)=3.5670;
x1(21,1)=21.2740; x1(21,2)=4.7658;
x1(22,1)=23.6375; x1(22,2)=3.3211;
x1(23,1)=23.8603; x1(23,2)=6.1293;
x1(24,1)=25.7806; x1(24,2)=1.8003;
x1(25,1)=26.4698; x1(25,2)=4.3295;
x1(26,1)=27.3746; x1(26,2)=3.1499;
x1(27,1)=27.6922; x1(27,2)=6.1123;
x1(28,1)=28.3321; x1(28,2)=3.6388;
x1(29,1)=29.3112; x1(29,2)=5.5035;
x1(30,1)=30.3822; x1(30,2)=2.5172;
x1(31,1)=31.8449; x1(31,2)=4.1858;
x1(32,1)=33.7120; x1(32,2)=5.0515;
x1(33,1)=33.9805; x1(33,2)=4.8947;
x1(34,1)=35.6319; x1(34,2)=5.7023;
x1(35,1)=35.9215; x1(35,2)=6.1456;
x1(36,1)=36.9147; x1(36,2)=3.8067;
x1(37,1)=37.9014; x1(37,2)=7.9138;
x1(38,1)=38.8244; x1(38,2)=7.3828;
x1(39,1)=40.8032; x1(39,2)=7.7581;
x1(40,1)=40.0112; x1(40,2)=8.0748;
x1(41,1)=41.5948; x1(41,2)=7.5525;
x1(42,1)=42.0983; x1(42,2)=5.4144;
x1(43,1)=44.3864; x1(43,2)=5.9879;
x1(44,1)=45.3002; x1(44,2)=7.9712;
x1(45,1)=46.9660; x1(45,2)=7.7468;
x1(46,1)=47.1053; x1(46,2)=5.5875;
x1(47,1)=47.8001; x1(47,2)=5.9673;
x1(48,1)=48.3976; x1(48,2)=7.1165;
x1(49,1)=50.2504; x1(49,2)=8.0479;
x1(50,1)=51.4667; x1(50,2)=8.6202;
x1(51,1)=49.7518; x1(51,2)=11.0474;
x1(52,1)=48.0198; x1(52,2)=9.7412;
x1(53,1)=47.8397; x1(53,2)=8.6673;
x1(54,1)=47.5073; x1(54,2)=9.6810;
x1(55,1)=46.5877; x1(55,2)=10.5484;
x1(56,1)=45.8399; x1(56,2)=8.6472;
x1(57,1)=44.6894; x1(57,2)=12.2699;
x1(58,1)=42.7355; x1(58,2)=13.1906;
x1(59,1)=42.2416; x1(59,2)=11.6802;
x1(60,1)=41.4626; x1(60,2)=9.1437;
x1(61,1)=39.3878; x1(61,2)=13.3631;
x1(62,1)=39.8096; x1(62,2)=12.6606;
x1(63,1)=38.1384; x1(63,2)=13.4300;
x1(64,1)=37.2636; x1(64,2)=10.7010;
x1(65,1)=35.4688; x1(65,2)=12.6869;
x1(66,1)=35.0976; x1(66,2)=12.7679;
x1(67,1)=34.8632; x1(67,2)=12.0533;
x1(68,1)=32.6704; x1(68,2)=15.7258;
x1(69,1)=32.3111; x1(69,2)=16.0957;
x1(70,1)=30.7838; x1(70,2)=14.5081;
x1(71,1)=30.2546; x1(71,2)=17.3737;
x1(72,1)=29.3982; x1(72,2)=13.6487;
x1(73,1)=27.7944; x1(73,2)=17.5663;
x1(74,1)=26.8273; x1(74,2)=17.3489;
x1(75,1)=26.3104; x1(75,2)=15.9892;
x1(76,1)=25.6752; x1(76,2)=17.1196;
x1(77,1)=23.7432; x1(77,2)=19.0045;
x1(78,1)=22.8505; x1(78,2)=17.6571;
x1(79,1)=22.1893; x1(79,2)=15.8862;
x1(80,1)=21.1315; x1(80,2)=16.5870;
x1(81,1)=20.4331; x1(81,2)=15.9183;
x1(82,1)=19.0226; x1(82,2)=17.5691;
x1(83,1)=18.5528; x1(83,2)=17.2806;
x1(84,1)=16.9787; x1(84,2)=17.6517;
x1(85,1)=15.3718; x1(85,2)=18.4702;
x1(86,1)=15.4013; x1(86,2)=15.8341;
x1(87,1)=14.9654; x1(87,2)=17.0939;
x1(88,1)=13.6133; x1(88,2)=18.5902;
x1(89,1)=12.4071; x1(89,2)=16.4305;
x1(90,1)=10.9699; x1(90,2)=18.6493;
x1(91,1)=9.2292; x1(91,2)=16.7441;
x1(92,1)=9.3297; x1(92,2)=18.4027;
x1(93,1)=7.7307; x1(93,2)=17.9292;
x1(94,1)=6.2801; x1(94,2)=17.6374;
x1(95,1)=6.1335; x1(95,2)=17.9167;
x1(96,1)=5.6460; x1(96,2)=19.7987;
x1(97,1)=4.3479; x1(97,2)=16.6548;
x1(98,1)=3.9989; x1(98,2)=16.7955;
x1(99,1)=2.9233; x1(99,2)=16.1468;
x1(100,1)=0.1177; x1(100,2)=16.4696;


x2(1,1)=18.0917; x2(1,2)=10.7906;
x2(2,1)=20.8946; x2(2,2)=7.9302;
x2(3,1)=27.0622; x2(3,2)=11.5962;
x2(4,1)=5.5730; x2(4,2)=9.4899;
x2(5,1)=26.7574; x2(5,2)=11.8636;
x2(6,1)=16.7292; x2(6,2)=9.3432;
x2(7,1)=19.0151; x2(7,2)=12.4156;
x2(8,1)=24.3078; x2(8,2)=11.8160;
x2(9,1)=22.4947; x2(9,2)=10.3850;
x2(10,1)=9.3145; x2(10,2)=9.4613;
x2(11,1)=9.9780; x2(11,2)=10.4605;
x2(12,1)=22.4415; x2(12,2)=9.6565;
x2(13,1)=13.5368; x2(13,2)=9.8577;
x2(14,1)=9.2123; x2(14,2)=11.2597;
x2(15,1)=16.8452; x2(15,2)=8.5662;
x2(16,1)=16.6143; x2(16,2)=8.5577;
x2(17,1)=12.0322; x2(17,2)=9.9863;
x2(18,1)=12.3410; x2(18,2)=8.6917;
x2(19,1)=20.7440; x2(19,2)=7.3164;
x2(20,1)=7.5386; x2(20,2)=11.5666;
x2(21,1)=26.8886; x2(21,2)=9.5218;
x2(22,1)=22.9919; x2(22,2)=7.1799;
x2(23,1)=17.3493; x2(23,2)=9.2897;
x2(24,1)=18.8619; x2(24,2)=9.7411;
x2(25,1)=13.5521; x2(25,2)=5.7984;
x2(26,1)=12.7381; x2(26,2)=7.5564;
x2(27,1)=21.2411; x2(27,2)=6.7318;
x2(28,1)=24.4092; x2(28,2)=11.9952;
x2(29,1)=26.6712; x2(29,2)=9.7442;
x2(30,1)=18.2293; x2(30,2)=12.0030;
x2(31,1)=22.6769; x2(31,2)=11.0067;
x2(32,1)=30.7391; x2(32,2)=10.8611;
x2(33,1)=32.5980; x2(33,2)=8.7771;
x2(34,1)=16.7562; x2(34,2)=10.5129;
x2(35,1)=32.4282; x2(35,2)=7.6007;
x2(36,1)=18.0425; x2(36,2)=8.6968;
x2(37,1)=14.7803; x2(37,2)=8.8055;
x2(38,1)=22.3809; x2(38,2)=10.3572;
x2(39,1)=18.6982; x2(39,2)=6.7692;
x2(40,1)=25.9816; x2(40,2)=7.5022;
x2(41,1)=22.9529; x2(41,2)=10.3560;
x2(42,1)=9.6995; x2(42,2)=12.6448;
x2(43,1)=0.8253; x2(43,2)=10.6597;
x2(44,1)=22.2435; x2(44,2)=9.2030;
x2(45,1)=12.9460; x2(45,2)=9.9126;
x2(46,1)=24.4483; x2(46,2)=10.1399;
x2(47,1)=28.4938; x2(47,2)=13.5242;
x2(48,1)=13.1255; x2(48,2)=9.8689;
x2(49,1)=25.0474; x2(49,2)=11.6899;
x2(50,1)=19.9509; x2(50,2)=9.8567;
x2(51,1)=15.4784; x2(51,2)=8.5583;
x2(52,1)=28.4445; x2(52,2)=10.5570;
x2(53,1)=15.9001; x2(53,2)=11.8933;
x2(54,1)=26.3668; x2(54,2)=7.0044;
x2(55,1)=28.5033; x2(55,2)=8.0366;
x2(56,1)=6.4663; x2(56,2)=9.6549;
x2(57,1)=36.6973; x2(57,2)=10.8341;
x2(58,1)=27.1367; x2(58,2)=8.3365;
x2(59,1)=25.3004; x2(59,2)=8.8306;
x2(60,1)=14.3970; x2(60,2)=11.1212;
x2(61,1)=17.4541; x2(61,2)=7.5968;
x2(62,1)=10.7689; x2(62,2)=11.8858;
x2(63,1)=11.3941; x2(63,2)=11.5540;
x2(64,1)=13.6303; x2(64,2)=8.8437;
x2(65,1)=22.5345; x2(65,2)=11.1880;
x2(66,1)=30.0558; x2(66,2)=12.3294;
x2(67,1)=27.0878; x2(67,2)=10.5662;
x2(68,1)=16.0525; x2(68,2)=11.5472;
x2(69,1)=5.9346; x2(69,2)=10.2316;
x2(70,1)=20.4220; x2(70,2)=10.0298;
x2(71,1)=14.5875; x2(71,2)=9.8690;
x2(72,1)=8.3235; x2(72,2)=8.2102;
x2(73,1)=21.6882; x2(73,2)=8.6548;
x2(74,1)=22.4873; x2(74,2)=9.9445;
x2(75,1)=24.1396; x2(75,2)=11.7790;
x2(76,1)=17.4024; x2(76,2)=8.9218;
x2(77,1)=16.4952; x2(77,2)=12.6580;
x2(78,1)=17.7652; x2(78,2)=8.5352;
x2(79,1)=17.9541; x2(79,2)=11.1611;
x2(80,1)=20.3055; x2(80,2)=7.6421;
x2(81,1)=29.1058; x2(81,2)=8.4386;
x2(82,1)=19.7172; x2(82,2)=9.1464;
x2(83,1)=22.9040; x2(83,2)=10.4004;
x2(84,1)=31.4804; x2(84,2)=13.3684;
x2(85,1)=16.9647; x2(85,2)=10.1980;
x2(86,1)=23.3807; x2(86,2)=9.2596;
x2(87,1)=20.7638; x2(87,2)=10.9202;
x2(88,1)=13.9053; x2(88,2)=7.8416;
x2(89,1)=3.7968; x2(89,2)=9.8266;
x2(90,1)=4.7264; x2(90,2)=10.8025;
x2(91,1)=16.9223; x2(91,2)=9.4675;
x2(92,1)=15.3344; x2(92,2)=9.6976;
x2(93,1)=16.5509; x2(93,2)=10.3232;
x2(94,1)=10.5063; x2(94,2)=10.7580;
x2(95,1)=20.3627; x2(95,2)=13.0198;
x2(96,1)=20.2478; x2(96,2)=13.5945;
x2(97,1)=14.6817; x2(97,2)=11.2545;
x2(98,1)=23.5466; x2(98,2)=7.5405;
x2(99,1)=9.1412; x2(99,2)=10.6535;
x2(100,1)=12.6591; x2(100,2)=8.9793;
x2(101,1)=17.4900; x2(101,2)=8.1205;
x2(102,1)=11.2094; x2(102,2)=11.2549;
x2(103,1)=26.3995; x2(103,2)=10.8712;
x2(104,1)=27.9677; x2(104,2)=9.6512;
x2(105,1)=20.4156; x2(105,2)=9.3194;
x2(106,1)=10.5738; x2(106,2)=8.4045;
x2(107,1)=22.9527; x2(107,2)=6.6876;
x2(108,1)=3.6603; x2(108,2)=7.0692;
x2(109,1)=15.7903; x2(109,2)=7.8872;
x2(110,1)=38.2151; x2(110,2)=9.2523;


for i=1:100 r1(i)=x1(i,1);end;
for i=1:100 r2(i)=x1(i,2);end;
for i=1:110 r3(i)=x2(i,1);end;
for i=1:110 r4(i)=x2(i,2);end;

figure(1);
plot(r1,r2,'*',r3,r4,'o');
hold on

g=0;h=0;x=1;
g(x)=100;h(x)=110;
flag=1;
z=1;
nr11=0;nt11=0;nr12=0;nt12=0;nr21=0;nt21=0;nr22=0;nt22=0;

index_1=0;
index_2=0;
index_3=0;
index_4=0;
while(flag<100)                          %两分剪辑，flag为结束标志
a=randint(1,1,[1,round(g(x)/2)]);        %随机产生考试集合参考集
b=randint(1,1,[1,round(h(x)/2)]);    
c=randperm(g(x));
d=randperm(h(x));
index_1=c(1:a);
index_2=c(a+1:g(x));
index_3=d(1:b);
index_4=d(b+1:h(x));
nt11=r1(index_1);
nt12=r2(index_1);
nr11=r1(index_2);
nr12=r2(index_2);
nt21=r3(index_3);
nt22=r4(index_3);
nr21=r3(index_4);
nr22=r4(index_4);
[m,n]=size(nt11);                      %求考试集和参考集的大小
[p,q]=size(nt21);
[l,k]=size(nr11);
[u,v]=size(nr21);
m=1;p=1;l=1;u=1;
t1=0;t2=0;
w=1;
for i=1:n                              %用参考集对考试集中第一类进行判断
    for j=1:k
        t1(j)=((nt12(i)-nr12(j))^2+(nt11(i)-nr11(j))^2)^(1/2);
        if(j>=2&&t1(j-1)<t1(j))
            t1(j)=t1(j-1);
        end
    end
    for t=1:v
        t2(t)=((nt12(i)-nr22(t))^2+(nt11(i)-nr21(t))^2)^(1/2);
        if(t>=2&&t2(t-1)<t2(t))
            t2(t)=t2(t-1);
        end
    end
    if(t1(j)<=t2(t))              %正确划分的留下
        r1(w)=nt11(i);
        r2(w)=nt12(i);
        w=w+1;
    end
end

for f=1:k                        %将参考集中的第一类放入新的样本  
    r1(w)=nr11(f);
    r2(w)=nr12(f);
    w=w+1;
end


for m=1:q                       %用参考集对考试集中第二类进行判断
    for p=1:k
        t3(p)=((nt22(m)-nr12(p))^2+(nt21(m)-nr11(p))^2)^(1/2);
        if(p>=2&&t3(p-1)<t3(p))
            t3(p)=t3(p-1);
        end
    end
    for l=1:v
        t4(l)=((nt22(m)-nr22(l))^2+(nt21(m)-nr21(l))^2)^(1/2);
        if(l>=2&&t4(l-1)<t4(l))
            t4(l)=t4(l-1);
        end
    end
    if(t4(l)<=t3(p))                  %正确划分的留下
        r3(u)=nt21(m);
        r4(u)=nt22(m);
        u=u+1;
    end
end 

for e=1:v                            %将参考集中的第二类放入新的样本
    r3(u)=nr21(e);
    r4(u)=nr22(e);
    u=u+1;
end

x=x+1;z=z+1;w=w-1;u=u-1;
t5=0;t6=0;
g(x)=w;h(x)=u;
t5=g(x-1)+h(x-1);
t6=g(x)+h(x);

if(t5==t6)                         %当两次剪辑样本不变时结束标志加1
    flag=flag+1;
end

end       
       
figure(2)                          %剪辑的最终结果
plot(r1,r2,'*',r3,r4,'o');
axis([0,60,0,20]) 


i=1;j=1;                           %压缩
t7=0;t8=0;t9=0;t10=0; 
grab1=0;grab2=0; 
flag1=0;                           %标号数组                                                  
a=randint(1,1,[1,w]);              %随机抽取分别属于两类的两个样本放入store袋中
b=randint(1,1,[1,u]);
store11(i)=r1(a);
store12(i)=r2(a);
store21(j)=r3(b);
store22(j)=r4(b);
index_5=randperm(w+u);             %将剪辑后的结果混合放进grab袋中
for n=index_5(1:w)
    grab1(n)=r1(i);
    grab2(n)=r2(i);
    flag1(n)=1;                    %第一类标号为1
    i=i+1;
end
for m=index_5(w+1:u+w)
    grab1(m)=r3(j);
    grab2(m)=r4(j);
    flag1(m)=0;                    %第二类标号为0
    j=j+1;
end
i=1;j=1;

for t=1:w+u                       %用store中的样本对grab中随机抽取的样本进行判断
  if(flag1(t)==1)                 %第一类的情况   
      for m=1:i
        t7(m)=((grab2(t)-store12(m))^2+(grab1(t)-store11(m))^2)^(1/2);
        if(m>=2&&t7(m-1)<t7(m))
            t7(m)=t7(m-1);
        end
      end
      for p=1:j
        t8(p)=((grab2(t)-store22(p))^2+(grab1(t)-store21(p))^2)^(1/2);
        if(p>=2&&t8(p-1)<t8(p))
            t8(p)=t8(p-1);
        end
      end
      if(t7(m)>=t8(p))            %错判的留下
        i=i+1;
        store11(i)=grab1(t);
        store12(i)=grab2(t);

      end
  else                           %第二类情况
      for k=1:i
        t9(k)=((grab2(t)-store12(k))^2+(grab1(t)-store11(k))^2)^(1/2);
        if(k>=2&&t9(k-1)<t9(k))
            t9(k)=t9(k-1);
        end
      end
      for l=1:j
        t10(l)=((grab2(t)-store22(l))^2+(grab1(t)-store21(l))^2)^(1/2);
        if(l>=2&&t10(l-1)<t10(l))
            t10(l)=t10(l-1);
        end
      end
      if(t10(l)>=t9(k))          %错判的留下
        j=j+1;
        store21(j)=grab1(t);
        store22(j)=grab2(t);
      end
  end
end
figure(3)                        %store中的即为最终结果
plot(store11,store12,'*',store21,store22,'o');
axis([0,60,0,20])