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

Visual C++

  1. //EigenvalueVectorRealSymmetryJacobi.cpp  
  2. //实对称阵特征值及特征向量雅可比法
  3. #include <iostream> //输入输出流头文件
  4. #include "Matrix.h" //矩阵类及相关函数头文件
  5. #include "EigenvalueVector.h" //特征值及特征向量头文件
  6. using namespace std; //名字空间
  7. void main() // 定义控制台应用程序的入口点
  8. {
  9. double a[3][3] = 
  10. {
  11. { 2.0, -1.0,  0.0},
  12. {-1.0,  2.0, -1.0},
  13. { 0.0, -1.0,  2.0}
  14. };
  15. matrix<double> v(3,3), da(&a[0][0],3,3);
  16. double eps = FLOATERROR;
  17.    cout.setf(ios::fixed); //输出数据为定点法
  18. cout.precision(6); //精度6位
  19. int k = EigenvalueVectorRealSymmetryJacobi(da, v, eps, 100);
  20. if(k > 0)
  21. {
  22. cout << "MATRIX A IS: " << endl;
  23. MatrixLinePrint(da);
  24. cout << endl;
  25. cout << "MATRIX V IS: " << endl;
  26. MatrixLinePrint(v);
  27. cout << endl;
  28. }
  29. double dmb[3][3] = 
  30. {
  31. {  3.0, -3.0, -2.0 },
  32. { -3.0,  8.0,  4.0 },
  33. {  -2.0, 4.0,  3.0 } 
  34. };
  35. }