韓亞秋， 賈云飛

(南京理工大學(xué)， 機(jī)械工程學(xué)院， 江蘇 南京 210094)

結(jié)合CPLD測時的改進(jìn)互相關(guān)算法在紅外測速系統(tǒng)中的應(yīng)用

韓亞秋， 賈云飛

(南京理工大學(xué)， 機(jī)械工程學(xué)院， 江蘇 南京 210094)

為了解決直接采用傳統(tǒng)互相關(guān)算法計算炮彈過靶時間差計算量大而不利于實時測量的問題， 提出了一種結(jié)合CPLD的改進(jìn)的互相關(guān)算法計算炮彈穿過兩靶的時間差的方法. 在分析炮彈穿過光幕時信號變化過程的基礎(chǔ)上， 設(shè)計CPLD測時模塊測出彈丸過靶的時間差的估計值， 根據(jù)CPLD測時結(jié)果篩選出參與互相關(guān)測時的相應(yīng)時間段的信號， 再進(jìn)行互相關(guān)分析從而得出炮彈過靶時間差， 計算出速度. 最后通過MATLAB對仿真信號進(jìn)行分析， 表明采用改進(jìn)的互相關(guān)算法不僅具有較少的計算量且有較高的精度.

CPLD測時； 互相關(guān)算法； MATLAB仿真； 紅外測速系統(tǒng)； 實時測量

0 引 言

傳統(tǒng)的測試彈丸飛行速度的方法是將區(qū)截裝置和測時儀結(jié)合起來， 其原理為測得兩個區(qū)截裝置的間距， 利用測時儀得到彈丸穿過兩個區(qū)截裝置的時間差， 從而計算出彈丸穿過兩區(qū)截裝置中間位置的平均速度[1].

對于一般輕武器的子彈速度測量， 區(qū)截裝置能很好地捕捉彈丸過靶信號并觸發(fā)測試儀測時， 目前已經(jīng)達(dá)到0.1%的精度， 但是對于炮彈大直徑的彈丸測量， 由于發(fā)射過程中有彈頭波的影響產(chǎn)生測時儀的誤觸發(fā)， 從而得不到彈丸的實際值， 因此提出了使用互相關(guān)算法進(jìn)行彈丸穿過兩個區(qū)截裝置時間差的計算. 目前而言該方法多采用采集儀采集數(shù)據(jù)， 然后將數(shù)據(jù)在MATLAB中進(jìn)行互相關(guān)分析， 但該方法直接將采集到的數(shù)據(jù)進(jìn)行互相關(guān)算法計算時間差， 需要參與運算的數(shù)據(jù)量較大， 無法實現(xiàn)實時測量. 因此， 提出了結(jié)合CPLD實現(xiàn)一種改進(jìn)的互相關(guān)算法計算， 先利用硬件編程語言設(shè)計CPLD測時功能， 由CPLD測時先估計彈丸過靶信號的時間差， 再由CPLD測時結(jié)果截取參與互相關(guān)運算的數(shù)據(jù)時間段， 該方法使參與互相關(guān)運算的數(shù)據(jù)量減少從而使計算量大大減少[2].

1 結(jié)合CPLD的改進(jìn)的互相關(guān)測時原理

圖 1 基于改進(jìn)的互相關(guān)算法紅外測速系統(tǒng)原理示意圖Fig.1 Schematic diagram of infrared velocity measurement system based on improved cross-correlation algorithm

基于改進(jìn)的互相關(guān)算法的紅外測速系統(tǒng)是雙區(qū)截裝置測速系統(tǒng)的一種實現(xiàn)方式， 我們也常將它稱為紅外測速靶， 測試原理如圖 1 所示. 光幕是由多個紅外二極管排成隊列形成的， 該光幕具有一定的寬度. 測出其兩靶之間的長度L和兩個觸發(fā)信號的時間差Δt， 利用式(1)確定彈丸在BC距離上彈丸的平均速度， 也可以認(rèn)為是兩光幕中間測量點D的瞬時飛行速度[3].

當(dāng)子彈穿過兩靶的紅外光幕時， 引起光通量的變化， 光電接收器件感受光通量變化并將其轉(zhuǎn)換為電信號， 分別產(chǎn)生兩個觸發(fā)信號. 信號調(diào)理電路實現(xiàn)濾波和放大功能， CPLD測出其彈丸穿過兩靶大致時間， 數(shù)據(jù)采集系統(tǒng)采集彈丸過靶的模擬信號并按照略大于CPLD的測時結(jié)果保存采集數(shù)據(jù)， 最后由主控模塊完成采集數(shù)據(jù)的互相關(guān)分析得到最終精確的測時結(jié)果.

為了減少互相關(guān)算法的計算量， 設(shè)計了結(jié)合CPLD測時和互相關(guān)算法測時相結(jié)合的方法. CPLD測時是基于硬件編程語言實現(xiàn)的， 當(dāng)?shù)谝粋€靶的模擬信號電壓達(dá)到一定閾值(固定閾值)即觸發(fā)CPLD計時器開始計時， 第二個靶的信號電壓達(dá)到固定閾值即停止計時， 然后將測時結(jié)果發(fā)送給數(shù)據(jù)分析模塊.

