MATLAB 寻找离散点拐点

最近用matlab拟合离散数据，曲线由两段圆弧和直线组成，如下图所示，现在想要通过代码寻找曲线的拐点也就是开始弯曲的位置，原理很简单，对离散数据求导，直线段的导数比较小，当曲线发生弯曲导数会突然变大，根据选取的阈值来控制选点的精度，阈值越小选取的点越精确。寻找拐点时分为两个部分分别寻找，上升区域和下降区域，m为数据点的个数。

matlab 寻找离散点拐点代码实现：

J= gradient(y,x);                                                          %求导
a=find(abs(J(1,1:floor(m/2)))>0.09);                                       %此处阈值设置为0.09
a1=find(abs(J(1,floor(m/2):end))>0.09)+floor(m/2)-1;
b=y(a(end));                                                               %上升趋势的最后一个点的纵坐标
b1=y(a1(1));                                                               %下降趋势的第一个点的纵坐标
c1=x(a1(1));                                                               %下降趋势的第一个点的横坐标