王 飛, 祝汝松, 張俊生

(1. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所, 四川 綿陽(yáng) 621000；2. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國(guó)家重點(diǎn)試驗(yàn)室, 四川 綿陽(yáng) 621000)

卡爾曼濾波在某跨聲速風(fēng)洞CTS機(jī)構(gòu)連續(xù)動(dòng)態(tài)軌跡捕獲試驗(yàn)技術(shù)中的應(yīng)用

王 飛1,2, 祝汝松1,2, 張俊生1,2

(1. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所, 四川 綿陽(yáng) 621000；2. 中國(guó)空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國(guó)家重點(diǎn)試驗(yàn)室, 四川 綿陽(yáng) 621000)

外掛物分離軌跡捕獲技術(shù)是世界上一種先進(jìn)而且難度較大的風(fēng)洞試驗(yàn)技術(shù)，國(guó)內(nèi)某跨聲速風(fēng)洞捕獲軌跡試驗(yàn)裝置(CTS)首次對(duì)連續(xù)動(dòng)態(tài)軌跡捕獲實(shí)驗(yàn)技術(shù)進(jìn)行了研究。本文首先介紹了該跨聲速風(fēng)洞CTS機(jī)構(gòu)動(dòng)態(tài)軌跡的預(yù)測(cè)方法，隨后對(duì)該機(jī)構(gòu)連續(xù)動(dòng)態(tài)軌跡同步控制方法的具體實(shí)現(xiàn)步驟進(jìn)行了較為詳細(xì)的描述，分別介紹了無(wú)緩沖和有緩沖的動(dòng)態(tài)軌跡連續(xù)同步控制過(guò)程及約束條件；在此基礎(chǔ)上，針對(duì)CTS機(jī)構(gòu)天平信號(hào)硬件濾波滯后和精度要求間矛盾的現(xiàn)實(shí)，引入泰勒中值定理，建立卡爾曼濾波(KF)模型，選取適當(dāng)參數(shù)，對(duì)天平信號(hào)進(jìn)行動(dòng)態(tài)實(shí)時(shí)濾波處理，并與傳統(tǒng)的位置控制方式進(jìn)行吹風(fēng)試驗(yàn)比較。試驗(yàn)結(jié)果表明：卡爾曼濾波較好地平衡了CTS機(jī)構(gòu)天平信號(hào)硬件濾波滯后和精度要求之間的矛盾，且分離軌跡的重復(fù)性得到了驗(yàn)證，機(jī)構(gòu)風(fēng)洞吹風(fēng)試驗(yàn)的效率也大為提高。

跨聲速風(fēng)洞；捕獲軌跡試驗(yàn)裝置；連續(xù)；動(dòng)態(tài)；卡爾曼濾波

0 引 言

外掛物從母機(jī)上分離(發(fā)射/投放)的軌跡特性模擬是一項(xiàng)非常重要的試驗(yàn)技術(shù)，對(duì)母機(jī)飛行安全性評(píng)估有著極其重要的意義，同時(shí)也有助于發(fā)射/投放類航空武器的研究和發(fā)展。

國(guó)外在20世紀(jì)60、70年代開(kāi)始對(duì)外掛物分離軌跡模擬試驗(yàn)設(shè)備(多自由度捕獲軌跡試驗(yàn)裝置，即CTS機(jī)構(gòu))及相關(guān)試驗(yàn)技術(shù)的進(jìn)行設(shè)計(jì)和研究，用于預(yù)測(cè)/模擬外掛物從母機(jī)上分離(發(fā)射/投放)的軌跡特性，這些特性包括外掛物的線位移、角位移、線速度、角速度、線加速度和角加速度以及氣動(dòng)載荷等隨分離時(shí)間的變化。CTS機(jī)構(gòu)運(yùn)行控制流程主要分為兩步：①分離軌跡預(yù)測(cè)生成；②使用適度控制策略和算法控制機(jī)構(gòu)按生成的軌跡運(yùn)動(dòng)。其運(yùn)動(dòng)控制具有高復(fù)雜性、多變量強(qiáng)耦合等特點(diǎn)，因此最初的軌跡預(yù)測(cè)/模擬是基于位置控制的，是一種“走-停-走”的控制過(guò)程，一條軌的預(yù)測(cè)/模擬運(yùn)行時(shí)間較長(zhǎng)，導(dǎo)致CTS機(jī)構(gòu)風(fēng)洞試驗(yàn)效率低下。為此，70年代末，美國(guó)Arnold工程發(fā)展中心成功實(shí)現(xiàn)了基于速度控制的連續(xù)動(dòng)態(tài)捕獲軌跡試驗(yàn)技術(shù)，從而使得同一軌跡的預(yù)測(cè)/模擬運(yùn)行時(shí)間大大縮短，提高了CTS機(jī)構(gòu)風(fēng)洞試驗(yàn)效率。

