[求助]求特征值的问题

huhust · 发表于 2005-12-25 21:58

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

下面是我求特征值时的出现的警告： *** FATAL ERROR 1310 from s_error_post. EIG function has generated error *** messages from lower-level routines: 为什么会出现不能求特征值的情况？ 我觉得，至少也应该有虚数特征值啊。 先谢谢各位了！

[此贴子已经被作者于2005-12-28 22:25:15编辑过]

Silence · 发表于 2005-12-26 08:58

这是用什么工具求特征值呢？？

huhust · 发表于 2005-12-26 10:18

用fortran编的

linqus · 发表于 2005-12-26 15:52

看起来像是直接调用imsl里面的函数了， 可能是你的矩阵不太符合其函数的格式要求所致， 比如数据格式等等方面的。

huhust · 发表于 2005-12-27 09:56

您帮我看看下面的程序，sh用0.01或者0.001，都可以得出结果，但0.0001时却不行。为什么？ use imslf90 use linear_operators implicit none integer, parameter :: n=7 double PRECISION sh double COMPLEX A(n,n) double COMPLEX dd(n)
 sh = 0.0001 A = 0.d0 A(1,2) = 1.d0 A(2,3) = 1.d0 A(3,2) = -1.d0; A(3,4) = -1.d0; A(3,6) = -1.d0 A(4,5) = 1.d0 A(5,6) = 1.d0 A(6,7) = 1.d0 A(7,2) = -301.d0+sh; A(7,4) = -301.d0*(1.0-sh); A(7,6) = -1.d0
! Compute the eigenvalues and eigenvectors. dd = EIG(A)
 write (*,*) dd
 end

linqus · 发表于 2005-12-30 14:16

的确挺奇怪的哈， 用matlab算过了，是可行的。 在看看。

账号		自动登录	找回密码
密码			我要加入

[经典算法] [求助]求特征值的问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（Silence）这是用什么工具求特征值呢？？

回复：（linqus）看起来像是直接调用imsl里面的函数了...