使用SymPy求解递推关系时rsolve返回None的问题咨询

• 系数精度问题：代码中a = 1/6是Python原生浮点数（循环小数的二进制近似值），不是SymPy可识别的精确有理数，rsolve对混有浮点系数的常系数线性递推鲁棒性不足，无法完成方程结构识别，直接返回None。
• 初值索引不匹配：rsolve求解六阶常系数线性递推时，默认求解域为n ≥ 6，要求提供的初值为y(0), y(1), y(2), y(3), y(4), y(5)共6个初始项。代码中传入的是y(1)到y(6)的初值，存在索引偏移，SymPy无法将给定初值对应到通解的系数求解环节，也会触发返回None的逻辑。
• 额外说明：该递推的特征方程存在有理根r=1，理论上属于SymPy可求解的常系数递推范畴，不存在方程本身不可解的问题。

• 替换精确系数：将浮点数写法的a = 1/6改为SymPy精确有理数形式a = sp.Rational(1,6)，保证递推方程的所有系数都是SymPy可识别的精确值，不混入Python原生浮点类型。
• 修正初值索引：顺着原递推逻辑补全缺失的y(0)项。根据递推式y(6) = a*(y(5)+y(4)+y(3)+y(2)+y(1)+y(0))，对比代码中y6的计算式y6 = a*(y5+y4+y3+y2+y1) +a，可推导出y0 = 1。随后将初值集合调整为y(0)到y(5)的6个项，匹配rsolve的索引规则。
• 备用兜底方案：如果调整参数后rsolve仍返回None，可手动拆分特征方程结构：原递推的特征方程为r**6 - a*(r**5 + r**4 + r**3 + r**2 + r +1) =0，可因式分解为(r-1)*(6*r**5 +5*r**4 +4*r**3 +3*r**2 +2*r +1)=0，提取所有特征根后写出通解形式，代入已知初值解线性方程组求待定系数，即可得到递推的闭式解。

import sympy as sp
# 替换为SymPy精确有理数
a = sp.Rational(1, 6)

y = sp.Function('y')
n = sp.symbols('n', integer=True)
f = y(n) - a * (y(n-1) + y(n-2) + y(n-3) + y(n-4) + y(n-5) + y(n-6))

# 按递推逻辑计算匹配索引的初值
y0 = 1
y1 = a
y2 = y1 * a + a
y3 = y2 * a + y1 * a + a
y4 = y3 * a + y2 * a + y1 * a + a
y5 = y4 * a + y3 * a + y2 * a + y1 * a + a

# 传入和rsolve索引规则匹配的初值集合
s = sp.rsolve(f, y(n), {
    y(0): y0,
    y(1): y1,
    y(2): y2,
    y(3): y3,
    y(4): y4,
    y(5): y5
})

print(sp.simplify(s))

内容的提问来源于stack exchange，提问作者peter.petrov