運(yùn)籌學(xué)實驗報告(編輯修改稿)
2025-08-16 06:26 本頁面
【文章內(nèi)容簡介】 60。vubr=x_f_b(2*nvars+2:3*nvars+1,j)。[xtemp,z,how]=lp(c,[A。c39。],[b。fvub],vlbr,vubr,x0,neqcstr,1)。ftemp=c39。*xtemp。ifstrcmp(how,39。ok39。)tempr0=round(xtemp)。tempr1=floor(xtemp)。tempr2=find(abs(xtemptempr0)pre)。ifisempty(tempr2)xall=[xall,xtemp]。fall=[fall,ftemp]。fvub=min([fvub,fall])。elseifftemp=fvubx_f_b=[x_f_b,[xtemp。ftemp。vlbr。vubr]]。endendend%%%%%%%%%%%%%%%%%%%%%i=i+1。end%thesecondwhilexint=x_f_b(1:nvars,:)。[m,mm]=size(xint)。j=j+1。ifjmmbreakend%theendbecausethebreak（因為中斷而結(jié)束）temp0=round(xint(:,j))。temp1=floor(xint(:,j))。temp2=find(abs(xint(:,j)temp0)pre)。mtemp=length(temp2)。end%theendofwhile（結(jié)束當(dāng)前）else%correspondthesecondif（符合第一個如果）x=xtemp。f=ftemp。end%theendofsecondif（第二個如果的結(jié)束）%%%%%%%%%%%%%%%%%%5if~isempty(fall)fmin=min(fall)。nmin=find(fall==fmin)。x=xall(:,nmin)。f=fmin。endelse%correspondthefirstif（符合第一個如果）x=nan*ones(1,nvars)。end LINDO程序例99頁第6題第二問在命令窗口鍵入以下內(nèi)容:max 11x14x2stx1+2x2=45x1+2x2=162x1x2=4endgin x1gin x2按solve鍵在reports window出現(xiàn): Global optimal solution found. Objective value: Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost X1 X2 Row Slack or Surplus Dual Price 1 2 3 4 非線性規(guī)劃程序名：unpfun1函數(shù)unpfun1函數(shù)執(zhí)行實例：（課本第四章152頁16題） min4x1+6x22x1x22x^2,取初始點(diǎn)x0=（1；1）在命令窗口鍵入以下內(nèi)容: f = 4*x(1)2*x(1)*x(2) x(2)^26*x(2)。 [x,fval]=fminunc(@unpfun1,x0)（調(diào)用無約束線性規(guī)劃函數(shù)）按運(yùn)行按鈕在solution report 窗口得到以下結(jié)果:Warning: Gradient must be provided for trustregion method。 using linesearch method instead. In fminunc at 265Optimization terminated: relative infinitynorm of gradient less than .x = * fval =
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1