如何 有效地 利用 contourplot中的数据点 得到 二维插值函数？

Question

下面是我的代码，但插值效果和计算效率都不好。

f[x_, y_] := Sin[x y]
(* 实际上，f是一个复杂得多的，含有数值积分的函数 *)

datTemp = 
  ContourPlot[f[x, y], {x, 0, 3}, {y, 0, 3}][[1, 1]] // 
   DeleteDuplicatesBy[First];

dat = {#[[1]], #[[2]], f[#[[1]], #[[2]]]} & /@ datTemp;

finterp = Interpolation[dat]

Plot3D[{f[x, y], finterp[x, y]}, {x, 0, 3}, {y, 0, 3}]

@namespace url(http://www.w3.org/1999/xhtml); /*背景图片*/

Answer 1

你对二维函数插值把x相同的数据点只保留一个，相当于y方向的变化都忽略了，那肯定插不准啊。。。

把DeleteDuplicatesBy[First]去掉就可以了，或者你要担心有重复点的话就改用DeleteDuplicates