clc, clear all, close all
%% WARUNKI POCZATKOWE
k= 0:9; % ilosc danych z tabeli
y= [-1.8 -1.2 0.5 1.6 1.0 0.1 -0.6 -0.8 -0.2 0.8];
s= [-1.8 -1.2 0.5 1.6 1.0 0.1 -0.6 -0.8 -0.2 0.8]% dane z tabeli
T= 0.005; % okres probkowania
l=1;
m=0;
P=1:(max(k)/T)+1; % prelokacja zmiennej P
u=0:1:9
%% I - INTERPOLACJA 3-PUNKTOWA
for i= 1:length(y)-2
d= m:T:2;
for h=1:length(d)
P(l)= 0.5*(2*y(i)-d(h)*(3*y(i)-4*y(i+1)+y(i+2))+(d(h)^2)*...
(y(i)-2*y(i+1)+y(i+2)));
l=l+1;
end
m=1+T;
end
%% GRAFICZNE PRZEDSTAWIENIE WYNIKU INTERPOLACJI 3-KROKOWEJ
t=0:T:9;
figure(1)
plot(t,P); hold on; grid on; plot(k,y,'ro');
title('Interpolacja 3-punktowa'); xlabel('Próbki'); ylabel('Amplituda');
legend('Funkcja interpolowana','Probki');
%% II - INTERPOLACJA Z UZYCIEM FUNKCJI MATLABA
y1 = interp1(k,y,t, 'linear'); % Interpolacja liniowa
y2 = interp1(k,y,t, 'cubic'); % Interpolacja wielomianowa
y3 = interp1(k,y,t, 'spline'); % Interpolacja pryey sklejanie funkcji
%% GRAFICZNE PRZEDSTAWIENIE WYNIKOW INTERPOLACJI Z FUNKCJI MATLABA
figure(2)
plot(t,y1, 'k'); hold on; grid on;
plot(t,y2, 'g'); plot(t,y3, 'm'); plot(k,y, 'ok');
title('Interpolacja z użyciem wbudowanych funkcji Matlaba');
xlabel('Próbki'); ylabel('Amplituda');
legend('Metoda Liniowa','Metoda Wielomianowa','Sklejanie Funckji','Probki')
figure(3)
plot(k,y,'ok')
hold on
plot(t,P);
plot(t,y2, 'g')
grid on
pierwszy\
t=0:1:9
p=[-1.8 -1.2 0.5 1.6 1.0 0.1 -0.6 -0.8 -0.2 0.8]
x=0:(1/3):9
y=interp1(t,p,x,'cubic');
figure(1)
plot(t,p,'oblack')
hold on
xlabel('czas [probki]')
ylabel('amlituda')
plot(x,y,'r');
y=interp1(t,p,x,'spline');
plot(x,y,'b');
y=interp1(t,p,x,'linear');
plot(x,y,'g');
legend('probki','cubic','spline','linear')
grid on