|  | 
 
| 之前,在论坛里帮忙别人解决了一个问题。找一条曲线的拐点。 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')
 
 | 
 |