雷達(dá)線性調(diào)頻脈沖壓縮地原理與其MATLAB仿真

1、WORD格式PAGE1 / NUMPAGES18LFM脈沖壓縮雷達(dá)仿真線性調(diào)頻（LFM）脈沖壓縮雷達(dá)仿真一雷達(dá)工作原理雷達(dá)是Radar（RAdioDetectionAndRanging）的音譯詞，意為“無(wú)線電檢測(cè)和測(cè)距”，即利用無(wú)線電波來(lái)檢測(cè)目標(biāo)并測(cè)定目標(biāo)的位置，這也是雷達(dá)設(shè)備在最初階段的功能。典型的雷達(dá)系統(tǒng)如圖1.1，它主要由發(fā)射機(jī)，天線，接收機(jī)，數(shù)據(jù)處理，定時(shí)控制，顯示等設(shè)備組成。利用雷達(dá)可以獲知目標(biāo)的有無(wú)，目標(biāo)斜距，目標(biāo)角位置，目標(biāo)相對(duì)速度等?，F(xiàn)代高分辨雷達(dá)擴(kuò)展了原始雷達(dá)概念，使它具有對(duì)運(yùn)動(dòng)目標(biāo)(飛機(jī)，導(dǎo)彈等)和區(qū)域目標(biāo)(地面等)成像和識(shí)別的能力。雷達(dá)的應(yīng)用越來(lái)越廣泛。圖1.1：簡(jiǎn)單脈沖

2、雷達(dá)系統(tǒng)框圖雷達(dá)發(fā)射機(jī)的任務(wù)是產(chǎn)生符合要求的雷達(dá)波形（RadarWaveform），然后經(jīng)饋線和收發(fā)開(kāi)關(guān)由發(fā)射天線輻射出去，遇到目標(biāo)后，電磁波一部分反射，經(jīng)接收天線和收發(fā)開(kāi)關(guān)由接收機(jī)接收，對(duì)雷達(dá)回波信號(hào)做適當(dāng)?shù)奶幚砭涂梢垣@知目標(biāo)的相關(guān)信息。假設(shè)理想點(diǎn)目標(biāo)與雷達(dá)的相對(duì)距離為R，為了探測(cè)這個(gè)目標(biāo)，雷達(dá)發(fā)射信號(hào)s(t),電磁波以光速C向四周傳播，經(jīng)過(guò)時(shí)間RC后電磁波到達(dá)目標(biāo)，照射到目標(biāo)上的電磁波可寫(xiě)成：Rs(t)C。電磁波與目標(biāo)相互作用，一部分電磁波被目標(biāo)散射，被反射的電磁波為Rs(t)C，其中為目標(biāo)的雷達(dá)散射截面（RadarCrossSection,簡(jiǎn)稱RCS），反映目標(biāo)對(duì)R電磁波的散射能力。再經(jīng)

3、過(guò)時(shí)間RC后，被雷達(dá)接收天線接收的信號(hào)為(2)stC。如果將雷達(dá)天線和目標(biāo)看作一個(gè)系統(tǒng)，便得到如圖1.2的等效，而且這是一個(gè)LTI（線性時(shí)不變）系統(tǒng)。圖1.2：雷達(dá)等效于LTI系統(tǒng)等效LTI系統(tǒng)的沖擊響應(yīng)可寫(xiě)成:Mh(t)(t)(1.1)iii1-1-LFM脈沖壓縮雷達(dá)仿真M表示目標(biāo)的個(gè)數(shù)，i是目標(biāo)散射特性，i是光速在雷達(dá)與目標(biāo)之間往返一次的時(shí)間，i2Ric(1.2)式中，R為第i個(gè)目標(biāo)與雷達(dá)的相對(duì)距離。i雷達(dá)發(fā)射信號(hào)s(t)經(jīng)過(guò)該LTI系統(tǒng)，得輸出信號(hào)(即雷達(dá)的回波信號(hào))()st：rMMs(t)s(t)*h(t)s(t)*(t)s(t)(1.3)riiiii1i1那么，怎樣從雷達(dá)回波信號(hào)()

