化工过程实验参数和模型参数拟合.pptVIP

第一章 案例计算示意 Y=a+bx 假设已有5组实验数据（x,y)如上表， 若要计算a和b,只要打开光盘，点击“cip”;找到“各章程序及操作文件”，再打开“第一章及一次拟合文件夹”，点击“DEM一次拟合.vbp”，打开计算程序如下： x 1 2 3 4 5 y 7 8 9 10 12 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句），可改变结果的小数点位数 “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 改为：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+bx^n 只要将原来一次拟合的程序中的一个语句加以修改即可，将原程序中： “x(i) = InputBox(”x(“ i ”)=“)””的后面增加一句： x(i) =(x(i))^n 即可，(x(i) = ax^3);例1-5，P19，具体数据如下表： x -3 -2 -1 2 4 y 14.3 8.3 4.7 8.3 22.7 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) = 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 已知下表实验数据，用二次拟合。 只要打开光盘，点击“cip”;找到“各章程序及操作文件”，再打开“第一章及二次拟合文件夹”，点击“DEM二次拟合.vbp”，打开计算程序如下。 x 1 2 3 4 5 6 y 6 16 32 54 82 116 二次拟合程序的应用 将Open “dem.dat” For Input As #1 改为 Open “ erci1.dat” For Input As #1,并将表上的数据利用上面介绍的记事本保存为erci1.