矩阵开方比较复杂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里面的函数替代实现求出的矩阵开方是有一定误差的:
输出如下,可以看出b和原始的a是有一定误差的
七次四项式 则最高次数为7次,
与2xy的2n次方的次数相同
你对这个回答的评价是?
版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。
点击添加站长微信