陳志高，柳銳，吳子豪，陳志平，吳自銀

（1.東華理工大學(xué) 自然資源部環(huán)鄱陽(yáng)湖區(qū)域礦山環(huán)境監(jiān)測(cè)與治理重點(diǎn)實(shí)驗(yàn)室，江西 南昌 330013；2.東華理工大學(xué) 測(cè)繪與空間信息工程學(xué)院，江西 南昌 330013；3.自然資源部海底科學(xué)重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 310012）

流速測(cè)量是開(kāi)發(fā)和利用水（海洋）資源、開(kāi)展水利（海洋）工程建設(shè)所面臨的基礎(chǔ)問(wèn)題之一［1］。聲學(xué)多普勒流速剖面儀（acoustic Doppler current profiler，ADCP）借助水跟蹤（water tracking，WT）獲得的相對(duì)流速剖面加上底跟蹤（bottom track‐ing，BT）獲得的船速即可實(shí)時(shí)獲取三維流速剖面，具有速度快、精度高、不干擾流場(chǎng)等優(yōu)勢(shì)，被國(guó)際海委會(huì)定為4 種先進(jìn)的海洋觀測(cè)儀器之一［2］，在流速觀測(cè)和流量測(cè)驗(yàn)中被廣泛采用［3-6］。然而，受復(fù)雜測(cè)量環(huán)境影響，ADCP 在實(shí)際應(yīng)用中仍存在底質(zhì)流動(dòng)或懸沙濃度較高時(shí)底跟蹤船速失效，進(jìn)而直接導(dǎo)致流速測(cè)量結(jié)果出現(xiàn)偏差或者數(shù)據(jù)丟失等問(wèn)題，嚴(yán)重影響了測(cè)量結(jié)果的可靠性，限制了ADCP 的作業(yè)環(huán)境。

針對(duì)底質(zhì)流動(dòng)或高濃度懸沙情況下ADCP 底跟蹤船速失效的情況，國(guó)內(nèi)外學(xué)者一般采用全球定位系統(tǒng)（global positioning system，GPS）前后歷元的GGA 定位信息計(jì)算出速度，并對(duì)底跟蹤船速進(jìn)行替換［7-9］。研究結(jié)果表明，采用GPS 船速替換方法可以較好地解決底跟蹤船速失效的問(wèn)題。然而，一方面，受岸邊樹(shù)木遮擋導(dǎo)致可用衛(wèi)星數(shù)目較少、風(fēng)浪導(dǎo)致船體姿態(tài)劇烈變化等影響，GPS 定位解易出現(xiàn)異常［10］，這將導(dǎo)致GPS 船速精度降低，進(jìn)而引起流速計(jì)算出現(xiàn)偏差；另一方面，底跟蹤船速精度通常優(yōu)于1 cm/s，高于利用GPS 定位解計(jì)算的船速精度［11］，直接的船速替換將有損底跟蹤數(shù)據(jù)有效時(shí)絕對(duì)流速的測(cè)量精度。因此，采用船速直接替換未能根本解決ADCP 船速基準(zhǔn)的準(zhǔn)確獲取問(wèn)題。

卡爾曼（Kalman）濾波在航跡估計(jì)、姿態(tài)融合解算等方面廣泛使用［12-13］。Kalman 濾波模型包括狀態(tài)方程和觀測(cè)方程，其核心思想為預(yù)測(cè)+測(cè)量反饋。在Kalman 濾波計(jì)算過(guò)程中，誤差協(xié)方差陣對(duì)濾波系統(tǒng)的精度起著重要的作用。因此，系統(tǒng)噪聲協(xié)方差陣Q和觀測(cè)誤差協(xié)方差R的確定尤為關(guān)鍵。在實(shí)際的測(cè)量過(guò)程中，噪聲的隨機(jī)特性主要是受測(cè)量?jī)x器以及周圍環(huán)境的影響，誤差協(xié)方差陣Q和觀測(cè)誤差協(xié)方差陣R不可能先驗(yàn)已知，因而單一固定的先驗(yàn)方差可能會(huì)嚴(yán)重地影響濾波精度。在針對(duì)不確定的誤差協(xié)方差陣Q和R的算法中，Sage-Husa 算法具有一定的代表性［14］。該算法依據(jù)最小均方誤差準(zhǔn)則，利用觀測(cè)數(shù)據(jù)在遞推求解的同時(shí)，通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器，實(shí)時(shí)估計(jì)和修正系統(tǒng)噪聲和量測(cè)噪聲的統(tǒng)計(jì)特性參數(shù)。該算法基本忽略了初值Q和R，這有助于實(shí)時(shí)估計(jì)噪聲的統(tǒng)計(jì)特性。動(dòng)態(tài)模型中的不確定性會(huì)影響濾波過(guò)程，并且GPS 異常解及底跟蹤失效等異常數(shù)據(jù)將直接影響船速濾波結(jié)果，進(jìn)而導(dǎo)致流量計(jì)算結(jié)果產(chǎn)生偏差。同時(shí)應(yīng)用Sage-Husa 算法時(shí)，容易出現(xiàn)R失去正定性或Q失去非負(fù)定性的問(wèn)題，使濾波發(fā)散。

