本發(fā)明涉及雷達探測技術(shù)領域，尤其涉及一種獲取實時脈沖多普勒譜的方法及裝置。

背景技術(shù)：

在雷達探測過程中，獲取到回波脈沖數(shù)據(jù)后往往還需要獲取回波脈沖數(shù)據(jù)的脈沖多普勒譜，如為提取目標的運動特征，則需要計算不斷積累的目標回波脈沖數(shù)據(jù)的脈沖多普勒譜。針對脈沖多普勒譜計算，目前通常是直接基于離散傅里葉變換(discretefouriertransform，dft)方法進行計算，即每次接收到回波脈沖數(shù)據(jù)時，對所有的回波脈沖數(shù)據(jù)進行dft變換以計算得到脈沖多普勒譜，但上述直接基于dft計算脈沖多普勒譜的方式存在以下問題：

(1)dft變換時要求計算點數(shù)需為2的冪次方，因而直接基于dft計算脈沖多普勒譜時，脈沖多普勒的計算點數(shù)必須為2的冪次方，無法實現(xiàn)任意點數(shù)的脈沖多普勒計算，使用不靈活；

(2)對于每次得到的回波脈沖數(shù)據(jù)均是基于全部脈沖回波數(shù)據(jù)進行計算，需要執(zhí)行大量的計算，計算復雜度高，且計算過程中存在大量的脈沖回波數(shù)據(jù)重復計算。

技術(shù)實現(xiàn)要素：

本發(fā)明要解決的技術(shù)問題就在于：針對現(xiàn)有技術(shù)存在的技術(shù)問題，本發(fā)明提供一種實現(xiàn)方法簡單、計算量小、復雜度低，且使用不受限制的獲取實時脈沖多普勒譜的方法及裝置。

為解決上述技術(shù)問題，本發(fā)明提出的技術(shù)方案為：

一種獲取實時脈沖多普勒譜的方法，步驟包括：

s01.系數(shù)矩陣確定：根據(jù)所需計算的點數(shù)確定用于計算脈沖多普勒譜的系數(shù)矩陣；

s02.多普勒譜計算：實時接收回波脈沖數(shù)據(jù)并根據(jù)所述系數(shù)矩陣計算回波脈沖數(shù)據(jù)的脈沖多普勒譜，且每次計算時，使用前一個脈沖回波數(shù)據(jù)所計算得到的脈沖多普勒譜，計算當前獲取的脈沖回波數(shù)據(jù)的脈沖多普勒譜。

作為本發(fā)明方法的進一步改進：所述步驟s02中計算脈沖多普勒譜時，具體預先確定參與計算脈沖多普勒的回波脈沖數(shù)量閾值n，當接收到的回波脈沖數(shù)據(jù)超過n時，由獲取的最新回波脈沖數(shù)據(jù)對參與計算脈沖多普勒的回波脈沖數(shù)據(jù)進行更新，對更新后的回波脈沖數(shù)據(jù)計算脈沖多普勒譜。

作為本發(fā)明方法的進一步改進，所述步驟s02的具體步驟為：

s201.定義用于緩存n個回波脈沖數(shù)據(jù)的緩沖矩陣b；

s202.每次接收回波脈沖數(shù)據(jù)時，判斷所接收到的回波脈沖數(shù)據(jù)是否超過回波脈沖數(shù)量閾值n，如果否，轉(zhuǎn)入執(zhí)行步驟s203，否則轉(zhuǎn)入執(zhí)行步驟s34；

s203.根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、所述緩沖矩陣b中緩存數(shù)據(jù)以及所述步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用當前回波脈沖數(shù)據(jù)更新所述緩沖矩陣b；

s204.根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、所述緩沖矩陣b中緩沖數(shù)據(jù)、當前回波脈沖數(shù)據(jù)以及所述步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用最新的回波脈沖數(shù)據(jù)更新所述緩沖矩陣b，以去除所述緩沖矩陣中第一個回波脈沖數(shù)據(jù)。

作為本發(fā)明方法的進一步改進：所述步驟s203中，具體采用下式計算當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜；

sn(m,:)＝sn-1(m,:)+bn(m,n)*c(n,:)；

