Matlab、R向量与矩阵操作 z

描述

Matlab

R

1

建立行向量v=[1 2 3 4]

v=[1 2 3 4]

v<-c(1,2,3,4)或v<-scan()，然后输入1 2 3 4，并按Enter

2

建立列向量v=[1 2 3 4]’

v=[1;2;3;4]

同上，R中不区分行列

3

建立矩阵A=(142536)

A=[1 2 3;4 5 6]

A<-matrix(c(1,2,3,4,5,6),

nrow=2,byrow=TRUE)

4

访问向量中元素，如第3个元素

v(3)

v[3]

5

访问矩阵中元素，如第二行，第三列

A(2,3)

A[2,3]

6

用访问向量元素的方式访问矩阵元素

A(6) **先数列后数行

A[6] **先数列后数行

7

生成等差递增序列(单位步长)，如2,3,4,5,6,7

2:7

2:7

8

生成等差递减序列(单位步长)，如7,6,5,4,3,2

7:-1:2

7:2

9

生成等差递增序列，

如2,5,8,11,14

2:3:14

seq(2,14,3)

10

7-9归纳起来，介于a和b之间，等差的n个数，可用函数

linspace(a,b,n)

seq(a,b,len=n)

11

构建以10为底的对数等距，即在10^a和10^b之间的n个数

logspace(a,b,n)

10^seq(a,b,len=n)

12

构建所有元素都为0的k维向量

zeros(k,1)列向量

zeros(1,k)行向量

rep(0,k)

13

构建所有元素都为j的k维向量

j*ones(k,1) 列向量

j*ones(1,k) 列向量

rep(j,k)

14

构建所有元素都为0的n×m矩阵

zeros(n,m)

matrix(0, nrow=n,ncol=m)

或matrix(0 ,n, m)

15

构建所有元素都为j的n×m矩阵

j*ones(n,m)

matrix(j, nrow=n, ncol=m)

或matrix(j, n, m)

16

n×n的单位矩阵

eye(n)

diag(n)

17

构建对角阵，使其对角线上的元素为向量v

diag(v)

diag(v, nrow=length(v))

18

提取矩阵A对角线元素

v=diag(A)

v<-diag(A)

19

矩阵横向连接(行数要相等)

[A1 A2]

cbind(A1,A2)

20

矩阵纵向连接(列数要相等)

[A1;A2]

rbind(A1,A2)

21

给定向量x和y，长度分别为m和n，生成m×n的矩阵X，其行为x；生成m×n的矩阵Y，其列为y；

[X,Y]=meshgrid(x, y)

m<-length(x); n=length(y);

X<-matrix(rep(x,each=n),nrow=n);

Y<-matrix(rep(y,m),nrow=n)

22

向量元素掉转顺序

v(end:-1:1)

rev(v)

23

访问矩阵A的列，如第二列

A(:, 2)

A[,2]#给出的结果是向量;

A[,2,drop=FALSE]#结果是n×1的矩阵

24

访问矩阵A的行，如第七行

A(7, :)

A[7, ]#给出的结果是向量;

A[7, ,drop=FALSE]#结果是1×m的矩阵

25

将矩阵A“拉直”为向量

A(:)

c(A)

26

提取子矩阵，如2-4行，3-6列

A(2:4,3:6)

A[2:4,3:6]

27

提取特殊元素形成矩阵，如一个3×2矩阵，包含行7,7,6，列2,1的元素

A([7 7 6],[2 1])

A[c(7,7,6),c(2,1)]

28

按行掉转矩阵元素顺序

fliplr(A)

t(apply(A,1,rev))

29

按列掉转矩阵元素顺序

flipud(A)

apply(A,2,rev)

30

截取向量v，如保留前10个元素

v=v(1:10)

v<-v[1:10]或

length(v)<-10

31

截取向量v元素，从a处开始至结尾

v=v(a:end)

v<-v[a:length(v)]

32

截取向量v中除k以外的所有元素

v=v(1:(k-1) (K+1):end)

v<-v[-k]

33

截取除k和j以外的所有元素

32中的方法

v<-v[c(-k, -j)]

34

对矩阵A中的元素重新排成n行m列的新矩阵(A中必须有nm个元素)

A=reshape(A, n, m )

dim(A)=c(n, m)

35

提取矩阵A的下三角矩阵

L=tril(A)

L<-A; L[upper.tri(A)]<-0

36

提取矩阵A的上三角矩阵

U=triu(A)

U<-A; U[lower.tri(A)]<-0

37

提取n×n的Hilbert矩阵

hilb(n)

Hilbert(n)

**需加载Matrix库

38

构建n维数组，如3×4×2

reshape(1:24, 3, 4, 2)或

reshape(1:24, [3 4 2])

array(1:24, c(3,2,4))