劉解放，劉思峰，吳利豐，方志耕（1.河南科技學(xué)院數(shù)學(xué)科學(xué)學(xué)院，河南　新鄉(xiāng)453003；.南京航空航天大學(xué)經(jīng)濟(jì)與管理學(xué)院，江蘇　南京10016）

分?jǐn)?shù)階反向累加離散灰色模型及其應(yīng)用研究

劉解放1，2，劉思峰2，吳利豐2，方志耕2
（1.河南科技學(xué)院數(shù)學(xué)科學(xué)學(xué)院，河南新鄉(xiāng)453003；
2.南京航空航天大學(xué)經(jīng)濟(jì)與管理學(xué)院，江蘇南京210016）

為了減小離散灰色預(yù)測模型解的擾動界，增加模型的解的穩(wěn)定性，并且充分利用系統(tǒng)的新信息，提出了分?jǐn)?shù)階反向累加離散灰色模型。通過數(shù)學(xué)推導(dǎo)，給出了模型具體的計(jì)算公式，并計(jì)算了其擾動界。通過比較一階反向累加離散灰色模型和分?jǐn)?shù)階反向累加離散灰色模型擾動界的大小，證明了分?jǐn)?shù)階反向累加離散灰色模型的優(yōu)越性。最后，將其應(yīng)用在某型號武器系統(tǒng)的儲存可靠性的預(yù)測中，計(jì)算結(jié)果優(yōu)于傳統(tǒng)灰色預(yù)測模型，預(yù)測精度有了大幅度提高，從而拓展了灰色預(yù)測模型的理論研究和實(shí)際應(yīng)用范圍。

灰色預(yù)測；分?jǐn)?shù)階；反向累加；可靠性

網(wǎng)址：w w w.sys-ele.co m

0 引 言

灰色系統(tǒng)理論主要研究“小樣本、貧信息”問題，近年來得到了學(xué)者們廣泛的關(guān)注和研究［1-8］?；疑到y(tǒng)理論最基本的方法是灰色生成，生成是建立灰色模型的基礎(chǔ)。累加生成是使灰色過程由灰變白的一種方法，其在灰色系統(tǒng)理論中占有極其重要的地位。通過累加生成可以看出灰量積累過程的發(fā)展態(tài)勢，使離亂的原始數(shù)據(jù)中蘊(yùn)含的積分特性或規(guī)律充分表現(xiàn)出來。

對于非負(fù)離散序列X（0），其一次累加生成序列X（1）是單調(diào)遞增的，因此，其擬合曲線（1）也應(yīng)該是單調(diào)遞增的才合理。但是，對于單調(diào)遞減的離散序列X（0），由于X（1）是單調(diào)遞增的，其模擬序列（1）也是單調(diào)遞增的，那么在對＾X（1）進(jìn)行累減生成還原時，就會產(chǎn)生不合理的計(jì)算誤差，基于此，文獻(xiàn)［9］提出了反向累加生成G O M（1，1）模型，由于反向累加生成序列為單調(diào)遞減序列，因此，在對于單調(diào)遞減序列進(jìn)行建模時能取得優(yōu)于G M（1，1）的計(jì)算效果。文獻(xiàn)［10-11］提出了參數(shù)的直接求解方法，簡化了計(jì)算過程，并將其應(yīng)用到了農(nóng)藥降解問題的建模中，取得了優(yōu)于原始G O M（1，1）的計(jì)算精度。在此基礎(chǔ)上，文獻(xiàn)［12］對于G O M（1，1）模型的初值和背景值進(jìn)行改進(jìn)，進(jìn)一步優(yōu)化了G O M（1，1）模型的參數(shù)求解方法，提高了模型的計(jì)算精度。文獻(xiàn)［13］給出了模型優(yōu)化的背景值，不僅提高了建模的精度，而且使得新模型對于遞減和遞增序列都適用，拓展了新模型的適用范圍。在文獻(xiàn)［14-15］離散G M（1，1）模型的基礎(chǔ)上，文獻(xiàn)［16］提出了離散G O M（1，1）模型，并給出了新模型的求解方法，證明了新模型具有白化指數(shù)律和線性不變性。上述對于G O M（1，1）模型的改進(jìn)，改善了該模型的計(jì)算效果，擴(kuò)大了其適用范圍，對于灰色理論的發(fā)展具有積極的理論和實(shí)際意義。

