Решение почти блок-диагональной линейной системы
x = slvblk(blokmat,b)
x = slvblk(blockmat,b,w)
x = slvblk(blokmat,b) возвращает решение (при наличии) линейной системы Ax = b, с матрицей A хранится в blokmat в сплайновой почти блок-диагональной форме. В настоящее время только команда spcol обеспечивает такое описание матрицы, типичной записью которой является значение некоторой производной (включая 0-ю производную, то есть значение) B-сплайна на некотором участке. Если линейная система переопределена (т.е. имеет больше уравнений, чем неизвестные, но имеет полный ранг), то возвращается решение наименьших квадратов.
Правая сторона b может содержать несколько столбцов и, как ожидается, будет содержать столько строк, сколько есть строк в матрице, описанной blokmat.
x = slvblk(blockmat,b,w) возвращает вектор x что минимизирует взвешенную сумму ((Ax-b) (j)) 2.
sp=spmak(knots,slvblk(spcol(knots,k,x,1),y.')) обеспечивает в sp B-форма сплайна s порядка k с узловой последовательностью knots который соответствует заданным данным (x,y), т.е. для которых(x) равняется y.
Команда bkbrk используется для получения существенных частей матрицы коэффициентов, описанной blokmat (в одной из двух доступных форм).
QR-факторизация выполняется для каждого диагонального блока, после того как она была дополнена уравнениями, не рассматриваемыми при факторизации предыдущего блока. Полученная факторизация затем используется для решения линейной системы путем обратной подстановки.