2015 2.13 matlab

                                                              数组运算

http://wenku.baidu.com/link?url=4WxXul51iVOqyNlyYapmMZmebo9fg-4J-GK_DfdSjw7JW1ELv17FvnmW9n5f7jib6n11AeJ9tQjYuoAZzYzeow2xGe0SZscY2Ua39FmHSv3

例程:>>A=magic(3)  A =

8     1     6      

3     5     7      

4     9     2

>>A(2:3,3:-1:1)

ans =

7     5     3      

2     9     4

>>A(:,end)

ans =      6      7      2

>>A(1,end-1)

ans =      1

>>A([2 1 3 3],[1 1 2 2 1])

ans =

3     3     5     5     3     

 8     8     1     1     8  

 4     4     9     9     4 

 4     4     9     9     4

A(ii,jj):其中ii和jj可以是一维向量、标量、“:”号或者“end”

大家对下标估计比较熟悉,由于在C语言中接触过,但是我这里需要强调的是,Matlab的下标是可以多行多列同时引用的,而像C语言等一次只能引用一个,比如

A(2:3,3:-1:1)表示引用数组中的2~3行,3~1列对应的元素

A(:,end)表示引用最后一列元素,“:”表示所有列或行,“end”表示最后一列或列,“end-n”表示倒数第n行或列

A(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列对应的元

1、向量的创建

1)直接输入:

行向量:a=[1,2,3,4,5]

列向量:a=[1;2;3;4;5]

2)用“:”生成向量

a=J:K生成的行向量是a=[J,J+1,…,K]

a=J:D:K生成行向量a=[J,J+D,…,J+m*D],m=fix((K-J)/D)

3)函数linspace用来生成数据按等差形式排列的行向量

x=linspace(X1,X2):在X1和X2间生成100个线性分布的数据,相邻的两个数据的差保持不变。构成等差数列。

x=linspace(X1,X2,n):在X1和X2间生成n个线性分布的数据,相邻的两个数据的差保持不变。构成等差数列。

4)函数logspace用来生成等比形式排列的行向量

X=logspace(x1,x2)在x1和x2之间生成50个对数等分数据的行向量。构成等比数列,数列的第一项x(1)=10x1,x(50)=10x2

X=logspace(x1,x2,n)在x1和x2之间生成n个对数等分数据的行向量。构成等比数列,数列的第一项x(1)=10x1,x(n)=10x2

注:向量的的转置:x=(0,5)’

2)函数eye,生成单位矩阵

eye(n) :生成n*n阶单位E

eye(m,n):生成m*n的矩阵E,对角线元素为1,其他为0

eye(size(A)):生成一个矩阵A大小相同的单位矩阵

eye(m,n,classname):对角线上生成的元素是1,数据类型用classname指定。其数据类型可以是:duoble、single、int8、uint8、int16、uint16、int32、uint32。

3)函数ones用ones生成全1的矩阵

ones(n) :生成n*n的全1矩阵

ones(m,n) :生成m*n的全1矩阵

ones(size(A)) :生成与矩阵A大小相同的全1矩阵

ones(m,n,p,…)生成m*n*p*….的全1的多维矩阵

ones(m,n,…,classname)制定数据类型为classname

4)函数zeros函数zeros生成全0矩阵

zeros(n):生成n*n的全0矩阵

zeros(m,n:)生成m*n的全0矩阵

zeros(size(A)):生成与矩阵A大小相同的全0矩阵

zeros (m,n,p,…)生成m*n*p*….的全0的多维矩阵

zeros (m,n,…,classname)指定数据类型为classname

5)函数rand函数rand用来生成[0,1]之间均匀分布的随机函数,其调用格式是:

Y=rand:生成一个随机数

Y=rand(n):生成n*n的随机矩阵

Y=rand(m,n):生成m*n的随机矩阵

Y=rand(size(A)):生成与矩阵A大小相同的随机矩阵

Y=rand(m,n,p,…):生成m*n*p*…的随机数多维数组

6)函数randn函数rand用来生成服从正态分布的随机函数,其调用格式是:

Y=randn:生成一个服从标准正态分布的随机数

Y=randn(n):生成n*n的服从标准正态分布的随机矩阵

Y=randn(m,n):生成m*n的服从标准正态分布的随机矩阵

Y=randn(size(A)):生成与矩阵A大小相同的服从标准正态分布的随机矩阵

Y=randn(m,n,p,…):生成m*n*p*…的服从标准正态分布的随机数多维数组

3、矩阵元素的提取与替换

1)单个元素的提取

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a(1,2)

b =

2

2)提取矩阵中某一行的元素,

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a(1,:)

b =

1     2     3

3)提取矩阵中某一列:

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a(:,1)

b =

1

3

4)提取矩阵中的多行元素

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a([1,2],:)

b =

1     2     3

3     4     5

5)提取矩阵中的多列元素

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a(:,[1,3])

b =

1     3

3     5

6)提取矩阵中多行多列交叉点上的元素

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入b=a([1,2],[1,3])

b =

1     3

3     5

7)单个元素的替换:

如:a=[1,2,3;3,4,5],运行后:

a =

1     2     3

3     4     5

输入:a(2,3)=-1

a =

1     2     3

3     4    -1

4、矩阵元素的重排和复制排列

1)矩阵元素的重排

B=reshape(A,m,n):返回的是一个m*n矩阵B,矩阵B的元素就是矩阵A的元素,若矩阵A的元素不是m*n个则提示错误。

B=reshape(A,m,n,p):返回的是一个多维的数组B,数组B中的元素个数和矩阵A中的元素个数相等

B=reshape(A,…,[],…):可以默认其中的一个维数

B=reshape(A,siz) :由向量siz指定数组B的维数,要求siz的各元素之积等于矩阵A的元素个数

2)矩阵的复制排列函数是repmat

B=repmat(A,n):返回B是一个n*n块大小的矩阵,每一块矩阵都是A

B=repmat(A,m,n):返回值是由m*n个块组成的大矩阵,每一个块都是矩阵A。

B=repmat(A,[m,n,p,…]):返回值B是一个多维数组形式的块,每一个块都是矩阵A

你可能感兴趣的:(2015 2.13 matlab)