一种求解无约束非线性规划的切线相交法
牛顿法与割线法求解非线性方程
牛顿法与割线法求解非线性方程在数学中，非线性方程是指方程中包含未知数的幂次大于等于2的项的方程。
求解非线性方程是数学中一个重要的问题，它在科学、工程和经济等领域中有着广泛的应用。
本文将介绍两种常用的非线性方程求解方法：牛顿法和割线法。
一、牛顿法牛顿法是一种迭代方法，用于求解非线性方程的根。
它基于泰勒级数展开的思想，通过不断迭代逼近方程的根。
牛顿法的基本思想是：选择一个初始值x0，然后通过迭代公式xn+1 = xn - f(xn)/f'(xn)，不断逼近方程的根。
具体步骤如下：1. 选择一个初始值x0；2. 计算函数f(x)在x0处的导数f'(x0)；3. 使用迭代公式xn+1 = xn - f(xn)/f'(xn)计算下一个近似解xn+1；4. 判断是否满足停止准则，如果满足，则输出近似解xn+1，算法结束；如果不满足，则将xn+1作为新的近似解，返回第2步继续迭代。
牛顿法的优点是收敛速度快，但缺点是对初始值的选择较为敏感，可能会陷入局部最优解。
二、割线法割线法也是一种迭代方法，用于求解非线性方程的根。
它与牛顿法类似，但是割线法不需要计算函数的导数。
割线法的基本思想是：选择两个初始值x0和x1，通过迭代公式xn+1 = xn - f(xn)(xn - xn-1)/(f(xn) - f(xn-1))，不断逼近方程的根。
具体步骤如下：1. 选择两个初始值x0和x1；2. 使用迭代公式xn+1 = xn - f(xn)(xn - xn-1)/(f(xn) - f(xn-1))计算下一个近似解xn+1；3. 判断是否满足停止准则，如果满足，则输出近似解xn+1，算法结束；如果不满足，则将xn+1作为新的近似解，返回第2步继续迭代。
割线法的优点是不需要计算函数的导数，但缺点是收敛速度相对较慢。
三、牛顿法与割线法的比较牛顿法和割线法都是求解非线性方程的有效方法，它们各有优缺点。
牛顿法的收敛速度较快，但对初始值的选择较为敏感；割线法不需要计算函数的导数，但收敛速度相对较慢。
6-3无约束非线性规划问题的求解
2. 搜索步长: k 取最优步长 , 即满足 f ( x k k d k ) min f ( x k d k ) 。
梯度法算法步骤：
1. 给定初始点 x 1 R n , 允许误差 0 , 令k 1 。 2. 计算搜索方向 d k f ( x k ) ;
3. 若 || d k || , 则停止计算， x k 为所求极值点； 否则，求最优步长 k
无约束非线性规划问题的求解方法分为解析法和直接 法两类。
解析法
需要计算函数的梯度，本节介绍最速下降法、共轭梯度 法、牛顿法、变尺度法等解析方法
直接法仅通过比较目标函数值的大小来移动迭代点。本章
主要介绍步长加速法等直接方法。
一般来说，无约束非线性规划问题的求解是通过一系 列一维搜索来实现。因此，如何选择搜索方向是解无约 束非线性规划问题的核心问题，搜索方向的不同选择， 形成不同的求解方法。
得最佳步长k
(f ( X k ))T f ( X k ) ( f ( X ( k ) ))T P ( k ) k T (f ( X k )) Af ( X k ) f ( X k )T H ( X k )f ( X k )
H( Xk ) A
若函数f(x)不是二次函数，可用上式得到k的近似值.
2 2 1 min f ( X ) 2 x x 例6-8 用最速下降法求解： ,2 初始点 X 0 1 1
1 精度 10 。
解 因为函数f(x)为二次函数，故可用上述公式计 算最佳步长k
4 x f ( X ) 1 2 x2
4 0 H(X ) 0 2
使得 f ( x k k d k ) min f ( x k d k )。
非线性方程的求解方法
非线性方程的求解方法一、引言在数学领域中，非线性方程是指未知量与其对自身的各次幂、指数以及任意函数相乘或相加得到的方程。
求解非线性方程是数学中一个重要而又具有挑战性的问题。
本文将介绍几种常见的非线性方程求解方法。
二、牛顿迭代法牛顿迭代法是一种经典的非线性方程求解方法，它利用方程的切线逼近根的位置。
设f(x)为非线性方程，在初始点x0附近取切线方程y=f'(x0)(x-x0)+f(x0)，令切线方程的值为0，则可得到切线方程的解为x1=x0-f(x0)/f'(x0)。
重复这个过程直到满足精确度要求或迭代次数达到指定次数。
三、二分法二分法是一种简单而又直观的非线性方程求解方法。
它利用了连续函数的中间值定理，即若f(a)和f(b)异号，则方程f(x)=0在[a, b]之间必有根。
根据中值定理，我们可以取中点c=(a+b)/2，然后比较f(a)和f(c)的符号，若同号，则根必然在右半区间，否则在左半区间。
重复这个过程直到满足精确度要求或迭代次数达到指定次数。
四、割线法割线法是一种基于切线逼近的非线性方程求解方法，它与牛顿迭代法相似。
由于牛顿迭代法需要求解导数，而割线法不需要。
设f(x)为非线性方程，在两个初始点x0和x1附近取一条直线，该直线通过点(x0,f(x0))和(x1, f(x1))，它的方程为y=f(x0)+(f(x1)-f(x0))/(x1-x0)*(x-x0)，令直线方程的值为0，则可得到直线方程的解为x2 = x1 - (f(x1)*(x1-x0))/(f(x1)-f(x0))重复这个过程直到满足精确度要求或迭代次数达到指定次数。
五、试位法试位法是一种迭代逼近的非线性方程求解方法。
它利用了函数值的变化率来逼近根的位置。
设f(x)为非线性方程，选取两个初始点x0和x1，然后计算f(x0)和f(x1)的乘积，如果结果为正，则根位于另一侧，否则根位于另一侧。
然后再选取一个新的点作为下一个迭代点，直到满足精确度要求或迭代次数达到指定次数。
无约束非线性优化算法
无约束非线性优化算法此类算法提供公式 2 的精确解。
但是，它们要耗费与 H 的几个分解成比例的时间。
因此，对于大型问题，需要一种不同方法。
文献（[42] 和 [50]）中提出了几种基于公式 2 的逼近和启发式方法建议。
Optimization Toolbox 求解器采用的逼近方法是将信赖域子问题限制在二维子空间 S 内（[39] 和[42]）。
一旦计算出子空间 S，即使需要完整的特征值/特征向量信息，求解公式 2 的工作量也不大（因为在子空间中，问题只是二维的）。
现在的主要工作已转移到子空间的确定上。
二维子空间 S 是借助下述预条件共轭梯度法确定的。
求解器将 S 定义为由 s1和 s2确定的线性空间，其中 s1是梯度 g 的方向，s2是近似牛顿方向，即下式的解或是负曲率的方向，以此种方式选择 S 背后的理念是强制全局收敛（通过最陡下降方向或负曲率方向）并实现快速局部收敛（通过牛顿步，如果它存在）。
现在，我们可以很容易地给出基于信赖域的无约束最小化的大致框架：1.构造二维信赖域子问题。
2.求解公式 2 以确定试探步 s。
3.如果 f(x + s) < f(x)，则 x = x + s。
4.调整Δ。
重复这四个步骤，直到收敛。
信赖域维度Δ 根据标准规则进行调整。
具体来说，它会在试探步不被接受（即f(x + s) ≥ f(x)）时减小。
有关这方面的讨论，请参阅[46] 和 [49]。
Optimization Toolbox 求解器用专用函数处理 f 的一些重要特例：非线性最小二乘、二次函数和线性最小二乘。
然而，其底层算法思路与一般情况相同。
这些特例将在后面的章节中讨论。
非线性优化与约束优化问题的求解方法
非线性优化与约束优化问题的求解方法非线性优化问题是在目标函数和约束条件中包含非线性项的优化问题。
约束优化问题是在目标函数中加入了一些约束条件的优化问题。
解决这些问题在实际应用中具有重要意义，因此研究非线性优化和约束优化问题的求解方法具有重要的理论和实际意义。
一、非线性优化问题的求解方法非线性优化问题的求解方法有很多，下面介绍几种常见的方法：1. 黄金分割法：黄金分割法是一种简单但有效的搜索方法，它通过不断缩小搜索范围来逼近最优解。
该方法适用于目标函数单峰且连续的情况。
2. 牛顿法：牛顿法利用目标函数的一阶和二阶导数信息来逼近最优解。
该方法收敛速度较快，但在计算高阶导数或者初始点选取不当时可能产生不稳定的结果。
3. 拟牛顿法：拟牛顿法是对牛顿法的改进，它通过逼近目标函数的Hessian矩阵来加快收敛速度。
拟牛顿法可以通过不同的更新策略来选择Broyden-Fletcher-Goldfarb-Shanno（BFGS）方法或者DFP方法。
4. 全局优化方法：全局优化方法适用于非凸优化问题，它通过遍历搜索空间来寻找全局最优解。
全局优化方法包括遗传算法、粒子群优化等。
二、约束优化问题的求解方法约束优化问题的求解方法也有很多，下面介绍几种常见的方法：1. 等式约束问题的拉格朗日乘子法：等式约束问题可以通过引入拉格朗日乘子来转化为无约束优化问题。
通过求解无约束优化问题的驻点，求得原始约束优化问题的解。
2. 不等式约束问题的罚函数法：不等式约束问题可以通过引入罚函数来转化为无约束优化问题。
罚函数法通过将违反约束条件的点处添加罚项，将约束优化问题转化为无约束问题。
3. 逐次二次规划法：逐次二次规划法是一种常用的求解约束优化问题的方法。
该方法通过依次处理逐个约束来逼近最优解，每次处理都会得到一个更小的问题，直至满足所有约束条件。
4. 内点法：内点法是一种有效的求解约束优化问题的方法。
该方法通过向可行域内部逼近，在整个迭代过程中都保持在可行域内部，从而避免了外点法需要不断向可行域逼近的过程。
非线性-无约束规划
6） 实用收敛性： ）
定义最优解集如下 S* = { x | x 具有某种性质 } 例：S*={x| x---g.opt} S*={x| x---l.opt} S*={x|∇f(x)=0} S*={x| f’(x)≤β} (β为给定实数,称为阈值） 当下列情况之一成立时 当下列情况之一成立时，称算法收敛具有该性质点 之一成立时， 1°∃x(k) ∈S*; ° 2°∀k，{X(k)}任意极限点∈S* ° 任意极限点∈ 任意极限点
* ak 为最优步长。 最优步长。 则称
根据单变量的驻点条件: 根据单变量的驻点条件 d f(xk+akPk)/dak=0 （当ak=ak* 时） 以及复合函数的求导法则可得： 以及复合函数的求导法则可得：
∇f ( x
k +1 T
) P =0
k
2） 缩小区间的非精确一维搜索
（1）单峰的概念 ） 若对任意λ 若对任意 1 ,λ2, α≤ 1º 若α2 ≤
停
11. 最优步长的一维搜索 1) 精确一维搜索（假定求目标函数极小值） 假定求目标函数极小值） * ak 是在给定 k和方向 是目标函数， 设f(X)是目标函数，如果 是在给定X 是目标函数 矢量P 通过f(x)=f(xk+akPk) 的极小化而产生 矢量 k下，通过
ak* = arg ak min f ( x k + ak P k )
∂ u ∂u ∂u ∂u = cos α + cos β + cos γ ∂ l ∂x ∂y ∂r
2. 海瑟矩阵
海瑟矩阵是对称形式：
∂2 f ( X ) ∂x12 ∂2 f ( X ) 2 H ( X ) = ∇ f ( X ) = ∂x2 ∂x1 ...... ∂2 f ( X ) ∂xn ∂x1
非线性-无约束规划
f(x) 为凹函数（严格凹函数）。
严格凸函数
凸函数
严格凹函数
x1
x 2
2.2 凸集、凸函数和凸规划(续)
定理： f(x) 为凸集 S 上的凸函数 S 上任 意有限点的凸组合的函数值不大于各点函 数值的凸组合。
△可行方_ 向：
设 x∈S，d∈Rn, d≠0, 若存在 0
_
使 x d S, (0, ) ，
称d 为该点的可行方向。
同时满足上述两个性质的方向称 下降可行方向。
迭代算法的停止标准
1）
|| X k1 X k || 1
或
||
X k 1 || X k
X ||
k
考虑（fs）
s.t. x∈S
常用一种线性搜索的方式构造{xk}序列来求解 迭代中从一点出发沿下降可行方向找一个
新的、更优的点。
△下降方向 ：
设 x _∈S，d ∈Rn,d≠0,若存在 ，0
使 在
_
_
x _f点(x的 下d )降 方f (x向),。 (0, )
，称d 为
4 常用的搜索算法结构
以及
4) 全局收敛： 对任意初始点x(1), 算法均收敛。
5) 局部收敛： 当x(1) 充分接近解x*时，算法才收敛。
2. 实用收敛性：
定义解集
S* = { x | x 具有某种性质 }
例：S*={x|x---g.opt} S*={x|x---l.opt}
S*={x| f(x)=0} S*={x|f′(x)≤β } (β为给定实数,称为阈值
xn2
无约束非线性规划求解方法及其实现
无约束非线性规划求解方法及其实现作者：杨玲指导老师：陈素根摘要:非线性规划是具有非线性约束条件或目标函数的数学规划，是运筹学的一个重要分支。
非线性规划属于最优化方法的一种，是线性规划的延伸。
非线性规划研究一个n元实函数在一组灯饰或不等式的约束条件下的极值问题，且目标函数和约束条件至少有一个是未知量的非线性函数。
目标函数和约束条件都是线性函数的情形则属于线性规划。
非线性规划是20世纪50年代才形成的一门新兴学科。
1951年H.W库恩和A.W塔克发表的关于最优性条件的论文是非线性规划正是诞生的一个重要标志。
在50年代还得出了可分离规划和二次规划的n种解法，它们大都是以G.B.丹齐克提出的解线性规划的单纯形法为基础的。
50年代末到60年代末出现了许多解线性规划问题的有效的算法，70年代又得到进一步的发展。
非线性规划在工程，管理，经济，科研，军事等发面都有广泛的应用，为最优设计提供了有力的工具。
20世纪80年代以来，随着计算机技术的快速发展，非线性规划在信赖域法、稀疏牛顿法、并行计算、内点法和有限存储法等领域取得了丰硕的成果，无约束非线性规划问题是非线性规划的一个重要内容，很多学者对非线性规划问题进行了深入且系统的研究，研究成果丰硕。
关键词最优化共轭梯度法非线性无约束1 引言1.1 无约束非线性规划问题是最基本的非线性规划问题，在1959~1963年幼三位数学家共同研究成功求解无约束问题的DFP变尺度法，该算法的研究成功是无约束优化算法的一个大飞跃，引起了一系列的理论工作，并陆续出现了许多新的算法。
20世纪80年代以来，随着计算机技术的快速发展，非线性规划在信赖域法、稀疏牛顿法、并行计算、内点法和有限存储法等领域取得了丰硕的成果。
无约束非线性规划问题是非线性规划的一个重要内容，很多学者对非线性规划问题进行了深入且系统的研究，研究成果丰硕。
1.2 本文主要研究无约束非线性规划问题，将文章分成四个部分，首先会具体介绍无约束非线性规划的相关概念，并在此基础上研究非线性规划的相关理论与基本算法问题，接着详细介绍无约束非线性规划的几种主要的求解方法，最后举例说明他在实际生活中的应用，并编程实现它。
凸优化之无约束优化（一维搜索方法：二分法、牛顿法、割线法）
凸优化之⽆约束优化（⼀维搜索⽅法：⼆分法、⽜顿法、割线法）1、⼆分法（⼀阶导）⼆分法是利⽤⽬标函数的⼀阶导数来连续压缩区间的⽅法，因此这⾥除了要求 f 在 [a0,b0] 为单峰函数外，还要去 f(x) 连续可微。
（1）确定初始区间的中点 x(0)=(a0+b0)/2 。
然后计算 f(x) 在 x(0) 处的⼀阶导数 f'(x(0))，如果 f'(x(0)) >0 , 说明极⼩点位于 x(0)的左侧，也就是所，极⼩点所在的区间压缩为[a0,x(0)]；反之，如果 f'(x(0)) <0，说明极⼩点位于x(0)的右侧，极⼩点所在的区间压缩为[x(0)，b0]；如果f'(x(0)) = 0，说明就是函数 f(x) 的极⼩点。
（2）根据新的区间构造x(1)，以此来推，直到f'(x(k)) = 0，停⽌。
可见经过N步迭代之后，整个区间的总压缩⽐为（1/2）N，这⽐黄⾦分割法和斐波那契数列法的总压缩⽐要⼩。
1 #ifndef _BINARYSECTION_H_2#define _BINARYSECTION_H_34 typedef float (* PtrOneVarFunc)(float x);5void BinarySectionMethod(float a, float b, PtrOneVarFunc fi, float epsilon);67#endif1 #include<iostream>2 #include<cmath>3 #include "BinarySection.h"45using namespace std;67void BinarySectionMethod(float a, float b, PtrOneVarFunc tangent, float epsilon)8 {9float a0,b0,middle;10int k;11 k = 1;12 a0 = a;13 b0 = b;14 middle = ( a0 + b0 )/2;1516while( abs(tangent(middle)) - epsilon > 0 )17 {18 #ifdef _DEBUG19 cout<<k++<<"th iteration:x="<<middle<<",f'("<<middle<<")="<<tangent(middle)<<endl;20#endif2122if( tangent(middle) > 0)23 {24 b0 = middle;25 }26else27 {28 a0 = middle;29 }30 middle =( a0+b0)/2;31 }3233 cout<<k<<"th iteration:x="<<middle<<",f'("<<middle<<")="<<tangent(middle)<<endl;34 }1 #include<iostream>2 #include "BinarySection.h"345float TangentFunctionofOneVariable(float x)6 {7return14*x-5;//7*x*x-5*x+2;8 }910int main()11 {12 BinarySectionMethod(-50, 50, TangentFunctionofOneVariable, 0.001);13return0;14 }1th iteration:x=0,f'(0)=-52th iteration:x=25,f'(25)=3453th iteration:x=12.5,f'(12.5)=1704th iteration:x=6.25,f'(6.25)=82.55th iteration:x=3.125,f'(3.125)=38.756th iteration:x=1.5625,f'(1.5625)=16.8757th iteration:x=0.78125,f'(0.78125)=5.93758th iteration:x=0.390625,f'(0.390625)=0.468759th iteration:x=0.195312,f'(0.195312)=-2.2656210th iteration:x=0.292969,f'(0.292969)=-0.89843811th iteration:x=0.341797,f'(0.341797)=-0.21484412th iteration:x=0.366211,f'(0.366211)=0.12695313th iteration:x=0.354004,f'(0.354004)=-0.043945314th iteration:x=0.360107,f'(0.360107)=0.041503915th iteration:x=0.357056,f'(0.357056)=-0.001220716th iteration:x=0.358582,f'(0.358582)=0.020141617th iteration:x=0.357819,f'(0.357819)=0.0094604518th iteration:x=0.357437,f'(0.357437)=0.0041198719th iteration:x=0.357246,f'(0.357246)=0.0014495820th iteration:x=0.357151,f'(0.357151)=0.0001144412、⽜顿法（⼆阶导）前提：f 在 [a0,b0] 为单峰函数，且[a0,b0] 在极⼩点附近，不能离的太远否则可能⽆法收敛。
求解非线性无约束优化问题的两种方法的分析
垫墼兰￡望叁兰堑圭兰垒篁塞１第一章预备知识§１．１共轭梯度方法§１．１．１引言共轭梯度法足最优化中最常用的方法之一。
它具有算法简便，存储需求小等优点，十分适合于大规模优化问题．在石油勘探，大气模拟，航天航空等领域出现的特大规模的优化问题是常常利用共轭梯度法求解。
在所有需要计算导数的优化方法中，最速下降是最简单的，但它速度太慢。
拟牛顿方法收敛速度很快，被广泛认为是非线性规划的最有效的方法。
但拟牛顿法需要存储矩阵以及通过求解线性方程组来计算搜索方向，这对于求解诸如上述问题等一些大规模问题几乎是不太可能办到的，共轭梯度法在算法的简便性，所需存储量等方面均与最速下降法差别不大，而收敛速度比最速下降法要快。
非线性共轭梯度法的收敛性分析的早期工作主要由Ｆｌｅｔｃｈｅｒ，Ｐｏｗｅｌｌ，Ｂｅａｌｅ等学者给出的，近年来，Ｎｏｃｅｄａｌ，Ｇｉｌｂｅｒｔ，Ｎａｚａｒｅｔｈ等学者在收敛性方面得到了不少的结果，使得共轭梯度法的研究由又热了起来．我国的学者也在共轭梯度法的理论研究中也取得了一定的成绩。
例如中科院应用数学所的韩继业，戴口等．§１．１．２共轭方向法共轭梯度法最本质的是共轭性质，共轭性是正交的一种推广。
定义１．１．２．１：设Ｗ∈咿×ｎ对称正定，ｄｌ，ｄ２，…，ｄ。
是咿中的一组非零向量，如果盯Ａｄｊ＝０，（ｉ≠Ｊ）．（１．１）则称ｄ１，ｄ２，…，ｄ。
是相互Ａ一共轭。
显然可见，如果ｄｌ，ｄ２，…，ｄ。
相互Ａ一共轭，则它们是线性无关的。
设Ｊ是单位阵则知，一共轭就是正交。
一般共轭方向法步骤如下：算法１．１．２．１：（一般共轭方向法）给出∞＋的初始点Ｘｌ，步ｌ：计算ｇｌ＝ｇ（Ｘ１）．步２：计算ｄｌ，使（ｆ｛’９ｌ＜０．步３：令女＝１．步４：计算口ｋ和Ｘｋ＋１，使得ｆ（ｘｋ－Ｆ‘１ｋｄｋ）。
Ｉ。
ｊ１１，‰十“呶），Ｘｋ＋１＝Ｘｋ＋ｖ。
ｋｄｋ．步５：计算以＋ｌ使得ｄ矗１Ｇｄｊ＝０，Ｊ＝１，２，…ｋ．步６：令ｋ：＝ｋ＋１，转步４．共轭方向法的一个基本性质是：只要执行精确线性搜索，就能得到二次终止性，这就足下面的共轭方向法基本定理。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
函数
优 化 方法
最 小 点
最小 值
迭代 次 数
０算 ，６ 计
转 Ｓｅ ２ ｔ ｐ
器
｝４ 次 值 。 ｏ５ ８．。。。 一 插 法 ‘ ８ ６４３４ 。 手三 ６ ０ ５ ４ ６ 。 ５ ２
（） （ ） ｏ 一 ｏ 相交 于点 （ｏ ｂ） 。 ６ ｏ）
（） ｆ ａｘ ＋
其 中
图 １ 算法 原 理 图
。若 － ｂ） Ｏ 则 ｂ 厂（ ｏ＝ ， 。
就 是最 优点 ； 厂（ ｏ＜ ， ａ ｂ，＝ ｏ若 （ｏ＞ ， ｂ ｂ； 若 ｂ） ０ 令 ＝ 。ｂ ａ； ｂ）Ｏ 令 ＝ 。
［，］ 。６ 上有一 阶连续 导数 。
就是 最优 点 ； 厂（ｏ ０ 令 ａ ｂ； （ ）０令 ａａ，＝０ 若 ｂ） ， ＝０ ＜ 若 ｂ ＞ ， ＝ｏｂ ｂ； ｏ 如果ｆ（ ＞ ， 作过 点 （叶 ６／ ， ｈ 的切 线 的平行 线 ）０则 （ ）２ｙ ） ｏ ） ， 中 ， ｐ （）＋ ６ ０ － ｂ （ ）／ 与 ， 其 ｙ ｂ ） ） ｂ ｆ（ ）叶６ ）２
个标 准 函数 的测 试表 明 了该 算 法的有 效性 。
关键 词 ： 单峰 函 数 ； 线 相 交 法 ； 约 束 非 线 性 规 划 切 无 中 图 分 类 号 ：Ｐ ０ Ｔ３１ 文献标识 码 ： Ａ 文 章 编 号 ：６ ２ ７ ０ （０ ０ ０ — ０ ５ ０ １７ — ８０ ２ １ ） ９ ０ ４ — ２
・
４ ・ ６
软 件 导 刊
表 １ ３种 函数 的 测 试 结 果与 对 比
２ １ 年 ００
Ｓ ｐ ： Ｉ － ， ｔ ４若 ａＩ 则输出 （ ）２停机 ｅ ｂ 叶６／ ， 若 ｆ ｎ ｆ ｓ 则输出 ｎ停机 ，（） ＜ ， ， 若 ｌ ６ ＜， 厂（）ｌ ｓ 则输出 ６ 停机 ，
线 方程为 ：
Ｓｅ２ 若 ｌ Ⅱ ＞ 贝 输出 ０ ｔ ： 厂（）ｌ ０ ４ ｐ 若Ｉ ６ ＜ － ）Ｉ ０则输出 ６ 尸（
ｓ ： 计算
作者 简介 ： 国建 （９ ４ ， ， 程 １６ 一） 男 陕西杨 凌人 ， 士 ， 博 西安 石油 大 学计 算机 学 院教授 ， 究 方 向为神 经 网络 、 化 计 算、 计算 方 法在 油藏描 述 中的应 研 进 软 用； 卢飞远 （９ ３ ） 男， 西成 阳人 ， 士 ， １８一 ， 陕 硕 西安石 油 大学计 算机 学院助 教 ， 究 方 向为 图像 处理 、 式 识 别 ； 华 （９ ９ ， ， 东 日照 研 模 房 １７ 一） 女 山 人， 士， 硕 西安石 油 大学计算 机 学院助教 ， 究方 向为神 经 网络 、 式识 别 。 研 模
０ 引言
本 文 首 先 介 绍 切 线 相 交 法 的 基 本 思 想 和 原 理 ， 导 其 算 法 推
ｆ＝ Ⅱ ａ－ 口 ｙｆ（）七 ）ａ ） ￣ ｆ（ （ｂ ６ ｂ 一，（ ） ｙ＝ ） ｆ（ ） ６ ｂ
解交 为ｏ 占 其点 ｎ ＝
点 （ 。 ） （ ６ ） 直 线 。 。 ）和 ｂ ）的
如果 ，（ ）Ｏ 则作 过点 （叶６／ ， ｈ 的切线 ｙ的平 行线 ａ＜ ， ｏ （ ）２ｙ ） ａ ｎ ｙｐ ｙｐ 厂（）＋７ｂ ｔ 口－ Ⅱ （＋ ） ２又 与 ｙ （）ｔ ） ａ ， 。 （（） ） ａ ｎ ｆ（）ａ ｂ ） ／ ６ ｂ ｂ
一
１ 算 法 基本 原 理
如图 １ 在 给定 的搜 索 区间 ［ ， ］ ， 点 处 的切 线为 ｙ ，ｂ ， ０ｂ 上 端 ａ ｙ 相交 于一点 ａ． ｏ且两 条切 线 的方程为 ：
重复上述迭代过程 ， 直到 ｆ － （ ａｌ 为精度控制参数 ） ｂ
并 输 出近试 最优值 ＝ 叶６ ／ 。 （ ）２
4函数测试与对比通过上述对3种函数的测试可以看出切线相交法在迭代次数和精确度上都优于平分法和黄精分割法甚至在复杂函数的迭代次数上和切线法2次插值法一样优于3次插值法而运算的时间却相当
第 ９ 第９ 卷 期
２ １ 年 ９月 ００
软 件 导 刊
Ｓｏｔ ｒ ｉｅ ｆｗａｅ Ｇｕｄ
Ｖｏ． ． １ ＮＯ９ ９ Ｓ ｐ． ｅ ２０１ ０
（） １
由 于 （ 口 ） （ （ ） 分 别 是 过 （ ） 两 条 切 线 ， 将 其 Ⅱ ） ，ｂ ６ ） １这 则
２ 算 法 步 骤
Ｓｅ １ 给 定 搜 索 区 间 ［，］设 定 精 度 控 制 参 数 ｅ 得 ２ （ ‘， 代 （， 到 两 切 人１ ｆ 。６ ＿并 到１ 得 这 条 ） 【 ｂ） Ⅱ 会 ） 可 一 ｂ） ｂ 一尸（ 可 ＝ ｆ
在求解 无 约束非 线性规 划 的单峰 函数 问题 时 ， 一 维探 索 若
（ ） 交 于 点 （ｏ ｂ） 。 ｂ相 ６ ｏ）
￣ ｅ ｏ ｏｂ＝
。若 （ｏ＝ ， ｂ ｂ）０ 则 。
区间 为 ［ ， ］ ｎ ６ 已经确 定 ， ， ）其 中／为实 值 连续 函数 ， 在 ｍ ， 且
一
种 求解无约束非线性规划 的切线相 交法
程 国建 ， 飞远 ， 卢 房 华
（ 西安石 油 大学 计 算机 学 院 ， 陕西 西安 ７ ０ ６ ） １０ ５
摘
一
要 ： 线相 交 法是基 于切线 法的 思想 并结合 黄金 分割 法 而应 用于 求解无 约束非 线性规 划 问题 的最优 解。提 出了 切 个 新的 求解 无约束 非线性 规 划的有 效算 法 ， 并在搜 索区间 已经 确定 和适 当的条件 下证 明 了该 算法 的收 敛性 ， 对几
。
其 中 ，ａ ｙ ｂ是 过
形成 过程 ； 次 ， 其 给出算 法 的一般求 解步 骤 ， 以及 算法 的收敛 性 证明； 再次 用一 些 函数 来测 试算 法 的正 确性 以及 和其他 经典 一 维搜 索算法对 比 ； 最后 对本 文进行 总结 以及 对算 法作进 一步 的
研 究
如果 － ａ）０， ａ 就是 最优点 ； 厂（ｏ＝ 则 ｏ