化工合成和分析2008年第二章Ch2-3.pptVIP

2.2 序贯模块法（接上堂课）; xk+1＝G(xk) 阻尼直接迭代法（加权直接迭代法）: xk+1＝qxk+(1-q)G(xk) q 阻尼因子，一般 0＜q1，可改善收敛的稳定性 q=0 直接迭代 xk+1＝G(xk) q=1 xk+1＝xk q0 外推直接迭代，加速收敛，但稳定性下降 q≥1 无意义 问题：q到底取多少合适? 韦格斯坦法来解决;韦格斯坦法（Wegstein）;k;k;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导出;这就是隐式方程的割线法，与Wegstein法相通的;提示：隐式方程更有普遍性，应用更广。但在流程模拟领域中，物流回路多用显式方程描述的，多用Wegstein法;有界 Wegstein 法;多维Wegstein法;令初始猜值为 xi1, 第二个初值可由直接迭代得到：;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;鄂樟坊三唾猴压砖藩妖陪努熬淬银梧伺正搀矿何灸茧辣胜羊算铀卷津恩错化