#include<iostream>
using namespace std;
int main()
{
double A[10]={1.0,1.0,0.5,1.0,1.0,0.25,0.5,0.25,2.0};//矩阵
double v[3];
double vk=-1;
double uk[3];
v[0]=1;
v[1]=1;
v[2]=1;
uk[0]=uk[1]=uk[2]=1;
int NUM=0; //迭代次数
while(NUM<10)
{
for(int i=0;i<3;i++)
{
double r=0;
for(int j=i*3;j<i*3+3;j++)
{
r+=A[j]*uk[j%3];
}
v[i]=r;
}
vk=0;
for(int i=0;i<3;i++)
{
if(vk<v[i])
{
vk=v[i];
}
}
for(int i=0;i<3;i++)
{
uk[i]=v[i]/vk;
cout<<uk[i]<<" ";
}
cout<<endl;
NUM++;
}
return 0;
}
幂法求矩阵的主特征值和特征向量
最新推荐文章于 2023-10-21 09:14:06 发布