唐榮桂王一帆盧曦

函數(shù)置換條件下P-Ⅲ參數(shù)的最小二乘計(jì)算

唐榮桂1王一帆2盧曦2

一、引言

眾所周知，因P-Ⅲ密度函數(shù)的原函數(shù)目前還無法求得，其參數(shù)只能應(yīng)用一些估計(jì)的方法獲取，如矩法、三點(diǎn)法、極大似然法、權(quán)函數(shù)法、線性矩法和適線法等，大部分方法誤差大，有的隨意性大，可信度低。在尊重樣本的基礎(chǔ)地位和樣本誤差，同時(shí)考慮到散點(diǎn)函數(shù)擬合的準(zhǔn)確度與唯一性的基礎(chǔ)上，筆者研究出在函數(shù)置換條件下利用最小二乘法原理進(jìn)行P-Ⅲ參數(shù)計(jì)算的方法。

二、原理與方法

1.目標(biāo)函數(shù)的建立

考察P-Ⅲ的密度函數(shù)：

令u=x-a0，對(duì)于其中的指數(shù)函數(shù)y=e-βu，取區(qū)間[a，b]，區(qū)間內(nèi)可近似用滿足兩端點(diǎn)連續(xù)（函數(shù)值相等）與平滑（一階導(dǎo)數(shù)值相等）的三次函數(shù)

進(jìn)行替代，三次函數(shù)的系數(shù)方程為：

對(duì)于發(fā)生在[a，b]內(nèi)事件的概率為：

如果把P-Ⅲ的密度函數(shù)劃分成很多小區(qū)間，各個(gè)區(qū)間的概率依次為ΔPj，就有如下累積概率公式：

對(duì)應(yīng)于P-Ⅲ的密度函數(shù)，樣本系列序號(hào)i從有限端（?。┑綗o限端（大）排列，給出如下樣本的有關(guān)計(jì)算公式：

這樣就可以根據(jù)公式（2）-（6）計(jì)算樣本的累積頻率Pei、模比系數(shù)Ki、有限端點(diǎn)坐標(biāo)a0、平移后坐標(biāo)ui（以0為起點(diǎn)）、部分概率ΔPi以及理論上的累積概率Pti。若要達(dá)到樣本點(diǎn)與理論曲線的最佳擬合，根據(jù)誤差的最小二乘原理，必然滿足如下目標(biāo)函數(shù)：

（7）式是多系數(shù)的復(fù)合函數(shù)方程式，沒有辦法利用求導(dǎo)數(shù)的辦法建立參數(shù)方程組，也就不可能用樣本統(tǒng)計(jì)計(jì)算的方法求解參數(shù)。下面將研究（7）式解的問題。

2.P-Ⅲ參數(shù)解的唯一性

P-Ⅲ的密度函數(shù)有3個(gè)部分組成，冪函數(shù)控制增率，指數(shù)函數(shù)控制減率，常數(shù)項(xiàng)是保證全體事件的概率為1的約束條件。因此，密度函數(shù)對(duì)于α始終是增函數(shù)，對(duì)于β又始終是減函數(shù)，都具有單調(diào)性，它們對(duì)分布函數(shù)的影響同樣具有單調(diào)性，因此，對(duì)于固定樣本，（7）式成立時(shí)兩個(gè)參變量必然只有唯一的平衡點(diǎn)，即（7）式只有唯一解。

3.簡單的求解方法——方向搜索法

如果利用計(jì)算機(jī)程序，依據(jù)（7）式，以不同的參數(shù)組合（排除不符合P-Ⅲ特征）進(jìn)行計(jì)算，尋找到最小值，即可獲得目標(biāo)函數(shù)的解。這里再給出簡單的解法——方向搜索法（只適用于具有單調(diào)性質(zhì)的參數(shù)），該方法可以減少大量的計(jì)算工作量，快速獲得高精度結(jié)果。具體步驟如下：

