关联向量机RVM程序

关联向量机程序 注释完全 步骤清晰 带有一个demo说明如何使用程序

clc;clear
X = [-10:20/150:10]‘;
Y = zeros（length（X）1）;
Y = 10*sin（X）+rand（length（X）1）;
Xt = [-10:20/160:10]‘;
Yt = 10*sin（Xt）;%+rand（length（Xt）1）;
%%
Xtrain = X;
Ytrain = Y;
Xtest = Xt;
Ytest = Yt;
%% 计算核矩阵
N = size（Xtrain1）;
delta = .3;
for i=1:N
    basis（i1） = 1;
    for j =1:N
        basis（ij+1） = kernel（Xtrain（i:）Xtrain（j:）delta）;
    end
end
%%
[mu model_idx] = SeqSparBayeLear（basisYtrain1000）;
%% test 
Nt = size（Xtest1）;
idx = find（model_idx==1）;
for i=1:Nt
    for j = 1:length（idx）
        if idx（j） == 1
            basis_t（ij） = 1;
        else
            basis_t（ij） = kernel（Xtrain（idx（j）-1:）Xtest（i:）delta）;
        end
    end
end
% 预测
y = basis_t*mu;
%%
% plot 训练结果
ytrain = basis（:idx）*mu;
figure;plot（ytrain）;hold on;plot（Ytrain‘r‘）
% plot 测试结果
figure;
plot（y）;hold on;plot（Ytest‘r‘）

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

     文件       2815  2011-06-22 16:51  SeqSparBayeLear.m

     文件        886  2011-07-20 16:53  demo.m

     文件         85  2011-06-22 15:04  kernel.m

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

                 3786                    3