资源简介
雷达中的运动补偿,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
- 上一篇:最小重量机器问题报告
- 下一篇:用汇编语言做的时钟,可以设置闹钟
相关资源
- 数字正交上变频器AD9857在高频雷达系
- ALOS-PALSAR雷达数据介绍
- 基于地震属性体的三维探地雷达快速
- PolSARpro Introduction PolSARpro操作说明
- ADS 24GHz雷达链路仿真,2019版ADS
- 机载雷达数据产品在油气输送管道工
- 雷达手册第三版中文 第1章 雷达概论
- 机载雷达导论.rar
- 倒车雷达算法文档
- 详解多功能相控阵雷达.pdf
- SNR-提取波的初至时
- 基于三维激光雷达的自动驾驶车辆高
- 雷达监测在滑坡信息化防治技术中的
- 微波(雷达)感应模块原理调试
- 基于共中心点道集约束的探地雷达波
- 雷达对抗原理_赵国庆_课后答案(部分
- 绝地求生雷达最新代码.sh
- 相控阵雷达系统设计
- 77GHz车载毫米波中远距雷达天线阵列设
- 雷达的真实回波数据仿真
- TI毫米波雷达原理
- 捷变频雷达信号代码
- 经典雷达资料
- 毫米波理论与技术 雷达,电子对抗,
- ISAR成像原理
- 雷达 导航与电子对抗
- 雷达原理和雷达对抗资料
- 雷达原理与对抗技术\\雷达原理与对抗
- 点云拼接和拟合算法
- 雷达成像原理Word完整版
评论
共有 条评论