①選定方程可能解的初值，如P-Ⅲ參數(shù)，可以利用其他方法給出的估計(jì)值作為初值；

②根據(jù)精度要求，設(shè)定參數(shù)的變化步長，可以一輪搜索，也可以多輪搜索；

③以一個(gè)參數(shù)固定，另一參數(shù)變化，通過（6）計(jì)算出數(shù)值，向數(shù)值減少方向逐步搜索到最小值點(diǎn)；

④再以該點(diǎn)為中心計(jì)算其他6個(gè)方位點(diǎn)上的值，如果沒有比其更小的值，該點(diǎn)對(duì)應(yīng)的參數(shù)即為所解。如果有比其還小的點(diǎn)，就沿著這兩個(gè)數(shù)值較小的點(diǎn)的減少方向繼續(xù)搜索，直至沒有更小的值為止。

三、應(yīng)用舉例

以江蘇省白馬湖地區(qū)（高良澗閘、運(yùn)東閘、阮橋閘三站平均）最大3日雨量的P-Ⅲ參數(shù)計(jì)算為例。這里主要介紹應(yīng)用

excel完成的計(jì)算過程。

1.制作誤差平方和計(jì)算模板（excel工作簿）

該簿由兩塊組成。

第一塊：坐標(biāo)平移計(jì)算和各個(gè)區(qū)間置換公式系數(shù)計(jì)算。樣本38年，變量u從0到38將密度曲線分成38個(gè)部分，建立38張工作表（1-38），每張表有一個(gè)固定部分和一個(gè)活動(dòng)部分。

固定部分，設(shè)立α、β輸入窗口，當(dāng)參數(shù)改變，密度函數(shù)的起點(diǎn)坐標(biāo)和各樣本點(diǎn)的距0值需要相應(yīng)計(jì)算。38張表的這部分都相同（直接復(fù)制），但表（2-38）的兩個(gè)輸入窗口改為調(diào)用表1相應(yīng)單元格的值。

活動(dòng)部分，目的是分別計(jì)算兩點(diǎn)間置換函數(shù)的系數(shù)，這里調(diào)用了excel行列式計(jì)算函數(shù)，同樣復(fù)制表（1）的活動(dòng)部分，后面各表只需在相同的列號(hào)欄依次向下移動(dòng)一行進(jìn)行粘貼即可。

第二塊：誤差平方和計(jì)算。工作簿中，專門設(shè)置一張計(jì)算誤差平方和的表格。在公式輸入時(shí)，α、β必須調(diào)用表1的指定單元格數(shù)據(jù)，Γ（α）可調(diào)用excel中的GAMMALN函數(shù)，ΔPi中各行的置換函數(shù)系數(shù)可依次調(diào)用表（1-38）中的相同單元格的數(shù)據(jù)（用活動(dòng)單元格，僅需改動(dòng)工作表號(hào)）。在表格的最末端有一個(gè)誤差平方總和欄，它是整個(gè)工作簿的輸出窗口，α、β變動(dòng)后的數(shù)字就取自該單元格。

2.方向搜索

根據(jù)誤差平方和最小來找出指定精度下的P-Ⅲ參數(shù)。

第一步：選擇參數(shù)初始值。利用武漢大學(xué)提供的《水文頻率分布曲線適線軟件》求得的樣本的參數(shù)為，本文計(jì)算中使用的是模比系數(shù)，x=1，利用公式（8）計(jì)算得：α=4，β=5.9。

第二步，選擇適合的步長進(jìn)行搜索。按照三位有效數(shù)字精度，分2次搜索。由于計(jì)算的誤差平方和隨著參數(shù)的變化規(guī)律明顯，只需要選擇代表性點(diǎn)，按照一定方向，很快就能取得結(jié)果。

粗搜索：以初始值為定點(diǎn)，先縱向搜索，找到最小值點(diǎn)，再向臨側(cè)逐步找最小值點(diǎn)。經(jīng)過初步搜索，可以確定最小值在3.25〈α〈3.35及4.85〈β〈4.95區(qū)間內(nèi)。

