met3Robaka.doc

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