24小时热门版块排行榜    

查看: 2492  |  回复: 13

NicoleLas

木虫 (正式写手)

引用回帖:
10楼: Originally posted by 月只蓝 at 2014-12-19 13:41:27
用quad系列函数,点除不能少,否则报错。
quad( @(t) exp( 90*( 1/298 - 1./(28*t + 52) ) ), 6.5, 9 )...

嗯,对嘀,谢谢您的提醒。
Be water, my friend.
11楼2014-12-19 13:43:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
11楼: Originally posted by NicoleLas at 2014-12-19 13:43:03
嗯,对嘀,谢谢您的提醒。...

共同进步。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
12楼2014-12-19 13:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

引用回帖:
12楼: Originally posted by 月只蓝 at 2014-12-19 13:58:13
共同进步。...

嗯,多谢前辈提携。

  奇怪,scipy+numpy算出来的结果怎么不对?
from scipy import integrate
from numpy import exp

def myIntegral(t):
    return exp(90*(1/298 - 1/(28*t + 52)))

result = integrate.quad( myIntegral, 6.5, 9 )
print result


或者:
from scipy import integrate
from numpy import exp
result = integrate.quad( lambda t: exp(90*(1/298 - 1/(28*t + 52))), 6.5, 9 )

结果:
(1.7862780518166446, 1.9831670215094072e-14)

是我Python忘了还是怎么,得做做功课了,看来。
Be water, my friend.
13楼2014-12-19 14:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

引用回帖:
12楼: Originally posted by 月只蓝 at 2014-12-19 13:58:13
共同进步。...

抱歉,打扰了,找到原因了,唉,果然是忘了。

from scipy import integrate
from numpy import exp
result = integrate.quad( lambda t: exp(90*(1.0/298 - 1.0/(28*t + 52))), 6.5, 9 )
print result

结果:
(2.4160828634341525, 2.6823908243868738e-14)
Be water, my friend.
14楼2014-12-19 14:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuyanlong 的主题更新
信息提示
请填处理意见