中国传媒大学《MATLAB基础与应用》课件-第11章优化问题的求解.pptVIP

【例11-3】计算下面函数在 上的极值。 解：首先绘制该二元函数在z=0上下两部分的曲面图， 下面编写求极值的主函数，这里由于不知道初始点，所以借助ginput()函数在等高线图形上取点作为初始点的横纵坐标，得到的结果如图： peaks_jizhi.m 【例11-4】Banana函数优化。 解：首先绘制出Banana函数的曲面图。 下面进行优化求解，为获得更好的效果，这里采用动态图形来显示每一步迭代的过程，这里首先介绍一个专门绘制“Banana”函数优化求解迭代点的函数bandemoutfcn()。 调用函数fminsearch()求解上述优化问题，在原绘图程序的基础上输入如下代码： plot3(-1.9,2,267.62,ko,markersize,12,linewidth,1,erasemode,none) text(-1.9,2.3,267.62,\fontname{隶书}\fontsize{16}开始,color,[0,0,0]) plot3(1,1,0,ko,markersize,12,linewidth,1,erasemode,none) text(1,1.3,100,\fontname{隶书}\fontsize{16}结束,color,[0,0,0]) set(gcf,Color,w) x0=[-1.9,2]; f=100*(x(2)-x(1)^2)^2+(1-x(1))^2; % 利用字符串定义函数表达式 OPTIONS=optimset(OutputFcn,@bandemoutfcn); % 设置优化控制参数OutputFcn [x,fval,exitflag,output]=fminsearch(f,x0,OPTIONS) example_11_4_1.m 下面再利用fminunc()函数求解上述优化问题，编写如下程序代码： x0=[-1.9,2]; OPTIONS=optimset(LargeScale,off,OutputFcn,@bandemoutfcn); % 设置优化参数 grad=[100*(4*x(1)^3-4*x(1)*x(2))+2*x(1)-2;100*(2*x(2)-2*x(1)^2)]; % 梯度表达式描述 f=100*(x(2)-x(1)^2)^2+(1-x(1))^2; % 函数表达式 OPTIONS=optimset(OPTIONS,HessUpdate,bfgs,gradobj,on,MaxFunEvals,200, InitialHessType,scaled-identity,LineSearchType,quadcubic); % 更新优化参数 [x,fval,exitflag,output]=fminunc({f,grad},x0,OPTIONS) % 优化问题求解 example_11_4_2.m 【例11-5】求解humps()函数在区间[0,1]内的极小值。 解：为显示该优化问题的每一步的迭代过程，先编写优化方法的输出函数myoutput.m并在优化参数中进行设置。 function stop = myoutput(x,optimvalues,state); % 优化方法的输出函数 history = []; stop = false; if state == iter history=[history;x]; end fplot(@humps,[0,2]) hold on plot(history,humps(history),r.,... MarkerSize,14) % 绘制迭代点 pause(1) % 暂停1秒 再编写如下主程序： options = optimset(OutputFcn, @myoutput); % 设置优化参数OutputFcn [x fval] = fminbnd(@humps,0,1,options) % 单变量优化求解 plot(x,fval,kv,markersize,8,MarkerFaceColor,k) % 绘制最优解 text(x-0.4,fval-7,[极小值点：(,num2str(x),,,num2str(fval),)]) % 添加标注 example_11_5.m 【例11-7】求解下面的多变量约束优化问题。 解：首先编写非线性约束条件的描述函数NonLinear_constr()，再编写如下主程序： x0=[1,1,1];sigma=10^5; lb=[1 0 0];ub=[3,inf,inf]; % 变量上下界 f=@(x,sigma)x(2)*sqrt(16+x(1)^2)+x(3)*sqrt(1+x(1)^2); opti