newknt

Новое распределение пропуска

Синтаксис

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 таким способом который

newknots(i)newknots(i+1)|Dkf|1/k=const, i=k:k+newl1