数值分析幂法c语言实现(共5页).doc

《数值分析幂法c语言实现(共5页).doc》由会员分享，可在线阅读，更多相关《数值分析幂法c语言实现(共5页).doc（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上1 实验目的： 1熟练掌握C语言程序设计，编程求解问题。 2.运用幂法求解住特征值和特征向量。2 实验内容： 例题：用幂法求 A= 的特征值和特征向量。完整代码以及截图如下:#include stdio.h#include math.h#define M 3void main()float fan(),max(),e1,e2,r1,r2;void au(),ex(),print_x(),std();static float aMM=1.0,1.0,0.5,1.0,1.0,0.25,0.5,0.25,2.0;static float u0M,u1M,maxn0,maxn

2、1;int i;printf(*n);printf(* 幂法 *n);printf(*求特征值与特征向量*n);printf(*nn);printf(input precision e1,e2:);scanf(%f,%f,&e1,&e2);printf(ninput u(%d):,M);for (i=0;ie1 | r2e2)printf(%4d,i+);print_x(u0);printf(n);ex(u0,u1);elsebreak; while (1);void au(a,u0,u1)float aM,u0,u1;int i,j;for (i=0;iM;i+)u1i=0;for (j=0

3、;jM;j+)u1i+=aij*u0j;void std(u)float u;int i;float t,max();t=max(u);for (i=0;iM;i+) ui=ui/t;float fan(u0,u1)float u0,u1;float max();int i;float uuM;for (i=0;iM;i+)uui=u0i-u1i;return max(uu);float max(u)float u;int i;float m;m=u0;for (i=0;im)m=ui;return m;void ex(u0,u1)float u0,u1;int i;for (i=0;iM;i+)u0i=u1i;void print_x(u)float u;int i; for (i=0;iM;i+)printf(%12.6f,ui);3. 运行结果： 专心-专注-专业