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

ulppknot

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

MATLAB中对矩阵元素的索引，包括常用的下标法和巧妙的逻辑索引，逻辑索引是用逻辑0或1去判断对应位置的元素是否被选择，所以矩阵和逻辑阵维数必须一致；是不是可能有点不明白，贴个例子你就懂啦：
9 S! q/ k5 h& R
3 n+ b5 z! M# M! Y3 i" T9 R) S>>A=[-1 1;3 -2]
7 z) J8 V5 H1 y, E* ~0 @A =
& {) o; X9 }# G# z2 |/ U; }% d! G   -1     1
    3    -2
' D1 |  {$ h- O- K, K0 x>> A > 0    %得到逻辑矩阵
ans =
    0     1
    1     0
1 b8 I  E3 Y8 I' _- ]>> A(A > 0)   %逻辑索引,逻辑1对应的元素被选中，按列选择
ans =
    3
    1
%当然，我知道你想得到的是酱紫的
ans=
! ^9 ?5 k  u: X; g9 ^7 E  [1 D    0     1
+ Q5 E* D2 e9 Z    3     0
1 v! s; f1 r* p& L4 ^* Y' h; m%这个可以这样得到:
: Q7 ?! i, h" ?2 `  i. p4 N4 i>>A.*(A > 0)   %此方法就是矩阵运算的好处，可以用来取满足条件的数，不必用循环,注意是点乘


, a7 l0 B) h3 I3 P! k7 c+ c另外，如果你要把矩阵中满足一定条件的数赋新值

4 A  R! r; ?9 a: R* t. y0 C" K比方说，一个名叫array数组里面，你要将里面大于1的都变成0，就不必用到循环：
; n8 U3 b( r, {" P
& e0 B) y- B3 G! E% ]array(array>1)=0;
, s3 O( R7 z" c3 ?4 k  y
! R2 h7 K- o% D1 q% H( G0 i把大于1小于3的变成0。

array(array>1 & array<3)=0;
1 u% T) o, d8 ^
: U; l2 T( [5 Z! U- x0 ^当然，还可以使用find，这个也很好用。
8 |& T5 X6 C( Y: X& F$ u
2 \: O: Q* E8 A" a6 }! ^' a但是，当你要把矩阵中满足一定条件的数做运算时，

3 m$ w+ ?4 O' U" K比如，把array数组里面大于50的数减去50，其他数不变，如果用上面的等式就不太行啦：

>>array=[51 34;57 23;45 89]
array =

) O2 o% f8 q) E% |    51    34
5 K9 ?1 F4 u. a/ O" e8 l' j: v% m    57    23
    45    89
>>array(array>50)-50
8 [7 `3 ~# H, \8 Z5 ]
ans =
9 d& c# |8 N$ k) [# k2 X     1
; c. g+ m3 n: S$ O4 m/ y! i     7
! M; u3 k) q4 a5 }4 s    39
, f/ ]5 a( L+ c9 K9 w
这显然不是我们想要的结果，但是怎么会没有解决的办法呢，MATLAB这么强大，看下面

>>array(array>50)=array(array>50)-50
5 n: b1 _5 h6 J3 |9 ]/ h( w5 P& harray =
+ v0 D1 G# F0 `- k8 U    1    34
4 f* K+ n9 c' x1 h2 l/ G; P    7    23
! A3 c% I% S4 J0 j  G" z    45   39
% E1 g, E8 C" r* \4 a
0 w/ a9 m0 X$ J5 x* L这个问题在我毕设的时候也困扰了我很久，如果用循环的话，未免low了点，这个方法就很实用。
' L0 [- Y' |# L
/ U! o+ D' a( M+ V* B/ Y7 A
* E+ m; Z! t1 P+ v2 s/ G; D
" L! n- ^6 h% z$ Q

xixihahaheihei

谢谢分享