|
之前,在论坛里帮忙别人解决了一个问题。找一条曲线的拐点。
http://www.ilovematlab.cn/thread-319551-1-1.html
现在,自己也用到了。哈。记录下来。
- function [gu_D,gu_ind] = fun_find_guaidian(D,meps)
- n_D=length(D);
- delt1 = diff(D);
- delt2 = diff(delt1);
- judm1 = delt1(1:(n_D-2)).*delt1(2:n_D-1)>=0;
- judm2 = abs(delt2)>meps;
- gu_ind = find(judm1 & judm2) +1;
- gu_D = D(gu_ind);
- end
复制代码
- % D= [1 1.1 1.2 1.1 1.05 1.0 0.997 0.995 0.8 0.7 0.6 0.5 0.67 0.72];
- % meps=0.04;
- % [gu_D,gu_ind] = fun_find_guaidian(D,meps)
- %
- % hold on
- % plot(D,'.-')
- % plot(gu_ind,gu_D,'ro')
复制代码
|
|