第10卷第4期(總第65期)煤礦開(kāi)采Vol. 10 No.4 ( Series No. 65 )2005年8月oal Mining TechnologyAugust 2005VC++在煤炭需求總量預(yù)測(cè)中的應(yīng)用李春睿，鄢麗(遼寧工程技術(shù)大學(xué)資源與環(huán)境工程學(xué)院,遼寧阜新123000 )[摘要]以煤炭行業(yè)實(shí)際數(shù)據(jù)為基礎(chǔ),結(jié)合數(shù)學(xué)中的一元非線性回歸預(yù)測(cè)理論,建立利用VC++6.0編程工具開(kāi)發(fā)的預(yù)測(cè)模型程序,利用該模型對(duì)我國(guó)近年煤炭需求總量作了預(yù)測(cè)。.[關(guān)鍵詞]煤炭經(jīng)濟(jì)預(yù)測(cè);回歸分析法;煤炭需求.[中圖分類號(hào)]0212.1[文獻(xiàn)標(biāo)識(shí)碼] A[文章編號(hào)] 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煤炭經(jīng)濟(jì)預(yù)測(cè)以往采用的方法主要有趨勢(shì)外推令y'=lny,x'=-,a=Inc,則有y'=a +法、時(shí)間序列預(yù)測(cè)法、灰色理論預(yù)測(cè)方法等，回歸分析是其中的一種。bx'回歸分析即假設(shè)-個(gè)系統(tǒng)的輸入和輸出之間存( 4)對(duì)數(shù)函數(shù)y=a +bln x變量替換公式:在著某種因果關(guān)系,它認(rèn)為輸入量的變化會(huì)引起系令x'=Inx ,則有y'=a+bx'統(tǒng)輸出量的變化,并且這種變化是符合某種規(guī)律(5)雙曲線函數(shù)一=a+變量替換公式:的。因而用回歸分析法建立的模型,不僅能用于經(jīng)濟(jì)預(yù)測(cè),而且能夠解釋系統(tǒng)內(nèi)部運(yùn)行的原因和各個(gè)令y'=+則有y' =a+bx'內(nèi)部因素之間的關(guān)系。下 面僅對(duì)回歸分析法中的一元非線性回歸法用VC ++ 6.0實(shí)現(xiàn)進(jìn)行探討。(6) S型曲線y=一the1非線性回歸模型的建立令y'=,x'=e~* ,則有y' =a+bx'非線性回歸在很多情況下,可以通過(guò)變量替換-般來(lái)說(shuō),可利用VC ++6.0強(qiáng)大的繪圖功能來(lái)轉(zhuǎn)化為線性回歸,從而可以為替換后的變量配一先作出x和y的散點(diǎn)圖,根據(jù)圖形的形狀大致判斷條直線方程,然后再經(jīng)反變換得到我們所需的曲線用哪種類型曲線最為合適。如果x和y的散點(diǎn)圖曲方程?！毕旅媪信e了-些最常用的函數(shù)及其變量替換線和以上幾種曲線圖形都類似,應(yīng)從中選取擬合誤公式。(曲線圖形略)差最小的曲線進(jìn)行模擬。為此對(duì)每-種曲線類型計(jì)( 1 )冪函數(shù)y=cx"變量替換公式:算它的剩余標(biāo)準(zhǔn)離差S :令y'=lny ,x'=Inx , a=Inc ,則有y'=a + bx's=-2(yx-站)(2)指數(shù)函數(shù)y=ce"變量替換公式:其中.γ,是自變量x取x,值時(shí),按相應(yīng)曲線類型計(jì)令y'=Iny ,a=Inc ,則有y' =a +bx'中國(guó)煤化工算擬MHCNMHG以合誤差。s取最小值(3)數(shù)函數(shù)y=cex變量替換公式:時(shí)所劉應(yīng)時(shí)四線機(jī)Er安遠(yuǎn)的曲線類型。與線性回[收稿日期] 2005 -03-16[作者簡(jiǎn)介] 李春睿( 1980- ),男，吉林九臺(tái)人，碩士研究生,聯(lián)系方式: E - Mail ( li100100100@ tom. com )。11總第65期煤礦開(kāi)采2005年第4期歸一樣,S可用來(lái)度量回歸方程的預(yù)測(cè)精度。對(duì)于2煤炭 需求量預(yù)測(cè)模型的實(shí)現(xiàn)非線性回歸,檢驗(yàn)所配方程有無(wú)實(shí)用價(jià)值,用量R作為衡量所配曲線擬合原始數(shù)據(jù)好壞的指標(biāo)。把1989 ~2004年的煤炭需求量X; ( Mt)以及若R越接近1,所配曲線的擬合效果就越好,工業(yè)用煤量Y;(Mt)作為的歷史數(shù)據(jù)作為樣本,否則不能選用該非線性回歸模型。共16組，如表1。1 - 2(y:-個(gè)尸3繪制散 點(diǎn)圖和程序設(shè)計(jì)R=采用非線性回歸進(jìn)行分析,定義2個(gè)單精度浮表1煤炭需求量及各 影響因素指標(biāo)的歷史數(shù)據(jù)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點(diǎn)型二維數(shù)組ary[ ][2], new[ ][2] 及單精de. Move'To(20 , 100) ;度浮點(diǎn)型一維數(shù)組P-ER[ ], 其中ary數(shù)組存放觀de. LineTo(212 , 100) ;測(cè)值X;和Y; , new數(shù)組存放對(duì)X;和Y;進(jìn)行變換以de. TextOut(212 ,92 ”X;" ,2) ;de. MoveTo(20 , 100) ;后的值,數(shù)組P-ER用來(lái)存放Y;的歷史數(shù)據(jù)的擬de. LineTo(20，10) ;合。建立相應(yīng)的數(shù)據(jù)錄入模塊,提供全屏幕錄入,de. Text0ut(20，10 " Y,",2) ;修改,自動(dòng)確認(rèn)和自動(dòng)清屏的功能。利用VC ++de. MoveTo(20，100) ;6. 0強(qiáng)大的繪圖功能作出X;和Y;的散點(diǎn)圖。具體for(i=1 ;i<16;i++ )做法:利用MFC AppWizard生成應(yīng)用程序Graph的de. LineTo(20+ (ary[i][0] -ary [0][0]) *10, 100- (ary[i][1] -ary[0][1 ]));基本框架后,只需對(duì)視圖類CGraphView加以修改,de. SletObjet( oldpen);即執(zhí)行View 菜單的ClassWizard 命令,添加On-newpen. DeleteObject( ) ; }Paint的實(shí)現(xiàn)代碼使程序在視圖區(qū)繪制出X;和Y;的程序首先對(duì)ary數(shù)組元素進(jìn)行變換處理,結(jié)果散點(diǎn)圖(圖1 )。放在new數(shù)組中,然后new數(shù)組中的new[ ][0]和new[][1]進(jìn)行-元線性回歸的處理,求得回歸系數(shù)后對(duì)系數(shù)進(jìn)行變換以求得非線性回歸系數(shù)。最后對(duì)ary[ ] [0]進(jìn)行擬合,求得數(shù)組P-ER的各元素值和各種參數(shù)的值(程序代碼略)。4運(yùn)行結(jié)果圖1 x;與Y散點(diǎn)表2運(yùn)行結(jié)果與樣本對(duì)比及誤差從圖1看出它與指數(shù)函數(shù)、對(duì)數(shù)函數(shù)和雙曲線實(shí)際值/10*t預(yù)測(cè)值/1041誤差/%函數(shù)類似。用上述3種曲線進(jìn)行擬合,根據(jù)程序計(jì)128532.2131515. 82. 32139248. 0140770. 01.09算結(jié)果發(fā)現(xiàn)用雙曲線函數(shù)的擬合誤差最小,故選擇126211. 3126233. 40.02雙曲線函數(shù)。程序的部分代碼為:程序運(yùn)行后得出實(shí)際值X; ,預(yù)測(cè)值x';,和誤void CGraphView : : OnPaint( )差P-ER,如表2所示,曲線類型為:.{inti;1/y= a+ b/xCPen newpen ;系數(shù)c=0CPen * oldpen ;中國(guó)煤化工CPaintDC de ( this) ;YHCNMHGnewpen. CreatePen( Ps - SOL標(biāo)準(zhǔn)離差S = 0.002632ID,1 , RCB(0,0 ,0)) ;相關(guān)系數(shù)R =0. 862172oldpen = de. SelectObject( &newpen) ;de. SetMapMode(0) ;(下轉(zhuǎn)51頁(yè))12徐坤:圍巖穩(wěn)定性綜合 控制技術(shù)的應(yīng)用2005年第4期長(zhǎng)錨固,以保證錨固效果,預(yù)緊力60~70kN。徑由大逐漸變小,前端孔徑8mm ,后端孔徑4mm。(6)幫部最上面1根錨桿都向上傾斜30° (高(3 )漿液選用化學(xué)漿液新型堵水材料馬麗散。幫最.上面1根錨桿向上傾斜約50~60°,第2根錨(4)注漿量每孔4~8桶，注漿時(shí)間每孔.桿向上傾斜約15 ~25°，低幫最上面1根錨桿可水20min左右, 注漿壓力不小于8MPa。平安設(shè))，最下面1根錨桿都向下傾斜30 ~45°安(5)根據(jù)具體的滲水淋水位置施工注漿錨桿設(shè),頂板兩側(cè)錨桿在垂直頂板方向再向外傾15~再根據(jù)堵水效果決定是否加打注漿孔及是否補(bǔ)注,20°，以便于施工和有效控制圍巖。通常以1根注漿錨桿控制2m2左右的面積計(jì)算注漿(7)幫頂破碎處,施工2~4根右旋全螺紋預(yù)錨桿的數(shù)量。拉力錨桿加強(qiáng)，錨桿長(zhǎng)度為2000mm,直徑為(6)注漿作業(yè)滯后迎頭10 ~20m ,在50m范18mm ,采用全長(zhǎng)錨固方式，每根錨桿采用2節(jié)圍內(nèi)完成。Z2550型中速樹(shù)脂藥卷錨固。4支護(hù)效果. (8)幫頂錨索滯后迎頭5 ~ 8m安裝齊全。3.3噴混凝土參數(shù)為了觀測(cè)錨桿支護(hù)效果,研究支護(hù)參數(shù)的合理對(duì)無(wú)水區(qū)的頂、幫,為防止煤巖體風(fēng)化,直接性,巷道每隔30m設(shè)置1個(gè)測(cè)站,對(duì)圍巖表面位.在圍巖表面噴射薄層混凝土,噴層厚30~50mm。移、頂板離層狀況、錨桿載荷進(jìn)行了觀測(cè)。采用原在頂、幫潮解區(qū)噴混凝土，噴層厚70 ~ 100mm。有支護(hù)形式的巷道圍巖表面位移-般在1m左右,混凝土配比為水泥:黃沙:石子= 1:2:2。噴射混凝頂板離層不斷增加,直到巷道修復(fù)為止,錨桿也隨土作業(yè)滯后迎頭10~ 20m ,在50m范圍內(nèi)完成。著巷道變形迅速失效。采用綜合控制技術(shù)施工的3.4煤巖體內(nèi)注漿 參數(shù)7,44風(fēng)巷圍巖表面位移最大值為400mm,頂板離(1)在頂、幫滲流區(qū)進(jìn)行煤巖層內(nèi)部注漿，層最大值180mm ,錨桿受力穩(wěn)定在70 ~ 110kN。從孔徑42mm,孔深3.5m,高幫/頂板/底幫注漿孔礦壓觀測(cè)結(jié)果結(jié)合直觀支護(hù)效果分析，上述圍巖穩(wěn)采用2/2/1型布置方式。定性綜合控制技術(shù)支護(hù)方案在桃園煤礦7,44風(fēng)巷. (2)注漿錨桿長(zhǎng)度3000mm,用12.7mm鋼管的應(yīng)用,基本能夠滿足安全生產(chǎn)需要,試驗(yàn)研究獲制成。鋼管底端1.5 ~2. 0m長(zhǎng)度內(nèi)錯(cuò)開(kāi)鉆孔,孔得了成功。[責(zé)任編輯:顏立新](上接12頁(yè))定了最佳的函數(shù)曲線圖象,合理的確定了最優(yōu)的預(yù)從上面結(jié)果R=0. 976217→1 ,可知預(yù)測(cè)模型測(cè)函數(shù)及其參數(shù)，因此該方法用于煤炭需求量的可行,利用該模型預(yù)測(cè)2005~2008年的煤炭需求中、短期預(yù)測(cè),具有很高的可信度。量如表3。[參考文獻(xiàn)]表3煤炭需求預(yù)測(cè)結(jié)果年份20062008[1]王端武,王浩，張爍,等.我國(guó)煤炭需求預(yù)測(cè)[J].中需求量/104t 167641. 1179465.4 172839. 6161357.2國(guó)煤炭, 1999 ,25 (4): 12-15.[2]郭嗣宗,陳剛.信息科學(xué)的軟計(jì)算方法[M].沈陽(yáng):東北5結(jié)論大學(xué)出版社, 2001.[3]樓順天，于衛(wèi)，閏華梁. MATLAB程序設(shè)計(jì)語(yǔ)言[ M].西用VC++6.0建立的煤炭需求預(yù)測(cè)模型,克服安:西安電子科技大學(xué)出版社, 1997.了傳統(tǒng)理論推理視圖界面的缺點(diǎn),經(jīng)過(guò)反復(fù)實(shí)驗(yàn)確[責(zé)任編輯:鄒正立]xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx(上接36頁(yè))(2)回收支柱和頂梁前,必須要先拴牢拉緊5結(jié)語(yǔ)后，方可回收,以防回下的鐵料竄下傷人。采用以上的管理方法. T314-'工作面安全實(shí)現(xiàn)(3)回柱時(shí),如遇到頂板來(lái)壓,要立即停止初次中國(guó)煤化工面收作,采出煤炭20作業(yè),待壓力穩(wěn)定后方可繼續(xù)回柱工作。余萬(wàn)YHCN M H沒(méi)有發(fā)生大的冒頂、倒(4)回柱不準(zhǔn)和打眼、裝藥等工序平行作業(yè)。棚、傷人等事故,取得了較好的經(jīng)濟(jì)效益。(5)回柱時(shí)，溜子道要填滿、閘死、蹬實(shí),[責(zé)任編輯:崔德仁]防止從溜子道濺出的煤塊傷人。51