國(guó)內(nèi)在外掛物分離捕獲軌跡試驗(yàn)技術(shù)領(lǐng)域的研究起步較晚，始于上世紀(jì)80、90年代，和國(guó)外的研究歷程相似，起初的CTS機(jī)構(gòu)運(yùn)行模式均基于位置控制。圖1為國(guó)內(nèi)某跨聲速風(fēng)洞CTS機(jī)構(gòu)示意圖，它首次對(duì)連續(xù)動(dòng)態(tài)軌跡捕獲試驗(yàn)技術(shù)進(jìn)行了系統(tǒng)性研究，既實(shí)現(xiàn)了位置控制，又實(shí)現(xiàn)了連續(xù)動(dòng)態(tài)控制。

圖1 國(guó)內(nèi)某跨聲速風(fēng)洞CTS機(jī)構(gòu)示意圖

在進(jìn)行上述跨聲速風(fēng)洞CTS機(jī)構(gòu)吹風(fēng)試驗(yàn)時(shí)，筆者發(fā)現(xiàn)現(xiàn)有的外掛物模型天平信號(hào)的硬件濾波對(duì)軌跡預(yù)測(cè)的準(zhǔn)確性和機(jī)構(gòu)運(yùn)行的實(shí)時(shí)性有一定的影響，同時(shí)這兩者又是矛盾的，為此在適當(dāng)提高硬件低通濾波截止頻率的基礎(chǔ)上，引入卡爾曼軟件濾波算法，以期改善矛盾，提高吹風(fēng)效率。

1 動(dòng)態(tài)軌跡預(yù)測(cè)

分離軌跡的預(yù)測(cè)生成是捕獲軌跡試驗(yàn)技術(shù)研究的基礎(chǔ)，圖2給出了CTS機(jī)構(gòu)控制原理。

圖2 CTS機(jī)構(gòu)控制原理圖

采集CTS機(jī)構(gòu)當(dāng)前軌跡點(diǎn)位置信息、六分量天平載荷信號(hào),綜合其它相關(guān)信息，通過(guò)適度算法，預(yù)測(cè)目標(biāo)位置值，由單軸(共6個(gè)軸)運(yùn)動(dòng)控制單元驅(qū)動(dòng)機(jī)構(gòu)各自由度同步運(yùn)行至目標(biāo)位置。

設(shè)f(k)為天平當(dāng)前軌跡點(diǎn)載荷信號(hào)，包含fx(k)、fy(k)、fz(k)、fα(k)、fβ(k)、fγ(k)6個(gè)分量；L(k)為機(jī)構(gòu)當(dāng)前軌跡點(diǎn)位置，包含x(k)、y(k)、z(k)、α(k)、β(k)、γ(k)6個(gè)分量；Vl(k+1)為目標(biāo)軌跡點(diǎn)理論速度，包含Vlx(k+1)、Vly(k+1)、Vlz(k+1)、Vlα(k+1)、Vlβ(k+1)、Vlγ(k+1)6個(gè)分量；L(k+1)為機(jī)構(gòu)目標(biāo)軌跡點(diǎn)的預(yù)測(cè)位置，包含x(k+1)、y(k+1)、z(k+1)、α(k+1)、β(k+1)、γ(k+1)6個(gè)分量，則可獲得機(jī)構(gòu)軌跡預(yù)測(cè)方程：

(1)

