求极值点

mst-tum

max 和min只能求出向量的最值点，有什么办法找出极值点吗?
也就是说要求x(i)<x(i+1)并且x(i)<x(i-1)的点
我试过用循环，不过结果不理想，每次得到的极值点都不同
谢拉！先

[ 本帖最后由 mst-tum 于 2007-5-4 20:08 编辑 ]

xjzuo

不妨将你的程序贴一下.
另:可以用现成的命令:fminunc, or fminsearch.
一般方法为：先画图，再在附近搜索极值点．

[ 本帖最后由 xjzuo 于 2007-5-4 20:17 编辑 ]

eight

原帖由 mst-tum 于 2007-5-4 20:01 发表

                               
登录/注册后可看大图

max 和min只能求出向量的最值点，有什么办法找出极值点吗?
也就是说要求x(i)

可参考 emd 程序中的 extr.m，即：

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  
    

复制代码

mst-tum

:lol ，谢谢拉。大家都好勤劳，假期还在。本来没指望这么快有回复的。
我先学习下，不懂再请教。

mst-tum

学习了以上的程序，还是有些地方不懂。
这里第27行的m指的是什么啊?
是不是写错了，应该是n吧?

eight

原帖由 mst-tum 于 2007-5-14 21:50 发表

                               
登录/注册后可看大图

学习了以上的程序，还是有些地方不懂。
这里第27行的m指的是什么啊?
是不是写错了，应该是n吧?

m = length(x);

其实仔细想想就猜到应该是原信号的程度

hyyly

居然不会用 ，到底返回什么值呢？

Happy99

hyyly 发表于 2010-11-5 08:18

                               
登录/注册后可看大图

...到底返回什么值呢？

什麼问题?
个人水平有限, 建议说清楚些

qibbxxt

这个不是用

1. imregionalmax

复制代码

和

1. imregionalmin

复制代码

就可以了吗？

hyyly

回复 Happy99 的帖子

我的意思是 eight  发的代码我不会使用，不过没关系，我找到更好的了。谢谢

zhouyang664

寻找极大值：
% from MATLAB-Wiki-FAQ
function index = localmax(x)
index = find( diff( sign( diff([0; x(:); 0]) ) ) < 0 );
LMax = x(index)
另外请参考副教授花如月的经典之作：
http://forum.vibunion.com/forum- ... A2%D0%F2%C1%D0.html