请教一个矩阵的特征值问题

Maier

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

Maier

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

iceyoghurt

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

orionsnow

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

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

iceyoghurt

回复 24# orionsnow

那就多谢了

orionsnow

回复  orionsnow

那就多谢了
iceyoghurt 发表于 2011-6-15 08:26

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

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

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

iceyoghurt

回复 26# orionsnow

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

orionsnow

回复  orionsnow

太感谢了，内牛满面
最近读了一篇paper，好像是说此类矩阵对round off erro ...
iceyoghurt 发表于 2011-6-28 15:51

gei wo kan kan

iceyoghurt

回复 28# orionsnow

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

orionsnow

回复  orionsnow
iceyoghurt 发表于 2011-6-28 15:56

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

你看看方程里头有没有可以调节计算精度的地方，把精度调高些看看。