化学二元反应蒸塔仿真mathematica实现程序马耀元.docxVIP

(*已知参数*) F = 100; xF = 0.5; xD = 0.95; xW = 0.05; q = 1; han = 6944; af = 2; DD = 50; W = 50; erro = 10^-8; time = 20; cur = 15; fun = NSolve[{yq*(q - 1) == q*xq - xF, yq*((af - 1)*xq + 1) == af*xq}, {xq, yq}]; Rmin = FromDigits[(xD - Evaluate[yq /. fun])/(Evaluate[yq /. fun] - Evaluate[xq /. fun])]; (*************************************控制部分*******************************************) (*可知F=D +W,所以w=0,newton法在此建立，即w-50=0,而w-time=LN-VN=L1+qF-Q/han-time=0=\ R*D+q*F-Q/han-50=0, 即： R*D+q*F-Q/han-W=0 \ ,其中，R,Q的值影响整个反应体系，具体分离目标和它俩的取值有关，利用newton-rapphson迭代法xn+1=xn-f \ (xn)/f(xn)求取合适的