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
大概看了下摘要,就是你说的那个意思。
你看看方程里头有没有可以调节计算精度的地方,把精度调高些看看。