shenyongjun 发表于 2008-6-10 16:48

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

% %% 程序的目的是寻找序列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 发表于 2008-6-10 17:33

早就有人编过类似程序:a=rand(1,100);
b=find(diff(sign(diff(a)))>0)+1
t=b(1)

shenyongjun 发表于 2008-6-10 17:47

回复 2楼 的帖子

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

bainhome 发表于 2008-6-10 17:57

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

shenyongjun 发表于 2008-6-10 21:11

回复 4楼 的帖子

学习了!:handshake
页: [1]
查看完整版本: 自己编的一个求时间序列中第一个局部极小值的程序,求高手指点