Контрольная работа: Аппроксимация функций 2

Решая эту систему, получим следующие зависимости для вычисления коэффициентов:

_____________________________________________________________

function [k]=Spl(aa,n,x,y);

c(1)=0;

b(1)=10*y(1,2)-10*y(1,1)-0.1*c(1);

for k=1:n-2

b(k+1)=0.2*c(k)+b(k);

c(k+1)=100*y(1,k+2)-100*y(1,k+1)-10*b(k+1);

end

j=floor(10*aa-109);

if j==6

j=5;

end

k=y(1,j)+b(j)*(aa-x(j))+c(j)*(aa-x(j))^2;

_____________________________________________________________

n=6;

clear yy; clear ff; clear x1; clear x1

for i=1:11

a=10.95+i*0.05;

ff(i)=Spline(a,n,x,y);

x3(i)=10.95+0.05*i;

end

for j=1:6

yy(j)=y(1,j);

x1(j)=x(j);

end

% Построение графика