DSP多普勒雷達測速測距

1、精選優(yōu)質文檔-傾情為你奉上DSP實驗課大作業(yè)設計一 實驗目的在DSP上實現(xiàn)線性調頻信號的脈沖壓縮、動目標顯示（MTI）和動目標檢測(MTD)，并將結果與MATLAB上的結果進行誤差仿真。二 實驗內容2.1 MATLAB仿真 設定帶寬、脈寬、采樣率、脈沖重復頻率，用MATLAB產生16個脈沖的LFM，每個脈沖有4個目標（靜止，低速，高速），依次做2.1.1 脈壓2.1.2 相鄰2脈沖做MTI，產生15個脈沖2.1.3 16個脈沖到齊后，做MTD，輸出16個多普勒通道2.2 DSP實現(xiàn)將MATLAB產生的信號，在visual dsp中做脈壓，MTI、MTD，并將結果與MATLAB作比較。三 實驗原

2、理3.1 脈沖壓縮原理及線性調頻信號雷達中的顯著矛盾是：雷達作用距離和距離分辨率之間的矛盾以及距離分辨率和速度分辨率之間的矛盾。雷達的距離分辨率取決于信號帶寬。在普通脈沖雷達中，雷達信號的時寬帶寬積為一常量（約為1），因此不能兼顧距離分辨率和速度分辨力兩項指標。脈沖壓縮（PC）采用寬脈沖發(fā)射以提高發(fā)射的平均功率，保證足夠的最大作用距離，而在接收時則采用相應的脈沖壓縮法獲得窄脈沖，以提高距離分辨率，因而能較好地解決作用距離和分辨能力之間的矛盾。一個理想的脈沖壓縮系統(tǒng)，應該是一個匹配濾波系統(tǒng)。它要求發(fā)射信號具有非線性的相位譜，并使其包絡接近矩形；要求壓縮網絡的頻率特性（包括幅頻特性和相頻特性）與發(fā)

3、射脈沖信號頻譜（包括幅度譜和相位譜）實現(xiàn)完全的匹配。脈沖壓縮按信號的調制規(guī)律（調頻或調相）分類，可分為以下四種：（1）線性調頻脈沖壓縮（2）非線性調頻脈沖壓縮（3）相位編碼脈沖壓縮（4）時間頻率編碼脈沖壓縮本實驗采用的是線性調頻脈沖壓縮。線性調頻信號是指頻率隨時間的變化而線性改變的信號。線性調頻可以同時保留連續(xù)信號和脈沖的特性，并且可以獲得較大的壓縮比，有著良好的距離分辨率和徑向速度分辨率，所以將線性調頻信號作為雷達系統(tǒng)中一種常用的脈沖壓縮信號。接收機輸入端的回波信號是經過調制的寬脈沖，所以在接收機中應該設置一個與發(fā)射信號頻率匹配的濾波器，使回波信號變成窄脈沖，同時實現(xiàn)了寬脈沖的能量和窄脈沖的

4、分辨能力。解決了雷達發(fā)射能量及分辨率之間的矛盾。匹配濾波器是指輸出信噪比最大準則下的最佳線性濾波器。根據匹配理論，匹配濾波器的傳輸特性：其中，為幅度歸一化常數(shù)，為信號的復共軛。傳輸特性還可用它的沖激響應來表示（時域表示）：。3.2 MTI（動目標顯示）原理在雷達終端的顯示器，當雜波與運動的目標在顯示器上同時出現(xiàn)時，觀察將會變得困難，特別是在強雜波時，弱的目標信號將會完全被淹沒。由于目標回波的多普勒頻移遠大于雜波的多普勒頻移，從而可在頻域上區(qū)分目標與雜波。MTI的設計主要是合理的選擇濾波器，動目標顯示濾波器（MTI）利用運動目標回波和雜波在頻譜上的區(qū)別，有效地抑制雜波而提取信號。由下圖可以看出回

