matlab在化工研究中的应用

安做瑾犬竽YOf序设计大作业题目的应用学：专：

化学工程学院化学工艺学：_________________学生姓名：阳开设课程：程设计与应用Matlab

在化工研究中的应用1概述化设与分、研与开、工、模预优和控领，常工程模化，用于建的或建过，既可过的、开和可于评估及运转情制改工模关必须建数模，就是质平热量平和量衡为基础并合动学而建立的模型方程。是美国MathWorks司于代中期推出的一套高性的数值计可视化软件。经过MathWorks司不断善，MATLAB已成自控制论、理统计字信号处理间列析态统仿真等领域的最用的和最受欢的软件。MATLAB有强大的扩展功能，使其在很多领域都能充发挥其强大的用。因此人们将称之为继C等高级语后的第四代计算机语言，利用其丰富的数，程人的程中出。MATLAB一性语言程算，数基为，的运和作并具有计、数和字统处理和在计算功具的功，计算结和程可视化具很强开性，不的用学都应的可发工自2070年代，基于和EISPACK软件的分，用开设计以，过30年完善，有数工具展89应用工具和能的应用研，等学领的应用研究等用。

程算学以及算发展，通过建过程理模型通过算机模拟，以便对化工过程进行设计和分析、型预测、优化和控制等。现阶段，用于数值计算的编程语言有VisualBasic6.0、CompacVisualFortan6.5、VisualC++和、MATLAB。尽管Fortran和程序执行效率很，可以分利用现有的代码资源，如NumericalRecipes和库，可提高运行速度，可以生成具有标准用户界面应用软件，数值方面用MATLAB言编程更具优势，Fortran和语言的部数、Recipes和的法数比有，对编能要求高编效较，MATLAB得算法齐，计算功能、图形可视化功能和号运算都十分大也持他向象高语的混编；以，MATLAB作为一方便、快捷的科学和工程计算语言具Fortran和等序无法比拟的优越。国内对的用究，最早为张文乐等，而在化学化领域的应用研起步较晚任佩林等首先用MATLAB编程计算热管固有频率的传递矩阵。本文首讨论其，而，国内化学工领域的应，其语的混合编程、其软的成、工的应用等面等，并以如应用提。2MATLAB在化工研究的应用21Delphi的合编程于MATLAB是一性的语言，运行效率，并且在其，序大都能其运行。利用其强语言的功能，其

的，充等以为化，以语言编的程序为化控制，并过的0LEvariant型MATLAB的工型库，MATLAB和混合编程的目的，实现了炼钢过程氧化期终态的准确，并提供详细源代码及说明。2.2与Basic混合编程付孝锦等引充分利用Basic的面设功能和的大算，通在中用的的方，开发了化学反方程式件。张云鹏等针对数值计算图形处理能力不足，结合的数学函数库，。值计形。3与程、基S上，用C，VB和语言混程，基于最小窗口MPCA非线性建，于，决批次线诊。际于某，早期预报和诊，和性的性。

在工算上应用化工中计问，数对通常数，计算值、、合解常分和方程、线性和非线性方等。解题的法编应的和1111等传统的计算语言具有一定的编程能力，调试过程中花费很多时间和精力。3.1.1线性最小二乘拟合用非线性最小二乘法从实验的二元物系气液平衡实验数据拟合程参数。温度：t=60C,⑴见1。表醇一数据液相组成Xmol

气相组成y

总压056650.8223466.95查得Antonie常数见2。表醇水Antonie常数i1实验实i1实验实验常数

甲醇

水C17.878638.10765C21473.111750.268C3230.0235.0解：由吉布斯自由能和活度系数知识推得，下二元物系实验数据的拟合中，可用式(1)-式⑷中的式⑷作为目标函数。E实验=xI*ln

Y

1实验+x2*lnY实验

⑴E计算=xi*lnY计算+x2*InY

计算

(

2)E实验-GE/(RT))计算=x*[InY实验+ln(xi+

A12X2)]+x2*[验+ln(x2+A⑶

21

*xm{1*[实验1+A12*x2)]+x2In实验2+A21*1)]}j1

⑷在求解Y1实验和Y实验将验据,y1,p的数组当变量待即将[体变量待。由平数公式=YPis/P得丫i=Ki，所以可以算到Y

Y

12121.13151.28071.4744]在进行拟合12和21，只要将标函数式写函数文，再调用这个最小拟合函，即得到结果A解非线性方程(组)

=0.6734

，A=1.0845例：简单蒸馏时，某刻釜残液量与低沸点组分x之间如下的系式：InFo/F=1/(o/x)+0)]对于一甲物系相对挥度a=2.5，开始时物系中含苯%，甲苯40%若蒸馏至原加料量一半，试求液中苯的含量。解：按意将值代入程，得运用函数运算功能将方写为f=并调用SOLVE函数求得此方程的确解为

在化工验据处中应化学程工专实是步解习掌握学程工科实研究方的个要的性节。业验同于础实，的仅证一原、观一现或是寻求个遍用的规律，是为了对性地解一个具有确工景的学程工。工验点，，处理。此算理据的数据处理程简单，数据处理个学程验步，是至要的个实验、分工件是化工实验数据的处理往简单，需通过复杂的数学计算，仅依靠手工计算则需要花费大量的时间，工验的理量大现很高因用算机处验据可大工作程艺一为专科

目有点程之定关，示某化所遵循规律。力、物间等限制，任验所完成验点都有，如何这些有实验点出各数的关，便实验数据的理问题。由化工过程的杂性，实验过中各参数之间关系往是非线的，处理或据拟的工作往，且过程比较琐。若能用算进行据理不仅理果准确很而且还省下多不浪费的人和时，大效率。321处理实验据的程序由于化工实验有很多，而且每一个实验数的处理的步、公式都不一样所很一程序描。是，每个都之处此一程序以描321.1

数据的输入数据的输入是利用提示输入函数put实现的，具体的方是（以输入温度为例—input(‘请入工作温度度；序中用到的组则都以矩的形输入。数理图在研工践常常量得许多散的验数和工数些通要各合方以得连续滑的线。二曲合。4计应用4.1值积分(1)值积基本理

其最用小求积数法多样如简的梯法、普生法牛顿特(Newton-Cotes)都常的。它的基思想是分间成区间i,xi+i],i=1,2,…,n，其中i=b。这样分就为。1

0

变步长辛生法基于变步辛普生法MATLAB出了函来求积分该函的调式为：其fname被积数名和b别是定积的下限和上限。tol来控制积分度，缺省时取。trace控制是否展积分过程，取0展现积分过程，取则现缺时取数积分值n为被积函数的用数。20顿－柯特斯基于牛顿柯特斯法，给出了函数来求定积分。函数的调用格式为：[I,n]=quad8('fname',a,b,tol,trace)其中数的quad函数是的缺省值取该数以地求出定积分的值，且一般情况下函数调用的步数明显小于函数，从而能以更高的效率求出所需的定积分值。30积函数由一个表格义在MATLAB中，对由表格形式义的函数关系的求积分问题用函数。其中向量义函关系Y=f(X)。(2)重定分的数值求解使用供的数就可以直接求出上述二重定积分的数值该函数的调用格式为：该函数求在[a,b]>[c,d]区上的二定积分。参数,trace的法函数quad全相同。4.数值微分在，没有