SAS矩阵运算

SAS矩阵运算 Proc iml; /*进入矩阵运算*/ X={1 2 3, 3 22 1, 3 1 2}; Print x; Y1=det(x); /*行列式*/ Print y1; Y2 =inv(X); /*矩阵求逆*/ Print y2; Y3=trace(X);/*迹*/ Print y3; Y4=diag(X);/*对角矩阵元素*/ Print y4; Y5=t(x); Y55=x`; Print y5 y55; Y6=x*t(x); Y66=t(x)*x; Print y6 y66;/*说明矩阵乘法不满足交换率*/ Xx={3 -1 2,4 0 5}; Yy={3 4, 6 -2, 4 3}; Xy=xx*yy; Print xy; A={3,-1,2}; La=sqrt(t(a)*a); Ea=a/la; Print a la ea; *解方程组AX=B,X=inv(A)*B a={321,-2 -14, 215}; B={10,3,20}; x=inv(a)*b; print x; *求矩阵的特征值与特征向量 A={13 -4 2,-4 13 -2, 2 -2 10}; E=eigval(a); V=eigvec(a); Print e ; Print v; Quit; /*退出矩阵运算*/ 单位矩阵 I3=i(3); Print i3; One矩阵(3行1列,全是1的元素的矩阵) J4=j(3,1,1); Print j4; MATLAB (X=[1 2 3;3 22 1;3 1 2]; a=[13 -42;-413 -2; 2 -210] [V D]=eig(a) I3=eye(3) J3=one(3,1,1)