孔新海，劉志斌，魏 勇

0 引言

灰色系統(tǒng)理論自誕生以來,已在許多領(lǐng)域得到了廣泛的應(yīng)用,特別是在社會(huì)、經(jīng)濟(jì)預(yù)測(cè)等方面應(yīng)用更為廣泛。作為灰色系統(tǒng)理論重要內(nèi)容之一的GM(1,1)模型,其應(yīng)用價(jià)值在越來越多的領(lǐng)域中得到體現(xiàn)。從GM(1,1)模型建模機(jī)理來看,對(duì)遞增離散點(diǎn)列x(0)，其一次累加生成序列x(1)是單調(diào)遞增的，我們常直接建模即可；對(duì)遞減離散點(diǎn)列x(0)，由于其一次累加序列x(1)也是單調(diào)遞增的，其模擬值也是遞增的，如果對(duì)進(jìn)行累減還原，常常被認(rèn)為會(huì)產(chǎn)生不合理的計(jì)算誤差，因此很多學(xué)者通過各種變換，如：反向累加生成、倒數(shù)變換、線性變換，把遞減序列變換成遞增序列，再建立GM(1,1)模型。其實(shí)，遞減序列也可以建立GM(1,1)模型，有時(shí)建模效果還非常好。本文在文獻(xiàn)[1][2]的基礎(chǔ)上，提出一種基于單調(diào)遞減序列的離散變換，以期提高序列的光滑度。

1 基本概念

定義 1設(shè)序列 X={x(1),x(2),…,x(n)}，稱 σx(k)=為序列X的級(jí)比；稱δx(k)=|1-σx(k)|為序列X的級(jí)比偏差；稱為序列X的光滑度。定義2若為兩單調(diào)遞減序列,且對(duì)?k=2,3,…,n,有 ρa(bǔ)(k)＜ ρb(k),則稱序列更光滑。

定理1對(duì)于單調(diào)遞減序列 x(k)＞0,非負(fù)變換y(k)=F(x(k))滿足不等式 ρx(k)＜ρy(k)的充要條件為F(x(k))可以表示為 F(x(k))=x(k)f(k)，其中序列 f(k)非負(fù)，且嚴(yán)格單調(diào)上升。

證明：先證充分性。由于x(k)、f(k)非負(fù)，且 f(k)嚴(yán)格單調(diào)上升，則對(duì)i=1,2,…,k-1，有

x(k)x(i)f(i)＜x(k)f(k)x(i)

兩邊對(duì)i求和，可得

即

從而

再證必要性。由

取k=2，得

x(2)/x(1)＜F(x(2))/F(x(1))

即有

F(x(1))/x(1)＜F(x(2))/x(2)

令 z(k)=x(k+1)，則

F(z(1))/z(1)＜F(z(2))/z(2)

顯然有

F(x(2))/x(2)＜F(x(3))/x(3)

2 基于單調(diào)遞減序列的離散變換

2.1 離散變換

定義3對(duì)于單調(diào)遞減的原始序列X={x(1),x(2),…,x(n)}，進(jìn)行離散變換y(k)=F(x(k))。

(1)若X為減速下降，即σx(1)＞σx(2)＞…＞σx(n-1)＞1，則變換算式為,d＞1,0＜s＜1。

(3)若X為加速下降，即1＜σx(1)＜σx(2)＜…＜σx(n-1)，則變換算式為：

令 φ(k)=(k+1)s-ks，當(dāng) σx(k)遞減，則 φ(k)也要遞減，即0＜s＜1；當(dāng)σx(k)遞增，則φ(k)也要遞增，即 s＞1。

令 ?(s)=(k+1)s-ks,則 d(k+1)s-ks=d?(s),所以 d 與 s的增減具有反向性。若d放大，s縮?。蝗鬱縮小，s放大。

性質(zhì)(1)為了保持變換序列與原始序列單調(diào)性一致，當(dāng) X為單調(diào)遞減序列時(shí)，對(duì)任意的k，d和 s滿足d(k+1)s-ks≤σx(k)；(2)該變換能夠縮小級(jí)比偏差；(3)該變換還原誤差不變。

證明：(1)當(dāng) X為單調(diào)遞減序列時(shí)，由于σy(k)=對(duì) 任 意 的 k,要 使y(k)≥ y(k+1),則,即 d(k+1)s-ks≤ σx(k)。

下面給出一種確定d和s的方法（級(jí)比偏差之和最小法）。

