实验六_无约束优化.pdfVIP

数学实验 无约束优化 林智慧 计03 2010011292 2013 春 数学实验实验报告 无约束优化 Lab6_无约束优化 实验目的： 1．掌握无约束规划的基本方法，对不同算法进行初步分析、比较。 2．练习建立和求解实际问题的无约束优化模型和非线性最小二乘拟合模型。 实验内容： 2.取不同的初值计算下列非线性规划，尽可能求出所有局部极小点，进而 找出全局极小点，并对不同算法 （搜索方向、搜索步长、数值梯度与分析梯度等）的 结果进行分析、比较。 （1） 2 2 5 2 min z (x x ) (1x ) [1x x (1x ) ] 1 2 1 1 2 1 实验过程： 设 2 2 5 2 f (x) (x x ) (1x ) [1  x x (1x ) ] 1 2 1 1 2 1 根据此函数的形式可以知道函数值恒大于或等于0，故求解f (x) 0 即求出 最优解，观察可得x1  0 或x1 1 或x2  0 时，函数值为0 。此函数梯度向量 和Hessian矩阵表达式相对比较复杂，通过求梯度的方法求解本题会十分麻 烦。 利用MATLAB工具求解 1、输出函数图象 程序如下： clear;clc; [x1,x2]=meshgrid(-5:0.01:5,-5:0.01:5); f=((x1.*x2).^2).*((1-x1).^2).*((1-x1-x2.*((1-x1).^5)).^2); figure(1); mesh(x1,x2,f) hold on; figure(2); contour(x1,x2,f,30) hold on; 与函数相对应，易于理解。 输出图像如下： 2013 春 数学实验实验报告 无约束优化 函数在(x1,x2)都在[-5,5]上取值的三维图像 与之对应的等高线示意图 从上面图中可知，我们所画范围较大，使细节不是很突出。远离解题目标。故 缩小范围到[-0.1,1.2]，画出图像如下： 2013 春 数学实验实验报告 无约束优化 从上图可以看出最优解的大致分布。 2013 春 数学实验实验报告 无约束优化 2、采用不同算法进行计算求解 定义基本函数f21 function y=f21(x) y=((x(1).*x(2)).^2).*((1-x(1)).^2).*((1-x(1)-x(2).*((1-x(1)).^5)).^2); end 取起始点为(0.5.0.1) 数值梯度：