matlab求曲线长度

曲线段在上的弧长为采用积分
所求弧长s=∫√(1+y'²)dx
matlab求出各点的导数，然后按照上式积分

clear
>> x=1:0.1:10;
>> y=rand(1,length(x));
>> dy=diff(y);
>> S=0.1*trapz((1+dy.^2).^0.5)

S =

    9.6474

>> plot(x,y,'o-')