请选择 进入手机版 | 继续访问电脑版
楼主: 815102137
15181 22

[编程问题求助] stata中如何得到一个现有矩阵的特征值和特征向量 [推广有奖]

  • 0关注
  • 4粉丝

讲师

71%

还不是VIP/贵宾

-

威望
0
论坛币
772 个
通用积分
13.8625
学术水平
13 点
热心指数
15 点
信用等级
10 点
经验
5281 点
帖子
423
精华
0
在线时间
520 小时
注册时间
2011-8-31
最后登录
2024-1-25

815102137 发表于 2014-8-25 10:19:25 |显示全部楼层 |坛友微信交流群
50论坛币
如题,在以前帖子中看到1.可将Excel中数据粘到stata的data editor 中,再执行命令“mkmat var*,mat(X)”即可得到矩阵X

2.计算手动输入的矩阵特征值和特征向量的命令为:
   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
我的是已经粘贴好的矩阵X,就没有执行手工输入矩阵的命令,直接执行下边的命令,结果出错了。请问如何解呀!求大神解救。

最佳答案

pkl 查看完整内容

好像要先讲stata中的矩阵转换成mata中的矩阵,现在mata中设置一个只有缺失值的矩阵,然后用st_view把stata中X矩阵的数据转换到mata中你设置的新矩阵中,然后就可以用求特征根的命令了。 思路是这样的,不知对你有没有帮助,自己查查help文档。 mata V=st_view(X, ., .) p=. l=. eigensystem(V, p,l)
关键词:Stata 特征向量 tata 特征值 System 如何
pkl 在职认证  发表于 2014-8-25 10:19:26 |显示全部楼层 |坛友微信交流群
好像要先讲stata中的矩阵转换成mata中的矩阵,现在mata中设置一个只有缺失值的矩阵,然后用st_view把stata中X矩阵的数据转换到mata中你设置的新矩阵中,然后就可以用求特征根的命令了。
思路是这样的,不知对你有没有帮助,自己查查help文档。
mata
V=st_view(X, ., .)
p=.
l=.
eigensystem(V, p,l)
已有 1 人评分论坛币 学术水平 热心指数 收起 理由
crystal8832 + 10 + 1 + 1 热心帮助其他会员

总评分: 论坛币 + 10  学术水平 + 1  热心指数 + 1   查看全部评分

使用道具

815102137 发表于 2014-8-25 21:33:16 |显示全部楼层 |坛友微信交流群
pkl 发表于 2014-8-25 10:19
好像要先讲stata中的矩阵转换成mata中的矩阵,现在mata中设置一个只有缺失值的矩阵,然后用st_view把stata中 ...
非常感谢你,我好好学习一下。

使用道具

815102137 发表于 2014-8-25 21:45:47 |显示全部楼层 |坛友微信交流群
pkl 发表于 2014-8-25 10:19
好像要先讲stata中的矩阵转换成mata中的矩阵,现在mata中设置一个只有缺失值的矩阵,然后用st_view把stata中 ...
执行V=st_view(X, ., .)时错误提示是“X not found”。之前已经按步骤创建了X矩阵啊。看了help文件,仍然没弄好。能麻烦您再指点一下吗?多谢了。

使用道具

pkl 在职认证  发表于 2014-8-25 22:31:56 |显示全部楼层 |坛友微信交流群
815102137 发表于 2014-8-25 21:45
执行V=st_view(X, ., .)时错误提示是“X not found”。之前已经按步骤创建了X矩阵啊。看了help文件,仍然 ...
不好意思,刚才弄错了,应该用v=st_matrix("x")把stata中矩阵转换到mata中,所以改写为:
v=st_matrix("x")
p=.
l=.
eigensystem(v, p,l)
p
l
结果为
p
                      -.723241249    .129430781 + .330324585i    .129430781 - .330324585i
  2 |                -.580525181                  .792993021                  .792993021
  3 |                -.284953867   -.007865144 - .383946807i   -.007865144 + .383946807i
  4 |                -.220148611    -.29553534 + .019452345i    -.29553534 - .019452345i
  5 |                 -.10123484    .008825556 + .100104836i    .008825556 - .100104836i
。。。
l
                               1                           2                           3
    +-------------------------------------------------------------------------------------
  1 |                 5.70476568   -.068718011 - 1.70348216i   -.068718011 + 1.70348216i
    +-------------------------------------------------------------------------------------
                               4                           5
     ---------------------------------------------------------+
  1      -.28366483 - .67491643i     -.28366483 + .67491643i  |
     ---------------------------------------------------------+
已有 1 人评分学术水平 热心指数 收起 理由
crystal8832 + 1 + 1 补偿

总评分: 学术水平 + 1  热心指数 + 1   查看全部评分

使用道具

815102137 发表于 2014-8-26 09:16:33 |显示全部楼层 |坛友微信交流群
pkl 发表于 2014-8-25 22:31
不好意思,刚才弄错了,应该用v=st_matrix("x")把stata中矩阵转换到mata中,所以改写为:
v=st_matrix(" ...
成功了!非常感谢你帮助我,谢谢!你太赞了!

使用道具

815102137 发表于 2014-8-27 08:22:01 |显示全部楼层 |坛友微信交流群
pkl 发表于 2014-8-25 22:31
不好意思,刚才弄错了,应该用v=st_matrix("x")把stata中矩阵转换到mata中,所以改写为:
v=st_matrix(" ...
能再问您一个问题吗?得到的特征矩阵怎样才能使得在mata模式之外也可用。退出mata模式后,我执行mat list p,显示的是“p not found”。

使用道具

pkl 在职认证  发表于 2014-8-27 12:58:50 |显示全部楼层 |坛友微信交流群
815102137 发表于 2014-8-27 08:22
能再问您一个问题吗?得到的特征矩阵怎样才能使得在mata模式之外也可用。退出mata模式后,我执行mat list ...
Mata是stata的一个编译程序,其实是相当于两个编译器,所以,stata中的变量在Mata中不能直接引用,Mata中的变量在stata中也不能直接引用,所以要转换,你看看Mata手册吧。

使用道具

815102137 发表于 2014-8-27 15:37:34 |显示全部楼层 |坛友微信交流群
pkl 发表于 2014-8-27 12:58
Mata是stata的一个编译程序,其实是相当于两个编译器,所以,stata中的变量在Mata中不能直接引用,Mata中 ...
好的,谢谢。

使用道具

815102137 发表于 2014-8-27 15:48:49 |显示全部楼层 |坛友微信交流群
这是我找到的,供后来读贴的人们使用吧。

mata模式下矩阵的输出.docx

18.86 KB

已有 1 人评分论坛币 收起 理由
crystal8832 + 10 补偿

总评分: 论坛币 + 10   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 00:30