大連海事大學(xué)matlab平時(shí)作業(yè)范本4

《大連海事大學(xué)matlab平時(shí)作業(yè)范本4》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、43廠(chǎng)址選擇問(wèn)題maxz=75xI2+75x2i+50xI3+50x3I+1OOx23+1OOx32+150y】]+240y12+210y21+120y22+160y31+220y32s.t.3yn+3y12+x12+X】3X21一兀3i<203y21+3y221X[2+X21+X23—兀32—&<3y31+3y32一X13一X23+X9+X32W24yii+y2i+y3i=7人2+『22+『32=13y2i+y?255所有變量均不小于01.源代碼:f=[07550750100501000150240021012001602200]?fA=[011-100-100330000

2、000;0-101010--10000330000;00-100-1110000000330;000000000000110000];b=[20;16;24;5];Aeq=[000000000100100100;000000000010010010];beq=[7;13];lb=zeros(18,1);[x,fval]=linprog(f,Azb,Aeq,beqzlb);2.程序運(yùn)行:00.00000.00001.000000.00000.00000.000007.00000.00000.00000.00008.00000fval=3.4850e+0032.結(jié)杲分析：B地向A

3、地運(yùn)送1萬(wàn)噸原料，ABC三地的建廠(chǎng)規(guī)模分別為7、5、8萬(wàn)噸。最小總費(fèi)用3485萬(wàn)元4.4生產(chǎn)計(jì)劃問(wèn)題max7=4?X

4、+3?x?s.t.2x]+x2<10“Xj+x2<8x2<7xpx2>01.源代碼：c=[-4-3];a=[21;1l;01];b=[1087];lb=zeros(2,1);[x,fval]=1inprog(c,a,b,[],[],lb);2.程序運(yùn)行：2.00006.0000fval=-26.00003.結(jié)果分析：甲機(jī)床2臺(tái)，乙機(jī)床6臺(tái)，最大總利潤(rùn)2萬(wàn)6000元5.4背包問(wèn)題一遺傳算法1.源代碼：1).主程序clearall;closeall;popsize

5、=50;chromlength=6;pm=0?5;a=[3,4,3,4,6,3];c=[12,12,9,16,30,12];b=12;pop=round(rand(popsize,chromlength));whilei<=200fit=fitness(pop,aAb,c);maxi=max(fit);[poplzpopnew]=selection(fit,pop);pop2=pmutation(popl,azb,pm);pop=pop2;i=i+lenddisp('最優(yōu)解為：，);disp(popnew(1,:));disp(*物品總價(jià)值為：J;disp(maxi);2).

6、適應(yīng)度子程序functionfit=fitness(pop,a,b,c);[pxzpy]=size(pop);total=zeros(pxz1);fit=zeros(pxA1);fori=l:pxforj=l:pytotal(i)=total(i)+pop(izj)*a(j);endiftotal(i)>bfit(i,1)=0;elseforj=l:pyfit(i)=fit(i)+pop(i,j)*c(j);endendend3).選擇子程序function[popl,popnew]=selection(fit,pop);[pxzpy]=size(pop);[fitlznum

7、]=sort(fit,1descend1);fori=l:10popnew(i,:)=pop(num(i)A:);endpopl=[popnew;popnew];4).變異子程序function[popl,popnew]=selection(fit,pop);[pxzpy]=size(pop);[fitlznum]=sort(fit,1descend1);fori=l:10popnew(i,:)=pop(num(i)A:);endpopl=[popnew;popnew];1.程序運(yùn)行：最優(yōu)解為：100011物品總價(jià)值等于：545.5指派問(wèn)題1.源代碼：f=[335332521

8、51646410];aeq=[l111000000000000;0000111100000000;0000000011110000;0000000000001111;1000100010001000;0100010001000100;0010001000100010;0001000100010001];beq=[11111111];lb=zeros(16,1);[xffval]=bintprog(f,[]Azaeq,beqzlb);1.程序運(yùn)行：x=0100000100101000fval=106.4基