因?yàn)?σy(k)=σx(k)dks-(k+1)s，所以

δy(k)=|1-σy(k)|=|1-σx(k)dks-(k+1)s|

(1)d(k+1)s-ks≤σx(k)；(2)d＞1；(3)s的范圍。

2.2 建立GM(1,1)模型

設(shè)X(0)={x(0)(1),x(0)(2),…,x(0)(n)}為單調(diào)遞減的原始序列，分析其級(jí)比變化趨勢(shì)，利用級(jí)比偏差之和最小確定d與 s，再進(jìn)行數(shù)據(jù)變換 y(k)=x(k)dks，對(duì)變換序列Y(0)={y(0)(k)}建立GM(1,1)微分方程模型

其中：

求解上述微分方程的時(shí)間響應(yīng)函數(shù)并離散化，則

還原得

3 實(shí)例分析

例1[3]某10mg血藥在相同時(shí)間間隔測(cè)得的濃度數(shù)據(jù)為

其級(jí)比序列為

這顯然是個(gè)增長速度減緩的序列。利用本文提出的離散變換y(k)=x(k)dks，d＞1,0＜s＜1，對(duì)原始數(shù)據(jù)進(jìn)行離散變換處理，然后再建立GM(1,1)模型。根據(jù)級(jí)比偏差之和最小,可求得d=1.466919,s=0.7264035,建立GM(1,1)模型，模擬效果如表1所示。

表1 模擬效果比較

從表1可以看出，經(jīng)本文方法變換之后建立的GM(1,1)模型模擬效果比原始模型的結(jié)果要好。

例2[4]在N2O5的熱分解反應(yīng)過程中,在等時(shí)間間隔里實(shí)驗(yàn)測(cè)得N2O5的分壓數(shù)據(jù)為X(0)={24.7,18.5,14.0,10.5,7.8,5.8,4.4,3.3},這里溫度T=318.2K ，單位：cm Hg，其級(jí)比序列為={1.3351,1.3214,1.3333,1.3462,1.3448,1.3182,1.3333}，這是個(gè)相對(duì)平穩(wěn)的序列。這里令s=1，取，建立GM(1,1)模型，模擬效果如表2所示。

表2 模擬效果比較

例3[5]某油田在遞減階段的年產(chǎn)量為X ={14.38,14,13.06,12.09,10.98,9.88,8.27}，單位：107t/a,其級(jí)比序列為={1.0271,1.0720,1.0802,1.1011,1.1113,1.1947這是個(gè)級(jí)比增長的遞減序列。利用本文提出的離散變換y(k)=x(k)dks(d＞1,s＞1)，對(duì)原始數(shù)據(jù)進(jìn)行離散處理，然后再建立GM(1,1)模型。根據(jù)級(jí)比偏差之和最小,可求得d=1.008491,s=2.057443,建立GM(1,1)模型，模擬效果如表3所示。

表3 模擬效果比較

4 結(jié)語

針對(duì)原始序列為單調(diào)遞減的情況，本文提出了一種提高光滑度的離散變換，這種變換能更好地適應(yīng)原始數(shù)據(jù)序列的變化趨勢(shì)。根據(jù)原始序列的增長速率來調(diào)整離散變換算式的參數(shù)，以縮小原始數(shù)據(jù)序列的級(jí)比偏差，提高GM(1,1)建模精度。文中所舉實(shí)例也說明了這種離散變換的可行性。本文只是針對(duì)原始序列為單調(diào)遞減的情況進(jìn)行了討論，其他情況將在其后作進(jìn)一步研究。

[1] 黃福勇.灰色系統(tǒng)建模的變換方法[J].系統(tǒng)工程理論與實(shí)踐,1994,14(11).

[2] Yong Wei,Yi Zhang.An Essential Characteristic of the Discrete Function Transformation to Increase the Smooth Degree of Data[J].The Journal of Grey System(UK),2007,18(3).

[3] 楊保華,張忠泉.灰色GRM(1,1)藥物動(dòng)力學(xué)模型[J].數(shù)理醫(yī)藥學(xué)雜志,2003,16(1).

[4] 安燕,吳啟勛.倒數(shù)累加生成灰色RGM(1,1)模型及其在化學(xué)動(dòng)力學(xué)上的應(yīng)用[J].計(jì)算機(jī)與應(yīng)用化學(xué),2006,23(2).

[5] 陳元千,郭二鵬.新型油田產(chǎn)量遞減模型的建立與應(yīng)用[J].中國海上油氣,2008,20(6).