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

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元 4.5 军事方案问题 1. 源代码： c=[log(0.9) log(0.92) log(0.8) log(0.84) log(0.85) log(0.88) log(0.75) log(0.8)]; a=[1 0 1 0 1 0 1 0;0 1 0 1 0 1 0 1;537.5 462.5 560 480 605 515 650 550]; b=[48 32 48000]; lb=zeros(8,1); [x,fval]=linprog(c,a,b,[],[],lb); 2. 程序运行： x = 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 48.0000 30.5455 fval = -20.6248 3. 结果分析: 向4号目标投放重型炸弹48枚，轻型炸弹31枚，剩余1枚轻型炸弹 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. 源代码： clear all; close all; a=0.99; t0=100; tf=3; t=t0; value=[3 3 5 3;3 2 5 2;1 5 1 6;4 6 4 10]; amount=size(value,1); num=randperm(amount); for i=1:amount if num(i)==1 pop(i,:)=[1 0 0 0]; elseif num(i)==2 pop(i,:)=[0 1 0 0]; elseif num(i)==3 pop(i,:)=[0 0 1 0]