这份名单插值是一个简单的函数,一个2D点映射到一个号码,如果 把每个{{x,y},z}
为f[x,y]=z
力数学对非结构化张格
{
{{1,3},9}, {{1,4},16},
{{2,4},8}, {{2,5},10}
}
我现在想要一个函数,插值/外推任何{x,y}
。
数学拒绝这样做:
Interpolation[{{{1,3},9}, {{1,4},16},{{2,4},8}, {{2,5},10}},
InterpolationOrder->1]
插值:: indim:本 坐标不处在一个结构性 张量积网格。
我明白为什么(数学想要一个“矩形”域),但 什么是强迫数学创建插值最简单的方法?
这不起作用:
f[1,3]=9; f[1,4]=16; f[2,4]=8; f[2,5]=10;
g=FunctionInterpolation[f[x,y],{x,1,2},{y,3,5}]
FunctionInterpolation :: nreal:
16附近{X,Y} = {1, - },这个函数不 计算为一个真正的数。 5 FunctionInterpolation :: nreal:
17在{x,y} = {1, - }附近,函数没有将 评估为实数。 5 FunctionInterpolation :: nreal:
18在{x,y} = {1, - }附近,函数没有将 评估为实数。 5 General :: stop:进一步输出 FunctionInterpolation :: nreal 将在此计算过程中被抑制。
即使你忽略上面的警告,评估摹给出错误
g[1.5,4] // FortranForm
f(1.5,4) + 0.*(-9.999999999999991*(f(1.4,4) - f(1.5,4)) +
- 0.10000000000000009*
- (9.999999999999991*
- (9.999999999999991*(f(1.4,4) - f(1.5,4)) +
- 4.999999999999996*(-f(1.4,4) + f(1.6,4))) +
- 0.5000000000000006*
- (-10.000000000000014*
- (-3.333333333333333*(f(1.3,4) - f(1.6,4)) -
- 4.999999999999996*(-f(1.4,4) + f(1.6,4))) -
- 9.999999999999991*
- (9.999999999999991*(f(1.4,4) - f(1.5,4)) +
- 4.999999999999996*(-f(1.4,4) + f(1.6,4))))))
其他“明显”的想法(内插值功能 本身)也不起作用。
靖岩崎,欢迎StackOverflow的数学界,并感谢您的贡献。 – 2011-05-07 06:35:16