我目前正在开发自己的内核用于分类，并希望将其纳入libsvm，替代libsvm提供的标准内核。

然而，我并不完全确定如何操作，显然也不希望犯任何错误。请注意，我的C++水平不是很好。我在libsvm的常见问题页面上发现了以下内容：

问：我想使用我自己的内核。有例子吗？在svm.cpp中，有两个用于内核评估的子程序：k_function()和kernel_function()。我应该修改哪一个？一个例子是“LIBSVM用于字符串数据”的LIBSVM工具。

我们有两个函数的原因如下。对于RBF内核exp(-g |xi – xj|^2)，如果我们先计算xi – xj然后计算范数平方，需要3n次操作。因此，我们考虑exp(-g (|xi|^2 – 2dot(xi,xj) +|xj|^2))，通过在开始时计算所有|xi|^2，可以将操作次数减少到2n。这是用于训练的。对于预测，我们不能这样做，因此需要一个使用3n次操作的常规子程序。最简单的添加你自己的内核的方法是，在这两个子程序中放入相同的代码，替换任何内核。

因此，我试图找到这两个子程序k_function()和kernel_function()。前者我在svm.cpp中找到了以下签名：

double Kernel::k_function(const svm_node *x, const svm_node *y,              const svm_parameter& param)

我理解正确吗，x和y分别存储我特征矩阵的一个观测值（=行）在数组中，我需要返回内核值k(x,y)？

另一方面，我完全找不到kernel_function()函数。Kernel类中有一个名为kernel_function的指针，并且有以下声明

double (Kernel::*kernel_function)(int i, int j) const; 

它在Kernel构造函数中被设置。在这种情况下，i和j是什么？我猜我需要设置这个指针，对吗？

一旦我重写了Kernel::k_function和Kernel::*kernel_function，我就完成了，libsvm将使用我的内核来比较两个观测值吗？

谢谢！

回答：

你不必为了使用自己的内核而修改LIBSVM的代码，你可以使用预计算内核选项（即，-t 4 training_set_file）。

因此，你可以根据需要在外部计算内核矩阵，将值存储在一个文件中，并将预计算的内核加载到LIBSVM中。在LIBSVM压缩包中的README文件中有如何操作的解释和示例（参见预计算内核部分的第236行）。