4、st提取出表征目標(biāo)特性的i(表征相對(duì)距離)和i(表征目標(biāo)反r射特性)呢？常用的方法是讓()st通過(guò)雷達(dá)發(fā)射信號(hào)s(t)的匹配濾波器，如圖1.3。r圖1.3：雷達(dá)回波信號(hào)處理st的匹配濾波器hr(t)為：()*h(t)s(t)(1.4)r于是，*s(t)s(t)*h(t)s(t)*s(t)*h(t)（1.5）orr對(duì)上式進(jìn)行傅立葉變換：*S(jw)S(jw)S(jw)H(jw)o2|S(jw)|H(jw)(1.6)如果選取合適的s(t)，使它的幅頻特性|S(jw)|為常數(shù)，那么1.6式可寫(xiě)為：S(jw)kH(jw)(1.7)oM其傅立葉反變換為：s(t)kh(t)k(t)(1.8)oiii1s(

5、t)中包含目標(biāo)的特征信息i和i。從so(t)中可以得到目標(biāo)的個(gè)數(shù)M和每個(gè)目標(biāo)相對(duì)o雷達(dá)的距離：cR（1.9）ii2這也是線性調(diào)頻（LFM）脈沖壓縮雷達(dá)的工作原理。-2-LFM脈沖壓縮雷達(dá)仿真二線性調(diào)頻（LFM）信號(hào)脈沖壓縮雷達(dá)能同時(shí)提高雷達(dá)的作用距離和距離分辨率。這種體制采用寬脈沖發(fā)射以提高發(fā)射的平均功率，保證足夠大的作用距離；而接受時(shí)采用相應(yīng)的脈沖壓縮算法獲得窄脈沖，以提高距離分辨率，較好的解決雷達(dá)作用距離與距離分辨率之間的矛盾。脈沖壓縮雷達(dá)最常見(jiàn)的調(diào)制信號(hào)是線性調(diào)頻（LinearFrequencyModulation）信號(hào),接收時(shí)采用匹配濾波器（MatchedFilter）壓縮脈沖。LFM

6、信號(hào)(也稱Chirp信號(hào))的數(shù)學(xué)表達(dá)式為：tK2j2(ftt)cs(t)rect()Te2(2.1)式中tf為載波頻率，rect()cT為矩形信號(hào)，t11trect()TT0,elsewise（2.2）KBTcTT，是調(diào)頻斜率，于是，信號(hào)的瞬時(shí)頻率為(22)fKtt，如圖2.1圖2.1典型的chirp信號(hào)（a）up-chirp(K0)（b）down-chirp(K0)將2.1式中的up-chirp信號(hào)重寫(xiě)為：j2fts(t)S(t)e（2.3）c式中，tS(t)rect()eTjKt2（2.4）是信號(hào)s(t)的復(fù)包絡(luò)。由傅立葉變換性質(zhì)，S(t)與s(t)具有相同的幅頻特性，只是中心頻率不同而以

7、，因此，Matlab仿真時(shí)，只需考慮S(t)。以下Matlab程序產(chǎn)生2.4式的chirp信號(hào)，并作出其時(shí)域波形和幅頻特性，如圖2.2。-3-LFM脈沖壓縮雷達(dá)仿真%demoofchirpsignalT=10e-6;%pulseduration10usB=30e6;%chirpfrequencymodulationbandwidth30MHzK=B/T;%chirpslopeFs=2*B;Ts=1/Fs;%samplingfrequencyandsamplespacingN=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.2);%generatechir

