VC++在煤炭需求总量预测中的应用

第10卷第4期(总第65期)煤矿开采Vol. 10 No.4 ( Series No. 65 )2005年8月oal Mining TechnologyAugust 2005VC++在煤炭需求总量预测中的应用李春睿，鄢丽(辽宁工程技术大学资源与环境工程学院,辽宁阜新123000 )[摘要]以煤炭行业实际数据为基础,结合数学中的一元非线性回归预测理论,建立利用VC++6.0编程工具开发的预测模型程序,利用该模型对我国近年煤炭需求总量作了预测。.[关键词]煤炭经济预测;回归分析法;煤炭需求.[中图分类号]0212.1[文献标识码] A[文章编号] 1006-6225 ( 2005 ) 04-0011-02Application of VC ++ in Coal Total Demand ForecastingLI Chun-rui , YAN Li( College of Resource and Environment Engineering , Liaoning Technical University , Fuxin 123000，China )Abstract : Based on the true data of coal industry ,a predict procedure by VC ++ 6. 0 has been built and studied , combining first-co-der non-linear regress analysis theory at mathematics. The practicability of the procedure in coal economic forecasting is proposed. Andthen a predict of the total need of coal these years in our country is made with the procedure. Finally , necessity of changing the proce-dure with change of coal need is shown.Key words : coal economic forecasting ; approach of regress analysis ; coal demand煤炭经济预测以往采用的方法主要有趋势外推令y'=lny,x'=-,a=Inc,则有y'=a +法、时间序列预测法、灰色理论预测方法等，回归分析是其中的一种。bx'回归分析即假设-个系统的输入和输出之间存( 4)对数函数y=a +bln x变量替换公式:在着某种因果关系,它认为输入量的变化会引起系令x'=Inx ,则有y'=a+bx'统输出量的变化,并且这种变化是符合某种规律(5)双曲线函数一=a+变量替换公式:的。因而用回归分析法建立的模型,不仅能用于经济预测,而且能够解释系统内部运行的原因和各个令y'=+则有y' =a+bx'内部因素之间的关系。下 面仅对回归分析法中的一元非线性回归法用VC ++ 6.0实现进行探讨。(6) S型曲线y=一the1非线性回归模型的建立令y'=,x'=e~* ,则有y' =a+bx'非线性回归在很多情况下,可以通过变量替换-般来说,可利用VC ++6.0强大的绘图功能来转化为线性回归,从而可以为替换后的变量配一先作出x和y的散点图,根据图形的形状大致判断条直线方程,然后再经反变换得到我们所需的曲线用哪种类型曲线最为合适。如果x和y的散点图曲方程。”下面列举了-些最常用的函数及其变量替换线和以上几种曲线图形都类似,应从中选取拟合误公式。(曲线图形略)差最小的曲线进行模拟。为此对每-种曲线类型计( 1 )幂函数y=cx"变量替换公式:算它的剩余标准离差S :令y'=lny ,x'=Inx , a=Inc ,则有y'=a + bx's=-2(yx-站)(2)指数函数y=ce"变量替换公式:其中.γ,是自变量x取x,值时,按相应曲线类型计令y'=Iny ,a=Inc ,则有y' =a +bx'中国煤化工算拟MHCNMHG以合误差。s取最小值(3)数函数y=cex变量替换公式:时所刘应时四线机Er安远的曲线类型。与线性回[收稿日期] 2005 -03-16[作者简介] 李春睿( 1980- ),男，吉林九台人，硕士研究生,联系方式: E - Mail ( li100100100@ tom. com )。11总第65期煤矿开采2005年第4期归一样,S可用来度量回归方程的预测精度。对于2煤炭 需求量预测模型的实现非线性回归,检验所配方程有无实用价值,用量R作为衡量所配曲线拟合原始数据好坏的指标。把1989 ~2004年的煤炭需求量X; ( Mt)以及若R越接近1,所配曲线的拟合效果就越好,工业用煤量Y;(Mt)作为的历史数据作为样本,否则不能选用该非线性回归模型。共16组，如表1。1 - 2(y:-个尸3绘制散 点图和程序设计R=采用非线性回归进行分析,定义2个单精度浮表1煤炭需求量及各 影响因素指标的历史数据MIt年份XY;19891033. 192 809. 12419931209. 195 993. 10219971392. 480 1216. 71020011262. 113 1136. 08019901055. 230810. 909199419981294. 9221149. 52420021366. 055 1271. 9541991 1108. 073892. 1261995 1376. 765 1175. 70719991263. 6531127. 5732001571. 8531480. 78319921159. 463 955. 3841996 1447. 344 1238. 85920001245. 374117. 30020041763. 4631750. 935点型二维数组ary[ ][2], new[ ][2] 及单精de. Move'To(20 , 100) ;度浮点型一维数组P-ER[ ], 其中ary数组存放观de. LineTo(212 , 100) ;测值X;和Y; , new数组存放对X;和Y;进行变换以de. TextOut(212 ,92 ”X;" ,2) ;de. MoveTo(20 , 100) ;后的值,数组P-ER用来存放Y;的历史数据的拟de. LineTo(20，10) ;合。建立相应的数据录入模块,提供全屏幕录入,de. Text0ut(20，10 " Y,",2) ;修改,自动确认和自动清屏的功能。利用VC ++de. MoveTo(20，100) ;6. 0强大的绘图功能作出X;和Y;的散点图。具体for(i=1 ;i<16;i++ )做法:利用MFC AppWizard生成应用程序Graph的de. LineTo(20+ (ary[i][0] -ary [0][0]) *10, 100- (ary[i][1] -ary[0][1 ]));基本框架后,只需对视图类CGraphView加以修改,de. SletObjet( oldpen);即执行View 菜单的ClassWizard 命令,添加On-newpen. DeleteObject( ) ; }Paint的实现代码使程序在视图区绘制出X;和Y;的程序首先对ary数组元素进行变换处理,结果散点图(图1 )。放在new数组中,然后new数组中的new[ ][0]和new[][1]进行-元线性回归的处理,求得回归系数后对系数进行变换以求得非线性回归系数。最后对ary[ ] [0]进行拟合,求得数组P-ER的各元素值和各种参数的值(程序代码略)。4运行结果图1 x;与Y散点表2运行结果与样本对比及误差从图1看出它与指数函数、对数函数和双曲线实际值/10*t预测值/1041误差/%函数类似。用上述3种曲线进行拟合,根据程序计128532.2131515. 82. 32139248. 0140770. 01.09算结果发现用双曲线函数的拟合误差最小,故选择126211. 3126233. 40.02双曲线函数。程序的部分代码为:程序运行后得出实际值X; ,预测值x';,和误void CGraphView : : OnPaint( )差P-ER,如表2所示,曲线类型为:.{inti;1/y= a+ b/xCPen newpen ;系数c=0CPen * oldpen ;中国煤化工CPaintDC de ( this) ;YHCNMHGnewpen. CreatePen( Ps - SOL标准离差S = 0.002632ID,1 , RCB(0,0 ,0)) ;相关系数R =0. 862172oldpen = de. SelectObject( &newpen) ;de. SetMapMode(0) ;(下转51页)12徐坤:围岩稳定性综合 控制技术的应用2005年第4期长锚固,以保证锚固效果,预紧力60~70kN。径由大逐渐变小,前端孔径8mm ,后端孔径4mm。(6)帮部最上面1根锚杆都向上倾斜30° (高(3 )浆液选用化学浆液新型堵水材料马丽散。帮最.上面1根锚杆向上倾斜约50~60°,第2根锚(4)注浆量每孔4~8桶，注浆时间每孔.杆向上倾斜约15 ~25°，低帮最上面1根锚杆可水20min左右, 注浆压力不小于8MPa。平安设)，最下面1根锚杆都向下倾斜30 ~45°安(5)根据具体的渗水淋水位置施工注浆锚杆设,顶板两侧锚杆在垂直顶板方向再向外倾15~再根据堵水效果决定是否加打注浆孔及是否补注,20°，以便于施工和有效控制围岩。通常以1根注浆锚杆控制2m2左右的面积计算注浆(7)帮顶破碎处,施工2~4根右旋全螺纹预锚杆的数量。拉力锚杆加强，锚杆长度为2000mm,直径为(6)注浆作业滞后迎头10 ~20m ,在50m范18mm ,采用全长锚固方式，每根锚杆采用2节围内完成。Z2550型中速树脂药卷锚固。4支护效果. (8)帮顶锚索滞后迎头5 ~ 8m安装齐全。3.3喷混凝土参数为了观测锚杆支护效果,研究支护参数的合理对无水区的顶、帮,为防止煤岩体风化,直接性,巷道每隔30m设置1个测站,对围岩表面位.在围岩表面喷射薄层混凝土,喷层厚30~50mm。移、顶板离层状况、锚杆载荷进行了观测。采用原在顶、帮潮解区喷混凝土，喷层厚70 ~ 100mm。有支护形式的巷道围岩表面位移-般在1m左右,混凝土配比为水泥:黄沙:石子= 1:2:2。喷射混凝顶板离层不断增加,直到巷道修复为止,锚杆也随土作业滞后迎头10~ 20m ,在50m范围内完成。着巷道变形迅速失效。采用综合控制技术施工的3.4煤岩体内注浆 参数7,44风巷围岩表面位移最大值为400mm,顶板离(1)在顶、帮渗流区进行煤岩层内部注浆，层最大值180mm ,锚杆受力稳定在70 ~ 110kN。从孔径42mm,孔深3.5m,高帮/顶板/底帮注浆孔矿压观测结果结合直观支护效果分析，上述围岩稳采用2/2/1型布置方式。定性综合控制技术支护方案在桃园煤矿7,44风巷. (2)注浆锚杆长度3000mm,用12.7mm钢管的应用,基本能够满足安全生产需要,试验研究获制成。钢管底端1.5 ~2. 0m长度内错开钻孔,孔得了成功。[责任编辑:颜立新](上接12页)定了最佳的函数曲线图象,合理的确定了最优的预从上面结果R=0. 976217→1 ,可知预测模型测函数及其参数，因此该方法用于煤炭需求量的可行,利用该模型预测2005~2008年的煤炭需求中、短期预测,具有很高的可信度。量如表3。[参考文献]表3煤炭需求预测结果年份20062008[1]王端武,王浩，张烁,等.我国煤炭需求预测[J].中需求量/104t 167641. 1179465.4 172839. 6161357.2国煤炭, 1999 ,25 (4): 12-15.[2]郭嗣宗,陈刚.信息科学的软计算方法[M].沈阳:东北5结论大学出版社, 2001.[3]楼顺天，于卫，闰华梁. MATLAB程序设计语言[ M].西用VC++6.0建立的煤炭需求预测模型,克服安:西安电子科技大学出版社, 1997.了传统理论推理视图界面的缺点,经过反复实验确[责任编辑:邹正立]xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx(上接36页)(2)回收支柱和顶梁前,必须要先拴牢拉紧5结语后，方可回收,以防回下的铁料窜下伤人。采用以上的管理方法. T314-'工作面安全实现(3)回柱时,如遇到顶板来压,要立即停止初次中国煤化工面收作,采出煤炭20作业,待压力稳定后方可继续回柱工作。余万YHCN M H没有发生大的冒顶、倒(4)回柱不准和打眼、装药等工序平行作业。棚、伤人等事故,取得了较好的经济效益。(5)回柱时，溜子道要填满、闸死、蹬实,[责任编辑:崔德仁]防止从溜子道溅出的煤块伤人。51