其中，sn(m,:)為第n個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，sn-1(m,:)為第n-1個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，bn(m,n)為第n個緩沖矩陣中第n個回波脈沖數(shù)據(jù)在第m個距離單元的值，c(n,:)為所述系數(shù)矩陣的第n行矢量；

計算得到第n個回波脈沖數(shù)據(jù)pn(m)在所有m個距離單元的脈沖多普勒譜sn(m,k)后，具體采用下式更新所述緩沖矩陣b；

bn(:,n)＝pn(m)；

其中bn(:,n)為第n個緩沖矩陣中第n個回波脈沖的全部距離單元矢量，m為距離單元的數(shù)量。

作為本發(fā)明方法的進一步改進：具體采用下式計算當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜；

sn(m,:)＝[sn-1(m,:)-bn-1(m,n-1)]*c^*(2,:)+pn(m)*c(n,:)；

其中，sn(m,:)為第n個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，sn-1(m,:)為第n-1個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，bn-1(m,n-1)為第n-1個緩沖矩陣中第n-1個回波脈沖數(shù)據(jù)在第m個距離單元的值，c*(2,:)為所述系數(shù)矩陣的共軛矩陣中第2行矢量，pn(m)為第n個回波脈沖數(shù)據(jù)，c(n,:)為所述系數(shù)矩陣的第n行矢量。

作為本發(fā)明方法的進一步改進：計算得到第n個回波脈沖數(shù)據(jù)pn(m)在所有m個距離單元的脈沖多普勒譜sn(m,k)后，具體采用下式更新所述緩沖矩陣b；

bn(:,1:(n-1))＝bn-1(:,2:n)；

bn(:,n)＝pn(m)；

其中，bn(:,1:(n-1))為第n個緩沖矩陣中第1個到第n-1個回波脈沖的全部距離單元矢量，bn(:,2:n)為第n個緩沖矩陣中第2個到第n個回波脈沖的全部距離單元矢量，bn(:,n)為第n個緩沖矩陣中第n個回波脈沖的全部距離單元矢量。

作為本發(fā)明方法的進一步改進：所述系數(shù)矩陣具體為采用下式計算到的復系數(shù)矩陣c、共軛矩陣c*；

ck×k(p,q)＝exp{-2πj*(p-1)(q-1)/k}；

其中，p＝1,2,…,k，q＝1,2,…,k，j為虛數(shù)單位，conj{·}表示復數(shù)的共軛運算，k為所需計算的點數(shù)。

一種獲取實時脈沖多普勒譜獲取裝置，包括：

系數(shù)矩陣確定模塊，用于根據(jù)所需計算的點數(shù)確定用于計算脈沖多普勒譜的系數(shù)矩陣；

多普勒譜計算模塊，用于實時接收回波脈沖數(shù)據(jù)并根據(jù)所述系數(shù)矩陣計算回波脈沖數(shù)據(jù)的脈沖多普勒譜，且每次計算時，使用前一個脈沖回波數(shù)據(jù)所計算得到的脈沖多普勒譜，計算當前獲取的脈沖回波數(shù)據(jù)的脈沖多普勒譜。

作為本發(fā)明裝置的進一步改進：所述多普勒譜計算模塊計算脈沖多普勒譜時，具體預先確定參與計算脈沖多普勒的回波脈沖數(shù)量閾值n，當接收到的回波脈沖數(shù)據(jù)超過n時，由獲取的最新回波脈沖數(shù)據(jù)對參與計算脈沖多普勒的回波脈沖數(shù)據(jù)進行更新后，對更新后的回波脈沖數(shù)據(jù)計算脈沖多普勒譜。

作為本發(fā)明裝置的進一步改進：所述多普勒譜計算模塊具體包括：

緩沖矩陣定義單元，用于定義用于緩存n個回波脈沖數(shù)據(jù)的緩沖矩陣b；

判斷單元，用于每次接收回波脈沖數(shù)據(jù)時，判斷所接收到的回波脈沖數(shù)據(jù)是否超過回波脈沖數(shù)量閾值n，如果否，轉(zhuǎn)入執(zhí)行第一計算單元，否則轉(zhuǎn)入執(zhí)行第二計算單元；

第一計算單元，用于根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、所述緩沖矩陣b中緩存數(shù)據(jù)以及所述系數(shù)矩陣確定模塊確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用當前回波脈沖數(shù)據(jù)更新所述緩沖矩陣b；

