qr方法求矩阵特征值c语言,QR方法求矩阵特征值

//QR方法求矩阵特征值

/******

数据

-1 2 1

2 4 -1

1 1 -6

******/

#include

#include

#include

using namespace std;

#include

#define N 3 //矩阵的维数

#define NUM 15 //QR分解次数

#define SNUM 13 //用来控制输出格式

void Print(long double A[N][N],long double Q[N][N],long double R[N][N]); //矩阵输出

void QR(long double A[N][N],long double Q[N][N],long double R[N][N]); //QR分解

void Multiplicate(long double A[N][N],long double R[N][N],long double Q[N][N]); //迭代，获得下一次矩阵A=QR

int main()

{

int i,j;

long double A[N][N];

long double Q[N][N]={0};

long double R[N][N]={0};

cout<

cout<

cout<

for(i=0;i

for(j=0;j

cin>>A[i][j];

cout<

cout<

cout<

for(i=1