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

[求助]求特征值的问题

<P><FONT color=#800080><FONT color=#000000>下面是我求特征值时的出现的警告：</FONT><br><br><U>*** FATAL    ERROR 1310 from s_error_post.EIG function has generated error<br>***          messages from lower-level routines:<br></U></FONT><br>为什么会出现不能求特征值的情况？<br>我觉得，至少也应该有虚数特征值啊。<br><br>先谢谢各位了！<br></P>
[此贴子已经被作者于2005-12-28 22:25:15编辑过]

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

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

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

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

<P>用fortran编的</P>

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

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

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

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

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

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

<P>的确挺奇怪的哈，<BR>用matlab算过了，是可行的。<BR>在看看。</P>
页: [1]
查看完整版本: [求助]求特征值的问题