应用于雷达系统匹配滤波器地matlab仿真

1、实用标准文案输入信号能量:E(s) =s2(t)dt :二(1.3)应用于雷达系统匹配滤波器的 matlab仿真.匹配滤波器原理在输入为确知加白噪声的情况下，所得输出信噪比最大的线性滤波器就是匹配滤波器，设一线性滤波器的输入信号为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)精彩文档输入、输出信号频谱函数:qQS( ) = j二s(t)etdtSo( -) = H( .)S

2、( .)(1.4)1 j tso(t)H( )S( )e d 2 二一输出噪声的平均功率：21: Eh二萩.').二 H2( )Pn( )d，(1.5)1 二j t,.H( -)S( )ej tod2 二-2.'-H( ) Pn( )d()SNRo =亍(1.6)利用Schwarz不等式得:1 SNRo - 一2 二产|s2d“aP1( )(1.7)上式取等号时，滤波器输出功率信噪比 SNR最大取等号条件:*.(1.8)H(.)=iSW 卜Pn()当滤波器输入功率谱密度是Pn(&) = No/2的白噪声时，MF的系统函数为:H( ) =kS*( )e,to, k =&#

3、163;(1.9)Nok为常数1, S (。)为输入函数频谱的复共腕，S ) = S(f),也是滤波器的传输函数H(©)02EvSNR =2Es(1.10)NoEs为输入信号s(t)的能量，白噪声n的功率谱为No/2SNR只输入信号s(t)的能量Es和白噪声功率谱密度有关。白噪声条件下，匹配滤波器的脉冲响应：* h(t) =ks (to-t)(1.11)如果输入信号为实函数，则与s(t)匹配的匹配滤波器的脉冲响应为：h(t) =ks(t。-t)(1.12)k为滤波器的相对放大量，一般k=1。匹配滤波器的输出信号：so(t)=s0(t)*h(t) -kR(t-to)(1.13)匹配滤波

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

5、t2)s(t) = rectq)e2式中fc为载波频率，rect(，)为矩形信号，ttrectGgT0 , elsewise(2.1)(2.2)在eqiiency 1K=B，是调频斜率，于是，信号的瞬时频率为 fc+Kt (-T2<t<T2),如图1图 1 典型的 chirp 信号(a) up-chirp(K>0) ( b) down-chirp(K<0)将2.1式中的up-chirp信号重写为：s(t) = S(t)ej26ct(2.3)式中，S(t) = rect (T；)ej71Kt(2.4 )是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的

6、幅频特性， 只是中心频率不同而以，因此， Matlab仿真时，只需考虑S(t)。通过MATLAB 仿真可得到信号时域和频域波形如下图所示：图2.LFM信号的时域波形和幅频特性三.线性调频信号的匹配滤波器 信号s(t)的匹配滤波器的时域脉冲响应为:一 一一 * 一h(t) =s (to -t)(3.1 )to是使滤波器物理可实现所附加的时延。理论分析时，可令 to = 0,重写3.1式,*h(t) =s )(3.2)将2.1式代入3.2式得：h(t) =rect(：t；)ejKt2 ej"fct(3.3 )S(t)*匹配滤波M)so班 ,图3.LFM信号的匹配滤波如图3, s(t)经过

7、系统h(t)得输出信号So(t),&(t)=s(t)*h(t)=j_s(u)h(t -u)du = h(u)s(t - u)duDO=ejKu-oOjoO2 Uj2-fcuj-K(t _u)2t -'U j2-fc (t _u)rect(：j：)ejerect(-)ej)du当0 土 ET时,T2s(o(t) = et J22j-Kte2-Ktudu2Cj21：KtuJ -Kt2 e二 e T2-j2JiKt t-T2sin 二 K(T-t)t j2- fct e二 Kt(3.4)当T Et E0时,t T2%(t)=ej - Kt2e-du工2一j2Ktu2 e-j2nKt

8、一%j2二 fcte(3.5)合并3.4和3.5两式:与(t) =Tsin 二 K(T t)t j2-%t e二 Ktt sin n KT (1-)t二 KTtt j2：fj rect( )e2T(3.6)3.6式即为LFMB冲信号经匹配滤波器得输出，它是一固定载频fc的信号。当 tWT时，包络近似为辛克(sinc )函数。So(t)=TSa(二 KTt)rect(')=TSa(二 Bt)rect(5)(3.7)4lSo(t)|图4.匹配滤波的输出信号1 1如图4,习惯上,当nBt=±n时，t = 士为其弟一苓点坐标；当nBt = 士时 t=±一B22B将此时的脉冲

9、宽度定义为压缩脉冲宽度。2B(3.8)LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度工之比通常称为压缩比D,D=T=TB(3.9)3.9式表明，压缩比也就是LFM信号的时宽频宽积。由(2.1 ) , (3.3 ) , (3.6)式，s(t),h(t),so(t)均为复信号形式，Matab仿真时，只需考 虑它们的复包络S(t),H(t),So(t) 即可。经MATLA昉真得线性调频信号经过匹配滤波器的波形信号如图5所示：华三一-dEWChirp signal after matched filter15-5O510Time in sec xBChirp signal after matched

10、fiber (Zoom) 2-4 mp/pn 莒 dEu-13 4I 51-0.5 Q 0 51Time in sec mE23图5.Chirp信号的匹配滤波图5中，时间轴进行了归一化，(t/(1/B)=txB)。图中反映出理论与仿真结 1果吻合良好。第一零点出现在 ±1即士一)处，此时相对幅度-13.4dB。压缩后 B11 一.的脉冲宽度近似为 一 士一)，此时相对幅度-4dB,这理论分析(图3.2) 一致。 B2B如果输入脉冲幅度为1,且匹配滤波器在通带内传输系数为 1,则输出脉冲幅度 为jkRn'斤B=JD,即输出脉冲峰值功率P0比输入脉冲峰值功率P增大了 D 倍。四.

