基于RBF 神經(jīng)網(wǎng)絡(luò)的新型AUV 自抗擾控制方法

烏云嘎，徐會(huì)希，姜志斌

(1.中國(guó)科學(xué)院沈陽(yáng)自動(dòng)化研究所機(jī)器人學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110016；2.中國(guó)科學(xué)院機(jī)器人與智能制造創(chuàng)新研究院，遼寧 沈陽(yáng) 110169；3.遼寧省水下機(jī)器人重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110169；4.中國(guó)科學(xué)院大學(xué)，北京 100049)

0 引 言

AUV 因其自主、智能的特性，在水下環(huán)境探測(cè)和資源勘探方面發(fā)揮了人工無(wú)法比擬的作用[1]。AUV 作業(yè)過(guò)程中，受到海浪和水流的干擾，所以在近海底自主觀測(cè)、軌跡跟蹤等要求精細(xì)控制的場(chǎng)景，必須考慮外界環(huán)境擾動(dòng)對(duì)AUV 穩(wěn)定性的影響。AUV 平臺(tái)具有多變量、非線性、強(qiáng)耦合、外界干擾多和內(nèi)部參數(shù)時(shí)變的特點(diǎn)，要想達(dá)到理想的控制效果，往往需要收斂速度快、抗干擾、魯棒能力強(qiáng)的控制器，這些要求給控制系統(tǒng)的設(shè)計(jì)帶來(lái)了很多困難。

基于對(duì)象信息的現(xiàn)代控制理論包含對(duì)模型不確定性的近似，這種近似性使得AUV 控制系統(tǒng)性能有待提高。針對(duì)控制器本身存在不足及其對(duì)擾動(dòng)抑制能力弱的問(wèn)題，韓京清等[3–5]提出了一種抗擾動(dòng)控制方法，對(duì)控制模型信息依賴少，具有天然的解耦性和優(yōu)秀的抗干擾能力，被廣泛應(yīng)用于機(jī)器人控制。針對(duì)非線性自抗擾控制參數(shù)比較多的問(wèn)題，高志強(qiáng)[6]提出了線性化版本，簡(jiǎn)化了工程實(shí)現(xiàn)和參數(shù)整定，更加便于工業(yè)推廣。

目前針對(duì)線性自抗擾控制器的理論研究主要集中在其時(shí)域和頻域特性的探索，以及與PID 的關(guān)系的證明上，總結(jié)出了很多控制器調(diào)參的指導(dǎo)方法。其中，JIN 等[7]證明二階線性ADRC 可以等效于增加了前饋和二階低通濾波器PID 控制，并提出了按照等效PID 參數(shù)對(duì)二階線性ADRC 調(diào)參的方法。崔文慶等[8]說(shuō)明在中低頻段線性ADRC 等效于PID 控制器串聯(lián)n階超前環(huán)節(jié)。除了根據(jù)頻域特性和類比PID 的調(diào)參方法，還有很多自適應(yīng)調(diào)參方法，如基于強(qiáng)化學(xué)習(xí)改進(jìn)的自抗擾控制算法在無(wú)人船中的應(yīng)用[9-10]。以上理論分析說(shuō)明ADRC 實(shí)際上是改進(jìn)版的PID，PID 是ADRC 的一種特殊情況。但是當(dāng)輸入階躍信號(hào)的幅值改變較大時(shí)，得到的過(guò)渡過(guò)程不同，原有的擴(kuò)張狀態(tài)觀測(cè)器對(duì)擾動(dòng)的觀測(cè)誤差將增大，進(jìn)而影響系統(tǒng)的穩(wěn)定性，需要重新調(diào)節(jié)控制器參數(shù)，才能取得好的控制效果。

為此，本文提出一種基于RBF 神經(jīng)網(wǎng)絡(luò)補(bǔ)償?shù)木€性自抗擾控制方法：利用神經(jīng)網(wǎng)絡(luò)的自學(xué)習(xí)和自適應(yīng)能力，將其用于擾動(dòng)觀測(cè)器的改進(jìn)，實(shí)現(xiàn)在不同階躍信號(hào)下對(duì)擾動(dòng)的準(zhǔn)確估計(jì)，以此改善ADRC 的控制性能，簡(jiǎn)化控制器調(diào)參過(guò)程，提高控制器的魯棒性。將改進(jìn)后的控制器應(yīng)用于AUV 的航向角控制，通過(guò)仿真實(shí)驗(yàn)說(shuō)明了方法的有效性，驗(yàn)證了所提方法對(duì)不同采樣頻率的魯棒性。仿真結(jié)果表明，在低采樣頻率下改進(jìn)后的控制方法有更好的控制效果，對(duì)實(shí)際工程具有參考意義。

