雷达信处理基本流程Word格式.docx
《雷达信处理基本流程Word格式.docx》由会员分享，可在线阅读，更多相关《雷达信处理基本流程Word格式.docx（23页珍藏版）》请在冰豆网上搜索。
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*f0*tau）;
%回波信号
120
l«
B:
-
•
;
——J
80
60
40
1
2
3
4
5
6
7
8
9
10
S=rectpuls（t-tp/2,tp）.*exp（i*pi*beita*（t-tp/2）.A2）;
%发射信号（参考信号）
脉压
J
丄
I
k
So=ifft（fft（Sb）.*conj（fft（S）））;
%figure（7）;
plot（t*c/2,
db（abs（So）/max（So））
）%归一化dB
gridon
-50
-100
1川X
iT
2」
f
I1r
1—m
-150
-200
-250
-300
-350
1000
2000
3000
4000
5000
6000
7000
二、去斜处理（宽带的匹配滤波）
去斜处理“有源相关”，通常用来处理极大带宽的LFM波形（如果直接采样的话因为频带很宽所以在高频的时候需要的采样率就很大，采样点数就很多，所以要经过去斜处理）
Stretch方法是针对线性调频信号而提出的，其方法是将输入信号与参考信号（经适当延迟的本振信号，延迟量通常由窄带信号测距结果估计出）混频，则每一个散射点就对应一个混频后的单频分量，对混频输出的信号进行DFT处理，即可获得目标的距离像，对参考信号的要求是应具有与输入信号相同的调频斜率。
去斜处理流程：
输入信号一►述~』频谱分析|一►输出信号
混频过程为回波信号在时域与参考信号的共轭相乘
混频后得到一个瞬时频率和目标距离成正比的单频信号，对其进行频谱分析即可得到目标的距离像；
去斜处理一般情况下可降低信号带宽；
%%%%%%%%%%%%%%%%%%%%去斜处理仿真程序%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
clearall;
closeall;
B=10e6;
%带宽10MHz
tp=10e-6;
%脉宽10us
k=B/tp;
%LFM系数
fs=50e6;
R0=3e3;
R1=2000;
R2=3500;
R=5000;
f0=60e6;
N=round（2*R/c*fs）;
fft_N=2Anextpow2（N）;
t=linspace（0,2*R/c,N）;
%%%%%%%%%%%%%%%%%%%%%%＜参考信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sref=exp（2i*pi*fO*t）.*exp（1i*pi*k*t.A2）;
X10
%%%%%%%%%%%%%%%%%%%%%%%%波信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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）.A2）.*exp（2j*pi*f0*（t-2*R2/c））;
Sb=Sb0+Sb1+Sb2;
f=fs*（0:
fft_N-1）/fft_N-fs/2;
%从-fs/2至Ufs/2
f=f*c*tp/2/B;
%瞬时频率对应的距离
sf=exp（-j*pi/k*f.A2）;
%滤波器传输函数
SSb=spectrum.*sf;
%从频域去距离扭曲，实现了压缩和去RVP
figure;
SSb=fftshift（SSb）;
SSb1=ifft（SSb）;
%消除了距离扭曲和RVP的时域信号subplot（211）;
plot（f,db（abs（SSb）/max（SSb）））xlabel（'
距离/m'
）;
gridonsubplot（212）;
plot（f,abs（SSb））xlabel（'
二、加窗
信号的截取产生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，在FFT分析中为了减少或消除频谱能量泄漏及栅栏效应可采用不同的截取函数对信号进行截短，截短函数称为窗函数，简称窗。
%%%%%%%%%%%%%%%%%%%%窄带加窗处理%%%%%%%%%%%%%%%%%%%%%%%%clc;
closeall
f0=10e9;
%载频
fs=100e6;
%采样频率
%光速
fft_N=2Anextpow2（N）;
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*f
0*tau）;
%%回波信号
%%时域加窗
sm=hamming（round（tp*fs））'
.*s（1:
round（tp*fs））;
%参考信号加窗
%%频域加窗
%找频谱的-4dB压缩点，窗函数严格与该压缩点之间的频谱对应
hamming1=[zeros（1855,1）'
hamming（387）'
zeros（1854,1）'
];
spectrum_sm=hamming1.*spectrum_s;
%%脉压
fft_ssb=ifft（fft（sb）.*conj（fft（s）））;
fft_smsb=ifft（fft（sb）.*conj（fft（sm,length（sb））））;
fft_spsb=ifft（fft（sb）.*conj（fftshift（spectrum_sm）））;
150
时域加窗后的参考信号频谱图
100
200
300
400
500
600
700
800
900
时域
加
窗前
的参
考
信号频
谱图
「
[
―L—
/
频率
f/Hz
x10
号
带加窗与否
-400
%%%%%%%%%%%%%%%%%%%%%%^斜加％%处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%带宽10MHztp=10e-6;
%脉宽10usu=B/tp;
%fs>
=2*B/tp*tauR0=3000;
%初始距离R=4500;
%距离波门
f0=60e6;
%%%%%%%%%%%%%%%%%%%%%%参考信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sref=exp（1i*pi*u*t.A2）;
%%%%%%%%%%%%%%%%%%%%%%%%回波信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Sb=rectpuls（t-2*R0/c,tp）.*exp（1j*pi*u*（t-2*R0/c）.A2）;
%%%%%%%%%%%%%%%%%%%%%%%%1频信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ssb=Sref.*conj（Sb）;
%%加窗
%w=hamming（502）'
%hamming=[zeros（749,1）'
w-min（w）,zeros（249,1）'
%hamming=abs（hamming）/max（hamming）;
hamming=[zeros（749,1）'
hamming（502）'
zeros（249,1）'
ssb0=hamming.*ssb;
spectrum_ssb0=fft（ssb0,fft_N）;
%一维距离像
spectrum_ssb=fft（ssb,fft_N）;
%%图6
plot（f,db（abs（fftshift（spectrum_ssb））/max（fftshift（spectrum_ssb））））
holdon
plot（f,db（abs（fftshift（spectrum_ssb0））/max（fftshift（spectrum_ssb
升级会员 