(9)--实验七 最优化方法数学实验.pptVIP

*（4）运算结果为：x=0.36600.36600.2679fval=0.5829（3）主程序为:x0=[0;0;0];A=[];b=[];Aeq=[111];beq=[1];vlb=[];vub=[];[x,fval]=fmincon(fun5,x0,A,b,Aeq,beq,vlb,vub,myfun5)↙*例13资金使用问题：设有400万元资金,要求4年内使用完,若在一年内使用资金x万元,则可得效益当年不用的资金可存入银行,年利率为10%.试制定出资金的使用计划,以使4年效益之和为最大.万元(效益不能再使用),解设变量表示第i年所使用的资金数,则有*（1）先建立M文件fun6.m,定义目标函数:functionf=fun6(x)f=-(sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4)));（2）再建立M文件mycon1.m定义非线性约束：function[g,ceq]=mycon1(x)g(1)=x(1)-400;g(2)=1.1*x(1)+x(2)-440;g(3)=1.21*x(1)+1.1*x(2)+x(3)-484;g(4)=1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4;ceq=0;*（3）主程序为:x0=[1;1;1;1];vlb=[0;0;0;0];vub=[];A=[];b=[];Aeq=[];beq=[];[x,fval]=fmincon(fun6,x0,A,b,Aeq,beq,vlb,vub,mycon1)↙得到*五、应用举例两料场分别向各工地运送多少吨水泥，使总的吨千米数最小．，某公司有6个建筑工地，每个工地的位置（用平面坐标a，b表示，单位：km）及水泥日用量d（单位：t）由表给出．目前有两个临时料场位于，，日储量各有20t．假设从料场到工地之间均有直线道路相连，试制订每天的供应计划，即从工地1234561.258.750.55.7537.251.250.754.7556.57.753547611为了进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量仍各为20吨，问应建在何处，节省的吨千米数有多大．*，因此目标函数为：解向工地从料场的运送量为记工地的位置为，水泥日用量，，料场位置为，日储量为约束条件*，这是一个线性规划问题，（1）当用临时料场时，决策变量为计算程序如下：a=[1.258.750.55.7537.25];b=[1.250.754.7556.57.75];d=[3547611];e=[2020];g=[51];h=[27];c=zeros(1,12);fori=1:6c(i)=sqrt((g(1)-a(i))^2+(g(2)-b(i))^2);c(i+6)=sqrt((h(1)-a(i))^2+(h(2)-b(i))^2);endA=[ones(1,6),zeros(1,6);zeros(1,6),ones(1,6)];Aeq=[eye(6),eye(6)];v1=zeros(1,12);[x,fval]=linprog(c,A,e,Aeq,d,v1,[])↙*x=3.00005.00000.00007.00000.00001.00000.00000.00004.00000.00006.000010.0000fval=136.2275即由料场A、B向6个工地运料方案如表所示：123456料场A350701料场B0040610总的吨千米数为136.2275．计算结果为：*（2）当为新建料场选址时，决策变量为和，这是一个非线性规划问题．计算程序如下：首先建立目标函数的Ｍ文件：functionf=liaoch(x)a=[1.258.750.55.7537.25];b=[1.250.754.7556.57.75];f1=0;fori=1:6s(i)=sqrt((x(13)-a(i))^2+(x(14)-b(i))^2);f1=s(i)*x(i)+f1;endf2=0;fori=7:1