你好,我对 python 很陌生,我正在尝试解决一组方程。我正在尝试将四积分的结果乘以浮点数,我的代码是:
from __future__ import division
import scipy.special as sp
from scipy import integrate
import math
Tb = 7.2
Tc = 9.3
t = Tb / Tc
n = 2*10**-6
L = 50*10**-6
W = 80*10**-9
a = 3*10**-2
s1 = W/ (2*n)
y1 = (L+(W/2)) / (2*n)
x0 = 0.015
r0 = 4*x0
s2 = (r0 / n)/1000000
print s2
y0 = (x0 / n)/1000000
def t1(t):
return t**-1*sp.kv(0, s2)
def t2(t):
return t**-1*sp.iv(0, s2)
print t2
Fk2 = (math.pi**-2) * integrate.quad(t1, s1, s2)
FI2 = (math.pi**-2) * integrate.quad(t2, s1, s2)
print Fk2 , FI2
我不断收到错误消息
25 print t2
---> 26 Fk2 = (math.pi**-2) * integrate.quad(t1, s1, s2)
27 FI2 = (math.pi**-2) * integrate.quad(t2, s1, s2)
28 print Fk2 , FI2
TypeError: can't multiply sequence by non-int of type 'float'
我很不确定要做什么,我尝试用 lambda 函数替换
def t1(t)
,但同样什么也没做。非常感谢任何帮助,并提前感谢您。
由于您是初学者,我首先要说的是,每当您遇到不理解的错误时,请将该行分成多个单独的步骤,以便您可以准确地找出错误所指的内容,如下所示:
temp = integrate.quad(t1, s1, s2)
print "temp:", temp
Fk2 = (math.pi**-2) * temp
然后你可以看到错误在于 temp 是一个元组,并且没有定义如何将元组乘以浮点数。
Fk2 = (math.pi**-2) * np.array(temp[:2])
您可以在文档中查看每个术语的含义。
integrate.quad 默认返回一个长度为 2 的元组,第一个条目是答案,第二个条目是答案误差的估计。要在进一步计算中使用结果,请尝试
Fk2 = (math.pi**-2) * integrate.quad(t1, s1, s2)[0]
FI2 = (math.pi**-2) * integrate.quad(t2, s1, s2)[0]
这仅选择答案并将其乘以浮点数。
试试这个:
Fk2 = [(math.pi**-2) * e for e in integrate.quad(t1, s1, s2)]
FI2 = [(math.pi**-2) * e for e in integrate.quad(t2, s1, s2)]