最優(yōu)化方法課堂記錄

星期二下午5、6在科A229

“離散隨機系統(tǒng)模擬”（2014年下半年開設(shè)——）地點：科A二樓

forifr在區(qū)間[0,0.5]課程“最優(yōu)化方法”：40學(xué)時（10次本年度記 第09周——無約束優(yōu)化方法線性規(guī) 第2周———隨機系統(tǒng)建 第3周———建模案 1997A零件參數(shù)的優(yōu)化設(shè) 1996B節(jié)水洗衣 2005BDVD租 建 寫 編 上屆記 10程序設(shè)計1：For While 1初始化s為空串WhileIfthen ，1 2 混合線性整數(shù)規(guī)劃（第函數(shù)0-1整數(shù)規(guī)劃（提供）—— cj變量最速下降tk用編程實現(xiàn)黃金分割法，并給出算例.輸入：1、一元函數(shù)2、a 3、b Fori=1:x(k+1)=x(k)+tk*ifsymsx1fun=4*(x1-2)^4+2*(x2-3)^2+4*x1*x2;H(1,1)=diff(fun,x1,2);H(1,2)=diff(diff(fun,x1,H(2,1)=diff(diff(fun,x2,1),x1);%H(1,2)熟悉求解函數(shù)：linprog（線性規(guī)劃），bintprog（0-1整數(shù)規(guī)劃 A=[23;14];b=[6;5];L=[0;0];functionf=fun1(x);cikit=M(:,[12]);？？find(t==[2%%n=10;r=fori=1:nr(i)=x/M;n=10000;fori=1:nx(i)f1=%f1=f2=1-functionX=rnd01r=rand;ifr<X=

X=ansColumns1through Columns6through ans ans >>ans= >>>> i=x(i)=x(i-while模擬時間未結(jié)束算法：法模擬產(chǎn)品的指標(biāo)Forc1= 課件中增加Dij的表達式.1.10.21997A1.10.31996B1.10.42005BDVD租 2， C求解優(yōu)化模型程序函 ---xii12, i輪的加水量N總加水量f(x1x2,xN)Excel文件式上行D105D106D107D108D109下行D110D109D108D107D106D2010D105上行D201D202D203D204D205下行D206D205D204D203D202第1行為L+線路編號2行為上行路線:且標(biāo)有“上行:”3行為下行路線:且標(biāo)有“下行 2 While未到文件末尾,工具箱函數(shù)（了解，熟悉一些函數(shù)）Toolbox 求解函數(shù)linprog（linearprogramming）面向?qū)ο缶幊淘O(shè)計模型（design文本文件操作函數(shù)輸入改變---文本文件，二進制文件（語言、引擎，Scriptlanguage）=’abc’height=500color=red height=axis這里假定自變量是x1,x2,…,xN的形式.課后3：面向?qū)ο螅琔ML標(biāo)準(zhǔn)建模語言/工具222 linprogfmincon 編程：語言（if，switch，for，while，常用函數(shù)sparse，spalloc，cell1 2 用法求解，并與lingo的求解結(jié)果對比——看差異。：難點——L形，水塘 索法（或法）來確定。x=x=x=[x,fval]=fmincon(...)[x,fval,exitflag]=fmincon(...)[x,fval,exitflag,output]=fmincon(...)[x,fval,exitflag,output,lambda]=fmincon(...)[x,fval,exitflag,output,lambda,grad]=fmincon(...)

Lingolinprogfminconfminbnd，fminsearchOptimizationToolbox（借助數(shù)學(xué)軟件：工具箱 fors=1:10沒有分配完，thenendfor Mfor rev(cust)<3 andd(id)>0，thenM(cust,rev(cust))=id;endendSyntaxx=bintprog(f)x=bintprog(f,A,b)x=bintprog(f,A,b,Aeq,x=bintprog(f,A,b,Aeq,beq,x=bintprog(f