自己编的一个求时间序列中第一个局部极小值的程序,求高手指点
% %% 程序的目的是寻找序列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 编辑 ] 早就有人编过类似程序:a=rand(1,100);
b=find(diff(sign(diff(a)))>0)+1
t=b(1)
回复 2楼 的帖子
我在计算过程中是先找到所有的极小值,然后返回第一个极小值的脚标。现在我是想让高手指点一下,看看能否不计算其他的极小值,而只计算到第一个极小值点即返回。这样不是可以节约时间吗?感觉楼上的方法好像也是把所有的极小值点都找出来了! 矩阵运算是MATLAB的推荐运算方式,恰恰节约而不是浪费时间。在这个问题上,你想当然了。该代码基本已经成为求极小(大)值的标准程序。不但达到还超过了你的要求,我认为是合适的。
假如只寻找第一个极小值,你的代码应该是唯一方式,无非就是循环移位,如果能够找到第n个极小值,也就是这种方式,估计谁指点也一样:victory:
回复 4楼 的帖子
学习了!:handshake
页:
[1]