化工过程实验参数及模型参数拟合_精品文档PPT资料.ppt
- 文档编号:15488645
- 上传时间:2022-11-02
- 格式:PPT
- 页数:31
- 大小:1.65MB
化工过程实验参数及模型参数拟合_精品文档PPT资料.ppt
《化工过程实验参数及模型参数拟合_精品文档PPT资料.ppt》由会员分享,可在线阅读,更多相关《化工过程实验参数及模型参数拟合_精品文档PPT资料.ppt(31页珍藏版)》请在冰豆网上搜索。
,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,由于数据太多,如果仍像刚才一样输入,就可能输错数据,也无法判断数据输入是否正确,可通过数据文件的形式输入。
具体方法如下:
先将原程序中:
Fori=1Tonx(i)=InputBox(x(&
i&
)=)y(i)=InputBox(y(&
)=)Nexti改为:
Open“shujv1.datForInputAs#1Fori=1TonInput#1,ax,ay/需在前面增加:
Dimax,ayx(i)=axy(i)=ayNextiClose#1,Y=a+bx,数据文件“shujv1.dat”可用“记事本”建立。
Y=a+bx,文件名为“shujv1.dat”,保存类型为“所有文件”;
可直接放在“vb”程序的目录下,这样可直接引用。
Y=a+bx,计算结果,文件存放形式,Y=a+bxn,只要将原来一次拟合的程序中的一个语句加以修改即可,将原程序中:
“x(i)=InputBox(”x(“例1-5,P19,具体数据如下表:
Y=aebx,将拟合公式两边取对数,可得线性形式:
lnY=lna+bx,程序修改如下:
原y(i)=InputBox(“y(”&
“)=”)后面,增加一句:
y(i)=log(y(i)原a0=(m*d-c*p)/(n*d-c2),增加一句:
a0=exp(a0),即可。
例1-6,P20,Y=aXb,将拟合公式两边取对数,可得线性形式:
lnY=lna+blnX,程序修改如下:
原X(i)=InputBox(“X(”&
X(i)=log(X(i)y(i)=InputBox(“y(”&
P28,第8题。
Y=aXb,Fori=1Tonx(i)=Exp(x(i)y(i)=Exp(y(i)eer=eer+(A00*x(i)A11-y(i)2sd=sd+Abs(A00*x(i)A11-y(i)Nexti,二次拟合程序的应用,y=a0+a1x+a2x2已知下表实验数据,用二次拟合。
只要打开光盘,点击“cip”;
找到“各章程序及操作文件”,再打开“第一章及二次拟合文件夹”,点击“DEM二次拟合.vbp”,打开计算程序如下。
二次拟合程序的应用,将Open“dem.dat”ForInputAs#1改为Open“erci1.dat”ForInputAs#1,并将表上的数据利用上面介绍的记事本保存为erci1.dat,运算可得结果。
二次拟合程序的应用,点击”开始计算”,输入”6”,回车,得下图结果,有关小数点的问题同一次拟合。
二次拟合程序的应用,利用二次拟合程序,通过适当的变换,可以拟合许多类型的实验数据,在程序中的修改工作量也非常小,可以非常方便使用,当然目前也有许多软件可以胜任该工作,但作者觉得,与其化很大的精力去找到拟合函数,还不如自己构建一个函数来得快,下面我们通过几个案例来说明。
y=a0+a1xn+a2xm,上面这个拟合公式,只要实验数据在不同的次方下有意义,就可以非常方便地利用光盘上的程序。
如拟合成:
(erci2.dat)y=a0+a1x0.5+a2x1.8数据为:
只要“x1(i)=xxx2(i)=xx2,x1(i)=xx0.5x2(i)=xx1.8,y=a0+a1x0.5+a2x1.8,Nu=c1Rec2Prc3,只要将上式两边进行对数运算,就可以得到线性表达式:
lnNu=lnc1+c2lnRe+c3lnPr以P11例1-3的数据为例,建立erci4.dat数据文件,程序的修改工作如下:
(dimxx1)Input#1,xx,YYx1(i)=xxx2(i)=xx2y(i)=YY,Input#1,xx,xx1,YYx1(i)=Log(xx)x2(i)=Log(xx1)y(i)=Log(YY),Nu=c1Rec2Prc3,C1=exp(-3.77255)=0.023;
C2=0.80003C3=0.30003注意必须对常数项数据进行变换处理。
y=a0+a1x1+a2x2,已知如下表实验数据:
Openerci5.datForInputAs#1Fori=1TomInput#1,xx1,xx2,YYx1(i)=xx1x2(i)=xx2y(i)=YYNexti,建立erci5.dat数据文件,修改部分见右边,计算结果如下:
单变量任意次拟合,拟合公式:
计算公式,单变量任意次拟合,已知下列数据,请分别用2次,3次,4次,5次拟合。
(假设已建立renci1.dat)调用光盘第一章的”任意次拟合中DEMN.frm”,将”m=7”,改为”m=6”,将Open“g:
dem.dat”ForInputAs#1改为Open“renci1.datForInputAs#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维变量拟合的程序作适当的修改,可以产生许多变型,希望同学们自己建立各种变型。
如,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 化工 过程 实验 参数 模型 拟合 精品 文档