弯弓射大雕 发表于 2008-3-16 15:53

广义特征值问题，请教

小弟用Fortran自带的imsl库中的GVCRG（N, A, n, B, n, ALPHA, bBETA, EVEC, N），求解广义特征值c，Ax=c*bx，因为矩阵a和b为三对角矩阵，大小为200*200，请教高手，a和b为稀疏矩阵吗，在计算特征值之前还要对a和b进行哪些处理，还是可以直接带到GVCRG中计算就可以，谢谢了。

欧阳中华 发表于 2008-3-16 17:22

.
    估计是矩阵太大，程序算不了吧，按一个字节8个BT，这个问题仅M阵和K阵就需要625K的内存...

弯弓射大雕 发表于 2008-3-16 17:50

谢谢，那有什么解决办法吗？

欧阳中华 发表于 2008-3-17 08:04

.
    在计算规模比较小的时候可以采用编辑器自带的特征值求解库函数，但由于库函数都是满阵运算，所以不能求解规模大的问题；当要求解大规模问题时，往往就需要有针对性的开发程序，还要采用一维压缩技术，技术难度要大很多。。。

弯弓射大雕 发表于 2008-3-17 09:22

我计算的矩阵是200*200的方阵的呀。

欧阳中华 发表于 2008-3-17 10:14

.
    关于结构有限元程序技术，有专门的书籍介绍，比如Bath等的著作，如果想学习的话，可以系统的看一下，不是一句二句就说清楚的... ..

弯弓射大雕 发表于 2008-3-17 13:15

谢了

VibrationMaster 发表于 2008-3-17 20:06

单次的200*200的方阵不算大，直接计算没问题。

弯弓射大雕 发表于 2008-3-18 10:03

楼主说的对，我程序中矩阵的定义应该为双精度类型，以前是单精度，所以出现了些问题，所以对矩阵计算我觉得用双精度好一些，与大家分享。
页: [1]
查看完整版本: 广义特征值问题，请教