Новое распределение пропуска
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 | аппроксимирован кусочной постоянной функцией, полученной локальным, дискретным, дифференцированием изменения Dk–1f. Новая последовательность пропуска выбрана, чтобы подразделить основной интервал кусочного полинома f таким способом который