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;
%%?計(jì)算核矩陣
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ù)測(cè)
y?=?basis_t*mu;
%%
%?plot?訓(xùn)練結(jié)果
ytrain?=?basis（:idx）*mu;
figure;plot（ytrain）;hold?on;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????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