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

ulppknot

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

) ]7 I3 f5 d5 ]% d! X5 {& |7 HMATLAB中对矩阵元素的索引，包括常用的下标法和巧妙的逻辑索引，逻辑索引是用逻辑0或1去判断对应位置的元素是否被选择，所以矩阵和逻辑阵维数必须一致；是不是可能有点不明白，贴个例子你就懂啦：
  O6 B5 ^# p8 h) Q" B5 S3 \$ L% s
>>A=[-1 1;3 -2]
A =
; w( d0 r  R9 j) ^   -1     1
' o: w3 }8 e& r1 S2 q: ?) }    3    -2
>> A > 0    %得到逻辑矩阵
ans =
6 g- V0 S, ]$ n5 L4 N    0     1
    1     0
>> A(A > 0)   %逻辑索引,逻辑1对应的元素被选中，按列选择
4 ^! R1 M  a& \+ _7 Rans =
    3
4 `7 c* t5 D# w1 Z4 ?    1
%当然，我知道你想得到的是酱紫的
: i; M7 N+ o3 \. Ians=
4 B4 r- q- b! Q, h$ Z    0     1
) X' M) e% u0 u# ]! A/ I    3     0
%这个可以这样得到:
/ X$ X9 g+ _; F9 e>>A.*(A > 0)   %此方法就是矩阵运算的好处，可以用来取满足条件的数，不必用循环,注意是点乘
8 P9 O, c& p4 p1 D3 v; t5 ]
! J& @' T+ R# o8 D4 B/ v# _
, e8 m" x/ H' _" q5 F+ Q5 s另外，如果你要把矩阵中满足一定条件的数赋新值
% C1 q% Y3 L2 A. p6 [0 w
比方说，一个名叫array数组里面，你要将里面大于1的都变成0，就不必用到循环：
, b, a' O: n1 T( M: w
array(array>1)=0;
/ `9 Z6 ?; z9 c$ a
3 `; g# _- S5 l* [( I, e把大于1小于3的变成0。
4 ~. H- Z$ n& W  a
7 U. W! m, [" r; A7 t2 R+ Sarray(array>1 & array<3)=0;

当然，还可以使用find，这个也很好用。
3 {4 [, F# `9 b0 K
但是，当你要把矩阵中满足一定条件的数做运算时，
. Y2 E( G3 z$ }
1 R# s! f+ ^7 S7 E$ Z6 u比如，把array数组里面大于50的数减去50，其他数不变，如果用上面的等式就不太行啦：

- i7 R8 W' Q7 d  P! v( X! l >>array=[51 34;57 23;45 89]
array =
( _* i  F6 I; a$ |& k1 S4 V0 b
    51    34
9 n7 [+ J- Z0 b* v) M+ C: g" C    57    23
    45    89
' ]  Q2 z" F1 T! N  v! @>>array(array>50)-50

5 Z& x3 K. A; E" L3 f- tans =
     1
' K! E0 K$ T* H4 q$ u0 f% C     7
3 X  Z9 T! Z. N/ Y# ]    39

7 M; V8 h5 b9 m这显然不是我们想要的结果，但是怎么会没有解决的办法呢，MATLAB这么强大，看下面
0 F$ \7 C3 p8 K) x
>>array(array>50)=array(array>50)-50
array =
    1    34
    7    23
8 a7 v2 f$ t/ J    45   39
6 Z$ t2 C0 k4 k. n& x
7 Q$ H& e7 `1 }# D/ K. S这个问题在我毕设的时候也困扰了我很久，如果用循环的话，未免low了点，这个方法就很实用。
5 t9 i9 S9 x# h5 f- j& t
# f4 z6 F% W* s" M$ t$ T- x' {) }0 S
- x/ J; W0 L0 B8 a3 ~9 S1 f

xixihahaheihei

谢谢分享