求助SVM工具箱（含程序）
我在网上下载了一个工具箱，里面有许多function函数，可是我不需要每个都去选择，只找出我需要的部分来编程。用工具箱主程序运行没问题，可是我在整理出我需要的部分后再运行会出问题。
我的程序如下：
clc
clear
C = 100;
e = 0.5;
%************* 构造训练样本 **************************
ker = struct('type','gauss','width',0.6);
n = 50;
rand('state',42);
X = linspace(-4,4,n);
Ys = (1-X+2*X.^2).*exp(-.5*X.^2);
f = 0.2;
Y = Ys+f*max(abs(Ys))*(2*rand(size(Ys))-1)/2;
figure;
plot(X,Ys,'b-',X,Y,'b*');
title('\epsilon-SVR');
hold on;
%************ 训练支持向量 ******************************
tic
options = optimset;
options.LargeScale = 'on';
options.Display = 'off';
n = length(Y);
Q = kernel(ker,X,X);
H = ;
f = ; % 符号不一样,决策函数就不一样,实际上是一回事!
A = [];
b = [];
Aeq = ;
beq = 0;
lb = zeros(2*n,1);
ub = C*ones(2*n,1);
a0 = zeros(2*n,1);
= quadprog(H,f,A,b,Aeq,beq,lb,ub,a0,options);
a = a(1:n)-a(n+1:end);
eXitflag
svm.ker = ker;
svm.x = X;
svm.y = Y;
svm.a = a';
t_train = toc
%*********** 寻找支持向量 *********************
a = svm.a;
epsilon = 1e-8;
i_sv = find(abs(a)>epsilon);
plot(X(i_sv),Y(i_sv),'ro');
%*********** 测试输出 *************************
tic
cathe = 10e+6;
nx = size(svm.x,2);
nt = size(Ys,2);
block = ceil(nx*nt/cathe);
num = ceil(nt/block);
for i = 1:block
if (i==block)
index = [(i-1)*num+1:nt];
else
index = (i-1)*num+;
end
end
tmp = a*kernel(ker,X,X(:,i_sv));
b = Y(i_sv)-tmp;
b = mean(b);
tmp = a*kernel(ker,X,Xt);
Yd = (tmp+b);
t_sim = toc
plot(X,Yd,'r--');
hold off;
调用的工具箱：(出问题的部分)
function Yd = svmSim(svm,Xt)
% ------------------------------------------------------------%
cathe = 10e+6; % kernel输出的元素个数的上限
nx = size(svm.x,2); % 训练样本数
nt = size(Xt,2); % 测试样本数
block = ceil(nx*nt/cathe); % 分块处理
num = ceil(nt/block); % 每块测试样本数
for i = 1:block
if (i==block)
index = [(i-1)*num+1:nt];
else
index = (i-1)*num+;
end
Yd(index) = svmSim_block(svm,Xt(:,index)); % 测试输出
end
% ------------------------------------------------------------%
function Yd = svmSim_block(svm,Xt)
% 输入参数:
% svm支持向量机(结构体变量)
% the following fields:
% type - 支持向量机类型{'svc_c','svc_nu','svm_one_class','svr_epsilon','svr_nu'}
% ker - 核参数
% type - linear :k(x,y) = x'*y
% poly :k(x,y) = (x'*y+c)^d
% gauss:k(x,y) = exp(-0.5*(norm(x-y)/s)^2)
% tanh :k(x,y) = tanh(g*x'*y+c)
% degree - Degree d of polynomial kernel (positive scalar).
% offset - Offset c of polynomial and tanh kernel (scalar, negative for tanh).
% width- Width s of Gauss kernel (positive scalar).
% gamma- Slope g of the tanh kernel (positive scalar).
% x - 训练样本
% y - 训练目标;
% a - 拉格朗日乘子
%
% Xt测试样本,d×n的矩阵,n为样本个数,d为样本维数
% 输出参数:
% Yd测试输出,1×n的矩阵,n为样本个数,值为+1或-1
% ------------------------------------------------------------%
type = svm.type;
ker = svm.ker;
X = svm.x;
Y = svm.y;
a = svm.a;
% ------------------------------------------------------------%
% 测试输出
epsilon = 1e-8; % 如果小于此值则认为是0
i_sv = find(abs(a)>epsilon); % 支持向量下标
switch type
case 'svc_c',
tmp = (a.*Y)*kernel(ker,X,X(:,i_sv)); % 行向量
b = Y(i_sv)-tmp;
b = mean(b);
tmp =(a.*Y)*kernel(ker,X,Xt);
tmp = tmp+b;
Yd = sign(tmp);
case 'svc_nu',
%------------------------------------%
% 与 'svc_c' 情况相同
tmp = (a.*Y)*kernel(ker,X,X(:,i_sv)); % 行向量
b = Y(i_sv)-tmp;
b = mean(b);
tmp =(a.*Y)*kernel(ker,X,Xt);
Yd = sign(tmp+b);
case 'svm_one_class',
n_sv = length(i_sv);
tmp1 = zeros(n_sv,1);
for i = 1:n_sv
index = i_sv(i);
tmp1(i) = kernel(ker,X(:,index),X(:,index));
end
tmp2 = 2*a*kernel(ker,X,X(:,i_sv)); % 行向量
tmp3 = sum(sum(a'*a.*kernel(ker,X,X)));
R_square = tmp1-tmp2'+tmp3;
R_square = mean(R_square); % 超球半径 R^2 (对所有支持向量求平均的结果)
nt = size(Xt,2); % 测试样本数
tmp4 = zeros(nt,1); % 列向量
for i = 1:nt
tmp4(i) = kernel(ker,Xt(:,i),Xt(:,i));
end
tmp5 = 2*a*kernel(ker,X,Xt); % 行向量
Yd = sign(tmp4-tmp5'+tmp3-R_square);
case 'svr_epsilon',
tmp = a*kernel(ker,X,X(:,i_sv)); % 行向量
b = Y(i_sv)-tmp; % 符号不一样,决策函数就不一样,实际上是一回事!
%b = Y(i_sv)+tmp;
b = mean(b);
tmp =a*kernel(ker,X,Xt); % 符号不一样,决策函数就不一样,实际上是一回事!
%tmp =-a*kernel(ker,X,Xt);
Yd = (tmp+b);
case 'svr_nu',
%------------------------------------%
% 与'svr_epsilon' 情况相同
tmp = a*kernel(ker,X,X(:,i_sv)); % 行向量
b = Y(i_sv)-tmp; % 符号不一样,决策函数就不一样,实际上是一回事!
%b = Y(i_sv)+tmp;
b = mean(b);
tmp =a*kernel(ker,X,Xt); % 符号不一样,决策函数就不一样,实际上是一回事!
%tmp =-a*kernel(ker,X,Xt);
Yd = (tmp+b);
otherwise,
end
问题是，用我的程序，MATLAB提示找不到Xt，我的程序确实没有定义Xt，可是工具箱里也同样没有定义Xt，为什么工具箱能运行，我这么就不行呢（整个工具箱出现Xt的地方也就只在上述程序里）
恳请大伙支招！！！谢谢！
[ 本帖最后由 ChaChing 于 2009-6-9 14:48 编辑 ] XT只是个接口
回复 沙发 tisang 的帖子
能否再详细些呢？怎么处理呢？
谢谢
你说的这个工具箱发给我一格吧
超球结构SVM工具箱,你能给我发一份吗?我研究一下。wangkai73@sina.comQQ:422574600,手机：15929327311（西安）
回复 板凳 paulliuqing 的帖子
你好，不知那个svm工具箱你那里有吗，可以发给我一份吗，谢了回复 楼主 paulliuqing 的帖子
楼主你好，请问你说的那个svm工具箱可以给我发一份吗，谢了，我QQ113304783 LZ是否可以分享一下svm的工具箱，在网上没下下来，，感激啊，，，邮箱：354773461@qq.com 谢谢楼主
页:
[1]