scatteredInterpolant
插值2-D或3-D分散数据
描述
使用scatteredInterpolant对二维或三维数据集进行插值分散的数据．scatteredInterpolant返回interpolant.F.对于给定的数据集。你可以评估F.在一组查询点，例如(yq xq)在2-D中，产生内插值vq = F (xq yq)．
使用griddedInterpolant用…进行内插网格数据．
创建
语法
描述
F.=散射interpolant.
F.= scatteredInterpolant (xgydF4y2Ba那ygydF4y2Ba那V.）xgydF4y2Ba和ygydF4y2Ba指定(x, y)样本点的坐标。V.是包含与点关联的示例值的矢量(x, y)．
F.= scatteredInterpolant (xgydF4y2Ba那ygydF4y2Ba那Z.那V.）
F.= scatteredInterpolant (___那方法那ExtrapolationMethod）方法和ExtrapolationMethod作为前两个语法中的任何一个中的最后两个输入参数。
方法可以是：“最近的”那“线性”,或“天然”．ExtrapolationMethod可以是：“最近的”那“线性”,或'没有'．
输入参数
属性
使用
使用scatteredInterpolant创建interpolant.那F.．然后你可以评估F.在特定的点使用下列任何语法:
VQ = F（PQ）
VQ = F（XQ，YQ）
Vq = F (Xq, Yq Zq)
Vq = F ({xq, yq})
Vq = F ({xq, yq zq})
VQ = F（PQ）指定矩阵中的查询点魁人党．在每一行魁人党包含查询点的坐标。VQ = F（XQ，YQ）和Vq = F (Xq, Yq Zq)将查询点指定为两个或三个大小相同的矩阵。Vq = F ({xq, yq})和Vq = F ({xq, yq zq})指定查询点为网格向量．使用此语法在要查询的大网格时要保存内存。
例子
二维插值
定义一些样本点，并计算在这些位置的三角函数值。这些点是插值的样本值。
t = linspace(3/4 *π,2 *π,50)';x = [3 * cos (t);2 * cos (t);0.7 * cost);y = [3 * sin (t);2 * sin (t);0.7 * sin (t)];v = repelem ([-0.5;1.5;2)、长度(t));
创建interpolant。
F = scatteredInterpolant (x, y, v);
在查询位置求插值(xq那yq）。
tq = linspace (3/4 * pi + 0.2, 2 *π- 0.2,40)';xq = (2.8 * cos (tq操作);1.7 * cos (tq操作);因为(tq)];yq = (2.8 * sin (tq操作);1.7 * sin (tq操作);sin (tq)];vq = F (xq, yq);
策划的结果。
Plot3（x，y，v，'。'，xq，yq，vq，'。')、网格在标题(线性插值的)包含('x'），Ylabel（'是'), zlabel (“价值”)传说(样本数据的那“插入查询数据”那“位置”那“最好”）
三维插值
为一组分散的采样点创建插值，然后在一组3-D查询点处评估内插。
定义200个随机点并采样一个三角函数。这些点是插值的样本值。
rng默认的；P = -2.5 + 5*rand([200 3]);v =罪(P(: 1)。^ 2 +(:, 2页)。^ 2 + P(: 3)。^ 2)。/ (P(: 1)。^ 2 + P (:, 2) ^ 2 + P(:, 3)。^ 2);
创建interpolant。
F = scatteredInterpolant (P, v);
在查询位置求插值(xq那yq那的zq）。
[xq, yq zq] = meshgrid (2:0.25:2);vq = F (xq, yq zq);
绘制结果的切片。
xslice =(闲置,1 2);yslice = (0, 2);zslice = (2,0);片(xq, yq zq、vq xslice, yslice, zslice)
更换样本值
中的元素替换价值属性，当您想要更改样本点上的值时。当计算新的插值函数时，可以立即得到结果，因为原来的三角剖分没有改变。
创建50个随机点并采样指数函数。这些点是插值的样本值。
rng (“默认”）x = -2.5 + 5 * rand（[50 1]）;Y = -2.5 + 5 * rand（[50 1]）;v = x。* exp（-x。^ 2-y。^ 2）;
创建interpolant。
F = scatteredInterpolant (x, y, v)
f =具有属性的散射interpolant：point：[50x2双]值：[50x1 double]方法：'linear'外推方法：'linear'
求插值（1.40,1.90）．
F (1.40, 1.90)
ans = 0.0069
改变插值的样本值，并在同一点重新计算插值。
x ^2 + y ^2;F.Values = vnew;F (1.40, 1.90)
ans = 5.6491
消除重复样本点
使用groupsummary消除重复的样本点，并在调用前控制它们的组合方式scatteredInterpolant．
创建一个200 × 3的样本点位置矩阵。在最后五行添加重复的点。
P = -2.5 + 5*rand(200,3);P (197:200:) = repmat (P(196年:),4,1);
在采样点上创建一个由随机值组成的向量。
V =兰德(大小(P, 1), 1);
如果你试图使用scatteredInterpolant对于重复的样本点，它抛出一个警告并对相应的值取平均值V.产生一个唯一的点。但是，你可以使用groupsummary在创建插值之前消除重复点。如果您想使用除平均之外的方法合并重复点，这是特别有用的。
使用groupsummary消除重复的样本点并保持最大值V.在重复的样本点位置。指定样本点矩阵作为分组变量，对应的值作为数据。
[V_unique, P_unique] = groupsummary (V, P, @max);
由于分组变量有三列，groupsummary返回唯一的组P_unique作为单元阵列。将单元阵列转换回矩阵。
p_unique = [p_unique {：}];
创建interpolant。由于样本点现在是唯一的，scatteredInterpolant不会抛出警告。
我= scatteredInterpolant (P_unique V_unique);
比较分散的数据插值方法
比较了几种不同的插值算法的结果scatteredInterpolant．
创建一个由50个分散点组成的样本数据集。点的数量被人为地缩小，以突出插值方法之间的差异。
x = -3 + 6 *兰特（50,1）;Y = -3 + 6 *兰特（50,1）;v = sin（x）。^ 4。* cos（y）;
创建插值和查询点网格。
F = scatteredInterpolant (x, y, v);[xq, yq] = meshgrid (3:0.1:3);
绘制结果使用“最近的”那“线性”,“天然”方法。每次插值方法改变时，您都需要查询插值器以获得更新的结果。
F.Method =.“最近的”；VQ1 = F（XQ，YQ）;Plot3（x，y，v，'莫')举行在网格（XQ，YQ，VQ1）标题（'最近的邻居')传说('样本点'那插值结果的那“位置”那“西北”）
F.Method =.“线性”；VQ2 = F（XQ，YQ）;图Plot3（x，y，v，'莫')举行在网格（XQ，YQ，VQ2）标题（“线性”)传说('样本点'那插值结果的那“位置”那“西北”）
F.Method =.“天然”；vq3 = F (xq, yq);图Plot3（x，y，v，'莫')举行在网格(xq yq vq3)标题(“天然的邻居”)传说('样本点'那插值结果的那“位置”那“西北”）
画出精确的解。
图Plot3（x，y，v，'莫')举行在网格(xq yq,罪(xq)。^4 .* cos(yq)) title('确切的解决方案')传说('样本点'那准确的表面的那“位置”那“西北”）
二维外推
使用最近的邻外推出在凸船外部的单点处查询插值。
定义一个由200个随机点组成的矩阵，并对一个指数函数进行抽样。这些点是插值的样本值。
rng (“默认”) P = -2.5 + 5*rand([200 2]);x = P (: 1);y = P (:, 2);v = x。* exp（-x。^ 2-y。^ 2）;
创建插值，指定线性插值和最近的邻外推。
f =散射interpolant（p，v，“线性”那“最近的”）
F = scatteredInterpolant with properties: Points: [200x2 double] Values: [200x1 double] Method: 'linear' ExtrapolationMethod: 'nearest'
求凸包外的插值。
VQ = F（3.0，-1.5）
VQ = 0.0029.
禁用外推并评估F.在同一点。
F.ExtrapolationMethod ='没有'；VQ = F（3.0，-1.5）
vq =南
更多关于
提示
求a的值比较快
scatteredInterpolant对象F.在许多不同的查询点集合，而不是使用函数单独计算插值griddata.或griddatan．例如：%快速创建插值F和评估多次f =散射interpolant（x，y，v）v1 = f（xq1，yq1）v2 = f（xq2，yq2）%使用网格数据单独计算插值较慢v1 = griddata（x，y，v，xq1，yq1）v2 = griddata（x，y，v，xq2，yq2）
要改变插值的样本值或插值方法，更新插值对象的属性更有效
F.比它创造一个新的scatteredInterpolant对象。当你更新价值或方法，输入数据的底层Delaunay三角剖分不会改变，因此可以快速计算新的结果。分散数据插值
scatteredInterpolant使用数据的Delaunay三角剖分，因此对样本点的缩放问题很敏感xgydF4y2Ba那ygydF4y2Ba那Z.,或P.．当出现这种情况时，您可以使用正常化重新归类数据并改进结果。看用不同的幅度对数据进行规范化为更多的信息。
算法
scatteredInterpolant使用分散的样本点的Delaunay三角测量来执行插值［1］．
参考文献
[1] Amidror,艾萨克。电子成像系统的散点数据插值方法:综述电子成像杂志．第11卷，第2期，2002年4月，157-176页。
你也可以从以下列表中选择一个网站:
