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

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

(10页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

19.90 积分

实用文档数值分析之幂法及反幂法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)的值。③求A的(谱范数)条件数和行列式：在(1)中用反幂法求矩阵A的按模最小特征值时，要用到Doolittle分解方法，在Doolittle分解完成后得到的两个矩阵分别为L和U，则A的行列式可由U阵求出，即：det(A)=det(U)。 求得det(A)不为0，因此A为非奇异的实对称矩阵，则：，和分别为模最大特征值与模最小特征值。2、程序源代码：#include#include#include#define N 501 //列#define M 5 //行#define R 2 //下带宽#define S 2 //上带宽#define K 39#define e 1.0e-12 //误差限float A[M][N]; //初始矩阵float u[N]; //初始向量float y[N],yy[N];float maximum,value1,value2,value_1,value_N,value_s,value_abs_max;const float b=0.16f,c=-0.064f;int max_sign,max_position;void Init_matrix_A() //初始化矩阵A{ int i; for(i=2;i

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。