sympy

Question

में काल्पनिक जड़ों पर ध्यान न दें मैं एक बहुपद हल करने के लिए sympy उपयोग कर रहा हूँ:

x = Symbol('x') 
y = solve(int(row["scaleA"])*x**3 + int(row["scaleB"])*x**2 + int(row["scaleC"])*x + int(row["scaleD"]), x) 

y संभव समाधान की एक सूची है। हालांकि, मुझे काल्पनिक लोगों को अनदेखा करने और केवल वास्तविक समाधानों का उपयोग करने की आवश्यकता है। साथ ही, मैं समाधान को एक मान के रूप में एक अभिव्यक्ति के रूप में नहीं चाहता हूं। अभी यह ऐसा लगता है:

[-2/3 - 55**(1/3)*(-1/2 - sqrt(3)*I/2)/3, -2/3 - 55**(1/3)*(-1/2 + sqrt(3)*I/2)/3, -55**(1/3)/3 - 2/3] 

मुझे अंतिम अभिव्यक्ति के मूल्य (-2.22756) की आवश्यकता है। क्या इसे सरल बनाने के लिए सिम्पी में कार्य हैं?

Answer 1

रूप Krastonov कहा था mpmath एक आसान विधि प्रदान की:

y = polyroots([int(row["scaleA"]), int(row["scaleB"]), int(row["scaleC"]), int(row["scaleD"])-value]) 
for root in y: 
    if "j" not in str(root): 
     value = root 

Answer 2

मैं अभिव्यक्ति का मूल्यांकन करने के लिए चरित्र "I" युक्त समाधानों को अनदेखा करने में कामयाब रहा और .evalf() का उपयोग किया। कोड है:

x = Symbol('x') 
    y = solve(int(row["scaleA"])*x**3 + int(row["scaleB"])*x**2 + int(row["scaleC"])*x + int(row["scaleD"]), x) 
    for root in y: 
     if "I" not in str(root): 
      print("This One:" + str(root.evalf())) 

Answer 3

आप x सेट करते हैं असली होने के लिए, SymPy केवल आप वास्तविक समाधान दे देंगे

x = Symbol('x', real=True) 
solve(..., x) 

Answer 4

solve() में विभिन्न प्रकार के समाधानों के लिए लगातार उत्पादन नहीं है, कृपया solveset(Eq,x,domain=S.Reals):

का उपयोग करें

from sympy import ImageSet, S 
x = Symbol('x') 
y = solveset(int(row["scaleA"])*x**3 + int(row["scaleB"])*x**2+int(row["scaleC"])*x + int(row["scaleD"]), x, domain=S.Reals) 

http://docs.sympy.org/latest/modules/solvers/solveset.html

Answer 5

यह ठीक बात की तरह है कि real_roots के लिए बना दिया है और अपने मामले जहां गुणांक पूर्णांक हैं करने के लिए विशेष रूप से लागू होता है है:

x = Symbol('x') 
eq = int(row["scaleA"])*x**3 + int(row["scaleB"])*x**2 + int(row["scaleC"])*x + int(row["scaleD"]) 
y = real_roots(eq, x) # gives [CRootOf(...), ...] 

CRootOf उदाहरणों में से मान हो सकता है आपको जो भी सटीकता चाहिए, उसका मूल्यांकन किया गया है और इसमें कोई काल्पनिक हिस्सा नहीं होना चाहिए। उदाहरण के लिए,

>>> [i.n(12) for i in real_roots(3*x**3 - 2*x**2 + 7*x - 9, x)] 
[1.07951904858] 

नोट: मुझे याद है के रूप में, वापस जड़ों कि यह मान्यताओं से मुलाकात की पुष्टि करने के लिए (में सक्षम नहीं था भेज देंगे हल यानी अगर वे इस धारणा है तो वे वापस आ रहे हैं के लिए गलत नहीं पाए गए)। इसके अलावा, अगर आप हल से अधिक सुसंगत आउटपुट चाहते हैं, तो @PyRick, ध्वज dict=True सेट करें।