数值分析幂法求矩阵特征值
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
#include
#include
#define N 3
void matrixx(double A[N][N],double x[N],double v[N])
{
for(int i=0;i
v[i]=0;
for(int j=0;j
}
} double slove(double v[N])
{
double max;
for(int i=0;i
return max; }
void main()
{ //data input double
double A[N][N]={-1.0,1.0,0,-4.0,3.0,0,1.0,0,2.0};
double x[N]={1,1,1};
double v[N]={0,0,0};
double u[N]={0,0,0};
double p[N]={0,0,0};
double e=1e-10,delta=1;
int k=0;
while(delta >= e)
{
for(int q=0;q
matrixx(A,x,v);
for(int i=0;i
delta=fabs(slove(v)-slove(p));
k++;
for(int l=0;l
cout << "迭代次数 " << k << endl;
cout << "矩阵的特征值 " << slove(v) << endl;
cout << "( " ;
for(int i=0;i
cout << ") " <<
endl;
}