雷达原理---脉冲雷达和连续波雷达
一、雷達分類
雷達按雷達信號的形式可分為脈沖雷達和連續波雷達。
脈沖雷達:
連續波雷達:
兩種情況:
- 大信號干擾使得接收機壓縮增益或者出現飽和,甚至造成接收機阻塞,通常可以通過將收發天線進行物理隔離來解決;
- 發射信號的邊帶噪聲將微弱的回波信號淹沒,對接收機的目標檢測造成影響。
邊帶(sideband):調制后的信號,在中心載頻的上下兩側各產生一個頻帶,稱作邊帶。邊帶的帶寬是由所使用的調制信號的帶寬和調制方式決定的。
解決方法:
??直接的信號泄露通常可以采用收發天線隔離和頻率分離相結合的方法得以解決。
??在多普勒導航器中,多普勒頻移可以提供足夠的頻率間隔,以保證發射信號不對接收機進行干擾。
脈沖體制雷達的優點:
??對于機載雷達,每個發射機不可避免的都會產生噪聲,并且調制到發射機的輸出,產生調制的邊帶噪聲,覆蓋了發射頻率左右很寬的頻帶,盡管這些邊帶噪聲功率很小,但是依然比來自目標的回波信號強很多個數量級。為了防止發射機邊帶噪聲干擾接收信號,必須將接收機與發射機隔離,采用獨立的發射機和接收機,并且發射機和接收機采用各種獨立的天線,從而實現發射機和接收機的隔離。地面和艦載連續波雷達就是如此。然而機載雷達因為空間受限,通常收發要共用一副天線,因此發射機的邊帶噪聲不可避免的通過天線進入接收機。脈沖體制雷達可以有效的避免出現發射機干擾接收機的問題。
- 單頻連續波雷達僅可用于測速,無法測距;
- 多頻連續波雷達能測距,并且能夠分辨出固定目標和活動目標;
- 調頻連續波雷達即可測距又可測速,但只適用于單個目標。
二、脈沖體制
脈沖多普勒雷達的工作原理可表述如下:當雷達發射一固定頻率的脈沖波對空掃描時,如遇到活動目標,回波的頻率與發射波的頻率出現頻率差,稱為多普勒頻率。根據多普勒頻率的大小,可測出目標對雷達的徑向相對運動速度;根據發射脈沖和接收的時間差,可以測出目標的距離。同時用頻率過濾方法檢測目標的多普勒頻率譜線,濾除干擾雜波的譜線,可使雷達從強雜波中分辨出目標信號。所以脈沖多普勒雷達比普通雷達的抗雜波干擾能力強,能探測出隱蔽在背景中的活動目標。
三、連續波體制
1. 單頻連續波
單頻連續波雷達可設發射信號st(t)=Acos(2πf0t+φ)s_t(t)=Acos(2πf_0t+\varphi)st?(t)=Acos(2πf0?t+φ)則接收信號為sr(t)=kAcos[2πf0(t?τ)+φ]=kAcos{2πf0[t?2(R0?vrt)c]+φ}=kAcos(2πf0t?4πR0λ+2πfdt+φ)=kAcos(2πf0t+2πfdt+φ?φ0)\begin{aligned}s_r(t)&=kAcos[2πf_0(t-\tau)+\varphi]\\ &=kAcos\{2πf_0[t-\frac{2(R_0-v_rt)}{c}]+\varphi\}\\ &=kAcos(2\pi f_0t-\frac{4\pi R_0}{\lambda}+2\pi f_dt+\varphi)\\ &=kAcos(2\pi f_0t+2\pi f_dt+\varphi-\varphi_0) \end{aligned}sr?(t)?=kAcos[2πf0?(t?τ)+φ]=kAcos{2πf0?[t?c2(R0??vr?t)?]+φ}=kAcos(2πf0?t?λ4πR0??+2πfd?t+φ)=kAcos(2πf0?t+2πfd?t+φ?φ0?)?
當目標與雷達之間存在相對速度時,接受到的回波信號的載頻相對于發射信號的載頻產生一個頻移,即多普勒頻移fd=2vrλf_d=\frac{2v_r}{\lambda}fd?=λ2vr??
根據多普勒頻移,便可求出目標速度,但無法完成測距。
多普勒信息的提取:發射信號與接收信號過相位檢波器(混頻后取低頻),輸出形式為urcos(2πfdt?φ0),測頻即可得到fdu_rcos(2\pi f_dt-\varphi_0),測頻即可得到f_dur?cos(2πfd?t?φ0?),測頻即可得到fd?
2. 調頻連續波
FMCWFMCWFMCW其基本原理為發射波是高頻連續波,其頻率隨時間按照三角波(鋸齒波可看作是三角波的特殊形式)規律變化,接收的回波頻率與發射的頻率變化規律相同,都是三角波規律,只是有一個時間差,利用這個微小時間差可計算出目標距離。
FMCWFMCWFMCW雷達的簡易框圖:
頻率合成器產生一個線性調頻信號(LFMLFMLFM),線性調頻信號通過發射天線輻射出去,遇到目標之后發生反射,反射的信號被接收天線接收,將發射信號與接收信號進行混頻,得到中頻信號(IFsignalIF signalIFsignal),中頻信號的頻率為發射信號與接收信號的頻率差。
2.1 鋸齒波調頻連續波
雷達發射的線性調頻信號(忽略幅度參數和噪聲)為:
st(t)=ej(2πf0t+πKt2)s_t(t)=e^{j(2πf_0t+πKt^2)}st?(t)=ej(2πf0?t+πKt2)
接收的回波信號為:
sr(t)=ej[2πf0(t?τ)+πK(t?τ)2]s_r(t)=e^{j[2πf_0(t-\tau)+πK(t-\tau)^2]}sr?(t)=ej[2πf0?(t?τ)+πK(t?τ)2]
2.1.1 對于靜止目標:τ=2R0c\tau=\frac{2R_0}{c}τ=c2R0??
得到的中頻信號為:
s0(t)=sr(t)st?(t)=ej[2πf0(t?τ)+πK(t?τ)2]e?j(2πf0t+πKt2)=e?j(2πf0τ+2πKtτ?πKτ2)=e?j2π(2f0R0c+2KR0ct?2KR02c2)\begin{aligned} s_0(t)& =s_r(t)s_t^*(t) \\ & =e^{j[2πf_0(t-\tau)+πK(t-\tau)^2]}e^{-j(2πf_0t+πKt^2)} \\ & =e^{-j(2πf_0\tau+2πKt\tau-πK\tau^2)} \\ &=e^{-j2π(\frac{2f_0R_0}{c}+\frac{2KR_0}{c}t-\frac{2KR_0^2}{c^2})} \end{aligned}s0?(t)?=sr?(t)st??(t)=ej[2πf0?(t?τ)+πK(t?τ)2]e?j(2πf0?t+πKt2)=e?j(2πf0?τ+2πKtτ?πKτ2)=e?j2π(c2f0?R0??+c2KR0??t?c22KR02??)?
中頻信號頻率fIF=12πdφ(t)dt=2KR0cf_{IF}=\frac{1}{2\pi}\frac{\mathrmozvdkddzhkzd\varphi(t)}{\mathrmozvdkddzhkzdt}=\frac{2KR_0}{c}fIF?=2π1?dtdφ(t)?=c2KR0??
因此我們可以對中頻信號進行傅里葉變換,找出幅度譜峰值位置,即可得出其頻率fIFf_{IF}fIF?, 從而得出目標的距離為R0=cfIF2KR_0=\frac{cf_{IF}}{2K}R0?=2KcfIF??
雷達的距離分辨率:雷達在屏幕上所能分辨的兩個目標物的最小實際距離。
??根據前面的分析,雷達對目標距離的檢測最后是通過FFT運算,轉變為對中頻信號譜峰的定位,那么可以想到,FFT運算的分辨率將會影響中頻信號頻譜的分辨率,從而影響雷達對目標距離的分辨率。
由數字信號處理知識可得,若信號觀察時間為TTT,則FFT運算之后頻域的分辨率為1T\frac{1}{T}T1?,設兩目標相隔 ΔR\Delta RΔR,那么此時兩目標對應的中頻信號頻率差為
ΔfIF=KΔτ=K2ΔRc\Delta f_{IF}=K\Delta \tau=K\frac{2\Delta R}{c}ΔfIF?=KΔτ=Kc2ΔR?
注意到此時的中頻信號來源于對單Chirp信號的采樣,觀察時間為TpT_pTp? , 從而有
ΔfIF=K2ΔRc>1Tp\Delta f_{IF}=K\frac{2\Delta R}{c}>\frac{1}{T_p}ΔfIF?=Kc2ΔR?>Tp?1?
代入K=BTpK=\frac{B}{T_p}K=Tp?B?,得
ΔR>c2B\Delta R>\frac{c}{2B}ΔR>2Bc?
雷達分辨率Rres=c2BR_{res}=\frac{c}{2B}Rres?=2Bc?
可見,雷達對目標距離的分辨率由發射信號的帶寬決定,增加帶寬,將會得到更好的距離分辨能力,但同時也會增加硬件成本及信號處理的難度。
2.1.2 對于運動目標:τ(t)=2(R0+vrt)c\tau(t)=\frac{2(R_0+v_rt)}{c}τ(t)=c2(R0?+vr?t)? (以徑向速度vrv_rvr?遠離雷達)
得到的中頻信號為:
s0(t)=sr(t)st?(t)=ej[2πf0[t?τ(t)]+πK[t?τ(t)]2]e?j(2πf0t+πKt2)=e?j[2πf0τ(t)+2πKtτ(t)?πKτ(t)2]=e?j2π(2f0R0c+2f0vrct+2KR0ct+2Kvrct2?2KR02c2?4KR0vrc2t?2Kvr2c2t2)=e?j[[2π(2f0vrc+2KR0c?4KR0vrc2)t+(2Kvrc?2Kvr2c2)t2]+(4πf0R0c?4πKR02c2)]=ej[2π(fbt+μbt2)+φ]\begin{aligned} s_0(t)& =s_r(t)s_t^*(t) \\ & =e^{j[2πf_0[t-\tau(t)]+πK[t-\tau(t)]^2]}e^{-j(2πf_0t+πKt^2)} \\ & =e^{-j[2πf_0\tau(t)+2πKt\tau(t)-πK\tau(t)^2]} \\ &=e^{-j2π(\frac{2f_0R_0}{c}+\frac{2f_0v_r}{c}t+\frac{2KR_0}{c}t+\frac{2Kv_r}{c}t^2-\frac{2KR_0^2}{c^2}-\frac{4KR_0v_r}{c^2}t-\frac{2Kv_r^2}{c^2}t^2)} \\ &=e^{-j[[2π(\frac{2f_0v_r}{c}+\frac{2KR_0}{c}-\frac{4KR_0v_r}{c^2})t+(\frac{2Kv_r}{c}-\frac{2Kv_r^2}{c^2})t^2]+(\frac{4\pi f_0R_0}{c}-\frac{4\pi KR_0^2}{c^2})]}\\ &=e^{j[2\pi(f_bt+\mu_b t^2)+\varphi]} \end{aligned}s0?(t)?=sr?(t)st??(t)=ej[2πf0?[t?τ(t)]+πK[t?τ(t)]2]e?j(2πf0?t+πKt2)=e?j[2πf0?τ(t)+2πKtτ(t)?πKτ(t)2]=e?j2π(c2f0?R0??+c2f0?vr??t+c2KR0??t+c2Kvr??t2?c22KR02???c24KR0?vr??t?c22Kvr2??t2)=e?j[[2π(c2f0?vr??+c2KR0???c24KR0?vr??)t+(c2Kvr???c22Kvr2??)t2]+(c4πf0?R0???c24πKR02??)]=ej[2π(fb?t+μb?t2)+φ]?
其中fb=?2f0vrc?2KR0c+4KR0vrc2μb=?2Kvrc+2Kvr2c2φ=?4πf0R0c+4πKR02c2\begin{aligned} 其中f_b&=-\frac{2f_0v_r}{c}-\frac{2KR_0}{c}+\frac{4KR_0v_r}{c^2}\\ \mu_b&=-\frac{2Kv_r}{c}+\frac{2Kv_r^2}{c^2}\\ \varphi&=-\frac{4\pi f_0R_0}{c}+\frac{4\pi KR_0^2}{c^2} \end{aligned}其中fb?μb?φ?=?c2f0?vr???c2KR0??+c24KR0?vr??=?c2Kvr??+c22Kvr2??=?c4πf0?R0??+c24πKR02???
中頻信號瞬時頻率fIF(t)=12πdφ(t)dt=fb+2μbtf_{IF}(t)=\frac{1}{2\pi}\frac{\mathrmozvdkddzhkzd\varphi(t)}{\mathrmozvdkddzhkzdt}=f_b+2\mu_btfIF?(t)=2π1?dtdφ(t)?=fb?+2μb?t
從嚴格的數學意義上來說,fIF(t)f_{IF}(t)fIF?(t)依然是一個線性調頻信號,
但是由于處理時間極短(通常在ms或者us量級),可忽略ttt的高次項t2t^2t2;同時忽略掉分母含c2c^2c2的項,上式可簡化為
fIF=?2vrf0c?2KR0c\begin{aligned} f_{IF}&=-\frac{2v_rf_0}{c}-\frac{2KR_0}{c}\\ \end{aligned}fIF??=?c2vr?f0???c2KR0???
這依然是一個單頻信號,但與靜止目標所不同的是,此單頻信號的頻率中既包含距離信息,也包含速度信息,通過FFT無法直接準確測量出目標的距離信息,這就是距離速度耦合現象。當目標的運動速度很大時,耦合現象就會更嚴重,對測距造成的誤差就會越大。通常采用的方法就是速度補償,先將目標的徑向速度vrv_{r}vr?求出來,再去進一步修正距離值。還有一種方法及時采用Chirp Sequence的發射波形,這種波形的調頻周期很短,可以近似忽略物體在這段時間的運動速度對測距造成的影響。
測距的本質:測距的本質就是計算延時,而延時與中頻信號頻率fIFf_{IF}fIF?成線性關系,所以通過計算中頻信號頻率fIFf_{IF}fIF?就可以計算出回波延時,進而計算出目標的距離。
如何測速
目標的徑向速度可通過在一段時間內的距離變化率獲得,也可以通過測量多普勒頻移獲得。
2.2 三角波調頻連續波
測速測距原理:
藍色為發射信號頻率,黃色為接收信號頻率,掃頻周期為TTT,掃頻帶寬為BBB,發射信號經過目標反射,回波信號會有延時,在三角形的頻率變化中,可以在上升沿和下降沿兩者上進行距離測量。
發射信號:
st(t)=Acos(2πf0t+πμt2)s_t(t)=Acos(2\pi f_0t+\pi \mu t^2)st?(t)=Acos(2πf0?t+πμt2)
接收信號:sr(t)=kAcos{2πf0[t?2(R0?vrt)c]+πμ[t?2(R0?vrt)c]2(運動目標以徑向速度vr靠近雷達)\begin{aligned}s_r(t)&=kAcos\{2\pi f_0[t-\frac{2(R_0-v_rt)}{c}]+\pi \mu[t-\frac{2(R_0-v_rt)}{c}]^2 (運動目標以徑向速度v_r靠近雷達) \end{aligned}sr?(t)?=kAcos{2πf0?[t?c2(R0??vr?t)?]+πμ[t?c2(R0??vr?t)?]2(運動目標以徑向速度vr?靠近雷達)?
發射信號頻率:
ft(t)=f0+μtf_t(t)=f_0+\mu tft?(t)=f0?+μt
接收信號頻率:
fr(t)=f0+2vrλ+μ[t?2(R0?vrt)c](1+2vrc)(忽略vrc項)=f0+fd+μ(t?2R0c)\begin{aligned}f_r(t)&=f_0+\frac{2v_r}{\lambda}+\mu[t-\frac{2(R_0-v_rt)}{c}](1+\frac{2v_r}{c})(忽略\frac{v_r}{c}項)\\ &=f_0+f_d+\mu(t-\frac{2R_0}{c}) \end{aligned}fr?(t)?=f0?+λ2vr??+μ[t?c2(R0??vr?t)?](1+c2vr??)(忽略cvr??項)=f0?+fd?+μ(t?c2R0??)?
發射信號頻率:
ft(t)=f0?μtf_t(t)=f_0-\mu tft?(t)=f0??μt
接收信號頻率:
fr(t)=f0+fd?μ(t?2R0c)f_r(t)=f_0+f_d-\mu(t-\frac{2R_0}{c})fr?(t)=f0?+fd??μ(t?c2R0??)
逆程:fbd=fr?ft=2μR0c+fdf_{bd}=f_r-f_t=\frac{2\mu R_0}{c}+f_dfbd?=fr??ft?=c2μR0??+fd?
頻率計讀數為平均值,即fb=fbu+fbd2=2μR0cf_b=\frac{f_{bu}+f_{bd}}{2}=\frac{2\mu R_0}{c}fb?=2fbu?+fbd??=c2μR0??
對于靜止目標,沒有多普勒頻移,即fbu=fbdf_{bu}=f_{bd}fbu?=fbd?
對于運動目標,中頻(差拍)信號頻率中也摻雜了fdf_dfd?,因此fbu≠fbdf_{bu}\neq f_{bd}fbu??=fbd?,且fbu=fb?fdf_{bu}=f_b-f_dfbu?=fb??fd?,fbd=fb+fdf_{bd}=f_b+f_dfbd?=fb?+fd?
其中fbf_bfb?為被探測物體靜止時的頻差,fdf_dfd?是被探測物體移動時的多普勒頻移.
那么,可以得到,與距離相關的差拍頻率fb=fbu+fbd2f_b=\frac{f_{bu}+f_{bd}}{2}fb?=2fbu?+fbd??,與速度相關的多普勒頻移fd=fbd?fbu2f_d=\frac{f_{bd}-f_{bu}}{2}fd?=2fbd??fbu??
從而得到待測目標的距離和速度:
R=cτ(t)2=cfb2μR=\frac{c\tau(t)}{2}=\frac{cf_b}{2\mu}R=2cτ(t)?=2μcfb??
vr=fdc2f0v_r=\frac{f_dc}{2f_0}vr?=2f0?fd?c?
(其中fbf_bfb?和fdf_dfd?是需要測量的頻差,μ\muμ是已知的調頻斜率μ=2BT\mu=\frac{2B}{T}μ=T2B?)
1. MATLAB繪制周期鋸齒波/三角波信號
%用sawtooth()函數生成周期鋸齒波/三角波信號 t=linspace(0,20,200); subplot(2,2,1); f_1=sawtooth(t); %周期為2π,width默認為1 plot(t,f_1); title('鋸齒波信號'); f_2=sawtooth(t,0.5); %周期為2π,width為0.5 subplot(2,2,2); plot(t,f_2); title('三角波信號');%用tripuls()函數生成單個三角形信號 fs=10e3; ts=-0.1:1/fs:0.1; subplot(2,2,3); plot(ts,tripuls(ts,0.05,1)); %斜率為1,寬度為0.05的三角波 title('鋸齒波脈沖'); subplot(2,2,4); plot(ts,tripuls(ts,0.05)); %寬度為0.05的三角波 title('三角形脈沖');總結
以上是生活随笔為你收集整理的雷达原理---脉冲雷达和连续波雷达的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: onchange事件的触发
- 下一篇: 用打怪游戏实例,深度理解python面向