大连海事大学作业汇编-matlab平日作业3.docVIP

4.3厂址选择问题 源代码： f=[0 75 50 75 0 100 50 100 0 150 240 0 210 120 0 160 220 0]; A=[0 1 1 -1 0 0 -1 0 0 3 3 0 0 0 0 0 0 0;0 -1 0 1 0 1 0 -1 0 0 0 0 3 3 0 0 0 0;0 0 -1 0 0 -1 1 1 0 0 0 0 0 0 0 3 3 0;0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0]; b=[20;16;24;5]; Aeq=[0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 1 0 0;0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 1 0]; beq=[7;13]; lb=zeros(18,1); [x,fval]=linprog(f,A,b,Aeq,beq,lb); 程序运行: x = 0 0.0000 0.0000 1.0000 0 0.0000 0.0000 0.0000 0 7.0000 0.0000 0 0.0000 5.0000 0 0.0000 8.0000 0 fval = 3.4850e+003 3. 结果分析： B地向A地运送1万吨原料，ABC三地的建厂规模分别为7、5、8万吨。最小总费用3485万元 4.4 生产计划问题 1. 源代码： c=[-4 -3]; a=[2 1;1 1;0 1]; b=[10 8 7]; lb=zeros(2,1); [x,fval]=linprog(c,a,b,[],[],lb); 2. 程序运行: x = 2.0000 6.0000 fval = -26.0000 3. 结果分析： 甲机床2台，乙机床6台，最大总利润2万6000元 5.4背包问题_遗传算法 1. 源代码： 1). 主程序 clear all; close all; popsize=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)); while i=200 fit=fitness(pop,a,b,c); maxi=max(fit); [pop1,popnew]=selection(fit,pop); pop2=pmutation(pop1,a,b,pm); pop=pop2; i=i+1 end disp(最优解为：); disp(popnew(1,:)); disp(物品总价值为：); disp(maxi); 2). 适应度子程序 function fit=fitness(pop,a,b,c); [px,py]=size(pop); total=zeros(px,1); fit=zeros(px,1); for i=1:px for j=1:py total(i)=total(i)+pop(i,j)*a(j); end if total(i)b fit(i,1)=0; else for j=1:py fit(i)=fit(i)+pop(i,j)*c(j); end end end 3). 选择子程序 function [pop1,popnew]=selection(fit,pop); [px,py]=size(pop); [fit1,num]=sort(fit,descend); for i=1:10 popnew(i,:)=pop(num(i),:); end pop1=[popnew;popnew]; 4). 变异子程序 function [pop1,popnew]=selection(fit,pop); [px,py]=size(pop); [fit1,num]=sort(fit,descend); for i=1:10 popnew(i,:)=pop(num(i),:); end pop1=[popnew;popnew]; 2. 程序运行： 最优解为： 1 0 0 0 1 1 物品总价值等于： 54 5.5 指派问题 1. 源代码： f=[3 3 5 3 3 2 5 2 1 5 1 6 4 6 4 10]; aeq=[1 1