雷达信号处理基本流程

基本雷达信号处理流程
一、脉冲压缩
窄带(或某些中等带宽)的匹配滤波: 相关处理，用FFT数字化执行，即快速卷积处理，可以在基带实现( 脉冲压缩) 快速卷积，频域的匹配滤波
脉宽越小，带宽越宽，距离分辨率越高；
脉宽越大，带宽越窄，雷达能量越小，探测距离越近；
D=BT (时宽带宽积)；
脉压流程：
频域：回波谱和参考函数共轭相乘
时域：相关
即输入信号的FFT乘上参考信号FFT的共轭再逆FFT;
Sc=ifft(fft(Sb).*conj(fft(S)));
匹
配
滤
波
器
Taskl
f0=10e9;% 载频tp=10e-6;% 脉冲宽度 B=10e6;% 信号带宽fs=100e6;% 采样率
R0=3000;% 目标初始距离 N=4096;c=3e8;tau=2*R0/c;beita=B/tp;t=(0:N-1)/fs;
Sb=rectpuls(t-tp/2-tau,tp).*exp(j*pi*beita*(t-tp/2-tau).A2).*exp(-2j*pi*fO*tau);% 回波信号
So=ifft(fft(Sb).*conj(fft(S)));%
脉压
S=rectpuls(t-tp/2,tp).*exp(i*pi*beita*(t-tp/2).A
2);
一
■
____
3
1
120
x 10
figure(7);
plot(t*c/2, db(abs(So)/max(So)) )% 归一化 dB
grid on
二、去斜处理（宽带的匹配滤波）
去斜处理“有源相关”通常用来处理极大带宽的LFM波形（如果直接采样的话因
为频带很宽所以在高频的时候需要的采样率就很大，采样点数就很多，所以要经过去斜处理）
Stretch方法是针对线性调频信号而提出的，其方法是将输入信号与参考信号（经适当延迟的本振信号，延迟量通常由窄带信号测距结果估计出）混频，则每一个散射点就对应一个混频后的单频分量，对混频输出的信号进行DFT处理，即可获得目标的距离像，对参考信号的要求是应具有与输入信号相同的调频斜率。

去斜处理流程：
一►风_J频谱分析|_►输出信号
输入信号
混频过程为回波信号在时域与参考信号的共轭相乘
混频后得到一个瞬时频率和目标距离成正比的单频信号，对其进行频谱分析即可得到目标的距离像；
去斜处理一般情况下可降低信号带宽；
%%%%%%%%%%%%%%%%%%%%%%%% 去斜处理仿真程
序 %%%%%%%%%%%%%%%%%%%%%%%%%
clc;clear all;close all;
B=10e6;% 带宽 10MHz
tp=10e-6;% 脉宽 10us
k=B/tp;%LFM 系数
fs=50e6;
R0=3e3;R1=2000;R2=3500;R=5000;
c=3e8;
f0=60e6;
N=rou nd(2*R/c*fs);
fft_N=2A nextpow2(N);
t=li nspace(0,2*R/c,N);
%%%%%%%%%%%%%%%%%%%%%%%%%% 参考信号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sref=exp (2i*pi*fO*t).*exp(1i*pi*k*t.A2);
%%%%%%%%%%%%%%%%%%%%%%%%%%% 回波信
号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sb0=exp(1j*pi*k*(t-2*R0/c).A2).*exp (2j *pi*f0*(t-2*R0/c));
Sb1=exp(1j*pi*k*(t-2*R1/c).A2).*exp (2j *pi*f0*(t-2*R1/c));
Sb2=exp(1j*pi*k*(t-2*R2/c).A 2).*exp (2j *pi*f0*(t-2*R2/c));
Sb=Sb0+Sb1+Sb2;
0 -1
0.5
度
幅
参考信号实部
1
1.5
2
2.5
仿真时间It 参考信号虚部
3
3.5
-5
X 10
1 1.5
2
2.5 3 度0 幅
-1 0
0.5
仿真时间/t 3.5 -5
x 10
100 度 50
幅
0 参考信号频谱
-1
1 1.5 -2.5
-2 -1.5
-0.5
0 0.5 频率f/Hz
2 2.5 7
x 10
%%%%%%%%%%%%%%%%%%%%%%%%%%% 混频信
号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%
SSb=Sref.*conj(Sb);% 去斜后时域信号
spectrum=fft(SSb,fft_N);% 去斜后频域信号
f=fs*(0:fft_N-1)/fft_N-fs/2;% 从-fs/2 至U fs/2
f=f*c*tp/2/B;%瞬时频率对应的距离
sf=exp(-j*pi/k*f.A 2);% 滤波器传输函数
SSb=spectrum.*sf;% 从频域去距离扭曲，实现了压缩和去 RVP figure;
-1
回波信号1实部
度
0 0.5
1.5
2 2.5 3
3.5
仿真时间t/S -5
x 10
度
0 回波信号2实部
0.5 1
1.5 2
2.5 仿真时间t/s 回波信号3实部
3
3.5
-1 x 10
-5
仿真时间t/s
-5
x 10
3
1
-1
-3
三个 目 标组成 的 回 波信号 的实 部
2
-2
0.5
1
1.5
时 间
2
t/s
2.5
3
3.5
-5
x 10
SSb=fftshift(SSb);
SSb 仁ifft(SSb)% 消除了距离扭曲和 RVP 的时域信号
subplot(211);
plot(f,db(abs(SSb)/max(SSb)))
xlabel('距离 /m');
grid on subplot(212); plot(f,abs(SSb)) xlabel('距离 /m');
grid on
二、加窗
信号的截取产生了能量泄漏，而用FFT 算法计算频谱又产生了栅栏效应，在FFT 分析中为了减少或消除频谱能量泄漏及栅栏效应可采用不同的截取函数对信号
jl j
丁乩 AF-H..丹Ki
Q 八
% J 3 i'
1
■■
"I 1
-20 -40 -60
-80
进行截短，截短函数称为窗函数，简称窗。

