24小时热门版块排行榜    

查看: 2690  |  回复: 17

houyanxu

金虫 (初入文坛)

[求助] 关于MATLAB,eig求特征值函数的问题 已有7人参与

矩阵
A =
     0    -2     2
    -1     0    -1
    -2    -2     0

用eig函数求特征值
[V,D]=eig(A)
结果为

V =

  -0.5773 - 0.0000i  -0.5773 + 0.0000i   0.5774         
   0.5774 + 0.0000i   0.5774 - 0.0000i  -0.5774         
   0.5774             0.5774            -0.5773         


D =

   1.0e-04 *

  -0.0540 + 0.0936i        0                  0         
        0            -0.0540 - 0.0936i        0         
        0                  0             0.1081         

但手算这个矩阵的特征值明显是 0, 0 , 0
请教一下各位~谢谢啦
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

因为太想学习,所以什么都学不好
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hlg154003231

金虫 (小有名气)

引用回帖:
11楼: Originally posted by houyanxu at 2016-01-06 17:03:09
什么意思,能具体说说么?这涉及计算特征值的算法么?...

就是A 原则上是不可能对角化的,,,但是MATLAB任然可以找到一个对角矩阵D和相应的矩阵V,,,使得A=V*D*V^-1(这里不可能严格相等),,,,,,,,这种分解在工程应用方面意义重大噢。。。。
13楼2016-01-06 19:19:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dk1013

木虫 (正式写手)

上面几位说的已经解释你的问题了。
首先你这个矩阵是奇异矩阵,条件数应该是Inf。上面给出的应该是数值近似。
但无论如何,这么大的条件数意味着你用matlab算出的结果基本上没什么用处。
你可以多看几位小数点(用format long命令)。
你会发现其实MATLAB给出的解(u,v)=eig(A)并不严格满足A*u=u*v(后面一些小数对不上)。换句话说就是解是错的。
你手算的没问题。

» 本帖已获得的红花(最新10朵)

到得还来别无事 庐山烟雨浙江潮
14楼2016-01-06 22:07:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jollage

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你手算怎么算的?一个矩阵的特征根都是零的唯一可能是零矩阵,但这个矩阵明显不是零矩阵。你可以很容易验算一下matlab算的对不对。

三个特征根的和是零,所以你看matlab求出来的特征根的和也是零。
2楼2016-01-05 22:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houyanxu

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by jollage at 2016-01-05 22:32:32
你手算怎么算的?一个矩阵的特征根都是零的唯一可能是零矩阵,但这个矩阵明显不是零矩阵。你可以很容易验算一下matlab算的对不对。

三个特征根的和是零,所以你看matlab求出来的特征根的和也是零。

你有手算么?

发自小木虫IOS客户端
因为太想学习,所以什么都学不好
3楼2016-01-05 23:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

getengqing

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
估计你算的是实数范围内的,matlab这里算的是复数范围内的···
一起交流学习/分享优秀资源
4楼2016-01-06 09:03:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
非对称矩阵特征根为负数。楼主确定没有弄错?

发自小木虫IOS客户端
数值计算
5楼2016-01-06 10:23:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houyanxu

金虫 (初入文坛)

引用回帖:
5楼: Originally posted by 独孤神宇 at 2016-01-06 10:23:51
非对称矩阵特征根为负数。楼主确定没有弄错?

关于MATLAB,eig求特征值函数的问题
求特征值步骤,请问我是哪里错了么?
因为太想学习,所以什么都学不好
6楼2016-01-06 13:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houyanxu

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by jollage at 2016-01-05 22:32:32
你手算怎么算的?一个矩阵的特征根都是零的唯一可能是零矩阵,但这个矩阵明显不是零矩阵。你可以很容易验算一下matlab算的对不对。

三个特征根的和是零,所以你看matlab求出来的特征根的和也是零。

关于MATLAB,eig求特征值函数的问题-1
请问我是哪里算错了么?
因为太想学习,所以什么都学不好
7楼2016-01-06 13:12:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxj126

木虫 (小有名气)


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 屏蔽内容, 欢迎常来 2016-01-07 09:39:35
本帖内容被屏蔽

8楼2016-01-06 13:55:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxj126

木虫 (小有名气)

【答案】应助回帖

引用回帖:
7楼: Originally posted by houyanxu at 2016-01-06 13:12:54
请问我是哪里算错了么?...

之前没看清。。回答错了。应该是matlab的算法不同吧
9楼2016-01-06 14:17:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlg154003231

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
A矩阵的特征值只有0,,,,但由于A矩阵无法对角化,所以MATLAB求解出来的是A矩阵的广义对角化分解;
即A*V=V*D;
10楼2016-01-06 16:23:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 houyanxu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 354求调剂 +5 Tyoumou 2026-03-18 8/400 2026-03-21 00:35 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
请填处理意见