Пользовательские произвольные постоянные имена для sympy dsolve?

Question

Есть ли способ обеспечить получение dsolve() для использования произвольных констант, названных другой последовательностью, чем C1, C2 и т. Д.?

cse() разрешает параметр symbols, который принимает бесконечный генератор имен, но я не вижу аналогичного параметра dsolve().

Я решаю схемы, полные символических емкостей, и константы dsolve запутывают.

Если это не так, есть ли быстрый способ заменить произвольные константы другими, такими как k_1, k_2 и т. Д.? Похоже, что dsolve() использует C1, C2 и т. Д. Для своих констант, в то время как мой код использует C_1, C_2 и т. Д. Для емкостей. Я мог бы изменить свое название повсюду, чтобы использовать нестандартный символ емкости, но я бы предпочел не делать этого.

---

С благодарностью @Marshmallow, я начал использовать эту процедуру, которая оборачивает dsolve позволяет мне изменить символы после факта. Там по-прежнему риск столкновения, но это достаточно хорошо сейчас:

def newdsolve(eq,*args,**kwds): 
    ss=kwds.get('symbols') 
    sln=dsolve(eq,*args,**kwds) 
    psyms=((newdsolve.csre.match(n.name),n) for n in eqVc.free_symbols) 
    if ss and isinstance(ss,str): 
     subsdict=dict([(n[1],'k_'+n[0].group(1)) for n in psyms if n[0]]) 
    elif ss: 
     subsdict=dict([(n[1],next(ss)) for n in psyms if n[0]]) 
    else: subsdict={} 
    return sln.subs(subsdict) 
newdsolve.csre=re.compile(r'C(\d+)') 

Answer 1

Я не вижу такой вариант в dsolve либо, но вы можете использовать subs заменять символы C1, C2 с вашим собственным. Функция sub_const заменяет их, например, k_1, k_2 и т. Д., Пока выражение не перестанет меняться. Другой символ можно использовать вместо k: это параметр.

from sympy import * 

def sub_const(expr, char): 
    i = 1 
    while True: 
     old = symbols('C{}'.format(i)) 
     new = symbols(char + '_{}'.format(i)) 
     new_expr = expr.subs(old, new) 
     if new_expr == expr: 
      return expr 
     expr = new_expr 
     i += 1 

x, f = symbols('x, f') 
print(sub_const(dsolve(Derivative(f(x), x, x) + 9*f(x), f(x)), 'k')) 

Выход:

Eq(f(x), k_1*sin(3*x) + k_2*cos(3*x))