基于上述問(wèn)題，本文結(jié)合自適應(yīng)抗差濾波原理提出一種改進(jìn)的Sage-Husa 自適應(yīng)濾波算法。通過(guò)狀態(tài)不符值統(tǒng)計(jì)量構(gòu)造自適應(yīng)因子ak并基于丹麥法構(gòu)建抗差等價(jià)權(quán)因子w，二者可分別抑制動(dòng)力學(xué)模型擾動(dòng)異常（波浪、潮汐對(duì)船速的擾動(dòng)）和觀測(cè)信息異常（GPS 和底跟蹤觀測(cè)異常值），提高濾波的自適應(yīng)性和魯棒性；為防止在更新時(shí)出現(xiàn)R失去正定性或Q失去非負(fù)定性導(dǎo)致濾波發(fā)散的問(wèn)題，對(duì)Q和R增加Kalman 濾波的約束條件。通過(guò)以上2 步，實(shí)現(xiàn)底跟蹤船速和GPS 船速的有機(jī)融合：既可在底質(zhì)流動(dòng)或高濃度懸沙情況下獲取穩(wěn)健的GPS 船速，又可保留正常情況下高精度的底跟蹤船速，進(jìn)而增強(qiáng)ADCP 船速基準(zhǔn)的準(zhǔn)確性與穩(wěn)定性，提高ADCP 走航觀測(cè)的普適性。

1 濾波模型

1.1 Kalman濾波

由ADCP 底跟蹤船速VBT與外部傳感器GPS 獲得的船速VGPS組成的觀測(cè)值矩陣可表示為：

式中：k為時(shí)間步長(zhǎng)；T 為矩陣轉(zhuǎn)置；下標(biāo)E和N分別表示東分量和北分量。

狀態(tài)向量為：

式中VE、VN分別是船速的東分量和北分量。

Kalman 濾波模型的觀測(cè)方程和狀態(tài)方程可以表示為：

式中：Hk為觀測(cè)矩陣；Φk為狀態(tài)轉(zhuǎn)移矩陣。相應(yīng)的，系統(tǒng)噪聲協(xié)方差陣Q和量測(cè)噪聲協(xié)方差陣R分別為：

式中：SP為真實(shí)的船舶動(dòng)態(tài)位置的方差；σGPS、σBT為測(cè)量誤差。

然而，在實(shí)際情況中，噪聲的隨機(jī)特性主要是受測(cè)量?jī)x器以及周圍環(huán)境的影響，Q陣和R陣不可能先驗(yàn)已知，因而單一固定的先驗(yàn)方差可能會(huì)嚴(yán)重地影響濾波精度。

1.2 Sage-Husa算法

在衍生于Kalman 濾波的自適應(yīng)算法中，Sage-Husa 算法是針對(duì)不確定的Q陣和R陣的經(jīng)典算法。該算法依據(jù)最小均方誤差準(zhǔn)則，利用觀測(cè)數(shù)據(jù)在遞推求解的同時(shí)，通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器，實(shí)時(shí)估計(jì)和修正系統(tǒng)噪聲和量測(cè)噪聲的統(tǒng)計(jì)特性參數(shù)，達(dá)到降低模型誤差、抑制濾波發(fā)散、提高濾波精度的目的［15］。具體的調(diào)節(jié)方式為：

式中：vk為觀測(cè)殘差向量；dk為加權(quán)因子；b為遺忘因子，一般取0.95～0.99。

然而，雖然通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器可以自適應(yīng)更新Q和R，提高Kalman 濾波的性能，但是動(dòng)態(tài)模型中的不確定性會(huì)影響濾波過(guò)程；同時(shí)，應(yīng)用Sage-Husa 算法時(shí),容易出現(xiàn)R失去正定性或Q失去非負(fù)定性的問(wèn)題。