文獻(xiàn)［17-18］把分?jǐn)?shù)階的思想引入到灰色建模中，取得了不錯的效果。為了減小灰色預(yù)測模型的擾動界，增加模型的穩(wěn)定性，并且充分利用系統(tǒng)的新信息，提出了分?jǐn)?shù)階反向累加離散灰色模型，并將其應(yīng)用在某型號武器系統(tǒng)的儲存可靠性預(yù)測中，計(jì)算結(jié)果顯示，新模型的計(jì)算精度優(yōu)于傳統(tǒng)灰色預(yù)測模型，從而拓展了灰色預(yù)測模型的理論研究和實(shí)際應(yīng)用范圍。

1 一階反向累加離散灰色模型及其擾動分析

定義1 設(shè)原始非負(fù)序列X（0）=｛x（0）（1），x（0）（2），…，x（0）（n）｝為原始序列，X（1）=｛x（1）（1），x（1）（2），…，x（1）（n）｝為X（0）的一階反向累加生成序列，其中

則稱

為G O M（1，1）模型［9］。

定理1［19-20］設(shè)A∈Cm×n，b∈Cm，A是矩陣A的廣義逆矩陣，當(dāng)A的列向量線性無關(guān)時，方程‖Ax-b‖有唯一解。

定理2［19-20］設(shè)A∈Cm×n，b∈Cm，A是矩陣A的廣義逆矩陣，B=A+F，c=b+k∈Cn，假設(shè)方程‖Bx-c‖2= min與‖Ax-b‖2=min的解分別為x+h和x，若rank（A）= rank（B）=n，而且當(dāng)‖A‖2‖F(xiàn)‖2＜1時，有

其中

稱x（1）（k+1）=β1x（1）（k）+β2（k=1，2，…，n-1）為一階反向累加離散灰色模型。

定理3 一階反向累加離散灰色模型

x（1）（k+1）=β1x（1）（k）+β2（k=1，2，…，n-1）的參數(shù)β1，β2可以利用最小二乘估計(jì)來計(jì)算，即

定義2 設(shè)原始非負(fù)序列為

一階反向累加序列為

其中

證明 略。

定理4 離散灰色模型的解可以按照最小二乘法‖Y-Bx‖2=min進(jìn)行求解，假設(shè)模型的解為x。下面考慮系統(tǒng)存在擾動的情況，為了證明簡便，不妨假設(shè)系統(tǒng)的第二項(xiàng)有擾動，即（0）（2）=x（0）（2）+ε，那么，解的擾動界為

由于B的列向量線性無關(guān)，因此，方程‖Y-Bx‖2= min有唯一解x=YB。

因?yàn)?/p>

又因?yàn)?/p>

所以

根據(jù)定理2，可以得到

限于篇幅，此處不再具體推導(dǎo)。

從L（x（0）（t））的計(jì)算公式可以看出，其擾動界是數(shù)據(jù)量t的增函數(shù)，即

則，系統(tǒng)的擾動界呈現(xiàn)出遞增的特性，那么，當(dāng)新信息發(fā)生擾動之后，其擾動界會大于老信息發(fā)生擾動之后產(chǎn)生的擾動界，說明模型對于新信息的變化更加敏感，這和灰色系統(tǒng)理論的“新信息優(yōu)先原理”比較相符。因此，與正向累加生成相比，反向累加生成具有一定的優(yōu)越性。

2 r階反向累加離散灰色模型及其擾動分析

為了進(jìn)一步分析分?jǐn)?shù)階反向累加離散灰色模型的性質(zhì)，增加模型解的穩(wěn)定性，減小灰色預(yù)測模型的擾動界，本節(jié)建立了r階分?jǐn)?shù)階反向累加離散灰色預(yù)測模型，并且對其擾動界進(jìn)行研究。

定理5 假設(shè)原始非負(fù)序列

則r次反向累加生成序列為

證明 利用數(shù)學(xué)歸納法進(jìn)行證明

當(dāng)r=1時

當(dāng)r=2時

假設(shè)當(dāng)r=p的時，等式成立，即