8、psignalsubplot(211)plot(t*1e6,real(St);xlabel(Timeinusec);title(Realpartofchirpsignal);gridon;axistight;subplot(212)freq=linspace(-Fs/2,Fs/2,N);plot(freq*1e-6,fftshift(abs(fft(St);xlabel(FrequencyinMHz);title(Magnitudespectrumofchirpsignal);gridon;axistight;仿真結(jié)果顯示：圖2.2：LFM信號(hào)的時(shí)域波形和幅頻特性-4-LFM脈沖壓縮雷達(dá)仿真三

9、LFM脈沖的匹配濾波信號(hào)s(t)的匹配濾波器的時(shí)域脈沖響應(yīng)為：*h(t)s(tt)（3.1）0t是使濾波器物理可實(shí)現(xiàn)所附加的時(shí)延。理論分析時(shí)，可令0t0，重寫(xiě)3.1式，0*h(t)s(t)（3.2）將2.1式代入3.2式得:t2j2ftjKth(t)rect()eecT(3.3)圖3.1：LFM信號(hào)的匹配濾波如圖3.1,s(t)經(jīng)過(guò)系統(tǒng)h(t)得輸出信號(hào)()st，os(t)s(t)*h(t)os(u)h(tu)duh(u)s(tu)duutu22()22()jKujfujKtujftuerect()ecerect()ecduTT當(dāng)0tT時(shí),T22jKtj2Ktus(t)eedu0tT2j2Kt

10、uTe222jKtjfteecj2KttT2(3.4)sinK(Tt)tKtj2ftec當(dāng)Tt0時(shí),tT22jKtj2Ktus(t)eedu0T2jKtuT2tejKt2j2ft2eecj2KtT2(3.5)sinK(Tt)tKtej2ftc-5-LFM脈沖壓縮雷達(dá)仿真合并3.4和3.5兩式：tsinKT(1)ttTj2fts(t)Trect()ec0KTt2T（3.6）3.6式即為L(zhǎng)FM脈沖信號(hào)經(jīng)匹配濾波器得輸出,它是一固定載頻fc的信號(hào)。當(dāng)tT時(shí)，包絡(luò)近似為辛克（sinc）函數(shù)。ttS0(t)TSa(KTt)rect()TSa(Bt)rect()2T2T（3.7）圖3.2：匹配濾波的輸出信

11、號(hào)如圖3.2，當(dāng)Bt時(shí)，t1B為其第一零點(diǎn)坐標(biāo)；當(dāng)Bt時(shí)，2t12B，習(xí)慣上，將此時(shí)的脈沖寬度定義為壓縮脈沖寬度。1122BB(3.8)LFM信號(hào)的壓縮前脈沖寬度T和壓縮后的脈沖寬度之比通常稱為壓縮比D，TDTB（3.9）3.9式表明，壓縮比也就是LFM信號(hào)的時(shí)寬頻寬積。由2.1,3.3,3.6式，s(t),h(t),so(t)均為復(fù)信號(hào)形式，Matab仿真時(shí)，只需考慮它們的復(fù)包絡(luò)S(t),H(t),So(t)。以下Matlab程序段仿真了圖3.1所示的過(guò)程，并將仿真結(jié)果和理論進(jìn)行對(duì)照。%demoofchirpsignalaftermatchedfilterT=10e-6;%pulsedura

12、tion10usB=30e6;%chirpfrequencymodulationbandwidth30MHzK=B/T;%chirpslopeFs=10*B;Ts=1/Fs;%samplingfrequencyandsamplespacingN=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.2);%chirpsignalHt=exp(-j*pi*K*t.2);%matchedfilterSot=conv(St,Ht);%chirpsignalaftermatchedfiltersubplot(211)L=2*N-1;-6-LFM脈沖壓縮雷達(dá)仿真t1=l

13、inspace(-T,T,L);Z=abs(Sot);Z=Z/max(Z);%normalizeZ=20*log10(Z+1e-6);Z1=abs(sinc(B.*t1);%sincfunctionZ1=20*log10(Z1+1e-6);t1=t1*B;plot(t1,Z,t1,Z1,r.);axis(-15,15,-50,inf);gridon;legend(emulational,sinc);xlabel(TimeinsectimesitB);ylabel(Amplitude,dB);title(Chirpsignalaftermatchedfilter);subplot(212)%zo

14、omN0=3*Fs/B;t2=-N0*Ts:Ts:N0*Ts;t2=B*t2;plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),r.);axis(-inf,inf,-50,inf);gridon;set(gca,Ytick,-13.4,-4,0,Xtick,-3,-2,-1,-0.5,0,0.5,1,2,3);xlabel(TimeinsectimesitB);ylabel(Amplitude,dB);title(Chirpsignalaftermatchedfilter(Zoom);仿真結(jié)果如圖3.3：圖3.3：Chirp信號(hào)的匹配濾波-7-LFM脈沖壓縮雷達(dá)仿真圖

15、3.3中，時(shí)間軸進(jìn)行了歸一化，（t/(1/B)tB）。圖中反映出理論與仿真結(jié)果吻合良好。第一零點(diǎn)出現(xiàn)在1（即1B）處，此時(shí)相對(duì)幅度-13.4dB。壓縮后的脈沖寬度近似為1B（12B），此時(shí)相對(duì)幅度-4dB,這理論分析（圖3.2）一致。上面只是對(duì)各個(gè)信號(hào)復(fù)包絡(luò)的仿真，實(shí)際雷達(dá)系統(tǒng)中，LFM脈沖的處理過(guò)程如圖3.4。圖3.4：LFM信號(hào)的接收處理過(guò)程雷達(dá)回波信號(hào)()st（1.4式）經(jīng)過(guò)正交解調(diào)后，得到基帶信號(hào)，再經(jīng)過(guò)匹配濾波脈沖壓縮后r就可以作出判決。正交解調(diào)原理如圖3.5，雷達(dá)回波信號(hào)經(jīng)正交解調(diào)后得兩路相互正交的信號(hào)I(t)和Q(t)。一種數(shù)字方法處理的的匹配濾波原理如圖3.6。圖3.5：正交解

16、調(diào)原理圖3.6：一種脈沖壓縮雷達(dá)的數(shù)字處理方式四：Matlab仿真結(jié)果（1）任務(wù)：對(duì)以下雷達(dá)系統(tǒng)仿真。雷達(dá)發(fā)射信號(hào)參數(shù)：幅度：1.0信號(hào)波形：線性調(diào)頻信號(hào)頻帶寬度：30兆赫茲（30MHz）脈沖寬度：10微妙（20us）9Hz）中心頻率：1GHz（10雷達(dá)接收方式：正交解調(diào)接收-8-LFM脈沖壓縮雷達(dá)仿真距離門(mén)：10Km15Km目標(biāo)：Tar1：10.5KmTar2：11KmTar3：12KmTar4：12Km5mTar5：13KmTar6：13Km2m（2）系統(tǒng)模型：結(jié)合以上分析，用Matlab仿真雷達(dá)發(fā)射信號(hào)，回波信號(hào)，和壓縮后的信號(hào)的復(fù)包絡(luò)特性，其載頻不予考慮（實(shí)際中需加調(diào)制和正交解調(diào)環(huán)節(jié)）

17、，仿真信號(hào)與系統(tǒng)模型如圖4.1。圖4.1：雷達(dá)仿真等效信號(hào)與系統(tǒng)模型（3）線性調(diào)頻脈沖壓縮雷達(dá)仿真程序LFM_radar仿真程序模擬產(chǎn)生理想點(diǎn)目標(biāo)的回波，并采用頻域相關(guān)方法（以便利用FFT）實(shí)現(xiàn)脈沖壓縮。函數(shù)LFM_radar的參數(shù)意義如下：T：chirp信號(hào)的持續(xù)脈寬；B：chirp信號(hào)的調(diào)頻帶寬；Rmin：觀測(cè)目標(biāo)距雷達(dá)的最近位置；Rmax：觀測(cè)目標(biāo)距雷達(dá)的最遠(yuǎn)位置；R：一維數(shù)組，數(shù)組值表示每個(gè)目標(biāo)相對(duì)雷達(dá)的斜距；RCS：一維數(shù)組，數(shù)組值表示每個(gè)目標(biāo)的雷達(dá)散射截面。在Matlab指令窗中鍵入：LFM_radar(10e-6,30e6,10000,15000,10500,11000,1200

18、0,12005,13000,13002,1,1,1,1,1,1)得到的仿真結(jié)果如圖4.2。（4）分辨率(Resolution)仿真改變兩目標(biāo)的相對(duì)位置，可以分析線性調(diào)頻脈沖壓縮雷達(dá)的分辨率。仿真程序默認(rèn)參數(shù)的距離分辨率為：R8C31062B230105m(4.1)圖4.3為分辨率仿真結(jié)果，可做如下解釋：(a)圖為單點(diǎn)目標(biāo)壓縮候的波形；(b)圖中，兩目標(biāo)相距2m，小于R，因而不能分辨；(c)圖中，兩目標(biāo)相距5m，等于R，實(shí)際上是兩目標(biāo)的輸出sinc包絡(luò)疊加，可以看到他們的副瓣相互抵消；(d)(h)圖中，兩目標(biāo)距離大于雷達(dá)的距離分辨率，可以觀察出，它們的主瓣變寬，直至能分辨出兩目標(biāo)。-9-LFM脈

19、沖壓縮雷達(dá)仿真圖4.2：仿真結(jié)果-10-LFM脈沖壓縮雷達(dá)仿真圖4.3：線性調(diào)頻脈沖壓縮雷達(dá)分辨率仿真-11-LFM脈沖壓縮雷達(dá)仿真附錄：LFM_radar.m%demoofLFMpulseradar%=functionLFM_radar(T,B,Rmin,Rmax,R,RCS)ifnargin=0T=10e-6;%pulseduration10usB=30e6;%chirpfrequencymodulationbandwidth30MHzRmin=10000;Rmax=15000;%rangebinR=10500,11000,12000,12008,13000,13002;%position

20、ofidealpointtargetsRCS=111111;%radarcrosssectionend%=%ParameterC=3e8;%propagationspeedK=B/T;%chirpslopeRwid=Rmax-Rmin;%receivewindowinmeterTwid=2*Rwid/C;%receivewindowinsecondFs=5*B;Ts=1/Fs;%samplingfrequencyandsamplingspacingNwid=ceil(Twid/Ts);%receivewindowinnumber%=%Gneratetheechot=linspace(2*Rmin/C,2*Rmax/C,Nwid);%receivewindow%openwindowwhent=2*Rmin/C%closewindowwhent=2*Rmax/CM=length(R);%numberoftargetstd=ones(M,1)*t-2*R/C*ones(1,Nwid);Srt=RCS*(exp(j*pi*K*td.2).*(abs(td)T/2);%radarechofrompointtargets%=%Digtalprocessing