雷达信号处理基本流程.doc

基本雷达信号处理流程一、脉冲压缩窄带（或某些中等带宽）的匹配滤波：相关处理，用FFT数字化执行，即快速卷积处理，可以在基带实现（脉冲压缩)快速卷积,频域的匹配滤波脉宽越小,带宽越宽，距离分辨率越高 ；脉宽越大,带宽越窄，雷达能量越小，探测距离越近；D=BT(时宽带宽积)；脉压流程：频域：回波谱和参考函数共轭相乘时域：相关即输入信号的FFT乘上参考信号FFT的共轭再逆FFT；Sc=ifft（fft（Sb）.*conj(fft（S)））;Task1f0=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）^2).*exp（-2j*pi＊f0*tau);％回波信号 S=rectpuls(t-tp/2，tp）exp(i*pi*beita＊(t—tp/2).^2）；%发射信号（参考信号）So=ifft(fft（Sb)。

＊conj（fft（S)）)；％脉压figure(7)；plot（t＊c/2，db(abs(So）/max（So)))％归一化dBgrid on二、 去斜处理(宽带的匹配滤波）去斜处理“有源相关”，通常用来处理极大带宽的LFM波形（如果直接采样的话因为频带很宽所以在高频的时候需要的采样率就很大，采样点数就很多，所以要经过去斜处理）Stretch方法是针对线性调频信号而提出的,其方法是将输入信号与参考信号(经适当延迟的本振信号,延迟量通常由窄带信号测距结果估计出）混频，则每一个散射点就对应一个混频后的单频分量,对混频输出的信号进行DFT处理，即可获得目标的距离像，对参考信号的要求是应具有与输入信号相同的调频斜率去斜处理流程：混频过程为回波信号在时域与参考信号的共轭相乘混频后得到一个瞬时频率和目标距离成正比的单频信号，对其进行频谱分析即可得到目标的距离像；去斜处理一般情况下可降低信号带宽；％%％%%%%％％％%％%%%%％％%％％%％% 去斜处理仿真程序 ％%%％%%％％%%％％%％％％%%％％％%％％％clc;clear all;close all;B=10e6;％带宽10MHztp=10e—6;％脉宽10usk=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=2^nextpow2(N)；t=linspace（0，2＊R/c,N); %%％%%%％%%％％％%％%％%%%％％％％%%％ 参考信号 %％%％%％％%%％％％％％％％%％%%%％%％％%%%％％Sref=exp（2i*pi*f0＊t)。

＊exp（1i＊pi＊k＊t^2)； %%％%％%％%%％％%％％％％%%%％％%%％%％% 回波信号 %%%％%％％％%%％％％%%%％％%％％％％％％%％%Sb0=exp(1j*pi＊k*（t—2＊R0/c）^2）.＊exp（2j＊pi＊f0＊(t—2*R0/c）);Sb1=exp（1j*pi＊k*（t-2*R1/c）exp(2j*pi*f0*(t-2＊R1/c）)；Sb2=exp(1j*pi*k＊（t-2＊R2/c)exp(2j＊pi*f0＊(t—2＊R2/c））；Sb=Sb0+Sb1+Sb2; %％%％％%%％％％％%％％%％％％％％％%%%％%% 混频信号 %％％%%%％％％%％%%％％%%%%%%％％%%％%%SSb=Srefconj(Sb)；%去斜后时域信号spectrum=fft(SSb，fft_N）;%去斜后频域信号f=fs*(0:fft_N-1）/fft_N—fs/2;％从—fs/2到fs/2f=f＊c*tp/2/B；%瞬时频率对应的距离sf=exp(—j*pi/k*f^2）；%滤波器传输函数SSb=spectrumsf；％从频域去距离扭曲，实现了压缩和去RVPfigure;SSb=fftshift（SSb);SSb1=ifft（SSb）;%消除了距离扭曲和RVP的时域信号subplot（211）；plot(f,db（abs(SSb)/max（SSb）)）xlabel（'距离/m');grid onsubplot（212）；plot(f，abs(SSb)）xlabel（’距离/m')；grid on三、加窗信号的截取产生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，在FFT分析中为了减少或消除频谱能量泄漏及栅栏效应可采用不同的截取函数对信号进行截短，截短函数称为窗函数，简称窗。

％%%％%％％%%%%%%％%％%%%％ 窄带加窗处理 ％%％%％％％％％％％%％％％%%%％％%%%％clc;clear all；close all；f0=10e9;％载频 B=10e6;%信号带宽 tp=10e-6;%脉冲宽度 fs=100e6；%采样频率 k=B/tp;%LFM系数，线性调频率 R0=3000；％初始距离c=3e8;％光速 R=6000；tau=2*R0/c；N=round（2＊R/c*fs）;fft_N=2^nextpow2（N）;t=(0:fft_N-1）/fs;s=rectpuls（t-tp/2，tp）exp(j＊pi*k*（t—tp/2）.^2); %% 发射信号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）^2).*exp(-2j＊pi*f0＊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)）)；%%%%%％%%％%%％％％%％％％％％%%%％％％% 去斜加窗处理 ％％%%%%%％%%%%%％％%%%％％%％%％%%%%B=10e6；％带宽10MHz tp=10e—6;％脉宽10us u=B/tp;%LFM系数fs=50e6；％fs>=2＊B/tp*tau R0=3000;％初始距离 R=4500;%距离波门c=3e8； f0=60e6；%载频N=round（2＊R/c＊fs）; fft_N=2^nextpow2(N）; t=linspace(0,2＊R/c，N）；f=fs＊（0:fft_N—1）/fft_N-fs/2；％从-fs/2到fs/2％%%%％％％％％％%％％%％%％%%％％％％％％% 参考信号 %％％％％%％%％％%%％%％％%%%％%％%％%%%％％％Sref=exp(1i＊pi*u＊t.^2);%％％％%％%%％%%%％％％%%％％%%％％％％％％％ 回波信号 %％%％％%％％％％%%％％％%％%%%%％％％％%%％Sb=rectpuls(t-2＊R0/c,tp）。

exp(1j＊pi*u＊（t—2*R0/c).^2);%%％%％%％%％%%%％%%％％%％％％%%％%%％％ 混频信号 %%％％％%%％％%%％％%％％%％％％％%%%％％%%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； %％图6plot(f，db（abs(fftshift(spectrum_ssb）)/max（fftshift（spectrum_ssb)）)）hold onplot（f，db（abs（fftshift(spectrum_ssb0））/max（fftshift（spectrum_ssb0)))，'r’）hold off二、检测1、脉冲多普勒（PD处理）多普勒效应:fd=2v/c*f0,v为镜像速度；慢时间维上的采样点做FFT可以测出目标的速度；使用复信号：频率正负可测量目标速度的方向;clc；clear all;close all；f0=10e9；％载频tp=10e—6;%脉冲宽度B=10e6;％带宽fs=100e6；％采样频率R0=3000；％初始距离c=3e8;%光速R=4500;%距离波门 gate=R+tp*c/2；％距离波门加脉宽对应距离N=round(2＊gate/c＊fs);%波门内采样点个数fft_N=2^nextpow2（N）；t=0：1/fs：tp;%信号长度echo_t=linspace(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。

^2)；%回波信号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）.^2—1i＊pi*2*f0＊round(2*R0/c*fs)+1i＊2＊pi*(2*f0＊v/c）*(m-1）＊Tr）+sqrt(0.1)＊（randn(1，N)+1i*randn（1，N)）；end%脉压fft_n=2^nextpow2（length（t）+N—1）；fft_s=fft(s，fft_n);for m=1：。