x=[];y=[];
while button==1
[xi,yi,button]=ginput(1);
%ginput(n)用鼠标选n个点,xi,yi分别为点的横坐标和纵坐标.
plot([xi],[yi],'r.','MarkerSize',10),hold on, %画所选的点.
[zi,dz1,dz2]=zxy6_3f(xi,yi); %计算函数值和梯度方向.
v=zi;
contour(X1,X2,Z,[v v],'-'), %在点所在的高度画一条等高线.
axis([xmin xmax ymin ymax]),
x=[x,xi];y=[y,yi];
H_line2=plot(x,y); %画已走的路径连线.
set(H_line2,'color','red','linewidth',2); %设置颜色和线宽.
xt=xi-dz1;yt=yi-dz2;
H_line=plot([xi xt],[yi yt],'k:','linewidth',1); %画最速下降方向路径.
end %若按左键button=1,继续循环.若按右键,button~=1,循环终止 . 】
zxy6_3f.m(模拟山谷的二次函数程序)
【 function [f,df1,df2]=zxy6_4f(x1,x2)
f=8*x1.*x1+9*x2.*x2-10*x1.*x2-12*x1-6*x2; %计算函数值.
if nargout > 1
df1=2*8*x1-10*x2-12*ones(size(x1)); %计算梯度向量.
df2=2*9*x2-10*x1-6*ones(size(x2));
end 】
6.3 .实验与观察(Ⅱ):Matlab优化工具箱简介
6.3.1多元函数无约束优化指令fminunc和fminsearch
1. 观察:运行香蕉函数的优化程序bandemo.m

2. 使用fminunc和fminsearch指令

◆ 观察:用inline生成函数.
【 f=inline('100*(x(2)-x(1)^2)^2+(1-x(1))^2'),
x=[2,2],y=f(x), %代入一个点计算看看效果 . 】
3. bandemo.m的简化和剖析
zxy6_4.m
【 clf; clear
%以下程序段是画香蕉函数图形.
xx = [-2:0.125:2]'; yy = [-1:0.125:3]'; [x,y]=meshgrid(xx',yy') ;
meshd = 100.*(y-x.*x).^2 + (1-x).^2; conts = exp(3:20);
xlabel('x1'),ylabel('x2'),title('Minimization of the Banana function')
contour(xx,yy,meshd,conts), hold on
plot(-1.9,2,'ro'),text(-1.9,2,'Start Point')
plot(1,1,'ro'),text(1,1,'Solution')

上一页下一页