1 AUV 模型

1.1 坐標(biāo)系、術(shù)語(yǔ)和符號(hào)定義

按照國(guó)際水池（ITTC）推薦的和造船與輪機(jī)工程學(xué)會(huì)（SNAME）術(shù)語(yǔ)公報(bào)的體系[11]，采用2 個(gè)右手坐標(biāo)系表述AUV 的運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)模型：一個(gè)是固定于地球上的定坐標(biāo)系E-ξηζ，原點(diǎn)取為水中任意一點(diǎn)；另一個(gè)是與AUV 固聯(lián)并隨之運(yùn)動(dòng)的動(dòng)坐標(biāo)系G-xyz，原點(diǎn)取為AUV 的浮心。AUV 的位置和姿態(tài)分別用動(dòng)系相對(duì)于定系的坐標(biāo)（ ξB， ηB， ζB）和姿態(tài)角（橫滾角 φ，縱傾角θ，航向角 ψ）表示，而AUV 運(yùn)動(dòng)的速度（u，v，w）、角速度（p，q，r）、所受力（X，Y，Z）和力矩（K，M，N）定義于動(dòng)坐標(biāo)系中，規(guī)定向量正方向沿對(duì)應(yīng)坐標(biāo)系正方向，所有旋轉(zhuǎn)變量采用右手螺旋規(guī)則確定方向。坐標(biāo)系和變量的定義如圖1 所示，各變量對(duì)應(yīng)的符號(hào)列于表1。

表1 變量對(duì)應(yīng)的符號(hào)Tab.1 Signs of parameters

圖1 坐標(biāo)系和變量示意Fig.1 Coordinate systems and parameters

1.2 AUV 運(yùn)動(dòng)學(xué)模型

將AUV 看做自由運(yùn)動(dòng)的剛體，根據(jù)剛體六自由度運(yùn)動(dòng)學(xué)方程，AUV 位姿和運(yùn)動(dòng)速度之間的轉(zhuǎn)換關(guān)系表達(dá)為：

式中：T表示速度矢量從動(dòng)系到定系的變換矩陣，R代表角速度矢量從動(dòng)系到定系的變換矩陣。

1.3 AUV 動(dòng)力學(xué)模型

將AUV 看做剛體，參考剛體動(dòng)力學(xué)公式推導(dǎo)AUV 動(dòng)力學(xué)模型。記v=[V?]T，根據(jù)剛體六自由度動(dòng)力學(xué)方程：

式中：MRB∈?6×6為 慣性矩陣；CRB∈?6×6為科氏力和向心力矩陣； τRB∈?6為剛體受到的所有外力和外力矩。對(duì)于水下航行器應(yīng)用環(huán)境，表達(dá)式[12]為：

式 中： τhyd=-MAv˙r-CA(vr)vr-D(vr)vr為 水 動(dòng) 力 項(xiàng)，τhs=-g(Θ)-g0為重力、浮力和載重產(chǎn)生的水靜力項(xiàng)，τdis表示來(lái)自水流、海浪等的環(huán)境干擾力、力矩，τ=為推進(jìn)器和槳、舵的操縱力。本文實(shí)驗(yàn)使用的模型為REMUS 100 AUV，其操縱力來(lái)自尾推進(jìn)器、垂直操縱舵和水平操縱舵。

1.4 AUV 航向角控制模型

由式（2）和式（3）得到AUV 狀態(tài)變量的加速度表達(dá)式為：

實(shí)際控制器按照解耦后的模型來(lái)設(shè)計(jì)，所以可以對(duì)模型進(jìn)行合理的簡(jiǎn)化。假設(shè)整個(gè)過(guò)程中推進(jìn)器的轉(zhuǎn)速保持不變，AUV 處于巡航狀態(tài)，當(dāng)控制航向角時(shí)，只有垂直舵動(dòng)作，水平舵微調(diào)用于保持縱傾角穩(wěn)定。基于上述假設(shè)，得到航向角控制模型：