1.3 改進(jìn)Sage-Husa算法

針對(duì)上述Sage-Husa算法的不足，本文根據(jù)自適應(yīng)抗差濾波理論，一方面，構(gòu)建自適應(yīng)因子ak和抗差等價(jià)權(quán)因子w，抑制動(dòng)力學(xué)模型擾動(dòng)異常和觀測(cè)信息異常，以提高濾波的自適應(yīng)性和魯棒性；另一方面，對(duì)Q和R矩陣增加約束條件，以保證Sage-Husa算法的穩(wěn)定性。本文改進(jìn)算法的基本流程如圖1所示。

圖1 改進(jìn)的Sage-Husa算法流程Fig.1 Flow chart of improved Sage-Husa algorithm

1.3.1 基本流程

依然假設(shè)量測(cè)噪聲Vk與系統(tǒng)噪聲Wk是統(tǒng)計(jì)特性已知的零均值高斯白噪聲且互不相關(guān)，Vk和Wk均與初始狀態(tài)X0無(wú)關(guān)，改進(jìn)算法的遞推方程如下。

預(yù)測(cè)方程為：

預(yù)測(cè)誤差協(xié)方差陣為：

觀測(cè)殘差向量為：

濾波增益矩陣為：

濾波解為：

協(xié)方差陣更新為：

系統(tǒng)噪聲協(xié)方差陣Q陣和量測(cè)噪聲協(xié)方差陣R陣更新為：

式中：ak為自適應(yīng)因子；Pt為觀測(cè)信息的抗差等價(jià)權(quán)矩陣；w為抗差等價(jià)權(quán)因子，其對(duì)角元素wi可采用觀測(cè)殘差判別統(tǒng)計(jì)量并基于丹麥法進(jìn)行構(gòu)造［16］。

與常規(guī)Sage-Husa 算法不同，本文改進(jìn)算法根據(jù)自適應(yīng)抗差濾波原理，首先基于丹麥法構(gòu)造抗差因子w，計(jì)算觀測(cè)信息的抗差等價(jià)權(quán)矩陣Pt，當(dāng)觀測(cè)信息中含有粗差時(shí)，抗差估計(jì)通過(guò)抗差因子w淘汰粗差；然后，根據(jù)狀態(tài)不符值統(tǒng)計(jì)量，基于三段函數(shù)法構(gòu)造自適應(yīng)因子ak，自適應(yīng)因子有助于減少預(yù)測(cè)狀態(tài)誤差模型中的不確定性［17］；最后，為了防止在更新時(shí)，Rk+1失去正定性或Qk+1失去非負(fù)定性的問(wèn)題，對(duì)Qk+1和Rk+1增加卡爾曼濾波的約束條件，將式（14）和（15）的等式右邊的第2 項(xiàng)處理結(jié)果的對(duì)角線元素取絕對(duì)值，非對(duì)角線元素取零［14］。

本文重點(diǎn)介紹自適應(yīng)因子ak和抗差等價(jià)權(quán)因子w的具體構(gòu)造過(guò)程。

1.3.2 自適應(yīng)因子ak的構(gòu)造

對(duì)于自適應(yīng)因子的構(gòu)造，楊元喜等［18］給出了4種構(gòu)造自適應(yīng)因子ak的誤差判別統(tǒng)計(jì)量，包括狀態(tài)不符值統(tǒng)計(jì)量、預(yù)測(cè)殘差統(tǒng)計(jì)量、方差分量比統(tǒng)計(jì)量和速度不符值統(tǒng)計(jì)量。為避免觀測(cè)信息不可靠時(shí)對(duì)ak構(gòu)造的影響，本文選擇基于狀態(tài)不符值統(tǒng)計(jì)量并采用三段函數(shù)法構(gòu)建ak，表示為［19-20］：

1.3.3 抗差等價(jià)權(quán)因子w的構(gòu)造

ADCP底跟蹤船速失效及GPS定位解易出現(xiàn)異常均會(huì)直接導(dǎo)致流速測(cè)量結(jié)果出現(xiàn)偏差，進(jìn)而影響流量計(jì)算結(jié)果的可靠性。因此需要運(yùn)用抗差估計(jì)原理抵制粗差干擾。丹麥法構(gòu)建抗差等價(jià)權(quán)因子可避免發(fā)生矩陣奇異，當(dāng)殘差的可信度處在模糊界限時(shí)，該方法能夠較有效地抑制粗差對(duì)濾波精度的影響。丹麥法抗差等價(jià)權(quán)函數(shù)為［21］：

