bench_bank32nh.m

本人编的一个程序

clear
load blfbank32nh
load bank32nhindex;
param.kernel = 'rbf'; param.kernelparam = 0.02; param.nsv = 150;
param.lamda = 1e-2; param.trace = 1;
trainnum = 5000;

index1 = bank32nhindex(1,1:trainnum);
index2 = bank32nhindex(1,trainnum+1:end);
samples = bank32nhsamples(index1,:);
labels = bank32nhlabels(index1,:);
testsamples = bank32nhsamples(index2,:);
testlabels = bank32nhlabels(index2,:);

[samples,A1,B1]   =  scaletrain(samples);
[labels,A2,B2]    =  scaletrain(labels);
mul               =  1/A2;
testsamples       =  scaletest(testsamples,A1,B1);
testlabels        =  scaletest(testlabels,A2,B2);

for i = 1:2
    param.epsilon = 0.02; param.delta = 0.3;
    [alpha,sv,beta,time,obj1] = ssvr_srstrain(samples,labels,param);
    time1(:,i) = time;
    for j = 1:param.nsv
        alpha = beta(1:j,j);
        result = ssvr_srstest(testsamples,samples,param,alpha,sv(1:j));
        testerror1(j,i) = sqrt(mean((mul*(result - testlabels)).^2));
    end
    
    param.epsilon = 0.02; param.delta = 0.3;
    [alpha,sv,beta,time,centers,obj2] = ssvr_crstrain(samples,labels,param);
    time2(:,i) = time;
    for j = 1:param.nsv
        alpha = beta(1:j,j);
        result = ssvr_crstest(testsamples,centers,param,alpha,sv(1:j));
        testerror2(j,i) = sqrt(mean((mul*(result - testlabels)).^2));
    end
end

plot(mean(testerror1,2),'b:');
hold on
plot(mean(testerror2,2),'r');
legend('SSVR-SRS','SSVR-CRS');
xlabel('Number of basis functions');
ylabel('Test error');