理 %%%%%%%%%%%%%%%%%%%%%%%%
clc;clear all;close all;
c=3e8;% 光速
tau=2*R0/c;
N=rou nd(2*R/c*fs);
fft_N=2A nextpow2(N);
t=(0:fft_N-1)/fs;
s=rectpuls(t-tp/2,tp).*exp(j*pi*k*(t-tp/2).A2); %% 发射信号
spectrum_s=fft(s,fft_N);% 参考信号频谱
spectrum_s=fftshift(spectrum_s);
sb=rectpuls(t-tp/2-tau,tp).*exp(j*pi*k*(t-tp/2-tau).A2).*exp(-2j*pi*fO*tau); %% 回波信号
%%时域加窗
sm=hammi ng(rou nd( tp*fs))'.*s(1:rou nd(tp*fs));% 参考信号加窗
%%频域加窗
%%%%%%%%%%%%%%%%%%%%%%%%
窄带加窗处
f0=10e9% 载频 B=10e6;%信号带宽 tp=10e-6;%脉冲宽度
fs=100e6;%采样频率
k=B/tp;%LFM 系数，线性调频率 R0=3000%初始距离
R=6000;
%找频谱的-4dB压缩点，窗函数严格与该压缩点之间的频谱对应
hammi ng仁[zeros(1855,1) ',hammi ng(387)', zeros(1854,1)'];
spectrum_sm=ham min g1.*spectrum_s;
%%脉压
fft_ssb=ifft(fft(sb).*conj(fft(s)));
fft_smsb=ifft(fft(sb).*conj(fft(sm,le ngth(sb))));
fft_spsb=ifft(fft(sb).*con j(fftshift(spectrum_sm)));
f/Hz
x 10
口湫凰点 %：009l7=y 凰点刖険 %：000£=0d neudva^=<sj%!ge09=sj 滋雲 IAId1%!dva=n snok 粵曹 %!g-eoL=d; zniAIOL 粵卑 %：9eOL=a
%%%%%%%%%%%%%%%%%%%%%%%%%%%% ®
咖皿魅¥%%%%%%%%%%%%%%%%%%%%%%%%%%%
oooz 0009 0009 ooot
^
000£000乙000 L 0
09£-
0017-
00£-
093-
09 k
00 k
09-
003-型
ef
c=3e8; f0=60e6;% 载频
N=rou nd(2*R/c*fs); fft_N=2A nextpow2(N);
t=li nspace(0,2*R/c,N);
f=fs*(0:fft_N-1)/fft_N-fs/2;% 从-fs/2 至U fs/2
%%%%%%%%%%%%%%%%%%%%%%%%%% 参考信
号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sref=exp(1i*pi*u*t.A2); %%%%%%%%%%%%%%%%%%%%%%%%%%%% 回波信号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sb=rectpuls(t-2*R0/c,tp).*exp(1j*pi*u*(t-2*R0/c).A2); %%%%%%%%%%%%%%%%%%%%%%%%%%%% 混频信号 %%%%%%%%%%%%%%%%%%%%%%%%%%%% ssb=Sref.*conj(Sb);
%%加窗
% w = hammi ng(502)';
% hammi ng=[zeros(749,1)', w-mi n( w),zeros(249,1)'];
% hammi ng=abs(hammi ng)/max(hammi ng);
hammi ng=[zeros(749,1) ',hammi ng(502) ：zeros(249,1)'];
ssb0=ham min g.*ssb;
spectrum_ssb0=fft(ssb0,fft_N); % 一维距离像
spectrum_ssb=fft(ssb,fft_N);
f=f*c*tp/2/B;%瞬时频率对应的距离
%%图6
plot(f,db(abs(fftshift(spectrum_ssb))/max(fftshift(spectrum_ssb))))
hold on
plot(f,db(abs(fftshift(spectrum_ssbO))/max(fftshift(spectrum_ssbO))),'r') hold off
去斜加窗与否 对比图
figure;
消除了 RVP
和距离扭曲的混频信号
400 度 幅
200
-4000 -3000 -2000 -1000
0 1000 2000 3000 4000
距离/m
二、检测
1、脉冲多普勒（PD 处理）
多普勒效应：fd=2v/c*f0 , v 为镜像速度； 慢时间维上的采样点做FFT 可以测出目标的速度; 使用复信号：频率正负可测量目标速度的方向；
clc;clear all;close all;
f0=10e9% 载频
tp=10e-6;%脉冲宽度
B=10e6;% 带宽
fs=100e6;%采样频率
R0=3000;%初始距离
消除了 RVP 和距离扭曲的混频信号
归
-3000 1000 2000 3000 度 幅 化 -150 -4000
-1000 距离/m
J
-'1 1
>
7";询
—1
-2000
4000
c=3e8;% 光速
R=4500;%距离波门
gate=R+tp*c/2;% 距离波门加脉宽对应距离
N=rou nd(2*gate/c*fs);% 波门内采样点个数
fft_N=2A nextpow2(N);
t=0:1/fs:tp;% 信号长度
echo_t=li nspace(0,2*gate/c,N);% 波门长度
tau=2*R0/c;
k=B/tp;%调频系数
Tr=100e-6;%脉冲重复周期
CPI=64;%总脉冲个数
v=60;%目标速度，朝向雷达
%发射信号
s=exp(i*pi*k*t.A2);
%回波信号
for m=1:CPI
sb(m,:)=rectpuls((echo_t-2*(R0-(m-1)*v*Tr)/c-tp/2)/(tp)).*exp(1i*pi*k*(echo_t-2*(R0-(m-1)*v*Tr) /c).A2-1i*pi*2*f0*rou nd(2*R0/c*fs)+1i*2*pi*(2*f0*v/c)*(m-1)*Tr)+sqrt(0.1)*(ra ndn (1,N)+1i*ra ndn( 1,N));
end
%脉压 fft_n=2An extpow2(le ngth(t)+N-1); fft_s=fft(s,fft_ n);
for m=1:1:CPI
fft_sb(m,:)=fft(sb(m,:),fft_ n); fft_ssb(m,:)=ifft(fft_sb(m,:).*conj(fft_s));
z(m,:)=abs(fft_ssb(m,(1:N))); z1(m,:)=z(m,:)/max(z(m,:));
z1(m,:)=20*log10(z1(m,:));
[maxval,maxpo]=max(z1(m,:));
end
%FFT
for fm=1:N
dop(:,fm)=fft(fft_ssb(:,fm));
a_dop(:,fm)=fftshift(abs(dop(:,fm))); end
%求极大值对应的坐标
[maxva,max_v]=max(a_dop(:,maxpo));
%PD测速
fd=(max_v-33)/CPI/Tr;
v_pd=fd*c/2/f0 %测速范围
fd_max=1/Tr/2;
v_max=fd_max*c/2/f0
%测速精度
det_fd=1/Tr/64;
det_v=det_fd*c/2/f0
figure;
mesh（echo_t*c/2,li nspace（-75,75,64）,a_dop）; axis tight;
xlabel（'距离：m'）;
ylabel（'速度：m/s'）;
title（'二维距离-多普勒平面'）;
二维距离-多普勒平面
斗
T m
v_pd =60.9375 v_max =75 det_v =2.3438
2、形心法测距测速
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 形心
法 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc;clear all;close all
f0=10e9;% 载频
tp=10e-6;%脉冲宽度
B=10e6;% 带宽
fs=100e6;%采样频率
R0=3000%初始距离 c=3e8;% 光速
N=4096; %此为培训期间数据，实际情况可以根据波门和信号宽度求岀N
t=(0:N-1)/fs;
sn r=20;
tau=2*R0/c;
k=B/tp;%调频系数
Tr=100e-6;%脉冲重复周期
CPI=64;%总脉冲个数
v=60;%目标速度，朝向雷达
sigmaf=1A2/(10A(s nr/10));
s=rectpuls(t-tp/2,tp).*exp(j*pi*k*(t-tp/2).A2);% 发射信号% figure;
% plot(t,real(s))
% xlabel('时间 /s');
% ylabel('幅度');
% title('发射信号实部');
% grid on for m=1:CPI
taum=2*(R0-m*Tr*v)/c;
sb=rectpuls(t-taum-tp/2).*exp(j*pi*k*(t-taum-tp/2).A2).*exp(-2j*pi*fO*taum);%
回波信号sb_n oise=sb+sqrt(sigmaf/2)*(ra ndn (1,N)+1i*ra ndn (1,N));% 加噪声的回波信号
fft_ssb=ifft(fft(sb).*conj(fft(s)));% 脉压处理
fft_ssb_s nr=ifft(fft(sb _n oise).*c on j(fft(s)));
Group(m,1:N)=fft_ssb;
Group_s nr(m,1:N)=fft_ssb_s nr;
end
figure;
imagesc(t*c/2,1:CPI,abs(fft_ssb))
figure;
imagesc(t*c/2,1:CPI,abs(fft_ssb_s nr))
for n=1:N
Group2=fft(Group(1:CPI,n));% 纵向做 FFT
Group2_2=fftshift(abs(Group2));
Group3( n,1:CPI)=Group2_2;
end
for n1=1:N
Group2_snr=fft(Group_snr(1:CPI,n1));% 纵向做 FFT
Group2_2_s nr=fftshift(abs(Group2_s nr));
Group3_s nr(n 1,1:CPI)=Group2_2_s nr;
end
figure;
mesh(abs(Group3))
figure;
mesh(abs(Group3_s nr))
[li ne,row]=fi nd(abs(Group3)==max(max(abs(Group3))));
[li ne_sn r,row_s nr]=fi nd(abs(Group3_s nr)==max(max(abs(Group3_s nr)))); Ran ge=t*c/2;
PRF=1/Tr;
fd=(-CPI/2:CPI/2-1)*PRF/CPI;
v=fd*c/2/f0;
for Ra=li ne-3:li ne+3
amp=abs(Group3(Ra,row));
C(Ra)=amp*Ra nge(Ra);
D(Ra)=sum(amp);
end
sum(C)/sum(D)
for V=row-3:row+3
in dex=abs(Group3(li ne,V));
E(V)=i ndex*v(V);
F(V)=sum(i ndex);
end
sum(E)/sum(F)
for Ra_sn r=li ne_sn r-3:l ine_sn r+3
amp_s nr=abs(Group3_s nr(Ra_s nr,row_s nr));
C_s nr(Ra_s nr)=amp_s nr*Ra nge(Ra_s nr);
D_s nr(Ra_s nr)=sum(amp_s nr);
end
sum(C_s nr)/sum(D_s nr)
for V_sn r=row_s nr-3:row_s nr+3
in dex_s nr=abs(Group3_s nr(l ine_sn r,V_s nr));
E_s nr(V_s nr)=i ndex_s nr*v(V_s nr);
F_sn r(V_s nr)=sum(i ndex_s nr);
end
sum(E_s nr)/sum(F_s nr)
结果：
ans =3.0000e+003
ans =60.3560
ans =3.0000e+003
ans =60.3533
3、信号检测
结果：
mean_ no ise =0.0010 + 0.0079i
var_ noise = 1.0123
pf1 =1.0000e-003
高斯白噪声取模值后的波形以及检测门限
点数
检测概率相对于SNR曲线987654321 率概测检
0 2 4 6 8 10 12 14 16 18 20
SNR/dB
4、单脉冲测角仿真
单脉冲跟踪雷达是通过比较来自两个或多个同时波束的信号获得目标角位置信息的一种雷达；
目前常用的单脉冲测角方法主要有幅度和差单脉冲测角和相位和差单脉冲测角。

幅度和差单脉冲测角通过比较两个相位中心重合但指向不同的波束得到目标角度信息；相位和差单脉冲测角则通过比较两个相位中心有一定距离但波束指向相互平行的波束得到目标角度信息。

相位和差单脉冲与幅度和差单脉冲的相似之处在于：目标角度坐标都是由一个和通道和两个差通道来提取的。

主要不同之处在是，幅度和差单脉冲产生的四个信号具有相同的相位但具有不同的幅度，而相位和差单脉冲信号具有相同的幅度但有不同的相位。

相位和差单脉冲对每个坐标系（方位和俯仰坐标）使用最少由两个阵元组成的阵列天线。

相位误差信号是由于不同天线阵元产生的信号之间的相位差来计算得出的。

%%%%%%%%%%%%%%%%%%%%%%%% 单脉冲测角仿
真 %%%%%%%%%%%%%%%%%%%%%%%%%%%%
f0=10e6;
R0=3e3;%目标距0号阵元的距离
d=10;%阵列接收天线之间的距离
theta0=0.2*pi/180;% 目标角度
%R0=R0+d*si n( theta0)/2;
c=3e8;% 光速
Iamda=c/f0;
tau0=2*R0/c;%到0号阵元的时延
theta=li nspace(-1*pi/180,1*pi/180,1000);
thetaP=0.15*pi/180;% 偏置角
N=4;%天线个数
%%幅度和差单脉冲测角 %相同相位不同幅度
%波束形成结果
Y=exp (2j *pi*fO*tauO)*exp(j*pi*(N-1)*d*si n( thetaO)/lamda).*(si n( N*pi*d*(si n( thetaO)-si
n( theta))/l amda)./si n( pi*d*(si n(thetaO)-si n(theta))/lamda));
thetaA=theta+thetaP;
thetaB=theta-thetaP;
%偏置波束A、B
Y_thetaA=exp (2j*pi*f0*tau0).*exp(j*pi*(N-1)*d*si n( thetaO)/lamda).*(si n( N*pi*d*(si n( theta0)-si n( thetaA))/lamda)./si n(pi*d*(s in (thetaO)-si n( thetaA))/lamda));
Y_thetaB=exp (2j*pi*f0*tau0).*exp(j*pi*(N-1)*d*s in (thetaO)/lamda).*(si n( N*pi*d*(si n( theta0)-si n( thetaB))/lamda)./si n(pi*d*(si n( theta0)-si n( thetaB))/lamda));
%差波束
Y_delta=Y_thetaA-Y_thetaB;
%和波束
Y_sigma=Y;
%复比
Y_AB=Y_delta./Y_sigma;
thetaAB=li nspace(0*pi/180,0.35*pi/180,1000);
波束形成配相结果
幅度和差单脉冲测角之和波束
度幅
-3 幅度和差单脉冲测角之差波束
5
5
a
a
-
-
1
.
5
5
°
0.1 0.15 0.2 0.25
角度/。

幅度和差单脉冲测角之和差比
0.3 0.35
6
4
2O
益增
压电
-
2
0.1 0.15 0.2
角度/。

0.25 0.3 0.35
%%相位和差单脉冲测角 %相同幅度不同相位
%第一种配相方法
%波束1
beam1=exp (2j *pi*fO*tauO)*[exp (2j *pi*O*d*(si n(thetaO)-si n(theta))/lamda)+exp (2j*pi*1*d*(si n(th etaO)-si n(theta))/lamda)];
%波束2
beam2=exp (2j*pi*f0*tau0)*exp (2j *pi*2*d*si n(thetaO)/lamda)*[exp (2j *pi*0*d*(si n(thetaO)-si n(the
ta))/lamda)+exp (2j *pi*1*d*(si n( theta0)-si n(theta))/lamda)];
beam_sigma=beam1+beam2;% 和口波束
beam_delta=beam1-beam2;% 差波束
beam_12=(beam1-beam2)./(beam1+beam2);% 和差比
3.9999
第一种配相法之和波束
角度/
第一种配相法之差波束
%%第二种配相方法
%波束1
beam1_2=exp (2j *pi*fO*tauO).*[exp (2j*pi*0*d*(si n(thetaO)-si n(theta))/lamda)+exp (:2j*pi*1*d*(s in( thetaO)-si n(theta))/lamda)];
%波束2
beam2_2=exp (2j *pi*fO*tauO).*exp (2j*pi*2*d*(si n(thetaO)-si
n(theta))/lamda).*[exp (2j*pi*0*d*(si n
(thetaO)-si n( theta))/lamda)+exp (2j*pi*1*d*(s in (thetaO)-s in (theta))/lamda)]; beam2_sigma=beam1_2+beam2_2;% 和波束
beam2_delta=beam1_2-beam2_2;% 差波束 beam12_2=(beam1_2-beam2_2)./(beam1_2+beam2_2);%
和差比
第二种配相法之和波束
角度/。

第二种配相法之差波束
角度/
第二种配相方法之和差比
角度/
Welcome To
Download !! !
欢迎您的下载，资料仅供参考！。