随着高新技术在航天领域的不断应用和发展,航天器的尺寸越来越大,结构越来越复杂。这种复杂的结构极易受到外界及结构本身扰动的影响而发生振动,强烈的振动将影响系统有效载荷的正常工作,导致系统性能下降甚至破坏。所以许多结构都需要采用附加约束阻尼层等措施来提高其结构阻尼。分析阻尼结构振动最主要的手段是复模态分析方法,即通过求解复特征值问题来确定系统的固有频率和损耗因子。本论文重点研究含阻尼结构的复模态分析算法,并编制高效的计算机程序来实现算法,为开发有自主知识产权的有限元求解器奠定基础。本文首先对子空间迭代法、乘幂法和逆幂法、Lanczos方法和QR方法等几种特征值的求解算法进行了深入分析和比较,给出了各种算法的优缺点及其适用范围。并据此提出了本论文重点研究的Lanczos-QR迭代方法。Lanczos-QR迭代方法结合了Lanczos算法快速高效的优点和QR方法稳定可靠的优点。针对大规模矩阵问题,采用了改进的截断Lanczos迭代方法对矩阵进行降阶,得到一个包含了原矩阵低阶特征值的三对角矩阵。针对实矩阵含有复特征值的问题,采用了双重步位移QR方法求解降阶后的实三对角矩阵的全部特征值和特征向量。本文利用FORTRAN语言实现了这一针对大型复杂结构模态分析的计算方法。并通过一个标准约束阻尼梁和一个复杂的卫星结构对算法和程序进行了验证,通过与实验结果和航天上应用广泛的NASTRAN软件分析的结果的比较表明,本算法及其程序具有较高的精度和较快的收敛速度。