非线性方程组求根算法概述.docVIP

 PAGE \* MERGEFORMAT 3 非线性方程组求根算法概述 -非线性方程组求根方法（5） 下面介绍一个求非线性方程组零点的迭代方法，它是一种类牛顿求根方法。为何说成是“类”呢，因为牛顿法除函数值外，还用到了梯度信息；而本法只用到函数值，是一种较简单且实用的算法。笔者30年前在 nth做访问研究时自己编制的，当时还真解决了実际问题，用之于7阶、12阶、18阶的复杂动力学系统，均得到了有意义的根值。 方程组 F=F(X)=0；X=[x1,x2,…….,xn]T Fi=Fi(x1,x2,…….xn)=0 ^^^^^^^^^^^^^^^^^^^^^ F1=F1(x1,x2,…….xn)=0 F2=F2(x1,x2,…….xn)=0 ① …………………………………….. Fn=Fn(x1,x2,…….xn)=0 改造后的方程组 改造的目的是使迭代中连续应用的JACOB矩阵尽量避免出现奇异的情况。 xi==xi+Fi=?i ^^^^^^^^^^^^^^^^^^ x1=?1(x1,x2,….,xn) x2= ?2(x1,x2,….,xn) ② ………………….. …….. xn=?n(x1,x2,….,xn) 迭代的过程 设迭代本次为(k),状态变量为x1(k),x2(k),….,xn(k),则下一次迭代为（k+1），状态变量是x1(k+1),x2(k+1),….,xn(k+1)。认为xi(k+1),i=1,2,…., n为方程组的根，则它们使原方程组两端恒等，即满足下式： xi(k+1)=?i(x1(k+1),x2(k+1),……,xn(k+1)), i=1,2,….,n 将上式右端展开成Taylor级数，且只取线性项，有 xi(k+1)=?ik+j=1j=n??i?xjxjk+1-xj(k), i=1,2,……,n ③ 其中 ??i?xj =JACOB(i,j), i=1,2,….,n; j=1,2,…..,n ④ 是雅可比矩阵。 将③式化成下列形式： ??1?x1-1???1?xn?????n?x1???n?xn-1 x1k+1??xn(k+1)=-?1k+j=1j=n??1?xjxj(k)?-?nk+j=1j=n??n?xjxj(k) ⑤ 上式可简化为 JACOBi,j-Inxjk+1= -?jk+JACOBi,jxj(k), i=1,2,…,n ⑥ 或写成 AX(k+1) = B ⑦ 其中， In是nxn阶单位矩阵: In = 1?0???0?1 矩阵A(n,n)及列向量B(n) 的各元分别是： A(i,j) = JACOBi,j-In, i=1,2,….,n; j=1,2,…???n ⑧ B(i)= -?ik+JACOBi,jxj(k), i=1,2,….,n; ⑨ 至于状态向量X(n)在相连的两个迭代时刻的分量则为： xi(k), i=1,2,….,n ; xi(k+1), i=1,2,….,n 上式中的A,B之各元都取k 时刻的值。 利用高斯消元法，解线性代数方程组⑦，即完成从X(k)→X(k+1)的迭代过程。达到收敛时，xi(k+1),i=1,2,…,n就是系统的根。 收敛标准 TOL (xi(k+1)-xi(k))xi(k+1≦TOL ， i=1,2,…,n ⑩ JACOB 的数值近似 ??i?xj = ?ix1,x2,…xj1+?,…,xn-?ix1,x2,…x