testrbf.m

来自「各种SVM分类算法」· M 代码 · 共 21 行

               clear;clc;close all;
               disp('-------------采用径向基内积函数的支持向量机的应用-------------');
               disp('输入样本矩阵(每行表示一个测试数据):');
               x=[randn(1,20)-3 randn(1,20)+3;randn(1,40)]'
               disp('样本所属类别:');
               y=[-ones(1,20) ones(1,20)]'
               [nsv,alpha,bias,T]=svm168(x,y,'rbf',10,4);
               fprintf('运行时间T=%g\n',T);
               fprintf('支持向量个数：%g\n',nsv);

               svcplot(x,y,'linear',alpha,bias);

               disp('测试数据:');
               X_test=[randn(1,20)-3 randn(1,20)+3;randn(1,40)]'
               disp('测试数据类别:');
               Y_test=[-ones(1,20) ones(1,20)]'
               disp('测试数据类别判决:');
               Output=kernel('rbf',X_test,x(find(alpha),:),4)*alpha(find(alpha))+bias
               disp('错误率:');
               number_error=sum(abs(Y_test-sign(Output)))/2;
               error=number_error/20