SVM对偶空间求解与直接求解效率比较

.m 资源 本代码利用对偶空间求解SVM模型与直接求解进行对比，实现简单的模式分类

clc;
clear all;
close all;

img=imread（‘1.jpg‘）;

img=double（img）;


figure
imshow（uint8（img））;
title（‘左键单击海洋部分，中键点击陆地部分，右键单击取天空部分‘）
[xxyybotton] = ginput（）;
xx=floor（xx）;
yy=floor（yy）;

%%%%%%
img0=rgb2hsv（img）;
img（::1）=img0（::2）;
img（::2）=img0（::3）;
%%%%%%
count=size（xx）;
for cou= 1:count
        trnx1（cou:）=[img（yy（cou）xx（cou）1）img（yy（cou）xx（cou）2）];
        if （botton（cou）==1）
            trny1（cou） = 1;
        else
            trny1（cou） = -1;
        end
end
trny1=trny1‘;

i=1;
for cou= 1:count
    if （botton（cou）~=1）
        trnx2（i:）=[img（yy（cou）xx（cou）1）img（yy（cou）xx（cou）2）];
        if （botton（cou）==2）
            trny2（i） = 1;
        end
        if（botton（cou）==3）
            trny2（i） = -1;
        end
        i=i+1;
    end
end
trny2=trny2‘;
figure
ker=‘linear‘;
[nsv1 alpha1 bias1] = svc（trnx1trny1ker10）;
svcplot（trnx1trny1keralpha1bias1）;
hold off
figure
[nsv2 alpha2 bias2] = svc（trnx2trny2ker10）;
svcplot（trnx2trny2keralpha2bias2）;
hold off
% 
% [MNk]=size（img）;
% img1=img;
% tic
% for i =1:M
%     for j= 1:N
%         testx=img（ij1:2）;
%         predicty=svcoutput（trnx1trny1testxkeralpha1bias1）;
%         
%         if（predicty==1）;
%             img1（ij:）=[25500];
%         else
%             predicty=svcoutput（trnx2trny2testxkeralpha2bias2）;
%             if（predicty==1）
%                 img1（ij:）=[255255255];
%             else
%                 img1（ij:）=[00255];
%             end
%         end
%     end
% end




% figure
% imshow（img1）;

testx=img（111:2）;
predicty0=svcoutput（trnx1trny1testxkeralpha1bias1）


for j=10:20
    imax=zeros（110）
    for i=1:10
         testx=img（floor（i*662/10）j1:2）;
         predicty=svcoutput（trnx1trny1testxkeralpha1bias1）
         if（predicty~=predicty0）
             imax（i）=0
         elseif（predicty==predicty0）
             imax（i）=1
         end
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2068  2014-04-06 10:13  svm-gunn\.m

     文件     173188  2014-04-08 01:36  svm-gunn\2 （3）.jpg

     文件        371  2012-03-06 00:30  svm-gunn\binomial.m

     文件        144  2012-03-06 00:30  svm-gunn\centrefig.m

     文件       1728  2012-03-06 00:30  svm-gunn\cmap.mat

     文件       1105  2012-03-06 00:30  svm-gunn\Contents.m

     文件       2808  2014-04-08 10:48  svm-gunn\line.c

     文件       4020  2014-04-08 10:49  svm-gunn\line.mexw32

     文件       1314  2014-04-08 11:21  svm-gunn\my_compare_svm.asv

     文件       1445  2014-04-08 11:40  svm-gunn\my_compare_svm.m

     文件        457  2012-03-06 00:30  svm-gunn\nobias.m

     文件      49152  2012-03-06 00:30  svm-gunn\qp.dll

     文件       2642  2012-03-06 00:30  svm-gunn\README

     文件       5113  2012-03-06 00:30  svm-gunn\road1.mat

     文件       5400  2012-03-06 00:30  svm-gunn\road2.mat

     文件        312  2012-03-06 00:30  svm-gunn\softmargin.m

     文件       2771  2012-03-06 00:30  svm-gunn\svc.m

     文件        837  2012-03-06 00:30  svm-gunn\svcerror.m

     文件       1228  2012-03-06 00:30  svm-gunn\svcinfo.m

     文件        973  2012-03-06 00:30  svm-gunn\svcoutput.m

     文件       3109  2012-03-06 00:30  svm-gunn\svcplot.m

     文件       1299  2012-03-06 00:30  svm-gunn\svdatanorm.m

     文件       2636  2014-03-19 17:13  svm-gunn\svkernel.m

     文件       2637  2012-03-06 00:30  svm-gunn\svkernel.m.bak

     文件        401  2012-03-06 00:30  svm-gunn\svtol.m

     文件     233118  2014-03-19 16:02  svm-gunn\TEST001.BMP

     文件     233118  2014-03-19 16:03  svm-gunn\TEST002.BMP

     文件       5386  2012-03-06 00:30  svm-gunn\uiclass.m

     文件      12592  2012-03-06 00:30  svm-gunn\uiclass.mat

     文件       5627  2012-03-06 00:30  svm-gunn\uiregress.m

............此处省略5个文件信息