使用Chebyshev节点和Lagrange多项式绘图报list与int减运算不支持错误求助

• 报错unsupported operand type(s) for -: 'list' and 'int'的核心原因：你传入的g是Python原生列表类型，原生列表不支持直接和数值做减法运算，只有numpy数组支持广播式的批量数值运算。另外你传入的节点参数d也是原生列表，运算时也存在同类隐患。
• 隐藏逻辑错误：lagrange函数中p += f[i]*L的缩进位置错误，放在了j循环内部，会导致每计算一个j项就重复累加L，不符合拉格朗日基函数的计算逻辑，就算类型报错修好也会得到错误结果。
• 你之前遇到的x, y, and format string must not be None报错，大概率是调用函数时没有正确接收返回值，绘图时传入了空值导致的。

import numpy as np
import matplotlib.pyplot as plt

def lagrange(x_nodes, z_target, y_nodes):
    d = len(x_nodes)
    # 输入转numpy数组，避免原生列表运算报错
    x_nodes = np.asarray(x_nodes)
    z_target = np.asarray(z_target)
    p = np.zeros_like(z_target)
    for i in range(d):
        L = np.ones_like(z_target)
        for j in range(d):
            if j != i:
                L *= (z_target - x_nodes[j]) / (x_nodes[i] - x_nodes[j])
        # 基函数计算完成后再加到结果上，缩进移到j循环外
        p += y_nodes[i] * L
    return p

def f(x):
    return np.cos(x)

# 插值节点x坐标
x_nodes = [0.1, 0.2, 4]
# 生成密集的目标计算点，方便绘图
z_target = np.linspace(-0.5, 4.5, 100)
# 插值节点对应的y值
y_nodes = f(np.asarray(x_nodes))

# 计算插值结果
interp_res = lagrange(x_nodes, z_target, y_nodes)

# 绘图逻辑
plt.plot(z_target, f(z_target), label='原始cos(x)')
plt.plot(z_target, interp_res, label='拉格朗日插值结果', linestyle='--')
plt.scatter(x_nodes, y_nodes, color='red', label='插值节点')
plt.legend()
plt.show()

• 将所有参与数值运算的列表都通过np.asarray()转为numpy数组，支持批量数值运算，解决类型不匹配报错
• 修正了累加插值结果的缩进位置，移到j循环外部，符合拉格朗日插值的计算逻辑
• 补充了完整的绘图逻辑，避免绘图时传入空值的旧报错
• 优化参数命名，明确区分插值节点和目标计算点，避免逻辑混淆

内容的提问来源于stack exchange，提问作者Ryan McAree