自己编的一个求时间序列中第一个局部极小值的程序，求高手指点

shenyongjun

% %% 程序的目的是寻找序列x中的第一个极小值的脚标
%%
%%% 程序有些繁琐，但是可以用
%%%


function tau=my_weizhi(x);

n=length(x);

x1=x(1);
x2=x(2);
x3=x(3);

nn=0;
for nnn=1:(n-3)
   
if x1 > x2 & x3 > x2
    nn=nn+1;
    tau1(nn)=x2;
end
    x1=x(nnn+1);
    x2=x(nnn+2);
    x3=x(nnn+3);
end

tau=find(x==tau1(1));

[ 本帖最后由 sigma665 于 2008-6-10 18:46 编辑 ]

bainhome

早就有人编过类似程序：

1. a=rand(1,100);
   
2. b=find(diff(sign(diff(a)))>0)+1
   
3. t=b(1)

复制代码

shenyongjun

我在计算过程中是先找到所有的极小值，然后返回第一个极小值的脚标。现在我是想让高手指点一下，看看能否不计算其他的极小值，而只计算到第一个极小值点即返回。这样不是可以节约时间吗？
感觉楼上的方法好像也是把所有的极小值点都找出来了！

bainhome

矩阵运算是MATLAB的推荐运算方式，恰恰节约而不是浪费时间。在这个问题上，你想当然了。该代码基本已经成为求极小(大)值的标准程序。不但达到还超过了你的要求，我认为是合适的。
假如只寻找第一个极小值，你的代码应该是唯一方式，无非就是循环移位，如果能够找到第n个极小值，也就是这种方式，估计谁指点也一样:victory:

shenyongjun

学习了！:handshake