Кубический сплайн сглаживания
Примечание
Для более простого, но менее гибкого метода, чтобы сгенерировать сплайны сглаживания, попробуйте приложение Curve Fitting или fit функция.
pp = csaps(x,y)(x,y) в ppform. Значение сплайна f на сайте данных x(j) аппроксимирует значение данных y(:,j) для j = 1:length(x).
Сплайн сглаживания f минимизирует
Здесь, n является количеством записей x и интеграл на самом маленьком интервале, содержащем все записи x. y j и x j отсылает к j th записи y и x, соответственно. D2f обозначает вторую производную функционального f.
Значения по умолчанию для ошибочных весов меры w j равняются 1. Значением по умолчанию для кусочной постоянной функции веса λ в мере по шероховатости является постоянная функция 1. По умолчанию, csaps выбирает значение для параметра сглаживания p на основе сайтов определенных данных x.
Чтобы оценить сглаживание шлицуют вне его основного интервала, необходимо сначала экстраполировать его. Используйте команду pp = fnxtr(pp) гарантировать, что вторая производная является нулем вне интервала, заполненного по условию сайты.
[___] = csaps({x1,...,xm}, обеспечивает ppform y,___)m- сглаживание продукта тензора варьируемой величины шлицует к данным по прямоугольной сетке, описанной {x1,...,xm}. Можно использовать этот синтаксис с любым из аргументов в предыдущих синтаксисах.
Подходящие сплайны сглаживания с помощью csaps функция с различными значениями для параметра сглаживания p. Используйте значения p между экстремальными значениями 0 и 1, чтобы видеть, как они влияют на форму и близость подходящего сплайна.
Загрузите набор данных титана.
[x, y] = titanium();
Когда p = 0, s0 подгонка прямой линии наименьших квадратов к данным. Когда p = 1S1 вариационный, или естественный, кубический сплайн interpolant.
Для 0 < p < 1, sp сплайн сглаживания, который является компромиссом между этими двумя экстремальными значениями: более сглаженный, чем interpolant s1 и ближе к данным, чем прямая линия s0.
p = 0.00009; s0 = csaps(x,y,0); sp = csaps(x,y,p); s1 = csaps(x,y,1);
figure fnplt(s0); hold on fnplt(sp); fnplt(s1); plot(x,y,'ko'); hold off title('Smoothing splines with different values for p'); legend('p = 0', ['p = ' num2str( p )], 'p = 1', 'Location', 'northwest')
Настройте параметр сглаживания, ошибочные веса меры и веса меры по шероховатости.
Создайте синусоиду с шумом.
x = linspace(0,2*pi,21); y = sin(x)+(rand(1,21)-.5)*.3;
Соответствуйте сплайну сглаживания к данным. Задайте параметр сглаживания p = 0.4 и ошибка измеряет веса w это варьируется через данные.
pp = csaps(x,y,0.4,[],[ones(1,10),repmat(5,1,10), 0]);
Функция возвращает сглаженную подгонку к зашумленным данным, которые намного ближе к данным в правильной половине из-за намного большего ошибочного веса меры там. Обратите внимание на то, что ошибка при взвешивании нуля для последней точки данных исключает эту точку из подгонки.
Теперь соответствуйте сплайну сглаживания с помощью тех же данных, сглаживая параметр и ошибочные веса меры, но с настроенными весами меры по шероховатости.
pp1 = csaps(x,y, [.4,ones(1,10),repmat(.2,1,10)], [], ...
[ones(1,10), repmat(5,1,10), 0]);Вес меры по шероховатости - только 0,2 в правильной половине интервала. Соответственно, подгонка более груба, но ближе на правой стороне данных (за исключением последней точки данных, которая проигнорирована).
Постройте оба соответствуют сравнению.
figure hold on fnplt(pp, 'b'); fnplt(pp1,'r--') plot(x,y,'ok') hold off ylim([-1.5 1.5]) title(['Cubic smoothing spline, with right half treated ',... 'differently']) legend('Larger error weight', 'Larger error and smaller roughness weight')
Соответствуйте сплайну сглаживания к двумерным данным, сгенерированным peaks функция с добавленным универсальным шумом. Используйте csaps получить новые, сглаживавшие точки данных и параметры сглаживания csaps определяет для подгонки.
Создайте сетку. В данном примере сетка является 51 61 регулярной координатной сеткой.
x = {linspace(-2,3,51),linspace(-3,3,61)};
[xx,yy] = ndgrid(x{1},x{2}); Сгенерируйте зашумленные данные с помощью peaks функциональные и случайные числа в интервале .
y = peaks(xx, yy);
noisy = y + (rand(size(y)) - 0.5);
figure
surf(xx,yy,noisy)
axis off
Соответствуйте данным. Используйте csaps получить сглаживавшие значения данных оценило по сетке x и параметр сглаживания значения по умолчанию используется в подгонке.
[sval,p] = csaps(x,noisy,[],x);
График подгонки показывает, что некоторая шероховатость остается. Обратите внимание на то, что необходимо транспонировать массив sval.
figure
surf(x{1},x{2},sval.')
axis off
Для несколько более сглаженного приближения задайте значение для p это немного меньше, чем csaps значение по умолчанию.
ssval = csaps(x,noisy,.996,x);
figure
surf(x{1},x{2},ssval.')
axis off
csaps реализация стандартной программы Фортрана SMOOTH из PGS.
Вычисление сплайна сглаживания требует решения линейной системы, матрица коэффициентов которой имеет форму p*A + (1-p)*B, с матрицами A и B в зависимости от сайтов данных x. Значение по умолчанию p делает p*trace(A) равный (1-p)*trace(B).