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

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

ID:45033266

大小:411.00 KB

页数:18页

时间:2019-11-08

数值分析报告之幂法及反幂法C语言程序实例_第1页
数值分析报告之幂法及反幂法C语言程序实例_第2页
数值分析报告之幂法及反幂法C语言程序实例_第3页
数值分析报告之幂法及反幂法C语言程序实例_第4页
数值分析报告之幂法及反幂法C语言程序实例_第5页
资源描述:

《数值分析报告之幂法及反幂法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)

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。