Skip to main content
 首页 » 编程设计

python之Python 中的 MATLAB interp2 函数

2025年05月04日7开发

我需要一个与 interp2 MATLAB 函数等效的 Python。我正在尝试制作这个 MATLAB example在 Python 中工作,但我不会。

import numpy as np 
 
from scipy.interpolate import interp2d 
from scipy.interpolate import RectBivariateSpline 
 
service = np.array(range(10, 31, 10)) 
years = np.array(range(1950, 1991, 10)) 
 
wage = np.array([[150.6970,199.5920,187.6250], 
    [179.3230, 195.0720, 250.2870], 
    [203.2120, 179.0920, 322.7670], 
    [226.5050, 153.7060, 426.7300], 
    [249.6330, 120.2810, 598.2430]]) 
 
ip = RectBivariateSpline(years, service, wage) 
 
print(ip(15, 1975)) 

但是我得到了这个我无法解决的错误(在 RectBivariateSpline 中):

    Traceback (most recent call last): 
  File "/Users/andrea/Documents/workspace/PythonProjects/pyArmBot/src/foo.py", line 15, in <module> 
    ip = RectBivariateSpline(years, service, wage) 
  File "/Library/Python/2.7/site-packages/scipy-0.10.1-py2.7-macosx-10.7-intel.egg/scipy/interpolate/fitpack2.py", line 728, in __init__ 
    kx,ky,s) 

请您参考如下方法:

我相信您需要使用以下行而不是您的打印行来调用 ip 以返回一个值。

print(ip.ev(15, 1975))

更新:

非常简单的调整。并盯着我的脸。您可以设置 kx 和 ky 的值,它们是双变量样条的度数 See documentation here

无论如何,只需将线调整为:

ip = RectBivariateSpline(years, service, wage, kx=2, ky=2)

如果你不想,你不需要调整kx。如果您将 kx 保留为 3,我怀疑它会大大改变插值。