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

uqHZau

1.Matlab中数组元素引用有三种方法

2 F0 j4 c4 V" K2 @0 O1.下标法(subscripts)
2.索引法(index)
2 B$ A& Y- d: p, E/ H3.布尔法(Boolean)
: j0 ]8 k$ D8 s
注意：在使用这三种方法之前，大家头脑一定要清晰的记住，Matlab中数组元素是按列存储(与Fortran一样)，比如说下面的二维数组
  T- G! P9 ?2 _  k- F$ P! ?
A=
     8     1     6
8 A6 e- H# `% o( ~9 I- K     3     5     7
     4     9     2
Matlab的存储顺序是8,3,4,1,5,9,6,7,2，也就是说先行后列，对于3维数组呢，就是先行后列再页对应个元素的索引和下标分别为

$ l3 }; x4 C& N$ E  L) Y, f5 G+ B  lElement Index Subscripts
4 H- v# L' U8 d8        1      (1,1)
3        2      (2,1)
# Q/ y: g  i4 `3 |9 l4        3      (3,1)
  B8 X1 p5 i# x: b% C1        4      (1,2)
# h- W0 n4 D" Y5        5      (2,2)
9        6      (3,2)
6        7      (1,3)
7        8      (2,3)
. C0 D5 ^6 p( I: w$ b2        9      (3,3)
; o0 J4 k  D6 }& i2 x# I
+ D& I, h; ?; _, V. H　　从上面的例子中已经很清晰的说明了下标和索引的区别了，也就是说Matlab为每个元素分配了一个唯一识别的ID(即index)
) m. O$ Z7 j( l  M/ V2 d- w

2.下标法引用

, G; I% q; x+ S7 [: NA(ii,jj)：其中ii和jj可以是一维向量、标量、“:”号或者“end”

* L$ s2 U$ t5 C7 h/ T' g大家对下标估计比较熟悉，由于在C语言中接触过，但是我这里需要强调的是，Matlab的下标是可以多行多列同时引用的，而像C语言等一次只能引用一个，比如

A(2:3,3:-1:1)表示引用数组中的2~3行，3~1列对应的元素
9 h' l7 r2 P8 K7 dA(:,end)表示引用最后一列元素，“:”表示所有列或行，“end”表示最后一列或列，“end-n”表示倒数第n行或列
7 ^0 U& a4 H* j. I/ hA(1,end-1)表示引用第1行倒数第2个元素
; t5 @3 H; H9 HA([2 1 3 3],[1 1 2 2 1])表示引用按两个向量引用指定的元素，即A中的第2,1,3,3行和第1,1,2,2,1列对应的元素
! I- q: {% E! M8 ~% l; Y
% z7 n) _5 E, F7 ?
举几个例子：
- @& w' V$ k4 }* l. l) {) h# O+ W: L
/ v0 ^0 ?& S1 h- J. t6 S>>A=magic(3)

A =
$ d9 N( i6 l6 |$ ^- N: l3 J     8     1     6
+ f, m. V5 b2 q7 I6 T# O% N9 L# F     3     5     7
% T& E" x  i# }" r' z     4     9     2

>>A(2:3,3:-1:1)
) f0 x: H  M! M8 M+ k. |
ans =
; B. |: Z( _0 W+ m, H     7     5     3
     2     9     4

. j, O, R! t" ]' Z>>A(:,end)
, w+ a; q  F0 O. A2 Q
ans =
     6
: X0 i7 M" T. Q5 K& }5 o7 ~" |     7
     2
: f/ Z0 d2 m( M$ r# @
>>A(1,end-1)
3 _; ?' v; V7 A" \
3 k2 I; u) G- uans =
) t+ [8 L8 i) [0 m- G     1

: ]6 I$ t* K1 Y>>A([2 1 3 3],[1 1 2 2 1])

. K% I1 F8 v! n7 P- k% _ans =
     3     3     5     5     3
2 R% i+ o/ q" Q! @9 _     8     8     1     1     8
     4     4     9     9     4
' g; x) q  b3 {& ^     4     4     9     9     4
3 d! u) {8 X$ j2 [
  N" k+ ~$ H6 y4 r/ `
3 L$ \1 _& a5 M/ _) j2.索引法引用(说白了索引就是存储顺序)

+ ~7 i/ Q" s: M( k. x. S! vA(index)：index可以是任意的数组，index的元素必须是正整数，且不大于numel(A)，返回的是一个尺寸与index一样的数组

下标和索引之间可以通过ind2sub和sub2ind函数相互转换，具体可以看帮助，很简单

[I,J] = ind2sub(siz,IND)
IND = sub2ind(siz,I,J)
还有使用A(：）就可以将数组A转换为列向量
: c. n; h. v7 C6 [9 h3 ~. d3 N
, |+ I) e" X# O; PA(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)
A([2 5 9;1 1 1;8 5 6])：返回的时侯是一个3*3的矩阵

3 Y2 {9 V3 l) i) f4 v6 h5 u>>A=magic(5)%括号中为索引值

A =
/ S2 l2 P, N: a; E7 V+ A; x    17 (1)     24 (6)     1 (11)     8 (16)    15 (21)
    23 (2)      5 (7)     7 (12)    14 (17)    16 (22)
8 P+ T: l7 X; _     4 (3)      6 (8)    13 (13)    20 (18)    22 (23)
    10 (4)     12 (9)    19 (14)    21 (19)     3 (24)
4 }) s8 Z) B+ I, M" ?$ P, u    11 (5)    18 (10)    25 (15)     2 (20)     9 (25)
7 m* c) o0 B& Q7 Q& }- b
9 H" [: p' r  Q% G6 z9 U3 C>>A(8)
$ ^$ k) L# ^3 n& [- Y' M# j+ \
( i; k* l, v$ G- p( P0 P7 [ans =
9 g/ E. c  }* T( ?3 W     6
2 I% h6 A5 L1 [% U  i( _3 R
>>A([1 10 5 2 2 1 3])
5 b2 G6 g: B9 \& e
ans =
' }. k. K9 p7 R1 P    17    18    11    23    23    17     4

: t: Z3 E) {1 ?7 m" g: l>>A([2 5 9;1 1 1;8 5 6])
& {1 }. m& `3 C$ f! t
9 W5 F# J  a4 g' z# e. oans =
1 N: v' @4 n& W" b    23    11    12
    17    17    17
& u6 P: q2 }4 X4 f     6    11    24
9 x# _8 |* L4 Y　　
& J2 S2 [8 y5 H$ Y3 n) v" `

2 u3 @, C9 M! J  k& S5 F! ^3.布尔法引用
" J. ~" f6 l5 f, Z& c; ~! M+ G
A(X)：X是一个有0和1组成布尔型数据，且size(A)=size(X)，对应位置为1则留下该数据，0则去掉，最后按A中的存储顺序，返回一个列向量
( Y4 u% V) \* m) {/ _4 s/ a
假如说A是3*3的数组

5 o% c6 b2 j5 a. e$ V( m# WA(logical([1 0 0;0 1 0;0 0 1]))：表示引用了数组A的对角线元素，注意必须使用logical将0/1数组转换为布尔型
6 _2 T# B3 u! Q+ ]3 @4 d
4 p7 s, m! ~  |3 n: h% u! A>>A=magic(3)%生成一个3*3的数组
7 Z1 \( E/ v8 K1 T  l6 K1 @
A=
     8     1     6
     3     5     7
, S/ t% O6 z- I6 K' Q     4     9     2

- f# W' \0 q8 E) _( T>>x=logical([1 1 0;0 1 1;1 0 1])%将double转化为boolean型数据

x =
8 R/ L0 |$ _- c     1     1     0
, K' n' @6 ]/ p4 t3 v/ A     0     1     1
7 t, V* S$ j% C# b     1     0     1
, y8 p+ }$ X' Y$ u
7 y' T! U3 S, p& n>>A(x)%引用对应位置为1的数据，返回列向量
. R+ ?' }0 Q" A/ g" o
ans =
     8
- L4 t7 ?4 @8 X' l. z# q     4
5 J" u( |( v: u6 ], @( Q     1
( R8 A' E2 R  {$ m$ J/ x     5
     7
     2
$ U; `0 |9 _9 T# Y( s6 e
>>x=A>5%是有了比较语句，返回布尔型数据，对应位置数据大于5的为1，否则为0
" \/ y2 b% \3 g3 t) ]" t1 f- T
x =
     1     0     1
4 `. u; r+ W# q) D0 Y     0     0     1
     0     1     0
/ K- w7 o! o# X! u
>>A(x)%返回大于A中大于5的元素，其实该命令可以一次性执行A(A>5)或者find(A>5)，前者返回具体元素，后者返回大于5的数据的索引值

ans =
     8
2 A% Y( W: w/ ?) l     9
+ O, M' ~! {: J( `- Q     6
     7
1 A8 B) y9 X7 ~2 ~
>>A(A>5)%一次性执行上面的命令
* U) N" z2 j+ M5 r( [! J. D
ans =
; S& `5 A- G6 t% s5 X     8
' U( C& m: [* M. [7 B5 y% D     9
     6
     7

>>indx=find(A>5)%查找A中对于5的元素，返回它们的索引(index)值，此时我们可以通过A(index)返回具体的元素

index =
     1
# O% l0 J+ V: x$ w! ~$ Y$ t     6
     7
. }  K* u: G2 N, m; l4 j' j     8
# P$ n, |! F( a/ R) L

xiaogegepcb

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