时域、空域对偶性

厚着脸皮要在同事公众号上写篇文章,尽量浅显、与专业相关,选了这个主题。


 一、时域与空域特性

以远场模型(平面波)为例,假设均匀线阵接收的为窄带信号,假设相邻振元间隔为d,入射角为

从空域坐标来看,相邻振元的间隔为:

等价到时间轴来看,采样点的间距为:,对应时间间隔为:

 

二、时、空域与采样定理

  A-空域角度理解

相邻振元的相位差为:

以干涉仪为例,如果存在相位模糊,有

k为非零整数,如果希望不出现相位模糊

对应扫描边界,则有

容易证明,同干涉仪一样,均匀线阵谱估计中的导向矢量,如果不满足上面的约束条件,同样会有多峰的问题。

  B-时域角度理解

前文提到,采样点对应的时间间隔为,即采样周期。空域均匀线阵对应时域均匀采样,采样频率:

入射信号的频率为:

如果采样无混叠,需要满足Nyquist采样定理:

该约束条件等价于:

可以看出均匀线阵的相位无模糊对应时域均匀采样的奈奎斯特定理。多说一句,如果是非均匀线阵、圆阵等形式,可以理解成对应维度的非均匀采样;从空域角度理解,非均匀阵列可以解决模糊问题,从时域角度理解,稀疏采样/非均匀采样可以突破奈奎斯特采样定理。

 

三、时、空域及MVDR算法

波束形成主要对感兴趣的方向进行增强/抑制,而谱估计更多是参数估计问题,前者操作多为主动,后者操作多为被动,MVDR算法对二者均适用。这里暂且抛开应用场景,仅从时、空角度理解MVDR的等价性。

接着上文的时域、空域思路,这里先从时域的角度来表述,为了简化均不考虑加窗情形。

  A-时域角度理解

对于N点均匀采样的信号,对其进行傅里叶变换:

的相关函数为:

容易证明有如下对应关系:

而相关函数对应的傅里叶变换为功率谱密度,可以求解功率谱密度:

  B-空域角度理解

 N个均匀线阵接收单元,对应的波束形成为:

即空域的波束形成可以理解为时域的傅里叶变换,

从而空域的功率谱密度可以等价为:

考虑到时域、空域具有等价性,空域的功率谱这么理解是合理的。

现在以常用的MVDR算法来理解这种等价性(MUSIC等谱估计算法的谱与波束形成完全不同,这里仅从时、空等价性上讨论“谱”的概念):

接收信号:

MVDR就是含有等式约束的最优化问题:

可以求解:

这个时候,如果将最优的w带入y,空域角度理解:y对应就是波束形成的结果。时域角度理解:y对应为傅里叶变换的结果。

通常MVDR的结果为的输出,根据上文分析可知,该结果从时域理解就是功率谱(差一个常数),所以从空域角度称作“谱”其实也是可以被接受的,对应功率(谱):

因为这是在空域,为了与时域功率谱的名字加以区别,可以称其为空间谱。

具体空间谱名称怎么由来,本文并没有考证。本文只是提供了一种理解“空间谱”名称的角度,至少MUSIC等算法的“谱”便与此不同,或许MUSIC等算法只是继承了“空间谱”这个名词也未可知。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/242792.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

开机重启,用户登录注销

开机重启,用户登录注销 用户管理 查询用户信息 用户组 在Linux下没有消息就是成功了(没有消息就是好消息)! 用户和组相关文件

Power of Two

https://leetcode.com/problems/power-of-two/ Given an integer, write a function to determine if it is a power of two. 数字 2^n 是大于0的&#xff0c;而且等于1左移n位得到的数字&#xff0c;所以2^n与2^n-1 相与运算得到0. bool isPowerOfTwo(int n) {if(n < 0)…

DCASE2013挑战赛介绍

简介2013 年起&#xff0c;为了评测现有的环境声音检测方法&#xff0c;电子与电气工程师学会音频和声学信号处理协会(Institute of Electrical and Electronics Engineers Audio and Acoustic Signal Process, IEEE AASP )开始举办声学场景和事件的检测与分类挑战赛(Detection…

DCASE挑战赛原始提案文件(详细信息)

本文是根据DCASE2013挑战赛的提案文件&#xff0c;加上个人的理解做了相应的翻译&#xff0c;可能有不对的地方&#xff0c;在之后的会慢慢改善。 背景在过去的十年里&#xff0c;人们对在代码公布和公共评估中提出方法的语音和音频处理社区的兴趣越来越浓厚。公共评估可以作为…

Number of 1 Bits

https://leetcode.com/problems/number-of-1-bits/ Write a function that takes an unsigned integer and returns the number of ’1’ bits it has (also known as the Hamming weight). For example, the 32-bit integer ’11’ has binary representation 0000000000000…

信号分析方法介绍

从开始的ICA&#xff0c;到稀疏表示&#xff0c;然后2001年发明NMF。 利用矩阵分解来解决实际问题的分析方法很多&#xff0c;如PCA(主成分分析)、ICA(独立成分分析)、SVD(奇异值分解)、VQ(矢量量化)等。在所有这些方法中&#xff0c;原始的大矩阵V被近似分解为低秩的VWH形式。…

文件目录类

文件目录类 enter表示查看下一行&#xff1b;空格键表示查看下一页

Counting Bits

https://leetcode.com/problems/counting-bits/ Given a non negative integer number num. For every numbers i in the range 0 ≤ i ≤ num calculate the number of 1’s in their binary representation and return them as an array. Example: For num 5 you should …

Detection and Classification of Acoustic Scenes and Events(DCASE2013详细介绍)

在DCASE 2013官网上&#xff0c;了解到在本篇文章中&#xff0c;详细描述了DCASE2013挑战赛的结果。下面加上个人的理解做了相应的翻译&#xff0c;可能有不对的地方&#xff0c;在之后的会慢慢改善。 摘要对于智能系统来说&#xff0c;使用音频形态是最好的&#xff0c;重点在…

MarkDown编辑器常用语法

分割线三个或更多-_*&#xff0c;必须单独一行&#xff0c;可含空格&#xff0c;如---、 ***、 ___。 效果如下&#xff1a; 如果上面的符号是紧挨着上文的文字&#xff0c;就是细线&#xff0c;否则就是上面显示的粗线。 图片链接![图片的名称](图片的本地连接或网络连接)&…

时间日期类和搜索查找类

时间日期类 搜索查找类

基于智能计算的降维技术研究与应用

《基于智能计算的降维技术研究与应用》皋军著&#xff08;2013&#xff09; 前言特征降维的需求随着社会信息化的发展&#xff0c;在具体的只能识别过程中需要处理的数据越来越多地呈现出高位特征&#xff0c;比如图像处理、文本分类、视频检索、计算机视觉、微阵列数据基因选择…

Add Binary

https://leetcode.com/problems/add-binary/ Given two binary strings, return their sum (also a binary string). For example, a “11” b “1” Return “100”. 计算过程类似Verilog的全加器。 char* addBinary(char* a, char* b) {int i;int l1, l2, l3, temp;…

SG仿真常用模块

workspace交互配合gateway in/out&#xff0c;实现信号仿真与workspace的互联。 滤波器可与FDATool同时使用&#xff0c;直接关联FDATool的参数&#xff0c;而不必输入FDATool的滤波器系数。 也可调用Workspace中的变量&#xff0c;假设滤波器系数设计并量化后&#xff0c;为My…

DCASE三次挑战赛概览

资料来源&#xff1a;http://www.cs.tut.fi/sgn/arg/dcase2017/dcase 所有对环境声音分类和检测方面感兴趣或研究工作的研究人员都应邀通过DCASE社区与同行分享想法、问题和意见。 这个研究领域正在迅速发展&#xff0c;在国际信号处理会议&#xff08;如ICASSP或EUSIPCO&#…

常用矩阵运算

作者&#xff1a;桂。 时间&#xff1a;2017-09-09 12:48:45 链接&#xff1a;http://www.cnblogs.com/xingshansi/p/7498175.html 一、复数相乘 可以表示为分块的形式&#xff1a; 二、范数 A-范数基本定义 p 0&#xff0c;0范数&#xff0c;对应非零元素个数&#xff1b; …

文件及目录实际案例

文件及目录实际案例 可以通过命名从jack用户切换到root用户。从低权限用户到高权限用户需要输入用户密码。 su - rootr(读权限)&#xff1a;代表是否能使用ls w(写权限)&#xff1a;代表在目录下创建和删除文件 x(执行权限)&#xff1a;代表可以进入到目录&#xff0c;例如cd

Power of Three

https://leetcode.com/problems/power-of-three/ Given an integer, write a function to determine if it is a power of three. Follow up: Could you do it without using any loop / recursion? 3的次方数没有显著的特点&#xff0c;最直接的方法就是不停地除以3&…