雷达脉冲压缩

第一篇：雷达脉冲压缩

雷达脉冲压缩

1、雷达工作原理

雷达是Radar（Radio Detection And Ranging）的音译词，意为“无线电检测和测距”，即利用无线电波来检测目标并测定目标的位置，这也是雷达设备在最初阶段的功能[1]。典型的雷达系统如图1.1，它主要由发射机，天线，接收机，数据处理，定时控制，显示等设备组成。利用雷达可以获知目标的有无，目标斜距，目标角位置，目标相对速度等。现代高分辨雷达扩展了原始雷达概念，使它具有对运动目标(飞机，导弹等)和区域目标(地面等)成像和识别的能力。雷达的应用越来越广泛。

图1.1 简单脉冲雷达系统框图

雷达发射机的任务是产生符合要求的雷达波形（Radar Waveform），然后经馈线和收发开关由发射天线辐射出去，遇到目标后，电磁波一部分反射，经接收天线和收发开关由接收机接收，对雷达回波信号做适当的处理就可以获知目标的相关信息。

假设理想点目标与雷达的相对距离为R，为了探测这个目标，雷达发射信号s(t),电磁波以光速C向四周传播，经过时间RC后电磁波到达目标，照射到目标上的电磁波可写成：R)。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为CRs(t)，其中为目标的雷达散射截面（Radar Cross Section ,简称RCS），反映目标对CR电磁波的散射能力[2]。再经过时间RC后，被雷达接收天线接收的信号为s(t2)。

Cs(t如果将雷达天线和目标看作一个系统，便得到如图1.2的等效，而且这是一个LTI（线性时不变）系统。

图1.2 雷达等效于LTI系统

等效LTI系统的冲击响应可写成:

h(t)(t)

(1.1)

iii1MM表示目标的个数，i是目标散射特性，i是光速在雷达与目标之间往返一次的时间，i2Ri

(1.2)c式中，Ri为第i个目标与雷达的相对距离。

雷达发射信号s(t)经过该LTI系统，得输出信号(即雷达的回波信号)sr(t)：

sr(t)s(t)*h(t)s(t)*i(ti)is(ti)

(1.3)

i1i1MM

图1.3 雷达回波信号处理

s(t)的匹配滤波器hr(t)为：

hr(t)s*(t)

(1.4)于是，so(t)sr(t)*hr(t)s(t)*s*(t)*h(t)

（1.5）对上式进行傅立叶变换：

So(jw)S(jw)S*(jw)H(jw)|S(jw)|H(jw)(1.6)如果选取合适的s(t)，使它的幅频特性|S(jw)|为常数，那么1.6式可写为：

So(jw)kH(jw)

(1.7)其傅立叶反变换为：

so(t)kh(t)ki1Mi(ti)

(1.8)so(t)中包含目标的特征信息i和i。从 so(t)中可以得到目标的个数M和每个目标相对雷达的距离：

Riic

（1.9）2这也是线性调频（LFM）脉冲压缩雷达的工作原理。

2、线性调频（LFM）信号

脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率，保证足够大的作用距离；而接受时采用相应的脉冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分辨率之间的矛盾。

脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation）信号,接收时采用匹配滤波器（Matched Filter）压缩脉冲。

LFM信号(也称Chirp 信号)的数学表达式为：

tj2(fctKt2)2s(t)rect()e

(2.1)

T式中fc为载波频率，rect()为矩形信号，tTtt1

1（2.2）

rect()TT0,elsewiseKB，是调频斜率，于是，信号的瞬时频率为fcKt(TtT)，如图2.1

22T

图2.1 典型的chirp信号（a）up-chirp(K>0)（b）down-chirp(K<0)将2.1式中的up-chirp信号重写为：

s(t)S(t)e式中，j2fct

（2.3）

tjKt2S(t)rect()e

（2.4）

T是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需考虑S(t)[3]。

3、FM脉冲的匹配滤波

信号s(t)的匹配滤波器的时域脉冲响应为：

h(t)s(t0t)

（3.1）

*t0是使滤波器物理可实现所附加的时延。理论分析时，可令t0＝0，重写3.1式，h(t)s*(t)

（3.2）将2.1式代入3.2式得:

h(t)rect()etTjKt2ej2fct

(3.3)

图3.1：LFM信号的匹配滤波

如图3.1,s(t)经过系统h(t)得输出信号so(t)，so(t)s(t)*h(t)

s(u)h(tu)duh(u)s(tu)du

当0tT时，uj2fcutuj2fc(tu)jKu2jK(tu)2erect()eerect()eduTTT22s0(t)tT2ejKtej2Ktuduej2KtuT2ej2fct

(3.4)

j2KttT ejKt2当Tt0时，sinK(Tt)tj2fcteKt2tT2s0(t)T2j2KtutT2j2fct

(3.5)

ejKt2eeejKtej2Ktuduj2KtT2合并3.4和3.5两式： sinK(Tt)tj2fcteKttsinKT(1)tTrect(t)ej2fct

（3.6）

s0(t)TKTt2T3.6式即为LFM脉冲信号经匹配滤波器得输出,它是一固定载频fc的信号。当tT时，包络近似为辛克（sinc）函数。

S0(t)TSa(KTt)rect(tt)TSa(Bt)rect()

（3.7）2T2T

图3.2：匹配滤波的输出信号

如图3.2，当Bt时，t11为其第一零点坐标；当Bt时，t，习

2B2B惯上，将此时的脉冲宽度定义为压缩脉冲宽度。

112

(3.8)2BBTTB

（3.9）LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度之比通常称为压缩比D，D3.9式表明，压缩比也就是LFM信号的时宽频宽积。

4、matlab程序与结果分析

（1）主要程序：

x(:,1:n)= 0.;

y(1:n)= 0.;

replica(1:n)= 0.;

replica = exp(i * pi *(b/taup).* t.^2);

for

j = 1:1:nmb

range = smb_range(j);

x(j,:)= smb_rcs(j).*exp(-i*2*pi*f0*2*range/c).* exp(i * pi *(b/taup).*(t +(2*range/c)).^2);

y = x(j,:)+ y;

end

rfft = fft(replica,nfft);

yfft = fft(y,nfft);

out= abs(ifft((rfft.* conj(yfft))))./(nfft);

s = taup * c /2;Npoints = ceil(rrec * nfft /s);

dist =linspace(0, rrec, Npoints);

（2）运行结果

雷达信号对三个目标进行探测，本设计中，分辨出相邻两点目标的最小距离是1.5米。故当两点距离小于1.5米时，雷达将分辨不出，即视为一点。

5、结论

脉冲压缩雷达既保持了窄脉冲的高距离分辨力，又能获得宽脉冲的强检测能力。脉冲压缩技术是大时宽带宽乘积信号经过匹配滤波器实现的, 不同的信号形式有不同的压缩性能, 其中线性调频脉冲信号的诸多优点使其称为脉冲压缩信号的首选,它也是最早、应用最广泛的脉冲压缩信号。脉冲压缩技术能在雷达发射功率受限的情况下, 提高目标的探测距离, 并且保持很高的分辨力, 是雷达反隐身、多目标分辨、抗干扰的重要手段, 在目前的雷达信号系统中有着广泛的应用。

6、参考文献

[1]丁鹭飞，耿富录·雷达原理（第三版）·西安：西安电子科技大学出版社，2006.[2]楼顺天，姚若玉，沈俊霞·MATLAB程序设计语言·西安：西安电子科技大学西电出版社，2007.[3]元春，苏广州，米红·宽带雷达信号产生技术[M]·北京：国防工业出版社，2002.附录：

nmb = 4;

%目标个数3个

rrec =200;

%测量目标最远距离 b = 100e6;

%调频信号带宽

smb_range = [10,30,40,100];

%三点目标的距离

最小分辨距离为s=c/2b=1.5m smb_rcs = [1 1 1 2];

%三点目标的横截面积 taup = 0.005e-3;

%信号持续脉宽 f0 = 5.6e9;

% 载频频率

c = 3e8;

% 信号传播的速度，即光速 fs = 2*b;

% 采样的频率 sampling_interval = 1/fs;

n = fix(taup/sampling_interval);

%总共点数（取整）nfft =n

% 采样点数 freqlimit = 0.5*fs;

freq = linspace(-freqlimit,freqlimit,n);% 频率采样间隔 = fs/n = 1/taup;t = linspace(-taup/2,taup/2,n);

%相邻点时间间隔 x(:,1:n)= 0.;

% x为矩阵 y(1:n)= 0.;

replica(1:n)= 0.;replica = exp(i * pi *(b/taup).* t.^2);%基带线性调频信号

for j = 1:1:nmb

%矩阵方法将接收信号叠加

range = smb_range(j);

x(j,:)= smb_rcs(j).*exp(-i*2*pi*f0*2*range/c).* exp(i * pi *(b/taup).*(t +(2*range/c)).^2);%接收信号

y = x(j,:)+ y;

%信号叠加 end rfft = fft(replica,nfft);yfft = fft(y,nfft);out= abs(ifft((rfft.* conj(yfft))))./(nfft);s = taup * c /2;Npoints = ceil(rrec * nfft /s);dist =linspace(0, rrec, Npoints);

%图片显示： figure subplot(311)plot(t,real(replica));axis tight;xlabel('Range in meters');ylabel('Amplitude in dB');title('线性调频信号');subplot(312)plot(t,real(y));axis tight;xlabel('Range in meters');ylabel('Amplitude in dB');title('压缩前雷达回波');subplot(313)plot(dist, out(1:Npoints))xlabel('Target relative position in meters')ylabel('压缩后雷达回波')

grid

第二篇：线性调频(LFM)脉冲压缩雷达仿真

线性调频（LFM）脉冲压缩雷达仿真

班级：信息五班

姓名：李辉 04141394 李港深 04141393 李少杰 04141395

分工：李辉

报告

李少杰李港深

程序及调试

概述：雷达工作原理

雷达是Radar（RAdio Detection And Ranging）的音译词，意为“无线电检测和测距”，即利用无线电波来检测目标并测定目标的位置，这也是雷达设备在最初阶段的功能。它是通过发射电磁波并接收回波信号，在后端经过信号处理将目标的各种特性分析出来的一个复杂的系统。其中，雷达回波中的可用信息包括目标斜距,角位置，相对速度以及目标的尺寸形状等。典型的雷达系统如图1.1，它主要由发射机，天线，接收机，数据处理，定时控制，显示等设备组成。利用雷达可以获知目标的有无，目标斜距，目标角位置，目标相对速度等。现代高分辨雷达扩展了原始雷达概念，使它具有对运动目标(飞机，导弹等)和区域目标(地面等)成像和识别的能力。雷达的应用越来越广泛。

图1.1：简单脉冲雷达系统框图

一．线性调频（LFM）脉冲压缩雷达原理

假设理想点目标与雷达的相对距离为R，为了探测这个目标，雷达发射信号s(t),电磁波以光速C向四周传播，经过时间RC后电磁波到达目标，照射到目标上的电磁波可写成：R)。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为CRs(t)，其中为目标的雷达散射截面（Radar Cross Section ,简称RCS），反映目标对CR电磁波的散射能力。再经过时间RC后，被雷达接收天线接收的信号为s(t2)。

