gevraagd is de concentratie in 1994 berekend door interpolatie met de lagrangeveelterm
Dit is de code van de functie
Code: Selecteer alles
function lagrange
X=[1978 1980 1982 1984 1986 1988 1990 1992];
f=[12; 12.7; 13; 15.2; 18.2; 19.8; 24.1; 28.1];
P=0;
x=1978:0.2:1994;
for j=1:numel(X)
L = 1;
for i=1:numel(X)
if (i~=j)
L=L.*(x-X(i))/(X(j)-X(i));
end
end
P=P+L*f(j);
end
plot(x,P)
end
Voor 1994 krijg ik een concentratie van -38.4
dit kan natuurlijk allemaal niet in werkelijkheid...
Anderzijds denk ik intuïtief wel dat eens je een punt zoekt buiten het interval van je spilpunten, je een grotere fout krijgt
Maar dit vind ik toch wel beetje extreem...
Dus kan zon grote (in werkelijkheid foute) afwijking wel?