用SAS中求一个矩阵的1/2次方，求解答~~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

假如u，v都是矩阵，矩阵u的平方等于矩阵v，已知v如何求u啊？？？我在iml的帮助里面找了很久都没找到，求SAS高人解答~~~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：求解答 高人解答 如何

没人知道吗？顶上去~~

我也不知，赐教

顶上~~~~~！！！！！

proc iml;
a={1 3,5 7};
b=sqrt(a);
print b;
quit;

这个是矩阵的每个元素开根号，不满足题目要求

矩阵开方比较复杂，sas里面没有现成的函数，理论上讲一般不是单值的，而且大多只对Hermite半正定矩阵定义算术根，即若A是Hermite半正定阵，其谱分解为A=Q*D*Q'，其中Q是酉阵，D是对角阵，Q'是Q的转置共轭，那么A^(1/k)定义成Q*D^(1/k)*Q'，D^(1/k)就是对D的每个对角元(是非负实数)取k次算术根，因此，可用sas里面的函数替代实现，求出的矩阵开方是有一定误差的：
proc iml;
a={7 10,15 24};
call svd(ua,qa,va,a);
sqrtm_a=ua*diag(sqrt(qa))*va`;
print sqrtm_a;
b=sqrtm_a*sqrtm_a;
print b;
quit;

输出如下，可以看出b和原始的a是有一定误差的。

                                                             sqrtm_a
                                                       1.7715868 1.4764294
                                                       2.4820666 4.4633641

                                                                b
                                                       6.8031158 9.2054646
                                                       15.475563 23.586215

版主，如果对题，麻烦赏点￥￥，因为我已经把链接给admin了，他要是早看到就以为没解决啊。。。。