误差理论与数据处理

  • 格式:doc
  • 大小:254.00 KB
  • 文档页数:7

误差理论与数据处理在喷雾雾滴粒径研究中的应用摘要:为提高雾滴粒径测量精确度以及分析雾滴粒径与其表面张力的关系,介绍了误差理论与数据处理的基本理论,并将其用于雾滴粒径测量数据分析和雾滴粒径与其表面张力线性回归分析,结果表明:误差理论与数据处理用于喷雾雾滴粒径研究是可行的,得出的实验结果较为理想。

关键词:误差理论数据处理雾滴粒径动态表面张力线性回归雾化性能是喷雾器的关键指标,研究雾滴对提高喷雾雾化性能有重要意义。

雾滴尺寸影响因素很多,一般通过实验分析,其特点是数据量大,误差存在可能性大。

利用误差理论与数据处理技术分析实验数据,能分析出误差产生原因、尽量减小误差,得到合理的实验结果,而且可以通过回归分析实验参数对雾滴尺寸的影响。

为此,本文介绍了误差分析与数据处理的基本理论并对雾滴粒径测量进行误差分析,对雾滴粒径与其动态表面张力关系进行回归分析。

1 误差分析与数据处理基本理论1.1 误差性质与处理方法(1)随机误差随机误差是指在一定试验条件下,以不可预知的规律变化着的误差,一般具有统计规律,大多服从正态分布,试验次数足够多时,随机误差会减小。

一般通过计算算术平均值、残余误差并对算术平均值进行校核,运用贝塞尔公式、别捷尔斯法、极差法或最大误差法计算测量标准差等来分析随机误差。

(2)系统误差系统误差是指由固定不变的或按确定规律变化的因素所造成的误差,通常由实验装置、环境、方法、人员等引起。

一般通过残余误差观察法、残余误差校核法、不同公式计算标准差比较法、计算数据比较法以及秩和检验法和t检验法等来分析发现系统误差。

减小或消除系统误差要从产生误差根源上消除或者用修正方法消除。

(3)粗大误差粗大误差是一种显然与事实不符的误差,没有一定规律,通常有实验人员粗心大意造成。

测量次数较大时,一般采用3σ准则(莱以特准则)来判别粗大误差,测量次数较少时,可采用罗曼诺夫斯基准则、格罗布斯准则或者狄克松准则来判别粗大误差,其中格罗布斯准则可靠性最高。

1.2 一元线性回归及其方差分析与显著性检验回归分析是处理变量间相关关系的一种数理统计方法。

一元线性回归分析通过试验,分析所得数据,找出两个变量之间关系的经验公式,并且这两个变量呈线性关系。

一元线性回归分析主要运用最小二乘原理计算回归系数。

