matlab离散数据拐点,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;matlab求导
b=y(a(end));                                                              %上升趋势的最后⼀个点的纵坐标
b1=y(a1(1));                                                              %下降趋势的第⼀个点的纵坐标
c1=x(a1(1));                                                              %下降趋势的第⼀个点的横坐标