5、波信號（含雜波與目標）：圖3.1 回波信號雜波主要集中在0頻附近，在頻率為0處，所設計的濾波器頻率響應應有凹口。以使地物雜波在通過MTI濾波器后將受到很大的抑制，對于每一個脈沖，其回波信號中的雜波分量都基本相同，而運動目標產生多普勒頻移。對連續(xù)脈沖對的回波進行相減，就可以完全對消雜波分量，而目標信號由于其相位的改變，通常不會完全被對消掉。常用的MTI濾波器為一次相消器（二脈沖對消）。在同一距離上等效一個有效的采樣間隔的離散時間序列。其時域方程為：，傳輸函數(shù)為：，在Z=1處有單零點，頻率響應為：。其在零頻有一凹口，如圖所示：圖3.3 多普勒濾波器頻響此濾波器能夠用來抑制噪聲。但同時把靜態(tài)目標也給

6、對消掉，也使得運動目標的譜分量有部分的衰減，因此用MTI一次相消器檢測不出靜止目標。3.3 MTD（動目標檢測）原理MTD與MTI同屬雷達信號處理的頻域處理范疇，但一般意義上來說，MTD是MTI的改進或更有效的頻域處理技術。廣義上地講，MTD處理又是多普勒處理的一種特殊形式與傳統(tǒng)的MTI相比，MTD改善了濾波器的頻率特性，使之更接近于最佳匹配線性濾波，以提高改善因子；能夠檢測強地物雜波的低速目標甚至切向飛行的大目標。MTD的核心是線性MTI加窄帶多普勒濾波器組。根據最佳線性濾波理論，在雜波背景下檢測運動目標回波，除了雜波抑制濾波器外，還應串接有對脈沖串信號匹配的濾波器，實際工作中，多普勒頻移不

7、能預知，需要采用一組相鄰且部分重疊的濾波器組覆蓋整個多普勒頻率范圍，這就是窄帶多普勒濾波器組。輸入N個脈沖的橫向濾波器組的實現(xiàn)是由N個輸出的橫向濾波器（N個脈沖和N-1根遲延線）經過各脈沖不同的加權并求和后形成的。其實現(xiàn)方法下圖所示： 圖3.4 橫向濾波器 圖3.5 N=16時濾波器頻響該濾波器的頻率覆蓋范圍為0到，為雷達工作重復頻率。氣象雜波(如云雨等)和箔條雜波受氣流和風力的影響，會相對雷達而動，其頻譜中心不是固定在0頻附近，而是在某個頻率區(qū)間內變化的，抑制此類雜波用普通的MTI濾波器是不行的，而MTD濾波器則可以抑制此類雜波。四 實驗步驟以下是該實驗中設定的幾個參數(shù)BandWidth=2

8、.0e6-帶寬TimeWidth=42.0e-6-脈寬Fs=2.0e6-采樣率PRT=240e-6-脈沖重復周期TargetDistance=3000 8025 8025 25600 -目標距離TargetVelocity=50 -100 0 283-目標速度%其中X=0; Y=8; Z=3;對應學號的后三位假設接收到的回波數(shù)是16個，噪聲為高斯隨機噪聲。4.1 在MATLAB中產生線性調頻信號。4.2 根據目標距離得出目標回波在時域的延遲量，根據目標速度得出多普勒相移，從而在MATLAB中產生4個目標的16個回波串（接收到的回波含噪聲）。4.3 由匹配濾波理論產生對應于目標回波的濾波系數(shù)（脈

9、壓系數(shù)），在時域中做線性卷積，實現(xiàn)匹配濾波（時域脈壓）；在頻域中做回波信號和脈壓系數(shù)的FFT，點乘后作逆FFT，實現(xiàn)頻域脈壓。兩者進行比較，討論其差別。4.4 對16個去暫態(tài)點后的脈沖串按接收順序進行排列，用一次相消器（一種濾波方式）實現(xiàn)MTI。4.5 做16通道的FFT，實現(xiàn)MTD。4.6 在DSP中對MATLAB產生的回波數(shù)據和脈壓系數(shù)進行處理，實現(xiàn)頻域脈壓。導入DSP的回波數(shù)據為時域數(shù)據，而脈壓系數(shù)為頻域數(shù)據。將導入DSP的時域回波數(shù)據進行一次FFT變換到頻域，然后將其與頻域脈壓系數(shù)進行點積，得到頻域脈壓結果。對該結果再做一次逆FFT，將頻域轉換成時域。在這一步中需要調用庫函數(shù)fft_f