那么，當(dāng)r=p+1時

綜上所述，根據(jù)數(shù)學(xué)歸納法即可得到

定義3 設(shè)序列X（0）、X（r）的定義同定理5，稱x（r）（k+1）= β1x（r）（k）+β2（k=1，2，…，n-1）為r階反向累加離散灰色模型。

為了分析r階反向累加離散灰色模型的穩(wěn)定性，以下分析其解的擾動情況。

定理6 假設(shè)r階反向累加離散灰色模型的解為x，根據(jù)反向累加生成的性質(zhì)和分析問題的方便，首先分析當(dāng)數(shù)據(jù)序列第n項(xiàng)發(fā)生擾動的情況，即＾x（0）（n）=x（0）（n）+ε，因?yàn)?/p>

因?yàn)?/p>

根據(jù)定理2，可以得到

當(dāng)數(shù)據(jù)序列的第n項(xiàng)發(fā)生擾動時，即＾x（0）（n）=x（0）（n）+ε，基于r階反向累加的離散灰色模型的擾動界為

從上面的證明過程可以推導(dǎo)出，當(dāng)數(shù)據(jù)序列的第t項(xiàng)發(fā)生擾動時，即＾x（0）（t）=x（0）（t）+ε，基于r階反向累加的離散灰色模型的擾動界為

限于篇幅，此處不再具體證明。

從第1節(jié)的計(jì)算結(jié)果可知，基于一階反向累加的離散灰色模型的擾動界為

因此，當(dāng)0＜r＜1時，基于r階反向累加的離散灰色模型的擾動界要小于基于一階反向累加的離散灰色模型的擾動界，從理論上來說，所建立的模型具有更好的穩(wěn)定性。

3 算例分析

某型號的武器裝備，其定期檢測數(shù)據(jù)如表1所示，運(yùn)用本文提出的分?jǐn)?shù)階反向累加離散灰色預(yù)測模型對系統(tǒng)可靠度的點(diǎn)估計(jì)進(jìn)行預(yù)測。其中，n表示檢測個數(shù)，k表示失敗個數(shù)，R表示可靠度的點(diǎn)估計(jì)。

表1 某型號武器裝備的試驗(yàn)信息

根據(jù)表1，可以計(jì)算得到某型號武器裝備的可靠度的點(diǎn)估計(jì)序列為

此處可靠度的點(diǎn)估計(jì)為沒有引入Bayes方法的計(jì)算結(jié)果。

為了和傳統(tǒng)的灰色預(yù)測模型進(jìn)行比較，本文選擇了G M（1，1）模型和D G M（1，1）模型作為比較的對象，不同模型的計(jì)算結(jié)果及相對誤差如表2所示，其中R D G M（1，1）表示本文提出的分?jǐn)?shù)階反向累加離散灰色預(yù)測模型，計(jì)算得到的平均絕對百分比誤差（mean absolute percent error，M A PE）如表2所示。根據(jù)遺傳算法求得最優(yōu)階數(shù)r=0.22，限于篇幅，具體的計(jì)算步驟不再逐步列舉。

表2 不同模型的計(jì)算結(jié)果和相對誤差

其中

從表2可以看出，當(dāng)原始數(shù)據(jù)的第3個數(shù)據(jù)出現(xiàn)擾動的時候，G M（1，1）模型和D G M（1，1）模型的預(yù)測相對誤差較大，而本文提出的模型，由于其擾動界較小，模型的解比較穩(wěn)定，計(jì)算結(jié)果的相對誤差較小，模型的模擬和預(yù)測精度較高。

4 結(jié) 論

本文提出了分?jǐn)?shù)階反向累加離散灰色模型，通過數(shù)學(xué)推導(dǎo)，得出了模型的求解公式，求出了其擾動界，通過和一階反向累加離散灰色模型進(jìn)行比較，發(fā)現(xiàn)其擾動界減小，從而提高了模型解的穩(wěn)定性。由于反向累加更加有效地利用了系統(tǒng)的新信息，所以模型的預(yù)測精度更高。最后，將分?jǐn)?shù)階反向累加離散灰色模型應(yīng)用在某型號武器裝備的儲存可靠度的預(yù)測中，取得了很好的模擬和預(yù)測效果，為此類武器系統(tǒng)的維修保障提供了理論上的支持，拓展了灰色預(yù)測模型的應(yīng)用領(lǐng)域，在理論和實(shí)際應(yīng)用上都有著重要的價值。