式(1)描述為：機(jī)構(gòu)目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置為當(dāng)前軌跡點(diǎn)位置、天平載荷及目標(biāo)軌跡點(diǎn)理論速度的函數(shù)方程，其函數(shù)法則Gk根據(jù)空氣動(dòng)力學(xué)推導(dǎo)獲得。

2 連續(xù)動(dòng)態(tài)軌跡同步控制方法

無(wú)論是最初的基于位置控制的捕獲軌跡試驗(yàn)技術(shù)，還是后來(lái)的基于速度控制的連續(xù)動(dòng)態(tài)捕獲軌跡試驗(yàn)技術(shù)，其控制流程均為“生成目標(biāo)軌跡點(diǎn)→機(jī)構(gòu)運(yùn)行至目標(biāo)軌跡點(diǎn)→生成下一目標(biāo)軌跡點(diǎn)→機(jī)構(gòu)運(yùn)行至下一目標(biāo)軌跡點(diǎn)→……”。

傳統(tǒng)的位置控制方法，機(jī)構(gòu)運(yùn)行至每一目標(biāo)軌跡點(diǎn)時(shí)速度均為零，在等待下一目標(biāo)軌跡點(diǎn)生成后，再運(yùn)行至下一目標(biāo)軌跡點(diǎn)，其運(yùn)動(dòng)控制流程是“走-停-走-……”的模式，一條軌跡的運(yùn)行時(shí)間較長(zhǎng)，效率低下；而連續(xù)動(dòng)態(tài)控制方法，機(jī)構(gòu)是以連續(xù)平滑的方式通過(guò)每一目標(biāo)軌跡點(diǎn)，在每一目標(biāo)軌跡點(diǎn)位置均不需等待，同一軌跡的運(yùn)行時(shí)間大大縮短，效率大為提高。

國(guó)內(nèi)某跨聲速風(fēng)洞CTS機(jī)構(gòu)連續(xù)動(dòng)態(tài)軌跡同步控制方法具體步驟描述如下：

(1)通過(guò)軌跡預(yù)測(cè)，實(shí)時(shí)獲得機(jī)構(gòu)目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置；

(2)以極小時(shí)間間隔ΔT為基準(zhǔn)，將當(dāng)前軌跡點(diǎn)位置至目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置的運(yùn)行軌跡分成N段；

(3)以每段實(shí)際起始位置及速度、目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置、基準(zhǔn)時(shí)間、速度極限、加速度極限等信息為輸入條件，推導(dǎo)出ΔT時(shí)刻后機(jī)構(gòu)目標(biāo)位置，即下一段起始位置；

(4)將ΔT時(shí)刻后機(jī)構(gòu)目標(biāo)位置作為定位目標(biāo)輸入運(yùn)動(dòng)控制系統(tǒng)，由運(yùn)動(dòng)控制系統(tǒng)完成機(jī)構(gòu)定位控制；

(5)重復(fù)上述步驟，保證機(jī)構(gòu)連續(xù)動(dòng)態(tài)地逼近目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置，直至該步長(zhǎng)運(yùn)行結(jié)束；

(6)重復(fù)(1)～(5)，直至機(jī)構(gòu)所有軌跡點(diǎn)實(shí)時(shí)預(yù)測(cè)并運(yùn)行結(jié)束。

CTS機(jī)構(gòu)在整個(gè)軌跡預(yù)測(cè)和運(yùn)行過(guò)程中，都以發(fā)揮機(jī)構(gòu)最大能力為原則，但在某些步長(zhǎng)內(nèi)，可能因某些特殊原因不能發(fā)揮機(jī)構(gòu)最大能力時(shí)，則需要進(jìn)行緩沖處理。

2.1 無(wú)緩沖的動(dòng)態(tài)軌跡連續(xù)同步控制

基準(zhǔn)時(shí)間ΔT為一極小值，根據(jù)實(shí)際情況選取。該步長(zhǎng)內(nèi)N段各點(diǎn)速度由控制系統(tǒng)實(shí)時(shí)采集獲得，各自由度速度極限值為Vmax，加速度極限值為amax。

