一種寬帶復(fù)相關(guān)流速測(cè)量方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種流速測(cè)量方法，尤其是涉及一種寬帶復(fù)相關(guān)流速測(cè)量方法。
【背景技術(shù)】
[0002] 隨著資源日益匱乏，開發(fā)海洋水域資源成為經(jīng)濟(jì)發(fā)展的急迫需求。海洋水域的流 速為海洋水域環(huán)境的重要參數(shù)，也是海洋水域活動(dòng)的重要依據(jù)?，F(xiàn)有的流速測(cè)量方法主要 是基于寬帶復(fù)相關(guān)技術(shù)，文獻(xiàn)《寬帶多普勒海流計(jì)測(cè)速方法研究》提出基于寬帶復(fù)相關(guān)算法 的流速研究。該文獻(xiàn)中，流速測(cè)量的具體過程為:首先利用脈沖信號(hào)發(fā)生器產(chǎn)生單頻脈沖信 號(hào)，采用Μ序列碼對(duì)該單頻脈沖信號(hào)進(jìn)行二進(jìn)制編碼得到子脈沖信號(hào)，將多個(gè)子脈沖信號(hào)依 次進(jìn)行拼接得到調(diào)制脈沖信號(hào);接收該調(diào)制脈沖信號(hào)經(jīng)海洋水域反射回的信號(hào)，并采集得 到該信號(hào)的回波信號(hào);最后采用寬帶復(fù)相關(guān)算法對(duì)該回波信號(hào)進(jìn)行處理，得到流速。
[0003] 該流速測(cè)量方法中，回波信號(hào)是一個(gè)含噪聲的隨機(jī)信號(hào)，回波信號(hào)的強(qiáng)度、噪聲干 擾和回波信號(hào)的傳播損耗等都影響流速測(cè)量的精度，復(fù)相關(guān)算法能夠?yàn)V除一定的白噪聲， 具有較好的抗干擾能力，由此該流速測(cè)量方法精度較高，但是復(fù)相關(guān)算法在對(duì)回波信號(hào)進(jìn) 行處理的過程中，得到的回波信號(hào)的復(fù)相關(guān)函數(shù)的相位信息具有周期性，流速的測(cè)量范圍 受到了限制，由此該方法量程較窄。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明所要解決的技術(shù)問題是提供一種測(cè)量精度高，且測(cè)量量程寬的寬帶復(fù)相關(guān) 流速測(cè)量方法。
[0005] 本發(fā)明解決上述技術(shù)問題所采用的技術(shù)方案為:一種寬帶復(fù)相關(guān)流速測(cè)量方法， 包括以下步驟：
[0006] ①獲取回波信號(hào)
[0007] ①-1利用脈沖信號(hào)發(fā)生器生成頻率為fQ的單頻脈沖信號(hào)，fQ= 1.5MHz;
[0008] ①-2采用碼元數(shù)為L(zhǎng)、碼元寬度為△ τ的Μ序列碼對(duì)單頻脈沖信號(hào)進(jìn)行二進(jìn)制編碼 得到子脈沖信號(hào)，碼元數(shù)L=15,碼元寬度Δτ為填充10個(gè)單頻脈沖信號(hào)周期，
將子脈沖信號(hào)的長(zhǎng)度記為T，T = LX Δ τ，符號(hào)X為乘運(yùn)算符號(hào)；
[0009] ①-3將多個(gè)子脈沖信號(hào)在待測(cè)流速的海洋水域中依次連續(xù)發(fā)射，依次連續(xù)發(fā)射的 多個(gè)子脈沖信號(hào)組成一個(gè)調(diào)制脈沖信號(hào)，將該調(diào)制脈沖信號(hào)記為s(t)，t為時(shí)域時(shí)間；
[0010] ①-4獲取該調(diào)制脈沖信號(hào)s(t)經(jīng)海洋水域反射回的信號(hào)，將該反射回的信號(hào)記為 8'（〇，根據(jù)8'（〇得到待處理的回波信號(hào)^(〇，^(〇 = 8'（〇+5(〇;其中5(〇為白噪聲信 號(hào)；
[0011] ②對(duì)回波信號(hào)rr(t)依次進(jìn)行前置放大處理、帶通濾波處理和低通濾波處理，得到 預(yù)采樣信號(hào)，將預(yù)采樣信號(hào)記為r (t);
[0012] ③采用交替采樣模塊對(duì)預(yù)采樣信號(hào)r(t)進(jìn)行四路時(shí)間交替采樣得到離散采樣序 列r(n)，n為離散序列號(hào)，1!=1，2,3，一；交替采樣模塊包括第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換 器、第三模數(shù)轉(zhuǎn)換器、第四模數(shù)轉(zhuǎn)換器和四選一選擇器，第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、 第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器和四選一選擇器的四個(gè)選擇端一一對(duì)應(yīng)連接，第一模數(shù) 轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器對(duì)預(yù)采樣信號(hào)r(t)進(jìn)行采樣， 四選一選擇器依次選擇第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn) 換器，第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器的采樣信號(hào)經(jīng) 四選一選擇器依次輸出按序拼接為離散采樣序列r(n);
[0013]④對(duì)離散采樣序列r(n)分別作復(fù)相關(guān)算法處理和互譜自相關(guān)算法處理；
[0014]復(fù)相關(guān)算法處理的具體過程為：
[0015] a.采用DDS模塊產(chǎn)生正交調(diào)制信號(hào)，采用該正交調(diào)制信號(hào)對(duì)離散采樣序列r(n)作 數(shù)字調(diào)制，得到兩路正交信號(hào)序列，分別記為I(n)和Q(n);采用I(n)和Q(n)構(gòu)造第一復(fù)數(shù)序 列X(n)，將I(n)作為第一復(fù)數(shù)序列X(n)的實(shí)部信號(hào)，Q(n)作為第一復(fù)數(shù)序列X(n)的虛部信 號(hào)，則第一復(fù)數(shù)序列表示X(n) = I(n)+jQ(n)，其中j表示復(fù)數(shù)的虛部單位；
[0016] b.將第一復(fù)數(shù)序列X(n)按照奇數(shù)和偶數(shù)離散序列號(hào)分為兩個(gè)復(fù)數(shù)序列，將離散序 列號(hào)為奇數(shù)的復(fù)數(shù)序列稱為第二復(fù)數(shù)序列，記為X(2m+1)，將離散序列號(hào)為偶數(shù)的復(fù)數(shù)序列 稱為第三復(fù)數(shù)序列，記為X(2m+2)，m = 0，1，2，3，…；
[0017] c.將第一復(fù)數(shù)序列X(n)的復(fù)相關(guān)函數(shù)記為RKn)，將第二復(fù)數(shù)序列X(2m+1)的復(fù)相 關(guān)函數(shù)記為R2 ( η )，將第三復(fù)數(shù)序列X ( 2m + 2 )的復(fù)相關(guān)函數(shù)記為R3 ( η )，
?其中Ν'表示第二復(fù) 數(shù)序列X(2m+1)的復(fù)相關(guān)長(zhǎng)度，Ν"表示第三復(fù)數(shù)序列X(2m + 2)的復(fù)相關(guān)長(zhǎng)度，且
,1為采樣時(shí)延值且其取值為1個(gè)Μ序列碼的編碼長(zhǎng)度，N = fs X (ti-to)，匕為 交替采樣模塊的采樣頻率，to為交替采樣模塊的采樣起始時(shí)刻，交替采樣模塊的采樣結(jié)束 時(shí)刻，X(2m+l+l) = I(2m+l+l)+jQ(2m+l+l)，X(2m+2+l) = I(2m+2+l)+jQ(2m+2+l)，I(2m+l+ 1)為I(2m+1)延時(shí)1的表達(dá)式，I(2m+2+l)為I(2m+2)延時(shí)1的表達(dá)式，Q(2m+1+1)為Q(2m+1)時(shí) 延1的表達(dá)式，Q(2m+2+l)為Q(2m+2)時(shí)延1的表達(dá)式;符號(hào)*為共輒運(yùn)算符，X\2m+1)為X(2m+ 1)的共輒復(fù)數(shù)序列，X$( 2m+2)為X (2m+2)的共輒復(fù)數(shù)序列；
[0018] d.將 R2(n)和 R3(n)代入公式1^(11)=1?2(11)+1?3(11)中得到 RKnM^RKn)的實(shí)部記為 ReUKn))，將將RKn)的虛部記為ImUKn));
[0019] e.將第一復(fù)數(shù)序列X(n)的復(fù)相關(guān)函數(shù)RKn)的相位角記為Θ，采用反正切函數(shù)求解

