MATLAB在化工數(shù)據(jù)處理中的應(yīng)用

1、化工數(shù)據(jù)處理與實驗設(shè)計期中論文Matlab在化工數(shù)據(jù)處理中的應(yīng)用姓名:陳奕含學(xué)號:2012115134班級:化學(xué)工程與工藝Matlab在化工數(shù)據(jù)處理中的應(yīng)用一、MATLA崎介Matlab軟件包最初就是1967年由ClereMaler用FORTRAN語言設(shè)計與編寫的。1984年Mathworks公司用C語言完成了Matlab的商業(yè)化版本并推向市場。經(jīng)過20余年的改進,Matlab已發(fā)展成為一個具有極高通用性的、帶有眾多實用工具的運算平臺,成為國際上廣泛認可的優(yōu)秀科學(xué)計算軟件。Matlab就是matrix&laboratory兩個詞的組合,意為矩陣工廠（矩陣實驗室）。就是由美國Mathwo

2、rks公司發(fā)布的主要面對科學(xué)計算、可視化以及交互式程序設(shè)計的高科技計算環(huán)境。它將數(shù)值分析、矩陣計算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模與仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中,為科學(xué)研究、工程設(shè)計以及必須進行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計語言（如C、Fortran）的編輯模式,代表了當(dāng)今國際科學(xué)計算軟件的先進水平。Matlab與Mathematica、Maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計算方面首屈一指。Matlab可以進行矩陣運算、繪制函數(shù)與數(shù)據(jù)、實現(xiàn)算法、創(chuàng)建用戶界面、連接其她編程語言的程序

3、等,主要應(yīng)用于工程計算、控制設(shè)計、信號處理與通訊、圖像處理、信號檢測、金融建模設(shè)計與分析等領(lǐng)域。Matlab的基本數(shù)據(jù)單位就是矩陣,它的指令表達式與數(shù)學(xué)、工程中常用的形式十分相似,故用Matlab來解算問題要比用C,FORTRAN等語言完成相同的事情簡捷得多,并且Matlab也吸收了像Maple等軟件的優(yōu)點,使Matlab成為一個強大的數(shù)學(xué)軟件。在新的版本中也加入了對C,FORTRAN,C+,JAVA的支持。其特點如下:1 .支持多平臺操作系統(tǒng)（Windows、Unix等）。2 .就是一種簡單易學(xué)的編程語言。3 .Matlab程序很容易維護。4 .編程效率很高。由于用戶程序可直接調(diào)用大量的Ma

4、tlab函數(shù),因此編程速度快。5 .用途廣泛?？捎糜跀?shù)值計算與符號計算、數(shù)據(jù)分析、工程與科學(xué)繪圖、圖形用戶界面設(shè)計、建模與仿真、控制系統(tǒng)設(shè)計、數(shù)字圖像信號處理以及財務(wù)工程6 .功能超強。包含600多個常用算法內(nèi)建函數(shù)，有眾多面向具體應(yīng)用的工具箱（如偏微分方程、最優(yōu)化方法、數(shù)理統(tǒng)計、樣條函數(shù)、神經(jīng)網(wǎng)絡(luò)工具箱等）與simulink仿真模塊。此外，其她產(chǎn)品延伸了Matlab的能力，包括數(shù)據(jù)采集與依靠Matlab語言編程產(chǎn)生獨立的C/C+代碼等等。其算法函數(shù)大多由國際知名專家完成，算法穩(wěn)定可靠、效率高。7 .具有開放式結(jié)構(gòu)，擴展功能強。Matlab的開放式結(jié)構(gòu)使Matlab產(chǎn)品族很容易針對特定的需求進

5、行擴充。8 .支持混合編輯技術(shù)。提供與其她面向?qū)ο蟮母呒壵Z言（如VC、VB與CVF等）進行混合編程的接口。9 .Matlab函數(shù)源代碼公開，有助于用戶學(xué)習(xí)與研究算法。10 .第三方公司Matlab軟件產(chǎn)品產(chǎn)品的強力支持，如femlab,可直接求解三維PDE問題。簡言之,Matlab具有非常強大的數(shù)值計算功能、圖形可視化功能與符號運算功能，且簡單易學(xué)，擴展性好,可以與其她面向?qū)ο蟮母呒壵Z言進行混合編程。二、軟件功能主要功能如下表所示：數(shù)值分析數(shù)值與符號計算工程與科學(xué)繪圖控制系統(tǒng)的設(shè)計與仿真數(shù)字圖像處理數(shù)字信號處理通訊系統(tǒng)設(shè)計與仿真財務(wù)與金融工程Matlab的主界面如圖1所示，由菜單欄、工具欄與幾

