雷达对抗原理.docx
《雷达对抗原理.docx》由会员分享，可在线阅读，更多相关《雷达对抗原理.docx（14页珍藏版）》请在冰点文库上搜索。
雷达对抗原理
雷达对抗原理大作业
学院：
电子工程学院
班级：
021215
姓名：
zkw
学号：
一、单脉冲和差测角仿真实验
单脉冲测角是雷达中常用的一种测角方法,它利用多个天线同时接收回波信号,通过比较回波信号的幅度或相位来获得目标的角位置信息。
确定目标准确空间位置是雷达的重要任务。
在大多数场合,雷达不仅要测定目标的距离,还需要测定目标的角度,即测定目标的角坐标,其中包括目标的方位角和高低角（仰角）。
雷达测角的性能包含测角范围、测角速度、测角精度和角分辨率等,准确度与测角误差是其中极为关键的指标。
1. 单脉冲和差测角原理
雷达测角的物理基础是电波在均匀介质中传播的直线性和雷达天线的方向性，分为振幅法和相位法两大类，其中振幅法测角又分为最大信号法和等信号法，等信号测角采用两个相同且彼此部分重叠的波束，其方向图如下图1所示，若目标处在两波束的交叠轴OA方向，则两波束收到的信号强度相等，否则一个波束收到的信号强度高于另一个，故常称OS为等信号轴。
当两个波束收到的回波信号相等时，等信号轴所指的方向即为目标方向。
若目标处在OB方向，波束2的回波比波束1的强，处在OC方向时，则与之相反，因此比较两个波束回波的强弱就可以判断目标偏离等信号轴的方向，并可用查表的方法估计出偏离等信号轴的大小。
2、天线电压方向性函数
雷达测角的基础是电波在均匀介质中传播的直线性和雷达天线的方向性。
天线的方向性可用其方向性函数或根据方向性函数画出的方向图表示。
但方向图的准确表达式往往很复杂，因而常采用简单函数来近似，常用的有余弦函数、高斯函数、辛克函数。
方向图的主要技术指标是半功率波束宽度θ以及副瓣电平。
3、设天线电压方向性函数为F（θ），等信号轴OA的指向为0，则波束1、2的方向性函数可分别写为
θk为与波束最大值方向的倾角
用等信号法测量时，波束1和波束2收到的回波信号为
θt为目标偏离等信号轴的角度
对信号进行和差法处理，可获得目标信号的差值及和值，即
在等信号轴附近差信号及和信号可近似表示为
归一化和差值为
Δ/∑由于正比于目标偏离θ0的角度θt，故可用它来判读的大小及方向。
4.单脉冲和差测角仿真
Matlab程序如下：
k=0.730;
d=0.160;
labda=2*pi/k;
theta_3db=1.2*labda/d;
theta_k=theta_3db/3;
theta=-2*theta_3db:
0.2:
2*theta_3db;
f1=exp（-1.3863*（theta-theta_k）.^2/theta_3db^2）;f2=exp（-1.3863*（theta+theta_k）.^2/theta_3db^2）;sigma=f1+f2;delta=f1-f2;
figure,subplot（221）,plot（theta,f1,'r-'）,gridon
holdon,plot（theta,f2）,xlabel（'角度\theta'）,ylabel（'两个响应'）;
subplot（222）,plot（theta,sigma）,xlabel（'角度\theta'）,ylabel（'和波束\Sigma'）,gridon
subplot（223）,plot（theta,delta）,xlabel（'角度\theta'）,ylabel（'差波束\Delta'）,gridon
subplot（224）,plot（theta,（delta./sigma））,gridon
xlabel（'角度\theta'）,ylabel（'\Delta/\Sigma'）
S
3.数据制表Δ\∑
由于最大单值测角范围为有限，因此只考虑在[-30°,30°]范围内的数据，本文中为简化内容，只讨论[-15°，15°]范围内的数值。
利用一次回归曲线拟合，得到对应的一次曲线。
k=0.730;
d=0.160;
labda=2*pi/k;
theta_3db=1.2*labda/d;
theta_k=theta_3db/3;
theta=-15:
1:
15;
f1=exp（-1.3863*（theta-theta_k）.^2/theta_3db^2）;
f2=exp（-1.3863*（theta+theta_k）.^2/theta_3db^2）;
sigma=f1+f2;
delta=f1-f2;
t=[-15:
1:
15]';
f=[-0.2115-0.1978-0.1840-0.1701-0.1562-0.1422-0.1281-0.1140-0.0999-0.0857-0.0715-0.0572-0.0429-0.0286-0.014300.01430.02860.04290.05720.07150.08570.09990.11400.12810.14220.15620.17010.18400.19780.2115]';
x=[ones（size（t））t];
a=x\f
f3=[ones（size（t））t]*a;
figure,plot（theta,（delta./sigma）,'r-'）
xlabel（'½Ç¶È\theta'）,ylabel（'\Delta/\Sigma'）;
gridon;
holdon;
plot（t,f3,'*'）
程序运行结果如下：
a=
0.0000
0.0142
得到直线方程为Δ\∑=0.0142θ
图中原始图像为红色直线，拟合结果为星号连点。
从程序运行结果可以看出，利用二次回归曲线拟合得到的星号直线和原始数据绘
制出来的直线基本上重合，可以看出，结果的误差是很小的。
θ
-20
-19
-18
-17
-16
-15
-14
-13
-12
Δ/∑
-0.4318
-0.4128
-0.3935
-0.3738
-0.3537
-0.3333
-0.3126
-0.2916
-0.2704
-11
-10
-9
-8
-7
-6
-5
-4
-3
-2
-0.2488
-0.227
-0.205
-0.1828
-0.1603
-0.1377
-0.115
-0.0922
-0.0692
-0.0462
-1
0
1
2
3
4
5
6
7
8
-0.0231
0
0.0231
0.0462
0.0692
0.0922
0.115
0.1377
0.1603
0.1828
9
10
11
12
13
14
15
16
17
18
0.205
0.227
0.2488
0.2704
0.2916
0.3126
0.3333
0.3537
0.3738
0.3935
19
20
0.4128
0.4318
由上表也可以看出计算得到的和曲线拟合出来的结果误差很小。
二、仿真线性调频信号
1.LFM脉冲的匹配滤波
脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。
这种体制采用宽脉冲发射以提高发射的平均功率，保证足够大的作用距离；而接受时采用相应的脉冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频信号,接收时采用匹配滤波器压缩脉冲。
从雷达回波信号
提取出表征目标特性的
（表征相对距离）和
（表征目标反射特性），常用的方法是让
通过雷达发射信号
的匹配滤波器，如图2.1所示
图2.1雷达回波信号匹配滤波处理
的匹配滤波器
为：
于是，
如图3.1,
经过系统
得输出信号
，
当
时,
（2.1）
当
时,
（2.2）
合并2.1和2.2两式：
（2.3）
2.3式即为LFM脉冲信号经匹配滤波器得输出,它是一固定载频
的信号。
当
时，包络近似为辛克（sinc）函数。
（2.4）
图2.2匹配滤波的输出信号
如图3.2，当
时，
为其第一零点坐标；当
时，
，习惯上，将此时的脉冲宽度定义为压缩脉冲宽度。
LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度
之比通常称为压缩比D，
上式表明，压缩比也就是LFM信号的时宽频宽积。
2.线性调频信号回波matlab仿真
closeall
clearall
eps=0.000001;
B=200.0e6;
T=10.e-6;
mu=B/T;
Fs=2*B;Ts=1/Fs;
N=T/Ts;
t=linspace（-T/2,T/2,N）;
fi=10e6;
St=exp（j*pi*mu*t.^2）;
Sr=exp（j*pi*mu*（t-1e-7）.^2）;
Ht=exp（-j*pi*mu*t.^2）;
Jt=exp（j*2*pi*（mu.*t.^2/2+fi.*t））;
S_J=St+Jt;
Sot=conv（St,Ht）;
Sot1=conv（Sr,Ht）;
Sot2=conv（S_J,Ht）;
L=2*N-1;
t1=linspace（-T,T,L）;
Z=abs（Sot）;Z=Z/max（Z）;
Z1=abs（Sot1）;Z1=Z1/max（Z1）;
Z2=abs（Sot2）;Z2=Z2/max（Z2）;
freq=linspace（-Fs/2,Fs/2,N）;
LFMFFT=fftshift（abs（fft（St）））;
figure
（1）,subplot（211）,plot（t*1e6,real（St））,grid;
xlabel（'时间/us'）,title（'实部'）,axis（[-11-11]）;
subplot（212）,plot（t*1e6,real（Sr）,'r'）,grid
axis（[-11-11]）,xlabel（'时间/us'）,title（'实部'）
figure
（2）,subplot（211）,plot（t*1e6,imag（St））,grid
axis（[-11-11]）,xlabel（'时间/us'）,title（'虚部'）
subplot（212）,plot（t*1e6,imag（Sr）,'r'）,grid
axis（[-11-11]）,xlabel（'时间/us'）,title（'虚部'）
figure（3）,plot（freq*1e-6,LFMFFT）,xlabel（'频率/MHz'）
grid,title（'LFM信号的频谱'）;
figure（4）,plot（t1*B,Z）,title（'LFM匹配滤波后的波形'）;
hold,plot（t1*B,Z1）,axis（[-305001]）,grid;
figure（5）,plot（t*1e6,real（S_J））,grid;
xlabel（'时间/us'）,ylabel（'幅度'）;
axis（[-1.51.5-1.51.5]）,title（'有杂波干扰的回波信号'）
figure（6）,plot（t1*B,Z2）,axis（[-20010001]）,grid;
title（'有杂波干扰的回波信号匹配滤波输出'）;
xlabel（'时间/us'）,ylabel（'幅度'）;
升级会员 