Maier 发表于 2011-6-14 09:21

1. matlab 求的ew 有时候排列顺序不是很一样,查一下。
2. 如果模差得不是很大,完全可以忽略不计。
3. 试试命令polyeig

Maier 发表于 2011-6-14 09:29

非常重要的是,你要说清楚你的矩阵是干什么用的,这样我们才能知道解的大概形式。
从你叙述上看,应该是个很小误差。
把左右特征向量相应乘一乘,应该是对的

iceyoghurt 发表于 2011-6-14 11:27

回复 21# Maier
谢谢回复{:5_335:}
1,左右特征值不一样,一个是实数,一个是复数,我在这里附图了。
2,我比较关心实数还是复数的问题
3,试了polyeig,结果差不多

orionsnow 发表于 2011-6-14 23:39

我从图上看觉得你这个是软件的误差,你有没有尝试用decomposition 把矩阵分离出来看看?

明天要是有时间我用mathematica 算算看。

iceyoghurt 发表于 2011-6-15 07:26

回复 24# orionsnow

那就多谢了{:4_298:}

orionsnow 发表于 2011-6-28 14:47

回复orionsnow

那就多谢了
iceyoghurt 发表于 2011-6-15 08:26 http://www.dolc.de/forum/images/common/back.gif

我算过了, 就是简单的求特征值命令。

和你的结果差不多,实部从0-1。5, 虚部+- 0。2。

你这个矩阵的 det = -1.14203*10^-20 好小,回头我再用更多的参数,或者矩阵分解看看。

iceyoghurt 发表于 2011-6-28 14:51

回复 26# orionsnow

太感谢了,内牛满面{:4_300:}
最近读了一篇paper,好像是说此类矩阵对round off error很敏感。。。

orionsnow 发表于 2011-6-28 14:53

回复orionsnow

太感谢了,内牛满面
最近读了一篇paper,好像是说此类矩阵对round off erro ...
iceyoghurt 发表于 2011-6-28 15:51 http://www.dolc.de/forum/images/common/back.gif

gei wo kan kan

iceyoghurt 发表于 2011-6-28 14:56

回复 28# orionsnow

http://www.jstor.org/stable/2157635?seq=1

orionsnow 发表于 2011-6-28 19:08

回复orionsnow
iceyoghurt 发表于 2011-6-28 15:56 http://www.dolc.de/forum/images/common/back.gif

大概看了下摘要,就是你说的那个意思。

你看看方程里头有没有可以调节计算精度的地方,把精度调高些看看。
页: 1 2 [3]
查看完整版本: 请教一个矩阵的特征值问题