请教取数组中次极大(小)的简单方法

lon91ong

数组中的极大(小)值和位置可以用maxval (minval) 和maxloc (minloc)得到, 那么次极大(小)要如何取呢?
我记得上面的几个函数是可以接受掩码参数的, 但是具体这个掩码怎么生成不太清楚

麻烦大侠能给出代码示例, 多谢了

mofangbao

mask就是一个逻辑运算 比如你找到了  最大值是50  数组是a ，a<50就是一个mask，这样再找最大值应该就能找到第二大的了。我也没试过 你试完了可以来反馈一下~

言深深

写程序处理，在找出来极大值之后可以通过赋值的方法去掉极大值，然后从剩下来的里面再挑。

lon91ong

mofangbao 发表于 2014-8-30 09:02
mask就是一个逻辑运算 比如你找到了  最大值是50  数组是a ，a

参考版主的意见, 做了一个测试代码, 如下:

1. 
   
2. program main
   
3.     implicit none
   
4.     integer :: n
   
5.     real(4), allocatable, dimension(:,:) :: A
   
6.     n=8
   
7.     allocate(A(n,n))
   
8.     A(1,:)=(/24.,64.,8.0,6.0,2.,6.,12.,26./)
   
9.     A(2,:)=(/12.,25.,5.0,2.0,6.,9.,12.,20./)
   
10.     A(3,:)=(/42.,86.,82.,74.,120.,45.,42.,80./)
    
11.     A(4,:)=(/38.,80.,42.,14.,18.,49.,21.,35./)
    
12.     A(5,:)=(/36.,78.,62.,54.,54.,40.,29.,96./)
    
13.     A(6,:)=(/96.,70.,52.,34.,38.,65.,12.,46./)
    
14.     A(7,:)=(/60.,120.,12.,2.,40.,8.,12.,24./)
    
15.     A(8,:)=(/80.,30.,24.,4.0,6.,16.,36.,56./)
    
16.     !call matrixOutput(n,A,6)
    
17.     print *,"Maxval of A:",maxval(A)
    
18.     print *,"Maxloc of A:",maxloc(A)
    
19.     print *,"Maxval2 of A:",maxval(A,A<120)
    
20.     print *,"Maxloc2 of A:",maxloc(A,A<120)
    
21. end program

复制代码

结果如下:

1. 
   
2. Maxval of A:   120.000000
   
3. Maxloc of A:           7           2
   
4. Maxval2 of A:   96.0000000
   
5. Maxloc2 of A:           6           1

复制代码

现在的问题是, 如果数组中有两个相同的极大值, 比如示例代码中有两个120, 我希望取到的第二个极大值是另外一个120
也就是说, 我希望是按位置取掩码, 而不是按值取掩码
请教大侠实现方法, 多谢了

mofangbao

那就用深深说的方法把最大值赋值为一个小值吧

lon91ong

mofangbao 发表于 2014-8-31 09:08
那就用深深说的方法把最大值赋值为一个小值吧

看来也只能如此了