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;plot（Ytrain‘r‘）
%?plot?測(cè)試結(jié)果
figure;
plot（y）;hold?on;plot（Ytest‘r‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????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