Переходы, т.е. f (x +) -f (x-)
jumps = fnjmp(f,x)
jumps = fnjmp(f,x)
как fnval(f,x)
кроме того, что он возвращает f перехода (x +) - f (x -) через x
(а не значение в x
) функции, f описывается как f
и что он работает только для одномерных функций.
Это функция для специалистов по сплайнам.
fnjmp(ppmak(1:4,1:3),1:4)
возвращает вектор [0,1,1,0]
начиная с pp
функция здесь 1 on [1.. 2], 2 на [2.. 3], и 3 на [3.. 4], следовательно, имеет нулевой переход на 1 и 4 и переход на 1 на обоих 2 и 3.
Если x
является cos([4:-1:0]*pi/4)
, затем fnjmp(fnder(spmak(x,1),3),x)
возвращает вектор [12 -24 24 -24 12]
(до округления). Это согласуется с тем фактом, что рассматриваемый сплайн является так называемым совершенным кубическим B-сплайном, т.е. имеет абсолютно постоянную третью производную (на своём базовом интервале). Измененная команда
fnjmp(fnder(fn2fm(spmak(x,1),'pp'),3),x)
возвращает вместо этого вектор [0 -24 24 -24 0]
, в соответствии с тем фактом, что, в отличие от B-формы, сплайн в ппформе не имеет разрыва ни в одной из его производных в конечных точках его основного интервала. Обратите внимание, что fnjmp(fnder(spmak(x,1),3),-x)
возвращает вектор [12,0,0,0,12]
начиная с -x
, хотя теоретически равен x
, отличается от x
путем округления, следовательно, третья производная B-сплайна, обеспечиваемая spmak(x,1)
не имеет перехода между -x(2)
, -x(3)
, и -x(4)
.