SVM分类与回归的

在matlab使用SVM工具箱进行分类与回归的小例子

%% Matlab神经网络43个案例分析

% 初识SVM分类与回归
% by 李洋（faruto）
% http://www.matlabsky.com
% Email:faruto@163.com
% http://weibo.com/faruto 
% http://blog.sina.com.cn/faruto
% 2013.01.01
%% A Little Clean Work
tic;
clear;
clc;
close all;
format compact;
%% 使用Libsvm进行分类的小例子
%{
一个班级里面有两个男生（男生1、男生2），两个女生（女生1、女生2），其中
男生1 身高：176cm 体重：70kg；
男生2 身高：180cm 体重：80kg；

女生1 身高：161cm 体重：45kg；
女生2 身高：163cm 体重：47kg；
如果我们将男生定义为1，女生定义为-1，并将上面的数据放入矩阵data中
并在label中存入男女生类别标签（1、-1）
%}
data = [176 70;
        180 80;
        161 45;
        163 47];
label = [1;1;-1;-1];
%{
这样上面的data矩阵就是一个属性矩阵，行数4代表有4个样本，列数2表示属性有两个，
label就是标签（1、-1表示有两个类别：男生、女生）。
%}

% 利用libsvm建立分类模型
model = svmtrain（labeldata）;

%{
此时该班级又转来一个新学生，其
身高190cm，体重85kg
我们想给出其标签（想知道其是男[1]还是女[-1]）
由于其标签我们不知道，我们假设其标签为-1（也可以假设为1） 
%}
testdata = [190 85];
testdatalabel = -1;

[predictlabelaccuracy] = svmpredict（testdatalabeltestdatamodel）;
predictlabel
if 1 == predictlabel
    disp（‘==该生为男生‘）;
end
if -1 == predictlabel
    disp（‘==该生为女生‘）;
end

% % 下面再使用libsvm工具箱本身带的测试数据heart_scale来进行一下测试
% 首先载入数据
load heart_scale;
data = heart_scale_inst;
label = heart_scale_label;

% 选取前200个数据作为训练集合，后70个数据作为测试集合
ind = 200;
traindata = data（1:ind:）;
trainlabel = label（1:ind:）;
testdata = data（ind+1:end:）;
testlabel = label（ind+1:end:）;

% 利用训练集合建立分类模型
model = svmtrain（trainlabeltraindata‘-s 0 -t 2 -c 1.2 -g 2.8‘）;

% 利用建立的模型看其在训练集合上的分类效果
[ptrainacctrain] = svmpredict（trainlabeltraindatamodel）;

% 预测测试集合标签
[ptestacctest] = svmpredict（testlabeltestdatamodel）;
%% 使用Libsvm进行回归的小例子
% 生成待回归的数据
x = （-1:0.1:1）‘;
y = -x.^2;

% 建模回归模型
model = svmtrain（yx‘-s 3 -t 2 -c 2.2 -g 2.8 -p 0.01‘）;

% 利用建立的模型看其在训练集合上的回归效果
[pymse] = svmpredict（yxmodel）;

scrsz = get（0‘ScreenSize‘）;
figure（‘Position‘[scrsz（3）*1/4 scrsz（4）*1/6  scrsz（3）*4/5 scrsz（4）]*3/4）;
plot（xy‘o‘）;
hold on;
plot（xpy‘r*‘）;
legend（‘原始数据‘‘回归数据‘）;
grid on;

% 进行预测
testx = 1.1;
display（‘真实数据‘）
testy = -testx.^2

[ptestytmse] = svmpredict（testytestxmodel）;
display（‘预测数据‘）;
ptesty

%% Record Time
toc;

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

     文件       2555  2013-08-18 10:51  Chapter_ClassifyRegressUsingLibsvm.m

     文件      28904  2005-03-22 20:37  heart_scale.mat

----------- ---------  ---------- -----  ----

                31459                    2