雷达管制的雷达管制的一般过程

在提供雷达管制服务之前，应建立和保持对航空器雷达识别(Radar Identification)，直到雷达服务终止。雷达识别是指将某一特定的雷达目标或雷达位置符号与某特定航空器相关联的过程。如果此后丢失雷达识别应该立刻重新识别航空器或者终止雷达服务，并且相应地通知飞行员，在适当时发布指令。雷达识别必须通过下面至少一种方法建立。

（一数则盯）一次监视雷达识别程序(PSR)

（二） 二次监视雷达识别程序(SSR)

一、一次监视雷达识别程序(PSR)

一次雷达可以分成三类：

机场监视雷达（ASR：Airport Surveillance Radar），它的作用距离为100海里，主要是塔台管制员或进近管制员使用。

航路监视雷达（ARSR:Air Route Surveillance Radar），设置在航管控制中心或相应的航路点上。它的探测范围在250海里以上，高度可达13000米。它的功率比机场监视雷达大，在航路上的各部雷达把整个航路覆盖，这样管制员就可以对航路飞行的飞机实施雷达间隔。

机场地面探测设备（ASD），它的功率小，作用距离一般为1英里，主要用于特别繁忙机场的地面监控，它可以监控在机场地面上运动的飞机和各种车辆，塔台管制员用来控制地面车辆和起降飞机的地面运行，保证安全。它主要的作用是在能见度低的时候提供飞机和车辆的位置信息，由于它的价格较高，机场通常没有这种设备。

使用一次监视雷达时，可以通过下列程序识别航空器：

a.航空器起飞后，其雷达目标在起飞跑道端2千米以内被发现；

b.当观察到仅有一个雷达目标，与雷达视频图上某点的相对位置与航空器的报告位置一致，其显示的航迹与航空器报告的航迹或计划航迹相符。

c.观察到仅有一个雷达目标，按照指示做不小于3 0度的识别转弯；（注：使用此识别时应注意避免航空器转弯时脱离管制区或雷达覆盖区，或使航空器低于最低安全高度，或最低引导高度）

d.通过雷达识别移交。

二、二次监视雷达识别程序(SSR)

二次雷达也叫做空管雷达信标系统（ATCRBS：Air Traffic Control Radar Beacon System）。它最初是在空战中为了使雷达分辨出敌我双方的飞机而发展的敌我识别系统，当把这个系统的基本原理和部件经过发展后用于民航的空中交通管制后，就成了二次雷达系统。

管制员从二次雷达上很容易知道飞机的编号、高度、方向等参数，使雷达由监视的工具变为空中管制的手段，二次雷达的出现是空中交通管制的最重大的技术进展，二次雷达要和一次薯和雷达一起工作，它的主天线安装在一次雷达的上方，和一次雷达同步旋转。二次雷达发射的脉冲是成对的，它的频率是1030MHz，每一对脉冲之间的时间间隔是固定的，这个间隔决定了二次雷达的模式。目前民航使用的是两种模式，一种间隔为8毫秒，称为A模式；另一种间隔21毫秒，称为C模式。

二次雷达系统的另一重要组成部分是飞机上装的应答机，应答机是一个在接受到相应的信号后能发出不同形式编码信号的无线电收发机，应答机在接收到地面二次雷达发出的询问信号后，进行相应回答。这些信号被地面的二次雷达天线接收，经过译码，就在一次雷达屏幕出现的显示这架飞机的亮点旁边显示出飞机的识别号码和高度，管制员就会很容易地了解飞机的位置和代号。为了使管制员在询问飞机的初期就能很快地把屏幕上的光点和所对应的飞机联系起来，机上应答机还具有识别功能，驾驶员在管制员要求时可以按下“识别”键，这时应答机发出一个特别位置识别脉冲（SPI），这个脉冲使地面站屏幕上的亮点变宽，以区别于屏幕上的其他亮点。使用二次监视雷达的，可以通过下盯橡列程序识别航空器：

a.观察到雷达显示有航空器按指示调定的应答机编码；

b.观察到仅有一个雷达目标，显示有航空器按指示使用特殊位置识别功能的信息；

c. 对于设定某一指定编码的航空器，观察其遵守指令的情况；

d.通过雷达识别移交。

