数学建模--最优化实验.ppt

实验六 最优化实验 最佳水槽断面问题 （矩形断面） 用宽 l = 24 cm的长方铁板折成一个断面为矩形的水槽，问怎样的折法可使水槽的断面面积达到最大 最佳水槽断面问题 （梯形断面） 将问题１推广等腰梯形的水槽，问怎样折法可使水槽断面面积达到最大? 最佳水槽断面问题 （对称五边形断面）将铁板折成如图所示的对称五边形，问怎样的折法可使水槽的断面面积达到最大？ 最佳水槽断面问题（五边） 最佳水槽断面问题（五边） 运行zxy6_6 [s,fval] = fmincon(zxy6_6S,x0,A,b,[],[],lb,ub) 求解 最佳水槽断面问题（多边和无限边） 优化变量数与最大断面面积的关系 断面形状 优化变量数 最大断面积 cm2 矩形断面 1 72 梯形断面 2 83.14 对称五边形 4 88.637 将铁板折成对称7边形、9边形，一般为对称2n+1边形 可以期望最大断面面积得到进一步的增加 随之而来是计算代价也随之增加。 最佳水槽断面问题（无限边） 最佳水槽断面问题（无限边） 最佳水槽断面问题（无限边） 微分法求最大和最小（高等数学） 微分法求最大和最小（高等数学） 微分法求最大和最小（高等数学） 运行zxy6_1.m syms x1 x2 %定义符号变量。 f=x1^3-x2^3+3*x1^2+3*x2^2-9*x1; % 函数z。 v=[x1 x2];df=jacobian(f,v) %计算雅可比。 [X,Y]=solve(df(1),df(2)) % 用指令solve求驻点。 运行zxy6_2画等值线图并将点标注在图上 微分法求最大和最小（高等数学） jacobian(f,v)： 计算函数的符号梯度、雅可比矩阵 如：若f(v),v=[v1 v2] 则df=[df/dv1 df/dv2] 如：若f(v)=[f1(v) f2(v)] ,v=[v1 v2] 则df=[df1/dv1 df1/dv1 df2/dv1 df2/dv2] 微分法求最大和最小（高等数学） solve指令： solve(eqn1,eqn2,…,eqnn) 求n个方程eqn1，eqn2，… ，eqnn所构成的方程组的根（符号解） 盲人下山与迭代寻优 一个盲人处于山上的某一点x0 要下到谷底，他应如何做？ 盲人下山与迭代寻优 Matlab优化工具箱简介 多元函数无约束优化指令fminunc、fminsearch的剖析 Matlab优化工具箱简介 观察： 在命令窗口键入bandemo 选择不同方法观察对香蕉函数的优化结果和过程。 Matlab优化工具箱简介 [x,fval,exitflag,output,grad,hessian]= fminunc(fun,x0,options,P1,…) 其中有些项可以缺省，如exitflag,output,grad,hessian,options,P1,P2,… 等等。 x0是初始点； fun是目标函数，可以用inline指令或建立M文件的方法生成目标函数； Matlab优化工具箱简介 参数的传递： 使用变量P1,P2,…可在目标函数和主程序之间需要传递某些参数 也可使用全局变量Gobal说明来进行传递。 输出： x,fval,exitflag,output,grad,hessian为输出信息， 最优点、最优函数值、算法结束的状态 （exitflag 0 算法收敛；=0达到最大步骤而停止；0算法不收敛）、算法结束后的某些信息（如迭代次数、所使用的优化方法等，可在命令窗口查看output的具体内容）、梯度值和海森矩阵，除x之外均可缺省。 Matlab优化工具箱简介 OPTIONS（控制参数） OPTIONS是一个数组，有多项内容 使用optimset对它进行修改和设置。 optimset(‘属性’，‘属性值’,…) 修改OPTIONS默认值，如默认属性‘LargeScale’＝‘on’， ( 使用“信赖域算法”)。 如果要使用其它方法，就要修改此项设置。 Matlab优化工具箱简介 内联函数inline inline(函数表达式,变量1,变量2, …)，不混淆的情形下变量项可以缺省。 f=inline(100*(x(2)-x(1)^2)^2+(1-x(1))^2), x=[2,2],y=f(x), 备注 Matlab优化工具箱简介 用M文件生成目标函数（套用如下格式） myfun.m func