行列式，矩阵，不能转换表达式

Question

我正在写一个简短的程序，它应该找到函数实部和虚部都为零的值。我不明白为什么在运行程序后会出现“无法将表达式转换为浮点”。（请原谅我写代码时的混乱！）我删去了符号a11-a88的定义，以节省阅读，但它们都是A型cmath.exp公司（bx），A1*cmath.exp公司（1.0j*b1*x）或1.0j*A2*cmath.exp公司（1.0j*b2*x）。我一直使用cmath函数而不是数学(cmath.exp公司不是exp，而且cmath.sqrt公司不是sqrt）。你知道吗

import sys
import math
from scipy import *
from numpy.linalg import *
from sympy import *
import numpy
from sympy.solvers import solve
import cmath
from scipy import optimize

plik=open('solution_e-.txt','w')

#I cut-off definitions of symbols a11-a88.

Det =((a77*a88+(-1.0)*a78*a87)*(a44*a55*a66+a45*a56*a64)+(a76*a88+(-1.0)*a78*a86)*(a44*a57*a65+a45*a54*a67))*(a11*(a22*a33+(-1.0)*a23*a32)+a21*(a13*a32+(-1.0)*a12*a33))+((a77*a88+(-1.0)*a78*a87)*(a34*a56*a65+a35*a54*a66)+(a76*a88+(-1.0)*a78*a86)*(a34*a55*a67+a35*a57*a64))*(a11*(a22*a43+(-1.0)*a23*a42)+a21*(a13*a42+(-1.0)*a12*a43))+((a77*a88+(-1.0)*a78*a87)*(a44*a56*a65+a45*a54*a66)+(a76*a88+(-1.0)*a78*a86)*(a44*a55*a67+a45*a57*a64))*(a11*(a23*a32+(-1.0)*a22*a33)+a21*(a12*a33+(-1.0)*a13*a32))+((a77*a88+(-1.0)*a78*a87)*(a34*a55*a66+a35*a56*a64)+(a76*a88+(-1.0)*a78*a86)*(a34*a57*a65+a35*a54*a67))*(a11*(a23*a42+(-1.0)*a22*a43)+a21*(a12*a43+(-1.0)*a13*a42))

equat = Det.real + Det.imag


for i in range (76500,76550,1):
    n=i/100000.0

    equat_lam = lambdify(x,equat)
    Solut = optimize.fsolve(equat_lam, n)
    plik.write(str(float(Solut))+'\n')
    print n



plik.close()

编辑：错误的完全回溯

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
C:\Anaconda\lib\site-packages\IPython\utils\py3compat.pyc in execfile(fname, glob, loc)
    195             else:
    196                 filename = fname
--> 197             exec compile(scripttext, filename, 'exec') in glob, loc
    198     else:
    199         def execfile(fname, *where):

C:\Users\Melania\Documents\doktorat\2017\analiza\Próbka I\poziomy_en\rozwiazanie_elektrony.py in <module>()
     33 print 'I defined other symbols'
     34 
---> 35 k1=(cmath.sqrt(2.0*(V1-x)*m))/hkr
     36 k2=(cmath.sqrt(2.0*(V2-x)*m))/hkr
     37 k3=(cmath.sqrt(2.0*x*m))/hkr

C:\Anaconda\lib\site-packages\sympy\core\expr.pyc in __complex__(self)
    210         result = self.evalf()
    211         re, im = result.as_real_imag()
--> 212         return complex(float(re), float(im))
    213 
    214     @_sympifyit('other', False)  # sympy >  other

C:\Anaconda\lib\site-packages\sympy\core\expr.pyc in __float__(self)
    205         if result.is_number and result.as_real_imag()[1]:
    206             raise TypeError("can't convert complex to float")
--> 207         raise TypeError("can't convert expression to float")
    208 
    209     def __complex__(self):

TypeError: can't convert expression to float

Answer 1

痕迹以

C:\Users\...
k1=(cmath.sqrt(2.0*(V1-x)*m))/hkr

最后你看到一个

TypeError: can't convert expression to float

由symphy的expr.__float__引发，被expr.__complex__调用，因此可以推断表达式2.0*(V1-x)*m不能转换为复数-通常发生这种情况是因为它包含一个自由符号。你知道吗

如果你想用数值计算平方根，你必须用一个数值来代替构成cmath.sqrt参数的所有符号，其中每个项，例如V1，都可以是一个包含大量符号的符号表达式。你知道吗

也就是说，如果你想找到函数实部和虚部都为零的值，显然你不应该写equat = Det.real + Det.imag