精搜索：按照給定精度的最小單位作為表格的計(jì)算單位，以粗搜索得出的最小值點(diǎn)為定點(diǎn)，先縱向后橫向進(jìn)行搜索。根據(jù)最小值點(diǎn)得出：α=3.32，β=4.92，據(jù)此換算出江蘇省白馬湖地區(qū)最大3日降水量的P-Ⅲ參數(shù)為：

3.成果驗(yàn)證

為了檢驗(yàn)本文方法，下面給出兩種方法擬合圖進(jìn)行對(duì)照，見圖1。

圖1 江蘇省白馬湖地區(qū)最大3日降水量頻率曲線對(duì)照?qǐng)D

計(jì)算機(jī)軟件計(jì)算得到的參數(shù)應(yīng)該優(yōu)于其他方法的估計(jì)值，本文方法的擬合度（0.983）優(yōu)于計(jì)算機(jī)軟件適線的擬合度（0.978）；點(diǎn)線配合情況看，本文方法做到了兼顧全部樣本，略優(yōu)于計(jì)算機(jī)軟件配線；以200年一遇的設(shè)計(jì)標(biāo)準(zhǔn)為例，軟件方法比本文方法偏小6.2%，這個(gè)差別對(duì)該區(qū)域的排澇設(shè)計(jì)帶來的影響還是很大的。

4.誤差評(píng)估

與直接進(jìn)行最小二乘法計(jì)算相比，本方法唯一的誤差來源是函數(shù)置換，為檢驗(yàn)其影響，取本文的最終結(jié)果β=4.92來驗(yàn)證。方法：模比系數(shù)以0.001為步長，用對(duì)應(yīng)的分段函數(shù)和指數(shù)函數(shù)來計(jì)算相對(duì)誤差，得出的結(jié)果是：負(fù)偏為主，平均誤差：-0.61‰，最大偏離-4.7‰，最大偏離區(qū)間在u36～u38（該區(qū)間平均偏離值-2.0‰），一般偏離數(shù)值與間距同步，若間距在0.03以下時(shí)，偏差小于十萬分之一。就本例的情況，總體影響小于千分之一，可以忽略。

四、關(guān)于本文方法與常規(guī)方法比較

目前評(píng)判P-Ⅲ參數(shù)優(yōu)劣，主要通過視覺觀察點(diǎn)線配合情況，按說幾率格紙繪制P-Ⅲ曲線是無法用簡單的數(shù)學(xué)解析式表達(dá)，依靠視覺調(diào)試（適線法）當(dāng)然會(huì)造成很大誤差，還可能出現(xiàn)與P-Ⅲ理論相悖的結(jié)果。另外現(xiàn)行方法對(duì)設(shè)計(jì)值推算也很麻煩，一般的查用表又很粗糙，對(duì)非專業(yè)人員來說，使用P-Ⅲ有一定難度。

本文方法是建立在P-Ⅲ理論基礎(chǔ)之上，實(shí)測樣本配線又是用公認(rèn)的最小二乘準(zhǔn)則，數(shù)學(xué)理念清晰，且結(jié)果唯一，不會(huì)造成歧義。為了滿足視覺效果，也可以繪制累積頻率點(diǎn)和擬合的Ki～Pti線加以參照。

利用本文方法編制簡單的小程序（需將置換函數(shù)的間距設(shè)計(jì)的很小），可以方便求出高精度的P-Ⅲ參數(shù)，還可以根據(jù)參數(shù)及設(shè)計(jì)標(biāo)準(zhǔn)輸出設(shè)計(jì)結(jié)果。作者相信，隨著技術(shù)的發(fā)展，P-Ⅲ也會(huì)和其他函數(shù)一樣成為各種計(jì)算工具的內(nèi)置函數(shù)，方便人們的應(yīng)用■

（作者單位：1.江蘇省洪澤湖水利工程管理處2231002.江蘇省灌溉總渠管理處223200）