11、雷达系统对线性调频信号的检测在实际实际雷达系统中，LFM永冲的处理过程如图6图6 LFM信号的接收处理过程雷达回波信号Sr(t)经过正交解调后，得到基带信号，再经过匹配滤波脉冲压缩后就可以作出判决。正交解调原理如图7,雷达回波信号经正交解调后得两路相互正交的信号I(t)和Q(t)。一种数字方法处理的的匹配滤波原理如图8。实用标准文案HnQXpFffrt)图7 正交解调原理图8 一种脉冲压缩雷达的数字处理方式以下各图为经过脉冲压缩输出的已加噪声的线性调频信号（模拟雷达回波信号）的matlab仿真结果：波形参数脉冲宽度T=10Ns,载频频率fc=10khz,脉冲宽 c度 B=30Mhz加噪线性调频

12、信号压缩前.SNR =3。IngRange in meters图9.SNR=30的脉冲压缩输入输出波形精彩文档实用标准文案精彩文档加噪线性调频信号压缩前,SNR =2058606264666670727476图10 SNR=20的脉冲压缩输入输出波形图11 SNR=0的脉冲压缩输入输出波形实用标准文案图13. SNR=-20的脉冲压缩输入输出波形精彩文档图12 SNR=-10的脉冲压缩输入输出波形实用标准文案精彩文档加噪线性调频信号压缩前.SNR =3058606264666870727476图14. SNR=-30的脉冲压缩输入输出波形S(t)信号中白噪声n为:n =sqrt(0.5* SN

13、R)* (randn (1, length(St) + j * randn(1,length (St)、仿真表明，线性调频信号经匹配滤波器后脉冲宽度被大大压缩，信噪比得到 了显著提高，但是雷达目标回波信号信号的匹配滤波仿真结果图9-14可以看出当信噪比小于零时随着信噪比的不断减小，所噪声对线性调频信号的干扰愈来愈 明显，当信噪比达到-30dB时已经有部分回波信号被淹没了，也就是说当信噪比 更小时即使是经过脉冲压缩，噪声仍能淹没有用信号。五.程序附录%demo of chirp signal T=10e-6;B=30e6;K=B/T;Fs=2*B;Ts=1/Fs;1 .线性频率调制信号(LFMD

14、仿真:%pulse duration10us%chirp frequency modulation bandwidth 30MHz%chirp slope%sampling frequency and sample spacingN=T/Ts;t=linspace(-T/2,T/2,N);St=exp(1i*pi*K*t，2);%generate chirp signalsubplot(211)plot(t*1e6,real(St);xlabel( 'Time in u sec' );title( 'Real part of chirp signal');gri

15、d on;axis tight ;subplot(212)freq=linspace(-Fs/2,Fs/2,N);plot(freq*1e-6,fftshift(abs(fft(St);xlabel( 'Frequency in MHz');title( 'Magnitude spectrum of chirp signal');grid on;axis tight ;2 LFM信号的匹配滤波仿真%demo of chirp signal after matched filterT=10e-6;B=30e6;30MHzK=B/T;%pulse duration1

16、0us %chirp%chirp slopefrequency modulation bandwidthFs=10*B;Ts=1/Fs;%sampling frequency and sample spacingN=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.A2);Ht=exp(-j*pi*K*t.A2);Sot=conv(St,Ht);subplot(211)L=2*N-1;t1=linspace(-T,T,L);Z=abs(Sot);Z=Z/max(Z);Z=20*log10(Z+1e-6);Z1=abs(sinc(B.*t1);%chirp

17、signal%matched filter%chirp signal after matched filter%normalize%sinc functionZ1=20*log10(Z1+1e-6);t1=t1*B;plot(t1,Z,t1,Z1,'r.');axis(-15,15,-50,inf);grid on;legend('emulational','sinc');xlabel('Time in sec timesitB');ylabel('Amplitude,dB');title('Chirp s

18、ignal after matched filter'); subplot(212)%zoomN0=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)；grid on；set(gca,'Ytick',-13.4,-4,0,'Xtick',-3,-2,-1,-0.5,0,0.5,123);xlabel('Time in sec timesitB');ylabel('Am

19、plitude,dB');title('Chirp signal after matched filter (Zoom)');3. LFM信号的雷达监测仿真% input('nPulse radar compression processing: n ');clear;close all ;T=10e-6;B=30e6;Rmin=8500;Rmax=11500;R=9000,10000,10200;RCS=1 1 1 ;C=3e8;K=B/T;Rwid=Rmax-Rmin;Twid=2*Rwid/C;Fs=5*B;Ts=1/Fs;Nwid=ceil(Twid/Ts);t=linspace(2*Rmin/C,2*Rmax/C,Nwid);M=length(R);td=ones(M,1)*t-2*R'/C*ones(1,Nwid);SNR=1,0.1,0.01,0.001,10,100,1000; for i=1:1:7Srt1=RCS*(exp(1i*pi*K*td.A2).*(abs(td)<T/2);n=sqrt(0.5*SNR(i)*(randn(size(Srt