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

ulppknot

Matlab的计算是以矩阵为单元，有着强大的矩阵处理功能，使用时应该多用矩阵或向量运算，把循环语句转变为向量运算会省很多时间，程序也更简洁易读。
- p# u( P  F3 p" [5 \# w7 @' G
MATLAB中对矩阵元素的索引，包括常用的下标法和巧妙的逻辑索引，逻辑索引是用逻辑0或1去判断对应位置的元素是否被选择，所以矩阵和逻辑阵维数必须一致；是不是可能有点不明白，贴个例子你就懂啦：
9 K; B/ Z* f# j% Q7 _) X' j. E; z
2 ]5 v& h* l" y: i3 j5 |>>A=[-1 1;3 -2]
* T, U4 P! e' `0 M# qA =
   -1     1
0 ?4 d( n6 F+ c    3    -2
>> A > 0    %得到逻辑矩阵
ans =
" A+ \9 @, q7 T    0     1
+ V6 W9 m, n. @8 ^0 L4 Q    1     0
>> A(A > 0)   %逻辑索引,逻辑1对应的元素被选中，按列选择
9 j* A, _6 @3 [$ pans =
0 t+ i% v/ {: `: Q1 u    3
    1
+ J3 c5 f( b# l0 @( B%当然，我知道你想得到的是酱紫的
ans=
    0     1
    3     0
1 {1 Q8 l5 a, d. `%这个可以这样得到:
0 E/ r8 a6 ~, |% k- e& [>>A.*(A > 0)   %此方法就是矩阵运算的好处，可以用来取满足条件的数，不必用循环,注意是点乘
- |7 }& d0 h; U9 j* R, R
/ x; k  p0 E" o' h& l, T
9 q- q* o; M3 C1 \另外，如果你要把矩阵中满足一定条件的数赋新值
0 Z! G& |7 [7 y  `6 v% \$ G1 k: J
比方说，一个名叫array数组里面，你要将里面大于1的都变成0，就不必用到循环：
8 I6 C8 h* c4 Y: u# J; [( ~
% w& y; n& u( \7 n+ O" X5 O7 i+ Farray(array>1)=0;
2 Y, P: X7 ^6 G& F: b0 E6 f0 b3 `
- h0 ^$ r5 z7 y' i" z把大于1小于3的变成0。

* D2 c. M( T& darray(array>1 & array<3)=0;

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

但是，当你要把矩阵中满足一定条件的数做运算时，

, K+ A1 ^: H/ w比如，把array数组里面大于50的数减去50，其他数不变，如果用上面的等式就不太行啦：
" }  r- t9 v( N" r" [
>>array=[51 34;57 23;45 89]
array =
* n/ F7 U9 Q5 `/ a1 ~
    51    34
; P) g, y- ~3 X: v; X6 b    57    23
    45    89
>>array(array>50)-50

ans =
. N7 i( g" K& I! C& ?: {  C; @     1
     7
    39
: O: |/ e& Q8 ~7 U+ w' F
6 c8 {6 b" K! T8 K# q0 n5 B% O# S这显然不是我们想要的结果，但是怎么会没有解决的办法呢，MATLAB这么强大，看下面

. w# Y5 C/ L; [. ?1 ~) q>>array(array>50)=array(array>50)-50
5 k  n8 v& a4 k1 o  u" c( |2 Uarray =
    1    34
    7    23
0 N0 P" Y; w0 O$ ]" I& [8 U; D    45   39

; W, n3 Y* V, R) F这个问题在我毕设的时候也困扰了我很久，如果用循环的话，未免low了点，这个方法就很实用。
3 K% T! l" i1 n0 f! B3 {


( o; j9 N( U& S6 B

xixihahaheihei

谢谢分享