形如bx b y +=∧0的线性回归方程,回归系数b b ,0的计算方法如下:xxxy l l b =,x b y b -=0式中:∑==Nt tx Nx 11,∑==Nt tyN y 1121121⎪⎭⎫ ⎝⎛-=∑∑==N t t Nt t xx x N x l ,21121⎪⎭⎫ ⎝⎛-=∑∑==N t t Nt t yy y N y lN y x y x l N t t N t t Nt t t xy ⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-=∑∑∑===111回归方程的方差分析是将由自变量取值不同所引起的误差和由其它因素所引起的误差从总变差中分解出来回归平方和: xy Nt t bl y y U =⎪⎭⎫ ⎝⎛-=∑=∧21 ,反映x 和y 线性关系对y 变差的作用;残余平方和: xy yy Nt t t bl l y y Q -=⎪⎭⎫ ⎝⎛-=∑=∧21,反映其它因素对y 变差的作用;总变差: ()yy Nt t l yy S =-=∑=21,Q U S +=回归方程显著性检验通常采用F 检验法,对于一元线性回归 回归平方和自由度:1=U υ;残余平方和自由度:2-=N Q υ 总离差平方和自由度:1-=N S υ F 检验统计量:)2/(1/-=N Q U F查F 分布表中三种不同显著性水平α(0.01、0.05、0.1)值,标记为)2,1(-N F α,将这个值与F 值比较,若)2,1(01.0-≥N F F ,则回归高度显著,若)2,1()2,1(01.005.0-≤≤-N F F N F ,则回归在0.05水平上显著,若)2,1()2,1(05.010.0-≤≤-N F F N F ,则回归在0.1水平上显著,若)2,1(1.0-<N F F ,则认为不显著。

2 喷雾雾滴粒径测量的误差分析实验采用扇形喷头在0.2MPa 的喷雾压力下喷雾,雾滴粒径测量采用济南微纳winner 激光粒度仪。

利用雾滴体积中径50D (以下简称D )来描述雾滴粒径,每组雾滴粒径测量重复8次,实验数据如表1:表1 雾滴尺寸测量结果与分析序号 i D (m μ)i v (m μ)2i v (2m μ)1 174.94 ﹣0.03 0.00092 174.98 ﹢0.01 0.00013 174.96 ﹣0.01 0.00014 175.00 ﹢0.03 0.00095 174.95 ﹣0.02 0.00046 174.97 0.00 0.0000 7 174.99 ﹢0.02 0.0004 8174.96﹣0.010.0001∑=81i iD=1399.75D =174.97∑=81i i v =﹣0.01∑=812i i v =0.0029实验数据分析: (1)求算术平均值D =n D i i ∑=81=1399.75/8≈174.97m μ(2)求残余误差残余误差D D v i i -=,计算值见表1 (3)校核算术平均值及其残余误差根据残余误差代数和校核规则,用第二种规则校核,即:n 为偶数时,A nv i i 281≤∑=,A 为D 末位数的一个单位04.001.04201.081=⨯=<=∑=A nv i i m μ,所以计算正确 (4)判断系统误差根据残余误差校核法,n =8,则42==nK ∑∑==-=∆8441i i i i v v =0.01m μ,差值较小,可认为测量无系统误差(5)求测量列单次测量标准差 根据贝塞尔公式:70029.0112=-=∑=n vni iσ=0.0204m μ (6)判别粗大误差由于测量次数较少,不适合用3σ判别准则,采用格罗布斯判别准则,将雾滴尺寸按大小顺序排列有:94.174)1(=D m μ 00.175)8(=D m μ03.094.17497.174)1(=-=-D D m μ 03.097.17400.175)8(=-=-D D m μ1.470.02040.03)1()8()1(==-==σDD g g查表03.2)05.0,8(0=g ,0)1(g g <,0)8(g g <,因此不含粗大误差 (7)求算术平均值的标准差80204.0==nDσσ≈0.0072m μ(8)求算术平均值的极限误差因测量次数较少,算术平均值的极限误差按t 分布计算71=-=n ν,取α=0.05, 查表得36.2=αt ,则算术平均值的极限误差D lim δ为:017.00072.036.2lim ±=⨯±=±=D t D σδαm μ(9)最后测量结果)017.097.174(lim ±=+=D D D δm μ其它各组喷雾雾滴尺寸数据按同样方法分析。

3 喷雾雾滴粒径与其动态表面张力关系的回归分析实验取20个不同动态表面张力的雾滴,测量其粒径,并分析雾滴粒径与其动态表面张力间的关系,各点雾滴动态表面张力和粒径数据如表2。

由数据可知,雾滴粒径与动态表面张力基本呈线性关系。

表2 实验数据表序号 r (m mN /)D (m )2r2DrD1 20.215 174.97 408.6462 30614.501 3537.0186 2 20.325 175.23 413.1056 30705.553 3561.5498 3 23.256 176.54 540.8415 31166.372 4105.6142 4 27.583 177.28 760.8219 31428.198 4889.9142 5 32.543 178.31 1059.0468 31794.456 5802.7423 6 33.283 177.48 1107.7581 31499.150 5907.06687 34.582 178.29 1195.9147 31787.324 6165.6248 8 37.293 179.23 1390.7678 32123.393 6684.02449 39.396 180.31 1552.0448 32511.696 7103.4928 10 42.256 180.12 1785.5695 32443.214 7611.1507 11 46.264 181.48 2140.3577 32934.990 8395.9907 12 46.583 180.87 2169.9759 32713.957 8425.4672 13 52.535 182.04 2759.9262 33138.562 9563.4714 14 55.795 183.56 3113.0820 33694.274 10241.7302 15 58.275 183.79 3395.9756 33778.764 10710.3623 16 60.076 184.24 3609.1258 33944.378 11068.4022 17 65.013 184.42 4226.6902 34010.736 11989.6975 18 70.026 186.32 4903.6407 34715.142 13047.2443 19 72.544 186.59 5262.6319 34815.828 13535.9850 20 75.785 187.78 5743.3662 35261.328 14230.9073 ∑913.6283618.8547539.2894655081.817166577.45673.1 求一元线性回归方程按回归方程回归系数的计算方法,各数据计算如表3:表3 回归系数计算表∑=Nt tr1=913.628∑=Nt tD1=3618.85N =20r =45.681 D =180.94∑=Nt tr12=47539.2894∑=Nt tD12=655081.817tNt tD r ∑=1=166577.4567N r N t t 21⎪⎭⎫ ⎝⎛∑==45.681N D N t t 21⎪⎭⎫ ⎝⎛∑==180.94N D r N t t N t t ⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛∑∑==11=8328.8728 N r r l N t t Nt t rr 2112⎪⎭⎫ ⎝⎛-=∑∑===5803.4833N D D l N t t Nt t DD2112⎪⎭⎫ ⎝⎛-=∑∑== =278.051N D r D r l N t t N t t Nt t rD⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-=∑∑∑===111 =1263.3223rrrDl l b ==0.2177 r b D b -=0=170.9952 r b b D 0+=∧=170.9952+0.2177r因此得回归方程:∧D =170.9952+0.2177r ,回归曲线如图1图1 雾滴粒径与其动态表面张力线性回归曲线4 小结本文总结了误差理论与数据处理所学基本知识,并将其应用于喷雾雾滴粒径的研究,对雾滴粒径测量进行了误差分析,包括随机误差、系统误差和粗大误差等,给出了雾滴最后测量结果为)017.097.174(lim ±=+=D D D δm μ;运用线性回归分析了雾滴粒径与其表面动态张力之间的关系,得到线性回归方程∧D =170.9952+0.2177r ,并对回归方程做了方差分析与显著性检验。