[请益] 积分 曲线长度

楼主: 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