chap9.cpp

jacobi 和QR 算法 用来求解特征值问题 自己编的 不知道是否好用

#include <iostream.h>
#include <math.h>
#include <iomanip.h>
#include <stdio.h>
#include <conio.h>
#include"chap9_Jacobi.h "
#include"chap9_QR.h "


void main()
{
//=====================================第一题===================================================

	cout<<"*****************************************************************************"<<endl;
	cout<<"                         程序说明（一）："<<endl<<endl;
	cout<<"第一题:编制用Jacobi方法计算矩阵的特征值和特征向量的程序"<<endl
		<<"   其中  A=[ 10  7   8   7"<<endl
		<<"             7   5   6   5"<<endl
		<<"             8   6   10  9"<<endl
		<<"             7   5   9   10]"<<endl;
    cout<<"*****************************************************************************"<<endl;
    cout<<endl;

	//建立eigenJacobi的一个对象m，矩阵的维数是四
	eigenJacobi m(4);
	m.no1();
	m.jacobi();
	m.output();

	cout<<endl;
	cout<<"按任意键继续..."<<endl;
	char c; 
    c=getch(); 

//======================================第二题===================================================
	
	cout<<"*****************************************************************************"<<endl;
	cout<<"                         程序说明（二）："<<endl<<endl;
	cout<<"第二题:编制用QR方法计算矩阵的特征值和特征向量的程序"<<endl
		<<"   其中  A=[ 2  -1   0   0"<<endl
		<<"             -1  2   -1  0"<<endl
		<<"             0   -1  2   -1"<<endl
		<<"             0   0   -1  2 ]"<<endl;
    cout<<"*****************************************************************************"<<endl;
    cout<<endl;

    //建立eigenQR的一个对象m，矩阵的维数也是四
	eigenQR n(4);
	n.no2();
	n.qr();
}