第二計算單元，用于根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、所述緩沖矩陣b中緩沖數(shù)據(jù)、當前回波脈沖數(shù)據(jù)以及所述步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用最新的回波脈沖數(shù)據(jù)更新所述緩沖矩陣b，以去除所述緩沖矩陣中第一個回波脈沖數(shù)據(jù)。

與現(xiàn)有技術(shù)相比，本發(fā)明的優(yōu)點在于：

1)本發(fā)明獲取實時脈沖多普勒譜的方法及裝置，通過根據(jù)所需計算的點數(shù)確定系數(shù)矩陣后計算脈沖多普勒譜，且每次計算脈沖多普勒譜時，僅使用前一個脈沖回波數(shù)據(jù)的脈沖多普勒譜計算當前脈沖回波數(shù)據(jù)的脈沖多普勒譜，即每次計算脈沖多普勒，只需利用前一次脈沖回波數(shù)據(jù)的脈沖多普勒對本次的脈沖多普勒譜進行更新，使得無需對全部的回波數(shù)據(jù)進行計算，可以減少大量的重復計算，計算過程簡單、計算復雜度低，因而多普勒譜的實時性以及獲取效率高；

2)本發(fā)明獲取實時脈沖多普勒譜的方法及裝置，通過設定參與脈沖多普勒計算的回波脈沖數(shù)量閾值，對于超過回波脈沖數(shù)量閾值的回波脈沖數(shù)據(jù)，則由獲取的最新回波脈沖數(shù)據(jù)對參與計算脈沖多普勒的回波脈沖數(shù)據(jù)進行更新，再進行脈沖多普勒譜的計算，使得能夠及時對當前脈沖回波進行更新，以去除舊的脈沖回波，保證能夠?qū)崟r反映目標信息的變換和更新；

3)本發(fā)明獲取實時脈沖多普勒譜的方法及裝置，通過定義用于緩存指定大小的回波脈沖數(shù)據(jù)的緩沖矩陣，基于緩沖矩陣計算回波脈沖多普勒譜，無需等待緩沖矩陣積滿，均能夠隨著得到的回波脈沖數(shù)據(jù)實時地計算更新脈沖多普勒譜，且無需復雜的完整傅里葉變換，僅需少量的加法、乘法運算進行更新，實時性以及獲取效率更高。

附圖說明

圖1是本實施例獲取實時脈沖多普勒譜的方法實現(xiàn)流程示意圖。

圖2是本實施例回波脈沖數(shù)據(jù)矩陣的原理示意圖。

圖3是本發(fā)明具體實施例中獲取實時脈沖多普勒譜的方法實現(xiàn)流程示意圖。

具體實施方式

以下結(jié)合說明書附圖和具體優(yōu)選的實施例對本發(fā)明作進一步描述，但并不因此而限制本發(fā)明的保護范圍。

如圖1、2所示，假定pn(m)表示第n個回波脈沖，每個回波脈沖包含m個距離單元，目標的位置落在第m個距離單元，本實施例獲取實時脈沖多普勒譜的方法步驟包括：

s01.系數(shù)矩陣確定：根據(jù)所需計算的點數(shù)確定用于計算脈沖多普勒譜的系數(shù)矩陣。

在具體實施例中，復系數(shù)矩陣c、共軛矩陣c^*分別為k×k的矩陣，計算公式分別為；

ck×k(p,q)＝exp{-2πj*(p-1)(q-1)/k}(1)

其中，p＝1,2,…,k，q＝1,2,…,k，j為虛數(shù)單位，conj{·}表示復數(shù)的共軛運算，k為所需計算的點數(shù)。

本實施例脈沖多普勒的點數(shù)k可根據(jù)實際需求任意設置，對于不同脈沖多普勒譜計算點數(shù)k，按照上述確定對應的復系數(shù)矩陣c、共軛矩陣c^*，再基于確定的復系數(shù)矩陣c、共軛矩陣c^*即可計算對應的脈沖多普勒譜，不受計算點數(shù)限制，可實現(xiàn)任意點數(shù)的脈沖多普勒譜計算。

