论文上说要用三次样条插值，但是相应的matlab程序用cubic出warning，说是换成了spline，程序里要迭代的，最后收敛。然后翻译成mma的代码，Interpolation没有cubic，只有spline和hermite。然后程序不能收敛，看了一下第一步的插值结果，相差0.001这么个量级。
interpolation怎么用三次样条插值呢，说明里没有cubic这个词。我需要插出来二阶可导。Interpolation的spline好像是二阶可导的啊，怎么会不收敛！！
程序比较大，所以就说说了。如果研究了一下觉得是我程序的问题，再来贴。

另外这个程序的matlab上用的interpn，是个不推荐使用要淘汰的指令。
另外做的是二元函数的插值。

interpolatingorder

无力了，不一样的。传来文件，已经影响我对mma的信心了。代码太长，希望有人能解决，很渺茫，放在这里了。里面SHOCK函数和这个没关系的，我已经注释掉了。
[无效] http://pan.baidu.com/s/1jGrb1gq

三次样条是要给边界的吧，matlab和mma自动给的边界是不是不同

我说啊，虽然我不知道你的程序最后是怎么改的，但是，你真的确定matlab和mathematica对高维矩阵各维度的排列是一样的？：
mathematica.stackexchange.com/a/10584/1871

我不太清楚
mathematica.stackexchange.com/a/10584/1871
上的关于mma转换到matlab的结果如何用到matlab上，我觉得这可能是适用于matlink的？
举个例子，我使用的是就是对于mma的单行或者单列要特殊处理一下，其他的没问题，包括赋值。

查到了B样条的一个性质，k次B样条的k-1次导数连续。B样条确实和平常的样条不一样，和matlab里平常的样条，不一样。mma只用B样条，窃以为这方法是足够的。