雷达信处理基本流程.docx
- 文档编号:1759272
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:23
- 大小:958.46KB
雷达信处理基本流程.docx
《雷达信处理基本流程.docx》由会员分享,可在线阅读,更多相关《雷达信处理基本流程.docx(23页珍藏版)》请在冰豆网上搜索。
雷达信处理基本流程
基本雷达信号处理流程
一、脉冲压缩
窄带(或某些中等带宽)的匹配滤波:
相关处理,用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*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);
1
%发射信号(参考信号)
1
脉压
J
丄
I
1
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;
c=3e8;
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('距离/m');
gridon
二、加窗
信号的截取产生了能量泄漏,而用FFT算法计算频谱又产生了栅栏效应,在FFT分析中为了减少或消除频谱能量泄漏及栅栏效应可采用不同的截取函数对信号进行截短,截短函数称为窗函数,简称窗。
%%%%%%%%%%%%%%%%%%%%窄带加窗处理%%%%%%%%%%%%%%%%%%%%%%%%clc;clearall;closeall
f0=10e9;%载频
fs=100e6;%采样频率
c=3e8;%光速
tau=2*R0/c;
N=round(2*R/c*fs);
fft_N=2Anextpow2(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*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
时域加窗后的参考信号频谱图
150
100
100
200
300
400
500
600
700
800
900
1000
时域
加
窗前
的参
考
信号频
谱图
「
[
―L—
/
J
1
2
3
4
6
7
8
9
10
频率
f/Hz
7
x10
号
f/Hz
考
-50
-100
-150
带加窗与否
-200
-250
-300
-350
-400
-100
-150
-200
-50
-250
-300
-350
%%%%%%%%%%%%%%%%%%%%%%^斜加%%处理%%%%%%%%%%%%%%%%%%%%%%%%%%%%
B=10e6;%带宽10MHztp=10e-6;%脉宽10usu=B/tp;%LFM系数
fs=50e6;%fs>=2*B/tp*tauR0=3000;%初始距离R=4500;%距离波门
c=3e8;f0=60e6;%载频
N=round(2*R/c*fs);fft_N=2Anextpow2(N);t=linspace(0,2*R/c,N);
f=fs*(0:
fft_N-1)/fft_N-fs/2;%从-fs/2至Ufs/2
%%%%%%%%%%%%%%%%%%%%%%参考信号%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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);
f=f*c*tp/2/B;%瞬时频率对应的距离
figure;%%图6
plot(f,db(abs(fftshift(spectrum_ssb))/max(fftshift(spectrum_ssb))))
holdon
plot(f,db(abs(fftshift(spectrum_ssb0))/max(fftshift(spectrum_ssb
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 雷达 处理 基本 流程