Новое распределение пропуска
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 таким способом который