大学数学实验七_无约束优化1.pdfVIP

大学数学实验七 无约束优化 实验报告 【实验目的】 1、掌握MATLAB 优化工具箱的基本用法，对不同算法进行初步分析、比较。 2 、练习用无约束优化方法建立和求解实际问题的模型 （包括非线性最小二 乘拟合）。 【实验内容】 2 取不同的初值计算下列非线性规划，尽可能求出所有局部极小点，进而找出全局极小点， 并对不同算法（搜索方向、搜索步长、数值梯度与分析梯度等）的结果进行分析、比较。 （1） 设 函数 是一个二变量的、乘积形式的函数，而且可以预见，其梯度向量和Hessian 矩 阵的表达式比较复杂，故若通过求梯度的方法求解本题会十分麻烦。 从这个函数表达式可以看出 恒 故若某个 能使 ，应该就是最优解。从表达式容易观察得，当 或 或 时， 。 下面用MATLAB 求解这个问题。 一、输出三维图像和等高线，直观观察最优解所在位置 首先，先让MATLAB 输出 的三维图像，直观地观察最优解的大致范围。程序如下。 上面的程序中，x 和y 分别表示 和 ，z 表示f (x , x ) 的值。自变量的范围均取在[-2, 2] 1 2 上，根据输出的图像进一步缩小范围。 输出的三维图像如下。 1 图1.1 三维图像 图1.2 等高线 7 第一幅图中，最大的函数值的数量级已到10 ，说明画图范围取得过大，函数值为0 的 电在该图中均呈现一平面，故需要缩小自变量的范围再绘图。修改程序如下。 2 输出的三维图像和等高线如下。 图1.3 三维图像 图1.4 等高线 从上面这幅等高线图可以看出，最优解大致分布在直线 或 或 附近。 3 二、对不同算法的结果进行分析比较 （数值方法计算梯度） 首先编制函数文件fun21 。 取起始点为(0.4, 0.2)进行计算。 1、搜索方向：BFGS ；搜索步长：混合二三次插值 输出结果如下。 x1 = 最优解x =0.4185 1 0.4185 0.0000 最优解x =0.0000 2 v1 = 5.9926e-015 最优值=5.9926e-015