北航数值分析大作业(一).docVIP

北京航空航天大学 《数值分析A》课程计算实习 《数值分析A》 计算实习题目一 姓名： 学号： 学院： 2014年11月 算法设计方案： A矩阵的存储与检索： 根据题目的设定，A矩阵为上半带宽与下半带宽均为2的带状矩阵，且除主对角线元素为变量ai外，其余元素均为常数b或c，故可以将A矩阵转存为新矩阵A[5][501]，在存入时原A矩阵的主对角元素存入新A矩阵的第三行，且各元素的列号保持不变。 求解λ1、λ501、λs逻辑关系图： 求解A与数μk=λ 可对A矩阵进行平移变换，对矩阵B=A+μkΙ采用反幂法求解按模最小的特征值βk，再经过反平移得到βk 求矩阵A的（谱范数）条件数cond(A)2和行列式 在采用反幂法求解λs的过程中要对矩阵A进行Doolittle分解（LU分解），det?(A)为U矩阵对角线元素的乘积。而根据公式cond(A)2=λmaxλmin可以得到A的条件数，其中λ 程序源代码 #includestdio.h #includemath.h #includestdlib.h #define N 501 #define M 5 #define e 1.0e-12 double a[M][N]={0.0},u[N]={0.0}; void read_A()/*设定A矩阵函数*/ { double b=0.16,c=-0.064; int i; dou