矩阵乘方和矩阵开方的命令

找了半天实在找不到矩阵乘方和矩阵开方的命令？找到的都是对数组的。矩阵乘方和矩阵开方的命令是？

正定阵与半正定阵的平方根，有确切的结论。一般的对称阵，也比较麻烦。

这里首先是一个数学问题。就一般的方阵及其平方根而言，它有没有平方根，有几个平方根，是一个比较麻烦的数学问题。

sungmoo 发表于 2012-9-9 18:30
这里首先是一个数学问题。就一般的方阵及其平方根而言，它有没有平方根，有几个平方根，是一个比较麻烦的数 ...

哦，这样子，其实我想做的是对一个对称的矩阵求开方

matrix a= (1,0.5,0.5,0.5, 0.5,0.5 \ 0.5,1,0.5,0.5, 0.5,0.5 \ 0.5,0.5,1,0.5, 0.5,0.5\0.5,0.5,0.5,1, 0.5,0.5\0.5,0.5,0.5, 0.5,1,0.5\0.5,0.5,0.5, 0.5,0.5,1)

mat list a
math b  = sqrt(a)   //a等于b*b(矩阵的乘法)
mat list b

问题是sqrt是对每个矩阵元素的开方，不是对矩阵的开方。

sungmoo，谢谢您的回答，有没有命令可以试试，先不管我这个特例有没有解？

针对stata中的矩阵  开根号的情形
help matpowersym  要在mata中运行 才行，并且不能出现虚数的情况，否则得不到结果。仅仅实验了对称矩阵
在stata运行如下代码
mata
a=(1,0,1\0,2,2\1,2,5)
b=matpowersym(a,0.5)
b
a
b*bend
如果直接把代码写入stata中，如果是循环语句碰到end会终止循环的，需要把这段代码单独保存do文件（my.do），在循环里面执行do文件:do my.do  即可，当然最好保存到stata的安装目录下

econfj 发表于 2012-9-9 19:45
哦，这样子，其实我想做的是对一个对称的矩阵求开方

matrix a= (1,0.5,0.5,0.5, 0.5,0.5 \ 0.5,1,0.5 ...

sqrt()也需要在mata下使用