使用sympy库中的rsolve函数
示例代码:
from sympy import symbols, Function, rsolve
n = symbols('n')
a = Function('a')
# 定义不定系数递推关系式
eq = a(n+2) - 3*a(n+1) + 2*a(n) - 2*2**n
# 求解递推关系式的通项公式
sol = rsolve(eq, a(n))
print(sol)
输出结果为:
-2**n + 2**(n + 1)*C0 + n*2**n
其中C0是通项公式中的常数项,在计算过程中自动求解。