Cs(t如果将雷达天线和目标看作一个系统，便得到如图1.2的等效，而且这是一个LTI（线性时不变）系统。

图1.2：雷达等效于LTI系统

等效LTI系统的冲击响应可写成:

h(t)

(t)

(1.1)

iii1MM表示目标的个数，i是目标散射特性，i是光速在雷达与目标之间往返一次的时间，i2Ri

(1.2)c式中，Ri为第i个目标与雷达的相对距离。

雷达发射信号s(t)经过该LTI系统，得输出信号(即雷达的回波信号)sr(t)：

sr(t)s(t)*h(t)s(t)*i(ti)is(ti)

(1.3)

i1i1MM那么，怎样从雷达回波信号sr(t)提取出表征目标特性的i(表征相对距离)和i(表征目标反射特性)呢？常用的方法是让sr(t)通过雷达发射信号s(t)的匹配滤波器，如图1.3。

图1.3：雷达回波信号处理

s(t)的匹配滤波器hr(t)为：

hr(t)s*(t)

(1.4)于是，so(t)sr(t)*hr(t)s(t)*s*(t)*h(t)

（1.5）对上式进行傅立叶变换：

So(jw)S(jw)S*(jw)H(jw)|S(jw)|H(jw)2

(1.6)如果选取合适的s(t)，使它的幅频特性|S(jw)|为常数，那么1.6式可写为：

So(jw)kH(jw)

(1.7)

式中，tjKt2S(t)rect()e

（2.4）

T是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需考虑S(t)。以下Matlab程序产生2.4式的chirp信号，并作出其时域波形和幅频特性，如图2.2。>> %%线性调频信号的程序 T=10e-6;B=30e6;K=B/T;

Fs=2*B;Ts=1/Fs;N=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.^2);线性调频信号 subplot(211)plot(t*1e6,real(St));xlabel('Time in u sec');title('Real part of chirp signal');grid 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.2：LFM信号的时域波形和幅频特性三． LFM脉冲的匹配滤波信号s(t)的匹配滤波器的时域脉冲响应为：

h(t)s*(t0t)

t0是使滤波器物理可实现所附加的时延。理论分析时，可令t0＝0，重写3.1式，h(t)s*(t)

将2.1式代入3.2式得:

h(t)rect(tjKt2fctT)eej2

图3.1：LFM信号的匹配滤波

如图3.1,s(t)经过系统h(t)得输出信号so(t)，3.1）3.2）

(3.3)

（

（so(t)s(t)*h(t)

s(u)h(tu)duh(u)s(tu)du

当0tT时，uj2fcujK(tu)2tuj2fc(tu)jKu2erect()eerect()eduTTT2s0(t)tT2ejKt2j2Ktuedu

ejKt2ej2KtuT2ej2fct

(3.4)

j2KttT2当Tt0时，sinK(Tt)tj2fcteKttT2s0(t)T2ejKtej2Ktuduej2KtutT2j2fct

(3.5)eTj2Kt22

ejKt2合并3.4和3.5两式： sinK(Tt)tj2fcteKttsinKT(1)tTrect(t)ej2fct

（3.6）

s0(t)TKTt2T3.6式即为LFM脉冲信号经匹配滤波器得输出,它是一固定载频fc的信号。当tT时，包络近似为辛克（sinc）函数。

S0(t)TSa(KTt)rect(tt)TSa(Bt)rect()

（3.7）2T2T

xlabel('Time in sec timesitB');ylabel('Amplitude,dB');title('Chirp signal after matched filter(Zoom)');仿真结果如图3.3：

图3.3：Chirp信号的匹配滤波

图3.3中，时间轴进行了归一化，（t/(1/B)tB）。图中反映出理论与仿真结果吻合良好。第一零点出现在1（即（11）处，此时相对幅度-13.4dB。压缩后的脉冲宽度近似为BB1），此时相对幅度-4dB,这理论分析（图3.2）一致。2B

上面只是对各个信号复包络的仿真，实际雷达系统中，LFM脉冲的处理过程如图3.4。

图3.4： LFM信号的接收处理过程

雷达回波信号sr(t)（1.4式）经过正交解调后，得到基带信号，再经过匹配滤波脉冲压缩后就可以作出判决。正交解调原理如图3.5，雷达回波信号经正交解调后得两路相互正交的信号I(t)和Q(t)。一种数字方法处理的的匹配滤波原理如图3.6。

图3.5：正交解调原理

图3.6：一种脉冲压缩雷达的数字处理方式

四：Matlab仿真结果

（1）任务：对以下雷达系统仿真。

雷达发射信号参数：幅度：1.0 信号波形：线性调频信号

频带宽度：30兆赫兹（30MHz）脉冲宽度：10微妙（20us）中心频率：1GHz（109Hz）雷达接收方式：正交解调接收

距离门：10Km~15Km 目标：

Tar1：10.5Km Tar2：11Km Tar3：12Km Tar4：12Km＋10m Tar5：13Km Tar6：13Km＋25m（2）系统模型：

结合以上分析，用Matlab仿真雷达发射信号，回波信号，和压缩后的信号的复包络特性，其载频不予考虑（实际中需加调制和正交解调环节），仿真信号与系统模型如图4.1。

图4.1：雷达仿真等效信号与系统模型

01112-

第三篇：雷达线性调频脉冲压缩的原理及其MATLAB仿真

LFM脉冲压缩雷达仿真

线性调频（LFM）脉冲压缩雷达仿真

宋萌瑞 201421020302 一．雷达工作原理

雷达是Radar（RAdio Detection And Ranging）的音译词，意为“无线电检测和测距”，即利用无线电波来检测目标并测定目标的位置，这也是雷达设备在最初阶段的功能。典型的雷达系统如图1.1，它主要由发射机，天线，接收机，数据处理，定时控制，显示等设备组成。利用雷达可以获知目标的有无，目标斜距，目标角位置，目标相对速度等。现代高分辨雷达扩展了原始雷达概念，使它具有对运动目标(飞机，导弹等)和区域目标(地面等)成像和识别的能力。雷达的应用越来越广泛。

图1.1：简单脉冲雷达系统框图

假设理想点目标与雷达的相对距离为R，为了探测这个目标，雷达发射信号s(t),电磁波以光速C向四周传播，经过时间RC后电磁波到达目标，照射到目标上的电磁波可写成：R)。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为CRs(t)，其中为目标的雷达散射截面（Radar Cross Section ,简称RCS），反映目标对CR电磁波的散射能力。再经过时间RC后，被雷达接收天线接收的信号为s(t2)。

