如何在Python中使用Z3查找Lagrange插值多项式？

Question

我正在使用Z3库和Python语言学习数学，但是我不知道如何找到基于拉格朗日插值多项式的多项式？ 例如，

I have
f(0) = y0
f(1) = y1
f(2) = y2
....

我想找到f（X）

我与Scipy配合良好，但是当更改为Z3库时，找不到更多文档了：(

Answer 1

我不知道如何在Z3的背景下回答您的问题。 Z3不包含用于拉格朗日插值的任何内置求解器。 您可以设置一个多项式等式的系统，其中的系数保持符号化。 然后，您可以询问Z3该系统是否可解决。 它会为系数生成一般的代数数，因此我不确定这是否适合您的应用程序。 当然，在python中，您可以编写：

a, b, c = Real('a b c')
def f(X):
    return a*X*X + b*X + c

solve(f(0) == y0, f(1) == y1, f(2) == y2)

其中y0 ， y1 ， y2是常数。