三、PSR与SSR主要性能比较 在下列情况下，管制员应当向航空器通报其位置：

（一）航空器第一次被识别时；

（二）航空器驾驶员要求提供服务时；

（三）航空器报告的位置与雷达管制员根据雷达观察到的位置有显著差别时；

（四）雷达引导后，如果现行指令使航空器偏离其原规定的航路，指示航空器恢复自主领航时；

（五）结束雷达服务前，如果观察到航空器偏离原规定的航路时。实施雷达管制（可能雷达引导、雷达间隔或者进行速度控制）

通过下列方式进行雷达识别的可不用通报航空器位置：

（一）根据航空器的位置报或起飞识别；

（二）通过使用指定的二次雷达编码或S模式，并且观察到航空器显示位置与其现行飞行计划一致；

（三）通过雷达识别的移交。

提供给航空器的位置情报，可以采用下列方式之一：

（一）相对于一个显著的地理位置；

（二）到一个重要点、航线上一个导航设备的方位和距离；

（三）距一个已知位置的方向和距离；

（二） 作最后进近的航空器距接地点的距离。 为提供不间断的雷达服务，应实施雷达管制移交。只有当认为航空器在接受管制员的雷达覆盖区内时，一名管制员才能试图将雷达识别移交个另一名雷达管制员。进行航空器雷达识别移交的方法如下:

a.两个雷达管制席相邻或者使用同一显示器时，移交方直接在接受方显示器上指出雷达位置指示符的名称；

b.两个雷达显示屏上都标有同一地理位置或导航设备，利用通信设备说明航空器距上述位置的相对方位和距离，必要时，应当指明航空器的航向；

c.利用自动化手段指定雷达位置指示符的名称；

d.当S模式覆盖有效时，将装有S模式航空器识别功能的航空器通知接受方；

f.移交方雷达管制员指示航空器变换编码或用特殊位置识别，接受方雷达管制员予以证实。 当雷达管制员要把某架航空器的管制权移交给非雷达管制员时，实施移交前，要确保该航空器与任何其它受管制航空器之间具有非雷达间隔。雷达管制服务中断或终止时应当立即通知航空器，但在下列情况可不必通报：

a.航空器改为目视飞行；

b.航空器已经着陆,或已经按指令转换到其他频率上；

c.航空器已经结束精密雷达进近。

CS

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%产生Stripmap SAR的回波

clear all

clc

thetaT=0%T平台波束斜视角

thetaT=thetaT*pi/180%rad

thetaR=0%R平台波束斜视角

thetaR=thetaR*pi/180

c=3e8%光速

fc=1.5e9%载频

lambda=c/fc%波长

%%测绘带区域

X0=200%方位向[-X0,X0]

Rtc=3000

Rrc=3000

Rc=(Rtc+Rrc)/2

R0=150%距离向[Rc-R0,Rc+R0]

%%距离向(Range),r/t domain

Tr=1.33e-6%LFM信号脉宽1.33us (200m)

Br=150e6%LFM信号带宽 150MHz

Kr=Br/Tr%调频斜率

Nr=1024

r=Rc+linspace(-R0,R0,Nr)

t=2*r/c%t域序列

dt=R0*4/c/Nr%采样周期

f=linspace(-1/2/dt,1/2/dt,Nr)%f域序列

%%方位向(Azimuth,Cross-Range),x/u domain

v=100%SAR 平台速度

Lsar=300%合成孔径长度

Na=512

x=linspace(-X0,X0,Na)%u域序列

u=x/v

du=2*X0/v/Na

fu=linspace(-1/2/du,1/2/du,Na)%fu域序列

ftdc=v*sin(thetaT)

ftdr=-(v*cos(thetaT))^2/lambda/Rtc

frdc=v*sin(thetaR)

frdr=-(v*cos(thetaR))^2/lambda/Rrc

fdc=ftdc+frdc%Doppler调频中心频率

fdr=ftdr+frdr%Doppler调频斜率

%％目标位置

Ntar=3%目标个数

Ptar=[Rrc,0,1 %距离向坐标,方位向坐标,sigma

Rrc+50,-50,1

Rrc+50,50,1]

%%产生回波

s_ut=zeros(Nr,Na)

U=ones(Nr,1)*u%扩充为矩阵

T=t'*ones(1,Na)

for i=1:1:Ntar

