《环境系统分析第六次作业》.docxVIP

环境系统分析 上机实验报告 上机主题 优化工具箱的使用（II） 有条件极值的求解 上机时间 班 级 学 号 姓 名 评 分 市政环境学院 环境科学与工程系 指导教师 一.上机实验目的 学会使用Matlab优化工具箱求解规划问题。 学会将待求函数线性化后再进行求解。 二.上机实验题目 已知费用函数形式为：Ci -200Qi°78 1000Qi°.78 i2.5 河流水质约束为：L2三6 试用Matlab优化工具箱求解 题目5所给数据的排放口最优规划问题 （提示：1.不考虑溶解氧的约束； L U」（L；-浦），式中L为各河段在满足水质约束情况下的 BOD排 放浓度，则i =土上，如：1=!5土互 Li0 150 可以直接将目标函数代入优化工具箱进行计算，也可以将目标函数分段 线性化之后，再代入优化工具箱进行计算（但要注意目标函数中的决策变量要 和约束方程中的变量匹配）；最好是两种方法都试一下，比较一下结果有何不 同？ .上机程序要点 0 1 「1.6378] 0 ..1.1753 和m= 0 0.8359 0.05263- 05867- 0.04854 0 0 0.03483 0.03125 0 U = 0.02477 0.02223 0.0400 。01738 0.01560 0.02807 给定的河流水质约束为L2 6 费用函数的形式为：Ci =200Qi°78 ?1000Q：78 i2.5 ,,T T T ,U L ￡(L°-m) 0.0485 3.5903 J0.03483L, +0.03125L2 3.9883 0.02477Li +0.02223L2 +0.0400La 4.2805 O.01738Li +0.01560L2 P.02807l_3 扣.05263L4 4.4950 分别将Q = 0.5 分别将Q = 0.5、 Q =0.3、Q=0.4、Q =0.5 代入可得 C1 C1 =116.5 582.4 茨 C2 =78.2 391 广 C3 = 97.9 489.3 ；5 C4 =116.5 582.4 ；5 函数直接代入优化工具箱 其程序为 Li[nit=[6 6 6 6]P ： 、约束条件 Lnex = Ln^(U)*(Liwit-in) 嫌i 入水质最大BOP 值 Ltnirc=[0 0 0 0]\ LawFLimit-sn; Lc^fmneim (0fuix Lib魏』U Ld叽[],[L Lmitx Lit削〕 、最就处理EOD n-l~Lc. /L; q=[0a5 0.3 0.4 0. 5]\ a=20Oq. 0, 78; b=1000*qt D. 78 h c-a+br *Tir 2, 5; CES^(e) 、满是约束荥件最低费用 得到的解为 LC = 94.5599 54.2078 413222 33.9916 C = 1.3353e+Q03 即排入口排入的 BOD含量分别为94.5599、54.2078、43.3222、 33.9916 mg / Lmg 此时有最费用1335.3万元。 ②将函数线性化后再代入优化工具箱 分三段对费用函数线性化，效率区间为0W%壬0.3 0.3壬听i2《0.85 0.85 4^3刍 程序如下： Limit=[6 6 6 ; Ltuax - iny(U)+(Lintit-ra) 冬输入水辰最大：量 LonFLimitTn; e=[0 0.3 0. 85 1]; a=[116.5 78.2 97. 9 口上曰; b=[5BX4 391 4S9. 3 582.41*; 5=甚丁05(土 4)； for rFl ：I ：4 k(n, 1)=0； 姑0 l)=b(n)A S*(e(2)*4.5-e(l)*4.5); Bw(tbi)=Q, CsG^ 】)=0; s(il 2)=3*As(p,l)/(e(2)-e(l))*3; Fl far mi=2：l ：3 As.tn, m)=b(n)/4J 5+^e(nH-l) 4. 5^e(in) 4. S); w)=_b (n)+c (m) /3- S+fcfnH-I j 3p 5-t(ni) 3- 5); k (nP re)=k，n, m-1 )+s (n, m)* ? (id) -? (m-1)); Cs (n,m)=-(€ (nH-1 -e .mJ *2/2*k(nJmJ; s(n, nrt-1 =3*(As(n, m)+Bs (.n^ nl+CsGi, n))/(e .nrt-l 3; end end %求解5 计算各段线性函数的斜率如下： 1.0e+003 * 0 0.0638 0.S642 1.0417 0 0.04 网 0.3788 1.2364 0 0.0536 0.4740 1.5473 0 0.06