2 RBF 神經(jīng)網(wǎng)絡(luò)改進(jìn)的自抗擾控制器

2.1 自抗擾控制技術(shù)

ADRC 將積分器串聯(lián)型作為標(biāo)準(zhǔn)型，將總擾動(dòng)擴(kuò)張為系統(tǒng)的一個(gè)狀態(tài)，并經(jīng)觀測(cè)器估計(jì)后補(bǔ)償?shù)娇刂戚斎胫屑右韵_(dá)到簡(jiǎn)化控制系統(tǒng)、提高控制抗干擾性的目的。

ADRC 由擴(kuò)張狀態(tài)觀測(cè)器（Extended State Observer，ESO）、跟蹤微分器（Tracking Differentiator，TD）和狀態(tài)誤差反饋控制律（State Error Feedback，SEF）組成，如圖2 所示。其中ESO 用于估計(jì)系統(tǒng)的狀態(tài)和總擾動(dòng)，生成補(bǔ)償量，TD 用于配置過(guò)渡過(guò)程和獲取微分信號(hào)，SEF 用于生成控制信號(hào)。

圖2 ADRC 控制器結(jié)構(gòu)Fig.2 Structure of ADRC controller

設(shè)二階對(duì)象：

式中：w為外部擾動(dòng)，f(y,y˙,w,t)為包含內(nèi)擾和外擾的總擾動(dòng)。取系統(tǒng)狀態(tài)變量：v1=y，v2=y˙，可將上式轉(zhuǎn)化為狀態(tài)方程：

將航向角控制模型除控制量一階項(xiàng)外的部分視為該通道的內(nèi)部擾動(dòng)，則式（5）變?yōu)椋?/p>

對(duì)比式（6），可對(duì)航向角控制應(yīng)用擴(kuò)張狀態(tài)觀測(cè)器，估計(jì)總擾動(dòng)并補(bǔ)償?shù)娇刂菩盘?hào)中。

2.1.1 跟蹤微分器

TD 的算法為：

其中fhan(x1,x2,r0,h0)為最速綜合函數(shù)，表達(dá)式為：

式中：r0，h0為 函數(shù)控制參量，r0稱為速度因子，根據(jù)過(guò)渡過(guò)程快慢的需要和系統(tǒng)的承受能力決定，h0稱為濾波因子，常取系統(tǒng)采樣步長(zhǎng)的若干倍。d=r0h20，a0=h0x2，x3=x1+a0。

TD 的狀態(tài)變量v1能 夠跟蹤輸入信號(hào)v0，配置過(guò)渡過(guò)程，生成快速、無(wú)超調(diào)的過(guò)渡信號(hào)，而v2作 為v0的近似微分能夠跟蹤輸入信號(hào)的微分。

2.1.2 擴(kuò)張狀態(tài)觀測(cè)器

ESO 是ADRC 的核心，其以控制信號(hào)δ和輸出信號(hào)y為輸入量，輸出對(duì)系統(tǒng)狀態(tài)和總擾動(dòng)的估計(jì)。

將總擾動(dòng)擴(kuò)張成新的狀態(tài)變量：

則可對(duì)系統(tǒng)式（7）建立線性擴(kuò)張狀態(tài)觀測(cè)器：

2.2 自適應(yīng)RBF 神經(jīng)網(wǎng)絡(luò)改進(jìn)的自抗擾控制器

根據(jù)擴(kuò)張狀態(tài)觀測(cè)器的表達(dá)式（13），觀測(cè)器帶寬是固定不變的，只對(duì)一定頻率范圍內(nèi)的干擾估計(jì)效果比較好，如果擾動(dòng)的頻率范圍較寬，觀測(cè)效果就會(huì)降低，需要重新調(diào)整觀測(cè)器帶寬。同時(shí)，由于AUV模型具有非線性、耦合、時(shí)變的特點(diǎn)，不同階躍輸入和采樣頻率也會(huì)引起內(nèi)部擾動(dòng)的變化，影響觀測(cè)器的效率。針對(duì)這一問(wèn)題，可以用RBF 神經(jīng)網(wǎng)絡(luò)改進(jìn)觀測(cè)器，使之具有自適應(yīng)特性。

