matlab reshape用法_MATLAB概率和统计

4.3 概率和统计

MATLAB不仅提供了强大的矩阵运算功能，在线性代数方面有广阔的应用，而且还能对大量的数据进行分析和统计，比如最常用的操作求平均值、最大值、标准差等，还有统计工具箱提供的二项分布、正态分布、泊松分布等，另外还有其他更为深入的统计功能。

4.3.1 基本分析函数

1．sum函数

sum函数用于求矩阵列方向元素或向量的和，调用格式如下。

(1)B=sum(A)：若A为向量，则返回所有元素的和；如A为矩阵，则分别对其各列所有元素求和并返回结果。

(2)B=sum(A,dim)：分别对矩阵A中第dim维的所有元素求和。

【例4-16】 sum函数应用示例。

>> A = pascal(6)

A =

1 1 1 1 1 1

1 2 3 4 5 6

1 3 6 10 15 21

1 4 10 20 35 56

1 5 15 35 70 126

1 6 21 56 126 252

>> A=A(:,1:4) % 创建演示矩阵

A =

1 1 1 1

1 2 3 4

1 3 6 10

1 4 10 20

1 5 15 35

1 6 21 56

>> B=sum(A) % 求各列的和

B =

6 21 56 126

>> C=sum(A') % 求转置后矩阵各列的和

C =

4 10 20 35 56 84

>> D=sum(A,1) % 求第1维方向也就是列方向各元素的和

D =

6 21 56 126

>> E=sum(A,2) % 求第2维方向也就是行方向各元素的和

E =

4

10

20

35

56

84

2．cumsum函数

cumsum函数用于求矩阵或向量的累积和，调用格式如下。

(1)B=cumsum(A)：若输入参数A为一个向量，则返回该向量所有元素的累积和；若A为矩阵，则返回该矩阵列方向各元素的累积和。

(2)B=cumsum(A,dim)：A为矩阵，dim为指定维数，若dim=1，则表示在列方向上求累积和；若dim=2，则表示在行方向上求累积和，如果A是多维矩阵那么dim可以为其他维数。

【例4-17】 cumsum函数应用示例。

>> cumsum(1:5)

ans =

1 3 6 10 15

>> A=magic(4)

A =

16 2 3 13

5 11 10 8

9 7 6 12

4 14 15 1

>> cumsum(A) % 列方向上求累积和

ans =

16 2 3 13

21 13 13 21

30 20 19 33

34 34 34 34

>> cumsum(A,1) % 列方向上求累积和

ans =

16 2 3 13

21 13 13 21

30 20 19 33

34 34 34 34

>> cumsum(A,2) % 行方向上求累积和

ans =

16 18 21 34

5 16 26 34

9 16 22 34

4 18 33 34

通过比较【例4-16】和【例4-17】可以看出cumsum函数和sum函数的区别：cumsum函数是累积和，结果中含有每一步计算的结果，而sum给出的则是最终求和的结果。

3．prod函数

prod函数用于求矩阵元素的积，其调用格式如下。

(1)B=prod(A)：若A为向量，则返回所有元素的积；若A为矩阵，则返回各列所有元素的积。

(2)B=prod(A,dim)：返回矩阵A中的第dim维方向的所有元素的积。

【例4-18】 prod函数应用示例。

>> prod(1:10) % 计算10的阶乘

ans =

3628800

>> M = magic(3)

M =

8 1 6

3 5 7

4 9 2

>> prod(M) % 列方向求积

ans =

96 45 84

>> prod(M,2) % 行方向求积

ans =

48

105

72

4．cumprod函数

cumprod函数用来求矩阵或向量的累积乘积，其调用格式如下。

(1)B=cumprod(A)：若输入参数A为一个向量，则返回该向量所有元素的累积乘积；若A为矩阵，则返回该矩阵列方向各元素的累积乘积。

(2)B=cumprod(A,dim)：A为矩阵，dim为指定维数，若dim=1，则表示在列方向上求累积乘积；若dim=2就代表在行方向上求累积乘积。

【例4-19】 cumprod函数应用示例。

>> cumprod(1:10)

ans =

Columns 1 through 7

1 2 6 24 120 720 5040

Columns 8 through 10

40320 362880 3628800>> A = [1 2 3; 4 5 6]

A =

1 2 3

4 5 6

>> cumprod(A,1) % 列方向上求累积乘积

ans =

1 2 3

4 10 18

>> cumprod(A,2) % 行方向上求累积乘积

ans =

1 2 6

4 20 120

通过比较【例4-18】和【例4-19】可以看出cumprod函数和prod函数的区别：cumprod函数是求累积乘积，结果中含有每一步计算的结果，而prod给出的则是最终乘积的结果。

5．sort函数

sort函数用于对矩阵元素按升序或者降序进行排序，其调用语法如下。

(1)B=sort(A)：对A进行默认的升序排序。输入参量A可以是向量、矩阵或字符串，若为向量，则对向量中的所有元素进行排序；若为矩阵，则对列方向各元素进行排序；若为字符串，则按其对应的ASCII码的大小进行排序。

(2)B=sort(A,dim)：对矩阵A中的第dim维进行升序排序。

(3)B = sort(...,mode)：按照指定升序或降序进行排序。mode可以是'ascend'(默认，升序)，或者是'descend'(降序)。

(4)[B,IX] = sort(A,...)：对A进行排序，并返回排序后各元素的下标值。

【例4-20】 sort函数使用示例。

>> A = [ 3 7 5;0 4 2 ]

A =

3 7 5

0 4 2

>> sort(A,1) % 列方向排序

ans =

0 4 2

3 7 5

>> sort(A,2) % 行方向排序

ans =

3 5 7

0 2 4

>> sort(A,2,'descend') % 行方向降序排序

ans =

7 5 3

4 2 0

>> [B,IX] = sort(A,2) % 排序并返回下标

B =

3 5 7

0 2 4

IX =

1 3 2

1 3 2

>> B=reshape(18:-1:1,3,3,2) % 创建高维矩阵

B(:,:,1) =

18 15 12

17 14 11

16 13 10

B(:,:,2) =

9 6 3

8 5 2

7 4 1

>> sort(B,1) % 列方向排序

ans(:,:,1) =

16 13 10

17 14 11

18 15 12

ans(:,:,2) =

7 4 1

8 5 2

9 6 3

>> sort(B,2) % 行方向排序

ans(:,:,1) =

12 15 18

11 14 17

10 13 16

ans(:,:,2) =

3 6 9

2 5 8

1 4 7

>> sort(B,3) % 页方向排序

ans(:,:,1) =

9 6 3

8 5 2

7 4 1

ans(:,:,2) =

18 15 12

17 14 11

16 13 10

6．sortrows函数

sortrows函数用于在保持各行相对元素不变的情况下，对各行整体进行升序排列。sortrows函数调用语法如下。

(1)B = sortrows(A)：按行对A进行升序排列。输入变量A必须是矩阵或者列向量。

(2)B = sortrows(A,column)：基于向量column指定的列对矩阵A进行排序。

(3)[B,index] = sortrows(A,...)：在对矩阵A进行排序的同时，返回下标索引。如果A是一个列向量，则B = A(index)；如果A是一个m×n的矩阵，则B = A(index,:)。

【例4-21】 sortrows函数使用示例。

>> rand('state',0) % 设定随机数种子，以便于读者验证

>> A = floor(rand(6,7) * 100); % 创建测试矩阵，floor函数用于取整，以便于观察

>> A(1:4,1)=95; A(5:6,1)=76;

>> A(2:4,2)=7; A(3,3)=73 % 修改部分数据，以体现函数用法

A =

95 45 92 41 13 1 84

95 7 73 89 20 74 52

95 7 73 5 19 44 20

95 7 40 35 60 93 67

76 61 93 81 27 46 83

76 79 91 0 19 41 1

>> B = sortrows(A) % 按照第1列元素大小对矩阵A进行排序

B =

76 61 93 81 27 46 83

76 79 91 0 19 41 1

95 7 40 35 60 93 67

95 7 73 5 19 44 20

95 7 73 89 20 74 52

95 45 92 41 13 1 84

通过比较可以看到，矩阵A中第1列具有相等的元素。sortrows函数在进行排序操作时，如果指定列中存在相等元素，则通过比较指定列右侧列中的元素来进行排序，右侧列中若还有相等元素，则按照右侧再下一列的元素进行排序。

>> C = sortrows(A,2) % 按照第2列的大小进行排序

C =

95 7 73 89 20 74 52

95 7 73 5 19 44 20

95 7 40 35 60 93 67

95 45 92 41 13 1 84

76 61 93 81 27 46 83

76 79 91 0 19 41 1

>> D = sortrows(A,[1 7]) % 按照第1列和第7列进行排序

D =

76 79 91 0 19 41 1

76 61 93 81 27 46 83

95 7 73 5 19 44 20

95 7 73 89 20 74 52

95 7 40 35 60 93 67

95 45 92 41 13 1 84

>> E = sortrows(A,[1 -7]) % 按照第1列和第7列进行排序

E =

76 61 93 81 27 46 83

76 79 91 0 19 41 1

95 45 92 41 13 1 84

95 7 40 35 60 93 67

95 7 73 89 20 74 52

95 7 73 5 19 44 20

矩阵D和E都是按照第1列和第7列进行排序，即首先按照第1列进行排序，如果第1列中存在相等元素，则按照第7列进行排序。参数[1 -7]中的负号的含义是按照降序排序。

>> F = sortrows(A, -4) % 按照第4列进行降序排序

F =

95 7 73 89 20 74 52

76 61 93 81 27 46 83

95 45 92 41 13 1 84

95 7 40 35 60 93 67

95 7 73 5 19 44 20

76 79 91 0 19 41 1