理想情況下,設(shè)紅外測速靶的原理參數(shù)基本一致,則當(dāng)子彈飛過兩光幕時引起的波形完全相似， 因此可以通過采用互相關(guān)算法得到兩個波形之間的關(guān)系， 從而得到彈丸過靶時的時間差. 設(shè)彈丸穿過Ⅰ靶和Ⅱ靶的信號波形分別為x1(t)，x2(t)， 則它們的互相關(guān)函數(shù)為

式中：τ代表相關(guān)函數(shù)自變量；T為兩個信號的脈寬. 通常模擬信號采集下來就成了數(shù)字離散信號[4-5]. 則式(2)的離散形式為

(3)

根據(jù)相關(guān)函數(shù)的定義可知， 當(dāng)τ取不同的值時，Rx1x2(τ)也取不同的值與之對應(yīng)[6]， 而Rx1x2(τ)的最大值時的τ值便是彈丸穿過兩靶的時間間隔.

圖 2 彈丸穿過兩光幕靶的信號示意圖Fig.2 Schematic diagram of projectile signal through the two light screens

如圖 2 所示為彈丸穿過兩靶的信號波形圖， 圖中第一條為彈丸穿過光幕靶的信號波形， 第二條為穿過光幕靶的信號波形， 假設(shè)采集的信號時間長度為0～T0， CPLD測得時間為T1， 那么在CPLD測得時間結(jié)果基礎(chǔ)上對信號長度進(jìn)行截取， 保留第一個靶信號對應(yīng)時間為0～(T0-T1)的信號長度， 第二個靶信號對應(yīng)T1～T0的信號長度， 將保留相應(yīng)時間長度的信號作為新的時間序列的信號數(shù)據(jù).

2 基于MATLAB的改進(jìn)互相關(guān)測時算法仿真

由改進(jìn)的互相關(guān)算法可知， 參與運算的數(shù)據(jù)點大大減少， 使互相關(guān)算法計算量減少， 但是它的可靠性， 精度和穩(wěn)定性需要進(jìn)一步分析. 采用MATLAB仿真改進(jìn)的互相關(guān)算法的計算結(jié)果可以反映其精度.

圖 3 較大直徑的彈丸過靶信號圖Fig.3 The signal of thelarger diameter projectile over the target

圖 3 為示波器采集的直徑相對較大的彈丸的過靶信號， 從圖3中可以看出彈丸過靶時前端有個明顯的較小幅值的尖峰， 尖峰常常是彈丸飛行時的彈頭波產(chǎn)生的， 由于彈頭波會使傳統(tǒng)的觸發(fā)測時儀的光幕靶所測的數(shù)據(jù)與彈丸的速度有較大差距. 下面通過MATLAB仿真彈丸過靶的信號并運用互相關(guān)算法進(jìn)行計算.

CPLD測時原理為： 當(dāng)?shù)谝煌ǖ佬盘柍^一定的幅值就觸發(fā)計時模塊， 第二通道信號超過幅值時停止計時. 在MATLAB仿真時采用固定閾值方法模擬CPLD測時， 取一定的時間差T1， 即當(dāng)?shù)谝煌ǖ赖姆抡娌ㄐ纬^固定閾值(如2 V)后記錄該點的時間點； 第二通道的信號也按此方法記錄時間點， 二者之差為第一次時間測量的T1； 由于實際兩個通道的信號幅值范圍不同， 為方便數(shù)據(jù)處理， 對信號幅值做進(jìn)一步處理， 設(shè)兩個通道的最大幅值分別為xmax，ymax， 則每個通道信號隨時間變化的幅值x(n),y(n)分別為

圖 4 彈丸過靶信號仿真波形及互相關(guān)函數(shù)結(jié)果Fig.4 Simulation waveforms and cross correlation function results of projectile.

從圖 4 互相關(guān)函數(shù)離散圖大致看出為300左右的時間單位， 由MARLAB峰值檢測得出互相關(guān)函數(shù)最大值為300. 對比采用改進(jìn)的互相關(guān)算法的仿真結(jié)果， 由固定閾值方法得到T1=299， 圖 5 所示為采用固定閾值讀出時間后截取有效信號的波形圖及互相關(guān)算法計算結(jié)果.

圖 5 中截取后的信號波形的數(shù)據(jù)點為100， 而原始的信號點為1 000， 將參與互相關(guān)計算的數(shù)據(jù)降低了1/10， 其互相關(guān)計算的結(jié)果在0值附近； 將改進(jìn)后互相關(guān)計算的結(jié)果采用二次峰值檢索方式檢索其峰值點對應(yīng)的時間點， 第一次檢索采用隔點檢索， 即檢索0,2,4……等偶數(shù)序列的數(shù)據(jù)值， 檢索到最大值后前后各取互相關(guān)函數(shù)離散序列10個數(shù)， 進(jìn)行二次檢索， 二次檢索逐個檢索， 最后檢索到最大值， 其檢索的結(jié)果如圖 6 所示， 圖6顯示二次檢索的結(jié)果時間差為1， 則總時間差為300， 與理論值相同.