RBF 神經(jīng)網(wǎng)絡(luò)是20 世紀(jì)80 年代末提出的一種神經(jīng)網(wǎng)絡(luò)，能以任意精度逼近任意連續(xù)函數(shù)，具有學(xué)習(xí)速度快、避免陷入局部最小值的特點(diǎn)，采用RBF 神經(jīng)網(wǎng)絡(luò)對(duì)系統(tǒng)式（7）中的總擾動(dòng)進(jìn)行自適應(yīng)逼近。

RBFNN 是一種3 層的神經(jīng)網(wǎng)絡(luò)，結(jié)構(gòu)如圖3 所示，包括輸入層、隱層和輸出層，輸入層到隱層之間為非線性映射，激活函數(shù)采用徑向基函數(shù)：

圖3 RBF 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.3 Structure of RBF Neural Network

徑向基函數(shù)取值僅與距中心的距離有關(guān)，且兩邊對(duì)稱。標(biāo)準(zhǔn)的徑向基函數(shù)一般取歐式距離，最常用的是高斯核函數(shù)。隱層到輸出層之間為線性映射，取徑向基函數(shù)輸出值的加權(quán)和，權(quán)值可以通過(guò)自適應(yīng)律調(diào)整，稱為自適應(yīng)徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)。

取RBF 神經(jīng)網(wǎng)絡(luò)的輸入為TD 設(shè)定值和ESO 觀測(cè)值的誤差：

采用高斯核RBF 神經(jīng)網(wǎng)絡(luò)的算法為：

此時(shí)可使神經(jīng)網(wǎng)絡(luò)的逼近誤差 σ最小，且控制系統(tǒng)穩(wěn)定。其中， α ≥0 為神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)率，b=[0 1]T，P為一個(gè)正定矩陣，滿足李雅普諾夫方程：

其中，Q是任意的 2×2正 定矩陣， Λ定義為：

式中，kp和kd對(duì)應(yīng)式（16）中取值。

神經(jīng)網(wǎng)絡(luò)的實(shí)際輸出為：

可以用f?(E) 代替z3補(bǔ)償控制信號(hào)，使抗擾動(dòng)有自適應(yīng)、自學(xué)習(xí)的性能。RBF 神經(jīng)網(wǎng)絡(luò)改進(jìn)后的自抗擾控制器ADRC+RBFNN 如圖4 所示。

圖4 基于RBFNN 改進(jìn)的ADRC 控制器結(jié)構(gòu)Fig.4 Structure of ADRC controller based on RBFNN

3 仿真實(shí)驗(yàn)

為了驗(yàn)證本文提出方法的有效性，設(shè)計(jì)仿真實(shí)驗(yàn)。針對(duì)AUV 姿態(tài)控制，分別建立PID、ADRC 和ADRC+RBFNN 控制器，進(jìn)行姿態(tài)控制仿真。在不同時(shí)刻輸入不同幅值的階躍信號(hào)，對(duì)比控制器的動(dòng)態(tài)響應(yīng)性能；然后設(shè)計(jì)控制端擾動(dòng)和反饋信號(hào)噪聲，研究控制器對(duì)干擾的觀測(cè)和補(bǔ)償能力。最后，比較采樣頻率對(duì)改進(jìn)前后控制器的影響，驗(yàn)證控制器的工程可實(shí)現(xiàn)性。

實(shí)驗(yàn)中的主要參數(shù)包括ADRC 的速度因子r0，最速綜合函數(shù)步長(zhǎng)h0，觀測(cè)器帶寬w0，控制器帶寬wc，補(bǔ)償因 子b0，PID 參數(shù)kp，Ti，Td，N，RBF 神經(jīng)網(wǎng) 絡(luò)中心矢量cij，基寬向量bj， 學(xué)習(xí)率 α。PID 控制和ADRC控制仿真參數(shù)如表2 所示。為了保證仿真結(jié)果的客觀性，PID 控制的參數(shù)按照文獻(xiàn)[8]中提出的等效PID 的參數(shù)整定，ADRC 按照文獻(xiàn)[16]中方法得到。

