Matlab中数组元素引用——三种方法

uqHZau

1.Matlab中数组元素引用有三种方法
* j0 Z  R7 |7 A8 `; D1 [5 c, q, f
1.下标法(subscripts)
4 P, G0 L7 n* @3 e8 ]2.索引法(index)
6 @3 F, x3 v1 H/ [" J+ }& ?6 p3.布尔法(Boolean)

注意：在使用这三种方法之前，大家头脑一定要清晰的记住，Matlab中数组元素是按列存储(与Fortran一样)，比如说下面的二维数组
  n7 t; T6 v) M) |0 h
A=
     8     1     6
& O9 a$ y. J. q# B; X2 W( K     3     5     7
     4     9     2
) @! @5 n  F, L$ j  q0 d7 }3 S, CMatlab的存储顺序是8,3,4,1,5,9,6,7,2，也就是说先行后列，对于3维数组呢，就是先行后列再页对应个元素的索引和下标分别为
5 V( H# _- O) F2 g+ `
9 p0 n, i& _! OElement Index Subscripts
8        1      (1,1)
3        2      (2,1)
, Y1 s/ n* X) i' w# `4        3      (3,1)
1        4      (1,2)
5        5      (2,2)
9        6      (3,2)
6        7      (1,3)
! Y; e; S* x; J: V7        8      (2,3)
# d  M6 ^- ]& K9 i9 r2        9      (3,3)
+ ^3 ^/ f2 V1 s7 B. e
& S% ~) u, A2 [/ K9 \: T　　从上面的例子中已经很清晰的说明了下标和索引的区别了，也就是说Matlab为每个元素分配了一个唯一识别的ID(即index)
( j" K: n  [. L1 x3 H
9 ^5 T0 a. z6 ?, B: V
9 c4 x7 e5 l: U+ W3 C  Z% @7 g2.下标法引用
- e% K2 {+ M3 G+ A3 B
A(ii,jj)：其中ii和jj可以是一维向量、标量、“:”号或者“end”
+ j4 P8 f5 D' _( [; e& Z5 N
/ i, A$ E) x; x大家对下标估计比较熟悉，由于在C语言中接触过，但是我这里需要强调的是，Matlab的下标是可以多行多列同时引用的，而像C语言等一次只能引用一个，比如
4 }0 d$ V* H# l8 y% f2 p
A(2:3,3:-1:1)表示引用数组中的2~3行，3~1列对应的元素
A(:,end)表示引用最后一列元素，“:”表示所有列或行，“end”表示最后一列或列，“end-n”表示倒数第n行或列
7 F1 C+ J& [" a) RA(1,end-1)表示引用第1行倒数第2个元素
A([2 1 3 3],[1 1 2 2 1])表示引用按两个向量引用指定的元素，即A中的第2,1,3,3行和第1,1,2,2,1列对应的元素


! J* O8 A$ f6 H6 n, d; p举几个例子：
9 F" F# L% n: }8 ?, ]; H
>>A=magic(3)
+ h- q6 I2 Y( \
7 c; C& P# G& K# t/ m; UA =
     8     1     6
& V4 V/ W9 h+ l2 i5 \     3     5     7
     4     9     2

3 g8 m! N& \4 R9 X# G% k* b" s$ x>>A(2:3,3:-1:1)
9 [8 T: T0 r, b$ p6 V, B5 c  ^
6 ]* ]3 W1 U) m% A' `ans =
$ }. I. N- l" c1 p8 Y     7     5     3
     2     9     4
. I1 O5 J+ s; e) q
9 w6 H: n6 v3 n. n4 p3 n  w! Y+ L>>A(:,end)

# Y1 I4 E% y4 v# h# c4 X( yans =
     6
  [% a* o# u" N! u     7
     2

) f+ l% h) j8 \% [! {; D) H>>A(1,end-1)
" C+ N/ h6 G0 D/ X8 x* O
ans =
8 u2 u& r- t+ r' F     1

>>A([2 1 3 3],[1 1 2 2 1])
# h3 H& Q/ Q1 Q5 j' r; J
ans =
  }4 p2 ]& ~+ i     3     3     5     5     3
     8     8     1     1     8
% b+ H& l& ?) C+ ?$ l* d; d/ p     4     4     9     9     4
2 y: x+ ?8 r5 |: i     4     4     9     9     4


, U" {- C, E% `* M# m' z# x2.索引法引用(说白了索引就是存储顺序)

