Contents

Johdattelua interpolaatioon

L2demo1

clc
close all
xdata=[-2 -1 3]
ydata=[1 -2 5];
plot(xdata,ydata,'*')
xdata =
    -2    -1     3
 A=[ones(3,1) [-2;-1;3] [-2;-1;3].^2]
 c=A\ydata'
A =
     1    -2     4
     1    -1     1
     1     3     9
c =
   -3.1000
   -0.1500
    0.9500
axis([-2.5 3.5 -5 10 ])
grid on
hold on
x=linspace(-2.1,3.1,100);
p=c(1)+c(2)*x+c(3)*x.^2;
length(p)
ans =
   100
plot(x,p)
plot(x,p,'r')   % Eri varilla.

Vektorien lineaarikombinaatio "riviajattelulla"

plot(x,p,'r')
V=[ones(100,1) x' (x.^2)'];
y=V*c;
plot(x,y,'g')
title('Interpolaatiopoly 3:n pisteen dataan')