《中科大《優(yōu)化設(shè)計》課程大作業(yè)之約束優(yōu)化實驗報告(共8頁)》由會員分享，可在線閱讀，更多相關(guān)《中科大《優(yōu)化設(shè)計》課程大作業(yè)之約束優(yōu)化實驗報告(共8頁)（8頁珍藏版）》請在裝配圖網(wǎng)上搜索。
1、精選優(yōu)質(zhì)文檔-傾情為你奉上約束優(yōu)化設(shè)計實驗報告力學(xué)系型號：聯(lián)想y470CPU：i5-2450M內(nèi)存：2GB系統(tǒng)：win7-64位求解問題：如上是以下三個約束方法共同需要求解的問題，預(yù)估結(jié)果：在（x1,x2,x3）（23,13,12）點(diǎn)附近存在極值。其中，每個方法對應(yīng)的初始條件分別為：（1）隨機(jī)試驗法設(shè)計變量范圍：隨機(jī)試驗點(diǎn)數(shù)：N=1000精度：eps=0.001（2）隨機(jī)方向法初始點(diǎn)：x0=(25,15,5)初始步長：a0=0.5精度：eps=0.001（3）線性規(guī)劃單純形法初始復(fù)合形：X=20 23 25 30;10 13 15 20;10 9 5 0頂點(diǎn)個數(shù)：n=4精度：eps=0.01計
2、算結(jié)果：約束方法所需時間迭代次數(shù)極值點(diǎn)極值隨機(jī)試驗法28.230（22.63,12.63,12.06）3445.50隨機(jī)方向法0.6549（22.67,12.67,11.98）3441.99線性規(guī)劃單純形法6.235（22.59,12.59,12.12）3445.61程序說明：主程序為main，運(yùn)行main后按提示即可得到相應(yīng)約束方法的求解結(jié)果。程序如下：1、 主程序clear;global kk;kk=0;disp(1.隨機(jī)試驗法);disp(2.隨機(jī)方向法);disp(3.線性規(guī)劃單純形法);while 1 n0=input(請輸入上面所想選擇約束優(yōu)化方法的編號（1、2、3）：); if
3、n0=1|n0=2|n0=3 break; end disp(此次輸入無效.);end disp( );disp();xx,yy=fmins(n0);fprintf(迭代次數(shù)為： %8.0fn, kk);disp(所求極值點(diǎn)的坐標(biāo)向量為：);fprintf( %16.5fn, xx);fprintf(所求函數(shù)的極值為： %16.5fn, yy);2、 調(diào)用函數(shù)function xx,yy=fmins(n0)if n0=1 tic;xx,yy=suijishiyan();toc;elseif n0=2 tic;xx,yy=suijifangxiang();toc;elseif n0=3 tic;
4、xx,yy=danchunxing();toc;endend3、 單純形法function xx,yy=danchunxing()clear;global kk;syms a b c;f=-a*b*c;g=-a+2*b+2*c=0;a+2*b+2*c=72;abs(a-b-10)=10;b=20;X=20 23 25 30;10 13 15 20;10 9 5 0;alpha=1.3;sita=0.5;gama=1;beta=0.7;var=a;b;c;eps=0.001;N=size(X);n=N(2);FX=zeros(1,n);while 1 for i=1:n FX(i)=double
5、(subs(f,var,X(:,i); end XS,IX=sort(FX); Xsorted=X(:,IX); px=sum(Xsorted(:,1:(n-1),2)/(n-1); Fpx=double(subs(f,var,px); SumF=0; for i=1:n SumF=SumF+(FX(IX(i)-Fpx)2; end SumF=sqrt(SumF/(n-1); if SumF0 bcon_1=0; else cof_alpha=0.7*(cof_alpha); end end fx2=double(subs(f,var,x2); if fx20 bcon_2=0; if fx3
6、XS(1) count=1; else count=2; end else bcon_2=0; count=3; end end if count=1 Xsorted(:,n)=x3; X=Xsorted; continue else Xsorted(:,n)=x2; X=Xsorted; continue end else if fx2XS(n-1) Xsorted(:,n)=x2; X=Xsorted; continue else if fx2XS(n) Xsorted(:,n)=x2; cof_beta=beta; bcon_3=1; while bcon_30 bcon_3=5; el
7、se cof_beta=cof_beta/2; bcon_3=bcon_3+1; end end if min(gx4)0 fx4=double(subs(f,var,x4); FNnew=double(subs(f,var,Xsorted(:,n); if fx4=0)& (x(1)+2*x(2)+2*x(3)=72)&(abs(x(1)-x(2)-10)=10)&(x(2)=20) f =-x(1)*x(2)*x(3); if f =0)& (x(1)+2*x(2)+2*x(3)=72)&(abs(x(1)-x(2)-10)=10)&(x(2)=20) f =-x(1)*x(2)*x(3)
8、; if f =0)& (x(1)+2*x(2)+2*x(3)=72)&(abs(x(1)-x(2)-10)=10)&(x(2)=20) f =-x(1)*x(2)*x(3); if f fl fl = f; break end end end epsl = abs(f0-fl)/f0); if epsl eps&x(1)=0)& (x10(i)+2*x20(j)+2*x30(k)=72)&(abs(x10(i)-x20(j)-10)=1e-3) z=x10(i)*x20(j)*x30(k); if zmaxz zmax=z; x1=x10(i);x2=x20(j);x3=x30(k); end end end endendxx=x1;x2;x3;yy=zmax;end專心-專注-專業(yè)
中科大《優(yōu)化設(shè)計》課程大作業(yè)之約束優(yōu)化實驗報告(共8頁)