rn=Ptar(i,1)xn=Ptar(i,2)sigma=Ptar(i,3)

rtn=rn+Rtc-Rrc

RT=sqrt(rtn^2+(rtn*tan(thetaT)+xn-v*U).^2)

RR=sqrt(rn^2+(rn*tan(thetaT)+xn-v*U).^2)

R=RT+RR

DT=T-R/c

phase=-pi*Kr*DT.^2-2*pi/lambda*R

s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)<Tr/2).*(abs(v*U-xn)<Lsar/2)

end

%方位向fft

s_kt=fftshift(fft(fftshift(s_ut).')).'

%CS变换

kc=4*pi/lambda

kc=kc*ones(1,Na)

kx=fu/v

p_kx0=-sqrt(kc.^2-kx.^2)%相位项泰勒展开的系数函数

p_kx1=2*kc/c/p_kx0

p_kx2=-2.*kx.^2/早基c^2./p_kx0.^3

C_kx=-(c*p_kx1/2+1)

Ks_r=1-2*Kr*Rc.*p_kx2

Ks_kx_r=Kr/pi./Ks_r

r0=Rc

s2_ut=exp(j*pi*C_kx.*ones(Nr,1)*Ks_kx_r.*(t'*ones(1,Na)-2*r0*(1+C_kx)/c).^2)%设计的线性调频信号

S_cs=s_kt.*s2_ut

%距离向fft

S_kw=fftshift(fft(fftshift(S_cs)))

%距离向匹配滤波

w=2*pi*f

rmc_r=exp(j.*w*2*C_kx*r0/c).*exp(j.*w.^2/4/贺樱pi/Kr/(1+C_kx))

rmc_r=rmc_r'陆拍谨*ones(1,Na)

S_rmc=S_kw.*rmc_r

%距离向ifft

S_kt=fftshift(ifft(fftshift(S_rmc)))

d_kxr=4*pi/c^2*Kr*C_kx*(1+C_kx).*(Rc-r0).^2%CS变换带来的相位误差

S_kt=S_kt.*exp(-j*d_kxr)%消除相位误差

%方位向匹配滤波

FU=ones(Nr,1)*fu

H_kx=exp(j*pi/fdr*(FU-fdc).^2)%方位向压缩因子

I_ut=S_kt.*H_kx

I_ut=fftshift(ifft(fftshift(I_ut.'))).'

subplot(221)

G=20*log10(abs(s_ut)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,-G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(a)原始信号')

subplot(222)

G=20*log10(abs(S_rmc)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,-G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(b)距离向匹配滤波后频谱')

subplot(223)

G=20*log10(abs(S_kt)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(c)消除相位误差后频谱')

subplot(224)

G=20*log10(abs(I_ut)+1e-6)

gm=max(max(G))

gn=gm-60%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(d)目标图象')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

RD

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%产生Stripmap SAR的回波

clear all

thetaT=0%T平台波束斜视角

thetaT=thetaT*pi/180%rad

thetaR=0%R平台波束斜视角

thetaR=thetaR*pi/180

c=3e8%光速

fc=1.5e9%载频

lambda=c/fc%波长

%%测绘带区域

X0=200%方位向[-X0,X0]

Rtc=3000

Rrc=3000

Rc=(Rtc+Rrc)/2

R0=150%距离向[Rc-R0,Rc+R0]

%%距离向(Range),r/t domain

Tr=1.5e-6%LFM信号脉宽 1.5us (200m)

Br=150e6%LFM信号带宽 150MHz

Kr=Br/Tr%调频斜率

Nr=512

r=Rc+linspace(-R0,R0,Nr)

t=2*r/c%t域序列

dt=R0*4/c/Nr%采样周期

f=linspace(-1/2/dt,1/2/dt,Nr)%f域序列

%%方位向(Azimuth,Cross-Range),x/u domain

v=100%SAR 平台速度

Lsar=300%合成孔径长度

Na=1024

x=linspace(-X0,X0,Na)%u域序列

u=x/v

du=2*X0/v/Na

fu=linspace(-1/2/du,1/2/du,Na)%fu域序列

ftdc=v*sin(thetaT)

ftdr=-(v*cos(thetaT))^2/lambda/Rtc

frdc=v*sin(thetaR)

frdr=-(v*cos(thetaR))^2/lambda/Rrc

fdc=ftdc+frdc%Doppler调频中心频率

fdr=ftdr+frdr%Doppler调频斜率

%％目标位置

Ntar=3%目标个数

Ptar=[Rrc,0,1 %距离向坐标,方位向坐标,sigma

Rrc+50,-50,1

Rrc+50,50,1]

%%产生回波

s_ut=zeros(Nr,Na)

U=ones(Nr,1)*u%扩充为矩阵

T=t'*ones(1,Na)

for i=1:1:Ntar

rn=Ptar(i,1)xn=Ptar(i,2)sigma=Ptar(i,3)

rtn=rn+Rtc-Rrc

RT=sqrt(rtn^2+(rtn*tan(thetaT)+xn-v*U).^2)

RR=sqrt(rn^2+(rn*tan(thetaT)+xn-v*U).^2)

R=RT+RR

DT=T-R/c

phase=pi*Kr*DT.^2-2*pi/lambda*R

s_ut=s_ut+sigma*exp(j*phase).*(abs(DT)<Tr/2).*(abs(v*U-xn)<Lsar/2)

end

%%距离压缩

p0_t=exp(j*pi*Kr*(t-2*Rc/c).^2).*(abs(t-2*Rc/c)<Tr/2)%距离向LFM信号

p0_f=fftshift(fft(fftshift(p0_t)))

s_uf=fftshift(fft(fftshift(s_ut)))%距离向FFT

src_uf=s_uf.*(conj(p0_f).'*ones(1,Na))%距离压缩

src_ut=fftshift(ifft(fftshift(src_uf)))%距离压缩后的信号

src_fut=fftshift(fft(fftshift(src_ut).')).'%距离多普勒域

%%二次距离压缩,距离迁移校正原理仿真

src_fuf=fftshift(fft(fftshift(src_uf).')).'%距离压缩后的二维频谱

F=f'*ones(1,Na)%扩充为矩阵

FU=ones(Nr,1)*fu

p0_2f=exp(j*pi/fc^2/fdr*(FU.*F).^2+j*pi*fdc^2/fc/fdr*F-j*pi/fc/fdr*FU.^2.*F)

s2rc_fuf=src_fuf.*p0_2f

s2rc_fut=fftshift(ifft(fftshift(s2rc_fuf)))%距离多普勒域

%%方位压缩

p0_2fu=exp(j*pi/fdr*(FU-fdc).^2)%方位向压缩因子

s2rcac_fut=s2rc_fut.*p0_2fu%方位压缩

s2rcac_fuf=fftshift(fft(fftshift(s2rcac_fut)))%距离方位压缩后的二维频谱

s2rcac_ut=fftshift(ifft(fftshift(s2rcac_fut).')).'%方位向IFFT

subplot(221)

G=20*log10(abs(s_ut)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,-G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(a)原始信号')

subplot(222)

G=20*log10(abs(src_fut)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(fu,r-Rc,-G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(b)距离多普勒域频谱')

subplot(223)

G=20*log10(abs(s2rc_fut)+1e-6)

gm=max(max(G))

gn=gm-40%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(fu,r-Rc,-G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(c)RMC后的RD域频谱')

subplot(224)

G=20*log10(abs(s2rcac_ut)+1e-6)

gm=max(max(G))

gn=gm-60%显示动态范围40dB

G=255/(gm-gn)*(G-gn).*(G>gn)

imagesc(x,r-Rc,G),colormap(gray)

grid on,axis tight,

xlabel('Azimuth')

ylabel('Range')

title('(d)目标图象')

1、将雷达装载在车辆新装在车赋正确的位置，如贴在车尾后方或者装载在车族燃顶；

2、组装各探头：将位于雷达外壳上游槐的6个探头，按照要求依次装载到车辆新装在车位，确保各探头能够正确发出探测信号；

3、调节探头角度：根据要求，将各探头的角度定位。

4、插上探头和雷达的接口线：将各探头的接口线插入到雷达的接口上，按要求紧固，避免松脱；

5、上电调试：将雷达与电源相连，然后上电；根据雷达的测试程序，对雷达进行调试，神穗友确保探头正常工作。

6、验证雷达：用有物体模拟目标进行测试，观察雷达的检测结果，验证雷达的正常性。

---