Octave / task_7

x = [-6:2:6];
y = (x .^ 3) / 10 + 2 * x .^ 2;
x1 = [-6:0.1:6];
y_linear = line_interp(x, y, x1);
y_square = square_interp(x, y, x1);

m_linear = interp1(x, y, x1);
m_cub = pchip(x, y, x1);
m_spline = spline(x, y, x1);

figure;
hold on
y1 = (x1 .^ 3) / 10 + 2 * x1 .^ 2;
plot(x1, y_linear, 'b'); #лин интерполяция
plot(x1, y_square, 'm'); #квадратичная интерполяция
y_point = [50.4, 25.6, 7.2, 0, 8.8, 38.4, 93.6]; #точки из файла
plot(x, y_point, 'g'); #график по точкам из файла
plot(x1, y1, 'r'); #график функции
xlabel('X')
ylabel('Y')
grid on
legend('лин. интерполяция',
'квадр. интерполяция',
'по точкам из файла',
'истинная функция');

figure;
hold on
plot(x1, m_linear, 'm-'); #лин интерполяция octave
plot(x1, m_cub, 'b--'); #кубическая интерполяция octave
plot(x1, m_spline, 'g--'); #интреполяция сплайнами octave
plot(x, y_point, 'ro'); #график по точкам из файла
plot(x1, y1, 'r-'); #график функции
xlabel('X')
ylabel('Y')
grid on
legend('лин. интерполяция octave',
'куб. интерполяция octave',
'интерполяция сплайнами octave',
'по точкам из файла',
'истинная функция');