Функция вычисления матрицы коэффициентов для линейной интерполяции по заданным координатам.
M = linterp(X, Y);
Х – входной массив координат по оси x для вычисления матрицы коэффициентов,
Y – входной массив координат по оси y для вычисления матрицы коэффициентов.
Функция вычисляет матрицу коэффициентов для линейной интерполяции по заданным координатам точек. Вычисленная матрица коэффициентов используется функцией interpol для линейной интерполяции аргумента x. Входные массивы X и Y должны быть одного размера. Входные массивы X, Y могут задаваться:
M = linterp (X, Y);
M = linterp ([x1,x2,x3,x4],[y1,y2,y3,y4]);
M = linterp ([x1,x2,x3,x4],Y);
M = linterp (X,[y1,y2,y3,y4]);
M = linterp ([-1.80,-1.60,-1.40,-1.20],Y);
M = linterp (X,[-1.80,-1.60,-1.40,-1.20]);
M = linterp ([-1.80,-1.60,-1.40,-1.20],[-1.40,-0.78,-0.53,-0.35]);
M – вычисленная матрица коэффициентов.
input Q;//вход - аргумент output n;//выход - результат кубической сплайн-интерполяции //массивы координат для вычисления матрицы const Qpts = [-2.00, -1.80, -1.60, -1.40, -1.20, -1.00, -0.80, -0.60, -0.40, -0.20, 0.00, 0.20, 0.40, 0.60, 0.80, 1.00, 1.20, 1.40, 1.60, 1.80, 2.00]; const Hp = [-1.40,-0.78,-0.53,-0.35,-0.17, 0.05, 0.38, 0.80, 0.94, 1.07, 1.15, 1.26, 1.33, 1.40, 1.47, 1.40, 1.45, 1.51, 1.59, 1.69, 1.74]; Mn = linterp(Qpts, Hp); n = interpol(Q, Mn);
В результате переменной n будет присвоено интерполированное значение напора (по массивам Hp и Qpts) в соответствии со значением расхода Q.