数值分析作业第一题.docVIP

数值分析作业第一题 算法设计方案 利用带状Dollittle分解，将A[501][501]转存到数组C[5][501]，以节省存储空间 1、计算和 首先使用幂法求出矩阵的按模最大的特征值λ0：如果λ00，则其必为按模最大值，因此λ501=λ0，然后采用原点平移法，平移量为λ501，使用幂法迭代求出矩阵A-λ501I的按模最大的特征值，其特征值按从小到大排列应为λ1-λ501、λ2-λ501、……、。因此A-λ501I的按模最大的特征值应为λ1-λ501。又因为λ501的值已求得，由此可直接求出λ1。 2、计算，其中λ1和λn分别是矩阵A的模最大和最小特征值，直接利用上面求得的结果直接计算。 矩阵A的行列式可先对矩阵A进行LU分解后，det(A)等于U所有对角线上元素的乘积。 二、源程序： #includemath.h #includestdio.h #includestdlib.h #includeiostream.h #define s 2 #define r 2 int Max(int v1,int v2); int Min(int v1,int v2); int maxt(int v1,int v2,int v3); void storage(double C[5][501],double b,double c); double mifa(doubl