开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 数学计算 > 查看源码
NLequation.h
资源名称：科学与工程数值计算算法配套源码vc++.zip　[点击查看]
上传用户：weigute
上传日期：2007-03-02
资源大小：1287k
文件大小：4k
源码类别：

数学计算

开发平台：

Visual C++

NLequation.h：源码内容
  1. //////////////////////////////////////////////////////////////////////
  2. // NLequation.h
  3. //
  4. // 求解线性方程组的类 CNLequation 的声明接口
  5. //
  6. // 周长发编制, 2002/8
  7. //////////////////////////////////////////////////////////////////////
  9. #if !defined(AFX_NLEQUATION_H__72143033_E39F_404E_9A3B_F870E603EEE2__INCLUDED_)
  10. #define AFX_NLEQUATION_H__72143033_E39F_404E_9A3B_F870E603EEE2__INCLUDED_
  12. #if _MSC_VER > 1000
  13. #pragma once
  14. #endif // _MSC_VER > 1000
  16. #include <math.h>
  18. // 需要调用复数和矩阵运算
  19. #include "Matrix.h"
  20. #include "Complex.h"
  22. // 类声明
  23. class CNLequation  
  24. {
  25. //
  26. // 公有接口函数
  27. //
  28. public:
  30. //
  31. // 构造与析构
  32. //
  34. CNLequation();
  35. virtual ~CNLequation();
  37. //
  38. // 虚函数：计算方程左端函数值，必须在引申类中覆盖该类函数
  39. //
  41. virtual double Func(double x)
  42. {
  43. return 0.0;
  44. }
  45. virtual double Func(int n, double x[])
  46. {
  47. return 0.0;
  48. }
  49. virtual void Func(double x, double y[])
  50. {
  51. }
  52. virtual double Func(double x, double y)
  53. {
  54. return 0.0;
  55. }
  56. virtual double Func(double x[], double y[])
  57. {
  58. return 0.0;
  59. }
  60.     virtual void FuncMJ(int n, double x[], double p[])
  61. {
  62. }
  64. //
  65. // 非线性方程求解算法
  66. //
  68. // 求非线性方程实根的对分法
  69. int GetRootBisect(int nNumRoots, double x[], double xStart, double xEnd, double dblStep, double eps = 0.000001);
  70. // 求非线性方程一个实根的牛顿法
  71. BOOL GetRootNewton(double* x, int nMaxIt = 60, double eps = 0.000001);
  72. // 求非线性方程一个实根的埃特金迭代法
  73. BOOL GetRootAitken(double* x, int nMaxIt = 60, double eps = 0.000001);
  74. // 求非线性方程一个实根的连分式解法
  75. BOOL GetRootPq(double* x, double eps = 0.000001);
  76. // 求实系数代数方程全部根的QR方法
  77. BOOL GetRootQr(int n, double dblCoef[], double xr[], double xi[], int nMaxIt = 60, double eps = 0.000001);
  78. // 求实系数代数方程全部根的牛顿下山法
  79. BOOL GetRootNewtonDownHill(int n, double dblCoef[], double xr[], double xi[]);
  80. // 求复系数代数方程全部根的牛顿下山法
  81. BOOL GetRootNewtonDownHill(int n, double ar[], double ai[], double xr[], double xi[]);
  82. // 求非线性方程一个实根的蒙特卡洛法
  83. void GetRootMonteCarlo(double* x, double xStart, int nControlB, double eps = 0.000001);
  84. // 求实函数或复函数方程一个复根的蒙特卡洛法
  85. void GetRootMonteCarlo(double* x, double* y, double xStart, int nControlB, double eps = 0.000001);
  87. //
  88. // 非线性方程组求解算法
  89. //
  91. // 求非线性方程组一组实根的梯度法
  92. BOOL GetRootsetGrad(int n, double x[], int nMaxIt = 500, double eps = 0.000001);
  93. // 求非线性方程组一组实根的拟牛顿法
  94. BOOL GetRootsetNewton(int n, double x[], double t, double h, int nMaxIt = 500, double eps = 0.000001);
  95. // 求非线性方程组最小二乘解的广义逆法
  96. BOOL GetRootsetGinv(int m, int n, double x[], double eps1 = 0.000001, double eps2 = 0.000001);
  97. // 求非线性方程组一组实根的蒙特卡洛法
  98. void GetRootsetMonteCarlo(int n, double x[], double xStart, int nControlB, double eps = 0.000001);
  100. //
  101. // 内部函数
  102. //
  103. private:
  104. void g60(double* t,double* x,double* y,double* x1,double* y1,double* dx,double* dy,double* p,double* q,int* k,int* it);
  105. void g90(double xr[],double xi[],double dblCoef[],double* x,double* y,double* p,double* q,double* w,int* k);
  106. void g65(double* x,double* y,double* x1,double* y1,double* dx,double* dy,double* dd,double* dc,double* c,int* k,int* is,int* it);
  107. void g60c(double* t,double* x,double* y,double* x1,double* y1,double* dx,double* dy,double* p,double* q,int* k,int* it);
  108. void g90c(double xr[],double xi[],double ar[],double ai[],double* x,double* y,double* p,double* w,int* k);
  109. void g65c(double* x,double* y,double* x1,double* y1,double* dx,double* dy,double* dd,double* dc,double* c,int* k,int* is,int* it);
  110. double rnd(double* r);
  111. };
  113. #endif // !defined(AFX_NLEQUATION_H__72143033_E39F_404E_9A3B_F870E603EEE2__INCLUDED_)