由機(jī)構(gòu)當(dāng)前點(diǎn)實(shí)際位置及速度、目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置、基準(zhǔn)時(shí)間ΔT推導(dǎo)下一點(diǎn)機(jī)構(gòu)位置，則有：

(2)

式(2)中0≤n

圖3 CTS機(jī)構(gòu)同步過(guò)程示意(2軸)

通過(guò)以下幾個(gè)約束條件，可以對(duì)機(jī)構(gòu)位置更新函數(shù)法則Fk,n進(jìn)行求解：(1)機(jī)構(gòu)不斷逼近目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置；(2)發(fā)揮機(jī)構(gòu)最大運(yùn)行能力；(3)基準(zhǔn)時(shí)間ΔT極小，一般以ms為單位；(4)機(jī)構(gòu)各自由度運(yùn)行速度及加速度均不能超出機(jī)構(gòu)運(yùn)行速度Vmax及加速度極限amax。

2.2 有緩沖的動(dòng)態(tài)軌跡連續(xù)同步控制

當(dāng)遇到特殊情況(如某自由度緊急反向)，需要預(yù)先對(duì)機(jī)構(gòu)進(jìn)行緩沖處理，緩沖完成后，機(jī)構(gòu)進(jìn)入目標(biāo)軌跡點(diǎn)預(yù)測(cè)位置的同步運(yùn)行。

CTS機(jī)構(gòu)緩沖同步過(guò)程如圖4所示(圖示僅兩個(gè)自由度)，L(k)至L(k+1)步長(zhǎng)內(nèi)各點(diǎn)位置推導(dǎo)函數(shù)法則仍為：

圖4 CTS機(jī)構(gòu)緩沖同步過(guò)程示意(2軸)

當(dāng)緩沖完成后，機(jī)構(gòu)即進(jìn)入目標(biāo)點(diǎn)預(yù)測(cè)位置的同步運(yùn)行，位置更新及函數(shù)法則Fk,n求解同2.1。

3 天平信號(hào)的實(shí)時(shí)濾波

由式(1)知：CTS機(jī)構(gòu)當(dāng)前位置和六分量天平載荷信號(hào)直接影響目標(biāo)點(diǎn)位置的預(yù)測(cè)，實(shí)際應(yīng)用中，通常采用硬件或軟件濾波器對(duì)天平信號(hào)進(jìn)行濾波處理，以保證機(jī)構(gòu)軌跡預(yù)測(cè)的準(zhǔn)確，圖5為加入天平信號(hào)濾波器的CTS機(jī)構(gòu)控制原理圖。

圖5 加入濾波器的CTS機(jī)構(gòu)控制原理圖

3.1 硬件濾波

為充分濾除高頻干擾對(duì)天平信號(hào)的影響，吹風(fēng)中使用低通濾波器對(duì)天平信號(hào)進(jìn)行濾波處理，其低通截止頻率分別為10Hz、3Hz和1Hz。

從圖6可知，在采用1Hz低通濾波時(shí)，濾波后的信號(hào)相對(duì)于3Hz和10Hz是有滯后的；而硬件10Hz濾波后的信號(hào)擾動(dòng)最大，3Hz擾動(dòng)為其次，1Hz擾動(dòng)最小。

圖6 某型天平硬件10Hz、3Hz、1Hz濾波頻率下的數(shù)據(jù)對(duì)比

Fig.6 The comparison of the test results using hardware filtering on 10Hz, 3Hz and 1Hz of a balance

3.2 卡爾曼濾波

實(shí)際應(yīng)用中，要求濾波滯后應(yīng)盡量小，同時(shí)信號(hào)擾動(dòng)應(yīng)盡量減少。然而濾波滯后和信號(hào)擾動(dòng)是相互矛盾的，為此，本文采取適當(dāng)提高硬件低通濾波截止頻率，通過(guò)后續(xù)軟件濾波的方法來(lái)緩解上述矛盾，將卡爾曼濾波算法引入到CTS機(jī)構(gòu)連續(xù)動(dòng)態(tài)控制中，獲得了很好的效果。

(1)卡爾曼濾波算法簡(jiǎn)介