s02.多普勒譜計算：實時接收回波脈沖數(shù)據(jù)并根據(jù)系數(shù)矩陣計算回波脈沖數(shù)據(jù)的脈沖多普勒譜，且每次計算時，使用前一個脈沖回波數(shù)據(jù)所計算得到的脈沖多普勒譜，計算當前獲取的脈沖回波數(shù)據(jù)的脈沖多普勒譜。

本實施例通過僅使用前一個脈沖回波數(shù)據(jù)的脈沖多普勒譜，來計算當前脈沖回波數(shù)據(jù)的脈沖多普勒譜，即每次計算脈沖多普勒，只需利用前一次脈沖回波數(shù)據(jù)的脈沖多普勒對本次的脈沖多普勒譜進行更新，使得無需對全部的回波數(shù)據(jù)進行計算，可以減少大量的重復計算，使得計算過程簡單、計算復雜度低。

如圖2所示，對于第m個距離單元，前n個脈沖回波數(shù)據(jù)與前n+1個脈沖回波數(shù)據(jù)之間的脈沖回波數(shù)據(jù)存在重復，若每次計算脈沖多普勒譜時均基于全部脈沖回波數(shù)據(jù)進行計算，則對應的脈沖多普勒存在重復計算部分。本實施例相對于前n個脈沖回波數(shù)據(jù)的脈沖多普勒譜，計算前n+1個脈沖回波數(shù)據(jù)的脈沖多普勒譜時，只需更新第n+1個脈沖回波數(shù)據(jù)引入的新的脈沖數(shù)據(jù)信息，無需計算全部的n+1脈沖回波數(shù)據(jù)的脈沖多普勒。

在具體實施例中，步驟s02中計算脈沖多普勒譜時，具體預先確定參與計算脈沖多普勒的回波脈沖數(shù)量閾值n，當接收到的回波脈沖數(shù)據(jù)超過n時，由獲取的最新回波脈沖數(shù)據(jù)對參與計算脈沖多普勒的回波脈沖數(shù)據(jù)進行更新后，對更新后的回波脈沖數(shù)據(jù)計算脈沖多普勒譜。

考慮到脈沖回波數(shù)據(jù)數(shù)量較多時，若基于全部的脈沖回波數(shù)據(jù)計算脈沖多普勒譜，則無法及時反映目標信息的變換和更新，本實施例通過預先確定參與脈沖多普勒計算的回波脈沖數(shù)量閾值n，使得當接收的脈沖回波數(shù)據(jù)達到閾值n時，能夠及時對當前n個脈沖回波進行更新，以去除舊的脈沖回波，再進行脈沖多普勒譜的計算，保證能夠?qū)崟r反映目標信息的變換和更新。

如圖3所示，在具體實施例中步驟s02的具體步驟為：

s201.定義用于緩存n個回波脈沖數(shù)據(jù)的緩沖矩陣b，緩沖矩陣b大小為m×n。

s202.每次接收回波脈沖數(shù)據(jù)時，判斷所接收到的回波脈沖數(shù)據(jù)是否超過回波脈沖數(shù)量閾值n，如果否，轉(zhuǎn)入執(zhí)行步驟s203，否則轉(zhuǎn)入執(zhí)行步驟s34；

s203.根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、緩沖矩陣b中緩存數(shù)據(jù)以及步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用當前回波脈沖數(shù)據(jù)更新緩沖矩陣b；

s204.根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、緩沖矩陣b中緩沖數(shù)據(jù)、當前回波脈沖數(shù)據(jù)以及步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用最新的回波脈沖數(shù)據(jù)更新緩沖矩陣b，以去除緩沖矩陣中第一個回波脈沖數(shù)據(jù)。

步驟s203中，所接收到的回波脈沖數(shù)據(jù)n未超過回波脈沖數(shù)量閾值n(1≤n≤n)，具體采用式(3)計算當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜；

sn(m,:)＝sn-1(m,:)+bn(m,n)*c(n,:)(3)

其中，sn(m,:)為第n個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，sn-1(m,:)為第n-1個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，bn(m,n)為第n個緩沖矩陣中第n個回波脈沖數(shù)據(jù)在第m個距離單元的值，c(n,:)為系數(shù)矩陣的第n行矢量；

計算得到第n個回波脈沖數(shù)據(jù)pn(m)在所有m個距離單元的脈沖多普勒譜sn(m,k)后，具體采用式(4)更新緩沖矩陣b；

bn(:,n)＝pn(m)(4)

