MATLAB優(yōu)化設(shè)計(jì)實(shí)驗(yàn)

1、1會(huì)計(jì)學(xué)MATLAB優(yōu)化設(shè)計(jì)實(shí)驗(yàn)優(yōu)化設(shè)計(jì)實(shí)驗(yàn)1.1 MATLAB環(huán)境簡(jiǎn)介1.2 數(shù)據(jù)表示1.3 矩陣運(yùn)算1.4 MATLAB的圖形系統(tǒng)簡(jiǎn)介 1.5 源文件（M文件） 啟動(dòng)MATLAB其窗口如右nCommand History中的命令，即可重復(fù)那條命令。)n注意：沒(méi)有邏輯型，但與C語(yǔ)言那樣，非0為真，0為假n例如：nx=1 2 3 4 5 6 7 8 9 0n y=1,2,3,4,5,6,7,8,9,0n z=1,2,3,4,5;2,3,4,5,6,7a的值na nnmmnnmmnnfhgxba,111kccc,.,21n（a與b的維數(shù)必須相同）nnbababa1112211nnababab11

2、122110009871000654000321Ann所以 AB相當(dāng)于A1B9 , 8 , 7 , 6 , 5 , 4 , 3 , 2 , 1 a 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 b12047150108310076521321321xxxxxxxxbAx nLine(x,y)n點(diǎn)擊工具欄上的A(insert text)和右箭頭(insert arrow)按鈕，在圖形中添加輔助說(shuō)明的箭頭和文字nX=5,6,7； y=8,12,6；z=5,4,3nPie(x); bar(x,z)n使用方法。語(yǔ)句后加分號(hào)，則在命令窗口中顯示其計(jì)算結(jié)果。1.5.2非函數(shù)文件是無(wú)函

3、數(shù)頭的M文件，由若干命令和注釋構(gòu)成。相當(dāng)于批處理文件或主程序如： %Filename is a sine.m x = 0:0.1:2*pi; y = sin(x); plot(x,y) %可包含漢字注釋1、建立新M-文件 在命令窗中選File/New/M-file命令，打開(kāi)編輯窗口可同時(shí)打開(kāi)多個(gè)M文件2、保存M文件在文本編輯器中選File/Save As 命令min. .tsxcT0,xbAx,)(,),.,(,)(32, 1nmArankxxxxxaATnnmij0,),.,(,),.,(2121bbbbbccccTmTn且n x,fval,exitflag,output,lambda =

4、linprog() 321645)(xxxxf20321xxx42423321xxx302321 xx0, 0, 0321xxxn -78.0000練習(xí)解下面線性規(guī)劃問(wèn)題某車間生產(chǎn)甲、乙兩種產(chǎn)品。生產(chǎn)甲種產(chǎn)品每件需要材料9Kg、3個(gè)工時(shí)、4KW電，可獲利60元。生產(chǎn)乙種產(chǎn)品每件需要材料4Kg、10個(gè)工時(shí)、5KW電，可獲利120元。若每天能供應(yīng)材料360Kg，有300個(gè)工時(shí)，能供200KW電，問(wèn)每天生產(chǎn)甲、乙兩種產(chǎn)品各多少件，才能獲得最大的利潤(rùn)。2.1線性規(guī)劃及其優(yōu)化函數(shù)線性規(guī)劃及其優(yōu)化函數(shù)提示：)(min)(maxxfxfn fval = 2.00003)sin()(xxf練習(xí) 的最小值提示：

5、 x,fval = fminsearch(fun,x0)的基本格式不能動(dòng)。因此x和x0應(yīng)為12向量，目標(biāo)函數(shù)可寫(xiě)成723)(222121xxxxxf2.2無(wú)約束非線性及其優(yōu)化函數(shù)無(wú)約束非線性及其優(yōu)化函數(shù)cxbAxxxfTT21)(或)2(),1 (xxfy fminunc()nxHessian矩陣值22212123)(xxxxxf應(yīng)用背景給定區(qū)間x1xx2，求函數(shù)f(x)的最小值。X可以是多元向量nx = 3nfval = -11)3()(2 xxfnnx,fval,exitflag = fmincon()nx,fval,exitflag,output = fmincon()noutput 返

6、回優(yōu)化算法信息的一個(gè)數(shù)據(jù)結(jié)構(gòu)ngrad 返回目標(biāo)函數(shù)在最優(yōu)解x點(diǎn)的梯度nhessian 返回目標(biāo)函數(shù)在最游解x點(diǎn)Hessian矩陣值nonlcon：函數(shù)由每一個(gè)輸入變量x計(jì)算出c和ceq，用于估計(jì)給定x點(diǎn)的非線性不等式 和非線性等式 。 和 見(jiàn)綜合應(yīng)用例。應(yīng)用背景約束條件如下，求解函數(shù)f（x）的最小值ubxlbbeqxAeqbxA*其中A、Aeq是矩陣，b、x、beq、lb、ub均是向量0)(xc0)(xceq)(xc)(xceq321)(xxxxfTx10,10,100 27220321xxx022321xxx2722321xxxbxA*nLargeScale 是采用large-scale算法還是medium-scale算法nnoptions=optimset(optimfun)noptions=optimset(oldoptions, param1,value1,.)n options=optimset(Display, iter, TolFun,1e-8) 12424(*)(22122211xxxxxexfx)(minxfx5 . 12121xxxx1021xx01x02x05 . 12121xxxx01021xxnC=1.5+x(1)*x(2)-x(1)-x(2);-x(