作者sam122094 (sam)
看板Python
标题[问题] matplotlib在数据下积分
时间Sun Mar 13 02:10:57 2016
就是小弟目前有一笔实验数据
但没有其方程式
用matplotlib做出图大概长这样
http://imgur.com/WhLEnjV
由於没有方程式，不晓得该如何积分＠＠
我主要是要找一条铅直线
使得左边曲线以下面积跟右边曲线以上面积一样大
但遇到积分就卡了
google大部分都是要先有f(x)才能积分做图的
请问有没有办法呢？
谢谢大家
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.4.209
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Python/M.1457806260.A.543.html
1F：推 stan622: 哈哈 大胆推测这是物冶作业XD03/13 02:16
2F：推 stan622: 可以试试scipy.optimize.curvefit把方程式逼近出来後，03/13 02:51
3F：→ stan622: 用个while回圈之类的把matano interface 找出来03/13 02:51
4F：推 stan622: 也可以考虑整条曲线下面积=matano interface 右边大矩形03/13 03:05
5F：→ stan622: 面积，来求答案03/13 03:05
终於有关键字了！Google 都得到bumpy,scipy 底下的一堆class，根本无法下手QQ，感谢
stan大，它的确是物冶哈哈
6F：推 PhysiAndMath: 数值积分+辛普森法 喂狗03/13 15:33
辛普森法不用有原先方程式吗？
※ 编辑: sam122094 (101.13.160.29), 03/13/2016 16:32:01
7F：推 orzguy: scipy.integrate的Integrating using Samples 03/13 22:50
8F：推 painkiller: numpy.trapz 把每个数据点都积一遍 03/17 02:15
9F：→ painkiller: 然後numpy.interp1d找出50%面积的x值 03/17 02:15
10F：→ painkiller: 当然数据量如果很大你就binary search去猜吧 03/17 02:16