《MATLAB圖解精餾塔理論塔板數程序代碼》由會員分享，可在線閱讀，更多相關《MATLAB圖解精餾塔理論塔板數程序代碼(6頁珍藏版)》請在人人文庫網上搜索。

1、MATLAB圖解精餾塔理論塔板數程序代碼function distillation %文件名“distillation”可以更改% 輸入計算所需參數q=1;%輸入進料熱狀況參數R=1.5;%輸入回流比xD=0.95;%輸入塔頂輕組分摩爾分數xW=0.04;%輸入塔底輕組分摩爾分數xF=0.52;%輸入進料輕組分摩爾分數%以下輸入相平衡數據x0=00.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.1;y0=00.0.0.0.0.0.0.0.0.0.。

2、0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.700330.0.0.0.762990.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.1;Yr=(x)R/(R+1).*x+xD/(R+1);%精餾段操作線fun=(x)(q-1)*(R/(R+1).*x+xD/(R+1)-(q*(x-xF)+(q-1)*xF);xQ=fzero(fun,0.5);%求操作點yQ=Yr(xQ);xOP=xW,xQ,xD;yOP=xW,yQ,xD;yfit=linspace(0,1,1001);xfit=interp1(y0,x0,yfit,pchip);%繪制圖形ho。

3、ld onbox onplot(0 1,0 1,k);xlabel(x)ylabel(y)plot(x0,y0,r)plot(xfit,yfit,r-)plot(xF,xF,b*)plot(xQ,yQ,bo)plot(xOP,yOP,b-)k=1;yn(1)=xD;xn(1)=interp1(y0,x0,yn(1),pchip);plot(xD,xn(1),yn(1),yn(1),b-)text(xn(1),yn(1),num2str(1),.HorizontalAlignment,center,VerticalAlignment,bottom)while xn(k)xWyn(k+1)=int。

4、erp1(xOP,yOP,xn(k);k=k+1;xn(k)=interp1(y0,x0,yn(k),pchip);plot(xn(k-1),xn(k-1),yn(k-1),yn(k),b-)plot(xn(k-1),xn(k),yn(k),yn(k),b-)text(xn(k),yn(k),num2str(k),.HorizontalAlignment,center,VerticalAlignment,bottom )endN=k;plot(xn(N),xn(N),yn(N),xn(N),b-)text(xn(N),yn(N),num2str(N),.HorizontalAlignment,。

5、center,VerticalAlignment,bottom )N_Feed=find(xnxF);N_Feed=N_Feed(1);text(0.5,0.5,strcat(所需理論板：,num2str(N),.strcat(進料板位置：,num2str(N_Feed),.HorizontalAlignment,left,VerticalAlignment,top)%以下代碼是為了去掉頂端和右邊坐標軸的刻度box offax2 = axes(Position,get(gca,Position),.XAxisLocation,top,.YAxisLocation,right,.Color,none,.XColor,k,YColor,k);set(ax2,YTick, );set(ax2,XTick, );box on運行結果如圖1所示：圖1 圖解苯-甲苯精餾塔理論塔板數和進料位置。

總結

以上是生活随笔為你收集整理的精馏塔matlab,MATLAB图解精馏塔理论塔板数程序代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。