RVM1 RVM源代码

RVM源代码，有例子说明，供初学者使用，不错的代码

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;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         886  2011-07-20 16:53  RVM1\demo.m
     文件          85  2011-06-22 15:05  RVM1\kernel.m
     文件        2815  2011-06-22 16:51  RVM1\SeqSparBayeLear.m