% skrypt33
disp('skrypt interpoluje 4 typy wielomianow do funkcji 2sin(x)+sin(2x)+sin(3x+pi)');
disp('dla punktow od 0 do 20
x=0:20;
y=2*sin(x)+sin(2*x)+sin(3*x+pi);
pkt=0:0.1:20;
pierwszy=interp1(x,y,pkt,'linear');
trzeci=interp1(x,y,pkt,'pchip');
sasiedzi=interp1(x,y,pkt,'nearest');
sklejane=interp1(x,y,pkt,'spline');
cla;
subplot(2,2,1);
hold on;
plot(pkt,pierwszy,'--b',x,y,'xr');
fplot('2*sin(x)+sin(2*x)+sin(3*x+pi)',[0 20],'-k');
title('pierwszy');
hold off;
subplot(2,2,2);
hold on;
plot(pkt,trzeci,'--g',x,y,'xr');
fplot('2*sin(x)+sin(2*x)+sin(3*x+pi)',[0 20],'-k');
title('trzeci');
hold off;
subplot(2,2,3);
hold on;
plot(pkt,sasiedzi,'--c',x,y,'xr');
fplot('2*sin(x)+sin(2*x)+sin(3*x+pi)',[0 20],'-k');
title('sasiedzi');
hold off;
subplot(2,2,4);
hold on;
plot(pkt,sklejane,'--m',x,y,'xr');
fplot('2*sin(x)+sin(2*x)+sin(3*x+pi)',[0 20],'-k');
title('sklejane');
hold off;
