如何求一维数据的局部极小（或极大）值

newsgae

求教：用matlab如何实现一组一维数据的局部极小（或极大）值及其对应横坐标的位置或定义域的范围，谢谢！

w89986581

help min
hlep max

yangch

[a,b]=max(y(m:n));
a是值 b是位置
好象是这样

eight

原帖由 newsgae 于 2007-5-11 13:45 发表

                               
登录/注册后可看大图

求教：用matlab如何实现一组一维数据的局部极小（或极大）值及其对应横坐标的位置或定义域的范围，谢谢！

这个问题之前刚讨论过，以后请多搜索一下论坛

1. 
   
2. d = diff(x);
   
3. n = length(d);
   
4. d1 = d(1:n-1);
   
5. d2 = d(2:n);
   
6. indmin = find(d1.*d2<0 & d1<0)+1;
   
7. indmax = find(d1.*d2<0 & d1>0)+1;
   
8. if any(d==0)
   
9.   
   
10.   imax = [];
    
11.   imin = [];
    
12.   
    
13.   bad = (d==0);
    
14.   dd = diff([0 bad 0]);
    
15.   debs = find(dd == 1);
    
16.   fins = find(dd == -1);
    
17.   if debs(1) == 1
    
18.     if length(debs) > 1
    
19.       debs = debs(2:end);
    
20.       fins = fins(2:end);
    
21.     else
    
22.       debs = [];
    
23.       fins = [];
    
24.     end
    
25.   end
    
26.   if length(debs) > 0
    
27.     if fins(end) == m
    
28.       if length(debs) > 1
    
29.         debs = debs(1:(end-1));
    
30.         fins = fins(1:(end-1));
    
31.       else
    
32.         debs = [];
    
33.         fins = [];
    
34.       end      
    
35.     end
    
36.   end
    
37.   lc = length(debs);
    
38.   if lc > 0
    
39.     for k = 1:lc
    
40.       if d(debs(k)-1) > 0
    
41.         if d(fins(k)) < 0
    
42.           imax = [imax round((fins(k)+debs(k))/2)];
    
43.         end
    
44.       else
    
45.         if d(fins(k)) > 0
    
46.           imin = [imin round((fins(k)+debs(k))/2)];
    
47.         end
    
48.       end
    
49.     end
    
50.   end
    
51.   
    
52.   if length(imax) > 0
    
53.     indmax = sort([indmax imax]);
    
54.   end
    
55.   if length(imin) > 0
    
56.     indmin = sort([indmin imin]);
    
57.   end
    
58.   
    
59. end  
    

复制代码

newsgae

谢谢，我搜索了，感觉这个是针对函数形式的

eight

原帖由 newsgae 于 2007-5-12 10:51 发表

                               
登录/注册后可看大图

谢谢，我搜索了，感觉这个是针对函数形式的

你的意思是对于数据不成立？一维数据是连续函数的离散表示，所以应该不存在这个问题

newsgae

好的，谢谢 。
我尝试了另外2个函数， imregionmax,imregionmin, 能够找到极值对应的位置，
谢谢啦。