CN112073873A - 一种无冗余阵元的一阶可调差分阵列的优化设计方法 - Google Patents
一种无冗余阵元的一阶可调差分阵列的优化设计方法 Download PDFInfo
- Publication number
- CN112073873A CN112073873A CN202010824086.5A CN202010824086A CN112073873A CN 112073873 A CN112073873 A CN 112073873A CN 202010824086 A CN202010824086 A CN 202010824086A CN 112073873 A CN112073873 A CN 112073873A
- Authority
- CN
- China
- Prior art keywords
- array
- angle
- order
- gamma
- theta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000013461 design Methods 0.000 title claims abstract description 15
- 230000004044 response Effects 0.000 claims description 29
- 230000005404 monopole Effects 0.000 claims description 10
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 230000001629 suppression Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/40—Details of arrangements for obtaining desired directional characteristic by combining a number of identical transducers covered by H04R1/40 but not provided for in any of its subgroups
- H04R2201/401—2D or 3D arrays of transducers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2430/00—Signal processing covered by H04R, not provided for in its groups
- H04R2430/20—Processing of the output signals of the acoustic transducers of an array for obtaining a desired directivity characteristic
- H04R2430/21—Direction finding using differential microphone array [DMA]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2430/00—Signal processing covered by H04R, not provided for in its groups
- H04R2430/20—Processing of the output signals of the acoustic transducers of an array for obtaining a desired directivity characteristic
- H04R2430/23—Direction finding using a sum-delay beam-former
Landscapes
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明公开了一种无冗余阵元的一阶可调差分阵列的优化设计方法,利用二阶泰勒近似来设计一阶可调差分阵列,先构造最小二乘意义下实现差分阵列波束频率不变性的目标残差代价函数,再分析残差函数极值所对应调向角度和阵列夹角,最后根据给定的优化调向空间确定最优阵型。本发明通过优化阵型来提升无冗余阵元一阶可调差分阵列的鲁邦性能,与常规均匀等边三角形阵列相比,在优化调向空间内,保证频率不变性能的同时提升白噪声抑制能力。
Description
技术领域
本发明属于麦克风阵列领域,特别涉及了一种阵型优化设计方法。
背景技术
麦克风差分阵列是一种具有小尺寸,超指向性,频率不变特性的传声器阵列,被广泛应用于各种便携可持设备和移动通信终端中。
在一些实际应用中,声源仅在一定的角度范围移动,如鲁棒可调波束形成器(见文献1:C.C.Lai,S.E.Nordholm and Y.H.Leung,Design of robust steerable broadbandbeamformers incorporating microphone gain and phase error characteristics,"inProc.ICASSP.578(2011),pp.101-104.)、可调波束形成器(见文献2:C.C.Lai,S.E.Nordholm and Y.H.Leung,A Study into the Design of Steerable MicrophoneArrays(Springer,Australia,2017).),而差分阵列设计中还没有无冗余阵元数目实现约束空间调向的相关研究方法。
在差分阵列设计中,最基本一阶可调差分阵列的构成需要至少3个非共线阵元,如3阵元等边三角形阵型、等腰直角三角形阵型(见文献3:G.W.Elko,“Steerable andvariable first-order differential microphone array,”United States Patent,No.:6041127,Mar.21,2000.)、3阵元等腰直角三角形阵型(见文献4:B.De Schuymer andH.Brouckxon,Steerable microphone array system with a first order directionalpattern,European Patent Application,Application No.:10151106.1,2011-8-24.)。文献3和4均针对360°全空间调向,利用特殊阵型结构,根据传统一阶泰勒近似进行分析和构造可调差分阵列,这些方法在微小尺寸或非常低的频率条件下,近似响应和实际响应不会出现较大误差,但对于非特殊结构阵型(非正交或非等边阵型)或高频条件,传统方法近似分析精度将不能满足要求,因为当阵列夹角或主瓣指向发生变化时,由二阶以上近似项产生的阵列响应虚部也随之变化且不能忽略不计,从而影响阵列响应形成不同的波束。
发明内容
为了解决上述背景技术提到的技术问题,本发明提出了一种无冗余阵元的一阶可调差分阵列的优化设计方法,保证了连续调向优化空间内波束的频率不变性,提高了阵列鲁棒性,减小了对实际阵列分析的近似误差。
为了实现上述技术目的,本发明的技术方案为:
一种无冗余阵元的一阶可调差分阵列的优化设计方法,包括以下步骤:
其中,θ为方位角,φ为俯仰角,E=[E1,E2,E3]T为导向矢量,上标T表示转置,对应三个阵元的阵列响应分别为E1=1,E2=exp(2jΩsinφcosθ),E3=exp[2jΩsinφcos(θ-γ)],为虚数单位,γ为阵列夹角,Ω=ωd/(2c),ω=2πf,f为频率,d为相邻阵元的距离,c为声音在空气中的传播速度;
W表示三个阵元的加权因子:
W=αWm+(1-α)Wd
其中,Wm为归一化单极子特征波束加权因子,Wd为偶极子加权因子:
Wm=[1,0,0]T
其中,θs为调向角度;
(2)构造最小二乘意义下残差代价函数:
针对J(γ,θs)分析不同阵列夹角γ条件下,残差代价函数极值对应的调向角度θs;
进一步地,阵列夹角γ≠60°,即阵列为非等边三角形阵型分布。
进一步地,在步骤(2)中,根据阵列响应和加权因子计算阵列性能指标——白噪声增益和指向性因数:
其中,WNG(γ,θs)为白噪声增益,DF(γ,θs)为指向性因数。
采用上述技术方案带来的有益效果:
本发明能够通过优化阵型来提升无冗余阵元一阶可调差分阵列的鲁邦性能,与常规均匀等边三角形阵列相比,在优化调向空间内,保证频率不变性能的同时提升白噪声抑制能力。
附图说明
图1是本发明中一阶可调差分阵列坐标图;
图2是本发明中一阶可调差分阵列实现框图;
图3是实例1中一阶、二阶泰勒近似响应与实际阵列响应对比图;
图4是实例2中半空间优化残差代价函数J(γ,θs)图;
图5是实例2中γ=240°阵列与等边三角形阵列差分波束对比图;
图6是实例2中半空间优化白噪声增益WNG(γ,θs)图;
图7是实例2中半空间优化指向性因数DF(γ,θs)图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明设计了一种无冗余阵元的一阶可调差分阵列的优化设计方法,
其阵型结构如图1所示,M2-M1和M3-M1麦克风对夹角为γ,阵元M1和M2(或M1和M3)间的距离为d,在该阵型结构基础上,图2进一步给出了本发明优化设计的实际实现框图,具体过程如下:
1、采用短时傅里叶变换,根据二阶泰勒近似在频域上先分别设计单极子和偶极子,然后根据波束图调节参数α合成差分阵列响应。
1)单极子由3阵元加权组合而成,其阵列响应:
其中,(g)T表示转置,θ为方位角,φ为俯仰角,E=[E1,E2,E3]T为导向矢量,对应三个阵元的阵列响应分别为E1=1,E2=exp(2jΩsinφcosθ),E3=exp[2jΩsinφcos(θ-γ)],式中为虚数单位,γ为阵列夹角,Ω=ωd/(2c),ω=2πf,f为频率,d为阵元M1和M2(或M1和M3)间的距离,c为声音在空气中的传播速度。
单极子加权因子Wm=[w1,w2,1-w1-w2]T,单极子无指向性,对任意空间角度幅度响应均为1,根据二阶泰勒近似合成归一化单极子响应波束,即选用原点处阵元1来构造单极子,即Wm=[1,0,0]T。
2)偶极子由全向麦克风对M2-M1、M3-M1两两相减构成,同样利用二阶近似,0和γ两个方向:
对上面两个非正交偶极子进行加权,并做归一化处理,可得调向偶极子响应:
3)根据1)和2)推出合成阵列响应:
其中W=αWm+(1-α)Wd为合成阵列加权因子,α为波束图调节参数。
2、利用最小二乘方法构造残差代价函数。
由J(γ,θs)残差代价函数可知同一俯仰角φ下残差J(γ,θs)仅由夹角γ和调向角度θs决定。此外,根据1中阵列响应和加权因子可推出阵列性能指标白噪声增益和指向性因数:
下面针对残差代价函数进行分析:
1)对目标残差代价函数求导:
推出极值点θs=γ/2+k·90°(k=0,1,2,3)或γ=60°。其中,阵列夹角γ=60°对应残差代价函数值J(γ,θs)|γ=60°≈3(1-α)2Ω2sin4φ/8,在θs∈[0,360°]调向空间残差值基本不发生变化。
调向空间起始角度可根据实际设计要求,结合函数对称性对阵列做相应的空间旋转。本发明中所述的一种无冗余阵元的一阶可调差分阵列,指的是非共线3阵元可调差分阵列,仿真实例均与等边正三角形对比,下面结合仿真实例对发明进行详细说明。
实例1
对阵列响应做二阶泰勒近似并且分析其虚部是本发明的关键。为此,举例说明二阶泰勒近似产生的虚部对阵列响应的影响。图3为Ω=π/18rad、阵列夹角γ=90°等腰三角形阵列分别利用传统一阶泰勒近似方法、本发明二阶泰勒近似方法与实际阵列响应对比波束图。图3中的(a)、(b)为超心形指向hypercardioid,对应一阶和二阶方向图最大近似误差分别为-19.1197dB和-37.9761dB;图3中(c)、(d)为心形指向cardioid方向图,对应一阶和二阶方向图最大近似误差分别为-22.6416dB和-41.4979dB。从图中可以看出,调向角度θs不同,实际阵列响应波束零陷深度不同,波束发生扭曲,若采用传统一阶泰勒近似实际响应产生的误差较大,分析结果将出现错误,而本发明方法保留了含有阵列夹角γ和调向角度θs产生的虚部分量,能够很好地近似实际阵列响应,便于更加精确的分析阵列实际响应和优化设计。
实例2
图4为本发明方法设计的阵列残差代价函数在不同优化夹角γ和不同调向角度θs条件下的变化,通过仿真观察得到在优化空间内 优化后的阵型残差代价函数相对于均匀阵型其数值更小,以图4中夹角γ=240°为例,做进一步说明,图5给出了本发明方法设计的阵列夹角γ=240°差分阵列与均匀等边三角形阵列差分波束对比图,取Ω=π/18rad,根据理论分析推出夹角γ=240°对应最优调向范围为θs∈[0,60°],图5中方向图调向角度分别为(a)θs=10°,(b)θs=20°,(c)θs=30°,(d)θs=40°,该范围内残差代价函数值更小,对应零陷深度为-24.9236dB,-31.4479dB,-42.4049dB和-31.4503dB,而当γ=60°时,对应心形指向方向图零陷深度为-22.0959dB,-22.4865dB,-22.4976dB和-22.4866dB,优化后心形波束具有更深的零陷。
在阵列鲁棒性提升方面,从图6和图7可以看出,不同阵型,即夹角γ=220°,240°,260°,280°,相对于均匀等边三角形阵型,在对应优化调向空间内,差分波束在保证阵列指向性基本不发生变化的前提下(最小值偏差0.0391dB),最大白噪声增益分别可提升4.8323,4.2550,3.3524和2.0078dB,从而有效提升阵列的鲁棒性。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (5)
1.一种无冗余阵元的一阶可调差分阵列的优化设计方法,其特征在于,包括以下步骤:
其中,θ为方位角,φ为俯仰角,E=[E1,E2,E3]T为导向矢量,上标T表示转置,对应三个阵元的阵列响应分别为E1=1,E2=exp(2jΩsinφcosθ),E3=exp[2jΩsinφcos(θ-γ)],为虚数单位,γ为阵列夹角,Ω=ωd/(2c),ω=2πf,f为频率,d为相邻阵元的距离,c为声音在空气中的传播速度;
W表示三个阵元的加权因子:
W=αWm+(1-α)Wd
其中,Wm为归一化单极子特征波束加权因子,Wd为偶极子加权因子:
Wm=[1,0,0]T
其中,θs为调向角度;
(2)构造最小二乘意义下残差代价函数:
针对J(γ,θs)分析不同阵列夹角γ条件下,残差代价函数极值对应的调向角度θs;
2.根据权利要求1所述无冗余阵元的一阶可调差分阵列的优化设计方法,其特征在于,阵列夹角γ≠60°,即阵列为非等边三角形阵型分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010824086.5A CN112073873B (zh) | 2020-08-17 | 2020-08-17 | 一种无冗余阵元的一阶可调差分阵列的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010824086.5A CN112073873B (zh) | 2020-08-17 | 2020-08-17 | 一种无冗余阵元的一阶可调差分阵列的优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112073873A true CN112073873A (zh) | 2020-12-11 |
CN112073873B CN112073873B (zh) | 2021-08-10 |
Family
ID=73661902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010824086.5A Active CN112073873B (zh) | 2020-08-17 | 2020-08-17 | 一种无冗余阵元的一阶可调差分阵列的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112073873B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022170541A1 (en) * | 2021-02-10 | 2022-08-18 | Northwestern Polytechnical University | First-order differential microphone array with steerable beamformer |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6041127A (en) * | 1997-04-03 | 2000-03-21 | Lucent Technologies Inc. | Steerable and variable first-order differential microphone array |
US20040001598A1 (en) * | 2002-06-05 | 2004-01-01 | Balan Radu Victor | System and method for adaptive multi-sensor arrays |
CN104464739A (zh) * | 2013-09-18 | 2015-03-25 | 华为技术有限公司 | 音频信号处理方法及装置、差分波束形成方法及装置 |
US20160134969A1 (en) * | 2012-12-04 | 2016-05-12 | Jingdong Chen | Low noise differential microphone arrays |
CN110133579A (zh) * | 2019-04-11 | 2019-08-16 | 南京航空航天大学 | 适用于球面麦克风阵列声源定向的球谐波阶数自适应选择方法 |
CN110244286A (zh) * | 2019-07-09 | 2019-09-17 | 西北工业大学 | 一种高增益无左右舷模糊的阵列设计方法 |
CN110415720A (zh) * | 2019-07-11 | 2019-11-05 | 湖北工业大学 | 四元差分麦克风阵列超指向性频率不变波束形成方法 |
CN110457796A (zh) * | 2019-07-29 | 2019-11-15 | 南京大学 | 无规平面差分传声器阵列阵元位置的优化方法 |
WO2020059977A1 (ko) * | 2018-09-21 | 2020-03-26 | 엘지전자 주식회사 | 연속적으로 스티어링 가능한 2차 디퍼런셜 마이크로폰 어레이 및 그것을 구성하는 방법 |
-
2020
- 2020-08-17 CN CN202010824086.5A patent/CN112073873B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6041127A (en) * | 1997-04-03 | 2000-03-21 | Lucent Technologies Inc. | Steerable and variable first-order differential microphone array |
US20040001598A1 (en) * | 2002-06-05 | 2004-01-01 | Balan Radu Victor | System and method for adaptive multi-sensor arrays |
US20160134969A1 (en) * | 2012-12-04 | 2016-05-12 | Jingdong Chen | Low noise differential microphone arrays |
CN104464739A (zh) * | 2013-09-18 | 2015-03-25 | 华为技术有限公司 | 音频信号处理方法及装置、差分波束形成方法及装置 |
WO2020059977A1 (ko) * | 2018-09-21 | 2020-03-26 | 엘지전자 주식회사 | 연속적으로 스티어링 가능한 2차 디퍼런셜 마이크로폰 어레이 및 그것을 구성하는 방법 |
CN110133579A (zh) * | 2019-04-11 | 2019-08-16 | 南京航空航天大学 | 适用于球面麦克风阵列声源定向的球谐波阶数自适应选择方法 |
CN110244286A (zh) * | 2019-07-09 | 2019-09-17 | 西北工业大学 | 一种高增益无左右舷模糊的阵列设计方法 |
CN110415720A (zh) * | 2019-07-11 | 2019-11-05 | 湖北工业大学 | 四元差分麦克风阵列超指向性频率不变波束形成方法 |
CN110457796A (zh) * | 2019-07-29 | 2019-11-15 | 南京大学 | 无规平面差分传声器阵列阵元位置的优化方法 |
Non-Patent Citations (3)
Title |
---|
2017 25TH EUROPEAN SIGNAL PROCESSING CONFERENCE (EUSIPCO): "Design and analysis of second-order steerable differential microphone arrays", 《2017 25TH EUROPEAN SIGNAL PROCESSING CONFERENCE (EUSIPCO)》 * |
XIAOGUANG WU;HUAWEI CHEN;JIANJIANG ZHOU;TIANWEN GUO: "Study of the Mainlobe Misorientation of the First-Order Steerable Differential Array in the Presence of Microphone Gain and Phase Errors", 《IEEE SIGNAL PROCESSING LETTERS》 * |
白雅瑞: "差分麦克风阵列波束形成中的零陷展宽技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022170541A1 (en) * | 2021-02-10 | 2022-08-18 | Northwestern Polytechnical University | First-order differential microphone array with steerable beamformer |
Also Published As
Publication number | Publication date |
---|---|
CN112073873B (zh) | 2021-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1429581B1 (en) | Method of broadband constant directivity beamforming for non linear and non axi-symmetric sensor arrays embedded in a obstacle | |
Zhang et al. | Pattern synthesis for arbitrary arrays via weight vector orthogonal decomposition | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN105182302A (zh) | 一种抗快速运动干扰稳健零陷展宽波束形成方法 | |
Yang et al. | On the design of flexible Kronecker product beamformers with linear microphone arrays | |
CN115038012B (zh) | 基于admm的麦克风阵列鲁棒频率不变波束形成方法 | |
CN103020363A (zh) | 一种通过改善阵列波束方向图旁瓣特性设计天线的方法 | |
Xu et al. | A modified differential beamforming and its application for DOA estimation of low frequency underwater signal | |
CN112073873B (zh) | 一种无冗余阵元的一阶可调差分阵列的优化设计方法 | |
CN107170441B (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
JP2002118414A (ja) | アレーアンテナの制御装置及び制御方法 | |
Cheng et al. | An analytical solution for weighted least-squares beampattern synthesis using adaptive array theory | |
Nnonyelu et al. | Cardioid microphones/hydrophones in a collocated and orthogonal triad—A steerable beamformer with no beam-pointing error | |
CN111817766B (zh) | 联合阵列天线单元方向图的波束形成方法 | |
Haneda et al. | Design of multipole loudspeaker array based on spherical harmonic expansion | |
CN115276746B (zh) | 基于交替方向惩罚的频率一致宽带波束形成方法及系统 | |
Luo et al. | Constrained maximum directivity beamformers based on uniform linear acoustic vector sensor arrays | |
Liu et al. | Constant beamwidth receiving beamforming based on template matching | |
CN115358075A (zh) | 一种阵列天线宽波束增益优化方法 | |
Huang et al. | Properties and limits of the minimum-norm differential beamformers with circular microphone arrays | |
Li et al. | A null broadening beamforming method of virtual antenna array | |
Li et al. | A novel algorithm for null broadending beamforming based on subspace projection and virtual antenna array | |
Jin et al. | Design of optimal linear differential microphone arrays based array geometry optimization | |
Feng et al. | Design of Frequency-Invariant Beamformers Using a Weighted Spatial-Response-Variance Measure | |
Vescovo | Beam scanning with null and excitation constraints for linear arrays of antennas |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |