MATLAB優(yōu)化設(shè)計實驗

1、1會計學(xué)MATLAB優(yōu)化設(shè)計實驗優(yōu)化設(shè)計實驗1.1 MATLAB環(huán)境簡介1.2 數(shù)據(jù)表示1.3 矩陣運算1.4 MATLAB的圖形系統(tǒng)簡介 1.5 源文件（M文件） 啟動MATLAB其窗口如右nCommand History中的命令，即可重復(fù)那條命令。)n注意：沒有邏輯型，但與C語言那樣，非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點擊工具欄上的A(insert text)和右箭頭(insert arrow)按鈕，在圖形中添加輔助說明的箭頭和文字nX=5,6,7； y=8,12,6；z=5,4,3nPie(x); bar(x,z)n使用方法。語句后加分號，則在命令窗口中顯示其計算結(jié)果。1.5.2非函數(shù)文件是無函

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命令，打開編輯窗口可同時打開多個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ī)劃問題某車間生產(chǎn)甲、乙兩種產(chǎn)品。生產(chǎn)甲種產(chǎn)品每件需要材料9Kg、3個工時、4KW電，可獲利60元。生產(chǎn)乙種產(chǎn)品每件需要材料4Kg、10個工時、5KW電，可獲利120元。若每天能供應(yīng)材料360Kg，有300個工時，能供200KW電，問每天生產(chǎn)甲、乙兩種產(chǎ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)的基本格式不能動。因此x和x0應(yīng)為12向量，目標(biāo)函數(shù)可寫成723)(222121xxxxxf2.2無約束非線性及其優(yōu)化函數(shù)無約束非線性及其優(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)化算法信息的一個數(shù)據(jù)結(jié)構(gòu)ngrad 返回目標(biāo)函數(shù)在最優(yōu)解x點的梯度nhessian 返回目標(biāo)函數(shù)在最游解x點Hessian矩陣值nonlcon：函數(shù)由每一個輸入變量x計算出c和ceq，用于估計給定x點的非線性不等式 和非線性等式 。 和 見綜合應(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(