【机器学习】libsvm 简单使用示例（C++）
libsvm简单使用demo
一、libsvm使用说明
二、svm.h源码
#ifndef _LIBSVM_H //如果没有定义 _LIBSVM_H 宏
#define _LIBSVM_H //则定义 _LIBSVM_H 宏，用于防止重复包含
#define LIBSVM_VERSION 317 //定义一个宏，表示 libsvm 的版本号
#ifdef __cplusplus //如果是 C++ 编译器
extern "C" { //则使用 C 语言的链接方式
#endif
extern int libsvm_version; //声明一个外部变量，表示 libsvm 的版本号
struct svm_node //定义一个结构体，表示一个特征节点
{
int index; //特征的索引
double value; //特征的值
};
struct svm_problem //定义一个结构体，表示一个 SVM 问题
{
int l; //训练集的行数
double *y; //训练集的标签向量
struct svm_node **x; //训练集的特征矩阵
};
enum { C_SVC, NU_SVC, ONE_CLASS, EPSILON_SVR, NU_SVR }; /* svm_type */ //定义一个枚举类型，表示 SVM 的类型
enum { LINEAR, POLY, RBF, SIGMOID, PRECOMPUTED }; /* kernel_type */ //定义一个枚举类型，表示核函数的类型
struct svm_parameter //定义一个结构体，表示 SVM 的参数
{
int svm_type; //SVM 的类型
int kernel_type; //核函数的类型
int degree; /* for poly */ //多项式核函数的次数
double gamma; /* for poly/rbf/sigmoid */ //核函数的系数
double coef0; /* for poly/sigmoid */ //核函数的常数项
/* these are for training only */ //以下是只用于训练的参数
double cache_size; /* in MB */ //内核缓存的大小（以 MB 为单位）
double eps; /* stopping criteria */ //停止准则的容差
double C; /* for C_SVC, EPSILON_SVR and NU_SVR */ //C-SVM 分类器或回归器的惩罚系数
int nr_weight; /* for C_SVC */ //不同类别的权重的个数
int *weight_label; /* for C_SVC */ //不同类别的权重的标签
double* weight; /* for C_SVC */ //不同类别的权重的值
double nu; /* for NU_SVC, ONE_CLASS, and NU_SVR */ //nu-SVM 分类器或回归器的参数
double p; /* for EPSILON_SVR */ //epsilon-SVR 回归器的损失函数的参数
int shrinking; /* use the shrinking heuristics */ //是否使用启发式收缩
int probability; /* do probability estimates */ //是否计算概率估计
};
//
// svm_model
//
struct svm_model //定义一个结构体，表示 SVM 的模型
{
struct svm_parameter param; /* parameter */ //SVM 的参数
int nr_class; /* number of classes, = 2 in regression/one class svm */ //类别的个数，回归或单类 SVM 为 2
int l; /* total #SV */ //支持向量的总数
struct svm_node **SV; /* SVs (SV[l]) */ //支持向量的矩阵
double **sv_coef; /* coefficients for SVs in decision functions (sv_coef[k-1][l]) */ //支持向量在决策函数中的系数
double *rho; /* constants in decision functions (rho[k*(k-1)/2]) */ //决策函数中的常数
double *probA; /* pariwise probability information */ //成对概率信息
double *probB;
int *sv_indices; /* sv_indices[0,...,nSV-1] are values in [1,...,num_traning_data] to indicate SVs in the training set */ //支持向量在训练集中的索引
/* for classification only */ //以下是只用于分类的信息
int *label; /* label of each class (label[k]) */ //每个类别的标签
int *nSV; /* number of SVs for each class (nSV[k]) */ //每个类别的支持向量的个数
/* nSV[0] + nSV[1] + ... + nSV[k-1] = l */
/* XXX */
int free_sv; /* 1 if svm_model is created by svm_load_model*/
/* 0 if svm_model is created by svm_train */ //表示 svm_model 是由 svm_load_model 还是 svm_train 创建的
};
struct svm_model *svm_train(const struct svm_problem *prob, const struct svm_parameter *param); //声明一个函数，用于训练 SVM 模型，参数为 SVM 问题和 SVM 参数，返回值为 SVM 模型的指针
void svm_cross_validation(const struct svm_problem *prob, const struct svm_parameter *param, int nr_fold, double *target); //声明一个函数，用于进行交叉验证，参数为 SVM 问题、SVM 参数、交叉验证的折数和目标向量，无返回值
int svm_save_model(const char *model_file_name, const struct svm_model *model); //声明一个函数，用于保存 SVM 模型到文件，参数为文件名和 SVM 模型的指针，返回值为 0 表示成功，-1 表示失败
struct svm_model *svm_load_model(const char *model_file_name); //声明一个函数，用于从文件加载 SVM 模型，参数为文件名，返回值为 SVM 模型的指针，如果失败则为 NULL
int svm_get_svm_type(const struct svm_model *model); //声明一个函数，用于获取 SVM 的类型，参数为 SVM 模型的指针，返回值为枚举类型的值
int svm_get_nr_class(const struct svm_model *model); //声明一个函数，用于获取类别的个数，参数为 SVM 模型的指针，返回值为整数
void svm_get_labels(const struct svm_model *model, int *label); //声明一个函数，用于获取每个类别的标签，参数为 SVM 模型的指针和标签向量，无返回值
void svm_get_sv_indices(const struct svm_model *model, int *sv_indices); //声明一个函数，用于获取支持向量在训练集中的索引，参数为 SVM 模型的指针和索引向量，无返回值
int svm_get_nr_sv(const struct svm_model *model); //声明一个函数，用于获取支持向量的个数，参数为 SVM 模型的指针，返回值为整数
double svm_get_svr_probability(const struct svm_model *model); //声明一个函数，用于获取 SVR 的概率估计，参数为 SVM 模型的指针，返回值为浮点数，如果失败则为 0
double svm_predict_values(const struct svm_model *model, const struct svm_node *x, double* dec_values); //声明一个函数，用于预测一个样本的决策值，参数为 SVM 模型的指针、样本的特征向量和决策值向量，返回值为预测的标签
double svm_predict(const struct svm_model *model, const struct svm_node *x); //声明一个函数，用于预测一个样本的标签，参数为 SVM 模型的指针和样本的特征向量，返回值为预测的标签
double svm_predict_probability(const struct svm_model *model, const struct svm_node *x, double* prob_estimates); //声明一个函数，用于预测一个样本的概率估计，参数为 SVM 模型的指针、样本的特征向量和概率估计向量，返回值为预测的标签
void svm_free_model_content(struct svm_model *model_ptr); //声明一个函数，用于释放 SVM 模型的内容，参数为 SVM 模型的指针，无返回值
void svm_free_and_destroy_model(struct svm_model **model_ptr_ptr); //声明一个函数，用于释放并销毁 SVM 模型，参数为 SVM 模型的指针的指针，无返回值
void svm_destroy_param(struct svm_parameter *param); //声明一个函数，用于销毁 SVM 参数，参数为 SVM 参数的指针，无返回值
const char *svm_check_parameter(const struct svm_problem *prob, const struct svm_parameter *param); //声明一个函数，用于检查 SVM 参数是否合法，参数为 SVM 问题和 SVM 参数的指针，返回值为一个字符串，如果为 NULL，则表示参数合法，否则表示参数有误
int svm_check_probability_model(const struct svm_model *model); //声明一个函数，用于检查 SVM 模型是否支持概率估计，参数为 SVM 模型的指针，返回值为一个整数，如果为 0，则表示不支持，否则表示支持
void svm_set_print_string_function(void (*print_func)(const char *)); //声明一个函数，用于设置 libsvm 的输出函数，参数为一个函数指针，该函数接受一个字符串作为参数，无返回值
#ifdef __cplusplus //如果是 C++ 编译器
}
#endif //则结束 C 语言的链接方式
#endif /* _LIBSVM_H */ //结束防止重复包含的条件三、示例demo
#include "svm.h" //引入 libsvm 的头文件
#include //引入字符处理的头文件
#include //引入标准库的头文件
#include //引入向量容器的头文件
#include //引入输入输出流的头文件
using namespace std; //使用标准命名空间
#define Malloc(type,n) (type *)malloc((n)*sizeof(type)) //定义一个宏，用于分配内存
struct svm_parameter param; // 由 parse_command_line 函数设置，用于存储 SVM 的参数
struct svm_problem prob; // 由 read_problem 函数设置，用于存储训练集的特征和标签
struct svm_model *model; //用于存储训练后的 SVM 模型
struct svm_node *x_space; //用于存储训练集的特征向量
vector generateLabels(int labelsSize) //定义一个函数，用于生成标签向量
{
std::vector labels; //创建一个 int 类型的向量
for (int i=0; i > generateData(int problemSize, int featureNum) //定义一个函数，用于生成特征数据
{
std::vector > data; //创建一个 double 类型的二维向量
for (int i=0; i featureSet; //创建一个 double 类型的向量，用于存储一行特征
for (int j=0; j > data = generateData(sizeOfProblem,elements); //调用 generateData 函数，生成特征数据
vector labels = generateLabels(sizeOfProblem); //调用 generateLabels 函数，生成标签数据
cout<<"data size = "< 输出结果：
feature pushed(第1行)
……
feature pushed(第300行)
data size = 30
labels size = 30
prob.y[0] = 1
prob.y[1] = 2
prob.y[2] = 1
prob.y[3] = 2
prob.y[4] = 1
prob.y[5] = 2
prob.y[6] = 1
prob.y[7] = 2
prob.y[8] = 1
prob.y[9] = 2
prob.y[10] = 1
prob.y[11] = 2
prob.y[12] = 1
prob.y[13] = 2
prob.y[14] = 1
prob.y[15] = 2
prob.y[16] = 1
prob.y[17] = 2
prob.y[18] = 1
prob.y[19] = 2
prob.y[20] = 1
prob.y[21] = 2
prob.y[22] = 1
prob.y[23] = 2
prob.y[24] = 1
prob.y[25] = 2
prob.y[26] = 1
prob.y[27] = 2
prob.y[28] = 1
prob.y[29] = 2
x_space[0].index = 1（第1个数据）
x_space[0].value = 0
x_space[1].index = 2
x_space[1].value = 1
x_space[2].index = 3
x_space[2].value = 2
x_space[3].index = 4
x_space[3].value = 3
x_space[4].index = 5
x_space[4].value = 4
x_space[5].index = 6
x_space[5].value = 5
x_space[6].index = 7
x_space[6].value = 6
x_space[7].index = 8
x_space[7].value = 7
x_space[8].index = 9
x_space[8].value = 8
x_space[9].index = 10
x_space[9].value = 9
x_space[10].index = -1
x_space[10].value = 0
（第2个数据）
x_space[11].index = 1
x_space[11].value = 0
x_space[12].index = 2
x_space[12].value = 1
x_space[13].index = 3
x_space[13].value = 2
x_space[14].index = 4
x_space[14].value = 3
x_space[15].index = 5
x_space[15].value = 4
x_space[16].index = 6
x_space[16].value = 5
x_space[17].index = 7
x_space[17].value = 6
x_space[18].index = 8
x_space[18].value = 7
x_space[19].index = 9
x_space[19].value = 8
x_space[20].index = 10
x_space[20].value = 9
x_space[21].index = -1
x_space[21].value = 0
……
（第30个数据）
x_space[319].index = 1
x_space[319].value = 0
x_space[320].index = 2
x_space[320].value = 1
x_space[321].index = 3
x_space[321].value = 2
x_space[322].index = 4
x_space[322].value = 3
x_space[323].index = 5
x_space[323].value = 4
x_space[324].index = 6
x_space[324].value = 5
x_space[325].index = 7
x_space[325].value = 6
x_space[326].index = 8
x_space[326].value = 7
x_space[327].index = 9
x_space[327].value = 8
x_space[328].index = 10
x_space[328].value = 9
x_space[329].index = -1
x_space[329].value = 0
line 0
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 1
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 2
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 3
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 4
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 5
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 6
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 7
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 8
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 9
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 10
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 11
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 12
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 13
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 14
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 15
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 16
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 17
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 18
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 19
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 20
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 21
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 22
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 23
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 24
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 25
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 26
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 27
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 28
1---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
line 29
2---1:0 2:1 3:2 4:3 5:4 6:5 7:6 8:7 9:8 10:9
all ok
*
optimization finished, #iter = 15
nu = 1.000000
obj = -30.000000, rho = 0.000000
nSV = 30, nBSV = 30
Total nSV = 30参考网址
https://github.com/cjlin1/libsvm
https://github.com/niosus/SVM_example