exponenta event banner

Постобработка сплайнов

С любым примером сплайна можно использовать следующие команды, например: cs, ch и sp примеры, построенные в разделе Интерполяция кубических сплайнов.

Сначала создайте сплайн, например:

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);