如何解决使用 SymPy 解决递归问题
我正在尝试使用 SymPy 解决这个重复问题,但它给了我 None
有什么想法吗?
import sympy as sp
a = 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) )
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
y6 = y5 * a + y4 * a + y3 * a + y2 * a + y1 * a + a
s = sp.rsolve(f,y(n),{
y(1) : y1,y(2) : y2,y(3) : y3,y(4) : y4,y(5) : y5,y(6) : y6
}
)
print(s)
解决方法
[编辑回复]
尝试重写您的重复项,使其引用有限数量的先前术语,而不是所有先前术语:
>>> from sympy import Rational
>>> a = Rational(1,6)
>>> rsolve(y(n) - (1 + a)*y(n-1),y(n),[a])
(7/6)**n/6
>>> [_.subs(n,i) for i in range(6)] == [y1,y2,y3,y4,y5,y6]
True
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。