[0020] f.將多普勒頻偏記為fd，采用公式
其中q為相干重復(fù)周期，即q = fsX T；
[0021] g.將復(fù)相關(guān)測(cè)量的流速記為V。 計(jì)算得到Vc=，c為水中的聲 速，c = 1500m/s;
[0022] 互譜自相關(guān)算法處理的具體過程為：
[0023] a.將離散采樣序列r(n)的頻譜記為r'（k)，離散采樣序列r(n)的自相關(guān)函數(shù)頻譜 記為R(k)，k表示頻域的序列號(hào)，k = 0，l，2,3,…；
[0024] b.采用 FFT 運(yùn)算計(jì)算得到頻譜 r'（k)，令 R(k)=r'（k)Xr'（k);
[0025] c.對(duì)自相關(guān)函數(shù)頻譜R(k)作插值處理后再進(jìn)行IFFT運(yùn)算，得到離散采樣序列r(n) 的自相關(guān)函數(shù)，將離散采樣序列r(n)的自相關(guān)函數(shù)記為R'（t);
[0026] d.通過離散采樣序列r(n)的自相關(guān)函數(shù)R'（t)確定第一旁瓣出現(xiàn)的時(shí)刻，將第一 旁瓣出現(xiàn)的時(shí)刻記為T1;將不發(fā)生頻偏的理想狀態(tài)下第一旁瓣出現(xiàn)的時(shí)刻記為T2,其中T 2 = Τ
[0027] e.將發(fā)生多普勒時(shí)延的差值記為ΔΤ，Δ 1 = ;
[0028] f.將自相關(guān)測(cè)量的流速記為vz，
計(jì)算得到自相關(guān)測(cè)量的流 速vz;
[0029] ⑤將流速修正參數(shù)記為α，a = [ ( Vz-Vc^/^Vmax]，Vmax為模糊速度，Vmax = c/4f〇T，符號(hào) "[]"為取整運(yùn)算符；
[0030] ⑥采用公式V = Vc+a2Vmax計(jì)算得到待測(cè)流速V。
[0031] 與現(xiàn)有技術(shù)相比，本發(fā)明的優(yōu)點(diǎn)在于采用時(shí)間交替采樣技術(shù)對(duì)回波信號(hào)進(jìn)行采樣 得到離散采樣序列，另外分別對(duì)離散采樣序列作復(fù)相關(guān)算法處理和互譜自相關(guān)算法處理， 時(shí)間交替采樣在保證采樣精度的前提下提高了采樣率，復(fù)相關(guān)算法處理過程中相應(yīng)增加了 復(fù)相關(guān)點(diǎn)數(shù)，并米用奇偶兩路并彳丁復(fù)相關(guān)運(yùn)算處理，在提尚復(fù)相關(guān)算法的測(cè)流精度的基礎(chǔ) 上，提高計(jì)算速度，保證流速測(cè)量的實(shí)時(shí)性，互譜自相關(guān)算法處理通過對(duì)時(shí)間交替采樣數(shù)據(jù) 作互譜插值延時(shí)估計(jì)，提高時(shí)延值的估計(jì)精度，減小流速大小區(qū)間的判斷誤差，實(shí)現(xiàn)寬量程 流速測(cè)量，采用該延時(shí)值粗測(cè)流速值，該粗測(cè)流速值不存在測(cè)量模糊問題，可作為流速值大 小區(qū)間的判斷標(biāo)準(zhǔn)，輔助復(fù)相關(guān)算法流速測(cè)量，克服復(fù)相關(guān)算法流速測(cè)量可能存在模糊值 和流速測(cè)量量程的局限性問題，修正流速測(cè)量量程，由此復(fù)相關(guān)算法處理和互譜自相關(guān)算 法處理相結(jié)合，通過采用時(shí)間交替采樣技術(shù)提高信號(hào)的采樣率，結(jié)合互譜自相關(guān)算法粗測(cè) 的流速值輔助復(fù)相關(guān)算法進(jìn)行流速測(cè)量，測(cè)量精度高，且測(cè)量量程寬。
【具體實(shí)施方式】
[0032] 以下結(jié)合實(shí)施例對(duì)本發(fā)明作進(jìn)一步詳細(xì)描述。
[0033] 實(shí)施例:一種寬帶復(fù)相關(guān)流速測(cè)量方法，包括以下步驟：
[0034]①獲取回波信號(hào)
[0035]①-1利用脈沖信號(hào)發(fā)生器生成頻率為fo的單頻脈沖信號(hào)，fo= 1.5MHz;
[0036] ①-2采用碼元數(shù)為L(zhǎng)、碼元寬度為△ τ的Μ序列碼對(duì)單頻脈沖信號(hào)進(jìn)行二進(jìn)制編碼
得到子脈沖信號(hào)，碼元數(shù)L=15,碼元寬度Δτ為填充10個(gè)單頻脈沖信號(hào)周期 將子脈沖信號(hào)的長(zhǎng)度記為T，T = LX Δ τ，符號(hào)X為乘運(yùn)算符號(hào)；
[0037] ①-3將多個(gè)(即至少兩個(gè))子脈沖信號(hào)在待測(cè)流速的海洋水域中依次連續(xù)發(fā)射，依 次連續(xù)發(fā)射的多個(gè)子脈沖信號(hào)組成一個(gè)調(diào)制脈沖信號(hào)，將該調(diào)制脈沖信號(hào)記為S(t)，t為時(shí) 域時(shí)間；
[0038]①-4獲取該調(diào)制脈沖信號(hào)s(t)經(jīng)海洋水域反射回的信號(hào)，將該反射回的信號(hào)記為 8'（〇，根據(jù)8'（〇得到待處理的回波信號(hào)^(〇，^(〇 = 8'（〇+5(〇;其中5(〇為白噪聲信 號(hào)；
[0039]②對(duì)回波信號(hào)rr(t)依次進(jìn)行前置放大處理、帶通濾波處理和低通濾波處理，得到 預(yù)采樣信號(hào)，將預(yù)采樣信號(hào)記為r (t);
[0040] ③采用交替采樣模塊對(duì)預(yù)采樣信號(hào)r(t)進(jìn)行四路時(shí)間交替采樣得到離散采樣序 列r(n)，n為離散序列號(hào)，1!=1，2,3，一；交替采樣模塊包括第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換 器、第三模數(shù)轉(zhuǎn)換器、第四模數(shù)轉(zhuǎn)換器和四選一選擇器，第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、 第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器和四選一選擇器的四個(gè)選擇端一一對(duì)應(yīng)連接，第一模數(shù) 轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器對(duì)預(yù)采樣信號(hào)r(t)進(jìn)行采樣， 四選一選擇器依次選擇第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn) 換器，第一模數(shù)轉(zhuǎn)換器、第二模數(shù)轉(zhuǎn)換器、第三模數(shù)轉(zhuǎn)換器和第四模數(shù)轉(zhuǎn)換器的采樣信號(hào)經(jīng) 四選一選擇器依次輸出按序拼接為離散采樣序列r(n);
[0041] ④對(duì)離散采樣序列r(n)分別作復(fù)相關(guān)算法處理和互譜自相關(guān)算法處理；
[0