SAR成像(一):线性调频信号(LFM)和脉冲压缩
文章仅为个人理解,如有不妥之处欢迎指正。
1、脉冲压缩目的和意义
在雷达系统中,距离分辨率是很重要的一项性能指标,它决定了雷达能分辨两个单元之间的最小距离,更通俗的,即雷达分辨相离较近物体的能力。而决定雷达分辨率的是发射信号带宽,所以我们想要提高距离分辨率,就需要发射更窄的脉冲信号,但是,窄脉冲就意味着发射信号能量小,导致探测距离短,所以,对于一般的脉冲信号(脉宽、时宽乘积约为1),比如矩形脉冲信号,带宽和时宽不能同时增大,因此距离分辨率和探测距离是一对矛盾。
脉冲压缩技术就能比较好的解决上述两个参量之间的矛盾,用宽脉冲发射信号,保证足够的探测距离,在接收端用相应的匹配滤波器,通过脉冲压缩技术得到窄脉冲,以此来提高距离分辨率。
2、线性调频信号
大时宽的宽频信号有很多形式,在雷达系统中最常用的是线性调频(Linear Frequency Modulation)脉冲信号。LFM脉冲信号的载频在脉冲宽度范围内线性变化,因此也有更宽的带宽。
LFM信号的复数表达式为:
s
(
t
)
=
r
e
c
t
(
t
τ
)
e
j
2
π
(
f
0
+
K
t
2
/
2
)
s(t)=rect(\frac{t}{\tau})e^{j2\pi(f_{0}+Kt^{2}/2)}
s(t)=rect(τt)ej2π(f0+Kt2/2)
其中,
f
0
f_{0}
f0为初始载频,
τ
\tau
τ为脉冲宽度,
K
K
K为发射线性调频信号的调频斜率,且
r
e
c
t
(
t
τ
)
=
{
1
,
∣
t
∣
≤
τ
/
2
0
,
∣
t
∣
≥
τ
/
2
rect(\frac{t}{\tau})=\begin{cases} 1, \text{ $|t|\leq \tau/2$} \\ 0, \text{ $|t|\geq\tau/2$} \end{cases}
rect(τt)={1, ∣t∣≤τ/20, ∣t∣≥τ/2
K
=
B
/
τ
K=B/\tau
K=B/τ
B
B
B为信号带宽。
信号的瞬时频率
f
(
t
)
f(t)
f(t)为:
f
(
t
)
=
1
2
π
d
d
t
[
2
π
(
f
0
t
+
K
t
2
/
2
)
]
=
f
0
+
K
t
f(t)=\frac{1}{2\pi}\frac{d}{dt}[2\pi(f_{0}t+Kt^2/2)]=f_{0}+Kt
f(t)=2π1dtd[2π(f0t+Kt2/2)]=f0+Kt
当
B
B
B=20MHz,
τ
\tau
τ=10
μ
s
\mu s
μs时,线性调频信号波形和频谱如下图所示:
还有一点不得不强调,LFM脉冲信号的带宽和时宽都是可以自己选择的,不像矩形脉冲信号那样带宽和时宽相互制约。
3、脉冲压缩
脉冲压缩的理论基础是匹配滤波。对雷达接收的回波信号进行脉冲压缩后有两个好处:增加信噪比;压缩信号宽度。
匹配滤波是一种最优滤波器,只要我们给予输入的是某一确知信号,并加上白噪声,那么就能让输出的信噪比达到最大。
将回波信号记为
f
(
t
)
f(t)
f(t),设滤波器的传输函数为
H
(
w
)
H(w)
H(w),则为了满足最大输出信噪比,
H
(
w
)
H(w)
H(w)的表达式为:
H
(
w
)
=
K
F
∗
e
−
j
w
t
0
H(w)=KF^*e^{-jwt_{0}}
H(w)=KF∗e−jwt0
上式中,
F
∗
F^*
F∗为
f
(
t
)
f(t)
f(t)傅里叶变换的共轭;
t
0
t_0
t0是物理器件的时间延时;
K
K
K为增益常数。
从时域角度看,匹配滤波器也可以用冲激函数
h
(
t
)
h(t)
h(t)来表示,则上式描述的滤波器冲击响应为:
h
(
t
)
=
K
f
∗
(
t
0
−
t
)
h(t)=Kf^*(t_0-t)
h(t)=Kf∗(t0−t)
脉冲压缩要在雷达后端接收机由数字信号处理器件完成,而在这个过程中,由于器件速度的限制,脉冲压缩的过程通常需要在零中频进行,此时,线性调频信号的表达式变为:
s
(
t
)
=
r
e
c
t
(
t
τ
)
e
j
π
K
t
2
s(t)=rect(\frac{t}{\tau})e^{j\pi Kt^{2}}
s(t)=rect(τt)ejπKt2少了初始频率
f
0
f_0
f0,从
h
(
t
)
h(t)
h(t)的表达式可以看出,如果不考虑器件延时
t
0
t_0
t0和增益常数
K
K
K,匹配滤波器的冲击响应是滤波器输入信号
f
(
t
)
f(t)
f(t)的共轭倒置,即:
h
(
t
)
=
f
∗
(
−
t
)
=
r
e
c
t
(
t
τ
)
e
−
j
π
K
t
2
h(t)=f^*(-t)=rect(\frac{t}{\tau})e^{-j\pi Kt^{2}}
h(t)=f∗(−t)=rect(τt)e−jπKt2滤波器输出可以由
f
(
t
)
、
h
(
t
)
f(t)、h(t)
f(t)、h(t)卷积得到,也可以由他们的傅里叶变换乘积再经过傅里叶反变换得到。在实际应用中,通过卷积直接做脉冲压缩的运算量很大。,更多是从频域角度来实现匹配滤波。
4、仿真
参数名称 | 参数值 |
---|---|
带宽(B) | 20MHz |
脉冲宽度(tao) | 10e-6s |
脉冲重复间隔 (T) | 1e-4s |
采样频率(fs) | 20MHz |
目标1距离(R1) | 5000m |
目标2距离(R2) | 8000m |
目标3距离(R3) | 8300m |
目标4距离(R4) | 9500m |
得到回波如下图:
对回波进行脉冲压缩后的波形如下图:
可以看到,由于目标2、目标3和目标3相距较近,他们的回波信号由接收机接收后会混在一起,难以分辨,但是经过脉冲压缩后,可以很容易的分辨出来。
附上仿真程序:
clear;clc;close all
T = 1e-4; %最大探测距离 c*T/2=15000m
tao = 10e-6;
B = 20e6;
fs = 2*B;
c = 3e8; %光速
K = B/tao;
R = [5000 8000 9300 8500];
rcs = [0.2 0.5 0.3 0.35];
delay = 2*R./c;
t=-tao/2:1/fs:tao/2-1/fs;
x_linear=exp(1j*pi*(K.*(t.^2))); %发射信号
figure
subplot(211);plot(t,real(x_linear));title('时域波形')
X=fftshift(fft(x_linear));
f=linspace(0,fs,length(t))-fs/2;
subplot(212);plot(f,abs(X));title('频域')
tt = 0:1/fs:T-1/fs;
echo = zeros(1,length(tt));
for i = 1:length(R)
% 回波信号
echo = echo + randn(1,length(tt))*0.2 + ...
rcs(i)*rectpuls(tt-delay(i)-tao/2,tao).*exp(1j*pi*K*(tt-delay(i)-tao/2).^2);
end
figure;
plot(tt,real(echo));xlabel('time');title('回波')
h=exp(1j*pi*(K.*(t.^2))); %匹配滤波器冲击响应
output=ifft(conj(fft(h,length(echo))).*fft(echo,length(echo))); %h和echo做fft的长度要一样
figure;
subplot(211);plot(0:1/fs:T-1/fs,abs(output));title('脉冲压缩输出波形');xlabel('time');
subplot(212);plot((0:1/fs:T-1/fs)*c/2,db(abs(output)/max(abs(output))),'r');
title('脉压结果分贝图');ylabel('dB');xlabel('m');
相关内容:
SAR成像(零):【总结】SAR成像原理和仿真实现
SAR成像(一):线性调频信号(LFM)和脉冲压缩
SAR成像(二):高方位向分辨率原理
SAR成像(三):快时间与慢时间
SAR成像(四):多普勒频移的计算
SAR成像(五):回波模型
SAR成像(六):距离徙动矫正
SAR成像(七):RD成像算法
参考:
[1]陈小利. 合成孔径雷达成像研究[D].南京大学,2015.
[2]https://blog.csdn.net/weixin_45858061/article/details/102986993
SAR成像(一):线性调频信号(LFM)和脉冲压缩的更多相关文章
- 信号基础知识----线性调频信号LFM //matlab命令:chirp
%关于线性调频信号(LFM) %参考书目:声呐技术,第二章P33 clc;close all;clear all;%参数----------------------------------f0=100 ...
- 【Matlab】线性调频信号LFM 仿真
[知识点] 生成序列 i = a:step:b 举例: i = 1:1:9 画图(子图) subplot(m,n,p)或者subplot(m n p) 总结起来就是,画一个m行n列的图. p表示在第p ...
- 【遥感专题系列】微波遥感(二、合成孔径雷达SAR基础)
目前使用最广的成像雷达系统就是合成孔径雷达(Synthetic Aperture Radar:SAR),SAR几乎成为了雷达的代名词.本文从应用角度介绍SAR系统的基本知识. 本文主要包括: SAR基 ...
- 【遥感专题系列】微波遥感(三、SAR图像特征)
SAR是主动式侧视雷达系统,且成像几何属于斜距投影类型.因此SAR图像与光学图像在成像机理.几何特征.辐射特征等方面都有较大的区别.在进行SAR图像处理和应用前,需要了解SAR图像的基本特征. 本文主 ...
- Computer Vision_33_SIFT:SAR-SIFT: A SIFT-LIKE ALGORITHM FOR SAR IMAGES——2015
此部分是计算机视觉部分,主要侧重在底层特征提取,视频分析,跟踪,目标检测和识别方面等方面.对于自己不太熟悉的领域比如摄像机标定和立体视觉,仅仅列出上google上引用次数比较多的文献.有一些刚刚出版的 ...
- HackRF实现ADS-B飞机信号跟踪定位
硬件平台:HackRF One软件平台:MAC运行环境搭建系统平台:OS X 10.11 EI Capitan文章特点:捕捉程序支持HackRF One且基于MAC平台验证通过有效. 1. 原理概述 ...
- 合成孔径雷达(Synthetic Aperture Radar, SAR)
合成孔径雷达(Synthetic Aperture Radar, SAR)是20世纪最先进的科技发明之一.SAR有机载与卫载系统之分,能提供地表地形.乃至行星等远距离目标区高解析度图像.目前无人飞行载 ...
- 数字麦克风PDM信号采集与STM32 I2S接口应用(二)
在使用STM32的数字麦克风I2S接口时,计算采样率让人头疼,芯片手册上没有明确的说法,而手册上的计算方法经过测试确和实验不符.借助搜索引擎,大部分资料都是来自于开发板卖家或开发板论坛,主要是咪头采集 ...
- linux系统编程之信号(五)
今天继续对信号进行学习,开始正入正题: sigaction函数: 安装信号之前我们已经学过一个函数:signal,它最早是在unix上出现的,它是对不可靠信号进行安装的,之后出现了可靠信号和实时信号, ...
- fMRI数据分析处理原理及方法(转载)
原文地址:http://www.cnblogs.com/minks/p/4889497.html 近年来,血氧水平依赖性磁共振脑功能成像(Blood oxygenation level-depende ...
随机推荐
- 【实战】yolov8 tensorrt模型加速部署
[实战]yolov8 tensorrt模型加速部署 TensorRT-Alpha基于tensorrt+cuda c++实现模型end2end的gpu加速,支持win10.linux,在2023年已经更 ...
- 解决node.js报错Invalid character in header content ["Content-Disposition"]
遇到这种报错一般在于下载文件时候,如果Content-Disposition设置文件名有中文会出现此种问题,解决方案如下: 把第二段代码改为第一段,即可~~ 'Content-Disposition' ...
- Map接口中的常用方法-Map集合遍历建找值方式
Map接口中的常用方法 Map接口中定义了很多方法,常用的如下: public v put(K key,v value):把指定的键与指定的值添加到Map集合中. public v remove(0b ...
- ssh 连接 wsl2
修改ssh相关config https://cloud.tencent.com/developer/article/1538305 其中,端口可以设置成 2222 以防万一 然后在wsl 中执行如下命 ...
- mysql02-配置文件my.ini
https://www.cnblogs.com/isme-zjh/p/11542106.html 1.mysql安装-连接 1.1安装 略 1.2连接-退出 连接 mysql –h 主机 ...
- 2.1.新建项目及项目目录和预览uni项目
目录结构 一个uni-app工程,默认包含如下目录及文件 static目录 使用注意 编译到任意平台时,static 目录下除不满足条件编译的文件,会直接复制到最终的打包目录,不会打包编译.非 sta ...
- .NET Core项目部署到Windows系统Docker
1.新建一个ASP.NET Core Web应用(模型-视图-控制器) 2. 项目启动Docker.Docker OS选择Windows 3.系统默认的dockerfile文件如下 #See http ...
- 3D建模零代码平台
近几年,随着国内外文化产业的迅猛发展,3D建模行业迎来黄金发展期. 尤其是在元宇宙时代及数字体验经济时代的大背景下,越来越多的实时.可交互的3D内容将出现在人们的生活中. 关于3D建模师而言,无疑,行 ...
- CSS3--点击弹出菜单动画
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <meta name ...
- Cesium之Quick Start
1. 引言 Cesium是一款三维地球和地图可视化开源JavaScript库,使用WebGL来进行硬件加速图形,使用时不需要任何插件支持,基于Apache2.0许可的开源程序,可以免费用于商业和非商业 ...