数值分析报告之幂法及反幂法C语言程序实例

《数值分析报告之幂法及反幂法C语言程序实例》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用文档数值分析之幂法及反幂法C语言程序实例1、算法设计方案：①求、和的值：：表示矩阵的按模最小特征值，为求得直接对待求矩阵A应用反幂法即可。、：已知矩阵A的特征值满足关系，要求、及时，可按如下方法求解：a．对矩阵A用幂法，求得按模最大的特征值。b．按平移量对矩阵A进行原点平移得矩阵，对矩阵B用反幂法求得B的按模最小特征值。c．则：，即为所求。②求和A的与数最接近的特征值（k=0，1，…39）：求矩阵A的特征值中与最接近的特征值的大小，采用原点平移的方法：先求矩阵B=A-I对应的按模最小特征值，则+即为矩阵A与最接近的特征值。重复以上过程39次即可求得（k=0，1，…39）的

2、值。③求A的（谱范数）条件数和行列式：实用文档在（1）中用反幂法求矩阵A的按模最小特征值时，要用到Doolittle分解方法，在Doolittle分解完成后得到的两个矩阵分别为L和U，则A的行列式可由U阵求出，即：det(A)=det(U)。求得det(A)不为0，因此A为非奇异的实对称矩阵，则：，和分别为模最大特征值与模最小特征值。2、程序源代码：#include#include#include#defineN501//列#defineM5//行#defineR2//下带宽#defineS2//上带宽#defineK39#

3、definee1.0e-12//误差限floatA[M][N];//初始矩阵floatu[N];//初始向量floaty[N],yy[N];floatmaximum,value1,value2,value_1,value_N,value_s,value_abs_max;constfloatb=0.16f,c=-0.064f;intmax_sign,max_position;实用文档voidInit_matrix_A()//初始化矩阵A{inti;for(i=2;i

4、N;i++){A[2][i]=(1.64-0.024*(i+1))*sin(0.2*(i+1))-0.64*exp(0.1/(i+1));}for(i=0;i

5、])){max_position=i;maximum=fabs(u[i]);实用文档}}if(u[max_position]<0)max_sign=-1;elsemax_sign=1;}voidGet_y()//单位化迭代向量{inti;for(i=0;i

6、][N-4]*y[N-4]+A[3][N-3]*y[N-3]+A[2][N-2]*y[N-2]+A[1][N-1]*y[N-1];u[N-1]=A[4][N-3]*y[N-3]+A[3][N-2]*y[N-2]+A[2][N-1]*y[N-1];实用文档for(i=2;i

7、ax_position];}voidCheck_value()//幂法第二迭代格迭代{Init_u();Get_max();Get_y();Get_u();Get_value();while(1){Get_max();实用文档Get_y();Get_u();Get_value();if(fabs((value2-value1)/value1)