请教一个5行5列的矩阵从输入到计算其特征值的stata命令

aufeyuanph 发表于 2009-8-23 15:25
Sungmoo您好，我按照您的方式确实可以生成5*5的矩阵，但是stata 显示为：“
p=.
unrecognized command:  p
r(199);

. eigensystem(x,p,l)
unrecognized command:  eigensystem‘
请问在创造您演示的矩阵后，如何顺利解决求解该矩阵的特征根问题，烦劳了！

前面的命令是mata状态下的。

先要输入mata进入该状态。

后要输入end退出该状态。

aufeyuanph 发表于 2009-8-23 15:27 mat x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)或者以该矩阵为例

*注意一步也不能少。同时，eigensystem是mata命令。

mata
x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)   //手工输入5*5矩阵
p=.
l=.
eigensystem(x,p,l)
p   //显示特征向量
l    //显示特征值
end

aufeyuanph 发表于 2009-8-23 15:27 mat x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)或者以该矩阵为例，也可以

*使用arlion版主的命令的话，即

mat x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)  
mat symeigen hl=x
mat list l             /*特征根*/
mat list h             /*特征向量*/

很好的学习例子

你们写的我都看不懂呢，楼上哪位可以帮我计算下，                                       
        1        1        1        0.9427        1
        1        1        1        1        1
        1        1        1        1        1.6375
        1.0608        0.684        1        1        1.0068
        1        1        0.6107        0.9933        1
的最大特征值和特征向量，谢谢

mat x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)

mat eigenvalues H Lamda = x

mat list H

H[1,5]
              c1          c2          c3          c4          c5
real   5.7047657  -.28366483  -.28366483  -.06871801  -.06871801

mat list Lamda

Lamda[1,5]
                 c1          c2          c3          c4          c5
complex           0   .67491643  -.67491643   1.7034822  -1.7034822

晕！太高了吧

arlionn的方法很好的！

我通过计算了，分享给大家，需要调用mata命令，在命令框中输入mata后，会显示直线和冒号“：”，之后在命令框中直接输入冒号后面的算式即可，下面是显示结果，不要输入那个冒号“：”，希望你成功。

. mata
------------------------------------------------- mata (type end to exit) ---------------------------------------
: A=(2,-1,1\0,3,-1\2,1,3)
: X=.
: L=.
: eigensystem(A, X, L)
: X
                              1                           2                           3
    +-------------------------------------------------------------------------------------+
  1 |                -.577350269   -.577350269 - 3.9705e-23i   -.577350269 + 3.9705e-23i  |
  2 |                 .577350269    .577350269 - 1.9848e-08i    .577350269 + 1.9848e-08i  |
  3 |                -.577350269                  .577350269                  .577350269  |
    +-------------------------------------------------------------------------------------+
: L
                     1                 2                 3
    +-------------------------------------------------------+
  1 |                4   2 - 3.4378e-08i   2 + 3.4378e-08i  |
    +-------------------------------------------------------+
: end
----------------------------------------------------------------------------------------------------------------------

arlionn 发表于 2010-5-5 08:16
mat x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)

mat eigenvalues H  ...

请问，得到特征向量的命令是什么样的呢？这个是特征值吧？