上世紀(jì)60年代匈牙利數(shù)學(xué)家魯?shù)婪颉·卡爾曼(Rudolf Emil Kalman)在他的博士論文和發(fā)表的論文“A New Approach to Linear Filtering and Prediction Problems”中提出了卡爾曼濾波理論，并導(dǎo)出了一套遞推估計(jì)算法，其基本思想是：采用信號(hào)與噪聲的狀態(tài)空間模型，利用前一時(shí)刻的估計(jì)值和現(xiàn)時(shí)刻的觀測(cè)值來(lái)更新對(duì)狀態(tài)變量的估計(jì)，求出現(xiàn)時(shí)刻的估計(jì)值。

系統(tǒng)離散狀態(tài)方程定義為：

式中：X(k)為系統(tǒng)狀態(tài)向量，k表示時(shí)間點(diǎn)；Φ(k)為狀態(tài)方程系數(shù)矩陣；U(k)為控制項(xiàng)；W(k)均值為0，方差為Q(k)的白噪聲；Z(k)為觀測(cè)向量；H(k)為觀測(cè)方程系數(shù)矩陣；V(k)為均值為0，方差為R(k)的白噪聲。

迭代過(guò)程如下：

系統(tǒng)時(shí)間傳播方程為：

(a)

(b)

在狀態(tài)值預(yù)測(cè)時(shí)，忽略B(k)U(k)項(xiàng)，采用下式作為系統(tǒng)時(shí)間傳播方程：

系統(tǒng)測(cè)量修正方程，即濾波增益為：

(c)

系統(tǒng)狀態(tài)更新方程為：

(d)

(e)

上述“T”表示轉(zhuǎn)置，P(k)為系統(tǒng)狀態(tài)向量協(xié)方差。推導(dǎo)濾波過(guò)程按式(a)～(e)循環(huán)進(jìn)行，直到濾波結(jié)束。

(2)濾波模型建立

在濾波器模型建立中，引入泰勒中值定理：若函數(shù)f(x)在開(kāi)區(qū)間(a,b)上有N+1階導(dǎo)數(shù)，則當(dāng)函數(shù)在(a,b)時(shí)，可展開(kāi)為一個(gè)關(guān)于(x-x0)的多項(xiàng)式和一個(gè)余項(xiàng)的和，其形式為：

式中：x,x0∈(a,b)，fn(x0)是f(x)的n階導(dǎo)數(shù)在x0處的值，Rn稱為拉格朗日余項(xiàng)。

離散系統(tǒng)中，x(k)代表函數(shù)x(t)在t=kT時(shí)刻的采樣值，其中T為x(t)采樣周期。根據(jù)泰勒展開(kāi)式，建立n階系統(tǒng)離散狀態(tài)模型。

系統(tǒng)狀態(tài)向量為：

獲得系統(tǒng)狀態(tài)系數(shù)(T表示采樣周期)矩陣為：

4 試驗(yàn)結(jié)果

圖7給出了風(fēng)洞馬赫數(shù)M=0.5，模型天平某方向載荷在硬件10Hz濾波和2階卡爾曼濾波的數(shù)據(jù)對(duì)比。

圖7 10Hz硬件濾波及2階卡爾曼濾波下的數(shù)據(jù)對(duì)比

Fig.7 The comparison of the test results using hardware filtering on 10Hz and two-rank Kalman filtering

圖8為風(fēng)洞馬赫數(shù)M=0.5，位置控制方式和連續(xù)動(dòng)態(tài)控制方式(加入卡爾曼濾波)下，上述試驗(yàn)?zāi)Ｐ虲TS吹風(fēng)試驗(yàn)的分離軌跡對(duì)比。因各車次試驗(yàn)工況不盡相同，導(dǎo)致了圖8中軌跡分離的偏差，該偏差是在機(jī)構(gòu)運(yùn)行允許范圍內(nèi)的。

(a) X軸分離軌跡對(duì)比

(b) Y軸分離軌跡對(duì)比

Fig.8 The comparison of the test results generated by the move-and-stop control process and the velocity control process