6、個功能窗口組成。常用的窗口有:命令窗口，用于Matlab命令的輸入與計算結(jié)構(gòu)的顯示，就是最常用的用戶交互窗口；歷史窗口，用于保存用戶已輸入過得命令；當(dāng)前目錄窗口，用于顯示當(dāng)前目錄下的文件；工作窗口，用于顯示與編輯內(nèi)存中的變量?！伴_始”按鈕提供了Matlab的常用功能菜單。圖1Matlab的主界面三、具體舉例-冷卻器的最優(yōu)設(shè)計利用Matlab進行數(shù)據(jù)處理已經(jīng)得到了越來越廣泛的應(yīng)用?；ぶ杏泻芏嗲蠼庾顑?yōu)值的問題，最優(yōu)化方法在化工生產(chǎn)中的應(yīng)用，主要涉及研究與開發(fā)中的實驗方案最優(yōu)化、化工數(shù)學(xué)模型的參數(shù)估計與辨識、化工過程優(yōu)化設(shè)計、工藝操作參數(shù)的優(yōu)化、過程優(yōu)化控制以及最優(yōu)生產(chǎn)調(diào)度等等。本文將舉例冷卻器的

7、最優(yōu)設(shè)計問題。例題:某煉油廠需將煤油從Ti=140C冷去至ijT2=40C,煤油處理量G=3X104kg/h,冷卻介質(zhì)為水，初始溫度為ti=30C,要求設(shè)計一冷卻器(逆流換熱)，并使該冷卻器的年度總費用J盡可能小。已知數(shù)據(jù)：(1)冷卻器每單位面積的總投資費用(材料費與制作費)Ja=200元/m2。(2)冷卻器年折舊率B=15%(包括維修費)。(3)冷卻器總傳熱系數(shù)K=836、8kJ/(m2hC)。(4)冷卻器每年運行時間8=8000ho(5)冷卻水單價Jw=0、04元/噸。(6)冷卻水比熱容cpw=4、184kJ/(kgC)(7)煤油比熱容cpc=2、092kJ/(kgC)。根據(jù)傳熱學(xué)的基本原

8、理分析可知，該冷卻系統(tǒng)存在一個最佳的冷流出口溫度使總費用最低(參見化工原理，冷流出口溫度高，則傳熱推動力小，完成規(guī)定熱負荷所需的傳熱面積打。投資費用就大，而需要的冷卻水量減小，操作費用小)。試用最優(yōu)化方法確定該最優(yōu)冷流體出口溫度，并計算在此最優(yōu)值下冷卻器的最小年費用、冷卻器的傳熱面積、每小時冷卻水用量以及年度投資費與年度操作費各占總費用的百分比。數(shù)學(xué)模型(1)目標函數(shù)J JA AJW1000式中,vA>B、Jw、8為已知參數(shù);A為傳熱面積，m2;w為冷卻水用量，kg/h(2)關(guān)于A與w數(shù)學(xué)模型一一熱平衡方程冷卻器的熱負荷為Q=GcpT1-T2)由熱平衡方程Q=wcpw(t2-t1)=Gc

9、pc(T1-T2)=KAAtm,WQw(3)cpw(t2t1)QKtm式中tm(Ti t2) (T2 t1)lnTH:(5)由(2)與(5)帶入(3)與(4)，然后在帶入(1),則(1)式僅有一個未知變量t20因此該問題屬于單變量最優(yōu)化問題。程序說明采用Nelder-Mead單純形法函數(shù)fminsearch()進行優(yōu)化，其中TotalFee()定義目標函數(shù)式(1)，函數(shù)Area_Water()則根據(jù)式(3)與(4)求A與w。程序清單functionCoolerOptDes%冷卻器的最優(yōu)化設(shè)計(OptimalDesignofaCooler)clearall;clcglobalT1T2Gt1JAb

10、etaKthetaJwCwCcQT1=140;T2=40;G=2e4;t1=30;JA=200;beta=0、15;K=836、8;theta=8000;Jw=0、04;Cw=4、184;Cc=2092;Q=G*Cc*(T1-T2);t0=50;t2=fminsearch(TotalFee,t0);fprintf('優(yōu)化結(jié)果:nn')fprintf('冷卻器最優(yōu)出口溫度為:%、2f%sn',t2,'')allFee=TotalFee(t2); fprintf(' 最小年費用為 :% 、3f 元 n' ,allFee)A w=Are

