Подходящие значения на рассеянных 2D сайтах со сплайнами сглаживания Тонкой Пластины

Сплайны продукта Tensor хороши для с координатной сеткой (двумерный и даже многомерный) данные. Для работы с рассеянными двумерными данными тулбокс обеспечивает сплайн сглаживания тонкой пластины. Предположим, что у вас есть значения определенных данных y(j) на сайтах данных, имеющий разброс x(:,j), j=1:N, в плоскости. Дать определенный пример,

n = 65; t = linspace(0,2*pi,n+1); 
x = [cos(t);sin(t)]; x(:,end) = [0;0];

обеспечивает 65 сайтов, а именно, 64 точки, равномерно распределенные на модульном круге, плюс центр того круга. Вот соответствующие значения данных, а именно, шумные значения очень хорошей функции.

y = (x(1,:)+.5).^2 + (x(2,:)+.5).^2;
noisy = y + (rand(size(y))-.5)/3;

Затем можно вычислить разумное приближение к этим данным

st = tpaps(x,noisy);

и постройте получившееся приближение наряду с зашумленными данными

fnplt(st); hold on
plot3(x(1,:),x(2,:),noisy,'wo','markerfacecolor','k')
hold off

и тем самым произведите следующее изображение:

Тонкая пластина, сглаживающая приближение сплайна к зашумленным данным