單變量函數(shù)一維尋優(yōu)搜索方法

1、設(shè)計(jì)任務(wù)：?jiǎn)巫兞亢瘮?shù)一維尋優(yōu)搜索方法第一部分(20%)第二部分(20%)第三部分(30%)第四部分(30%)報(bào)告成績(jī)：指導(dǎo)教師：預(yù)習(xí)報(bào)告(20%)：平時(shí)成績(jī)(10%)：驗(yàn)收答辯成績(jī)(50%)：設(shè)計(jì)報(bào)告成績(jī)(20%)：總評(píng)成績(jī)：源代碼：1. GUI模塊screen=get(0,'Screensize');w=screen(3);h=screen(4);figure('color',1,1,1,'position',0.2*h,0.2*h,0.5*w,0.3*h,'Name','單變量函數(shù)','NumberTi

2、tle','off','MenuBar','none');hcount=uimenu(gcf,'label','&count'); hmenu1=uimenu(hcount,'label','全局搜索','Callback','six'); hmenu2=uimenu(hcount,'label','二分法','Callback','five'); hmenu3=uimen

3、u(hcount,'label','黃金分割','Callback','golden'); hmenu4=uimenu(hcount,'label','FABONONACI','Callback','four'); hplot=uimenu(gcf,'label','&plot'); hmenu1=uimenu(hplot,'label','誤差'); hmenu1=uimenu(hplot,&#

4、39;label','時(shí)間'); uimenu(gcf,'label','&quit','call','close(gcf)');2. 黃金分割函數(shù)模塊tic;a=0;b=1;e=1e-10;a1=b-0.618*(b-a);a2=a+0.618*(b-a);while b-a>e y1=8*a13-2*a12-7*a1+3; y2=8*a23-2*a22-7*a2+3; if y1>y2 a=a1; a1=a2; y1=y2; a2=a+0.618*(b-a); plot(a2,y1,

5、'bh') axis(0.6 0.65 -0.25 0.18); title('黃金分割法'); else b=a2; a2=a1; y2=y1; a1=b-0.618*(b-a); plot(a1,y2,'bh') axis(0.6 0.65 -0.25 0.18); title('黃金分割法'); endendxmin=(a+b)/2ymin=8*xmin3-2*xmin2-7*xmin+3 t=toc; disp('t='); disp(t); per=(xmin-0.6298)/0.62983.fabonn

6、aci法模塊tic;f0=1;f1=1;f=f0+f1;n=2;a=0;b=1;while 1/f>1e-10 f0=f1;f1=f;f=f0+f1; n=n+1;endfun=inline('8*x3-2*x2-7*x+3','x');for k=1:n-2 t1=b+f1/f*(a-b);t2=a+f1/f*(b-a); if fun(t1)<fun(t2) b=t2; else a=t1; end f=f1;f1=f0;f0=f-f1;endif fun(t1)<fun(t2) t=t1,y=fun(t1)else t=t2,y=fun(

7、t2)endplot(t,y,'r*')axis(0.6 0.65 -0.25 0.18);title('fabonnaci法'); t=toc; disp('用時(shí)t='); disp(t); per=(t-0.6298)/0.62984.二分法tic;a=0;b=1;e=1e-10;x0=(a+b)/2;fun=inline('8*x3-2*x2-7*x+3','x');f=8,-2,-7,3;y=polyval(polyder(f),x0);while b-a>e if y>0 b=x0; x0=(

8、a+b)/2; y=polyval(polyder(f),x0); elseif y<0 a=x0; x0=(a+b)/2; y=polyval(polyder(f),x0); else x0=(a+b)/2 endendx0y=fun(x0)plot(x0,y,'g*')axis(0.6 0.65 -0.25 0.18);title('二分法'); t=toc; disp('t='); disp(t); per=(x0-0.6298)/0.62985.全局搜索法模塊tic;h=0.000001; x=0:h:1; f=8*x.3-2*x.2-7*x