clear all
clf

X=linspace(0,2*pi,90);
R=rand(size(X));
idx=find(R>0.6);
x=X(idx)';

t=sin(x)+0.3*randn(size(x));
echo on
%%%%%%%> oto nasze dane
plot(X,sin(X),x,t,'o');
legend('oryginalna funkcja', 'zaszumione dane pochodzace z tej funkcji',0);
pause

%%%%%%%> teraz zbudujemy siec RBF

net = rbf(1, 5, 1, 'gaussian');
options(1, :) = foptions;
options(2, :) = foptions;
pause

%%%%%%%> i ja przecwiczymy na naszych danych
net = rbftrain(net, options, x, t);
pause

%%%%%%%> teraz przepuscimy dane przez siec i narysujemy wynik
Y = rbffwd(net, X');
plot(X,sin(X),x,t,'o',X,Y);
legend('oryginalna funkcja', 'zaszumione dane pochodzace z tej funkcji','interpolacja siecia RBF',0);
pause
%%%%%%%> teraz zbadamy jak nasza siec ekstrapoluje
X=linspace(0,3*pi,100);
Y = rbffwd(net, X');
plot(X,sin(X),x,t,'o',X,Y);
legend('oryginalna funkcja', 'zaszumione dane pochodzace z tej funkcji','extrapolacja siecia RBF',0);

echo off