表1給出了風(fēng)洞馬赫數(shù)M=0.5，位置控制方式和連續(xù)動(dòng)態(tài)控制方式(加入卡爾曼濾波)下，上述試驗(yàn)?zāi)Ｐ筒煌嚧蜟TS吹風(fēng)試驗(yàn)的風(fēng)洞運(yùn)行時(shí)間和耗氣量數(shù)據(jù)。

表1 連續(xù)動(dòng)態(tài)控制方式和位置控制方式風(fēng)洞運(yùn)行數(shù)據(jù)

5 結(jié) 論

大量地面和模型吹風(fēng)試驗(yàn)表明：卡爾曼濾波算法在CTS機(jī)構(gòu)中的應(yīng)用有效可行，它改善了硬件(主要是1Hz低通)濾波滯后，同時(shí)滿足系統(tǒng)對(duì)數(shù)據(jù)精度的要求；卡爾曼濾波器的引入對(duì)CTS機(jī)構(gòu)分離軌跡的重復(fù)性并無(wú)不利影響；CTS機(jī)構(gòu)連續(xù)動(dòng)態(tài)控制方式大大提高了風(fēng)洞的吹風(fēng)效率，節(jié)約成本50%以上，其效果已在大量吹風(fēng)試驗(yàn)中得到充分體現(xiàn)。

[1] Hill D W Jr. Development of a velocity control algorithm for controlling a 6-dof trajectory model support[R]. Propulsion Wind Tunnel Facility. Calspan Field Services, Inc./AEDC Division, Arnold Air Force Station, Tennessee.

[2] Roy G. Wind tunnel captive trajectory testing using continuous velocity control[R]. Butler Instrumentation Engineer Propulsion Wind Tunnel Facility, Arnold Air Force Station, Tennessee.

[3] Gargon F, Ph Taravel, Raffin J C. Recent developments in captive trajectory systems of the ONERA modane wind tunnels[C]. 39th Aerospace Sciences Meeting & Exhibit.

[4] Andrew Garrell, Jay Nowakowski, Roman Paryz. Design & development of a new captive trajectory simulation model support system[C]. 35th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, 1999.

[5] Zhu Rusong, Wang Fei, Li Ping, et al. The continuous dynamic simulation test of the captive trajectory system used in 2.4m transonic wind tunnel[C]. CCDC, 2011

[6] Zhou D H, Wang Q L. Strong tracking filtering of nonlinear systems with colored noise[J]. J of Beijing Institute of Technology, 1997, 17(3): 321-326.

[7] Paulo S R Diniz. Adaptive filtering: algorithms and practical implementation[M]. 劉郁林, 景曉軍, 等譯. 北京：電子工業(yè)出版社.

