matlab 图像 局部极值,[转载]matlab 图像局部求极值

求极值的具体语句：

极值：

data是你的数据，

find(diff(sign(diff(data)))==-2)+1

找到极大值的位置

find(diff(sign(diff(data)))==2)+1

找到极小值的位置

data(find(diff(sign(diff(data)))==-2)+1)和

data(find(diff(sign(diff(data)))==2)+1)

返回的是极大值和极小值

要想diff(sign(diff))出现-2，只有是sign的结果为-1,1排列的时候，-1代表，前一次的diff运算前面的数值比后面大，1代表，后一次diff运算后面的数值比前面大，因而，x-a,x,x+b，a,b均为正数，则x为极小值。

同理判断极大值。

find查找到极值的索引后还要+1，是因为，差分运算(diff)后的数据的长度比原始数据少1，见help。

%方法一

x=0:0.01:20;

y=2*sin(x/2)+cos(2*x)/2;

% indmax=find(diff(sign(diff(y)))<0)+1;%极大值点

% indmin=find(diff(sign(diff(y)))>0)+1;%极小值点

ind=find(diff(sign(diff(y)))~=0)+1;%极值点

plot(x,y,x(ind),y(ind),'ro')

%方法二

x=0:0.01:20;

y=2*sin(x/2)+cos(2*x)/2;

indmax=intersect(find(diff(y)>0)+1,find(diff(y)<0));%极大值点

indmin=intersect(find(diff(y)<0)+1,find(diff(y)>0));%极小值点

plot(x,y,x(indmax),y(indmax),'ro',x(indmin),y(indmin),'go')

%方法三

x=0:0.01:20;

y=2*sin(x/2)+cos(2*x)/2;

yf=y(2:end-1)-y(1:end-2); %前向差分

yb=y(2:end-1)-y(3:end); %后向差分

I=find(yf.*yb>=0)+1; %寻找极值点(不考虑首尾两点)

plot(x,y,'b',x(I),y(I),'r*');