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

Visual C++

  1. // OrdinaryDifferentialEguation.h 求解常微分方程(组)头文件
  2. // Ver 1.0.0.0
  3. // 版权所有(C) 何渝, 2002
  4. // 最后修改: 2002.5.31.
  5. #ifndef _ORDINARYDIFFERENTIALEGUATION_H //避免多次编译
  6. #define _ORDINARYDIFFERENTIALEGUATION_H
  7. #include <LinearEquation.h> //线性方程组头文件
  8. #include <valarray> //模板类数组的标准头文件
  9. #include <Matrix.h> //模板类矩阵头文件
  10. #include <comm.h> //公共头文件
  11. using namespace std; //名字空间
  12. //定步长全区间积分的欧拉法
  13. template <class _Ty>
  14. void ODE_EulerContentStep(_Ty t, valarray<_Ty>& y,  
  15. _Ty h, int k, matrix<_Ty>& z);
  16. //变步长积分欧拉法
  17. template <class _Ty>
  18. void ODE_EulerVariationalStep(_Ty t, _Ty h, valarray<_Ty>& y, _Ty eps);
  19. //全区间定步长积分维梯法
  20. template <class _Ty>
  21. void ODE_WittyContentStep(_Ty t, valarray<_Ty>& y,  
  22. _Ty h, int k, matrix<_Ty>& z);
  23. //全区间定步长积分龙格-库塔法
  24. template <class _Ty>
  25. void ODE_RungeKuttaContentStep(_Ty t, valarray<_Ty>& y,  
  26. _Ty h, int k, matrix<_Ty>& z);
  27. //变步长积分龙格-库塔法
  28. template <class _Ty>
  29. void ODE_RungeKuttaVariationalStep(_Ty t, _Ty h, 
  30. valarray<_Ty>& y, _Ty eps);
  31. //变步长积分一步基尔法
  32. template <class _Ty>
  33. void ODE_GillVariationalStep(_Ty t, _Ty h, valarray<_Ty>& y,
  34. _Ty eps, valarray<_Ty>& q);
  35. //全区间变步长积分基尔法
  36. template <class _Ty>
  37. void ODE_GillVariationalStepInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  38. _Ty eps, int k, matrix<_Ty>& z);
  39. //变步长积分默森法
  40. template <class _Ty>
  41. void ODE_MersonVariationalStep(_Ty t, _Ty h, int n, valarray<_Ty>& y,
  42. _Ty eps, int k, matrix<_Ty>& z);
  43. //积分一步连分式法
  44. template <class _Ty>
  45. void ODE_Fraction(_Ty t, _Ty h, valarray<_Ty>& y, _Ty eps);
  46. //全区间积分连分式法
  47. template <class _Ty>
  48. void ODE_FractionInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  49. _Ty eps, int k, matrix<_Ty>& z);
  50. //全区间积分双边法
  51. template <class _Ty>
  52. void ODE_BothSidesInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  53. _Ty eps, int k, matrix<_Ty>& z);
  54. //全区间积分阿当姆斯预报-校正法
  55. template <class _Ty>
  56. void ODE_AdamsInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  57. _Ty eps, int k, matrix<_Ty>& z);
  58. //全区间积分哈明法
  59. template <class _Ty>
  60. void ODE_HammingInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  61. _Ty eps, int k, matrix<_Ty>& z);
  62. //积分一步特雷纳法
  63. template <class _Ty>
  64. void ODE_Treanor(_Ty t, _Ty h, valarray<_Ty>& y);
  65. //全区间积分特雷纳法
  66. template <class _Ty>
  67. void ODE_TreanorInterval(_Ty t, _Ty h, valarray<_Ty>& y,
  68. int k, matrix<_Ty>& z);
  69. //积分刚性方程组吉尔法
  70. template <class _Ty>
  71. int ODE_Gear(_Ty a, _Ty b, _Ty hmin, _Ty hmax, _Ty h, _Ty eps, 
  72.  valarray<_Ty>& y0, int k, valarray<_Ty>& t, matrix<_Ty>& z);
  73. //二阶微分方程边值问题数值解法
  74. template <class _Ty>
  75. void ODE_LinearBoundaryValude(_Ty a, _Ty b, _Ty ya, _Ty yb, int n, 
  76. valarray<_Ty>& y);
  77. #include "OrdinaryDifferentialEguation.inl" //类及相关函数的定义头文件
  78. #endif // _ORDINARYDIFFERENTIALEGUATION_H