矿井主要通风机性能曲线的最优多项式拟合-西安科技大学学报.PDF

大 第31卷第1期 西安科技 学 学报 V01．3lNo．1 201 1年01月 OFXI’ANUNIVERSITYOF ANDTECHNOLOGYJurL201l JOURNAL SCIENCE 文章编号：1672—9315(2011)Ol一0014—05 矿井主要通风机性能曲线的最优多项式拟合 吴奉亮1’2 摘要：为确定主要通风机性能曲线多项式拟合函数的最优阶数，将回归系数的显著性作为最优 阶数的判定依据。分析了基于最小二乘法的任意阶数主扇性能曲线拟合模型；指出了回归多项 式的各项系数应与风压显著相关才是合理的，在合理基础之上阶数最大的即为最优；采用F检验 法计算了回归多项式各系数的显著性。对于有m个样点的测定数据，通过m—1次拟合，可确定 多项式的最优阶数。编制了相应的程序，实现了最佳阶数回归函数求解与曲线可视化，结合一包 含13个测点的实例显示了这种方法的优越性。 关键词：主扇性能曲线；最优多项式；显著性检验；拟合 722 中图分类号：TD 文献标志码：A 0引言 矿井主扇性能曲线一般是用测定的原始数据和曲线图给出，但在风机选型…、稳定性分析心1及各种 通风网络分析口‘41中需知道其函数表达式，因此常将曲线的测点作为输入数据，由软件求解拟合或插值函 数。对于设计阶段用到的出厂性能参数一般只有曲线图，从图上直接取点采用插值的方法求解，可达到 较高的精度。新风机安装完成后以及生产过程中每隔一定的年限都需进行主扇性能测定，对于这种有测 定原始数据的情况，函数的求解常用的是多项式最dx--乘拟合，其中多项式阶数的选择比较重要，不少学 者对此进行了研究。文献[5]认为拟合值的最大相对误差小于5％时目标函数阶数就是合理的，并用实例 认为5次多项式拟合效果最佳；文献[7]采用延拓逼近算法避免高阶拟合时的振荡现象。已有研究主要 是针对个体主扇性能曲线的拟合分析，得到一固定阶数的拟合方法。对于不同的主扇、或同一个主扇(轴 流式)不同叶片角下曲线的最优拟合阶数不一定都是相同的。本文仍将风机曲线的拟合函数设定为多项 式，但不事先确定阶数，通过对回归结果进行显著性检验来决定合理阶数，结合编制的可视化程序进行拟 合效果分析。 1任意阶数多项式曲线拟合原理 数为k时，函数表达式可记为 I P(q)=bo+∑blq‘． (1) (2)式，式中∑表示执行m个数据点的和。 收稿日期：2010一lO一19 基金项目：陕西省教育厅科学研究计划项目(2010JK683)；西安科技大学博士启动金项目(2010DJ003) 通讯作者：吴奉亮(1977一)，男，四川新都人，博士。讲师，主要从矿井通风与安全方面的教学与研究工作 第1期 吴奉亮：矿井主要通风机性能曲线的最优多项式拟合 15 mb。+(∑q)b。+(∑q2)6：+…+(∑qk)6l=∑p (2) … ； ； i i i 设向量b=(b。，b：，…，b。)，bi的值可以通过求解线性方程组 Ab=夏 (3) 得到，其中系数矩阵A的元素可表示为 口，：∑qt。一挲，(i√：1，2，…，．|})． (4)