式中：wi為第i個(gè)觀測(cè)向量的等價(jià)權(quán)因子；vi為第i個(gè)觀測(cè)向量的標(biāo)準(zhǔn)化殘差；v（ki）和分別為殘差的第i個(gè)分量和其標(biāo)準(zhǔn)差矩陣的第i個(gè)分量；c2一般取值為1.5～2.0。

2 船速濾波實(shí)驗(yàn)結(jié)果及分析

2.1 實(shí)驗(yàn)數(shù)據(jù)

為驗(yàn)證各種濾波方法的有效性，在長(zhǎng)江口水域（圖2）開(kāi)展了走航式ADCP 流速測(cè)量實(shí)驗(yàn)。實(shí)驗(yàn)水域?yàn)槌毕佣?，底質(zhì)穩(wěn)定；流速測(cè)量采用RDI公司生產(chǎn)的300 kHz 瑞江系列4 波束ADCP，Leica SR530用于導(dǎo)航定位，THALES 3011 GPS 羅經(jīng)為其提供外部方位（定向精度優(yōu)于0.5°），并采用實(shí)時(shí)動(dòng)態(tài)量測(cè)技術(shù)（real time kinematic，RTK）獲取船速，其精度優(yōu)于±0.05 m/s。實(shí)驗(yàn)設(shè)計(jì)的走航斷面CD如圖2所示，其長(zhǎng)度約為1 200 m；對(duì)CD 斷面進(jìn)行往返觀測(cè)，共2個(gè)測(cè)回，4個(gè)測(cè)次。

圖2 位于長(zhǎng)江口的實(shí)驗(yàn)區(qū)域Fig.2 Located in the experimental area of the Yangtze River estuary

為測(cè)試不同底質(zhì)流動(dòng)速度對(duì)觀測(cè)流速及船速濾波結(jié)果的影響，實(shí)驗(yàn)將模擬的底質(zhì)流動(dòng)速度疊加到實(shí)測(cè)的底跟蹤船速中。模擬的底質(zhì)流動(dòng)速度包括移動(dòng)偏差（moving bed，MB）與隨機(jī)噪聲2 部分。移動(dòng)偏差的強(qiáng)度f(wàn)MB由以下函數(shù)確定［22-23］：

式中：t為觀測(cè)時(shí)刻；T為觀測(cè)總時(shí)間；P為每個(gè)觀測(cè)點(diǎn)相對(duì)于斷面中心點(diǎn)的位置；e 為自然常數(shù)；F為最大強(qiáng)度；fMB為每個(gè)歷元的強(qiáng)度。由式（20）可以生成一個(gè)在斷面邊緣附近為0，然后沿?cái)嗝婧桔E逐漸增大，直至斷面中心達(dá)到最大值的移動(dòng)偏差序列。最后，在每個(gè)歷元中加入N（0，5 cm/s）的高斯白噪聲（隨機(jī)噪聲），即可得到最終模擬的底質(zhì)流動(dòng)速度。如圖3 所示，實(shí)驗(yàn)共模擬了移動(dòng)偏差fMB為5 cm/s（MB5），10 cm/s（MB10）和15 cm/s（MB15）的3 種不同強(qiáng)度底質(zhì)流動(dòng)速度。

圖3 實(shí)測(cè)流速與模擬的3種底質(zhì)流動(dòng)速度Fig.3 Measured water velocity and simulated three kinds of moving bed velocity

2.2 抗差性能比較

由于測(cè)量船在走航過(guò)程中存在障礙物遮擋或GPS 衛(wèi)星星座幾何形狀較差的情況，因此GPS（RTK）定位信息可能存在粗差。如圖4 所示，為了展現(xiàn)不同的濾波方法的抗差性能，實(shí)驗(yàn)在fMB=5 cm/s的條件下，對(duì)第1測(cè)次的GPS數(shù)據(jù)在131～133 s、151～155 s 之間在東、北分量加入大小為0.5 m/s 的粗差。由圖4 可知，Kalman 濾波、Sage-Husa 算法受粗差影響較大，導(dǎo)致濾波后船速曲線產(chǎn)生明顯“隆起”，所以需要在濾波前事先將其剔除；而本文改進(jìn)算法得到的濾波結(jié)果受粗差影響很小，所以在濾波前可以不作濾波剔除工作。

圖4 加入粗差后不同濾波方法的抗差結(jié)果Fig.4 Results of different methods after adding gross er‐ror

