现代设计方法上机.docVIP

现代设计方法上机 现代设计方法上机 PAGE 现代设计方法上机 现代设计方法上机 一、黄金有哪些信誉好的足球投注网站法 1.程序框图 2.参考程序： function [xstar,minvalue,a,b]=golden(f,a,b,e) % golden 为黄金有哪些信誉好的足球投注网站方法，本程序只针对一元函数 % [ xstar,fminvalue]=golden(f,a,b,eps)为输入函数 f为函数 eps为允许误差 (a,b)为初始点，xstar为极小值点，fminvalue为极值; % 昆明理工大学机电工程学院机械工程系 刘孝保 2012年11月 a1=a+*(b-a); a2=a+*(b-a); y1=subs(f,x,a1); y2=subs(f,x,a2); flag=true; k=0; disp([黄金有哪些信誉好的足球投注网站过程如下：]); while(flag) disp([k= num2str(k) ;a= num2str(a) ;a1= num2str(a1) ;a2= num2str(a2) ;b= num2str(b) ;y1= num2str(subs(f,x,a1)) ;y2= num2str(subs(f,x,a2)) ]); if(y1y2) b=a2; a2=a1; y2=y1; a1=a+*(b-a); y1=subs(f,x,a1); else a=a1; a1=a2; y1=y2; a2=a+*(b-a); y2=subs(f,x,a2); end length=b-a; if(lengthe) flag=false; end k=k+1; end xstar=(b+a)/2; minvalue=subs(f,x,xstar); disp([******极小值点为：x= mat2str(xstar) ；******极小值为fmin= num2str(minvalue) 。******]); 3.自己尝试利用求f=x^2+2x在区间【-3,5】的极小值。 二、最速下降法 1.最速下降法程序框图 2.参考程序 function [xstar,fminvalue]=fsxsteep(f,e,a,b) % fsxsteep 为最速下降法，本程序只针对二元函数 % [ xstar,fminvalue]=fsxsteep(f,e,a,b)为输入函数 f为函数 e为允许误差 (a,b)为初始点，xstar为绩效指点，fminvalue为极值; % 昆明理工大学机电工程学院机设系 刘孝保 disp([ ================================================================]); disp([ ******************2.无约束优化之最速下降法************************]); disp([ ================================================================]); x1=a;x2=b; xx0=[x1,x2]; k=0; fx1=diff(f,x1); %对x1求偏导数 fx2=diff(f,x2); %对x2求偏导数 g=[fx1 fx2]; %梯度 g1=subs(g); %把符号变量转为数值 d=-g1; fig=figure(); plot(x1,x2,o)%画初始点 text(x1,x2,[ ( num2str(x1) , num2str(x2) )])%画初始点坐标 hold on; axis equal ff=f-subs(f); hhh=ezplot(ff);%画目标函数等值线 set(hhh,linewidth,3); hold on; pause(1); alpha=YJFC(f,x1,x2,d);%求第一步最优步长 disp([初始迭代方向为：( mat2str(d) )。]); disp([初始迭代的最优步长为alpha= num2str(alpha) 。]); while (abs(norm(g1))=e) disp([*******************第 num2str(k) 次迭代*********************]); xx1=subs(xx0+alpha*d); x1=xx1(1); x2=xx1(2); xx1=subs(xx1); plot(