cardiag.asv

标准svm3.0程序 内有mex优化程序

clear all
close all
load data\car.txt;
load discar;
%if(0)
%end
value1=discar(:,13);
standdata=discar(:,2:12);
[Samples1,Labels,testS1,testL]=createset(value1,standdata);
standdata=[Samples1,Labels];
[row column]=size(standdata);
for i=1:column
    r=[];
    for j=1:row      
       a=standdata(j,i);
       x=standdata(:,i);
        for m=1:length(x)
            r(j,m)=dissim(a,x(m));
        end
 end
    eval(['ssr' num2str(i) '=r;']);
 end

%%%%%reduct%%%%%%%%%%%
n=[];
sig=[];
x=0;
base=ones(row);
r=eval(['ssr' num2str(column)]);
entropyd=entropy(r);
attrinu=column-1;
for j=attrinu:-1:1
    for i=1:attrinu
       r1=eval(['ssr' num2str(i)]);
        sig(i)=entropyd+entropy(min(r1,base))-entropy(min(min(r1,r),base));
    end
    [x1,n1]=max(sig);
    x=[x;x1];
    len=length(x);
    if abs(x(len)-x(len-1))>0.001
        base1=eval(['ssr' num2str(n1)]);
        base=min(base,base1);
        n=[n;n1];
    else
        break
    end
end
n
pause
%%%%%%%%%%%%%%%%分割训练集和测试集%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
value1=car(:,13);
standdata=car(:,2:12);
for i=1:11
    maxd=max(standdata(:,i));
    mind=min(standdata(:,i));
    standdata(:,i)=(standdata(:,i)-mind)/(maxd-mind);
end
[Samples1,Labels,testS1,testL]=createset(value1,standdata(:,n));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Degree=3;
C=1;
Gamma=1;
Coeff=1;
%Degree ---(default: 3)
%                 Gamma  ---(default: 1)
%                 Coeff  ---(default: 1)
%     C: Cost of the constrain violation  (default: 1)
[AlphaY1, SVs1, Bias1, Parameters1, nSV1, nLabel1] =PolySVC(Samples1', Labels'%,Degree, C, Gamma, Coeff);
[ClassRate1, DecisionValue1, Ns1, ConfMatrix1, PreLabels1]= SVMTest(testS1',testL', AlphaY1, SVs1, Bias1,Parameters1, nSV1, nLabel1);
ClassRate1