［1］Deng J L.Control problem of grey system［J］.Systems and Control Letters，1982，1（5）：288-294.

［2］Liu S F，Lin Y.Grey systemstheory and applications［M］.Berlin：Springer-Verlag，2010：107-148.

［3］Dang Y G，Liu S F，Liu B.The G M models thatx（1）（n）be taken as initial value［J］.Chinese Journal of Management Science，2005（1）：132-134.（黨耀國，劉思峰，劉斌.x（1）（n）為初始條件G M模型［J］.中國管理科學(xué)，2005（1）：132-134.）

［4］W ei Z，H e J M.Generalized G M（1，1）model and its application in forecasting of fuel production［J］.Applied Mathematical Modelling，2013，37（9）：6234-6243.

［5］Cui J，Liu S F，Zeng B，et al.Anovelgrey forecasting m odel and its optimization［J］.Applied Mathematical Modelling，2013，37（6）：4399-4406.

［6］Lin J，Lian R J.Design of a grey-prediction self-organizing fuzzy controller for active suspension systems［J］.A pplied Soft Computing，2013，13（10）：4162-4173.

［7］Li D C，Chang C J，Chen C C，et al.Forecasting short-term electricity consumption using the adaptive grey-based approach—an Asian case［J］.O mega，2012，40（6）：767-773.

［8］Samvedi A，Jain V.Agrey approach for forecasting in a supply chain during intermittent disruptions［J］.Engineering Applications of Artificial Intelligence，2013，26（3）：1044-1051.

［9］Song Z M，Deng J L.The accumulated generating operation in opposite direction and its use in grey model G O M（1，1）［J］.Systems Engineering，2001，19（1）：66-69.（宋中民，鄧聚龍.反向累加生成及灰色G O M（1，1）模型［J］.系統(tǒng)工程，2001，19（1）：66-69.）

［10］Liu J Y，Yang T X，Wang S L.Directed method for computing parameters of G O M（1，1）［J］.Journal of Jilin University （Engineering and Technology Edition），2003，33（2）：75-79.）（劉金英，楊天行，王淑玲.反向G O M（1，1）模型參數(shù)的直接求解方法［J］.吉林大學(xué)學(xué)報（工學(xué)版），2003，33（2）：75-79.）

［11］Liu J Y，Yang T X，W ang S L.Grey opposite-direction m odelG O M（1，1）and its application［J］.Journal of Jilin A gricultural University，2002，24（4）：91-94.（劉金英，楊天行，王淑玲.反向G O M（1，1）模型及其在農(nóng)藥降解分析中的應(yīng)用［J］.吉林農(nóng)業(yè)大學(xué)學(xué)報，2002，24（4）：91-94.）

［12］Yang Z，Ren P，Dang Y G.Grey opposite-direction accu m ulated generating and optimization of G O M（1，1）model［J］.Systems Engineering-Theory&Practice，2009，29（8）：160-164.（楊知，任鵬，黨耀國.反向累加生成與灰色G O M（1，1）模型的優(yōu)化［J］.系統(tǒng)工程理論與實(shí)踐，2009，29（8）：160-164.）

［13］Lian ZW，Dang Y G，Wang Z X.Properties of accumulated generating operation inopposite-direction and optimization of G O M（1，1）m odel［J］.Systems Engineering-Theory&Practice，2013，33（9）：2306-2312.（練鄭偉，黨耀國，王正新.反向累加生成的特性G O M（1，1）模型的優(yōu)化［J］.系統(tǒng)工程理論與實(shí)踐，2013，33（9）：2306-2312.）

［14］Xie NM，Liu S F.Discrete G M（1，1）and mechanism of grey forecasting m odel［J］.Systems Engineering-T heory&Practice，2005，25（1）：93-99.（謝乃明，劉思峰.離散G M（1，1）模型與灰色預(yù)測模型建模機(jī)理［J］.系統(tǒng)工程理論與實(shí)踐，2005，25（1）：93-99.）