其中bn(:,n)為，m為距離單元的數(shù)量。

如圖2、3所示，對于其第m個距離單元，當?shù)玫降趎(1≤n≤n)個回波脈沖數(shù)據(jù)pn(m)時，對于第m個距離單元，采用上述公式(3)計算得到該距離單元的脈沖多普勒譜sn(m,k)，依次對全部m個距離單元的回波脈沖數(shù)據(jù)進行計算，得到第n(1≤n≤n)個回波脈沖時刻的全部距離單元的脈沖多普勒譜sn(m,k)，并按照公式(4)更新緩沖矩陣b，完成n個內(nèi)的回波脈沖數(shù)據(jù)的多普勒譜計算。

步驟s204中，所接收到的回波脈沖數(shù)據(jù)n超過回波脈沖數(shù)量閾值n(n≥n+1)，具體采用式(5)計算當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜；

sn(m,:)＝[sn-1(m,:)-bn-1(m,n-1)]*c^*(2,:)+pn(m)*c(n,:)(5)

其中，sn(m,:)為第n個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，sn-1(m,:)為第n-1個回波脈沖數(shù)據(jù)在第m個距離單元的脈沖多普勒譜，bn-1(m,n-1)為第n-1個緩沖矩陣中第-1個回波脈沖數(shù)據(jù)在第m個距離單元的值，c*(2,:)為系數(shù)矩陣的共軛矩陣中第2行矢量，pn(m)為第n個回波脈沖數(shù)據(jù)，c(n,:)為系數(shù)矩陣的第n行矢量。

計算得到第n個回波脈沖數(shù)據(jù)pn(m)在所有m個距離單元的脈沖多普勒譜sn(m,k)后，具體采用式(6)、(7)更新緩沖矩陣b；

bn(:,1:(n-1))＝bn-1(:,2:n)(6)

bn(:,n)＝pn(m)(7)

其中，bn(:,1:(n-1))為第n個緩沖矩陣中第1個到第n-1個回波脈沖的全部距離單元矢量，bn(:,2:n)為第n個緩沖矩陣中第2個到第n個回波脈沖的全部距離單元矢量，bn(:,n)為第n個緩沖矩陣中第n個回波脈沖的全部距離單元矢量。由公式(6)、(7)使得去除緩沖矩陣b的第一個回波脈沖數(shù)據(jù)矢量b(:,1)，并將pn(m)寫入緩沖矩陣，從而用最新的回波脈沖數(shù)據(jù)pn(m)更新緩沖矩陣b。

由于第n-1個回波脈沖的脈沖多普勒譜sn-1(m,k)計算時：

bn(m,:)＝bn-1(m,:)-b1(m,1)+pn(:)(8)

則由第n-1個回波脈沖的脈沖多普勒譜sn-1(m,k)計算第n個回波脈沖的脈沖多普勒譜sn(m,k)為：

即第n個回波脈沖的脈沖多普勒譜sn(m,k)可采用上述公式(5)計算得到。

如圖2、3所示，對于其第m個距離單元，當?shù)玫降趎≥n+1個回波脈沖數(shù)據(jù)pn(m)時，采用上述公式(5)計算得到該距離單元的脈沖多普勒譜sn(m,k)，依次對全部m個距離單元的回波脈沖數(shù)據(jù)進行計算后，得到第n(n≥n+1)個回波脈沖時刻的全部距離單元的脈沖多普勒譜sn(m,k)，此時緩沖矩陣b已滿，則按照公式(6)、(7)更新緩沖矩陣b，完成超過n的回波脈沖數(shù)據(jù)的多普勒譜計算。即無需等待緩沖矩陣積滿即可計算得到脈沖多普勒譜，無論緩沖矩陣是否積滿均能夠隨著得到的回波脈沖數(shù)據(jù)實時地計算更新脈沖多普勒譜，實時性好，且無需復雜的完整傅里葉變換，僅需少量的加法和乘法運算，計算量少、計算復雜度低，且計算效率高。