11、a_Water(t2); fprintf(' 冷卻器傳熱面積為:%3fmA2n',A)fprintf( ' 每小時冷卻水用量為:% 、1fkg/hn' ,w)fee1=JA*A*beta; fee2=Jw*theta*w/1000;fprintf( ' 年度投資費為 :% 、 1f 元,:% 、 2f%sn',fee1,fee1/allFee*100,'%' )fprintf( ' 年度操作費為:% 、 1f 元,:% 、 2f%sn',fee2,fee2/allFee*100,'%' )%func

12、tionJ=TotalFee(t2)globalT1T2Gt1JAbetaKthetaJwCwCcQAw=Area_Water(t2);J=JA*A*beta+Jw*theta*w/1000;%functionAw=Area_Water(t2)globalT1T2Gt1JAbetaKthetaJwCwCcQvar1=T1-t2;var2=T2-t1;dtm=(vaU-var2)/log(van/var2)A=Q/(K*dtm);w=Q/Cw/(t2-t1);Editor-C:U5ersC¥HDocumentsMATlABCooIerOptDes.mFileEditTextGoCell

13、ToolsDebugDesktopIVindowHelp修？©IB日益st"日。Ik西工密瑁i-1.0JTT小口翼I噫喔m1 HfunctionCoolerOptBes¥；令卻罌的最優(yōu)化設(shè)計(OpftinalDesignofaCooler)2 -clearall；clc3 -globalT112GtlJAbeta£thetaJwCvCcQ4 -11=140;T2=4Q;G=2e4;11=305 TA=2i0;beta-0IS;6 -theta=8000;JwQ,口；Cy=4.194;:-Cc=2092;C=C*<c*(Tl-I2);7 -tQNQ;

14、9 -t2=finnEeatch(®TotalFeejtOi;10 -fpnntf(J優(yōu)化結(jié)果；11 -年七Ltrtf冷卻器最優(yōu)出口溫度為士為.二f如n)12 -allFee=IotalFee(t2''fijrintf最小年斐用為:立3f1元F',allFeQ13 -0w=Arta_Wa,t*r(t2);fprintf冷卻器傳熟面積為：S3fn'2njA)It-fprlritf每小時冷卻水用聿為：Ifke/hn")15-feslsjAM*beia；fe*2=JvHheta*w/l000；16-加Hntf年度垠資需為：*If元，占總雷用；%.2

15、f»sW,feel,feel/allFee*100/%J)IT-空口田£仁年度操作題為：,.If元*占總矍用士%.2f%s'rC,fe#2Jf#e2/allF#e*100,)IS9619 functionJ=TotalFee(t2)20 -jlatalTI12GtlJAbetaKthetaJvCwCcQ21 -Av=Ar*a_Water(t2);J=JA*A*beta+Jw*thtta*v/lCOO;22 %23 functionAw=Area_Water(t2)24 -globalTI12GtlJAbelaKthetaJwCwCcQ25 -v-arl=U-t2：v

16、ar2=I2-t1：dtn=tvarl-var2)/log'varl/var2)26 -;v=&/C-r/(t2-t11；27CohfiiriAndWindowONewtoMATLAB?WatchthisVidedsee口erntJEorrejdGettingStarted.dtm=24.0715dtm24.0714優(yōu)化結(jié)果：冷卻器最優(yōu)出口溫度為；92.49VdtJTL=24.0714最小年事用為：L1菸脂二元dtm=24.0714冷卻器傳熱面租為：2Q7715T09115(hw"2每小時冷卻水用量為：1CK3511業(yè)"h年度股資費為二6231452"元，占總?cè)跤枚?4.8騏年度操作費為:5121125.2元，占總?cè)茫?5.11黑»1計算結(jié)果冷卻器最優(yōu)出口溫度為92、49,最小年費用為11352、578元,冷卻器傳熱面積為207、715m2,每小時冷卻水用量為16003、5kg/h。總費用中，年度投資費用為6231、5元,占總費用的54、89%;年度操作費為5121、1元,占總費用的45、11%。四、未來展望Matlab具有非常強大的數(shù)值計算功能、圖形可視化功能與符號運算功