表2 PID 和ADRC 控制器參數(shù)設(shè)置[2]Tab.2 Parameters setting of PID and ADRC controllers

為了得到準(zhǔn)確的結(jié)果，仿真中改進(jìn)方法的參數(shù)r0、h0、w0、wc、b0保持與ADRC 控制器參數(shù)一致。為了簡(jiǎn)化調(diào)參過(guò)程，高斯函數(shù)中心矢量cij取為零向量，基寬向量bj取 為單位向量，學(xué)習(xí)率 α =1.0。RBFNN隱含層有5 個(gè)節(jié)點(diǎn)，即采用2–5–1 結(jié)構(gòu)。

3.1 不同階躍輸入姿態(tài)控制仿真實(shí)驗(yàn)

為了比較不同階躍輸入下控制器的響應(yīng)性能，設(shè)計(jì)了AUV 航向角跟蹤控制仿真實(shí)驗(yàn)。仿真時(shí)間總長(zhǎng)150 s，分別于0 s、20 s、40 s、90 s 時(shí)輸入幅值為1 0?、-23.5?、 3 0?、 - 40?的階躍輸入信號(hào)，輸入幅值根據(jù)控制超調(diào)量10%得到，以便比較。

3.1.1 不含外部擾動(dòng)時(shí)的仿真實(shí)驗(yàn)結(jié)果

仿真中采用的參數(shù)在階躍輸入為1 0?時(shí)整定，不含外部擾動(dòng)時(shí)的航向角跟蹤控制仿真實(shí)驗(yàn)的結(jié)果如圖5所示。ADRC 和ADRC+RBFNN 控制方法對(duì)v1和v2的估計(jì)誤差分別如圖6 和圖7 所示。

圖5 無(wú)外擾時(shí)航向角跟蹤控制仿真結(jié)果Fig.5 Simulation results of yaw trajectory tracking without external disturbances

圖7 不含外部擾動(dòng)時(shí) v2估計(jì)誤差Fig.7 Estimation error of v2 without external disturbances

根據(jù)仿真結(jié)果可以發(fā)現(xiàn)，基于ADRC 的方法可以取得與PID 近似的效果，且在20 s 之前對(duì)信號(hào)的觀測(cè)誤差比較小，說(shuō)明控制器此時(shí)可以準(zhǔn)確的觀測(cè)模型不確定部分和各通道耦合部分。

階躍輸入幅值和過(guò)渡過(guò)程的不同會(huì)影響ADRC 的控制效果，結(jié)果是當(dāng)階躍值增加到2 3.5?時(shí)，傳統(tǒng)ADRC控制存在13.3%的超調(diào)量，而ADRC+RBFNN 控制超調(diào)量很小，對(duì)v1和v2的 觀測(cè)誤差較小。當(dāng)階躍值為30?時(shí)，改進(jìn)后方法超調(diào)量為17.2%，觀測(cè)誤差不存在振蕩，而且能夠更快穩(wěn)定。當(dāng)階躍值為 40?時(shí)，改進(jìn)后方法的控制效果也具有明顯的優(yōu)勢(shì)。

仿真結(jié)果說(shuō)明ADRC+RBFNN 方法在不同階躍信號(hào)輸入時(shí)，對(duì)內(nèi)部擾動(dòng)和系統(tǒng)狀態(tài)估計(jì)具有自適應(yīng)能力，能夠保證更大范圍的穩(wěn)定性，簡(jiǎn)化了控制器參數(shù)整定過(guò)程。

3.1.2 含外部擾動(dòng)時(shí)的仿真實(shí)驗(yàn)結(jié)果

根據(jù)文獻(xiàn)[17]，外部擾動(dòng)w包括控制端擾動(dòng)和反饋信號(hào)噪聲2 種。本次實(shí)驗(yàn)中控制端擾動(dòng)是幅值為 6?，周期為10 s 的正弦信號(hào)，反饋噪聲是均值為 0.5?的高斯白噪聲信號(hào)。含外部擾動(dòng)時(shí)的航向角根蹤控制仿真實(shí)驗(yàn)結(jié)果如圖8 所示，改進(jìn)前后控制方法對(duì)v1和v2的估計(jì)誤差分別如圖9 和圖10 所示，圖11 為兩者對(duì)控制端擾動(dòng)的估計(jì)誤差對(duì)比。