Cs(t如果将雷达天线和目标看作一个系统，便得到如图1.2的等效，而且这是一个LTI（线性时不变）系统。

图1.2：雷达等效于LTI系统

等效LTI系统的冲击响应可写成:

LFM脉冲压缩雷达仿真

h(t)(t)

(1.1)

iii1MM表示目标的个数，i是目标散射特性，i是光速在雷达与目标之间往返一次的时间，i2Ri

(1.2)c式中，Ri为第i个目标与雷达的相对距离。

雷达发射信号s(t)经过该LTI系统，得输出信号(即雷达的回波信号)sr(t)：

sr(t)s(t)*h(t)s(t)*i(ti)is(ti)

(1.3)

图1.3：雷达回波信号处理

s(t)的匹配滤波器hr(t)为：

hr(t)s*(t)

(1.4)于是，so(t)sr(t)*hr(t)s(t)*s*(t)*h(t)

（1.5）对上式进行傅立叶变换：

So(jw)S(jw)S*(jw)H(jw)|S(jw)|H(jw)2

(1.6)如果选取合适的s(t)，使它的幅频特性|S(jw)|为常数，那么1.6式可写为：

So(jw)kH(jw)

(1.7)其傅立叶反变换为：

so(t)kh(t)ki1Mi(ti)

(1.8)so(t)中包含目标的特征信息i和i。从 so(t)中可以得到目标的个数M和每个目标相对

LFM脉冲压缩雷达仿真

雷达的距离：

Riic

（1.9）2

这也是线性调频（LFM）脉冲压缩雷达的工作原理。二．线性调频（LFM）信号

脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation）信号,接收时采用匹配滤波器（Matched Filter）压缩脉冲。

LFM信号(也称Chirp 信号)的数学表达式为：

tj2(fctKt2)2s(t)rect()e

(2.1)

T式中fc为载波频率，rect()为矩形信号，tTtt11

（2.2）

rect()TT0,elsewiseKB，是调频斜率，于是，信号的瞬时频率为fcKt(TtT)，如图2.1

22T

图2.1 典型的chirp信号（a）up-chirp(K>0)（b）down-chirp(K<0)将2.1式中的up-chirp信号重写为：

s(t)S(t)e式中，j2fct

（2.3）

tjKt2S(t)rect()e

（2.4）

LFM脉冲压缩雷达仿真

是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需考虑S(t)。以下Matlab程序产生2.4式的chirp信号，并作出其时域波形和幅频特性，如图2.2。

仿真结果显示：

图2.2：LFM信号的时域波形和幅频特性

三． LFM脉冲的匹配滤波

信号s(t)的匹配滤波器的时域脉冲响应为：

h(t)s*(t0t)

（3.1）

t0是使滤波器物理可实现所附加的时延。理论分析时，可令t0＝0，重写3.1式，h(t)s*(t)

（3.2）将2.1式代入3.2式得:

LFM脉冲压缩雷达仿真

tjKt2ej2fct

(3.3)

h(t)rect()eT

图3.1：LFM信号的匹配滤波

如图3.1,s(t)经过系统h(t)得输出信号so(t)，so(t)s(t)*h(t)

s(u)h(tu)duh(u)s(tu)du

当0tT时，uj2fcujK(tu)2tuj2fc(tu)jKu2erect()eerect()eduTTT2s0(t)tT2ejKtej2Ktuduej2KtuT2ej2fct

(3.4)

j2KttT22

ejKt2当Tt0时，sinK(Tt)tj2fcteKttT2s0(t)T2ejKtej2Ktuduej2KtutT2j2fct

(3.5)ej2KtT22

ejKt2合并3.4和3.5两式： sinK(Tt)tj2fcteKttsinKT(1)tTrect(t)ej2fct

（3.6）

s0(t)TKTt2T3.6式即为LFM脉冲信号经匹配滤波器得输出,它是一固定载频fc的信号。当tT时，包络近似为辛克（sinc）函数。

S0(t)TSa(KTt)rect(tt)TSa(Bt)rect()

（3.7）2T2T-567

LFM脉冲压缩雷达仿真

（3）线性调频脉冲压缩雷达仿真程序LFM_radar

仿真程序模拟产生理想点目标的回波，并采用频域相关方法（以便利用FFT）实现脉冲压缩。函数LFM_radar的参数意义如下：

T：chirp信号的持续脉宽； B：chirp信号的调频带宽；

Rmin：观测目标距雷达的最近位置； Rmax：观测目标距雷达的最远位置；

R：一维数组，数组值表示每个目标相对雷达的斜距； RCS：一维数组，数组值表示每个目标的雷达散射截面。得到的仿真结果如图4.2。（4）分辨率(Resolution)仿真

改变两目标的相对位置，可以分析线性调频脉冲压缩雷达的分辨率。仿真程序默认参数的距离分辨率为：

C3108R5m

(4.1)

2B230106图4.3为分辨率仿真结果，可做如下解释：(a)图为单点目标压缩候的波形；

(b)图中，两目标相距2m，小于R，因而不能分辨；

(c)图中，两目标相距5m，等于R，实际上是两目标的输出sinc包络叠加，可以看到他们的副瓣相互抵消；

(d)－(h)图中，两目标距离大于雷达的距离分辨率，可以观察出，它们的主瓣变宽，直至能分辨出两目标。

LFM脉冲压缩雷达仿真

图4.2：仿真结果

LFM脉冲压缩雷达仿真

图4.3：线性调频脉冲压缩雷达分辨率仿真

第四篇：脉冲多普勒雷达测速仿真汇总

脉冲多普勒雷达测速仿真

任务书

雷达进行PD测速主要是利用了目标回波中携带的多普勒信息，在频域实现目标和杂波的分离，它可以把位于特定距离上、具有特定多普勒频移的目标回波检测出来，而把其他的杂波和干扰滤除。因此要求雷达必须具备很强的抑制杂波的能力，能在较强的杂波背景中分辨出运动目标的回波。

如今，不管是在军用还是民用上，雷达都在发挥着它很早重要的作用，与早期雷达采用距离微分方法测速相比，基于脉冲多普勒理论的雷达测速技术具有实时性好、精度高等优点。特别是现代相控阵技术在雷达领域的应用，实现了波束的无惯性扫描和工作方式的快速切换，更便于应用脉冲多普勒技术进行雷达测速。

本篇课程设计目的在于介绍脉冲多普勒雷达测速的原理，并对这种技术进行介绍和仿真。

脉冲多普勒雷达测速仿真

摘要

脉冲多普勒(PD)雷达以其卓越的杂波抑制性能受到世人瞩目。现代飞行器性能的改进和导航手段的加强，使其能在低空和超低空飞行，因此防御低空入侵己成重要问题，由此要求机载雷达，包括预警机雷达和机载火控雷达具有下视能力，即要求能在强的地杂波背景中发现微弱的目标信号，所以现代的预警机雷达和机载火控雷达皆采用PD体制。脉冲多普勒雷达包含了连续波雷达和脉冲雷达两方面的优点，它具有较高的速度分辨能力，从而可以更有效地解决抑制极强的地杂波干扰问题;此外，脉冲多普勒雷达能够同时敏感地测定距离和速度信息;能够利用多普勒处理技术实现高分辨率的合成孔径图像;而且亦具有良好的抗消极干扰能力和抗积极干扰能力。

本文介绍了脉冲多普勒雷达测速的原理，信号处理。并用matlab简单的仿真了雷达系统对信号的处理.关键词：脉冲多普勒雷达恒虚警脉冲压缩线性调频

Abstact

Pulse Doppler(PD)radar is famous for it`s outsdanding clutter suppression.Modern aircraft`s function and GPS has been strengthen.now.it makes the aircraft can fly lower and lower.So.nowadays,Defensing.Low altitude invasion has been an important problem.so we require airborne radar.Early warning radar and airborne fire control radar have the ability to look down.That is to say.The radar is be required the ability to find Weak target signal in the strong Groung clutter.So.The modern airborne early warning radar and airborne fire control radar use the PD system.Pulse Doppler(PD)radar concludes two adervantages of Continuous wave radar and impulse radar.It has a higher velocity resolution.thus it can effectively.soveing the problem of strong ground clutter.what`s more.Pulse Dppler(PD)radar can Sensitive text the Distance and speed on the same time.Itcan use Doppler processing technology to realise Synthetic aperture images with high resolution.This article sinply introduced principle of pulse Doppler radar and signal

脉冲多普勒雷达测速仿真

processing.Using matlab to simulation The signal processing of radar system.Linear frequency modulation.Keywords:Pulse Doppler(PD)radar.Constant false alarm rate.pulse compression.脉冲多普勒雷达测速仿真

一．脉冲多普雷达简介··········································1 1,多普勒效应··············································1

二、多普勒测速原理············································2

三、多普勒雷达简介············································4

四、多普勒雷达工作原理········································6

五、PD雷达信号处理仿真·······································8 5.1、正交双通道处理·······································9 5.2、脉冲压缩·············································10 5.3、线性调频信号的脉冲压缩·······························12 5.4、巴克码信号的脉冲压缩·································14 5.5、恒虚警处理···········································14 5.5.1、单元平均恒虚警处理(CA-CFAR)····················16 5.5.2、平均选大恒虚警处理(GO-CFAR)····················16 5.6、动目标检测(MTD)模型··································19

六、总结与展望················································20 参考文献······················································21

脉冲多普勒雷达测速仿真

二、脉冲多普雷达简介

1,多普勒效应

多普勒效应是指当发射源和接收者之间有相对径向运动时，接收到的信号将发生变化。这一物理现象首先由奥地利物理学家多普勒于1842年发现的。多普勒雷达的工作原理是以多普勒效应为基础的，多普勒效应不仅仅适用于声波，它也适用于所有类型的波，包括光波、电磁波。

当无线电波在行进的过程中碰到物体时，该无线电波会反弹，而且反弹回来的波其频率以及振幅会随着碰撞物体的运动状态而变化。若无线电波是固定不动的，那么反弹回来的电波频率是不会发生改变的，然而，若物体是朝着无线电发射的方向前进，此时所反弹回来的电波会被压缩，因此该电波的频率也会随之增加。反之，若远离波源运动时，接收到的频率较波源的实际频率降低。频率升高或者降低的数值为多普勒频率，多普勒效应的数学推理如下: 设某一波源的频率为f，波长为入，它们与波的传播速度的关系为: f  v

收

（1.1）

若波源以速度W向静止的接物体运动时，在接收物体处所接收到的频率为: f'v'vf

vw

(1.2)若波源以速度W背向静止的接收物体运动时，在接收物体处所接收到的频率为: vvf“f”vw到的频率为:

(1.3)同理，当波源不动时，若物体以速度W向波源运动时，在接收物体处所接收f'v'vwf v

(1.4)当波源不动时，若物体以速度W背向波源运动时，在接收物体处所接收到的频率为：

脉冲多普勒雷达测速仿真

vvwf“f ”v

(1.5)由上述的数学分析可知，若波源与接收物体之间存在相对运动时，接收的频率就不同与发射的频率。即两者之间的距离缩短(相向运动)时，接收频率高于发射频率;两者之间的距离增大时(反向运动)，接收频率低于发射频率。这种由于相对运动引起的频率变化，称为“多普勒效应”。

二、多普勒测速原理

假若有一个运动目标，利用多普勒效应，在天线场的有效范围内，如图2.1所示:

图2.1多普勒运动信号模型

天线1发射电磁波到运动物体上，发射频率为f0，速度为光速C，运动物体的速度为V，则运动物体实际接收到的频率为: f1cvf0 v

（2.1）

同理，由多普勒效应知，电磁波从运动物体反射回来，天线2所接收到的频率为: f2vf1

cv

（2.2）

所以，天线2接收到的频率与发射频率之间的关系如下:

脉冲多普勒雷达测速仿真

f2vcvcvf0f0

cvvcv

（2.3）

式(2.3)说明:由发射天线1发射的频率为fo的电磁波，相对天线以径向速度V作离向动，由于多普勒效应，接收天线2接收到的电磁波的频率将是发生两次多普勒效应后的频率，即变为f 2。将相对运动所引起的接收频率与发射频率之间的差频称为多普勒频率，用fa表示，则: fdf0f22v2vf0

（2.4）

由式((2.4)分析得知，多普勒频率与相对天线的径向速度成正比，只要能测出fd就可以求出V，这就是多普勒测速雷达测速的基本原理。

根据前面的分析与推导，可以对多普勒频率作出这样的定义: 2vrfdf0（2.5）c式中:fd：多普勒频率

Vr：运动目标的速度 C：光速 f0：发射波频率由式(2.5)我们可以得到: 2vrvrf0 c

（2.6）

从式(2.6)可以看到其它变量都是己知的，只要我们测出fd了，就可以计算出被测目标的速度。下面来分析当雷达与目标有相对运动时，雷达接收信号的特

脉冲多普勒雷达测速仿真

征。在这里我们设目标为理想“点”目标，即目标尺寸远小于雷达分辨单元。

三、多普勒雷达简介

脉冲多普勒雷达的应用和发展是雷达理论与技术的发展，特别是新型电子器件和数字信号处理技术发展的结果。目标环境变化，以及下视，下射，反低空突防任务的迫切需要，以及为了获得优良的抑制严重的杂波干扰的性能，对脉冲多普勒雷达信号源的稳定度，寄生调制提出了极高的要求，对其天线旁瓣电平提出甚严的限制，同时还要求具有高速运算能力和大存储空间的处理机。一部实际的脉冲多普勒雷达几乎包含了现代雷达的主要波形、测角体制和先进的信号处理技术，采用了各种现代优化的设计思想。因此，研制脉冲多普勒雷达是一个较为艰巨的工作。那么为了更好且更经济地研究其技术、战术性能和抗干扰性能，采用计算机仿真方式进行研究是最为合理地，因此，对脉冲多普勒雷达的计算机仿真成为一项迫切的任务。

关于PD雷达的精确定义，1970年M.I.斯科尔尼克曾以三点特征进行描述： 1)具有足够高的脉冲重复频率(PRF)，以致不论杂波或所观测的目标都没有速度模糊;2)能实现对脉冲串频谱单根谱线的多普勒滤波，即频域滤波;3)由于PRF很高，通常对所观测的目标产生距离模糊。

上述定义仅适用于高PRF的PD雷达，随着PD雷达技术的发展，对PD雷达的定义有所延伸。一些PD雷达的设计允许工作在低、中高频上来满足不同的需要。目前关于PD雷达的定义为:能实现对信号脉冲串频谱的单根谱线滤波，具有对目标进行速度分辨力的雷达称为PD雷达。

根据PD雷达的系统重复频率不同，可以对其进行分类，令雷达与目标的距离为RT，雷达与目标的相对速度为从，最大不模糊距离为R man，最大不模糊速度为Vmah，则可将PD雷达分为高、中、低PRF三类脉冲雷达:

脉冲多普勒雷达测速仿真

图3.1 脉冲多普勒雷达分类

可以看到，低PRF雷达无距离模糊，有速度模糊;中PFR既有距离模糊，也有速度模糊，高PRF无速度模糊，但有距离模糊。对于PD雷达，选择合适的PRF是一个十分重要的问题，需要考虑实际的应用和设备的复杂性。

PD雷达采用相参技术从相参到非相参，这是雷达技术的飞跃。相参，是从光学领域引申而来的术语，在雷达中的含义是指目标回波信号与发射信号之间应保持严格的相位关系，并用以提取目标的有关信息。早期的普通脉冲雷达都是非相参结构，它们只是利用发射脉冲与回波脉冲之间的时间差提取目标的距离信息，当与目标相等距离上，存在着大量的干扰物体的反射波(地物和海浪杂波，云雨气象杂波以及投放物金属箔条杂波等)时，微弱的目标回波淹没在其中，而使雷达失效。非相参的普通脉冲体制的雷达，无论架设在地面，舰船或是飞机上都无法对付这些来自低空的袭击。在光学中早己为人们所熟知的相参概念，此刻却为雷达技术带来了生机。当用时间差无法区别目标与背景干扰时，速度差却很容易把它们区别开来。

对于地面假设的雷达，与周围地物不存在相对运动，因此从这些地物来的雷达回波信号的频率与反射信号的频率相同;而从运动目标来的回波信号频率，则有别于发射信号，而且这种差别额大小正比于目标的速度。对于机载雷达，地物背景与雷达虽然存在着相对运动，但是这种相对运动的速度和方向是可知的，并有别于机载雷达和空中运动目标的相对速度。以雷达载机同运动目标和地物的速度差为基础发现运动目标的技术途径是可行的。但是要具体实现这种发射信号和

脉冲多普勒雷达测速仿真

回波信号的有效相参并不十分容易。它要求发射信号具有很高的稳定性，否则微小的回波频率变化是无法察觉的。主振放大式发射机就是目前广泛用来获得高功率、高稳定的发射信号的主要设备，尽管其成本高，构造复杂，但这是雷达实现全相参所必须的。雷达实现了全相参，它的性能和应用领域可能大为扩展。

四、多普勒雷达工作原理

由于目标和干扰物相对于雷达的径向速度不同，回波信号也有不同的多普勒频率，可用频域滤波的方法选出目标的多普勒频率谱线，滤除干扰杂波的谱线，使雷达从强杂波中分离和检测出目标信号。为实现这一目的，一方面PD雷达的发射脉冲信号必须有稳定的相参性能，通常采用主振放大式发射机，另一方面在接收机的信号处理中，把每一脉冲重复周期分成若干个距离门，再用多普勒频率范围内的窄带滤波器组对信号和杂波进行过滤。下面先扼要讲述PD雷达信号工作的基本原理，其中信号处理模块、几种典型的干扰，数据处理模块中的解模糊算法，滤波算法是本论文重点讲述的内容，典型PD雷达原理方框图4.1为例:图中只画出了一个接收通道。为了连续跟踪一个目标，还包括速度跟踪、距离和角度跟踪设施。PD雷达可以把位于特定距离上，具有特定多普勒频移的目标回波检测出来，把其它的杂波和干扰滤除。

主振放大式相参发射机:这种发射机信号质量高，可以满足PD雷达对发射信号的频率稳定度和频谱纯度的要求。影响PD雷达性能的主要因素是信号的短期不稳定性，PD雷达的频率稳定度要达到很高的程度才能达到相参测速的要求。距离门电路完成距离选通，雷达可以有多个距离通道。

单边带滤波器是一个带宽近似等于脉冲重复频率的带通滤波器，主要作用是从回波频谱中只滤出一根谱线和它所对应的多普勒频移分量，使后面的各种滤波处理在单根谱线上进行，通常单边带滤波设置在中频，选取回波频谱的中心谱线。

主杂波抑制电路对干扰很强的主瓣杂波进行抑制;机载PD雷达的高度杂波是由地面的垂直反射所形成的杂波，通常可以采用单独的固定频率抑制滤波器(零频率滤波器)来滤除它。

脉冲多普勒雷达测速仿真

信号处理机:PD雷同常规脉冲雷达的主要区别在于PD雷达利用了目标回波中携带的多普勒信息，在频域内实现目标和杂波的分离，从而可以从很强的地物杂波背景中检测出动目标，并对其精确测速。简而言之，PD雷达信号处理机的核心。就是一个窄带滤波器组，它滤除了各种干扰及杂波，保留了所需的目标信号。由于雷达从天线到视频的全部过程都是线性的，运动目标信息在接收机中无畸变地保存到视频级，因此多普勒滤波器可以在视频(零中频)实现。由于视频滤波比中频滤波易于实现，数字处理在视频进行量化也比在中频容易，所以多普勒滤波将变得更加简单可靠。这是因为，数字式滤波器具有高可靠性、抗干扰性和灵活性等优点，精度也好。基本的数字滤波器的原理是利用快速傅立叶变换求取信号频谱，从而为下一步的处理提供回波信号的频谱分布信息。

数据处理机:数据处理机的主要功能，一个是对工作方式、扫描图形的产生、脉冲重复频率的选择、杂波频潜的预测进行控制;另一个功能是对数据相关和滤波、距离跟踪、角度跟踪、雷达罩和天线角度误差修正进行数据处理;再一个功能是对系统进行性能监视和机内自检。

PD雷达必须采取恒虚警率((CFAR)处理技术以便防止干扰增大时虚警概率过高，努力使得当噪声、杂波和干扰功率或其他参数发生变化时，输出端的虚警概率保持恒定;速度、距离和角度跟踪系统实现对目标的有效跟踪。

PD雷达系统中的很多理论己经是比较成熟的理论，在本文中主要进行雷达回波信号、信号处理中相关算法、数据处理相关算法、及PD雷达几种主要干扰技术的建模仿真，本文主要做PD雷达测距、测速的研究。

多普勒雷达工作原理框图如图4.1:

脉冲多普勒雷达测速仿真

图4.1 典型脉冲多普勒雷达工作原理框图

五、PD雷达信号处理仿真

PD雷达信号处理模块的主要框图:

图5,1 雷达信号处理框图

信号处理器主要是用来监测目标并利用一定的方法来抑制各种杂波和人为干扰所产生的不希望有的信号，处理后的视频输出信号再与某个门限比较。若信号超过监测门限，便判断为“发现目标”，然后把目标信号输送下一级处理器进行监测，以便测量出目标的距离，角度，径向速度和其它一些目标特性。

脉冲多普勒雷达测速仿真

下面对于其主要模块进行仿真。

5.1、正交双通道处理

为将中频信号转换为视频信号进行数字处理，采用低通滤波法进行相干检波，低通滤波法是一种完全仿照传统的模拟正交采样的实现方法，只是将频移放在了A/D变换之后。这样混频和滤波都是由数字系统来实现的，其原理框图如图5.2所示。

图5.2 正交双通道处理框图

首先将中频输出信号进行数字混频，这样就将正频谱的中心移到了零频，时域信号也应的分解为实部和虚部，再让混频后的信号经过低通滤波器，滤除高频分量，即可得到所需的基带正交双路信号。

脉冲多普勒雷达测速仿真

图5.3 对信号滤波处理后

5.2、脉冲压缩

匹配滤波器是以输出最大信噪比为准则的最佳线性滤波器。当滤波器的频率响应H(f)为信号频谱S(f)的复共扼时，称之为信号的匹配滤波。根据匹配滤波器理论，在白噪声背景下，滤波器输出端信号噪声功率比的最大峰值为2ElNo，即当噪声功率谱密度给定后，决定雷达检测能力的是信号能量E。使用简单矩形脉冲信号，可以用增大脉冲宽度了来提高信号能量，但是这样使得距离分辨力降低。因此，提高雷达的探测能力和保证必需的距离分辨力这对矛盾很难解决，有必要去寻找和采用较为复杂的信号形式。

匹配滤波器输出信号是波形的自相关函数，它是信号功率谱的傅里叶反变换。由于脉宽和频谱宽度之间成反比关系，因此B越大，距离分辨力越好，距离分辨力取决于所用信号的带宽B。对于简单矩形脉冲，信号带宽B与其脉冲宽度:满足Bt约等于1的关系，因此用宽脉冲时必然降低其距离分辨力。如果在宽脉冲内采用附加的幅度或相位调制，以增加信号带宽B，那么当接收时用匹配滤波器进行处理，可将长脉冲压缩到1/B宽度，这样既可使雷达用长的脉冲去获得大的能量，同时又可以得到短脉冲所具备的距离分辨力。这种信号称为脉冲压缩信号或大时宽带宽信号，这种方法称为脉冲压缩。因为脉冲内有附加调制，其脉宽:和带宽B的乘积大于1，压缩比，一般采用BT约等于1.脉冲多普勒雷达测速仿真

幅度调制会降低发射信号的平均功率，不能最大程度地利用发射管的效能，使雷达威力下降。相位调制是常用的宽带信号产生方法，常用的相位调制信号有线性调频信号和二相编码信号，下面具体分析对这两种信号的脉冲压缩。脉冲压缩滤波器是雷达视频信号的匹配滤波器，是在时间轴上进行的，可以针对一个脉冲重复周期或整个脉冲压缩处理。

匹配滤波器的实现是脉冲压缩处理的核心部分，主要分别可以通过时域卷积或频域相乘实现。所谓匹配滤波器就是这样一种最佳线形滤波器，在输入为确知信号加白噪声的情况下，所得输出信噪比达到最大。匹配滤波器的冲激响应函数为:

要实现数字脉冲压缩，需先对模拟的线性调频信号进行A/D变换，经过采样量化后的信号简记为s(n)。

频域处理公式为y(n)=IFFT(S(f)H(f))=IFFT(FFT(s(n))FFT(h(n)))。通常情况下，基于FFT算法特有的高速性，频域实现比直接时域卷积运算要快。用快速傅立叶变换实现的脉冲压缩频域处理过程如图5.4。

图5.4 脉冲压缩频域处理

脉冲多普勒雷达测速仿真

图5.5 脉冲压缩频域处理

5.3、线性调频信号的脉冲压缩

线性调频信号的数字脉冲压缩处理可在时域，也可在频域进行。由于高速A/D变换器、大规模集成电路技术以及FFT的应用，使宽带信号的实时处理成为可能，采用DSP处理的频域数字脉冲压缩具有处理速度高、工作稳定、重复性好、灵活性高的优点。所以，在仿真中对于线性调频信号在频域进行压缩。

线性调频信号脉冲压缩后的旁瓣较高，最大旁瓣为13.2dB，一般不能满足系统要求，常采取窗函数抑制旁瓣，可根据对旁瓣和主瓣宽度的要求选择窗函数，其中海明窗较为常用。加窗处理同样可在时域进行也可在频域进行，对于时宽带宽积较大的信号时域和频域加窗效果接近，但时宽带宽积较小的信号时域加窗的效果优于频域加窗，所以仿真中采用时域加窗频域压缩的方法。图5.6为线性调频信号脉冲压缩处理的流程框图。

脉冲多普勒雷达测速仿真

图5.6 线性调频信号脉压处理框图

由于线性信号存在距离速度祸合，当目标回波中存在多普勒频率时，经脉冲压缩后会产生距离误差。多普勒频率大于零时脉冲压缩后的主瓣向左偏移，即距离变小，多普勒频率小于零时脉冲压缩后的主瓣向右偏移，即距离变大。在精密跟踪雷达中需进行误差补偿，用于补偿由于目标的运动引起的测距误差。

误差补偿模型仅适用于单目标环境，而且仅是在雷达工作在跟踪状态时才做误差补偿，雷达工作在搜索状态时不做误差补偿处理。

图5.7 线性调频信号脉冲压缩处理结果

脉冲多普勒雷达测速仿真

5.4、巴克码信号的脉冲压缩

巴克码脉冲压缩器采用和线性调频脉冲压缩相同的方法，时域加窗频域压缩。当不存在多普勒失谐时，巴克码经过相关处理后输出的主瓣高度等于码长N，旁瓣高度为1，即主副瓣比为N。当采用最大长度13时，主副瓣比约为22.3dB需要抑制旁瓣以提高目标检测能力。如果回波信号有多普勒频移，主瓣幅度下降，旁瓣幅度上升，应设法克服由多普勒失谐带来的影响。

图5.8 巴克码处理过的信号

5.5、恒虚警处理

PD雷达用于机载下视或类似的条件下，由不同的地物回波所形成的杂波的强度和分布情况极为复杂。为了在复杂的杂波环境中检测出所关心的运动目标回波，要求PD雷达必须采用某种恒虚警率((Const False Alarm Rate, CFAR)处理技术，以便在杂波环境变化时，防止雷达的虚警概率发生太大的变化，同时保证一定的检测概率。同常规脉冲雷达不同，PD雷达一般是在经过多普勒滤波后在

脉冲多普勒雷达测速仿真

频域进行检测，因此PD雷达所采用的CFAR处理方式也有所不同。

参量法CFAR处理是基于数理统计中的参数估计理论，它的前提是要知道杂波幅度分布的类型，未知的只是分布函数中的参数。对于不同的杂波分布类型，实现CFAR处理的方法也不同，但他们的基本思想是一样的，即在假设检测单元附近各背景单元的杂波输出为独立同分布的前提下，通过处理各背景单元的输出而对杂波分布中的参数进行估计，然后构造一个分布不依赖杂波强度的新的统计量去与门限比较，从而达到CFAR的目的。

图5.9 恒虚警处理结果

雷达杂波可以看作是来自大量独立反射体的回波信号的总和。这些回波信号的幅度和相位都是随机的，因此合成后杂波的幅度是一个随机变量，其统计特性随杂波环境的不同而变化。若粗糙地面的散射单元是均匀的，则根据中心极限定理，合成后的杂波近似服从高斯分布，经过包络检波后，杂波包络幅度服从瑞利分布，其概率密度函数为

脉冲多普勒雷达测速仿真

xx2exp[]x0

2f(x)220,其他（5.1）

瑞利分布是描述雷达杂波时适用范围最广的一种分布。对于工作在无杂波区的高重复频率的PD雷达而言，由于接收机噪声一般服从高斯分布，所以噪声包络幅度也可以用瑞利分布来描述。下面针对瑞利分布杂波介绍CFAR技术。5.5.1、单元平均恒虚警处理(CA-CFAR)瑞利分布杂波幅度的数字特征为

E(x)2,E(x2)22

（5.2）

其中E(X)为均值，E(Xz)表示杂波的功率。当方差变化时，X的概率密度函数也将随之变化，对于给定的门限UT，相应的虚警概率

pfaxJT2x2VTexp[]dexp[]x2

2（5.3）

222也发生变化。若定义一个新的随机变量Y，令

YX

（5.4）

则容易求得Y的概率密度函数为

fY(y)yexp(y22),y0

（5.5）

CA-CFAR检测在杂波边缘中要引起虚警率的上升，而在多目标环境中将导致检测性能下降，这些不足促进了对其它CFAR方案的寻求。新的CFAR方案应该能够区分干扰和主目标，平均选大GO(Greatest o介CFAR是一种满足上述要求的修正方案。

5.5.2、平均选大恒虚警处理(GO-CFAR)

脉冲多普勒雷达测速仿真

平均选大在杂波边缘环境中能保持较好的控制虚警性能。图5.10是单元平均选大CFAR.模型，通过两侧L/2个距离单元的数据平均值估计杂波功率，用杂波功率对所检测的距离单元数据进行归一化并乘以门限乘子后作为检测门限，与检测门限比较，超过门限判断为有目标，低于门限判断为无目标。

图5.10 GO-CFAR实现框图

呈现在脉冲多普勒接收机输出端的杂波分布特征是各种杂波的混合，在距离一多普勒图的范围内这个混合杂波的变化可能很大，并且强点状杂乱回波的位置也是随机的。因此，必须对邻近测试单元的多个单元组成的当地参考窗口进行检查，用这种方法来估计测试单元干扰功率。根据PD雷达工作频率的不同，在距离一多普勒图中设置参考窗口的方式也不同。对于低重复频率，需要在距离维中设置参考窗口:对于高重复频率，需要在多普勒频率维设置窗口;然而对于中重复频率的雷达，则需要设置两维窗口。

以上分析的是在瑞利杂波情况下的恒虚警处理过程和相关电路。雷达在实际工作环境中所接收的杂波可能并不是瑞利分布，那么通过上面处理步骤得到的恒虚警效果可能不明显。某种具体的恒虚警电路一般只适用于对应的特定杂波情况，对于不同的杂波，应该对恒虚警处理电路做必要的调整才能满足要求。总的原则是根据杂波的具体模型，通过参数估计和适当的变换构造出新的随机变量，使它的概率密度函数与输入杂波强度无关。

脉冲多普勒雷达测速仿真

5.6、动目标检测(MTD)模型

由于MTI的杂波抑制能力有限，为了弥补MTI的缺陷，为此在MTI后串接一窄带滤波器组覆盖整个频率的范围。当不同的目标有不同的速度时，应的不同多普勒频率对应不同的窄带滤波器输出，因此可以利用MTD测出多普勒频移，确定目标的速度，这样，滤波器组实现了速度分辨和精确测量的作用。具有N个输出的横向滤波器(N个重复周期和N-1根延迟线)，经过各重复周期的不同加权并求和后，即可实现上面所述的窄带滤波器组。其原理结构框图如5.11。

图5.11 MTD横向滤波器结构

当满足采样定理的条件时，一个时间函数的取样序列经过DFT处理之后，输出为该信号的频谱取样。我们可以将每根谱线看成一个窄带滤波器的输出，这就是DFT的滤波特性。因而利用DFT的滤波特性可以形成窄带多普勒滤波器组，其振幅特性可以表示为|fg|。

|Hk(f)||sin[N(fTrKN)]|,K0,1,2...,N1（5.6）

sin[(fTrKN)]

脉冲多普勒雷达测速仿真

图5.12即N=8, fr =lOKHz时DFT等效横向滤波器振幅特性。

图5.12 DFT等效横向滤波器的幅频响应

现在用来产生数字滤波器组的最普遍的方法是FFTo FFT是离散傅立叶变换的快速算法，通常采用每个脉冲同一个距离单元的一组数据FFT来实现滤波器组。

其实使用数字方法计算离散信号的频谱，每个固定频率分量的输出就相当于中心频率在此固定频率的窄带滤波器输出。

脉冲多普勒雷达测速仿真

图5.13 MTD求模处理结果

六、总结与展望

本文简单的对脉冲多普勒雷达做了简单介绍，对脉冲多普勒雷达系统中的信号处理部分，有源干扰及数据处理中的典型算法进行建模仿真的探索，脉冲多普勒雷达系统仿真是一个内容复杂、涉及广泛的研究课题，采用的仿真软件不同，仿真方法不同，实现的仿真系统的功能、复杂程度也会有所不同。本文采用功能强大的计算仿真软件MATLAB对其进行仿真，主要目的在于使仿真模块的程序易读易懂，结构简单，能够基本反映出该体制雷达主要处理过程的工作原理。

对雷达系统仿真是一个非常有价值的研究领域，随着计算机仿真技术的不断提高，仿真领域将不断扩大。由于本人学识有限，时间有限，本文还存在很多的缺点不足，例如，本文仅尝试对雷达信号处理部分及部分干扰的建模，系统的跟踪模块未能完善，雷达数据处理只涉及到算法的研究，并未将其与前面的信号处理组合到一个完整的系统，希望本文些许的分析可以对以后的学友提供一些帮助，更希望以后的研究能在此基础上不断完善，如果能够结合到实际系统中，那么其仿真系统的意义将更加深广。

脉冲多普勒雷达测速仿真

本人水平有限，论文中错误及不当之处，恳请各位老师批评指正。参考文献

[1]刘丽华，董天临连续波多普勒测速雷达射频前端电路设计与仿真.电子科技.[2]张武娟.雷达测速在列车运行中的研究与应用【D】:[硕士论文].长沙:中南大学，2008.[3]周晶晶.雷达测速仪数字信号处理系统的设计【D】:【硕士论文】.西安:西北工业大学，2006.[4]汪源源.多普勒信号分析技术的比较研究.复旦学报(自然科学版)，1996，2:53-61。[5]张伟.基于ARM的雷达信号处理系统的研究【D】:【硕士论文】.南京:南京理工大学，2008 [6]刘峰，初韵辉，许家栋等.一种PO雷达半实物仿真系统的研究.西北工业大学学报.2002.[7]韦伟.脉冲多普勒雷达信号处理关键技术研究.南京理工大学硕士学位论文.2004.

第五篇：大作业-雷达线性调频脉冲压缩的原理及其MATLAB仿真(DOC)

LFM脉冲压缩雷达仿真

西安电子科技大学

线性调频（LFM）脉冲压缩雷达仿真

概述：雷达工作原理

图1.1：简单脉冲雷达系统框图

一．线性调频（LFM）脉冲压缩雷达原理

假设理想点目标与雷达的相对距离为R，为了探测这个目标，雷达发射信号s(t),电磁波以光速C向四周传播，经过时间RC后电磁波到达目标，照射到目标上的电磁波可写成：R)。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为CRs(t)，其中为目标的雷达散射截面（Radar Cross Section ,简称RCS），反映目标对CR电磁波的散射能力。再经过时间RC后，被雷达接收天线接收的信号为s(t2)。

Cs(t如果将雷达天线和目标看作一个系统，便得到如图1.2的等效，而且这是一个LTI（线性时不变）系统。

LFM脉冲压缩雷达仿真

西安电子科技大学

雷达的距离：

Riic

（1.9）2

这也是线性调频（LFM）脉冲压缩雷达的工作原理。二．线性调频（LFM）信号

脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率，保证足够大的作用距离；而接收时采用相应的脉冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分辨率之间的矛盾。

脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation）信号,接收时采用匹配滤波器（Matched Filter）压缩脉冲。

LFM信号(也称Chirp 信号)的数学表达式为：

tj2(fctKt2)2s(t)rect()e

(2.1)

T式中fc为载波频率，rect()为矩形信号，tTtt11

（2.2）

rect()TT0,elsewiseKB，是调频斜率，于是，信号的瞬时频率为fcKt(TtT)，如图2.1

22T

图2.1 典型的chirp信号（a）up-chirp(K>0)（b）down-chirp(K<0)将2.1式中的up-chirp信号重写为：

s(t)S(t)e式中，j2fct

（2.3）

tjKt2S(t)rect()e

（2.4）

LFM脉冲压缩雷达仿真

西安电子科技大学

是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需考虑S(t)。以下Matlab程序产生2.4式的chirp信号，并作出其时域波形和幅频特性，如图2.2。>> %%线性调频信号的程序 T=10e-6;B=30e6;K=B/T;

图2.2：LFM信号的时域波形和幅频特性

LFM脉冲压缩雷达仿真

西安电子科技大学

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));

辛克函数（理论波形）Z1=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 signal after matched filter');subplot(212)

N0=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,1,2,3]);xlabel('Time in sec timesitB');ylabel('Amplitude,dB');title('Chirp signal after matched filter(Zoom)');仿真结果如图3.3：

图3.3：Chirp信号的匹配滤波

图3.3中，时间轴进行了归一化，（t/(1/B)tB）。图中反映出理论与仿真结果吻合良好。

LFM脉冲压缩雷达仿真

西安电子科技大学

第一零点出现在1（即（11）处，此时相对幅度-13.4dB。压缩后的脉冲宽度近似为BB1），此时相对幅度-4dB,这理论分析（图3.2）一致。2B

上面只是对各个信号复包络的仿真，实际雷达系统中，LFM脉冲的处理过程如图3.4。

图3.4： LFM信号的接收处理过程

图3.5：正交解调原理

图3.6：一种脉冲压缩雷达的数字处理方式

四：Matlab仿真结果

（1）任务：对以下雷达系统仿真。

雷达发射信号参数：幅度：1.0 信号波形：线性调频信号

频带宽度：30兆赫兹（30MHz）脉冲宽度：10微妙（20us）中心频率：1GHz（109Hz）雷达接收方式：正交解调接收

距离门：10Km~15Km 目标：

LFM脉冲压缩雷达仿真

西安电子科技大学

Tar1：10.5Km Tar2：11Km Tar3：12Km Tar4：12Km＋10m Tar5：13Km Tar6：13Km＋25m（2）系统模型：

图4.1：雷达仿真等效信号与系统模型（3）线性调频脉冲压缩雷达仿真程序LFM_radar

仿真程序模拟产生理想点目标的回波，并采用频域相关方法（以便利用FFT）实现脉冲压缩。函数LFM_radar的参数意义如下： T：chirp信号的持续脉宽； B：chirp信号的调频带宽；

Rmin：观测目标距雷达的最近位置； Rmax：观测目标距雷达的最远位置；

R：一维数组，数组值表示每个目标相对雷达的斜距； RCS：一维数组，数组值表示每个目标的雷达散射截面。在Matlab指令窗中键入：

LFM_radar(10e-6,30e6,10000,15000,[10500,11000,12000,12010,13000,13025],[1,1,1,1,1,1])得到的仿真结果如图4.2。

图4.2：仿真