2.3 濾波精度分析

圖5 為在模擬移動(dòng)偏差fMB為5 cm/s 的條件下，由標(biāo)準(zhǔn)Kalman 濾波、Sage-Husa 算法和本文改進(jìn)Sage-Husa 算法這3 種濾波方法獲得的不同測(cè)次的船速結(jié)果。由圖5 可知，與Kalman 濾波和Sage-Hu‐sa 算法所得結(jié)果相比，本文方法的濾波結(jié)果更接近用于參考的GPS 船速，表明在存在底質(zhì)流動(dòng)的情況下，采用本文改進(jìn)濾波方法可以提高船速精度。

圖5 不同濾波方法得到的船速Fig.5 Ship velocity obtained by different filtering methods

圖6 所示為在fMB為5 cm/s 的條件下，由3 種濾波方法（下標(biāo)AKF表示Kalman濾波，下標(biāo)SHA表示Sage-Husa 算法，下標(biāo)ISH 表示本文改進(jìn)方法）獲得的船速與GPS 船速之間的東、北分量的均方根誤差（root mean squared error，RMSE）計(jì)算公式為：

圖6 不同濾波方法的船速東、北分量誤差Fig.6 Error of ship velocity east and north components with different filtering methods

式中：Vfit為濾波后的船速；Vdata為實(shí)測(cè)的GPS 船速；n為該測(cè)次的總歷元。

從圖6 可以看出，本文方法與Kalman 濾波和Sage-Husa算法相比，其東、北分量誤差曲線更接近0；其RMSE 在4 個(gè)測(cè)次取均值，可以求得本文方法在東分量的RMSE 為7.6 cm/s，而Kalman 濾波和Sage-Husa 算法分別為10.0 cm/s 和9.0 cm/s，分別相對(duì)提高了2.4 cm/s 和1.4 cm/s；本文方法在北分量的RMSE 為8.6 cm/s，而Kalman 濾波和Sage-Husa算法分別為12.1 cm/s和9.3 cm/s，分別提高了3.5 cm/s和0.7 cm/s。

為進(jìn)一步驗(yàn)證不同底質(zhì)流動(dòng)情況下3 種濾波方法的效果，實(shí)驗(yàn)采用矢量標(biāo)準(zhǔn)差σ和均方根誤差RMSE 作為精度指標(biāo)進(jìn)行評(píng)價(jià)。其中，，σe和σn分別為東和北分量的標(biāo)準(zhǔn)差；RMSE計(jì)算公式與式（21）一致。表1和表2分別為3種濾波方法在不同底質(zhì)流動(dòng)下的矢量標(biāo)準(zhǔn)差σ和均方根誤差RMSE計(jì)算結(jié)果。由表1和表2可知：

表1 3種濾波方法下的船速矢量標(biāo)準(zhǔn)差σTable1 Standard deviation σ of boat velocity under 3 filtering methods cm/s

表2 3種濾波方法下的船速RMSETable 2 Root mean square error of velocity under 3 filtering methods cm/s

1）Sage-Husa 算法在底質(zhì)流動(dòng)MB5、MB10、MB15下的σ對(duì)比Kalman濾波分別提高了1.6 cm/s、1.4 cm/s 和1.2 cm/s，RMSE 則分別提高了1.1 cm/s、1.0 cm/s 和0.7 cm/s。究其原因，受波浪、風(fēng)、泥沙濃度等影響ADCP 走航測(cè)量環(huán)境較為復(fù)雜，單一固定的Q陣和R陣會(huì)影響濾波效果，而Sage-Husa算法通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器可以自適應(yīng)更新Q和R，提高Kalman濾波的性能。

2）本文改進(jìn)方法在底質(zhì)流動(dòng)MB5、MB10、MB15下的σ相對(duì)Sage-Husa算法又分別提高了1.0 cm/s、1.2 cm/s 和1.5 cm/s，RMSE 則分別提高了1.0 cm/s、1.0 cm/s 和1.3 cm/s。本文改進(jìn)方法相對(duì)Sage-Husa算法精度得到了進(jìn)一步提升，主要原因是本文改進(jìn)方法不僅通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器更新Q和R，還通過(guò)構(gòu)造自適應(yīng)因子控制動(dòng)力學(xué)模型誤差的影響。

2.4 流量計(jì)算

由ADCP 流量測(cè)量原理可知，走航式ADCP 測(cè)量所得流量是由各個(gè)ping（歷元）獲得的微斷面流量累加得到，斷面總流量計(jì)算公式為［24］：

式中：Qm為實(shí)測(cè)流量；Vx/y，WT為ADCP 水跟蹤獲得的相對(duì)流速的東/北分量；Vx y,BV為GPS 或底跟蹤獲得的船速的東/北分量；N為微斷面總個(gè)數(shù)；i為第微斷面序號(hào)；hi為第i個(gè)微斷面的深度；Δti為相鄰微斷面i和i-1之間的采樣時(shí)間間隔。

圖7所示為模擬底質(zhì)流動(dòng)速度MB為5 cm/s的條件下，利用實(shí)測(cè)GPS 船速與3 種濾波方法得到的船速計(jì)算得到的絕對(duì)流速及其對(duì)應(yīng)的流量計(jì)算結(jié)果。圖7 的各個(gè)子圖中，從左至右分別為以實(shí)測(cè)GPS 速度、傳統(tǒng)Kalman濾波、Sage-Husa濾波以及本文改進(jìn)方法濾波后的速度作為船速參考計(jì)算的絕對(duì)流速，相應(yīng)的流量計(jì)算結(jié)果分別記為QGPS、QAKF、QSHA及QISH。流量單位為m3/s，負(fù)號(hào)表示漲潮，正號(hào)表示落潮。為便于比較，表3 還列出了fMB為10 cm/s 和15 cm/s 的條件下，由GPS 船速與3 種濾波方法船速作為參考計(jì)算出來(lái)的流量。由圖7及表3可得出以下結(jié)論：

表3 3種濾波方法下計(jì)算的斷面流量Table 3 Section discharge calculation under three filtering methods m3·s-1

圖7 不同濾波方法的流量計(jì)算結(jié)果Fig.7 Discharge calculation of different filtering methods

1）隨著底質(zhì)流動(dòng)速度f(wàn)MB的增大，計(jì)算得到的斷面流量結(jié)果變小。究其原因，由于底質(zhì)流動(dòng)速度與實(shí)際流速方向基本一致，受底質(zhì)流動(dòng)影響，實(shí)測(cè)的底跟蹤船速產(chǎn)生和底質(zhì)流動(dòng)速度一樣的矢量偏差，而計(jì)算流量的式（22）實(shí)際上是由船速矢量和流速矢量相乘轉(zhuǎn)化而來(lái)，因此該公式計(jì)算的流量結(jié)果必然偏?。?/p>

2）隨著底質(zhì)流動(dòng)速度的增大，本文改進(jìn)算法相對(duì)于傳統(tǒng)Kalman 濾波方法其流量計(jì)算精度明顯提高。在fMB為5 cm/s的條件下，本文方法的精度相對(duì)于Kalman 濾波和Sage-Husa 算法分別提升了0.9%和0.6%，在fMB為10 cm/s 的條件下分別提升了3.0%和0.4%，在fMB為15 cm/s 的條件下則分別提升了6.0%和1.2%。盡管本文方法對(duì)船速精度的提高并不明顯（約為2 cm/s），但對(duì)于年均流量為2.8×104m3/s的長(zhǎng)江口這種大型感潮河段來(lái)說(shuō)，本文方法1%～6%的流量相對(duì)精度的提高顯得必要且有意義。

3 結(jié)論

1）傳統(tǒng)Kalman濾波結(jié)果受粗差影響較大，需事先將粗差剔除；本文改進(jìn)算法引入的抗差因子具有較強(qiáng)的抗差性，濾波前可不作粗差剔除工作。

2）本文改進(jìn)算法相較Kalman 濾波方法和Sage-Husa算法，其船速濾波精度提高了1.0～2.7 cm/s。主要原因是本文改進(jìn)方法不僅通過(guò)時(shí)變?cè)肼暯y(tǒng)計(jì)估值器更新誤差協(xié)方差陣Q和R，還通過(guò)構(gòu)造自適應(yīng)因子控制動(dòng)力學(xué)模型誤差的影響。

3）本文方法的流量計(jì)算精度相對(duì)于傳統(tǒng)Kal‐man濾波和Sage-Husa算法在底質(zhì)流動(dòng)速度為5 cm/s、10 cm/s 與15 cm/s 的情況下分別提升了0.9%～3.0%和0.4%～1.2%。盡管本文方法對(duì)船速精度的提高并不明顯，但對(duì)于年均流量為2.8×104m3/s的長(zhǎng)江口這種大型感潮河段來(lái)說(shuō)，本文方法1%～6%的流量相對(duì)精度的提高顯得必要且有意義。