(Не рекомендуется) Интерполяция разрозненных данных
TriScatteredInterp не рекомендуется. Использовать scatteredInterpolant вместо этого.
TriScatteredInterp используется для интерполяции рассеянного набора данных, находящегося в 2-D или 3-D пространстве. Разрозненный набор данных, определенный местоположениями X и соответствующие значения V может быть интерполирован с помощью триангуляции Делоне X. Это приводит к образованию поверхности формы V = F(X). Поверхность можно вычислить в любом месте запроса. QX, использование QV = F(QX), где QX лежит в пределах выпуклого корпуса X. interpolant F всегда проходит через точки данных, указанные образцом.
F = TriScatteredInterp создает пустой интерполятор рассеянных данных.
F = TriScatteredInterp( создает интерполятор, который подходит к поверхности формы Q,V)V = F(Q) к разрозненным данным в (Q, V). Q является матрицей размера mptsоколо-ndim, где mpts - количество точек и ndim - размер пространства, в котором находятся точки (ndim равно 2 или 3). Вектор столбца V определяет значения в Q, где длина V равняется mpts.
F = TriScatteredInterp( использует указанное DT,V)DelaunayTri объект DT в качестве основы для вычисления интерполятора. DT является триангуляцией Делоне рассредоточенных местоположений данных, DT.X. Матрица DT.X имеет размер mptsоколо-ndim, где mpts - количество точек и ndim - размер пространства, в котором находятся точки, 2 <= ndim <= 3. V - вектор столбца, определяющий значения в DT.X, где длина V равняется mpts.
F = TriScatteredInterp(___, определяет метод интерполяции method)method используется для интерполяции данных. Можно использовать любую из предыдущих комбинаций входных аргументов.
Чтобы оценить интерполятор, выразите утверждение в форме Монжа Vq = F(Xq), Vq = F(Xq,Yq), или Vq = F(Xq,Yq,Zq) где Vq - значение интерполятора в местоположении запроса и Xq, Yq, и Zq - векторы местоположений точек.
Создайте набор данных из 2-D случайных рассеянных точек.
rng default
x = rand(100,1)*4-2;
y = rand(100,1)*4-2;
z = x.*exp(-x.^2-y.^2);Создайте интерполятор для данных.
F = TriScatteredInterp(x,y,z);
Оценка интерполятора в местоположениях точек запроса (qx,qy). Соответствующие интерполированные значения в этих местоположениях возвращаются в qz.
ti = -2:.15:2; [qx,qy] = meshgrid(ti,ti); qz = F(qx,qy);
Постройте график данных и интерполированной поверхности.
mesh(qx,qy,qz) hold on plot3(x,y,z,'o')
delaunayTriangulation | interp1 | interp2 | interp3 | meshgrid | scatteredInterpolant | triangulation