0043算法笔记【随机化算法】解非线性方程组.docxVIP

问题描述? ? ?求解下面的非线性方程组? ??其中，x1,x2,…,xn是实变量，fi是未知量x1,x2,…,xn的非线性实函数。要求确定上述方程组在指定求根范围内的一组解。? ??问题分析? ? ?解决这类问题有多种数值方法，如：牛顿法、拟牛顿法、粒子群算法等。最常用的有线性化方法和求函数极小值方法。为了求解所给的非线性方程组，构造一目标函数? ???式中，x=(x1,x2,……xn)。易知，上式取得极小值点即是所求非线性方程组的一组解。? ??求解思路? ? 在指定求根区域D内，选定一个随机点x0作为随机搜索的出发点。在算法的搜索过程中，假设第j步随机搜索得到的随机搜索点为xj。在第j+1步，计算出下一步的随机搜索增量dxj。从当前点xj依dxj得到第j+1步的随机搜索点。当x时，取为所求非线性方程组的近似解。否则进行下一步新的随机搜索过程。? ??题外话：笔者在读王晓东《算法设计与分析》中这一节时，发现书上所给的代码似乎有些问题。在这里指出，如果提得不对，还请大侠们拍砖。书中给出的代码具体如下：[cpp]?view plain?copybool?nonlinear(double?*x0,double?*dx0,?double?*x,double?a0,????????double?epsilon,?double?k,int?n,?int?steps,int?m)??{???????static?randomnumber?rad;??????bool?success;??????double?*dx,?*r;??????dx=new?double[n+1];??????r=new?double[n+1];??????int?mm=0;??????int?j=0;??????double?a=a0;??????for(int?i=1;i=n;i++)??????{??????????x[i]=x0[i];??????????dx[i]=dx0[i];??????}??????double?fx=f(x,n);??????double?min=fx;??????while?((minepsilon)(jsteps))??????{??????????if(fxmin)??????????{??????????????min=fx;??????????????a*=k;??????????????success=true;???????????}??????????else??????????{??????????????mm++;??????????????if?(mmm)a/=k;??????????????success=false;??????????}??????????for(int?i=1;i=n;i++)??????????????r[i]=2.0*rnd.fRandom()-1;??????????if(success)??????????????for(int?i=1;in;i++)??????????????????dx[i]=a*r[i];??????????else???????????for?(int?i=1;in;i++)??????????????dx[i]=a*r[i]-dx[i];??????????for?(int?i=1;in;i++)??????????????x[i]+=dx[i];??????????fx=f(x,n);????????}?????if(fx=epsilon)?????????return?ture;?????else?????????return?false;??}??? ? ??问题1：while ((minepsilon)(jsteps))句中，循环控制变量jsteps，但在循环体中，j没有自加，while循环永远不会结束。? ? ?问题2：在循环体结束时，有赋值语句：fx=f(x,n); ? 笔者发现如果此时fxepsilon时，程序while循环体会继续执行，到while开始几句是会更新min的值，可循环不会跳出，而会继续往下执行，更新数组x的值，重新计算目标函数f(x)。这样即时程序找到方程的解，也不会停下来。? ? 笔者纠正后具体程序如下：[cpp]?view plain?copy//随机化算法?解线性方程组??#include?stdafx.h??#include?RandomNumber.h??#include?iostream??using?namespace?std;????bool?NonLinear(double?*x0,double?*dx0,double?*x,double?a0,??????????????????????double