环境系统分析第六次作业

1、共享知识分享快乐环境系统分析上机实验报告上机主题优化工具箱的使用（ii）有条件极值的求解上机时间班级学号姓名评分市政环境学院卑微如蝼蚁、坚强似大象共享知识分享快乐环境科学与工程系指导教师一上机实验目的1.学会使用matlab优化工具箱求解规划问题。2.学会将待求函数线性化后再进行求解。二上机实验题目已知费用函数形式为：c=200q0.78+1000q0.78hiiii2.5vvvvv河流水质约束为：l062试用matlab优化工具箱求解题目5所给数据的排放口最优规划问题（提示：1.不考虑溶解氧的约束;2.lu-1(l0-m),式中l为各河段在满足水质约束情况下的bod排2l放浓度，则h=li0

2、-li，如：hii01=150-li；1503可以直接将目标函数代入优化工具箱进行计算，也可以将目标函数分段线性化之后，再代入优化工具箱进行计算（但要注意目标函数中的决策变量要和约束方程中的变量匹配）；最好是两种方法都试一下，比较一下结果有何不同？卑微如蝼蚁、坚强似大象共享知识分享快乐0.017380.015600.028070.052630.5867三上机程序要点0000.048541.63780.034830.03125001.1753u=和m=0.024770.022230.040000.8359给定的河流水质约束为l062ul(l-m)费用函数的形式为：c=200q0.78+1000q

3、0.78hiii0i2.50.04854l0.03483l1+0.03125l120.02477l1+0.02223l20.01738l1+0.01560l2+0.0400l3+0.02807l3+0.05263l43.59033.98834.28054.4950分别将q=0.5、q=0.3、q=0.4、q=0.5代入可得c=116.5+582.4h2.5c=78.2+391h2.51122c=97.9+489.3h2.5c=116.5+582.4h2.53344函数直接代入优化工具箱其程序为卑微如蝼蚁、坚强似大象共享知识分享快乐得到的解为即排入口排入的bod含量分别为94.5599、54.2

4、078、43.3222、33.9916mg/lmg此时有最费用1335.3万元。卑微如蝼蚁、坚强似大象共享知识分享快乐将函数线性化后再代入优化工具箱分三段对费用函数线性化，效率区间为0h0.30.3h0.850.85h1i1i2i3程序如下：计算各段线性函数的斜率如下：即为下表：0h0.30.3h0.850.85h1i1i2i3卑微如蝼蚁、坚强似大象共享知识分享快乐c=116.5+582.4h211c=78.2+391h222c=97.9+489.3h233c=116.5+582.4h244s=63.811s=42.821s=53.631s=63.841s=564.212s=378.822s=

5、47432s=564.242s=1841.713s=1236.423s=1547.333s=1841.743minz=cul(l-m)代入目标函数可得4ii=1=409.1+0.319l+2.821l+9.2085l111213+0.214l+1.894l+6.182l212223+0.268l+2.37l+7.7365l313233+0.319l+2.821l+9.2085l4142430s.t0.0485l0.0348l1+0.0312l120.0248l1+0.0222l20.0174l1+0.01560l2+0.0400l3+0.0281l3+0.0526l44.59034.98835

6、.28055.4950s.t即0.0106l0.0070l1+0.0063l120.0047l1+0.0042l20.0032l1+0.0028l2+0.0076l3+0.0051l3+0.0096l41111又因为l=200-l=200-(l+l+l)iiri1i2i3s.t1111l+l+ll11+l12+l131213l+l+l1213l11+l12+l13105.4401+0.897(l+l+l)236.2161212223+0.897(l+l+l)+1.6144(l+l+l)489.1625212223313233+0.897(l+l+l)+1.6144(l+l+l)+3.0268(l+l+l)991.6286212223313233414243将以上条件带入优化工具箱程序为卑微如蝼蚁、坚强似大象共享知识分享快乐最后可得解即最低费用为1270.7万元。显然可得，将其线性化后再带入优化工具箱，其最优费用明显下降。四上机实验心得卑微如蝼蚁、坚强似大象共享知识分享快乐优化工具箱使用得当可以简便地算出有条件极值的问题，若与lingo进行对照，lingo编写程序更为简捷，但是局限性较大，只能用于求解优化问题，而matlab求解问题面较广。但在求解优化问题方面，个人感觉lingo使用起来比matlab更简单。同时matlab在help函数的时候，显示的全是英文，对于