矩阵开方比较复杂，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了，他要是早看到就以为没解决啊。。。。