[8] 黃敘輝, 羅新福, 于志松. FL-24風(fēng)洞新型捕獲軌跡系統(tǒng)設(shè)計(jì)與發(fā)展[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2008, 26(2): 145-149.

Huang X H, Luo X F, Yu Z S. Design and development of a new captive trajectory simulation system in FL-24 wind tunnel[J]. Acta Aerodynamica Sinica, 2008, 26(2):145-149.

[9] 黃敘輝, 龐旭東, 宋斌. 1.2m跨超聲速風(fēng)洞新型捕獲軌跡系統(tǒng)研制[J]. 實(shí)驗(yàn)流體力學(xué), 2008, 22(2): 95-98.

Huang X H, Pang X D, Song B. Development of a new captive trajectory simulation system in the 1.2m transonic wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2008, 22(2): 95-98.

[10] 高勤, 李志強(qiáng), 都學(xué)新. 一種新型自適應(yīng)卡爾曼濾波算法[J]. 現(xiàn)代雷達(dá), 2001, 12(6): 29-34.

Gao Q, Li Z Q, Du X X. A new adaptive Kalman filtering algoritnm[J]. Modern Radar, 2001, 12(6): 29-34.

[11] 范文兵, 劉春風(fēng), 張淑貞. 一種強(qiáng)跟蹤擴(kuò)展卡爾曼濾波器的改進(jìn)算法[J]. 控制與決策, 2006, 36(1): 52-56.

Fan W B, Lin C F, Zhang S Z. Improved method of strong tracking extended Kalman filter[J]. Control and Decision, 2006, 36(1): 52-56.

[12] 虞旦, 韋巍, 張遠(yuǎn)輝. 一種基于卡爾曼預(yù)測(cè)的動(dòng)態(tài)目標(biāo)跟蹤算法研究[J]. 光電工程, 2009, 21(1): 73-76.

Yu D, Wei W, Zhang Y H. Dynamic target tracking with Kalman filter as predictor[J]. Opto-Electronic Engineering, 2009, 21(1): 73-76.

(編輯：張巧蕓)

更 正 啟 事

本刊2014年第5期第112頁(yè)表2“優(yōu)秀會(huì)議論文簡(jiǎn)介”中序號(hào)1的論文信息有誤。原論文題名和作者“A cylindrical converging shock tube for Richtmyer-Meshkov instability study”和“Juchun Ding, Minghu Wang, Zhigang Zhai, Ting Si, Xisheng Luo”，應(yīng)更正為“Experimental study on the interaction of a planar shock wave with a three-dimensional gaseous cylinder”和“Mojun Chen, Zhigang Zhai, Ting Si, Xisheng Luo”。由于本刊工作失誤對(duì)相關(guān)作者和讀者造成了困擾，特此致歉，同時(shí)向發(fā)現(xiàn)該錯(cuò)誤的中國(guó)科學(xué)技術(shù)大學(xué)陳模軍表示由衷的感謝。

Application of Kalman filtering in continuous dynamic trajectory simulation technique of the captive trajectory system used in a transonic wind tunnel

Wang Fei1,2, Zhu Rusong1,2, Zhang Junsheng1,2

(1. Facility Design & Instrumentation Institute, China Aerodynamics Research and Development Center, Mianyang Sichuan 621000, China; 2. The State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang Sichuan 621000, China)

The flying trajectory simulation of the store after being separated from the mother plane is an advanced and difficult technique in wind tunnel test. In this paper, the continuous dynamic trajectory simulation technique is studied in the Captive Trajectory System (CTS) used in a transonic wind tunnel. This paper initially introduces the dynamic trajectory predicting principle. Then, the specific steps of the coordinated control method of continuous dynamic trajectory simulation are amply described. The algorithm and constraint condition of the coordinated control and the buffering coordinated control are also discussed. Finally, the Taylor mean value theorem is applied to build a Kalman filtering model to resolve the conflict between filtering delay and accuracy requirements. After choosing appropriate parameter, the Kalman filtering is used to process the balance signal dynamically and in real time, and the comparison between the continuous dynamic control process and the move-and-stop control process is carried out. The wind tunnel CTS test result shows that the Kalman filtering can reduce the signal fluctuation greatly, meet the real-time requirements of CTS, and dynamic separation trajectory of the continuous dynamic control process metches well with the one of the move-and-stop cortrol, while the wind tunnel CTS test efficiency has also been considerably improved.

transonic wind tunnel;captive trajectory system;continuous;dynamic;Kalman filtering

1672-9897(2015)01-0103-06

10.11729/syltlx20120029

2013-03-03;

2014-04-10

WangF,ZhuRS,ZhangJS.ApplicationofKalmanfilteringincontinuousdynamictrajectorysimulationtechniqueofthecaptivetrajectorysystemusedinatransonicwindtunnel.JournalofExperimentsinFluidMechanics, 2015, 29(1): 103-108. 王 飛, 祝汝松, 張俊生. 卡爾曼濾波在某跨聲速風(fēng)洞CTS機(jī)構(gòu)連續(xù)動(dòng)態(tài)軌跡捕獲試驗(yàn)技術(shù)中的應(yīng)用. 實(shí)驗(yàn)流體力學(xué), 2015, 29(1): 103-108.

V211.73

A

王 飛(1982-)，男，碩士，江蘇阜寧人，主要從事風(fēng)洞測(cè)控系統(tǒng)設(shè)計(jì)、開(kāi)發(fā)與研究。通信地址：四川省綿陽(yáng)市二環(huán)路南段6號(hào)設(shè)備設(shè)計(jì)及測(cè)試技術(shù)研究所403室(621000)。E-mail：wfei_fancy@163.com