10、lp32.asm。該子程序可實現(xiàn)8192點復數(shù)的FFT功能。由于C語言中無法實現(xiàn)復數(shù)運算，因此，對8192個復數(shù)按照實部虛部交替的順序進行重排列，用長度為16384的數(shù)組來存放時域回波數(shù)據。頻域相乘后，做乘積結果的逆FFT，得到脈壓結果。做逆FFT，仍需調用庫函數(shù)fft_flp32.asm，此時要通過FFT子程序實現(xiàn)逆FFT的功能，要對頻域的數(shù)據進行處理，才能達到這一目的。4.7 對脈壓后的數(shù)據按照脈沖號重排，相鄰序列的數(shù)據相減（滑動對消），實現(xiàn)MTI。4.8 調用子程序fft_16.asm，做16通道FFT，實現(xiàn)MTD。入口參數(shù)為16通道的脈壓數(shù)據。五 實驗結果及討論本部分將詳細分析實驗得到

11、的數(shù)據、圖像、誤差、產生速度模糊的原因以及脈壓頻域、時域執(zhí)行周期。5.1MATLAB中時域脈壓和頻域脈壓結果的比較: 圖5.1 總回波信號，時域脈壓，頻域脈壓，時域脈壓和頻域脈壓的差別從圖中可以看出，時頻域脈壓結果差別很小，相對誤差停留在10-14數(shù)量級上。其實，二者結果應該是相等的。由于MATLAB與DSP軟件處理平臺及編程函數(shù)的計算精度和中間結果處理的差異，才出現(xiàn)以上很小的誤差。由數(shù)字信號處理理論，M點時域離散信號與N點時域離散信號做線性卷積后的信號長度為N+M-1。設為數(shù)字信號，長度分別為M和N，如果的長度，則?；夭〝?shù)據echo長為7680點，脈壓系數(shù)coeff長為84點，N+M-1=7

12、680+84-1=7763，為提高FFT速度，取L=8192，對回波數(shù)據、脈壓系數(shù)的FFT乘積后的結果做IFFT就可以實現(xiàn)脈沖壓縮。經統(tǒng)計，在DSP上頻域脈壓處理（從FFT到IFFT）的時鐘數(shù)為。設兩個卷積的序列長度分別為M 和N，且N>M，則卷積所用的復乘數(shù)為N×M，而N點FFT的運算量為次復數(shù)乘法運算和次復數(shù)加法運算。當序列長度N、M很大時，頻域的總運算量將大大低于時域的運算量。在本實驗中，回波數(shù)據的長度N=7680，脈壓系數(shù)長度M=84，時域脈壓的運算量為N×M=7680×84+83=次復數(shù)加法運算和次復數(shù)乘法運算。由于L>N+M-1，將x(n)

13、，h(n)分別補L-N和L-M個零變成兩個長度均為L（8192）的序列，分別做L點FFT，點積后再做逆FFT，可以實現(xiàn)頻域脈壓。頻域脈壓的8192點FFT和IFFT的運算量為2×0.5×8192×=次復數(shù)乘法運算和0.5×2×8192×=次復數(shù)加法運算。可以看出，用頻域脈壓節(jié)省了很大的的運算量，當序列長度很大時，用FFT進行運算的優(yōu)勢更明顯。5.2 MATLAB與DSP的結果比較:5.2 MATLAB與DSP處理結果比較由圖5.2可以看出，MATLAB與DSP之間處理數(shù)據的誤差很小，脈壓結果絕對誤差范圍在數(shù)量級上，MTI、MTD絕對誤

14、差范圍分別在和上。由于兩種工具處理數(shù)據的精度不一樣，MATLAB用CPU處理數(shù)據，DSP用DSP內核處理數(shù)據。所以結果會有一定誤差，這個誤差在可以接受的范圍內。六、實驗思考題及心得6.1 脈沖壓縮分析脈沖壓縮提高了雷達分辨率，從信號幅度上分析，第一個目標由于距離比較近，雷達發(fā)射與接收是分開進行的，當回波到達是雷達還處于發(fā)射狀態(tài)，雷達處于閉鎖期不能接收回波，從而第一個目標的回波幅度很小。與此同時，由于閉鎖期的影響，回波信號的帶寬減小，而脈壓后的信號時寬帶寬積是定值，故時寬變大，導致旁瓣變大。第二與第四個目標的距離相同，在時域上兩個信號的回波進行疊加，所以回波幅度比較強，但是由于兩個目標回波相位上

