文章仅为个人理解,如有不妥之处欢迎指正。

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π1​dtd​[2π(f0​t+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)和脉冲压缩的更多相关文章

  1. 信号基础知识----线性调频信号LFM //matlab命令:chirp

    %关于线性调频信号(LFM) %参考书目:声呐技术,第二章P33 clc;close all;clear all;%参数----------------------------------f0=100 ...

  2. 【Matlab】线性调频信号LFM 仿真

    [知识点] 生成序列 i = a:step:b 举例: i = 1:1:9 画图(子图) subplot(m,n,p)或者subplot(m n p) 总结起来就是,画一个m行n列的图. p表示在第p ...

  3. 【遥感专题系列】微波遥感(二、合成孔径雷达SAR基础)

    目前使用最广的成像雷达系统就是合成孔径雷达(Synthetic Aperture Radar:SAR),SAR几乎成为了雷达的代名词.本文从应用角度介绍SAR系统的基本知识. 本文主要包括: SAR基 ...

  4. 【遥感专题系列】微波遥感(三、SAR图像特征)

    SAR是主动式侧视雷达系统,且成像几何属于斜距投影类型.因此SAR图像与光学图像在成像机理.几何特征.辐射特征等方面都有较大的区别.在进行SAR图像处理和应用前,需要了解SAR图像的基本特征. 本文主 ...

  5. Computer Vision_33_SIFT:SAR-SIFT: A SIFT-LIKE ALGORITHM FOR SAR IMAGES——2015

    此部分是计算机视觉部分,主要侧重在底层特征提取,视频分析,跟踪,目标检测和识别方面等方面.对于自己不太熟悉的领域比如摄像机标定和立体视觉,仅仅列出上google上引用次数比较多的文献.有一些刚刚出版的 ...

  6. HackRF实现ADS-B飞机信号跟踪定位

    硬件平台:HackRF One软件平台:MAC运行环境搭建系统平台:OS X 10.11 EI Capitan文章特点:捕捉程序支持HackRF One且基于MAC平台验证通过有效. 1. 原理概述 ...

  7. 合成孔径雷达(Synthetic Aperture Radar, SAR)

    合成孔径雷达(Synthetic Aperture Radar, SAR)是20世纪最先进的科技发明之一.SAR有机载与卫载系统之分,能提供地表地形.乃至行星等远距离目标区高解析度图像.目前无人飞行载 ...

  8. 数字麦克风PDM信号采集与STM32 I2S接口应用(二)

    在使用STM32的数字麦克风I2S接口时,计算采样率让人头疼,芯片手册上没有明确的说法,而手册上的计算方法经过测试确和实验不符.借助搜索引擎,大部分资料都是来自于开发板卖家或开发板论坛,主要是咪头采集 ...

  9. linux系统编程之信号(五)

    今天继续对信号进行学习,开始正入正题: sigaction函数: 安装信号之前我们已经学过一个函数:signal,它最早是在unix上出现的,它是对不可靠信号进行安装的,之后出现了可靠信号和实时信号, ...

  10. fMRI数据分析处理原理及方法(转载)

    原文地址:http://www.cnblogs.com/minks/p/4889497.html 近年来,血氧水平依赖性磁共振脑功能成像(Blood oxygenation level-depende ...

随机推荐

  1. 【实战】yolov8 tensorrt模型加速部署

    [实战]yolov8 tensorrt模型加速部署 TensorRT-Alpha基于tensorrt+cuda c++实现模型end2end的gpu加速,支持win10.linux,在2023年已经更 ...

  2. 解决node.js报错Invalid character in header content ["Content-Disposition"]

    遇到这种报错一般在于下载文件时候,如果Content-Disposition设置文件名有中文会出现此种问题,解决方案如下: 把第二段代码改为第一段,即可~~ 'Content-Disposition' ...

  3. Map接口中的常用方法-Map集合遍历建找值方式

    Map接口中的常用方法 Map接口中定义了很多方法,常用的如下: public v put(K key,v value):把指定的键与指定的值添加到Map集合中. public v remove(0b ...

  4. ssh 连接 wsl2

    修改ssh相关config https://cloud.tencent.com/developer/article/1538305 其中,端口可以设置成 2222 以防万一 然后在wsl 中执行如下命 ...

  5. mysql02-配置文件my.ini

     https://www.cnblogs.com/isme-zjh/p/11542106.html   1.mysql安装-连接   1.1安装 略   1.2连接-退出 连接 mysql –h 主机 ...

  6. 2.1.新建项目及项目目录和预览uni项目

    目录结构 一个uni-app工程,默认包含如下目录及文件 static目录 使用注意 编译到任意平台时,static 目录下除不满足条件编译的文件,会直接复制到最终的打包目录,不会打包编译.非 sta ...

  7. .NET Core项目部署到Windows系统Docker

    1.新建一个ASP.NET Core Web应用(模型-视图-控制器) 2. 项目启动Docker.Docker OS选择Windows 3.系统默认的dockerfile文件如下 #See http ...

  8. 3D建模零代码平台

    近几年,随着国内外文化产业的迅猛发展,3D建模行业迎来黄金发展期. 尤其是在元宇宙时代及数字体验经济时代的大背景下,越来越多的实时.可交互的3D内容将出现在人们的生活中. 关于3D建模师而言,无疑,行 ...

  9. CSS3--点击弹出菜单动画

    <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <meta name ...

  10. Cesium之Quick Start

    1. 引言 Cesium是一款三维地球和地图可视化开源JavaScript库,使用WebGL来进行硬件加速图形,使用时不需要任何插件支持,基于Apache2.0许可的开源程序,可以免费用于商业和非商业 ...