圖8 含外部擾動(dòng)時(shí)航向角跟蹤控制仿真結(jié)果Fig.8 Simulation results of yaw trajectory tracking with external disturbances

圖9 含外部擾動(dòng)時(shí) v1估計(jì)誤差Fig.9 Estimation error of v1 with external disturbances

圖10 含外部擾動(dòng)時(shí) v2估計(jì)誤差Fig.10 Estimation error of v2 with external disturbances

圖11 航向角控制擾動(dòng)估計(jì)仿真結(jié)果Fig.11 Simulation results of disturbances estimated value of controllers

由圖8 可知，ADRC 控制相對(duì)PID 具有更好的抗擾動(dòng)效果，可以更快穩(wěn)定，但是隨著階躍值的增加，原始ADRC 控制發(fā)生一定程度的振蕩。圖9 和圖10說(shuō)明此時(shí)ADRC 對(duì)總擾動(dòng)的估計(jì)出現(xiàn)較大誤差，圖11說(shuō)明其對(duì)控制端擾動(dòng)的估計(jì)出現(xiàn)較大的誤差。反觀ADRC+RBFNN 方法，雖然也會(huì)受到階躍變化的影響，但是總體超調(diào)量小于原始ADRC 方法，且對(duì)擾動(dòng)和狀態(tài)量的估計(jì)誤差能夠迅速收斂。

仿真結(jié)果說(shuō)明，相比ADRC 控制，本文改進(jìn)方法能夠更好適應(yīng)不同階躍輸入的變化，對(duì)總擾動(dòng)的估計(jì)更加準(zhǔn)確，提高了控制器的穩(wěn)定性和抗干擾能力。

3.2 不同采樣頻率姿態(tài)控制仿真實(shí)驗(yàn)

由于采樣頻率會(huì)影響擴(kuò)張狀態(tài)觀測(cè)器的效果，而AUV 屬于大慣性系統(tǒng)，實(shí)際平臺(tái)中的采樣頻率比較低，有明顯的離散特性。為了比較不同采樣頻率對(duì)改進(jìn)前后控制方法的影響，設(shè)計(jì)仿真實(shí)驗(yàn)。選擇采樣頻率分別為50 Hz、20 Hz、10 Hz，對(duì)比控制器的效果，仿真結(jié)果如圖12 所示。

圖12 不同采樣頻率航向角跟蹤控制仿真結(jié)果Fig.12 Simulation results of yaw trajectory tracking control at different sampling frequencies

仿真結(jié)果說(shuō)明改進(jìn)后的ADRC+RBFNN 控制器在低采樣頻率下依然能夠取得理想的控制效果，進(jìn)而可以通過(guò)嵌入式系統(tǒng)硬件和工業(yè)軟件組件的方式在實(shí)際AUV 平臺(tái)中實(shí)現(xiàn)，具有很大的工程應(yīng)用潛力。

4 結(jié) 語(yǔ)

針對(duì)AUV 作業(yè)過(guò)程中存在的環(huán)境干擾，為了提高在不同階躍輸入和采樣頻率下的擾動(dòng)觀測(cè)能力和抗干擾能力，本文提出一種基于RBF 神經(jīng)網(wǎng)絡(luò)補(bǔ)償?shù)男滦妥钥箶_控制方法。在線性ADRC 控制的基礎(chǔ)上，將自適應(yīng)RBF 神經(jīng)網(wǎng)絡(luò)應(yīng)用于改進(jìn)擴(kuò)張狀態(tài)觀測(cè)器，提升了擾動(dòng)估計(jì)的快速性、準(zhǔn)確性。仿真實(shí)驗(yàn)證明，改進(jìn)后的ADRC 在不同階躍輸入下對(duì)擾動(dòng)和系統(tǒng)狀態(tài)的觀測(cè)更加準(zhǔn)確，控制器跟蹤航向角能夠更快達(dá)到穩(wěn)定值，且適用于低采樣頻率。

后續(xù)工作將重點(diǎn)研究本文方法在AUV 平臺(tái)中的應(yīng)用問(wèn)題，通過(guò)水池和開(kāi)放水域?qū)嶒?yàn)驗(yàn)證本文提出方法的有效性。