如圖3所示，在具體實施例中，首先指定所需計算的脈沖多普勒譜的點數(shù)k以及開辟m×n大小的緩沖矩陣b，并設定參與計算脈沖多普勒的回波脈沖數(shù)量閾值n；由脈沖多普勒譜的點數(shù)k按照公式(1)、(2)計算得到復系數(shù)矩陣c、共軛矩陣c^*；判斷當前回波脈沖的個數(shù)n是否滿足n>n，如果否，按照式(3)更新計算各個距離單元的脈沖多普勒譜s(m,k)，然后按照公式(4)更新緩沖矩陣b；如果是，按照公式(5)更新計算各個距離單元的脈沖多普勒譜s(m,k)，然后按照公式(6)、(7)更新緩沖矩陣b，最終輸出脈沖多普勒譜s(m,k)，返回執(zhí)行下一次計算，直至退出計算。

本實施例獲取實時脈沖多普勒譜獲取裝置，包括：

系數(shù)矩陣確定模塊，用于根據(jù)所需計算的點數(shù)確定用于計算脈沖多普勒譜的系數(shù)矩陣；

多普勒譜計算模塊，用于實時接收回波脈沖數(shù)據(jù)并根據(jù)系數(shù)矩陣計算回波脈沖數(shù)據(jù)的脈沖多普勒譜，且每次計算時，使用前一個脈沖回波數(shù)據(jù)所計算得到的脈沖多普勒譜，計算當前獲取的脈沖回波數(shù)據(jù)的脈沖多普勒譜。

本實施例中，多普勒譜計算模塊計算脈沖多普勒譜時，具體預先確定參與計算脈沖多普勒的回波脈沖數(shù)量閾值n，當接收到的回波脈沖數(shù)據(jù)超過n時，由獲取的最新回波脈沖數(shù)據(jù)對參與計算脈沖多普勒的回波脈沖數(shù)據(jù)進行更新后，對更新后的回波脈沖數(shù)據(jù)計算脈沖多普勒譜。

本實施例中，多普勒譜計算模塊具體包括：

緩沖矩陣定義單元，用于定義用于緩存n個回波脈沖數(shù)據(jù)的緩沖矩陣b；

判斷單元，用于每次接收回波脈沖數(shù)據(jù)時，判斷所接收到的回波脈沖數(shù)據(jù)是否超過回波脈沖數(shù)量閾值n，如果否，轉(zhuǎn)入執(zhí)行第一計算單元，否則轉(zhuǎn)入執(zhí)行第二計算單元；

第一計算單元，用于根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、緩沖矩陣b中緩存數(shù)據(jù)以及步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用當前回波脈沖數(shù)據(jù)更新緩沖矩陣b；

第二計算單元，用于根據(jù)前一個回波脈沖數(shù)據(jù)所計算得到的脈沖多普勒譜、緩沖矩陣b中緩沖數(shù)據(jù)、當前回波脈沖數(shù)據(jù)以及步驟s01確定的系數(shù)矩陣，計算得到當前回波脈沖數(shù)據(jù)在各距離單元的脈沖多普勒譜，并使用最新的回波脈沖數(shù)據(jù)更新緩沖矩陣b，以去除緩沖矩陣中第一個回波脈沖數(shù)據(jù)。

本實施例中，系數(shù)矩陣確定模塊指定所需計算的脈沖多普勒譜的點數(shù)k，由脈沖多普勒譜的點數(shù)k按照上述公式(1)、(2)計算得到復系數(shù)矩陣c、共軛矩陣c^*；多普勒譜計算模塊中，由判斷單元判斷當前回波脈沖的個數(shù)n是否滿足n>n，如果否，按照上述式(3)更新計算各個距離單元的脈沖多普勒譜s(m,k)，然后按照公式(4)更新緩沖矩陣b；如果是，按照上述公式(5)更新計算各個距離單元的脈沖多普勒譜s(m,k)，然后按照上述公式(6)、(7)更新緩沖矩陣b，最終輸出脈沖多普勒譜s(m,k)，返回執(zhí)行下一次計算，直至退出計算。

上述只是本發(fā)明的較佳實施例，并非對本發(fā)明作任何形式上的限制。雖然本發(fā)明已以較佳實施例揭露如上，然而并非用以限定本發(fā)明。因此，凡是未脫離本發(fā)明技術(shù)方案的內(nèi)容，依據(jù)本發(fā)明技術(shù)實質(zhì)對以上實施例所做的任何簡單修改、等同變化及修飾，均應落在本發(fā)明技術(shù)方案保護的范圍內(nèi)。