怎样快速找出任意一组数列中递增数的位置

sigma665

比如：
a=[1 2 3 2 1 2 3]
递增数的位置应该是 1 2 3   5 6 7

下面是我写的一个函数
    function p=posi(a) %找出数列递增项
        NN=length(a);
        for ii=2:NN
            if a(ii)-a(ii-1)>1e-6 %递增
                p(ii)=ii;
                p(ii-1)=ii-1;
            elseif abs(a(ii)-a(ii-1))<1e-6 %相等
                p(ii)=0;
                p(ii-1)=0;
            else
                p(ii)=0;
                p(ii-1)=0;               
            end         
        end        
    end

用这个算出来的结果： 1 2  5 6 7
这个会覆盖掉一部分。

另外，这还是用循环算的，效率较低。

求解答。

qibbxxt

用diff和find试一试吧

sigma665

qibbxxt 发表于 2012-1-13 09:28
用diff和find试一试吧

diff我也想过，但是diff会少一个数

sigma665

1. function p=posi(a) %找出数列递增项
   
2.         NN=length(a);
   
3.         for ii=2:NN
   
4.             if a(ii)-a(ii-1)>1e-6 %递增
   
5.                 p(ii)=ii;
   
6.                 p(ii-1)=ii-1;
   
7.             elseif abs(a(ii)-a(ii-1))<1e-6 %相等
   
8.                 p(ii)=0;               
   
9.             else
   
10.                 p(ii)=0;                        
    
11.             end         
    
12.         end        
    
13.     end

复制代码

不用循环还是不会。。

rocwoods

sigma665 发表于 2012-1-14 13:54
不用循环还是不会。。

indmax=find(diff(sign(diff(a)))<0)+1;%极大值点索引坐标
indmin=find(diff(sign(diff(a)))>0)+1;%极小值点索引坐标
求出来后每个极小值到极大值索引之间就是你要的递增索引项。

qibbxxt

sigma665 发表于 2012-1-13 15:54
diff我也想过，但是diff会少一个数

按照吴兄的思想，来实现：

1. a = [1 2 3 1 2 3];
   
2. b =find(bitor(imregionalmin(a),imregionalmax(a)));
   
3. c = arrayfun(@(x)b(x):b(x + 1),1:2:length(b),'uni',0);
   
4. celldisp(c);

复制代码

1. 
   
2. c{1} =
   
3. 
   
4.      1     2     3
   
5. 
   
6. 
   
7. 
   
8. c{2} =
   
9. 
   
10.      4     5     6

复制代码

bainhome

bitor,imregionalmin,imregionalmax这些命令比较奇怪，从来没用过。
但这段命令的计算结果似乎算的是递减数列，例如下面这段代码：

1. clc
   
2. clear
   
3. data=round(100*rand(1,15))
   
4. b =find(bitor(imregionalmin(data),imregionalmax(data)))
   
5. c = arrayfun(@(x) b(x):b(x + 1),1:2:length(b),'uni',0);
   
6. celldisp(c);

复制代码

1. data =
   
2.     49    44    45    31    51    51    82    79    64    38    81    53    35    94    88
   
3. b =
   
4.      1     2     3     4     7    10    11    13    14    15
   
5. c{1} =
   
6.      1     2
   
7. c{2} =
   
8.      3     4
   
9. c{3} =
   
10.      7     8     9    10
    
11. c{4} =
    
12.     11    12    13
    
13. c{5} =
    
14.     14    15

复制代码

另外向量长度变大后，会出现错误：

1. ??? Attempted to access b(18); index out of bounds because numel(b)=17.
   
2. 
   
3. Error in ==> @(x)b(x):b(x+1)
   
4. 
   
5. 
   
6. Error in ==> celldispxxxxx at 18
   
7. c = arrayfun(@(x) b(x):b(x + 1),1:2:length(b),'uni',0);

复制代码

取极值之间数据计算结果需要考虑极大值、极小值序列索引首个序号大小不同的情况。

qibbxxt

bainhome 发表于 2012-1-16 11:27
bitor,imregionalmin,imregionalmax这些命令比较奇怪，从来没用过。
但这段命令的计算结果似乎算的是递减数 ...

恩，谢谢老兄提出，的确有这些问题，我改进一下

1. clear;clc;close all
   
2. a = [49,44,45,31,51,51,82,79,64,38,81,53,35,94,88];
   
3. bmin = find(imregionalmin(a));
   
4. bmax = find(imregionalmax(a));
   
5. c = cell2mat(arrayfun(@(x)find(bmax > x,1),bmin,'uni',0));
   
6. d = arrayfun(@(x)bmin(x):bmax(c(x)),1:length(c),'uni',0);
   
7. celldisp(d);

复制代码

sigma665

问题是解决了，不过这个程序看不太懂啊。。

比如c的作用是什么？？？