• 大小: 94KB
    文件类型: .zip
    金币: 2
    下载: 2 次
    发布日期: 2021-06-03
  • 语言: 其他
  • 标签: 雷达  

资源简介

雷达中的运动补偿,keystone变换的CZT实现,利用仿真数据的验证

资源截图

代码片段和文件信息

%% 模拟产生匀速运动目标  
% 先做多普勒模糊补偿,后做CZT补偿距离徙动    
clear all;clc; close all;

 %% 雷达参数
 fs=15e6;                           % 基带采样率
 B=10e6;                            % 基带带宽
 ts=1/fs;
 tp=200e-6;                         % 脉宽
 kr=B/tp;                           % 调频斜率
 Tr=1e-3;                           % 脉冲重复周期
 PRF=1/Tr;
 c=3e8;
 fc=3e8;                            % 载频
 lamda=c/fc;
 N=200;                             % 脉冲数
 %% 目标参数
 N01=0:N-1;
 R0=5e4;                            % 初始距离 m
 v=499;                             % 径向运动速度m/s
 R=R0-v*N01*Tr;                     % 目标与雷达径向距离 m
 fd=2*v/lamda;                      % 多普勒频率  Hz
 
 %%回波信号
 t0=0:ts:tp-ts;                             %发射序列
 Tt=exp(j*pi*kr*(t0-tp/2).^2);              %发射波形
 tt=tp:ts:(Tr-ts);                          %采样序列
 n_fft=2^nextpow2(length(tt)-1);            %频域脉压FFT点数
 st=zeros(length(tt)N);                    %回波波形
for m=1:N
     tm(m)=2*R(m)/c-floor(2*R(m)/c/Tr)*Tr;
     st(ceil(tm(m)*fs):ceil(tm(m)*fs)+tp*fs-1m)=exp(j*pi*kr*(t0-tm(m)+ceil(tm(m)*fs)*ts-tp/2).^2 -j*2*pi*2*R(m)/lamda);   %+j*2*pi*fd*t0
end
% SNR=-18; %输入信噪比
% st=noisegen(stSNR);
%figure;plot(real(st(:1)));
 
%% 脉压
fft_ht=fft(Ttn_fft);                       %匹配函数
fft_st=fft(stn_fft);
Array_ones=ones(N1);
fft_ht_N=Array_ones*fft_ht;
xt=fft_st.*conj(fft_ht_N.‘);                %共轭相乘
% yt=ifft(xt);                                %返回到时域,脉压
% figure;contour(abs(yt)); 
% % xlabel(‘距离单元‘‘fontsize‘mesh_TextFontSize);ylabel(‘脉冲数‘‘fontsize‘mesh_TextFontSize)
% % set(gca‘fontname‘‘Helvetica‘‘fontsize‘mesh_NumFontSize);
% % set(gcf‘windowstyle‘‘normal‘);
% % set(gcf‘Position‘[500300mesh_widthmesh_height]);
% title(‘脉压结果‘)

%% MTD
% inte_xt=fft(yt.‘);                          %脉间FFT
% figure;mesh(abs(inte_xt));
% % xlabel(‘距离单元‘‘fontsize‘mesh_TextFontSize);ylabel(‘多普勒单元‘‘fontsize‘mesh_TextFontSize);
% % zlabel(‘幅度‘‘fontsize‘mesh_TextFontSize)
% % set(gca‘fontname‘‘Helvetica‘‘fontsize‘mesh_NumFontSize);
% % grid on;
% % set(gcf‘windowstyle‘‘normal‘);
% % set(gcf‘Position‘[500300mesh_widthmesh_height]);
% title(‘未补偿‘)

%% keystone距离徙动补偿
N01=0:N-1;
k1=floor(fd/PRF);     % 多普勒模糊数,速度模糊部分的补偿
%%%% 频谱正半部分补偿
for l=0:min(floor(n_fft*B*1.15/fs/2)n_fft/2)-1      % 1:min(floor(n_fft*B*1.1/fs/2)n_fft/2)-1  选择频谱主要能量 1.1*n_fft*B/fs/2
    phi=(fc+l*fs/n_fft)/fc;                         % (fc+l*fs/n_fft)/fc 为对应快时间频率(l)处的系数,大小=1左右
    f_amb_xt(l+1:)=xt(l+1:).*exp(-j*2*pi*k1*phi*N01);   %先做多普勒模糊补偿        
    w=exp(-j*2*pi*phi/N);
    Czt_KT(l+1:) = czt(f_amb_xt(l+1:)Nw);       % 慢时间维CZT实现距离徙动补偿 
end
%%%% 频谱付半部分补偿
for l=max(-floor(n_fft*B*1.15/fs/2)-n_fft/2):-1     % max(-floor(n_fft*B*1.1/fs/2)-n_fft/2):-1  选择频谱主要能量 1.1*n_fft*B/fs/2
    phi=(fc+l*fs/n_fft)/fc;               
    f_amb_xt(l+1+n_fft:)=xt(l+1+n_fft:).*exp(-j*2*pi*k1*phi*N01);     %先做多普勒模糊补偿    
    w=exp(-j*2*pi*phi/N);
    Czt_KT(l+1+n_fft:) = czt(f_a

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      102504  2017-04-15 05:39  CZT仿真.docx
     文件         817  2017-04-15 05:39  tryCZT2.m
     文件        3450  2017-04-15 05:39  keystone1.m

评论

共有 条评论