一种基于拟Newton法解非线性方程组的数值算法课程.docVIP

一种基于拟Newton法解非线性方程组的数值算法 (吴志华 中国计量学院信息工程学院 310018)关键字:非线性方程组 牛顿法1?理论背景我们先考虑线性方程,线性方程组的解便不难得出了与线性方程相比,非线性方程问题无论是从理论上还是从计算公式上,都要复杂得多,计算方程的根既无一定章程可寻也无直接法可言,求解高次方程组的根,求解含有指数和正弦函数的超越方程的零点,牛顿(I . Newton)提出了方程求根的一种迭代方法,被后人称为牛顿算法,人们一直用牛顿算法,改善牛顿算法,不断推广算法的应用范围,可以说是数值计算方面的最有影响的计算方法对于方程式,如果是线性函数,则它的求根是容易的,其基本思想是将非线性方程式逐步归结为某种线性方程来求解2?主要理论考虑方程组 其中均为多元函数, 就可写成 当,且中至少有一个是自变量 的非线性函数时,则称方程组为非线性方程组求根的直接推广,实际上只要把单变量函数看成向量函数则可将单变量方程求根方法推广到方程组的一个近似根 将函数的分量在用多元函数泰勒展开,并取其线性部分,则可表示为 令上式右端为零,得到线性方程组 其中 称为为雅可比(Jacobi)矩阵,并记解为,则得 这就是解非线性方程组的牛顿法3?算法 牛顿法主要思想是用 进行迭代 的雅可比矩阵,再求过求出它的逆,当它达到某个精度时即停止迭代 具体算法如下:首先宏定义方程组 ,确定步长和精度; 求的雅可比矩阵; 可用求出雅可比矩阵;求雅可比矩阵的逆; 将右乘一个单位矩阵,通过单位矩阵变换实现求 的逆,用指针来存贮雅可比矩阵与其逆的相乘;用来迭代;当精度大于时,重复执行2——5步,直到精度小于或等于停止迭代,就是最后的迭代结果其中4?代码#include iostream.h #include stdlib.h #include math.h #include conio.h #define f0(x1,x2) (x1+2*x2-3) #define f1(x1,x2) (2*x1*x1+x2*x2-5) #define x_ 0.000001 #define matrixNum 2 double *matrixF2(double *x); int y=0; void main() { int i,j,n; double p,*x; double *b; double *matrixF; //矩阵F double *matrixF_; //矩阵F的雅可比矩阵的逆 b=(double *)malloc(matrixNum); matrixF=(double *)malloc(matrixNum); matrixF_=(double *)malloc(matrixNum*matrixNum); cout请输入初值:; for(i=0;imatrixNum;i++) cin*(x+i); do { p=0.0; for(i=0;imatrixNum;i++) *(b+i)=0; *matrixF=f0(*x,*(x+1)); *(matrixF+1)=f1(*x,*(x+1)); matrixF_=matrixF2(x); for(i=0;imatrixNum;i++) { for(j=0;jmatrixNum;j++) *(b+i)+=*(matrixF_+i*matrixNum+j)*(*(matrixF+j)); *(x+i)=*(x+i)-*(b+i); cout*(x+i) ; } coutendl; for(i=0;imatrixNum;i++) p+=pow(*(b+i),2); y++; }while(sqrt(p)x_); cout停止迭代,最终迭代结果为*x,*(x+1)endl; delete [] matrixF; delete [] matrixF_; getch(); } double *matrixF2(double *x) { int i,j; double t; double *matrixF1; //矩阵F的雅可比矩阵