统计建模与R软件第二章作业

1、

x<-matrix(c(1,2,3),3,1)
y<-matrix(c(4,5,6),3,1)
e<-matrix(c(1,1,1),3,1)
z<-2*x+y+e
crossprod(x,y)###x,y的内积///////x%*%y
tcrossprod(x,y)###x,y的外积///////x%o%y

2、

A<-matrix(c(1:20),4,5,byrow=FALSE)
B<-matrix(c(1:20),4,5,byrow=TRUE)
c<-A+B
D<-tcrossprod(A,B)
E<-A*B
F<-A[c(-4,-5),c(-4,-5)]###保留A的前三行，前三列。
G<-B[,-3]###不包含B的第三列。

3、

X<-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2))
结果：[1] 1 1 1 1 1 2 2 2 3 3 3 3 4 4

4、生成一个5阶的HilbertJ矩阵n<-5

H<-matrix(0,5,5)
for(i in 1:n){
for(j in 1:n){
H[i,j] = 1/(i+j-1)
}
}
det(H)###求行列式
solve(H)###求矩阵逆
ev<-eigen(H)###求特征值和特征向量
错误形式：
H<-matrix(0,5,5)
i=1
j=1
if(j<=5){