Новое распределение пропуска
newknots = newknt(f,newl)
newknt(f)
[...,distfn] = newknt(...)
newknots = newknt(f,newl) возвращает последовательность узла, внутренние узлы которой сокращают основной интервал f в части newl таким способом, которым определенная кусочная линейная монотонная функция, связанная с высокой производной f, является equidistributed.
Намерение состоит в том, чтобы выбрать последовательность узла, подходящую для прекрасного приближения функционального g, грубое приближение которого в f принято, чтобы содержать достаточно информации о g, чтобы сделать это выполнимым.
newknt(f) использование для newl его значение по умолчанию, а именно, количество полиномиальных частей в f.
[...,distfn] = newknt(...) также возвращается, в distfn, ppform той кусочной линейной монотонной функции, являющейся equidistributed.
Если ошибка в приближении наименьших квадратов sp к некоторым данным x,y сплайном порядка, k кажется неровным, вы можете попробовать за больше справедливого распределения узлов при помощи
spap2(newknt(sp),k,x,y);
Для другого примера смотрите Решение Нелинейного ОДУ с Пограничным слоем Словосочетанием.
Это - стандартная программа Фортрана NEWNOT в PGS. С k порядок функции кусочного полинома f в pp функция |Dkf | аппроксимирована кусочной постоянной функцией, полученной локальным, дискретным, дифференцированием изменения D k –1f. Новая последовательность пропуска выбрана, чтобы подразделить основной интервал кусочного полинома f таким способом который