【正文】 1138 1315 130 1260 2398 1404 188 2615 110 177 2215 116 1461 1094 1403 1269 115 1227 1088 105 1670 41 1344 1162 1813 107 1261 1229 633 106 563 514 1104 1651 28 1076 1567 1033 1321 1030 1152 2145 499 ]。 aa=size(p)。 n=aa(1)。 sigma=2。 T=[ones(50,1)。ones(50,1)]。 for i=1:n for j=1:n K(i,j)=T(i)*p(i,:)*p(j,:)39。*T(j)。 end end f=ones(n,1) aeq=T39。 beq=0。 la=zeros(1,n)。 lb=ones(1,n)。 [x,fval]=quadprog(K,f,[],[],aeq,beq,la,lb)。 w=0。 for j=1:n w=w+x(j)*p(j,:)*T(j)。 end cc=find(x)。 zb=cc(1)。 if T(zb)0 b1=1p(zb,:)*w39。 else if T(zb)0 b1=1p(zb,:)*w39。 end end t_test=p。 for j=1:n g(j)=t_test(j,:)*w39。+b1。 end for j=1:n if g(j)0 class(j)=1。 else if g(j)0 class(j)=1。 else class(j)=0。 end end end 結(jié)果 class39。 ans = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2. clear all p=[ 520 561 99 409