圖 5 截取后的波形及互相關(guān)函數(shù)離散圖Fig.5 Waveform after interception and cross correlation function

圖 6 互相關(guān)二次峰值檢索的結(jié)果Fig.6 The results of the second peak search

基于MATLAB的仿真結(jié)果與理論值如表 1、 表 2 所示， 時間長度以單位時間個數(shù)來算， 從仿真結(jié)果將二次計算的測時結(jié)果與理論值對比發(fā)現(xiàn)， 二次計算的結(jié)果與理論值結(jié)果吻合， 且誤差在正負(fù)1個時間單位.

表 1 1 000個時間單位的仿真結(jié)果對比

表 2 2 000個時間單位的仿真結(jié)果對比

3 結(jié) 論

結(jié)合CPLD測時與互相關(guān)算法計算時差的方法經(jīng)過仿真計算可以看出： 該方法從減少計算點數(shù)出發(fā)， 能大幅度減少參與互相關(guān)運算的數(shù)據(jù)量， 且仿真結(jié)果與理論值吻合， 具有較高的精度； 從實現(xiàn)該方法的角度看， CPLD測時模塊實現(xiàn)起來簡單， 互相關(guān)算法計算量的減少會更利于實現(xiàn)實時測量， 因此該方法對于紅外測速系統(tǒng)實現(xiàn)互相關(guān)實時測量有重要意義.

[1] 智小軍， 倪晉平， 田會， 等. 相關(guān)分析技術(shù)在光幕靶速度測量中的應(yīng)用[J]. 探測與控制學(xué)報， 2005， 27(1)： 44-46. Zhi Xiaojun， Ni Jinping， Tian Hui， et al. Application of cross-correlation analysis techniques in velocity measurement of light screen[J]. Journal of Detection & Control， 2005， 27(1)： 44-46. (in Chinese)

[2] 范磊. 基于S3C6410的粒子圖像測速技術(shù)研究[D]. 南京： 南京理工大學(xué)， 2014.

[3] 劉繼平. 無線式高精度紅外測速系統(tǒng)研究[D]. 南京： 南京理工大學(xué)， 2015.

[4] 胡廣書. 數(shù)字信號處理理論、 算法與實現(xiàn)[M] . 北京： 清華大學(xué)出版社， 1997.

[5] 董濤， 倪晉平， 宋玉貴， 等. 高射頻武器彈丸連發(fā)速度測量系統(tǒng)[J]. 彈道學(xué)報， 2010， 22(1)： 33-36. Dongtao， Ni JinPing， Song Yugui， et al. Velocity measurement system of high-frequ-ency string-firing bullet[J]. 2010， 22(1)： 33-36. (in Chinese)

[6] 陳興旺. 大口徑彈丸光幕靶測速技術(shù)[J]. 四川兵工學(xué)報， 2012， 33(12)： 7-10. Chen Xingwang. Light screen of velocity measurement technology of large caliber projectile[J]. Journal of Sichuan Ordnance， 2012， 33(12)： 7-10. (in Chinese)

[7] 趙萍， 劉洪利. 一種基于速度檢測的快速互相關(guān)算法[J]. 信息技術(shù)， 2010， 10(3)： 5-7-10. Zhao Ping， Liu Hongli. A fast algorithm for cross-correlation based on velocity measureing[J]. Information Technology， 2010， 10(3)： 5-7-10. (in Chinese)

Research on the Application of Improved Cross Correlation Algorithm in Combined with CPLD Time Measurement in Infrared Velocity Measurement System

HAN Yaqiu, JIA Yunfei

(College of Mechanical Engineering, Nanjing University of Science and Technology, Jiangsu 210094, China)

Thetraditional cross-correlation algorithm has a large amount of computation in calculating time difference of projectile passing through two targets,it is not conducive to real-time measurement.In order to solve the problem， an improved cross correlation algorithm combinedwith CPLD time measurement was proposed to calculate the time difference.Based on the analysis of thesignal changes of the projectile passing through the targets,the estimated value of the time difference was measured by the CPLD test module.It was proposed that the cross correlation analysis the effective signals after screening based on the the estimated valueto obtain the time difference of the projectile, and then calculating speed.Finally, through the analysis of the MATLAB simulation signal, it shows that the improved algorithm has less computation and higher precision.

CPLD measurement; cross correlation algorithm; MATLAB simulation; infrared velocity measurement system; real-time measurement

1671-7449(2017)01-0029-05

2016-10-15

韓亞秋(1990-)， 女， 碩士生， 主要從事嵌入式智能儀器設(shè)備開發(fā)， 武器外彈道測試技術(shù)等研究.

TH89

A

10.3969/j.issn.1671-7449.2017.01.005