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;
Pri...
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)
本文档为【SAS矩阵运算】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。