[請益] 積分 曲線長度

作者: totemist (totem)   2016-03-28 23:17:45
各位大大們好
我想請教該如何使用積分quad/quadl來計算曲線長度的問題
本身data有7000多條曲線,每一條都在X軸上取相同等距的點(X = 0, 6, 12,..., 180)
自己原本是用 L = sqrt((x1 - x0)^2 + (y1 - y0)^2);
但運算速度非常的慢(算一次7000多條曲線個別的長度,大約3~4小時)
上網查了一下,quad 和 quadl 好像要有一條線的函式才能使用。
而我自己這7000多條曲線都只有"點"的資訊,並不是函式。
這樣的情況下我能往哪裡找尋資料呢?
謝謝
作者: sunev (Veritas)   2016-03-28 23:19:00
把for去掉應該就可以了
作者: ksr310135 (馨右小小)   2016-03-28 23:56:00
積函式用quad,積點用trapz另外內插法interp1可以將點轉換成函式
作者: totemist (totem)   2016-03-29 01:21:00
明天立馬來試試看,謝謝。
作者: celestialgod (天)   2016-03-29 09:55:00
trapz可以直接積一個矩陣,每一個column是一條曲線

Links booklink

Contact Us: admin [ a t ] ucptt.com