HessenbergTransform.cpp
上传用户:fxromeo
上传日期:2010-04-08
资源大小:89k
文件大小:1k
开发平台:

Visual C++

  1. //HessenbergTransform.cpp  
  2. //约化一般实矩阵为赫申伯格阵的初等相似变换法
  3. #include <iostream> //输入输出流头文件
  4. #include "Matrix.h" //矩阵类及相关函数头文件
  5. #include "EigenvalueVector.h" //特征值及特征向量头文件
  6. using namespace std; //名字空间
  7. void main() // 定义控制台应用程序的入口点
  8. {
  9. double a[5][5] = 
  10. {
  11. {  1.0,   6.0, -3.0, -1.0,  7.0},
  12. {  8.0, -15.0, 18.0,  5.0,  4.0},
  13. { -2.0,  11.0,  9.0, 15.0, 20.0},
  14. {-13.0,   2.0, 21.0, 30.0, -6.0},
  15. { 17.0,  22.0, -5.0,  3.0,  6.0} 
  16. };
  17. matrix<double> da(&a[0][0],5,5);
  18. cout.setf(ios::fixed); //输出数据为定点法
  19. cout.precision(6); //精度6位
  20. int iRet = HessenbergTransform(da);
  21. if(iRet > 0)
  22. {
  23. cout << "MATRIX A IS: " << endl;
  24. MatrixLinePrint(da);
  25. cout << endl;
  26. }
  27. }