SAR成像基础知识急救箱(一)卷积 相关 滤波器那些事儿

1 相关与卷积

1.1 相关

自相关定义:

A(τ)=∫μ(t)μ(t+τ)dt

自相关的涵义是一个函数和平移后的自己的乘积的积分,注意自相关是平移量的函数。直观上理解:如果平移量为0,则对应上式的结果最大。对于如下式所示的信号:

μ 1 (t)={10  for ∣ t∣≤T/2otherwise

对于T=1  ,对应的图形和自相关函数分别为:

可以这样理解:积分代表面积,所以函数平移之后与原来的自己再相乘再积分,即是看乘积函数与水平轴围成的面积。对于上式,可以想象的是,若平移量较大,大于1,则乘积为零,面积为零,即相关函数值为零;同样的,若平移量为零,则在范围为(?0.5,0.5) 内乘积始终为1,面积为1,相关函数的值为1。

互相关是自相关的扩充,定义为,

C μv (τ)=∫μ ? (t)v(t+τ)dt

其中上标表示复共轭,对于实数,可以忽略。

几何上,相关解释为滑动内积,通常,相关是对等式右边两个信号的比较。

1.2 卷积

卷积的定义为:

f(τ)=∫f 1 (t)f 2 (τ?t)dt

对比上式,卷积与互相关的区别很明显,卷积是先将一个函数对折之后再平移,再与另一个函数相乘,再积分。它也是平移量的函数。显然若函数是是对称函数,则对折与否并不影响结果。

卷积有性质:线性和、结合律和交换律。

另外,还有二维卷积、离散卷积。对于离散卷积需要注意的是其输出结果存在边缘效应。

上面说到对于对称实信号,卷积和相关相同,但是对于下式来说,卷积和相关就体现出了差异,

μ 2 (t)={10 for 0<t<Totherwise

函数图形、自相关函数图形、卷积的图形分别为:

自相关函数图形和卷积函数图形的区别是“对折”造成的。

相关的物理意义还是比较好理解的,相关程度;卷积的物理意义是什么呢?牵强一点可以理解与一个函数和另一个函数对折之后的版本

的相关程度。应该不止如此吧。。。。。。

另外卷积是可交换的,相关不可:

C μv (t)=C ? μv (?t)

2 滤波器

2.1 滤波器的冲激响应和传递函数

因为:

f(t)=∫ +∞ ?∞ f(τ)δ(t?τ)dτ

现有线性时不变系统H ,且:

Hδ(t)=h(t)

则,for any input function x(t), the output function y(t) can be expressed as:

y(t) =Hx(t)=H∫ +∞ ?∞ x(τ)δ(t?τ)dτ=∫ +∞ ?∞ x(τ)h(t?τ)dτ=x(t)?h(t)

这是一个很有意思的结论,后面我们会用到。

若把滤波器看成是一个线性时不变系统,则h(t) 即为滤波器的冲激响应。

滤波器的传递函数定义为H(f) :

h(t)=∫ +∞ ?∞ H(f)e j2πft df

H(f)=∣ H(f)∣e j?(f)

其中:

  • ∣ H(f)∣ 表示滤波器的幅频特性,即当某一幅度为1,频率为f的信号通过滤波器后输出信号的幅度变化;
  • ?(f) 表示滤波器的相频特性,即当某一幅度为1,频率为f且相位为零的信号通过滤波器后输出信号的相位变化;

考虑到时域的卷积对应频域的乘积, 继续转换上式:

y(t) =x(t)?h(t)=∫ +∞ ?∞ x(τ)h(t?τ)dτ=∫ +∞ ?∞ X(f)H(f)e j2πft df

2.2 匹配滤波器

匹配滤波器在SAR成像中是非常关键的一个存在。直观的理解是:把回波信号的已知形式稍作转换后放到一个盒子里,当回波信号通过盒子时候,如果跟盒子里面的信号配上了,或者说比较像,盒子就尖叫一声,这样我们根据不同时间的尖叫声就可以知道对应不同距离的目标的回波信号啦。这里的数学原理是如下。

假设发射信号是μ(t) ,回波信号是r(t)=σμ(t?τ)

数学上描述两个信号是否配上,用相关:

y(t)=∫μ ? (s)r(s+t)ds=σ∫μ ? (s)μ(s+t?τ)ds=σA(t?τ)

显然,当t=τ ,输出值最大,注意这里的τ 为回波时延。如果有N个目标则输出中就有N个尖峰。若r(t)=∑ N i=1 σ i μ(t?τ i ) ,则对应的输出为:

y(t) =∫μ ? (s)r(s+t)ds=∑ i=1 N ∫μ ? (s)μ(s+t?τ i )ds=∑ i=1 N σ i A(t?τ i )

稍微作形式上的转换:

y(t)=∫μ ? (s)r(s+t)ds=∫r(s)μ ? (?(t?s))ds=∫r(s)h(t?s)ds

注意这是一个输入信号和h(t)=μ ? (?t) 的卷积。那么可想而知,这里的h(t)=μ ? (?t) 即为匹配滤波器的冲激响应。考虑到时域的卷积对应频域的乘积,若Y(f),R(f),H(f) 分别为y(t),r(t),h(t) 的傅里叶变换,则:

Y(f)=R(f)H(f)

其中:

H(f) =∫h(t)e ?j2πft dt=∫μ ? (?t)e ?j2πft dt=U ? (f)

注意这里f前并没有负号,可以这样理解:频域的复共轭对应时域的复共轭加上时间反折。

前面说过根据匹配滤波器的输出(尖峰的位置)我们可以进行目标定位,这并没有错。但匹配滤波器的输出尖峰的位置既可能对准目标回波信号输入的前沿也可能对准其零频位置,这跟匹配滤波器的具体实现方式有关。

值得注意的是以上的讨论都只考虑了连续信号的情况,实际处理中一般都是采样后得到的离散信号。这样带来的一个问题是,卷积会存在弃置区(因为边缘效应)。

时间: 2024-10-12 02:21:18

SAR成像基础知识急救箱(一)卷积 相关 滤波器那些事儿的相关文章

【搜索引擎基础知识3】搜索引擎相关开源项目及网站

部分内容转自:http://blog.csdn.net/hguisu/article/details/8024799 一. 开源项目 1.Lucene全文检索系统 http://lucene.apache.org和 http://www.lucene.com.cn/ Lucene是apache软件基金会4 jakarta项目组的一个子项目,是一个开放源代码的全文检索引擎工具包,即它不是一个完整的全文检索引擎,而是一个全文检索引擎的架构,提供了完整的查询引擎和索引引擎,部分文本分析引擎(英文与德文

计算机学科基础知识和主流编程语言相关内容的总结

这个仓库主要用来对计算机学科的基础理论知识和主流编程语言相关内容进行总结,每个目录对应一个大的方向. 主要内容 主要包含的计算机基础有以下几个大方向: DataBase DataStructure Linux_OS Network Algorithm Math 编程语言目前有: Python C++ 此外,还有几个目录: Coding:一些 C++ 的代码 OnlineJudge:一些 OJ 的解题代码 Offer:主要分享一些和找工作相关的内容 Others:一些无法归类,但是可能有用的东西

Oracle数据库基础知识_字符串操作相关2

6.LPAD,RPAD 作用:左/右边的字符串填充一些特定的字符语法: LPAD(string , n, [pad_String])          string:可是字符或者参数                 n:字符的长度,是返回的字符串的数量,如果这个数量比原字符串的长度要短,lpad函数将会把字符串截取成从左到右的n个字符; pad_string:是个可选参数,这个字符串是要粘贴到string的左边,如果这个参数未写,lpad函数将会在string的左边粘贴空格. 例如 SELECT

dubbo源码学习(一)基础知识及使用的相关技术

初学dubbo的源码,只做尝试性的去学习,做为自己学习的一个记录,各位看官如果觉得写的有错误或理解的不对,请在留言区告诉我,互相学习.本人能力有限,有大神进入 时请指点. Dubbo是Alibaba开源的分布式服务框架,它最大的特点是按照分层的方式来架构,使用这种方式可以使各个层之间解耦合(或者最大限度地松耦合),我们可以非常容易地通过Dubbo来构建分布式服务,并根据自己实际业务应用场景来选择合适的集群容错模式,这个对于很多应用都是迫切希望的,只需要通过简单的配置就能够实现分布式服务调用,也就

POSIX 消息队列基础知识复习,以及相关例程

1.1        Posix消息队列 1.1.1       消息队列的创建和删除 1.1.1.1     mq_open( ) #include<mqueue.h> mqd_tmq_open( const char *name, int flag ) mqd_t mq_open( const char *name, int flag, mode_t mode, mq_attrattr ) 创建或获取一个消息队列.成功返回消息队列描述符:失败返回-1. 参数name指定与消息队列相关联的名

基础知识:Web开发相关技术

Web应用需要遵循哪些规范 ? ? 主要包括目录结构,文件存放方式等方面的规范,比如Java Web程序的所有文件需要包含在一个文件夹中,这个文件夹就是我们的应用程序文件夹,在eclipse下就是建立一个工程的那个目录,应用程序文件夹下还需要包括一个WEB-INF文件夹,WEB-INF文件夹还需要包括Classes文件夹和lib文件夹,以及web.xml描述文件 ? ? 第三方的jar文件需要放在WEB-INF文件夹下的lib文件夹下 ? ? web.xml文件里定义了应用程序可以访问的资源和这

dubbo源码学习一:基础知识及使用的相关技术

Dubbo是Alibaba开源的分布式服务框架,它最大的特点是按照分层的方式来架构,使用这种方式可以使各个层之间解耦合(或者最大限度地松耦合),我们可以非常容易地通过Dubbo来构建分布式服务,并根据自己实际业务应用场景来选择合适的集群容错模式,这个对于很多应用都是迫切希望的,只需要通过简单的配置就能够实现分布式服务调用,也就是说服务提供方(Provider)发布的服务可以天然就是集群服务. Dubbo的产生背景.最初的需求.架构设计 等可以详细看官方的文档:http://dubbo.apach

图像增强相关基础知识

图像增强处理-1 图像增强是图像处理中一个重要的内容,在图像生成,传输或变换的过程中,由于多种因素的影响,造成图像质量下降,图像模糊,特征淹没,给分析和识别带来困难.因此,按特定的需要将图像中感兴趣的特征友选择地突出,衰减不需要的特征,提高图像的可懂度是图像增强的主要内容.图像增强不考虑图像降质的原因,而且改善后的图像也不一定逼近原图像,这是它与图像复原本质的区别.图像增强的主要目的有两个:一是改善图像的视觉效果,提高图像的清晰度:二是将图像转换成一种更适合人类或机器进行分析处理的形式,一边从图

linux内存基础知识和相关调优方案

内存是计算机中重要的部件之一,它是与CPU进行沟通的桥梁.计算机中所有程序的运行都是在内存中进行的,因此内存的性能对计算机的影响非常大.内存作用是用于暂时存放CPU中的运算数据,以及与硬盘等外部存储器交换的数据.只要计算机在运行中,CPU就会把需要运算的数据调到内存中进行运算,当运算完成后CPU再将结果传送出来,内存的运行也决定了计算机的稳定运行.对于整个操作系统来说,内存可能是最麻烦的的设备.而其性能的好坏直接影响着整个操作系统. 我们知道CPU是不能与硬盘打交道的,只有数据被载入到内存中才可