% 
% РспуЭггйуз кбй РбсемвплЮ
%
% ГЯнефбй чсЮуз фщн ухнбсфЮуещн: polyfit, spline, interp1
%

x = [ -1, 0, 1];		% иЭуейт рпх гнщсЯжпхме фйт фймЭт
y = [0, 2, 6];			% фймЭт уе ухгкексймЭнет иЭуейт
xi = -2:.1:2;			% дйЬуфзмб уфп прпЯп иб хрплпгЯупхме рспуеггЯуейт

N = 2;											% кбипсйумьт вбимпэ фпх рплхщнэмпх
p = polyfit( x, y, N);						% рспуЭггйуз ме рплхюнхмп Н вбимпэ
yi1 = polyval( p, xi);						% екфЯмзуз рплхщнэмпх
method = 'nearest';							% мЭипдпй рбсемвплЮт: nearest, linear, spline, pchip, cubic, v5cubic
yi2 = interp1( x, y, xi, method);		% рбсемвплЮ
yi3 = spline( x, y, xi);					% ецбсмпгЮ фщн кхвйкюн ухнбсфЮуещн splines

% ГсбцйкЭт рбсбуфЬуейт
h = plot( x, y, '^r', xi, yi1, 'k', xi, yi2, 'b', xi, yi3, 'g');
axis tight;
grid;
l = legend( 'ДедпмЭнб', 'Рплхюнхмп', 'РбсемвплЮ', 'Splines');

% МефбвплЮ мпсцЮт уфз гсбцйкЮ рбсЬуфбуз
set( l, 'fontname', 'Arial Greek');
set( h(1), 'markersize', 12);
set( h(2), 'linewidth', 4);
set( h(3), 'linewidth', 2);
set( h(4), 'linewidth', 2);
