北航数值分析大作业1(学硕).docxVIP

北京航空航天大学数值分析（1）学号： SY1307145姓名： 龙 安 林学院：机械工程及自动化学院 2013 年 10 月 20 日算法设计方案开始带状存储矩阵A，并输入允许误差ε求λ1与λ501用带平移量的幂法子函数求按模最大特征值（此时平移量为零），设为λ1（即暂存于λ1中）；以λ1为平移量，用带平移量的幂法子函数求特征值序列中另一端特征值，设为λ501（即暂存于λ501中）；如果特征值λ1大于零且特征值λ501小于零，或者，λ1与λ501均大于零，则交换λ1与λ501的值；否则，不更改λ1与λ501的值；完成λ1与λ501的求解；求λs如果，λ1与λ501均大于零，则λs=λ1；如果，λ1与λ501均小于零，则λs=λ501；否则，调用带平移量的反幂法子函数（此时平移量为零），求按模最小特征值λs；求λik令平移量；循环调用带平移量的反幂法子函数（先平移，后LU分解），求λik；求A的（谱范数）条件数cond（A）2设按模最大特征值为；条件数cond（A）2=； 求行列式detA调用三角分解函数，对A进行LU分解；；；结束程序源代码#include stdafx.h#includestdio.h#includemath.h//声明用具求解相应特征值的调用函数double p_method(double pm_mat[][501],double sub_epsilo