［15］Xie N M，Liu S F.Research on the affine properties of discrete grey model［J］.Control and Decision，2008，23（2）：200-203.（謝乃明，劉思峰.離散灰色模型的仿射特性研究［J］.控制與決策，2008，23（2）：200-203.）

［16］Wang F X.Discrete G O M（1，1）model and its application［J］. Journal of Xinjiang University（Natural Science Edition），2012，29（1）：75-77.（王豐效.離散G O M（1，1）模型及其應(yīng)用［J］.新疆大學(xué)學(xué)報（自然科學(xué)版），2012，29（1）：75-77.）

［17］W u L F，Liu S F，Yao L G，et al.Grey system model with the fractional order accumulation［J］.Communicationsin Nonlinear Science and Numerical Simulation，2013，18（7）：1775-1785.

［18］W u L F，Liu S F，Chen D，et al.Using gray model with fractional order accu m ulation to predict gas emission［J］.Natural Hazards，2014，71（3）：2231-2236.

［19］Sun J G.M atrix disturbance analysis ［M］.Beijing：Science Press，1987：355-356.（孫繼廣.矩陣擾動分析［M］.北京：科學(xué)出版社，1987：355-356.）

［20］Stewart GM.O n the perturbation of pseudo-inverses，projections and linear square problems［J］.SIA M Review，1977（19）：634-662.

Fractional order reverse accumulative discrete grey model and its application

LIU Jie-fang1，2，LIU Si-feng2，W U Li-feng2，F(xiàn) A N G Zhi-geng2
（1.School of M athematical Science，H enan Institute of Science and Technology，Xinxiang 453003，China；2.College of Economics and M anagement，N anjing University of Aeronautics and Astronautics，N anjing 210016，China）

In order to reduce the perturbation bounds of the solution of discrete grey forecasting m odel，increase the stability ofthe solution ofthe m odel，and make full use of new information ofthe system，a fractional reverse accu m ulative discrete grey m odelis put forward.Its calculation form ula is deduced and its perturbation bounds are calculated.The specific calculation form ula is given by mathematical deduction，and the perturbation bounds are calculated.The superiority of the fractional order reverse accu m ulative discrete grey m odelis proved by co m paring the perturbation bounds ofit with the first-order reverse accu m ulative discrete grey m odel.Finally，it is applied in storage reliability prediction of a certain type of weapon system，and the calculation resultis better than the traditional grey prediction m odel.The prediction accuracy has im proved dramatically.The proposed model expands the scope of the theoretical research and practical appl ication of grey forecasting model.

grey prediction；fractional order；reverse accu m ulation；reliability

N 941.5

A

10.3969/j.issn.1001-506 X.2016.03.37

1001-506 X（2016）03-0719-06

2014-10-30；

2015-11-05；網(wǎng)絡(luò)優(yōu)先出版日期：2015-11-20。

網(wǎng)絡(luò)優(yōu)先出版地址：http：∥w w w.cnki.net/kcms/detail/11.2422.T N.20151120.1801.012.html

歐盟第7研究框架瑪麗·居里國際人才引進(jìn)計(jì)劃（FP7-PIIF-G A-2013-629051）；國家自然科學(xué)基金（91324003，71171113，71271226）；國家自然科學(xué)基金與英國皇家學(xué)會國際合作交流項(xiàng)目（71111130211）；國家社會科學(xué)基金重點(diǎn)項(xiàng)目（12AZD102）資助課題

劉解放（1980-），男，講師，博士，主要研究方向?yàn)榛疑到y(tǒng)理論、復(fù)雜裝備研制管理。

E-mail：liujf101@126.com

劉思峰（1955-），男，教授，博士研究生導(dǎo)師，主要研究方向?yàn)榛疑到y(tǒng)理論、系統(tǒng)工程。

E-mail：sfliu@nuaa.edu.cn

吳利豐（1983-），男，博士研究生，主要研究方向?yàn)榛疑到y(tǒng)理論、復(fù)雜裝備研制管理。

E-mail：wlf6666@126.com

方志耕（1962-），男，教授，博士研究生導(dǎo)師，主要研究方向?yàn)橘|(zhì)量管理、系統(tǒng)工程。

E-mail：zhigengfang@163.com