15、存在角度，第四個目標存在多普勒頻移，在不同的脈沖重復周期內兩個信號的角度不同，矢量相加的結果則信號幅度也不同，所以幅度存在著波動。 6.2 MTI分析根據一次相消得頻率響應，可以知道不同頻率的幅頻響應的差異，當處在中間位置時響應最大，通過濾波器后信號幅度也就最強。根據頻率響應的不同也就出現(xiàn)MTI結果圖中的三個梯度的變化。在MTI中不同速度目標幅度有變化是因為有不同的多普勒頻移，可以說目標的幅度被多普勒頻移信號調制了。可以從另外一個角度理解，當采用連續(xù)波雷達測速時經過差頻處理后獲得的信號是一個以多普勒頻移為中心頻率的單載頻信號，當采用脈沖信號時，相當于是以PRF為采樣頻率對單載頻信號進行一個采樣

16、，幅度自然就會產生變化。6.3 MTD速度/多普勒通道的含義MTD中不同的多普勒通道是匹配與不同多普勒頻移的匹配濾波器。當目標速度產生的多普勒移處在哪個通道的范圍內是哪個通道就有信號輸出，其他通道不匹配則無輸出，測速分辨力，假設目標速度都是正的時候，當信號處在第n通道時（通道數(shù)為116），其速度為。 實際上目標的速度有正負的時候多普勒頻移時在0頻兩側的，測速范圍為。在一個重復周期內，根據目標處在的采樣點數(shù)n可以得出目標的距離為。距離分辨力 。² 第一個目標X=37,Y=2,則：距離=2775m，速度= 49.76m/s² 第二個目標X=107,Y=-4,則：距離=8025m

17、，速度=99.52m/s² 第三個目標X=107，Y=0,則：距離=8025m，速度= 0m/s² 第四個目標X=341，Y=-5，則：距離=25575m，速度= 124.4m/s與設定的參數(shù)進行對比：TargetDistance=2850 8025 8025 25600-目標距離TargetVelocity=50 -100 0 283- -目標速度對比可知，第一、二、三個目標在測量的分辨力內達到了準確探測目標的作用，而第四個目標能準備探測到距離，但速度出現(xiàn)了模糊。6.4 速度模糊及解模糊6.4.1 速度模糊及其原因脈沖多普勒雷達以脈沖周期為的間隔采集目標物回波信號，在兩次

18、采樣之間目標物相對雷達運動造成信號中相位的變化超過±180°時，雷達無法確認其準確的相位變化量，也無法確定其準確的多普勒速度的現(xiàn)象。因此，限制了測速及測距的范圍。雷達可測的有效的最大距離為=0.5×240×300=36000米，可測的有效的最大速度為=398.1250米/秒，但本實驗中為了方便觀察負向宿舍，用fftshift將零頻搬移到中間，故可測的速度范圍是-199199m/s。而第四個目標的速度為283m/s,超出了最大可測的速度范圍，故出現(xiàn)了速度模糊。6.4.2 解速度模糊工程上通常用多重頻的方法來解模糊，即用兩個不同的脈沖重復頻率的信號來測量目標的速度。假設這兩個不同的重頻分別為PTF1和PTF2，測得的目標的多普勒頻率為Fd1,Fd2，目標實際的多普勒頻率為Fd，則m=0,1,2. n=0,1,2.可知，F(xiàn)d與m,n都是線性關系，在坐標圖中分別將其表示出來，兩直線交點即為該目標的真實多普勒頻率。6.5 頻譜泄漏由于進行DFT運算是存在頻譜泄漏的問題，所以進過MTD處理之后第四個目標存在速度泄漏的問題。泄漏的原因是在進行DFT處理的過程中在時域是采用截斷函數(shù)，若用N點DFT對一復數(shù)序列進行變換，則只有當該序列頻率等于DFT頻率間隔的1/(NTR)的整數(shù)倍時，在fr 內才會