Cтепенная интерполяция X по матрице M с вычислением производной
y = interpolwithderi(x,M,dYdx);
x - входная скалярная величина;
M - матрица коэффициентов для кубической сплайн-интерполяции;
dYdx - выходное значение, текущая производная текущего выходного значения функции.
interpolwithderi – функция получения интерполированного значения y и его производной dYdx по входному значению x и матрице интерполяционных коэффициентов M. Матрицу коэффициентов можно получить например с помощью функции cspline. В случае, если нет необходимости получать производную интерполированной функции, можно использовать функции языка программирования interpol, linear_interpn.
Функция возвращает переменную типа double, полученное путем сплайн-интерполяции по заданным аргументам.
input u; output y,dydu; initialization Xs = linspace(0,100,1000); Ys = sin(Xs); M = cspline(Xs,Ys); end; y = interpolwithderi(u,M,dydu);
Текст программы для блока типа 0985.dita#. Происходит получение коэффициентов интерполяционной матрицы, после чего по полученной матрице и входному значению формируются значение функции и ее производной.