【问题标题】:Python: Simpson's rule, negative answer for positive area under the curve?Python:辛普森规则,曲线下正面积的否定答案?
【发布时间】:2016-08-16 16:05:07
【问题描述】:

我有

y1=[  9.49110000e-004   4.74145420e-004   1.41847155e-008   3.33228420e-028
       3.76352289e-081   4.48206815e-193   0.00000000e+000   0.00000000e+000
       0.00000000e+000   0.00000000e+000   0.00000000e+000] 

x=[ 112.  111.97667396  111.90666665  111.78989038  111.62619837  
     111.41538384  111.15717866  110.85125168  110.49720654  110.09457901
    109.64283388]

所有y 值都是正数,因此曲线下的面积应该是正数。但是当我尝试使用辛普森法则整合时

  from scipy.integrate import simps
  b= simps(y1, x)

  print b

我得到-2.45630795891e-05 作为答案。我究竟做错了什么?

【问题讨论】:

    标签: python scipy integrate simpsons-rule


    【解决方案1】:

    x 正在减少,因此您正在从右向左积分。这会翻转结果的符号。

    【讨论】:

    • 我把y1x都颠倒了,还是一样的答案。
    • @HuShu:嗯,when I reverse them,答案是肯定的。
    • 哇!谢谢。我不知道这是如何/为什么发生的,但是在我重新启动内核并颠倒顺序后,它给了我一个正数。
    猜你喜欢
    • 2016-02-16
    • 2016-05-05
    • 2014-12-17
    • 2013-12-21
    • 2015-12-15
    • 2017-12-04
    • 2015-02-13
    • 2020-01-19
    • 2019-04-25
    相关资源
    最近更新 更多