匹配滤波器设计仿真

1、输入信号能量: ：2 E(S)= s (t)dt ：(1.3)雷达系统匹配滤波器的仿真匹配滤波器原理在输入为确知加白噪声的情况下，所得输出信噪比最大的线性滤波器就是匹配滤波器，设一线性滤波器的输入信号为x(t):(1.1)x(t) =s(t) n(t)其中：s(t)为确知信号，n为均值为零的平稳白噪声，其功率谱密度为No/2。设线性滤波器系统的冲击响应为h(t)，其频率响应为H()，其输出响应:(1.2)y(t) =s°(t)n°(t)精选资料，欢迎下载输入、输出信号频谱函数:OOSC O= s(t) tdt-SaSOCO = HC)SC)(1.4)1 ： - j tSo(

2、t)H( JS( Je d 2兀输出噪声的平均功率：(1.5)21 严Eno(t)Pno(-)d 二2 兀-CO1 閃小ALOH )S)eaod<S"" - 1 ©H®)2Pn®)d®(1.6)利用SChWarZ不等式得:SNRo乞丄 P1C-)(1.7)上式取等号时，滤波器输出功率信噪比 SNR)最大取等号条件:H(J =:S C )Pn( Je jto(1.8)当滤波器输入功率谱密度是Pn(J=No2的白噪声时，MF的系统函数为:HC kSC )j to,k 土No(1.9)k为常数1，S*(')为输入函数频谱的复共

3、轭,S*( J = S(-'),也是滤波器的传输函数H CJSNR)2 EsNo(1.10)ES为输入信号s(t)的能量，白噪声n的功率谱为No/2SNR)只输入信号s(t)的能量ES和白噪声功率谱密度有关。白噪声条件下，匹配滤波器的脉冲响应：h(t) =ks*(t°-t)(1.11)如果输入信号为实函数，则与s(t)匹配的匹配滤波器的脉冲响应为:h(t) =ks(t° -t)(1.12)k为滤波器的相对放大量，一般k=1。匹配滤波器的输出信号：s°(t) =s°(t)* h(t) = kR(t 1°)(1.13)匹配滤波器的输出波形是

4、输入信号的自相关函数的k倍，因此匹配滤波器可以看成是一个计算输入信号自相关函数的相关器，通常k =1线性调频信号(LFM脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。 这种体制采用宽脉 冲发射以提高发射的平均功率，保证足够大的作用距离；而接受时采用相应的脉 冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分 辨率之间的矛盾。脉冲压缩雷达最常见的调制信号是线性调频(Lin ear Freque ncyModulation )信号,接收时采用匹配滤波器(MatChed FiIter )压缩脉冲LFM信号(也称ChirP信号)的数学表达式为:s(t) = rect(*)eK

5、2 j2二(fct -t )(2.1)frepencL 式中fc为载波频率，rect(*)为矩形信号,(2.2)K =，是调频斜率，于是，信号的瞬时频率为fc Kt (-T2 _t - T2) ,如图1(2.3)(2.4)图 1 典型的 ChirP 信号(a) UP-ChirP(K>0)( b) down-chirp(K<0)将2.1式中的UP-ChirP信号重写为:s(t) = S(t)ej2 二fct式中,S(t) = rect (*)ej Kt是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性， 只是中心频率不同而以，因此， MatIab仿真时，只

6、需考虑S(t)。通过MATLAB 仿真可得到信号时域和频域波形如下图所示：Real Part Of ChirP Signal图2.LFM信号的时域波形和幅频特性三线性调频信号的匹配滤波器 信号S(t)的匹配滤波器的时域脉冲响应为：h(t) =s* (to t)(3.1 )to是使滤波器物理可实现所附加的时延。理论分析时，可令to = 0,重写3.1式,h(t) =s*(-t)(3.2 )将2.1式代入3.2式得:(3.3 )匹配滤波h(t)图3.LFM信号的匹配滤波如图3, s(t)经过系统h(t)得输出信号So(t),sb(t)=s(t)* h(t)Oano=s(u)h(t-u)du = h

7、(u)s(t-u)du-OO-OO2rect(U)ej2K ej")2rect(u)ej2呱T2s°(t) = ey-2_j2 二 KtUJ "Kt2 e=e 2j*te j 二KtUdU一 -j2JiKt t-T2Sin 二 K(T - t)t j2-fct eejf二 Kt(3.4)t T 2的=ej“：Kt2e»KtUdU合并3.4和3.5两式:s3(t) =T一j2Ktu#Kt2 ej2兀 Kt 一2Sin : K(T t)t j2-fcte(3.5): KttSZKTF二 KTtrect(2F)严(3.6)3.6式即为LFM脉冲信号经匹配滤波器

8、得输出,它是一固定载频fc的信号 t乞T时，包络近似为辛克(SinC )函数。S0(t) =TSa(二 KTt)rect() =TSa(二 Bt)rect()(3.7 )图4.匹配滤波的输出信号如图4,习惯上,1JF1当二Bt=二时，t =为其第一零点坐标；当二Bt 时，t -B22B将此时的脉冲宽度定义为压缩脉冲宽度。鎳12B B(3.8)LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度之比通常称为压缩比D,(3.9)D=- =TB3.9式表明，压缩比也就是LFM信号的时宽频宽积。由(2.1),(3.3), (3.6 )式，s(t),h(t),so(t)均为复信号形式，Matab仿真时，只需考

9、虑它们的复包络 S(t),H(t),So(t)即可。经MATLA筋真得线性调频信号经过匹配滤波器的 波形信号如图5所示:8pnluv-13 4ChlrP SIgnaI after matched filtere ImulationalSinCO -2-4 ¾>pn<Ev15-IO-5O510Tjnrle in SeC XSChirP Signai after matched filter (ZoOm)-S-3-2-1-0.5 O 0 5123Time in SeC XS图5.Chirp信号的匹配滤波图5中，时间轴进行了归一化，（t心B）=txB）°图中反映出理论

10、与仿真结1果吻合良好。第一零点出现在_1 （即一）处，此时相对幅度-13.4dB。压缩后B11的脉冲宽度近似为 一（± ），此时相对幅度-4dB,这理论分析（图3.2 ） 一致。B2B如果输入脉冲幅度为1,且匹配滤波器在通带内传输系数为 1,则输出脉冲幅度 为.kT2 -.TB 一 D ,即输出脉冲峰值功率PO比输入脉冲峰值功率P增大了 D 倍。四雷达系统对线性调频信号的检测在实际实际雷达系统中，LFM脉冲的处理过程如图6图6 LFM信号的接收处理过程雷达回波信号Sr（t）经过正交解调后，得到基带信号，再经过匹配滤波脉冲压缩后就可以作出判决。正交解调原理如图7,雷达回波信号经正交解调

11、后得两路相互正交的信号（t）和Q（t）。一种数字方法处理的的匹配滤波原理如图&, PrfflJtt)SmC27=PrlOIt)图7正交解调原理图8 一种脉冲压缩雷达的数字处理方式以下各图为经过脉冲压缩输出的已加噪声的线性调频信号（模拟雷达回波信号）的matlab仿真结果：波形参数脉冲宽度T=10s ,载频频率fc=10khz,脉冲宽度 B=30MhzRange in meters1j图9.SNR=30的脉冲压缩输入输出波形图10 SNR=20的脉冲压缩输入输出波形图11 SNR=O的脉冲压缩输入输出波形图12 SNR=-Io的脉冲压缩输入输出波形Range in metersX-IQ*

12、图13. SNR=-20的脉冲压缩输入输出波形图14. SNR=-30的脉冲压缩输入输出波形S(t)信号中白噪声n为：n = Sqrt(0.5* SNR) * (randn (1, length(St) j * randn( 1,length (St) >仿真表明，线性调频信号经匹配滤波器后脉冲宽度被大大压缩，信噪比得到了显著提高，但是雷达目标回波信号信号的匹配滤波仿真结果图9-14可以看出当信噪比小于零时随着信噪比的不断减小，所噪声对线性调频信号的干扰愈来愈 明显，当信噪比达到-30dB时已经有部分回波信号被淹没了，也就是说当信噪比 更小时即使是经过脉冲压缩，噪声仍能淹没有用信号。五程

13、序附录1.线性频率调制信号(LFM仿真:%demo Of ChirP Sig nalT=10e-6;B=30e6;K=B/T;Fs=2*B;TS=1/Fs;%pulse durati On 10us%chirp freque ncy modulati On ban dwidth 30MHz%chirp slope%sampli ng freque ncy and SamPIe SPaCi ngN=TTs;t=li nspace(-T2,T2,N);St=exp(1i*pi*K*t.2);%ge nerate ChirP Sig nalsubplot(211)plot(t*1e6,real(St

14、);xlabel('Time in U sec');title('Real Part of ChirP Sig nal');grid on; axis tight;subplot(212)freq=li nspace(-Fs2,Fs2,N);plot(freq*1e-6,fftshift(abs(fft(St);xlabel('Freque ncy in MHz');title('Mag nitude SPeCtrUm of ChirP Sig nal');grid on; axis tight;2 LFM信号的匹配滤波仿真%de

15、mo of ChirP Sig nal after matched filterT=10e-6;B=30e6;30MHz%pulse durati On 10us%chirpfreque ncymodulati Onban dwidthK=B/T;%chirp slopeFs=10*B;Ts=1/Fs;N=T/Ts;%sampli ng freque ncy and SamPIe SPaCi ngt=li nspace(-T2,T2,N); St=exp(j*pi*K*t.2); Ht=exp(-j*pi*K*t.2); Sot=co nv(St,Ht);subplot(211) L=2*N-

16、1;%chirp Sig nal%matched filter%chirp Sig nal after matched filtert1=li nspace(-T,T,L); Z=abs(Sot)Z=Zmax(Z); Z=20*log10(Z+1e-6); Z1=abs(si nc(B.*t1); Z1=20*log10(Z1+1e-6); t1=t1*B;%no rmalize%sinc fun Cti Onplot(t1,Z,t1,Z1,'r.'); axis(-15,15,-50,i nf);grid on; lege nd('emulati On al'

17、,'s in c'); xlabel('Time in SeC timesitB'); ylabel('Amplitude,dB');title('Chirp Sig nal after matched filter'); subplot(212)%zoomN0=3*FsB;t2=-N0*Ts:Ts:N0*Ts;t2=B*t2;plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');axis(-i nf,in f,-50,i nf);grid on;set(gca,'Ytic

18、k',-134,-4,0,'Xtick',-3,-2,-1,-0.5,0,05,1,2,3);XIabeI('Time in SeC timesitB');ylabel('Amplitude,dB');titIe('Chirp Sig nal after matched filter (Zoom)');3. LFM信号的雷达监测仿真% in put(' nPulse radar COmPreSSi On PrOCeSS ing: n ');clear;close all;T=10e-6;B=30e6;Rmi

19、 n=8500;RmaX=11500;R=9000,10000,10200;RCS=1 1 1 ;C=3e8;K=BT;RWid=RmaX-Rmi n;TWid=2*RwidC;Fs=5*B;TS=1/Fs;NWid=Ceil(TwidTs);t=li nspace(2*Rmi nC,2*RmaxC,Nwid);M=le ngth(R);td=o nes(M,1)*t-2*R7C* On es(1,Nwid);SNR=1,0.1,001,0001,10,100,1000;for i=1:1:7SrtI=RCS*(exp(1i*pi*K*td.2).*(abs(td)<T2);n=sqrt(0.5*SNR(i)*(ra ndn (SiZe(Srt1)+1i*ra ndn (SiZe(Srt1);Srt=Srt1+ n;%Digt