帶權(quán)最小二乘多項式擬合的泛克立格法

帶權(quán)最小二乘多項式擬合的泛克立格法

克立格法是一種基于隨機過程的統(tǒng)計預(yù)測方法。理論上，克立格估計具有最小估計散射差異的統(tǒng)計特征。此外，克立格應(yīng)用中最重要的是半變異函數(shù)的確定，即協(xié)方差函數(shù)的確定。然而，由于觀測值較低，在實踐中難以應(yīng)用更精確的協(xié)方差函數(shù)?？肆⒏窆浪闶欠窬哂凶钚』烙嫴町惖慕y(tǒng)計特征？從統(tǒng)計學(xué)的角度來看，在某些情況下，它是最小差分的最優(yōu)計算公式，但從數(shù)學(xué)的角度來看，它只是一種相對特殊的近似方法，即通過旋轉(zhuǎn)曲線生成的基函數(shù)。這已經(jīng)被許多科學(xué)家證明。在這項工作中，我們將重點關(guān)注克立格估算的分析和解釋。1系數(shù)矩陣估計的可恢復(fù)算法設(shè)以節(jié)點(xi,yi)為中心生成的旋轉(zhuǎn)面為以此為基礎(chǔ)可建立其線性組合再設(shè)數(shù)據(jù)點為zi=z(xi,yi),i=1,2…,n,在數(shù)據(jù)點上可建立如下插值條件寫成矩陣形式為式中:適當(dāng)?shù)剡x擇節(jié)點和基函數(shù)可保證系數(shù)矩陣非奇異,因此由(3)式可得將(4)式代入(2)式可建立未觀測點(xp,yp)的插值函數(shù)為如將基函數(shù)取為各向同性的協(xié)方差函數(shù)則有它形式上與已知數(shù)學(xué)期望的簡單克立格估計相同.2克立格估計設(shè)z(x)是具有常數(shù)數(shù)學(xué)期望的平穩(wěn)隨機過程,我們的目的是通過觀測值zi=z(xi),i=1,2,…,n的加權(quán)線性組合最優(yōu)地確定z(xp),式中wi為待定的克立格權(quán),利用條件極值可推得確定系數(shù)wi的方程組為式中:μ為拉格朗日乘數(shù),引入矩陣記號(8)式可改寫為權(quán)系數(shù)與拉格朗日乘數(shù)可分別求解.由(9c)式有因為ITC-1I只是一個數(shù),故有(13)式再回頭代入(11)式可得再將(14)式代入(8)式可得記M可理解為觀測值的廣義平均值,當(dāng)觀測值獨立C=E或觀測值弱相關(guān)C≈E,則有即為簡單平均值.在(15)式兩端同時減去M可得(18)式與(7)式相對比可知帶一個未知均值的克立格估計可解釋為觀測值與觀測值的廣義均值之差被旋轉(zhuǎn)曲面生成的基函數(shù)所逼近.3拉格朗日乘數(shù)估計設(shè)平穩(wěn)隨機過程z(x)=m(x)+s(x)由二部分組成,第一部分m(x)為確定性漂移,第二部分s(x)為平穩(wěn)隨機信號.一般漂移可表示為多項式形式式中:fT(x)=[f0(x),f1(x),…,fk(x)]=[1,x,y,xy…],a=(a0,a1,…,ak)T.令zi,i=1,2,…,n為觀測值,z(xp)為待估值,并可由觀測值的加權(quán)線性組合求出由克立格估計準(zhǔn)則應(yīng)有限于篇幅,略去推導(dǎo)過程可得求解克立格權(quán)與拉格朗日乘數(shù)的方程組為式中:μ=(μ0,μ1,…,μk)T,fT=[f0(x),f1(x),…,fk(x)]T,F為由漂移部分組成的系數(shù)陣.C非奇異,由(22)式可得(23)式代入(22)式,易知然后將(24)式再次回代至(23)式,可知克立格權(quán)系數(shù)矢量為(25)式代入(20)式有記它就是(19)式的系統(tǒng)參數(shù),因此(26)式可改寫為或(29)式寫成元素形式為z(xp)-fpTa=(30)式仍與(7)式類似,因此可以說在除去系統(tǒng)部分后,泛克立格與使用旋轉(zhuǎn)面作為基函數(shù)的函數(shù)逼近類似.此外與(22)式直接解泛克立格方程組相比,本文(27),(28)式先求系統(tǒng)漂移部分再求隨機信號部分,既降低了求逆矩陣階數(shù),又提高了計算效率,并且也易于區(qū)分漂移部分與隨機部分的大小量級.4協(xié)方差參數(shù)值本節(jié)我們基于一定的數(shù)學(xué)分析觀點來分析協(xié)方差函數(shù)參數(shù)的選取.為討論方便以一維為例,(5)式略去下標(biāo)p可改寫為(31)式表明預(yù)測值只是函數(shù)k(x)移位后的線性組合.取k(x)為最常用的Gauss函數(shù)Gauss函數(shù)衰減很快,它的局部支持區(qū)間受到了參數(shù)d的強烈控制.不同參數(shù)的Gauss函數(shù)可示于圖1,從圖1可看出,當(dāng)d較小時,Gauss函數(shù)類似于不連續(xù)的脈沖函數(shù)δ(x).作為理想化的一個特例,我們考查ai≡1,xi=i時不同參數(shù)的Gauss函數(shù)之線性組合構(gòu)成的函數(shù)形態(tài).d=1.0,d=0.5和d=0.1時Gauss函數(shù)移位后的線性組合可示于圖2.實際上ai≡1隱含著采樣值(或觀測值)為一水平直線的采樣,理想的逼近結(jié)果也應(yīng)為一條直線.但事實并非如此,圖2c中的曲線類似于一把梳子,這是由于d=0.1時Gauss函數(shù)的支撐區(qū)間太窄所致.圖2b中的曲線已較圖2a好得多,但仍有許多連續(xù)的極大值和極小值,類似于一條正余弦曲線.只有圖2a達(dá)到了較理想的效果.這些事實說明協(xié)方差函數(shù)的逼近效果的確與其參數(shù)有很大關(guān)系,并且可以肯定地說合理的參數(shù)是存在的.首先我們假定信噪比或噪聲部分已被確定或已另外考慮,此處我們只考慮連續(xù)的信號部分.準(zhǔn)則1協(xié)方差參數(shù)應(yīng)使得證明:(1)d=min意味著Gauss函數(shù)具有較小的支撐區(qū)間,由此形成的矩陣帶寬將比較小.(2)在假定z(x)平穩(wěn)的情況下,半節(jié)點的預(yù)測值與節(jié)點上的預(yù)測值應(yīng)有相同的數(shù)學(xué)期望,即在E{z(x)}=C(C為常數(shù))或limC※0)的情況下,由于ai為觀測值的線性組合,有故(35)式代入(34)式有(36)式稍加變形也即(33)式,此處(33)式亦可從代數(shù)上理解為一條水平直線被核函數(shù)移位生成的基函數(shù)所逼近時,其半節(jié)點處的值應(yīng)與節(jié)點處的值相等.證畢.為求出合理的協(xié)方差參數(shù)值,引入比值函數(shù)式中考慮到Gauss函數(shù)衰減很快,對區(qū)間[-4,4]以外的Gauss函數(shù)已視其為零略去不計.為更好地了解比值函數(shù)隨d的變化規(guī)律,圖3畫出了它的變化曲線.從圖3可得出兩條主要的結(jié)論:(1)合理的參數(shù)應(yīng)位于區(qū)間[0.8,1.7],它可以使半節(jié)點與節(jié)點處的協(xié)方差函數(shù)線性組合大致相等;(2)d<0.5的參數(shù)應(yīng)盡量避免使用,因為它使半節(jié)點處的協(xié)方差函數(shù)的線性組合遠(yuǎn)小于節(jié)點處協(xié)方差函數(shù)的線性組合.準(zhǔn)則2Gauss函數(shù)的參數(shù)可比照B樣條函數(shù)確定.說明:B樣條函數(shù)在函數(shù)擬合或逼近中具有許多優(yōu)點因而獲得了非常廣泛的應(yīng)用.n次B樣條在數(shù)學(xué)上可方便地由矩形脈沖經(jīng)n-1次卷積得到式中:*表示卷積,π(x)是矩形脈沖函數(shù)特別地,對3次B樣條有除有限支撐外,B樣條與Gauss函數(shù)比較相似,都為鐘形函數(shù).Martheon和Dubrule曾討論過克立格估計和樣條間的等效性和相似性,特別是Unser等研究了Gauss函數(shù)與B樣條函數(shù)之間的相似關(guān)系,并建立了如下逼近公式取n=3,3次B樣條與Gauss函數(shù)已對比列于圖4.(41)式Gauss函數(shù)的參數(shù)轉(zhuǎn)化為常用的參數(shù)形式有不同階數(shù)的樣條對應(yīng)的Gauss函數(shù)參數(shù)可列于表1.一般3～7次的樣條均可取得較好的逼近效果,否則階次過低或過高的逼近結(jié)果將不是太粗糙就是太平滑或支撐區(qū)間過寬.因此由類比可知Gauss函數(shù)的參數(shù)亦應(yīng)與此相應(yīng),宜取為[0.82,1.16].這與我們前述的準(zhǔn)則1的估計是相一致的.本節(jié)以上推導(dǎo)是以格網(wǎng)為單位長度導(dǎo)出的,如格網(wǎng)非單位長度時,可通過下列變換化格網(wǎng)為單位長度式中:Δx,Δy分別為格網(wǎng)長度,x,y和x′,y′分別為變換前后的坐標(biāo)在實際應(yīng)用時,數(shù)據(jù)一般是非格網(wǎng)散布在一定區(qū)域.設(shè)S為區(qū)域面積,N為域內(nèi)數(shù)據(jù)個數(shù),在假定x和y方向無異向性時等效的格網(wǎng)長度可直觀地由下式求出:試算證明用這種方法確定散亂數(shù)據(jù)的等效格網(wǎng)長度效果比較好.5gauss函數(shù)的海淡設(shè)理論模型為:f(x)=sin2πx,并令采樣間距Δx=0.0625,在0<x<1區(qū)間共有17個數(shù)據(jù)點,以這些數(shù)據(jù)點為基礎(chǔ),用d=0.1,d=0.5和d=1.0(以Δx=0.0625為單位長度)時的Gauss函數(shù)去逼近該理論模型,圖5畫出了相應(yīng)的逼近曲線示意圖.觀察圖5可以發(fā)現(xiàn)d=0.1和d=0.5時Gauss函數(shù)的逼近曲面產(chǎn)生了許多不應(yīng)有的低谷和高峰,與原理論模型相差較大,這是由于d=0.1和d=0.5時Gauss函數(shù)支撐區(qū)間太窄,在半節(jié)點上達(dá)不到應(yīng)有的預(yù)期結(jié)果,在半節(jié)點上出現(xiàn)了許多不應(yīng)有的極小值.而d=1.0時Gauss函數(shù)產(chǎn)生了非常好的逼近效果,幾乎與理論模型一致.6協(xié)方差函數(shù)的函數(shù)應(yīng)該說平穩(wěn)過程有著確定的譜S(w),通過Wiener＿Chintsch定理它又有著確定的協(xié)方差函數(shù)C(x),式中F和F-1分別為Fourier正演和反演.但是從另一方面來看,在一定條件下克立格可理解為由旋轉(zhuǎn)曲面生成的基函數(shù)進行的函數(shù)逼近,該基函數(shù)并不一定要取為被逼近過程的協(xié)方差函數(shù),具有一定的可選擇性,并且,由于矩陣求逆的數(shù)值穩(wěn)定性等原因這種逼近還與觀測數(shù)據(jù)的取樣密度存在一定的內(nèi)在聯(lián)系.對于同樣一個隨機過程,采樣密度變了,基函數(shù)(即“協(xié)方差函數(shù)”)中的參數(shù)亦應(yīng)相應(yīng)調(diào)整,否則將導(dǎo)致數(shù)值不穩(wěn)定性出現(xiàn)(如矩陣求逆失敗).:(1)在一定程度上與采樣數(shù)據(jù)的密度(或間隔)有一定聯(lián)系.(2)以一維Gauss函數(shù)為例,設(shè)Δx為數(shù)據(jù)間隔,則其最優(yōu)參數(shù)應(yīng)取d=[0.8Δx,1.7Δ