interp2

插值2-d中meshgrid格式栅格数据

描述

VQ= interp2(X,YVXQ,Y染色体长臂返回内插使用线性内插在特定的查询点的两个变量的函数的值。结果总是通过功能的原始采样。Xÿ包含样品点的坐标。V包含在每个采样点对应的函数值。XQYQ包含查询点的坐标。

VQ= interp2(VXQ,Y染色体长臂假设采样点的默认网格。默认的网格点覆盖的矩形区域,X = 1:n的Y = 1:米,其中[M,N] =尺寸(V)。当你想节省内存,并且不关心两点之间的绝对距离此语法。

VQ= interp2(V返回上通过将采样值之间的间隔一次在每个维度形成的精制网格中的内插值。

VQ= interp2(Vķ返回上通过重复地减半的间隔形成一精制网格内插值ķ次在每个维度。这导致2 ^ k-1个样本值之间内插点。

VQ= interp2(___方法指定的替代的内插方法:“线性”“最近”'立方体''makima', 要么“花”。默认的方法是“线性”

VQ= interp2(___方法extrapval还规定extrapval,分配给所有查询标量值横亘该采样点的域的外部。

如果省略extrapval对于样本点的域之外的查询参数,则基于该方法争论interp2返回下列之一:

  • 对于外推值“花”'makima'方法

  • 为NaN对于其它的内插方法的值

例子

全部收缩

粗样品功能。

[X,Y] = meshgrid(-3:3);V =峰(X,Y);

画出粗采样。

图冲浪(X,Y,V)标题(“原始采样”);

创建0.25间距查询网格。

[Xq的,Y染色体长臂] = meshgrid(-3:0.25:3);

内插在查询点。

VQ = interp2(X,Y,V,Xq的,Y染色体长臂);

绘制的结果。

图冲浪(Xq的,YQ VQ);标题(“线性插值使用更细的栅格”);

粗略采样峰功能。

[X,Y] = meshgrid(-3:3);V =峰(7);

画出粗采样。

图冲浪(X,Y,V)标题(“原始采样”);

创建0.25间距查询网格。

[Xq的,Y染色体长臂] = meshgrid(-3:0.25:3);

内插在查询点,并指定立方插值。

VQ = interp2(X,Y,V,Xq的,YQ'立方体');

绘制的结果。

图冲浪(Xq的,YQ VQ);标题(“立方插值在更精细的网格”);

加载一些图像数据到工作空间。

加载flujet.mat色彩表灰色

隔离图像的小区域,并将其转换为单精度。

V =单(X(200:300,1:25));

显示图像的区域。

于imagesc(V);轴标题(“原始图像”

插入通过反复把在每个维度上的精制网格五倍的点之间的间隔内插的值。

VQ = interp2(V,5);

显示结果。

于imagesc(VQ);轴标题('线性插值'

粗采样函数过的范围内,[-2,2]在这两个维度。

[X,Y] = meshgrid(-2:0.75:2);R = SQRT(X. ^ 2 + Y ^ 2)+ EPS;V = SIN(R)./(R);

画出粗采样。

图冲浪(X,Y,V)XLIM([ -  4 4])ylim([ -  4 4])标题(“原始采样”

创建超出的域名查询网格Xÿ

[Xq的,Y染色体长臂] = meshgrid(-3:0.2:3);

域内执行三次插值Xÿ,并将其赋值会超出零所有查询。

VQ = interp2(X,Y,V,Xq的,YQ'立方体',0);

绘制的结果。

图冲浪(Xq的,YQ VQ)标题(“三次插值与X和Y的VQ = 0外域”);

输入参数

全部收缩

样品网格点,指定为实数矩阵或向量。样品网格点必须是唯一的。

  • 如果Xÿ是矩阵,然后它们包含的坐标满格(在meshgrid格式)。使用meshgrid函数来创建Xÿ矩阵在一起。这两个矩阵必须是相同的大小。

  • 如果Xÿ被矢量,则它们被视为网格矢量。在两种载体中的值必须严格单调,无论是增加或减少。

注意

在未来的版本中,interp2将不接受用于样品和查询网格的行和列向量的混合组合。相反,必须通过构建满格meshgrid。或者,如果你有一个大的数据集,你可以使用griddedInterpolant代替interp2

例:[X,Y] = meshgrid(1:30,-10:10)

数据类型:|

样本值,指定为实数或复数矩阵。对于尺寸要求V取决于大小Xÿ

  • 如果Xÿ是代表一个完整的网格矩阵(在meshgrid格式),然后V必须大小相同Xÿ

  • 如果Xÿ是格矢量,然后V必须是含有基质冗长)长度(X)列。

如果V包含复杂的数字,然后interp2插值单独的实部和虚部。

例:兰特(10,10)

数据类型:|
复数支持:万博1manbetx

查询点,指定为实标量,矢量,矩阵或阵列。

  • 如果XQYQ是标量,那么他们是一个单一的查询点的坐标。

  • 如果XQYQ是不同取向的矢量,然后XQYQ被视为网格向量。

  • 如果XQYQ具有相同的大小和方向的矢量,然后XQYQ被视为散点在2-d的空间。

  • 如果XQYQ是矩阵,然后它们代表任一的查询点的满格(在meshgrid格式)或散点。

  • 如果XQYQ是N d阵列中,2-d的空间然后它们代表散射点。

注意

在未来的版本中,interp2将不接受用于样品和查询网格的行和列向量的混合组合。相反,必须通过构建满格meshgrid。或者,如果你有一个大的数据集,你可以使用griddedInterpolant代替interp2

例:[Xq的,Y染色体长臂] = meshgrid((1:0.1:10),( - 5:0.1:0))

数据类型:|

细化系数,指定作为一个真正的,非负整数标量。此值指定的次数反复划分细化网格的间隔在每个维度的数目。这导致2 ^ k-1个样本值之间内插点。

如果ķ0, 然后VQ是相同的V

interp2(V,1)是相同的interp2(V)

内插值的下图显示了放置(红色)之间的九个取样值(黑色),用于k = 2时

例:interp2(V,2)

数据类型:|

插值方法,规定为在该表中的选项之一。

方法 描述 连续性 注释
“线性” 在查询点内插的值在每个相应维度相邻网格点是基于所述值的线性内插。这是默认的插值方法。 C0
  • 需要在每一维的至少两个网格点

  • 需要的不仅仅是更多的内存“最近”

“最近” 在查询点的插补值是在最近的采样网格点的值。 间断
  • 要求在每个维度上的两个网格点。

  • 适度的内存需求最快的计算

'立方体' 在查询点内插的值在每个相应维度相邻网格点是基于三次插值的值的。内插基于三次卷积。 C1
  • 电网必须有统一的间距在每个维度,但间距不必对所有尺寸相同

  • 要求在每个维度上至少四个点

  • 需要更多的内存和计算时间比“线性”

'makima' 改性阿克玛三次Hermite内插。在查询点的插补值最多三个使用每个相应维度相邻网格点的值进行评价是基于与多项式的分段函数。所述阿克玛式进行修饰以避免过冲。 C1
  • 需要在每个维度上至少2分

  • 产生比起伏少“花”

  • 计算时间一般小于“花”,但对内存的要求是相似的

“花” 在查询点内插的值在每个相应维度相邻网格点是基于三次插值的值的。内插是使用不-A-结结束条件基于三次样条。 C2
  • 需要在每个维度四点

  • 需要更多的内存和计算时间比'立方体'

函数值的外部域Xÿ,指定为实数或复数标量。interp2返回外域的所有点这个常数的值Xÿ

例:

例:5 + 1I

数据类型:|
复数支持:万博1manbetx

输出参数

全部收缩

内插的值,返回为实数或复数标量,矢量或矩阵。的大小和形状VQ取决于你使用的语法,在某些情况下,输入参数的大小和值。

语法 特殊情况 VQ的大小
interp2(X,Y,V,Xq的,Y染色体长臂)
interp2(V,Xq的,Y染色体长臂)
这些语法包括的和变化方法要么extrapval
XQYQ是标量 纯量 大小(VQ)= [1 1]当你通过XQYQ作为标量。
同上 XQYQ具有相同的大小和方向的矢量 相同尺寸和方向的矢量XQYQ 如果大小(Xq的)= [100 1]
大小(Y染色体长臂)= [100 1]
然后大小(VQ)= [100 1]
同上 XQYQ是混合取向的矢量 矩阵,其中行的数量是长度(Y染色体长臂)和列的数目是长度(Xq的) 如果大小(Xq的)= [1 100]
大小(Y染色体长臂)= [50 1]
然后大小(VQ)= [50 100]
同上 XQYQ具有相同的尺寸的矩阵或阵列 矩阵或相同尺寸的阵列XQYQ 如果大小(Xq的)= [50 25]
大小(Y染色体长臂)= [50 25]
然后大小(VQ)= [50 25]
interp2(V,k)的
该语法的和变化包括方法要么extrapval
没有

矩阵,其中行数是:
2 ^ K *(大小(V,1)-1)+1

和列数为:
2 ^ K *(大小(V,2)-1)+1

如果大小(V)= [10 20]
k = 2时
然后大小(VQ)= [37 77]

更多关于

全部收缩

严格单调

一组值是始终增加或减少,而不会逆转的。例如,序列,A = [2 4 6 8]是严格单调和提高。序列,B = [2 4 4 6 8]是不是严格单调的,因为在之间的值没有变化B(2)B(3)。序列,C = [2 4 6 8 6]包含之间的逆转C(4)C(5),所以它是不是在所有单调。

全网格(在meshgrid格式)

对于interp2中,完整的网格是一对矩阵,其元素代表点的超过一个矩形区域的网格。一个矩阵包含X坐标 - ,另一个矩阵包含ÿ坐标 - 。在值X- 矩阵是严格单调并沿行增加。沿其列中的值是恒定的。在值ÿ- 矩阵是严格单调和沿着列增加。沿其行的值是恒定的。使用meshgrid函数来创建一个完整的网格,可以传递给interp2

例如,下面的代码创建该区域的满格,-1≤X≤3和1≤ÿ≤4:

[X,Y] = meshgrid(-1:3,(1:4))
X = -1 0 1 2 3 -1 0 1 2 3 -1 0 1 2 3 -1 0 1 2 3 Y = 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4

格矢量来表示的网格比全网格更紧凑的格式。两种格式和样本值的矩阵之间的关系V

网格矢量

对于interp2,网格矢量由一对矢量的定义的X- 和ÿ坐标 - 在网格中。行向量定义X坐标 - ,和列向量定义ÿ坐标 - 。

例如,下面的代码创建网格向量指定的区域中,-1≤X≤3和1≤ÿ≤4:

X = -1:3;Y =(1:4)';

散点

对于interp2散点由一对限定散在2-d空间中的点的集合阵列的。一个阵列包含X坐标 - 和另一个包含ÿ坐标 - 。

例如,下面的代码指定的点,(2,7),(5,3),(4,1),和(10,9):

X = [2 5;4 10];Y = [7 3;1 9];

扩展功能

R2006a前推出