0

我正在编写一个简短的程序，它应该找到函数的实部和虚部都为零的值。我不明白为什么我在运行程序后得到“无法将表达式转换为浮点数”。（请原谅我写代码时的乱七八糟！）我切断了符号a11-a88的定义，以节省您的阅读时间，但它们都是A类型cmath.exp(b x), 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
4

1 回答 1

0

跟踪开始于

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

最后你会看到一个

TypeError: can't convert expression to float

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

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

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

于 2018-11-03T01:26:40.193 回答