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

用戶名（公共郵箱mathzyforstudc@）下周開(kāi)始：2：10開(kāi)始上課；5：10結(jié)束最優(yōu)化方計(jì)算方上課記 2009-3-08第1 2009-3-15第2 2009-3-22第3 2009-3-29第4 2009-4-05第5 2009-4-19第6 2009-4-26第7 課后練習(xí) 綜合應(yīng)用 綜合應(yīng)用 綜合應(yīng)用 上課記2009-3-081參 等，最優(yōu)化理論及其方法，電子科技大 概述：最優(yōu)化建模步基本概念：梯度程序設(shè)計(jì)基無(wú)約束優(yōu)化方法一維搜索：黃金分割

fx黃金分割

fxf(x為單峰函數(shù)最速下降法

fx注意：方法的適用范圍；特黃金分割法算法（偽代碼初始化a,b,,la(1)(ba),ra(ba),flf(l),frf

(ba) fl

//定區(qū)間,定al,lr,flfr,ra(ba),frf(r)else//fl

br,rl,frfl,la(1)(ba),flf2009-3-152內(nèi)容安排線性規(guī)劃建模與整數(shù)規(guī)劃建模與無(wú)約束最優(yōu)化方（練習(xí)）某廠生產(chǎn)甲乙兩種口味的飲料,每百箱甲飲料需用原料6千克,工人名,可獲利10萬(wàn)元;每百箱乙飲料需用原料5千克,工人20名,可獲利9萬(wàn)元.今廠共有原料60千克,工人150名,又由于其他條件所限甲飲料產(chǎn)量不超過(guò)8百箱.若投資0.8萬(wàn)元可增加原料1千克,問(wèn)應(yīng)否作這項(xiàng)投資若每百箱甲飲料獲利可增加1萬(wàn)元,問(wèn)應(yīng)否改變生產(chǎn)計(jì)劃. 和Lingo軟件求解.決策:甲乙產(chǎn)量x1x2(箱),整數(shù) 10*s1+9*s2s1=x1/100max=10*s1+9*s2;s1=x1/100;

10*s1+9*s2-s1=x1/100max=10*s1+9*s2-s1=x1/100;隨機(jī)跳躍在可行域內(nèi)產(chǎn)生N個(gè)可行for產(chǎn)生1個(gè)可行解rand用來(lái)產(chǎn)生區(qū)間(0,1)上均勻分布的隨機(jī)數(shù)for語(yǔ)for2009-3-2231部分：最優(yōu)化方法最速下降法注：最速下降法例子推導(dǎo)要修改（符號(hào)問(wèn)題2部分：優(yōu)化建模與求解技術(shù)course_soft.rar（示例程序操作平臺(tái) 命令行輸入mathexpguicmm打開(kāi)程序L（1000位顧客訂單矩陣M(10003局部變量C（1000個(gè)元素）,C（i）第i位顧客已經(jīng)分配得到DVD數(shù)Left(j)表示第jDVD剩余的數(shù)fork=1:10for如果Li行有元素等于k，且現(xiàn)有數(shù)量C(i)<3Left(j)>0，則將第jDVD配給第i顧客：M(i,C(i)+1)=j;2009-3-294第1部分：有約束優(yōu)化方第2部分：建模與求解綜合外點(diǎn)法實(shí)現(xiàn)總結(jié)（測(cè)試函數(shù)test_outer語(yǔ)法方函數(shù)句柄的使用 函數(shù)值計(jì)算，可以調(diào)用feval（第1個(gè)參數(shù)為函數(shù)文件名或函數(shù)句柄，2個(gè)參數(shù)為自變量取值）實(shí)現(xiàn)無(wú)約束規(guī)劃求解：調(diào)用了fminsearch（系統(tǒng)函數(shù)2009-4-055動(dòng)態(tài)規(guī)劃建模階段如何劃狀態(tài)的定最優(yōu)指標(biāo)函數(shù)（包含指標(biāo)函數(shù)，狀態(tài)轉(zhuǎn)移方程明確前向還是后向：決定了狀態(tài)轉(zhuǎn)移方程、最優(yōu)指標(biāo)函數(shù)的形邊界條件（最后考慮2009-4-196隨機(jī)系統(tǒng)模隨機(jī)模擬基模擬（隨機(jī)試驗(yàn)法系統(tǒng)模擬（計(jì)算機(jī)模擬、計(jì)算機(jī)仿真隨機(jī)系統(tǒng)的特點(diǎn)隨機(jī)因素—— 量（服從什么分布—— 量建模設(shè)計(jì)模隨機(jī)的模擬（產(chǎn)生隨機(jī)數(shù)normcdfx2normpdfx22f 2

x2 x2 2 2

x3基本思想：產(chǎn)生若干個(gè)可行點(diǎn)，找出這些可行點(diǎn)中的最小值點(diǎn)fori=1:nifx1^2+x2^2+x3^2>=14x2>=2x1,thenifoptvaluef(x);//保留函數(shù)值optx=x;//當(dāng)前最優(yōu)可行點(diǎn)建模與求解綜合應(yīng)用題：電子管題：修理廠理發(fā)店系統(tǒng)模2009-4-267課后練習(xí)綜合應(yīng)用對(duì)于高速問(wèn)題問(wèn)題，如果約束條件改0x1x2x3x4x530，該如何修改隨機(jī)跳躍法程改進(jìn)DVD租賃問(wèn)題的算法設(shè)計(jì)及其求解程序