《化工过程实验参数及模型参数拟合》由会员分享，可在线阅读，更多相关《化工过程实验参数及模型参数拟合（31页珍藏版）》请在金锄头文库上搜索。
1、第一章,案例计算示意,Y=a+bx,假设已有5组实验数据（x,y)如上表， 若要计算a和b,只要打开光盘，点击“cip”;找到“各章程序及操作文件”，再打开“第一章及一次拟合文件夹”，点击“DEM一次拟合.vbp”，打开计算程序如下：,Y=a+bx,运行该程序，出现界面如下： 将实验次数改为“5”，点击“开始计算” 出现以下界面：,Y=a+bx,根据提示，分别输入5组实验数据，可得结果，如下图：,Y=a+bx,根据计算结果可知，a=5.6,b=1.2,若x=1代入，则得y=6.8,有0.2的偏差。整体的平均绝对偏差为0.24。注意可以通过改变程序中以下语句（共有4句），可改变结果的小数点位数
2、“A00 = Int(a0 * 100000 + 0.5) / 100000” 原设置为“5”位小数点，若改为“6”，只需将该语句改为： “A00 = Int(a0 * 1000000 + 0.5) / 1000000”即可，其他3句也仿照。,Y=a+bx,若要计算P4中的例1-1，由于数据太多，如果仍像刚才一样输入，就可能输错数据，也无法判断数据输入是否正确，可通过数据文件的形式输入。具体方法如下： 先将原程序中:For i = 1 To n x(i) = InputBox(“x(“ & i & “)=“) y(i) = InputBox(“y(“ & i & “)=“) Next i 改为
3、：Open “ shujv1.dat“ For Input As #1 For i = 1 To n Input #1, ax, ay /需在前面增加：Dim ax,ay x(i) = ax y(i) = ay Next i Close #1,Y=a+bx,数据文件“ shujv1.dat” 可用“记事本”建立。,Y=a+bx,文件名为“ shujv1.dat”，保存类型为“所有文件”；可直接放在“vb”程序的目录下，这样可直接引用。,Y=a+bx,计算结果,文件存放形式,Y=a+bxn,只要将原来一次拟合的程序中的一个语句加以修改即可，将原程序中： “x(i) = InputBox(”x(“
4、 例1-5，P19，具体数据如下表：,Y=aebx,将拟合公式两边取对数，可得线性形式： lnY=lna+bx，程序修改如下： 原y(i) = InputBox(“y(” & i & “)=”)后面，增加一句 ： y(i) = log(y(i) 原a0 = (m * d - c * p) / (n * d - c 2) ,增加一句 ： a0=exp(a0),即可。例1-6，P20,Y=aXb,将拟合公式两边取对数，可得线性形式： lnY=lna+blnX，程序修改如下： 原X(i) = InputBox(“X(” & i & “)=”)后面，增加一句 ： X(i) = log(X(i) y(i
5、) = InputBox(“y(” & i & “)=”)后面，增加一句 ： y(i) = log(y(i) 原a0 = (m * d - c * p) / (n * d - c 2) ,增加一句 ： a0=exp(a0),即可。P28,第8题。,Y=aXb,For i = 1 To n x(i) = Exp(x(i) y(i) = Exp(y(i) eer = eer + (A00 * x(i) A11 - y(i) 2 sd = sd + Abs(A00 * x(i) A11 - y(i) Next i,二次拟合程序的应用,y=a0+a1x+a2x2 已知下表实验数据，用二次拟合。 只要打
6、开光盘，点击“cip”;找到“各章程序及操作文件”，再打开“第一章及二次拟合文件夹”，点击“DEM二次拟合.vbp”，打开计算程序如下。,二次拟合程序的应用,将Open “dem.dat” For Input As #1 改为 Open “ erci1.dat” For Input As #1,并将表上的数据利用上面介绍的记事本保存为erci1.dat，运算可得结果。,二次拟合程序的应用,点击”开始计算”，输入”6”，回车，得下图结果，有关小数点的问题同一次拟合。,二次拟合程序的应用,利用二次拟合程序，通过适当的变换，可以拟合许多类型的实验数据，在程序中的修改工作量也非常小，可以非常方便使用，
7、当然目前也有许多软件可以胜任该工作，但作者觉得，与其化很大的精力去找到拟合函数，还不如自己构建一个函数来得快，下面我们通过几个案例来说明。,y=a0+a1xn+a2xm,上面这个拟合公式，只要实验数据在不同的次方下有意义，就可以非常方便地利用光盘上的程序。如拟合成：(erci2.dat) y=a0+a1x 0.5+a2x1.8 数据为： 只要“x1(i) = xx x2(i) = xx 2,x1(i) = xx0.5 x2(i) = xx 1.8,y=a0+a1x 0.5+a2x1.8,Nu=c1Rec2Prc3,只要将上式两边进行对数运算，就可以得到线性表达式：lnNu=lnc1+c2lnR
8、e+c3lnPr 以P11例1-3的数据为例，建立erci4.dat数据文件，程序的修改工作如下：（dim xx1) Input #1, xx, YY x1(i) = xx x2(i) = xx 2 y(i) = YY,Input #1, xx, xx1, YY x1(i) = Log(xx) x2(i) = Log(xx1) y(i) = Log(YY),Nu=c1Rec2Prc3,C1=exp(-3.77255)=0.023; C2=0.80003 C3=0.30003 注意必须对常数项数据进行变换处理。,y=a0+a1x1 +a2x2,已知如下表实验数据：,Open “erci5.dat
9、“ For Input As #1 For i = 1 To m Input #1, xx1, xx2, YY x1(i) = xx1 x2(i) = xx2 y(i) = YY Next i,建立 erci5.dat 数据文件，修改部分见右边，计算结果如下：,单变量任意次拟合,拟合公式： 计算公式,单变量任意次拟合,已知下列数据，请分别用2次，3次，4次，5次拟合。(假设已建立renci1.dat） 调用光盘第一章的”任意次拟合中DEMN.frm”,将”m=7”,改为”m=6”, 将Open “g:dem.dat” For Input As #1 改为Open “renci1.dat“ Fo
10、r Input As #1,数据按下式构建：1+2*x+3*x2+0.01x3+0.0001*x4+0.0002*x5,单变量任意次拟合,点击运行，分别输入2，输入3，输入4，输入5，可以得到左图的计算结果，当输入5时，得到的结果和构建的公式相同。,1+2*x+3*x2+0.01x3+0.0001*x4+0.0002*x5,单变量任意次拟合,该程序没有太多的变型，只能将x变成sinx或cosx或logx等简单的变化，在化工应用中，某些特殊的场合，强制要求A0=0，则可调用强制过原点任意次拟合程序，具体操作和“任意次拟合”相仿。,N维变量拟合,拟合公式： 计算公式：,N维变量拟合,人为构建3变量数据，公式如下：y=1+2x1+x2+x3 调用光盘第一章的” N维变量拟合多变量拟合.vbp”,N维变量拟合,如对N维变量拟合的程序作适当的修改，可以产生许多变型，希望同学们自己建立各种变型。如,
