[發(fā)明專利]一種基于經(jīng)驗(yàn)?zāi)B(tài)分解的自動氣象站數(shù)據(jù)質(zhì)量控制方法

有效

• 申請?zhí)枺?/td>: 201310159628.1
• 申請日: 2013-05-03
• 公開（公告）號: CN103278867A
• 公開（公告）日: 2013-09-04
• 發(fā)明（設(shè)計(jì)）人: 葉小嶺;張穎超;熊雄;張齊東;胡凱;孫寧
• 申請（專利權(quán)）人: 南京信息工程大學(xué)
• 主分類號: G01W1/00
• 分類號: G01W1/00
• 代理公司: 南京經(jīng)緯專利商標(biāo)代理有限公司 32200
• 代理人: 朱小兵

1.一種基于經(jīng)驗(yàn)?zāi)B(tài)分解的自動氣象站數(shù)據(jù)質(zhì)量控制方法，其特征在于：包括以下步驟：

步驟1.采集自動氣象站溫度觀測數(shù)據(jù)x(t)，其中t為采樣時(shí)間，且t為大于0的自然數(shù)；

步驟2.對x(t)進(jìn)行基本質(zhì)量控制，得到基本控制后的數(shù)據(jù)

步驟3.根據(jù)公式對基本控制后的數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解，得到本征模分量IMF_i與的趨勢項(xiàng)r；其中k為大于1的自然數(shù)，i∈k；

步驟4.計(jì)算本征模分量IMF_i的能量E_i，根據(jù)E_i繪制能量柱狀圖，并觀察能量柱狀圖，得到本征模分量能量的階躍點(diǎn)M，根據(jù)該階躍點(diǎn)確定存在噪聲的本征模分量IMF_j；其中M為自然數(shù)，且1<M≤k；j為自然數(shù)，且1≤j≤M-1；

步驟5.對基本控制后的數(shù)據(jù)與本征模分量IMF_j進(jìn)行時(shí)間序列相關(guān)性分析得到相關(guān)系數(shù)v_j，根據(jù)相關(guān)系數(shù)v_j計(jì)算賦權(quán)本征模分量NEWIMF_j和賦權(quán)后的趨勢項(xiàng)NEWr；

步驟6.根據(jù)賦權(quán)本征模分量NEWIMF_j和賦權(quán)后的趨勢項(xiàng)NEWr進(jìn)行數(shù)據(jù)重構(gòu)得到去噪聲后的數(shù)據(jù)，完成去噪質(zhì)量控制。

2.根據(jù)權(quán)利要求1所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的自動氣象站數(shù)據(jù)質(zhì)量控制方法，其特征在于：步驟2中所述基本質(zhì)量控制包括界限值檢查、氣候極值檢查、內(nèi)部一致性檢查、時(shí)變一致性檢查、空間一致性檢查、特殊天氣事件檢查。

3.根據(jù)權(quán)利要求1所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的自動氣象站數(shù)據(jù)質(zhì)量控制方法，其特征在于：步驟4中所述計(jì)算本征模分量的能量E_i運(yùn)用如下公式計(jì)算：

Ei=Σi=1kIMFi(t)2,]]>

其中IMF_i(t)為第i個(gè)本征模分量的時(shí)間序列，t為采樣時(shí)間點(diǎn)。

4.根據(jù)權(quán)利要求1所述的基于經(jīng)驗(yàn)?zāi)B(tài)分解的自動氣象站數(shù)據(jù)質(zhì)量控制方法，其特征在于：步驟5中所述計(jì)算賦權(quán)本征模分量運(yùn)用如下公式計(jì)算：

NEWIMF_j=v_j·IMF_j；j=1,2…M-1

所述賦權(quán)后的趨勢項(xiàng)運(yùn)用如下公式計(jì)算：

NEWr=Σp=MkIMFp+r.]]>