如何将scipy.interpolate.interp2d生成的2D插值函数转为SymPy符号形式

报错原因 #

• 你的代码存在两个核心逻辑错误：
  1. scipy.interpolate.interp2d 是纯数值实现的接口，仅支持传入数值类型的输入，传入SymPy符号对象会触发类型转换失败的报错。
  2. sympy.lambdify 的功能是将SymPy符号表达式转换为可执行的数值函数，你的用法刚好颠倒了调用逻辑，自然无法得到预期结果。

最优实现方案（适配Q4单元场景） #

你明确说明插值由4个角点生成，对应Q4有限元单元的双线性形函数，本身存在显式的符号表达式，完全不需要通过Scipy插值对象中转，可以直接构造符号插值函数：

import sympy as sp

# 定义自然坐标符号变量
r, s = sp.symbols('r s')

# 替换为你四个角点对应的实际Z值，也可以保留符号形式参与后续计算
z_bottom_left = 1.0
z_bottom_right = 2.0
z_top_right = 3.0
z_top_left = 4.0

# Q4单元标准双线性形函数（自然坐标范围[-1, 1]）
N1 = (1 - r) * (1 - s) / 4  # 左下节点形函数
N2 = (1 + r) * (1 - s) / 4  # 右下节点形函数
N3 = (1 + r) * (1 + s) / 4  # 右上节点形函数
N4 = (1 - r) * (1 + s) / 4  # 左上节点形函数

# 构造完整符号插值函数
f_symbolic = N1*z_bottom_left + N2*z_bottom_right + N3*z_top_right + N4*z_top_left

# 直接执行符号微分操作示例
# 对r求一阶偏导
df_dr = sp.diff(f_symbolic, r)
# 对s求一阶偏导
df_ds = sp.diff(f_symbolic, s)
# 代入数值计算示例
f_value = f_symbolic.subs({r:0, s:0}).evalf()

通用场景方案（从interp2d提取参数转符号） #

如果你需要适配其他插值类型的通用转换逻辑，可以提取interp2d对象内部存储的节点、系数参数，手动拼接符号表达式：
线性插值的interp2d对象的tck属性存储了x轴节点、y轴节点、插值系数三个数组，提取后按照对应插值算法的规则即可拼接出对应的SymPy表达式。

内容的提问来源于stack exchange，提问作者rada-dev