大连海事大学Matlab平时作业范本6.docxVIP

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.5455fval = -20.62483. 结果分析:向4号目标投放重型炸弹48枚，轻型炸弹31枚，剩余1枚轻型炸弹4.6检验员聘用问题源代码:f=[40 36]; a=[-5 -3;1 0;0 1];b=[-45 9 15]; lb=zeros(2,1); [x,fval]=linprog(f,a,b,[ ],[ ],lb); 程序运行:x = 9.0000 0.0000fval = 3605.4 背包问题_模拟退火1. 模拟退火源代码：clearall; closeall;k=[12;12;9;16;30;12]; k=-k; d=[3;4;3;4;6;3]; restriction=12;num=6;sol_new=ones(1,num); E_current=inf;E_best=inf;sol_current=sol_new;sol_best=sol_new;a=0.95;t0=97;tf=3;t=t0;p=1;while t=tffor r=1:1000tmp=ceil(rand.*num);sol_new(1,tmp)=~sol_new(1,tmp); while 1 q=(sol_new*d=restriction); if ~q p=~p; tmp=find(sol_new==1); if p sol_new(1,tmp)=0;elsesol_new(1,tmp(end))=0;endelsebreakendendE_new=sol_new*k; ifE_newE_currentE_current=E_new;sol_current=sol_new;ifE_newE_best′E_best=E_new;sol_best=sol_new;endelseif randexp(-(E_new-E_current)./t) E_current=E_new;sol_current=sol_new;elsesol_new=sol_current;endendend t=t.*a; enddisp(最优解为：)sol_bestdisp(物品总价值等于：)val=-E_best;disp(val)disp(背包中物体总重量是：)disp(sol_best*d)2．程序运行：最优解为：sol_best = 1 0 0 0 1 1物品总价值等于： 54背包中物品总重量是：125.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:amountifnum(i)==1pop(i,:)=[1 0 0 0];elseifnum(i)==2pop(i,:)=[0 1 0 0];elseifnum(i)==3pop(i,:)=[0 0 1 0];elseifnum(i)==4pop(i,:)=[0 0 0 1];endendE_best=inf;pop_best=pop;while t=tf ind1=0;ind2=0;while(ind1==ind2) ind1=ceil(rand.*amount); ind2=ceil(rand.*amount);end tmp1=num(ind1);num(ind1)=num(ind2);num(ind2)=tmp1;for i=1:amountifnum(i)==1