( W1 d$ [! q- cA(index)：index可以是任意的数组，index的元素必须是正整数，且不大于numel(A)，返回的是一个尺寸与index一样的数组

下标和索引之间可以通过ind2sub和sub2ind函数相互转换，具体可以看帮助，很简单
! J6 x* \" C3 ~- i( D4 q+ t, Y( W$ _
[I,J] = ind2sub(siz,IND)
IND = sub2ind(siz,I,J)
还有使用A(：）就可以将数组A转换为列向量

* O& I/ B# C/ A+ D0 t' kA(8)：表示引用A的第8个元素
B=A([1 10 5 2 2 1 3])：表示依次引用A的第1,10,5,2,2,1,3个元素，返回与index尺寸相同的数组，也就是说size(B)=size(index)
/ K6 M$ o; G$ I  E, o' m. RA([2 5 9;1 1 1;8 5 6])：返回的时侯是一个3*3的矩阵

>>A=magic(5)%括号中为索引值
7 j4 M" v  z% o/ q) P/ L
" V3 T& M5 f3 w% j/ ~+ O1 hA =
    17 (1)     24 (6)     1 (11)     8 (16)    15 (21)
- G" q! R+ Q( h8 I; C! H7 K    23 (2)      5 (7)     7 (12)    14 (17)    16 (22)
     4 (3)      6 (8)    13 (13)    20 (18)    22 (23)
    10 (4)     12 (9)    19 (14)    21 (19)     3 (24)
    11 (5)    18 (10)    25 (15)     2 (20)     9 (25)
  \( B" M& G! C- s' m
>>A(8)
3 u7 A' V2 O/ `* [* m& ~$ y
  f7 u" K* _6 F( u5 y& Aans =
     6

7 z9 E4 X4 G! p7 f' v' O>>A([1 10 5 2 2 1 3])
& N" ]% P0 X& h" _3 s( ^# j
ans =
: e) v1 _/ o2 [& e% t    17    18    11    23    23    17     4

- K% d- ?8 A& w* s9 p4 o>>A([2 5 9;1 1 1;8 5 6])

ans =
; s* T2 W& z& B4 z0 G2 Z0 ~    23    11    12
    17    17    17
% R( \, @& ~* A! M: n6 S8 p+ v     6    11    24
　　


3.布尔法引用

A(X)：X是一个有0和1组成布尔型数据，且size(A)=size(X)，对应位置为1则留下该数据，0则去掉，最后按A中的存储顺序，返回一个列向量

假如说A是3*3的数组

A(logical([1 0 0;0 1 0;0 0 1]))：表示引用了数组A的对角线元素，注意必须使用logical将0/1数组转换为布尔型

' u9 }) ]# s  |0 {3 [. p+ Z>>A=magic(3)%生成一个3*3的数组
, \8 Q  P' G1 T6 ~  Z$ {2 {$ t
A=
0 X& h1 Y" d9 r/ e$ L7 X- X# U     8     1     6
3 [% l( ]# U4 D' [     3     5     7
     4     9     2

3 @4 _, `' R( |>>x=logical([1 1 0;0 1 1;1 0 1])%将double转化为boolean型数据

x =
     1     1     0
     0     1     1
     1     0     1

. {8 @8 F: ?2 W/ ?- ]# _& A>>A(x)%引用对应位置为1的数据，返回列向量
' A* I1 l% |, L- j
4 z; x0 B, ?/ |: d& L4 Mans =
3 l: G9 i% p( B- \1 x  j     8
$ [7 S* w" [* P& H8 _$ m     4
     1
* H5 W' p; {" d! r     5
     7
4 P0 X9 ?" ~8 `# f2 D% o5 J& l     2

* J2 p, A# m0 S  `; U6 k+ Q( c>>x=A>5%是有了比较语句，返回布尔型数据，对应位置数据大于5的为1，否则为0

8 v/ V8 P/ r$ o0 }4 A& Z) hx =
     1     0     1
     0     0     1
6 Y. m) g4 q0 C) p     0     1     0

>>A(x)%返回大于A中大于5的元素，其实该命令可以一次性执行A(A>5)或者find(A>5)，前者返回具体元素，后者返回大于5的数据的索引值

8 i! ?' y: J7 b4 u$ Q3 O% d2 g8 Zans =
     8
     9
     6
( Z, \1 g$ T, d8 V2 B     7

- I% X0 G( ?0 Q>>A(A>5)%一次性执行上面的命令
0 m- `; ?6 t& ?) p
4 t4 _1 ?" o. Uans =
     8
. u3 f/ v4 Z; [$ `3 `     9
* y! e9 R$ v; @+ M' l     6
; _- _! z7 `( @! Q# |1 a+ Q     7

) u2 J3 f5 J2 g& z( o>>indx=find(A>5)%查找A中对于5的元素，返回它们的索引(index)值，此时我们可以通过A(index)返回具体的元素
9 h' W- f( t1 m% @/ e
6 _( T5 H7 i5 J2 L  Nindex =
     1
     6
7 B9 @5 H" @/ r3 s     7
; p+ Z; a; V8 o( K  Q1 q     8

xiaogegepcb

Matlab中数组元素引用——三种方法