Можно использовать следующие команды с любым сплайном в качестве примера, такие как cs
, ch
и примеры sp
, созданные в разделе Cubic Spline Interpolation.
Сначала создайте сплайн, например:
sp = spmak(1:6,0:2)
fnplt(sp)
Чтобы получить значение в a
, используйте синтаксис fnval(f,a)
, например:
fnval(sp,4)
Создать вторую производную сплайна:
DDf = fnder(fnder(sp))
Альтернативный способ создать вторую производную:
DDf = fnder(sp,2);
Получить определенный интеграл сплайна на интервале [a
.. b
], в этом примере от 2 до 5:
diff(fnval(fnint(sp),[2;5]))
Чтобы вычислить различие между двумя сплайнами, используйте форму fncmb(sp1,'-',sp2)
, например:
fncmb(sp,'-',DDf);