RD算法及其距离徙动校正(原理及matlab仿真)
RD算法流程图:
图1 RD算法流程图
因为是在word里写的,这里直接上截图了,不然mathtype公式显示不出来。
对于RD算法的原理这里就不展开讨论了,重点分析距离徙动校正的原理。
仿真程序,F_s=2.5B,距离徙动最大2.92m,对应7.3个距离单元。PRF取1.2倍多普勒带宽,目标点相对于参考点的位置为(-50,50,0),则仿真结果如下:观察图2,可以看出经过距离徙动校正后,距离多普勒域图像由弯曲变直;观察图3,可以看出RD算法比简单的二维脉压成像结果分辨率高,即经过距离徙动校正能提高分辨率。
图2-1 距离多普勒域成像结果
图2-2 放大后的距离多普勒域成像结果
图3-1 二维脉压成像结果
图3-2 RD算法成像结果
程序
%仿真具有3个以上距离徙动回波信号,进行RD成像
close all;clear;%基本参数
c = 3*10^8; %波速
fc = 5.3*10^9; %载频5.3GHz
lamda = c/fc; %波长
v = 150; %雷达速度150m/s
B = 150*10^6; %发射信号带宽150MHz
Daz = 2; %方位天线尺寸2m
R0 = 20*10^3; %参考点斜距20km
Tr = 1.5*10^-6; %脉冲持续时间
Kr = B/Tr; %线性调频信号调频斜率
H = 1000; %雷达高度
Yc = sqrt(R0^2-H^2); %成像中线
center = [0,Yc,0]; %场景中心坐标
CR = Daz/2; %横向分辨率
thetaaz = lamda/Daz; %波束宽度
Fs = 2.5*B; %采样频率
Dsar = lamda*R0/(2*CR); %合成孔径长度
Tsar = Dsar/v; %一个合成孔径时间
Ka = -2*v^2/(lamda*R0); %多普勒调频斜率
Ba = abs(Ka*Tsar); %多普勒带宽
PRF = 1.2*Ba; %脉冲重复频率
X = 150;Y = 150; %场景大小
Rmin = sqrt((Yc-Y/2)^2+H^2); %场景最小距离
Rmax = sqrt((Yc+Y/2)^2+H^2+(Dsar/2)^2); %场景最大距离
Nfast = ceil(((2*(Rmax-Rmin)/c+Tr)*Fs)); %快时间维采样点数
tf = 2*R0/c+(-(Nfast/2):(Nfast/2)-1)/Fs; %快时间采样序列
Nslow = ceil((X+Dsar)/v*PRF); %慢时间维(方位向)采样点数
ts = (-(Nslow/2):(Nslow/2)-1)/PRF; %慢时间采样序列
pos = [-50,50,0,3]; %目标相对中心点位置,[x,y,z,rcs],第四列为后向反射系数
disp('目标位置(方位,斜距,高度):');
Rpos(1:3) = pos(1:3)+center %目标绝对位置
Rpos(:,4) = pos(:,4);%回波信号
signal = zeros(Nfast,Nslow);
Xs = ts.*v-Rpos(1);
Ys = 0-Rpos(2);
Zs = H-Rpos(3);
sigma = Rpos(4); %反射系数
R = sqrt(Xs.^2+Ys^2+Zs^2); %斜距
tau = 2*R/c; %时延
Tfast = tf'*ones(1,Nslow)-ones(Nfast,1)*tau; %列是快时间 行是慢时间
Phase = pi*Kr.*Tfast.^2-(4*pi/lamda)*ones(Nfast,1)*R; %相位延迟
signal = signal+sigma*exp(j*Phase).*(abs(Tfast)<=Tr/2).*(ones(Nfast,1)*(abs(Xs)<=Dsar/2)); %回波
S = fftshift(fft(fftshift(signal)));%二维脉压
hf = exp(j*pi*Kr*(tf-2*R0/c).^2).*(abs(tf-2*R0/c)<=Tr/2); %距离向参考函数
Hf = (fftshift(fft(fftshift(hf))).')*ones(1,Nslow);
ComF = S.*conj(Hf); %距离向匹配滤波
Sr = fftshift(ifft(fftshift(ComF))); %距离向IFFT
Coms = fftshift(fft(fftshift(Sr.'))).'; %方位向FFT
hs = exp(j*pi*Ka*ts.^2).*(abs(ts)<Tsar/2); %方位向参考函数
Hs = ones(Nfast,1)*fftshift(fft(fftshift(hs)));
ComS = Coms.*conj(Hs); %方位向匹配滤波
Saz = fftshift(ifft(fftshift(ComS.'))).';%RD算法,sinc插值
Coms_rcmc = zeros(Nfast,Nslow);
N = 6; %插值点数
Rp = sqrt(sum((Rpos(2:3)-[0,H]).^2)); %目标到雷达的最近距离
h = waitbar(0,'Sinc插值中......'); %生成一个进度条
for m = 1:Nslow %慢时间for n = N/2+1:Nfast %快时间%距离偏移量deltaR = (lamda/v)^2*(Rp+(n-Nfast/2)*c/2/Fs)*((m-Nslow/2)/Nslow*PRF)^2/8;DU = deltaR/(c/2/Fs); %偏移距离单元deltaDU = DU-floor(DU); %偏移距离单元小数部分for k = -N/2:N/2-1if (n+floor(DU)+k)>Nfast %超出边界Coms_rcmc(n,m) = Coms_rcmc(n,m)+Coms(Nfast,m)*sinc(DU-k);elseComs_rcmc(n,m) = Coms_rcmc(n,m)+Coms(n+floor(DU)+k,m)*sinc(deltaDU-k);endendendwaitbar(m/Nslow);
end
close(h); %关闭进度条
ComS_rcmc = Coms_rcmc.*conj(Hs); %方位向压缩
Saz_rcmc = fftshift(ifft(fftshift(ComS_rcmc.'))).';%画图
rf = c*tf/2; %距离
az = v*ts; %方位
faz = (-Nslow/2:Nslow/2-1)/Nslow*PRF; %多普勒频率
figure(1);
[f,Rf] = meshgrid(faz,rf);
subplot(121);
mesh(f,Rf,abs(Coms));view(0,90);
title('(a) 未RCMC');
xlabel('多普勒/Hz');ylabel('斜距R/m');
subplot(122);
mesh(f,Rf,abs(Coms_rcmc));view(0,90);
title('(b) RCMC');
xlabel('多普勒/Hz');ylabel('斜距R/m');
figure(2);
[Az,Rf] = meshgrid(az,rf);
mesh(Az,Rf,abs(Saz));view(0,90);
title('二维脉压成像结果');
xlabel('方位x/m');ylabel('斜距R/m');
figure(3);
mesh(Az,Rf,abs(Saz_rcmc));view(0,90);
title('RD算法成像结果');
xlabel('方位x/m');ylabel('斜距R/m');
RD算法及其距离徙动校正(原理及matlab仿真)相关推荐
- matlab 脉冲压缩算法,线性调频脉冲压缩原理及其MATLAB仿真+程序
摘 要:脉冲压缩技术在现代雷达系统和超带宽通信系统等领域具有越来越重要的应用.本文首先阐述了线性调频脉冲压缩的基本原理,并提出了增大信号等效带宽可以提高雷达距离分辨率,然后分析了匹配滤波器的特性,给出 ...
- 光栅原理及其c语言仿真,光栅原理及MATLAB仿真
<光栅原理及MATLAB仿真>由会员分享,可在线阅读,更多相关<光栅原理及MATLAB仿真(8页珍藏版)>请在人人文库网上搜索. 1.取样光纤光栅的原理及基于MATLAB的反射 ...
- matlab 光栅 傅里叶,光栅原理及MATLAB仿真汇编.doc
光栅原理及MATLAB仿真汇编 取样光纤光栅的原理及基于MATLAB的反射谱仿真 取样光纤光栅其实与相称光纤光栅基本上一致,不同的地方在于,相移光栅是在均匀布拉格光栅的某一点处引入相移,导致在反射谱中 ...
- SAR成像系列:【9】合成孔径雷达(SAR)成像算法-波数域(omega-K)成像算法[也叫距离徙动(RM)算法](附Matlab代码)
波数域()成像算法作为本系列的最后一种成像算法介绍.关于SAR成像的其他的各种改进算法就不一一列举了.在实际成像中,万变不离其踪,最主要的是关注成像的几何模型,再根据指标选择不同的基础成像算法,然后进 ...
- MIMO大规模天线阵列原理与matlab仿真(含GUI)
一 设计实现流程 1.MIMO系统原理 传统的无线通信系统是采用一个发送天线和一个接收天线的通信系统,即单输入单输出(SISO)天线系统.但由于单天线系统的信道容量较低,不能满足4G.5G数据传输需求 ...
- 数字PID控制算法原理及Matlab仿真
引言 最近碰到一个项目需要用到PID控制算法,于是在网上找了一些资料学习了一下,发现网上对于PID算法的Matlab仿真方面的内容比较少,所以我就把我自己所学习到的内容分享给大家.本次博文主要介绍了位 ...
- 基于FMCW的测距原理及matlab仿真
FMCW是什么 FMCW(Frequency Modulated Continuous Wave),即调频连续. FMCW实现测距 利用发射信号与接收信号进行混频,得到包含目标距离和速度信息的中频信号 ...
- 光栅原理及其c语言仿真,取样光栅原理及MATLAB仿真.doc
取样光纤光栅的原理及基于MATLAB的反射谱仿真 张睿 摘要 文章主要运用了基于耦合模理论的传输矩阵法来分析取样光栅的原理,并利用MATLAB模拟和分析了取样光栅长度.调制折射率强度.取样光栅节点的长 ...
- 使用EM算法估计GMM参数的原理及matlab实现
相关数学概念 协方差矩阵 多维高斯分布 其中k=n,即x的维度. GMM的原理 GMM,高斯混合模型,是一种聚类算法. 1.GMM概念: -将k个高斯模型混合在一起,每个点出现的概率是几个高斯混合的结 ...
最新文章
- pyspark dataframe基本用法
- [Java] Welcome to Artifactory 使用说明
- python下载图片的命令_网上的图片不知道怎么批量下载?python教你怎么把网站上面的图片都爬下来...
- C语言 数组排序 – 选择法排序 - C语言零基础入门教程
- db2 语句包括不必要的列表_DB2的一些常用SQL写法(转)
- 序列化模块和hashlib模块
- Java基础之StringBuffer和StringBuilder的区别
- mysql基本sql语句总结(一)
- SSD,单次多框检测器Single Shot Multibox Detector,超越YOLO和Fast-RCNN
- python 气泡图 聚类_用SPSS制作气泡图
- 拉格朗日插值一个公式概括
- QQ、微信头像照片模糊怎么办?如何处理高清?
- phython基础之数据类型和变量
- Java项目:JSP会议-会议室管理系统
- 服务器vga转hdmi显示器不亮,如何排除HDMI转VGA的常见故障_排除故障的四种方法
- cnn和rnn可以结合使用吗,rnn和cnn优点缺点对比
- 白炽灯和led哪个护眼?分享真正适合孩子的护眼台灯
- 基于抽象语法树的Online Judge语义查重系统设计与实现(毕设)
- GBase 8a集群目录结构
- ciscn_2019_c_1 1
热门文章
- 屏蔽不显示VS警告warning
- 自由软件之父:苹果微软罪恶滔天
- 迁移gogs直接拷贝其gogs-repositories导致hook脚本中的相关路径与实际不相符
- csv文件和excel文件
- Tsar学习笔记---tsar tsar2db mysql 集成调测篇
- springboot查看请求ip、ip归属地、设备信息、本地MAC地址、操作系统信息
- 美泰推出首款获CarbonNeutral(R)认证的Matchbox(R) Tesla Roadster压铸模型车,使用9
- java利用栈进行进制转换
- 怎样填写备注是什么意思_手机号备注234824是什么意思什么梗 数字梗意思大全...
- 314-Leetcode 堆排序