[問題] matplotlib在數據下積分

[sam122094 (sam)]

就是小弟目前有一筆實驗數據

但沒有其方程式

用matplotlib做出圖大概長這樣

http://imgur.com/WhLEnjV

由於沒有方程式，不曉得該如何積分＠＠

我主要是要找一條鉛直線

使得左邊曲線以下面積跟右邊曲線以上面積一樣大

但遇到積分就卡了

google大部分都是要先有f(x)才能積分做圖的

請問有沒有辦法呢？

謝謝大家

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.4.209
※ 文章網址: https://www.ptt.cc/bbs/Python/M.1457806260.A.543.html
推 stan622: 哈哈 大膽推測這是物冶作業XD
推 stan622: 可以試試scipy.optimize.curvefit把方程式逼近出來後，
→ stan622: 用個while迴圈之類的把matano interface 找出來
推 stan622: 也可以考慮整條曲線下面積=matano interface 右邊大矩形
→ stan622: 面積，來求答案
終於有關鍵字了！Google 都得到bumpy,scipy 底下的一堆class，根本無法下手QQ，感謝
stan大，它的確是物冶哈哈

推 PhysiAndMath: 數值積分+辛普森法  餵狗
辛普森法不用有原先方程式嗎？
※ 編輯: sam122094 (101.13.160.29), 03/13/2016 16:32:01
推 orzguy: scipy.integrate的Integrating using Samples
推 painkiller: numpy.trapz  把每個數據點都積一遍
→ painkiller: 然後numpy.interp1d找出50%面積的x值
→ painkiller: 當然數據量如果很大你就binary search去猜吧