clear;
clc;
%p0=[-10:0.1:10];  %輸入數(shù)據(jù)
p0=20*rand(1,100)-10;
o=0.2;   %噪聲正態(tài)分布的標(biāo)準(zhǔn)差
lC=size(p0,2);
z=o*randn(1,lC);   %輸出時伴隨的噪聲
for i=1:lC;
    t0(i)=sin(p0(i))/p0(i)+z(i);    %輸出數(shù)據(jù)
end
    
%LC=75;  %學(xué)習(xí)樣本數(shù) 
%P=p0(1:LC);      
%T=t0(1:LC);     
X=p0';            %訓(xùn)練用輸入數(shù)據(jù)
Y=t0';            %訓(xùn)練用輸出數(shù)據(jù)
ker='bspline';
C=1;
loss='einsensitive';
e=0.3;
global p1;
p1=1;

[nsv,beta,bias]=svr(X,Y,ker,C,loss,e)
x=[-10:0.1:10];
y=sin(x)./x;
%tstX=x';
%targetY=y';
%tstY=svroutput(X,tstX,ker,beta,bias);
%err=svrerror(X,tstX,targetY,ker,beta,bias,loss,e)  %計(jì)算預(yù)測誤差

figure;
plot(x,y) 
hold on;
svrplot(X,Y,ker,beta,bias,e)
hold off;