学习一下MATLAB中矩阵的逻辑索引方法

ulppknot

Matlab的计算是以矩阵为单元，有着强大的矩阵处理功能，使用时应该多用矩阵或向量运算，把循环语句转变为向量运算会省很多时间，程序也更简洁易读。

MATLAB中对矩阵元素的索引，包括常用的下标法和巧妙的逻辑索引，逻辑索引是用逻辑0或1去判断对应位置的元素是否被选择，所以矩阵和逻辑阵维数必须一致；是不是可能有点不明白，贴个例子你就懂啦：

>>A=[-1 1;3 -2]
! y3 s2 }$ _0 O; NA =
( l4 F7 i8 H" Q5 ^# ]   -1     1
    3    -2
& j, O9 A; P8 v. n>> A > 0    %得到逻辑矩阵
ans =
    0     1
* N8 o. m$ P' [$ O2 }5 M3 \    1     0
>> A(A > 0)   %逻辑索引,逻辑1对应的元素被选中，按列选择
ans =
    3
, P/ w% Q' L& b3 M+ Q3 Z' y/ x    1
%当然，我知道你想得到的是酱紫的
7 H) F! \- P( U2 wans=
. K, R$ a5 G1 a! n  k9 [1 P, N$ r    0     1
    3     0
%这个可以这样得到:
+ g4 B* ]: A, b; b>>A.*(A > 0)   %此方法就是矩阵运算的好处，可以用来取满足条件的数，不必用循环,注意是点乘


另外，如果你要把矩阵中满足一定条件的数赋新值
8 k5 R4 k4 u1 a9 J& a* y
比方说，一个名叫array数组里面，你要将里面大于1的都变成0，就不必用到循环：
) z" n) }" c& k" V5 c
% M# g& r* \9 x; xarray(array>1)=0;
: `+ P+ o# J  T) w
把大于1小于3的变成0。

array(array>1 & array<3)=0;
3 \& i+ L! Z' s7 g# Y
$ O  m( J- L: K3 m当然，还可以使用find，这个也很好用。
  d, k% ?5 o# K0 ^: x
5 @' W& }! d4 ]: i9 v7 r. d5 m5 t但是，当你要把矩阵中满足一定条件的数做运算时，

比如，把array数组里面大于50的数减去50，其他数不变，如果用上面的等式就不太行啦：
# \/ o- \3 X9 G9 i5 ~
>>array=[51 34;57 23;45 89]
8 a5 K, Q" z/ j- d1 y, warray =
1 b  s% \5 y7 z1 ]% [  f2 s
4 Q4 G  h  N0 p* ?. ]% V    51    34
    57    23
& _. o( E* L) j  {, B- [    45    89
" h: R& q( A$ ]$ D>>array(array>50)-50
9 ^, Y+ N% e4 G4 O  V$ p7 [# `2 N
& ?3 U4 ]- Z8 D9 U( yans =
     1
     7
8 e" _2 O2 m+ A' }7 Z    39
7 w* j( F0 o  |
. i8 C6 Z3 B2 |, ?, |) o这显然不是我们想要的结果，但是怎么会没有解决的办法呢，MATLAB这么强大，看下面
8 o4 J0 k3 h& b( c- k
>>array(array>50)=array(array>50)-50
array =
    1    34
, t% A3 w9 N% s' F/ [/ e# J$ E1 a    7    23
    45   39

这个问题在我毕设的时候也困扰了我很久，如果用循环的话，未免low了点，这个方法就很实用。
  ^5 b1 k5 Q# K5 n  l
7 m4 Y1 \! y$ g" _7 ?
, T- W' ^& ~  |  {

xixihahaheihei

谢谢分享