1、） %归一化 dB-50-100-150-200-25030004000500060007000-300-3502000二、去斜处理（宽带的匹配滤波）去斜处理“有源相关”，通常用来处理极大带宽的 LFM 波形（如果直接采样的话 因为频带很宽所以在高频的时候需要的采样率就很大， 采样点数就很多， 所以要 经过去斜处理）Stretch 方法是针对线性调频信号而提出的， 其方法是将输入信号与参考信号 （经 适当延迟的本振信号， 延迟量通常由窄带信号测距结果估计出） 混频 ，则每一个 散射点就对应一个混频后的单频分量，对混频输出的信号进行 DFT 处理 ，即可 获得目标的距离像，对参考信号的要求是应具
2、有与输入信号相同的调频斜率。去斜处理流程 ：输入信号 频谱分析 输出信号混频过程为回波信号在时域与参考信号的共轭相乘 混频后得到一个瞬时频率和目标距离成正比的 单频信号 ，对其进行频谱分析即可 得到目标的距离像；去斜处理一般情况下可降低信号带宽；%去%斜%处%理 仿真程序 % clc;clear all ;close all ;B=10e6; %带宽 10MHztp=10e-6; %脉宽 10usk=B/tp; %LFM 系数fs=50e6;R0=3e3;R1=2000;R2=3500;R=5000;f0=60e6;N=round（2*R/c*fs）; fft_N=2nextpow2（N）;
3、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）.2）.*exp（2j*pi*f0*（t-2*R1/c）;Sb2=exp（1j*pi*k*（t-2*R2/c）.2）.*exp（2j*pi*f0*（t-2*R2/c）;Sb=Sb0+Sb1+Sb2;f=fs*（0:fft_N-1）/fft_N-fs/2; %从 -fs/2
4、到 fs/2f=f*c*tp/2/B; %瞬时频率对应的距离sf=exp（-j*pi/k*f.2）; %滤波器传输函数SSb=spectrum.*sf; %从频域去距离扭曲，实现了压缩和去 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（ 三、加窗信号的截取产生了能量泄漏， 而用 FFT算法计算频谱又产生了栅栏效应，
5、 在 FFT 分析中为了减少或消除频谱能量泄漏及栅栏效应可采用不同的截取函数对信号 进行截短，截短函数称为窗函数，简称窗。%窄%带%加%窗 处理 %clc;close all f0=10e9; %载频 fs=100e6; %采样频率 c=3e8; %光速 tau=2*R0/c;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
6、（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（
7、sm,length（sb）; fft_spsb=ifft（fft（sb）.*conj（fftshift（spectrum_sm）;归 一 化 发 射 信 号 频 谱23456789频 率 f/Hzx 1050150时 域 加 窗 后 的 参 考 信 号 频 谱 图2003004005006007008009001000时域加窗前的参考信号频谱图01-400窄 带 加 窗 与 否 对 比 图距 离 /m窄 带 频 域 加 窗 与 否 对 比 图%去%斜%加%窗%处 理 % %带宽 10MHz tp=10e-6; %脉宽 10us u=B/tp; %fs=2*B/tp*tau R0=3000; %
8、初始距离 R=4500; %距离波门 f0=60e6; %载频%参%考%信%号 % 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（
9、249,1）ssb0=hamming.*ssb;spectrum_ssb0=fft（ssb0,fft_N）; %一维距离像spectrum_ssb=fft（ssb,fft_N）; %图 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-120度幅化归-4000消 除 了 RVP 和 距 离 扭 曲 的 混 频 信 号-3000-2000-1000检
10、测1、脉冲多普勒（ PD处理） 多普勒效应： fd=2v/c*f0 ，v 为镜像速度； 慢时间维上的采样点做 FFT可以测出目标的速度； 使用复信号：频率正负可测量目标速度的方向； clc;f0=10e9; %脉冲宽度 %带宽fs=100e6; %采样频率R0=3000; %初始距离 %光速R=4500;gate=R+tp*c/2; %距离波门加脉宽对应距离N=round（2*gate/c*fs）; %波门内采样点个数fft_N=2nextpow2（N）;t=0:1/fs:tp; %信号长度echo_t=linspace（0,2*gate/c,N）; %波门长度 %调频系数Tr=100e-6;
11、 %脉冲重复周期CPI=64; %总脉冲个数v=60; %目标速度，朝向雷达%发射信号s=exp（i*pi*k*t.2）;%回波信号for m=1:CPIsb（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=2nextpow2（length（
12、t）+N-1）;fft_s=fft（s,fft_n）;1: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,:）;）=20*log10（z1（m,:maxval,maxpo=max（z1（m,:%FFTfor fm=1:Ndop（:,fm）=fft（fft_ssb（:,fm）;a_dop（:,fm）=fftshift（abs（dop（:,fm）;%求极大值对应的坐标maxva,max_v=max（a
13、_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/f0mesh（echo_t*c/2,linspace（-75,75,64）,a_dop）;axis tight ;xlabel（ 距离： mylabel（ 速度： m/s title（ 二维距离 - 多普勒平面 v_pd =60.9375 v_max =75 det_v =2.34382、形心法测距测速%形%心%法% % % clc;clo
14、se all %载频 tp=10e-6; %带宽 fs=100e6;N=4096; %此为培训期间数据，实际情况可以根据波门和信号宽度求出 Nsnr=20;sigmaf=12/（10（snr/10）; %发射信号% figure;% plot（t,real（s）% xlabel（ 时间 /s）;% ylabel（ 幅度 % title（ 发射信号实部 % grid onCPI taum=2*（R0-m*Tr*v）/c;sb=rectpuls（t-taum-tp/2）.*exp（j*pi*k*（t-taum-tp/2）.2）.*exp（-2j*pi*f0%加噪声的回波*taum）; %回波信号
15、sb_noise=sb+sqrt（sigmaf/2）*（randn（1,N）+1i*randn（1,N）; %脉压处理fft_ssb_snr=ifft（fft（sb_noise）.*conj（fft（s）;Group（m,1:N）=fft_ssb;Group_snr（m,1:N）=fft_ssb_snr; imagesc（t*c/2,1:CPI,abs（fft_ssb）imagesc（t*c/2,1:CPI,abs（fft_ssb_snr）for n=1:Group2=fft（Group（1:CPI,n）; %纵向做 FFTGroup2_2=fftshift（abs（Group2）;Group
16、3（n,1:CPI）=Group2_2;for n1=1:Group2_snr=fft（Group_snr（1:CPI,n1）;Group2_2_snr=fftshift（abs（Group2_snr）;Group3_snr（n1,1:CPI）=Group2_2_snr;mesh（abs（Group3）mesh（abs（Group3_snr）line,row=find（abs（Group3）=max（max（abs（Group3）;line_snr,row_snr=find（abs（Group3_snr）=max（max（abs（Group3_snr）; Range=t*c/2;PRF=1/T
17、r;fd=（-CPI/2:CPI/2-1）*PRF/CPI;v=fd*c/2/f0;for Ra=line-3:line+3 amp=abs（Group3（Ra,row）; C（Ra）=amp*Range（Ra）; D（Ra）=sum（amp）;sum（C）/sum（D）for V=row-3:row+3 index=abs（Group3（line,V）; E（V）=index*v（V）; F（V）=sum（index）;sum（E）/sum（F）for Ra_snr=line_snr-3:line_snr+3 amp_snr=abs（Group3_snr（Ra_snr,row_snr）; C
18、_snr（Ra_snr）=amp_snr*Range（Ra_snr）; D_snr（Ra_snr）=sum（amp_snr）;sum（C_snr）/sum（D_snr）for V_snr=row_snr-3:row_snr+3 index_snr=abs（Group3_snr（line_snr,V_snr）; E_snr（V_snr）=index_snr*v（V_snr）; F_snr（V_snr）=sum（index_snr）;sum（E_snr）/sum（F_snr）结果：ans =3.0000e+003ans =60.3560ans =60.35333、信号检测mean_noise =0
19、.0010 + 0.0079ivar_noise = 1.0123pf1 =1.0000e-003高斯白噪声取模值后的波形以及检测门限3.5检测概 率相对于 SNR曲线9876543200000000率概测检4 6 8 10 12 14 16 18 20SNR/dB4、单脉冲测角仿真单脉冲跟踪雷达是通过比较来自两个或多个同时波束的信号获得目标角位置信 息的一种雷达；目前常用的单脉冲测角方法主要有幅度和差单脉冲测角和相位和差单脉冲测角。 幅度和差单脉冲测角 通过比较两个 相位中心重合 但指向不同 的波束得到目标角 度信息；相位和差单脉冲测角 则通过比较两个 相位中心有一定距离 但波束指向相 互平
20、行 的波束得到目标角度信息。相位和差单脉冲与幅度和差单脉冲的相似之处在于： 目标角度坐标都是由一个和 通道和两个差通道来提取的。 主要不同之处在是， 幅度和差单脉冲产生的四个信 号具有相同的相位但具有不同的幅度， 而相位和差单脉冲信号具有相同的幅度但 有不同的相位。 相位和差单脉冲对每个坐标系 （方位和俯仰坐标） 使用最少由两 个阵元组成的阵列天线。 相位误差信号是由于不同天线阵元产生的信号之间的相 位差来计算得出的。%单%脉%冲%测 角仿真 %f0=10e6; %目标距 0 号阵元的距离 d=10; %阵列接收天线之间的距离 theta0=0.2*pi/180; %目标角度%R0=R0+d*sin（theta0）/2;lamda=c/f0; tau0=2*R0/c; %到 0 号阵元的时延theta=linspace（-1*pi/180,1*pi/180,1000）;thetaP=0.15*pi/180; %偏置角N=4; %天线个数% 幅度和差单脉冲测角%相同相位不同幅度%波束形
