本文整理汇总了Python中sympy.Poly.subs方法的典型用法代码示例。如果您正苦于以下问题：Python Poly.subs方法的具体用法？Python Poly.subs怎么用？Python Poly.subs使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类sympy.Poly的用法示例。
在下文中一共展示了Poly.subs方法的3个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: ratint_ratpart
# 需要导入模块: from sympy import Poly [as 别名]
# 或者: from sympy.Poly import subs [as 别名]
def ratint_ratpart(f, g, x):
"""Horowitz-Ostrogradsky algorithm.
Given a field K and polynomials f and g in K[x], such that f and g
are coprime and deg(f) < deg(g), returns fractions A and B in K(x),
such that f/g = A' + B and B has square-free denominator.
"""
f, g = Poly(f, x), Poly(g, x)
u = poly_gcd(g, g.diff())
v = poly_div(g, u)[0]
n = u.degree - 1
m = v.degree - 1
d = g.degree
A_coeff = [ Symbol('a' + str(n-i), dummy=True) for i in xrange(0, n+1) ]
B_coeff = [ Symbol('b' + str(m-i), dummy=True) for i in xrange(0, m+1) ]
symbols = A_coeff + B_coeff
A = Poly(zip(A_coeff, xrange(n, -1, -1)), x)
B = Poly(zip(B_coeff, xrange(m, -1, -1)), x)
H = f - A.diff()*v + A*poly_div(u.diff()*v, u)[0] - B*u
result = solve(H.coeffs, symbols)
A = A.subs(result)
B = B.subs(result)
rat_part = Poly.cancel((A, u), x)
log_part = Poly.cancel((B, v), x)
return rat_part, log_part示例2: test_subs
# 需要导入模块: from sympy import Poly [as 别名]
# 或者: from sympy.Poly import subs [as 别名]
def test_subs():
p = Poly(t*x*y**2 + x*y + t**2, x, y)
assert p.subs(x, 2) == \
Poly(((2*t, 2, t**2), ((2,), (1,), (0,))), y)
assert p.subs(y, 2) == \
Poly(((2 + 4*t, t**2), ((1,), (0,))), x)
assert p.subs(x, y) == \
Poly(((t, 1, t**2), ((3,), (2,), (0,))), y)
assert p.subs(y, x) == \
Poly(((t, 1, t**2), ((3,), (2,), (0,))), x)
assert p.subs(x, z) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), z, y)
assert p.subs(y, z) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), x, z)
assert p.subs(t, z) == \
Poly(((z, 1, z**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(z, t) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(t, x) == \
Poly(((1, 1, 1), ((2, 2), (2, 0), (1, 1))), x, y)
assert p.subs(t, y) == \
Poly(((1, 1, 1), ((1, 3), (1, 1), (0, 2))), x, y)
assert p.subs(t, sin(1)) == \
Poly(((sin(1), 1, sin(1)**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(t, sin(x)) == \
x*y**2*sin(x) + x*y + sin(x)**2
assert p.subs(x, sin(x)) == \
t*y**2*sin(x) + y*sin(x) + t**2示例3: test_subs
# 需要导入模块: from sympy import Poly [as 别名]
# 或者: from sympy.Poly import subs [as 别名]
def test_subs():
p = Poly(t*x*y**2 + x*y + t**2, x, y)
assert p.subs(x, 2) == \
Poly(((2*t, 2, t**2), ((2,), (1,), (0,))), y)
assert p.subs(y, 2) == \
Poly(((2 + 4*t, t**2), ((1,), (0,))), x)
assert p.subs(x, y) == \
Poly(((t, 1, t**2), ((3,), (2,), (0,))), y)
assert p.subs(y, x) == \
Poly(((t, 1, t**2), ((3,), (2,), (0,))), x)
assert p.subs(x, z) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), z, y)
assert p.subs(y, z) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), x, z)
assert p.subs(t, z) == \
Poly(((z, 1, z**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(z, t) == \
Poly(((t, 1, t**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(t, x) == \
Poly(((1, 1, 1), ((2, 2), (2, 0), (1, 1))), x, y)
assert p.subs(t, y) == \
Poly(((1, 1, 1), ((1, 3), (1, 1), (0, 2))), x, y)
assert p.subs(t, sin(1)) == \
Poly(((sin(1), 1, sin(1)**2), ((1, 2), (1, 1), (0, 0))), x, y)
assert p.subs(t, sin(x)) == \
x*y**2*sin(x) + x*y + sin(x)**2
assert p.subs(x, sin(x)) == \
t*y**2*sin(x) + y*sin(x) + t**2
coeffs = [ Symbol('a' + str(i)) for i in [3,2,1,0] ]
values = [ 4, 0, 0, 1 ]
f = Poly(zip(coeffs, [3,2,1,0]), x)
assert f.subs(dict(zip(coeffs, values))) == Poly(4*x**3+1, x)