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

ulppknot

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

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

: g0 L9 k5 N4 F% W% S( A, p# k>>A=[-1 1;3 -2]
2 a8 P! I4 u  t1 A( u  B* ^A =
! p1 T" V% ?' B( \& f5 C1 r: C   -1     1
2 y6 m+ [) g" @5 j  s    3    -2
>> A > 0    %得到逻辑矩阵
ans =
    0     1
% b: m. A1 x6 P8 S5 b    1     0
0 z" j4 g+ p3 Q- A/ V>> A(A > 0)   %逻辑索引,逻辑1对应的元素被选中，按列选择
ans =
    3
    1
%当然，我知道你想得到的是酱紫的
/ F8 R' z$ p+ P/ B4 c/ Lans=
% h8 ~! Y7 K9 ]! M: y8 `, R' b    0     1
( |5 J# D) ^' d    3     0
1 Y3 J# [6 P$ y( M+ y  f%这个可以这样得到:
>>A.*(A > 0)   %此方法就是矩阵运算的好处，可以用来取满足条件的数，不必用循环,注意是点乘
  B8 _2 N# F: D- b

. F& q! m2 N* _% d另外，如果你要把矩阵中满足一定条件的数赋新值

比方说，一个名叫array数组里面，你要将里面大于1的都变成0，就不必用到循环：

array(array>1)=0;
; Z& T/ B7 k& u' G' f9 h& c2 ?, o* G
把大于1小于3的变成0。
$ p6 G+ V' a/ E8 [7 ^; V
5 S' u/ H3 c7 Zarray(array>1 & array<3)=0;

当然，还可以使用find，这个也很好用。

3 {5 d0 f9 z6 A, p( a2 l3 m  T但是，当你要把矩阵中满足一定条件的数做运算时，

8 l, r$ b  G' E比如，把array数组里面大于50的数减去50，其他数不变，如果用上面的等式就不太行啦：
& q) ]6 E, S' E' v
3 S9 v2 n2 V# F5 `/ I >>array=[51 34;57 23;45 89]
array =
0 s  f6 z. T% z
    51    34
    57    23
    45    89
>>array(array>50)-50
4 F1 ?8 J) ?8 P+ `0 ]" P; i7 x
, l2 P: {- [$ O' X5 x) rans =
9 ~8 H2 T& ?& f1 H  K. w     1
     7
4 B  D* B: H" Z& ~$ |8 o4 [. [+ r    39

! s) D1 W2 y6 p1 c6 f9 w0 b这显然不是我们想要的结果，但是怎么会没有解决的办法呢，MATLAB这么强大，看下面
5 d" W( a2 P; _( w6 M' l4 y
( ~; O& \  k! M% p) D+ t/ q# j/ J>>array(array>50)=array(array>50)-50
' N8 H; \  d$ t* E  c% Narray =
" o& l/ D! J2 M5 g6 e- @    1    34
: Q" }! l- q# l! W* D! Z$ Y$ J    7    23
/ I/ F7 m( P3 M# g( F6 l. @1 ~; [    45   39
" O$ M' x" M# C& H9 [
这个问题在我毕设的时候也困扰了我很久，如果用循环的话，未免low了点，这个方法就很实用。
1 `" V, a) L$ G
7 O/ l9 x. m2 P$ ^; Y

xixihahaheihei

谢谢分享