作者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/m.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