化工合成与分析年Ch.pptVIP

function out=weg1dk=1; e1=0.0001; e2=0.0001;x(1)=100;x(2)=x(1)^(1/3)+6;for i=1:100, k=k+1; G(k-1)=x(k-1)^(1/3)+6; G(k)=x(k)^(1/3)+6; S(k)=(G(k)-G(k-1))/(x(k)-x(k-1)); q(k)=S(k)/(S(k)-1); x(k+1)=q(k)*x(k)+(1-q(k))*G(k); if abs(x(k+1)-x(k))=e1, if abs((x(k+1)-x(k))/x(k+1))=e2, break; end endendk[S q G]x 直接迭代与韦格斯坦 Wegstein 令猜值为X11＝3X21＝10X31＝6 令猜值为X11＝3X21＝10X31＝6 function out=weg3d clc; k=1; e1=0.0001; e2=0.0001; x(1,1)=3; x(2,1)=10; x(3,1)=6; x(1,2)=(4-x(3,1)^0.5)^3/x(2,1); x(2,2)=(81-x(1,1)^2-x(3,1)^2)^0.5; x(3,2)=(33-x(1,1)^0.5)/x(2,1); for i=1:100, k=k+1; G(1,k-1)=(4-x(3,k-1)^0.5)^3/x(2,k-1); G(2,k-1)=(81-x(1,k-1)^2-x(3,k-1)^2)^0.5; G(3,k-1)=(33-x(1,k-1)^0.5)/x(2,k-1); G(1,k)=(4-x(3,k)^0.5)^3/x(2,k); G(2,k)=(81-x(1,k)^2-x(3,k)^2)^0.5; G(3,k)=(33-x(1,k)^0.5)/x(2,k); S(1,k)=(G(1,k)-G(1,k-1))/(x(1,k)-x(1,k-1)); S(2,k)=(G(2,k)-G(2,k-1))/(x(2,k)-x(2,k-1)); S(3,k)=(G(3,k)-G(3,k-1))/(x(3,k)-x(3,k-1)); q(1,k)=S(1,k)/(S(1,k)-1); q(2,k)=S(2,k)/(S(2,k)-1); q(3,k)=S(3,k)/(S(3,k)-1); x(1,k+1)=q(1,k)*x(1,k)+(1-q(1,k))*G(1,k); x(2,k+1)=q(2,k)*x(2,k)+(1-q(2,k))*G(2,k); x(3,k+1)=q(3,k)*x(3,k)+(1-q(3,k))*G(3,k); if ((x(1,k)-x(1,k-1))^2+(x(2,k)-x(2,k-1))^2+(x(3,k)-x(3,k-1))^2)^0.5=e1, if ((x(1,k)-x(1,k-1))^2+(x(2,k)-x(2,k-1))^2+(x(3,k)-x(3,k-1))^2)^0.5/(x(1,k)^2+x(2,k)^2+x(3,k)^2)^0.5=e2, break; end end end k G S q x 三、过程系统参数优化过程系统模型与最优化模型联解得到一组使工况目标函数最佳的决策变量（优化变量），从而实现最佳工况 面向方程法总结 原理：将描述整个过程系统的数学方程式联立求解，从而得出模拟计算结果 优点：解算快速有效，对设计、优化问题灵活方便，效率较高 缺点：形成通用软件比较困难；不能利用现有大量丰富的单元模块；缺乏实际流程的直观联系；计算失败之后难于诊断错误所在；对初值要求比较苛刻；计算技术难度较大 过程系统稳态模拟三种方法的比较 10月17日作业 2.6（1）用Wegstein法作上节题2.4(1)，迭代至y4 ，看是否收敛？ （2）与直接迭代法比较，Weigstein迭代对收敛有什么好处？ 2.7 用Wegstein法求解右侧方程组，令初值为x11=2，x21=2，计算到 x3，写出过程。 2.8上机练习，13周作：分别用直接迭代法、Wegstein法求解右侧方程组：（令初值为x11=3，x21=3, x31=3，计算到 x3） 推荐软件Matlab，最多迭代100步，程序和x中间结果用WORD编排并打印好后，各班收齐后统一上交。可以用课件中的程序修改成自己的程序，绝对误差 、相对误差均不大于0