1、实验2 R基础(二)1. (习题2.1)建立一个R文件,在文件中输入变量x = (1,2,3)T, y = (4,5,6)T,并作以下运算(1) 计算z = 2x + y + e,其中 e = (1,1,1)T;(2) 计算x与y的内积;(3) 计算x与y的外积。解:源代码:(1)x-c(1,2,3)y-c(4,5,6)e-c(1,1,1)z=2*x+y+ez1=crossprod(x,y) #z1为x与y的内积或者x%*%yz2=tcrossprod(x,y) #z2为x与y的外积或者x%o%yz;z1;z22. (习题2.2)将1,2,20构成两个45阶的矩阵,其中矩阵A是按列输入,矩阵B
2、 是按行输入,并作如下运算(1) C = A + B;(相对应的数相加)(2) D = ABT;(3) E = (eij )45,其中 eij = aij bij;(相对应的数相乘)(4) F是由A的前3行和前3列构成的矩阵;(5) G是由矩阵B的各列构成的矩阵,但不含B的第3列。解:源代码:A-matrix(1:20,c(4,5);A # A-matrix(1:20,nrow=4,byrow=F);A/ A-matrix(1:20,nrow=4);A # A-matrix(1:20,ncol=5);A # A-matrix(1:20,ncol=5,byrow=F);AB-matrix(1:20,nrow=4,byrow=TRUE);BC=A+B;CD-A%*%t(B);DE=A*B;EF-A1:3,1:3;FG-B,-3;G3. (习题2.3)构造一个向量x,向量是由5个1,3个2,4个3和2个4构成,注意用到 rep()函数。解:源代码:x-c(rep(1,times=5),rep(2,times=3),rep(3,times=4),rep(4,times=2); xx-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2);x运行截图: