Delay and sum beamforming

作者:桂。

时间:2017-06-03  15:40:33

链接:http://www.cnblogs.com/xingshansi/p/6937576.html



前言

本文主要记录麦克风阵列的几个基本知识点,并简单介绍基本的波束合成方法:Delay and sum (DSB).

一、栅瓣效应

类似干涉仪中的phase-wrapping问题:传感器的距离超过半波长的距离。

  A-Beam Pattern

以线性阵列为例

输出为

对应的增益,也就是output,取频率为1KHz的数据:

有时候也用极坐标表示

  B-Frequency Response

考虑如下场景

实际只有一个波峰,却在不同角度长出多个一样的,这就是栅瓣

相位是角度、频率的共同作用,对应的栅瓣从二维平面更容易理解

从每一个来看,grating lobe就是栅瓣:

二、Delay and sum

DSB核心的思想就是延迟相加,对于栅瓣的DOA估计需要解模糊。假设估计出DOA,进行补偿之后(delay)通过叠加(sum)就可以实现目标信号的增强:

关于delay的操作,可以在时域进行也可以在频域进行,本质上是小数延迟的问题(fractional delays).例如对于下面的dual-channel问题:

需要延迟 是0.5149ms,也就是24.717samples,而不是整数。

小数延迟实现的思路有很多种,例如可以通过线性相位(延迟)的FIR构造,对应就是sinc函数,通常构造完成之后会做一个加窗的平滑处理:

double delay = 0.25;               // Fractional delay amount
int filterLength = 11;             // Number of FIR filter taps (should be odd)
int centreTap = filterLength / 2;  // Position of centre FIR tap

for (int t=0 ; t<filterLength ; t++)
{
   // Calculated shifted x position
   double x = t - delay;

   // Calculate sinc function value
   double sinc = Math.sin(Math.PI * (x-centreTap)) / (Math.PI * (x-centreTap));

   // Calculate (Hamming) windowing function value
   double window = 0.54 - 0.46 * Math.cos(2.0 * Math.PI * (x+0.5) / filterLength);

   // Calculate tap weight
   double tapWeight = window * sinc;

   // Output information
   System.out.printf("%3d % f % f % f\n", t, sinc, window, tapWeight);
}

  其他如:重采样技术(Resample)、Farrow-filter技术,都可以实现这样的功能。

参考

  • http://www.labbookpages.co.uk/audio/beamforming/delaySum.html#pattern
时间: 2024-10-07 07:27:22

Delay and sum beamforming的相关文章

Pascal小游戏 打飞机

一个经典的打飞机游戏(1)Pascal代码 十分经典,有一种街机的感觉 奇葩青年的又一控制台神作. uses crt; type list=record         ty,ax:integer;         end;      xy=record         bx,by:integer;         end;      l1=array[1..4,1..5]of char;      l2=array[1..5,1..6]of char;      l3=array[1..8,1

翻译:使用红外传感器与Arduino进行简单动作与手势检测

译注:昨天看 Adruino 的 Twitter 推了这篇项目,第一眼就觉得非常有趣,翻译给大家看看.文中的红外传感器比较高级,和淘宝上5块钱的那种只能输出0和1的不一样, TPA81 是可以输出温度的,还是8个连续点. MLX90614 可以输出一点的温度,还卖将近三十,可以想象 TPA81 的价格--而且淘宝上并没有卖 TPA81 的,可见国内还没有仿的.文章也给出了另一个解决方案,使用8个 HC-SR505 来替代,这就便宜了很多.我想了一下,如果要只是实现像文章那样的增减手势,3个 HC

Verilog之delay的两种用法(inter/intra)

verilog语言中有两种延迟方式:inter-delay和intra-delay,关于inter和intra.这两个英文前缀都有“内部,之间”的意思,但又有所不同.inter表达不同事物之间,intra表达同类事物之间,两者具体的含义请细细体会:).以阻塞式赋值为例(block assignment): 1.inter-delay的表达式为:  #delay-value a=b+c: 先说说阻塞式赋值语句执行的一般过程: block assignment方程的RHS先估值(evaluation

复习图---Delay Constrained Maximum Capacity Path(SPFA+二分)

Delay Constrained Maximum Capacity Path Time Limit:10000MS     Memory Limit:65535KB     64bit IO Format:%I64d & %I64u Submit Status Description Consider an undirected graph with N vertices, numbered from 1 to N, and M edges. The vertex numbered with

【HDU 1839】 Delay Constrained Maximum Capacity Path(二分+最短路)

[HDU 1839] Delay Constrained Maximum Capacity Path(二分+最短路) Delay Constrained Maximum Capacity Path Time Limit: 10000/10000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Others) Total Submission(s): 1515    Accepted Submission(s): 481 Problem

hdu 1839 Delay Constrained Maximum Capacity Path(spfa+二分)

Delay Constrained Maximum Capacity Path Time Limit: 10000/10000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Others)Total Submission(s): 1790    Accepted Submission(s): 577 Problem Description Consider an undirected graph with N vertices, nu

LeetCode OJ - Sum Root to Leaf Numbers

这道题也很简单,只要把二叉树按照宽度优先的策略遍历一遍,就可以解决问题,采用递归方法越是简单. 下面是AC代码: 1 /** 2 * Sum Root to Leaf Numbers 3 * 采用递归的方法,宽度遍历 4 */ 5 int result=0; 6 public int sumNumbers(TreeNode root){ 7 8 bFSearch(root,0); 9 return result; 10 } 11 private void bFSearch(TreeNode ro

129. Sum Root to Leaf Numbers

Given a binary tree containing digits from 0-9 only, each root-to-leaf path could represent a number. An example is the root-to-leaf path 1->2->3 which represents the number 123. Find the total sum of all root-to-leaf numbers. For example, 1 / 2 3 T

Leetcode 494 Target Sum 动态规划 背包+滚动数据

这是一道水题,作为没有货的水货楼主如是说. 题意:已知一个数组nums {a1,a2,a3,.....,an}(其中0<ai <=1000(1<=k<=n, n<=20))和一个数S c1a1c2a2c3a3......cnan = S, 其中ci(1<=i<=n)可以在加号和减号之中任选. 求有多少种{c1,c2,c3,...,cn}的排列能使上述等式成立. 例如: 输入:nums is [1, 1, 1, 1, 1], S is 3. 输出 : 5符合要求5种