CN110720226A - 用于混合n个信息信号的方法和装置 - Google Patents
用于混合n个信息信号的方法和装置 Download PDFInfo
- Publication number
- CN110720226A CN110720226A CN201880037021.7A CN201880037021A CN110720226A CN 110720226 A CN110720226 A CN 110720226A CN 201880037021 A CN201880037021 A CN 201880037021A CN 110720226 A CN110720226 A CN 110720226A
- Authority
- CN
- China
- Prior art keywords
- components
- complex
- combined
- component
- spectral values
- 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
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
- H04R1/00—Details of transducers, loudspeakers or microphones
- H04R1/20—Arrangements for obtaining desired frequency or directional characteristics
- H04R1/32—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only
- H04R1/40—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers
- H04R1/406—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers microphones
-
- 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/03—Synergistic effects of band splitting and sub-band processing
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S2420/00—Techniques used stereophonic systems covered by H04S but not provided for in its groups
- H04S2420/07—Synergistic effects of band splitting and sub-band processing
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S5/00—Pseudo-stereo systems, e.g. in which additional channel signals are derived from monophonic signals by means of phase shifting, time delay or reverberation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
- H04S7/30—Control circuits for electronic adaptation of the sound field
Landscapes
- Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Complex Calculations (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
- Electrophonic Musical Instruments (AREA)
Abstract
本发明涉及一种用于混合N个信息时间信号(s1(t)、s2(t)、...)的方法和装置,信息时间信号分别首先在频域中转换成N个复数信息信号(v1(f,t1)、v2(f,t1)、...)之一,其中,N是大于1的整数。对此进行以下步骤。在N个复数信息信号的频率一致的频谱值中分别转换成第一和第二分量(208)。在频率一致的N个频谱值中的N个第一分量结合成第一组合分量(210)。在频率一致的N个频谱值中的N个第二分量结合成第二组合分量(212)。第一和第二组合分量结合成结果频谱值(214)。上述步骤也针对N个复数信息信号的频率一致的其他频谱值进行以便产生其他的结果频谱值(216、220)。由此获得的结果频谱值被结合成复数的输出信息信号(m(f,t1))。
Description
技术领域
本发明涉及一种用于混合N个信息时间信号的方法和装置,信息时间信号分别首先从时域转换成N个复数信息信号之一的频域,其中N是大于1的整数。这样的方法或装置例如用于对麦克风信号进行内插或外推。
背景技术
从EP 2994094B1中已知一种方法和装置,其中通过混合,麦克风信号从至少两个麦克风信号生成内插或外插的信号。
已知的方法涉及以下应用,其中,麦克风处于声场中,在该声场中,它们将各自麦克风位置处的声场测量变量(例如,声压)转换为麦克风信号,并且需要估计麦克风位置之外的声场测量变量的值,确切地说是从麦克风位置内插或外推的位置。
在已知方法中,内插或外推信号类似于在内插或外推位置处的声场测量变量。已知的方法使用能量相关的复数频谱值的加权以及加权的复数频谱值的总和,该总和包括补偿能量误差的校正。在已知方法中,通过校正,内插或外插信号具有在其平均能量上仅微不足道地偏离内插或外插位置处的声场测量变量的值的特性,并且即使通过多于一个声源由声波引起声场也保留了该特性。已知方法中的加权因子是从内插或外推的“虚拟”位置的计算表示中的系数得出的。
在已知方法中,内插或外推信号的相位与内插或外推位置处的声场测量变量的相位不匹配。在已知方法中从单个声源发出的直接声场的情况就是这种情况。在声场由多个声源的声波产生的情况下,在已知方法中,内插或外推信号在其相位上与内插或外推位置处的声场测量变量的偏差甚至更大。此外,在已知方法中,不可能进行超过麦克风的一倍距离以外的外推。提及的麦克风信号和内插或外推信号是复数值信号,其以通常的方式参照频率来描述变量的状态,在当前情况下是声场测量的变量。
通常将内插或外推位置计算为表示为矢量的位置的混合,特别是作为矢量的系数加权总和,其附加条件是系数总和等于1。由于附加条件,内插或外推的维数比位置数少1。例如,在这种情况下,在两个位置的情况下,通过该位置描述直线上的一维内插或外推位置;在三个位置的情况下,通过位置描述在平面上的二维内插或外推位置;或者在4个位置的情况下,在空间中进行三维内插或外推。
发明内容
本发明的目的是,提供系数作为控制参数。
需要指出的是,对于从单个声源出发的直接声场的情况得出在内插或外推位置处的关于声场测量变量的相位的有利结论,因为相位和空间位置存在物理上合法的关系,因此可以近似为假设平面波前的线性函数。
应当指出,在弥散声场的情况下,可得出在内插或外推位置处的关于声场变量的能量的有利结论,因为能量和空间位置在物理上是合法的关系,其在假设时间平均的情况下可以近似为常数。
在许多实际应用中,存在一个声场,该声场是来自多个声源的声波或直接声和扩散声叠加而成的。
附图说明
根据若干实施例在下面附图说明中详细阐述本发明。其中示出:
图1示出了根据本发明如何实现(N=)两个复数信息信号的混合;
图2示出了根据本发明的混合方法的流程图;
图3示出了根据第一实施例如何使(N=)两个复数信息信号的(N=)两个在频率上一致的矢量混合;
图4示出了根据第二实施例如何使(N=)2个复数信息信号的两个在频率上一致的矢量混合;
图5示出了用于混合(N=)两个信息信号的混合装置的实施例,该信息信号首先分别从时域转变到频域;
图6示出了用于混合(N=)三个信息信号的混合装置的实施例,该信息信号首先分别从时域转变到频域;
图7示出了从三个第一分量中推导出组合分量的实施例。
具体实施方式
首先根据图1详细阐述根据本发明的混合方法。对此,从两个信息信号、例如两个麦克风信号开始,麦克风信号彼此混合,例如用于麦克风信号的内插或外推。
然后,可以借助虚拟的麦克风的麦克风信号在内插或外推中均衡由混合物产生的结果信号,该虚拟麦克风的麦克风信号位于两个麦克风的连接线上的两个麦克风之间的位置。在外推时,此时可借助虚拟的麦克风的麦克风信号均衡结果信号,虚拟的麦克风的麦克风信号此时在两个麦克风的连接线上在两个麦克风之外。
两个麦克风信号作为时间的函数在图1中用s1(t)和s2(t)表示。首先借助从时域到频域的转换来转换这些信号。为此,将位于由W1表示的时间间隔中的时间信号转换为频域。此转换例如可以通过傅立叶变换进行。由此产生变换后的复数信息信号v1(f,t1)或v2(f,t1)作为频率f的函数。
然后,在图1中由附图标记100示意性地表示的混合方法中使两个变换后的复数信息信号的频率一致的复数频谱值v1(f1,t1)或v2(f1,t1)混合,以获得结果频谱值m(f1,t1)。然后,稍后详细阐述的该方法针对依次连续的频率一致的频谱值v1(f2,t1)或v2(f2,t1)进行重复。在图1中示意性地用附图标记101表示该重复的方法,且产生结果矢量m(f2,t1)。始终再次重复该混合方法,以便获得复数的输出信息信号m(f,t1)作为频率的函数。
此处需要提及的是,通过图1中的方框100和101示出的混合方法因此可依次重复进行,或可彼此同时并行地进行,从而在混合方法的控制的系统周期中可产生复数的输出信息信号m(f,t1)。
在例如借助逆傅立叶变换使复数的输出信息信号m(f,t1)从频域逆变换到时域中之后,然后在时间间隔W1中获得混合时间信号sc(t)。
然后可在接下来的时间间隔中重复此时描述的方法,例如在图1中通过W2示出。
图2在流程图中示出了,如何混合两个频率一致的复数频谱值。在开始该方法之后,在方框202中,首先在方框204中将N(N在该情况下等于二)个麦克风信号s1(t)和s2(t)从时域转换到频域中。由此产生N(=2)个转换的复数信息信号v1(f,t1)或v2(f,t1)。然后在方框206中选出N(=2)个复数信息信号的N(等于2)个频率一致的频谱值。这例如是图1中的第一频率值f1的频谱值v1(f1,t1)或v2(f1,t1)。
根据本发明此时在方框208中(在图1中也表示为步骤A)将(N=)两个复数频谱值转换成第一分量和第二分量。这此后根据图3a进一步阐述。在方框210中(也表示为步骤B),(N=)两个复数频谱值的第一分量结合成第一组合分量。这此后根据图3b进一步阐述。在下一方框212中(在图2中也表示为步骤C),(N=)两个复数频谱值的第二分量结合成第二组合分量。这此后根据图3c进一步阐述。然后在方框214中(在图2中也表示为步骤D),将第一组合分量和第二组合分量结合以便获得结果频谱值。这此后根据图3d进一步阐述。因此,以这种方式从两个频谱值v1(f1,t1)或v2(f1,t1)中推导出结果频谱值m(f1,t1)。此时针对下一频率值f2的下一(N=)两个频率一致的频谱值v1(f2,t1)或v2(f2,t1)重复该方法。这在图2中通过方框216和222来显示。在方框216,确定是否所有频谱值已经通过该方法处理。在方框222,选出两个复数信息信号的下一(N=)两个频率一致的频谱值并且继续传递至方框208。然后采用应用于频谱值v1(f2,t1)和v2(f2,t1)的根据本发明的方法来获得结果频谱值m(f2,t1)。
由此为(N=)两个复数信息信号的所有频率一致的频谱值执行该方法,直至在方框218中获得复数的输出信息信号m(f,t1)。然后在方框220中通过从频域到时域中的逆变换将复数的输出信息信号转换成混合的时间信号sc(t)。
此处还适用的是,在流程图的另一实施例中可同时彼此并行地执行方框206至214来直接获得复数的输出信息信号m(f,t1)。
在图3中详细阐述如何实施图2中的方框208至214的方法。图3a示出了在复平面中作为矢量OP1或OP5绘出的两个频率一致的频谱值v1(f1,t1)和v2(f1,t1),其中,O是复平面的原点。在方框208中将频谱值v1(f1,t1)(=OP1)转换成第一分量OP3和第二分量OP4。如此选择第一分量OP3和第二分量OP4,使得它们在分量OP3和OP4的复数值相加时得到频谱值OP1。在方框208中进一步将频谱值v2(f1,t1)(=OP5)转换成第一分量OP7和第二分量OP8。第一分量OP7和第二分量OP8如此选择,使得它们在分量OP7和OP8的复数值相加时得到频谱值OP5。
第一分量OP3和OP7以及第二分量OP4和OP8的端点位于一个圆K上。这在本发明的该实施例中意味着,第一和第二分量的幅度或矢量长度彼此相等。圆K的半径与两个频谱值v1(f1,t1)和v2(f1,t1)的绝对值(数值)相关。尤其适用以下关系:
第一能量值E1(f1,t1)等于:ABS(v1(f1,t1))2
第二能量值E2(f1,t1)等于:ABS(v2(f1,t1))2
圆K的半径R此时等于:SQRT{(E1+E2)/2}。
因此,对该能量值的数学平均值进行开根号是半径的大小。
在第一实施例中确定半径意味着使用以下假设,即,声场由两个直接声场的叠加组成,其中,两个假定的直接声场相等,导致在内插或外推位置(英语:location)处声场变量的估计值会很宽。
图3b示出了,在方框210(步骤B)中第一分量OP3和OP7如何结合成第一组合分量OP9。在此,第一组合分量的端点P9如下确定。
对此首先在该部分中描述对混合、例如内插或外推的系数的确定。已知的是,从给定位置内插或外推的位置可以用数学方法表示,例如通过线性组合来表示,这将在下文中使用。
如果混合是内插或外推,则线性组合的系数之和等于1。一维线性内插或外推位置L的两个位置L1和L2中给定位置之一的数学表示为
L=L1*c1+L2*c2
其中,c1和c2是系数,且
c1+c2=1。
如果对于L1和L2使用输出相应麦克风信号s1(t)和s2(t)的麦克风的麦克风位置,则c1和c2是根据本发明内插或外推的系数。
第一分量OP3和OP7的内插在图3b中产生第一组合分量OP9。对此,点P9将圆区段P3–K-P7分成两个部分,即,
(圆区段P3–P9的弧长)/(圆区段P9–P7的弧长)=c1/c2。
图3c示出了,在方框212中(步骤C)第二分量OP4和OP8如何结合成第二组合分量OP10。在此,第二组合分量的端点P10如下确定。
如上结合图3b所述,在图3c中也将圆区段P4–K–P8分成两个部分,确切地说通过点P10分成两个部分。对此点P10将圆区段P4–K–P8分成两个部分,即,
(圆区段P4–P10的弧长)/(圆区段P10–P8的弧长)=c1/c2。
图3d示出了在方框214中(步骤D)第一组合分量OP9和第二组合分量OP10如何结合成结果频谱值OP11。这通过使组合分量OP9和OP10的复数值相加实现。
因此,依次重复地或彼此并行地实施上面根据图3描述步骤,并且与图2相结合以便获得在该方法的第一实施例中的复数的输出信息信号。
对此还额外地提出,必须始终针对新的复数频谱值对,例如针对v1(f2,t1)或v2(f2,t1)重新进行半径计算。
在上述方法中进行混合,这引起两个信息时间信号内插。这是因为c1和c2均为正且小于1。上述方法也可引起外推。在这种情况下,两个系数之一c1或c2为负,另一个系数大于1,其中仍然适用:c1+c2=1。这意味着点P9和P10还始终位于圆上,但是位于区段P3–K-P7或P4–K–P8之外。
在根据图4进一步扩展的第二实施例中,两个复数信息信号的混合如下进行。图4a再次示出了在复平面中作为矢量OP1或OP5绘出的两个频率一致的频谱值v1(f1,t1)和v2(f1,t1),其中,O是复平面的原点。在方框208中将频谱值v1(f1,t1)(=OP1)转换成第一分量OP3和第二分量OP4。如此选择第一分量OP3和第二分量OP4,使得它们在分量OP3和OP4的复数值相加时得到频谱值OP1。在方框208中进一步将频谱值v2(f1,t1)(=OP5)转换成第一分量OP7和第二分量OP8。第一分量OP7和第二分量OP8如此选择,使得它们在分量OP7和OP8的复数值相加时得到频谱值OP5。
第一分量OP3和OP7的端点位于一个圆K’上。这在本发明的该实施例中意味着,第一分量OP3和OP7的幅度或矢量长度彼此相等。第二分量OP4和OP8的端点位于一个圆K”上。这在本发明的该实施例中意味着,第二分量OP4和OP8的幅度或矢量长度彼此相等。
两个圆K’和K”的半径此时彼此不等,但是又与两个频谱值v1(f1,t1)和v2(f1,t1)的绝对值(数值)相关。
在第二实施例中假设,两个假设的直接声场中的其中一个占主导地位,从而使得对于在声场中占主导地位的直接声场部分而言,尽可能准确地估计在内插或外推位置处对声场变量的值。特别地,以下公式适用于半径的计算:
EA=(E1+E2)/2+Ed
EB=(E1+E2)/2-Ed
Ed应大于零。另一方面,Ed不允许过大,因为此时就不可能再将两个频谱值之一划分为分量。具有较小矢量长度的一个在这种情况下为图4a中的OP5,并且通过示例的方式示出了Ed的最大值的边界情况,其中还是能够进行划分的,并且可以识别出,频谱值OP5与其分量OP7和OP8是共线的。
圆K′的半径R′此时等于SQRT(EA)。
圆K″的半径R″此时等于SQRT(EB)。
图4b示出了,在方框210(步骤B)中第一分量OP3和OP7如何结合成第一组合分量OP9。在此,第一组合分量的端点P9又以与上面根据图3b所述的相同的方式被确定。
点P9将圆区段P3–K’–P7分成两个部分,即,
(圆区段P3–P9的弧长)/(圆区段P9–P7的弧长)=c1/c2。
图4c示出了,在方框212中(步骤C)第二分量OP4和OP8如何结合成第二组合分量OP10。在此,第二组合分量的端点P10如下确定。
如上结合图4b所述,在图4c中也将圆区段P4–K″–P8分成两个部分,确切地说通过点P10分成两个部分。对此点P10将圆区段P4–K″–P8分成两个部分,即,(圆区段P4–P10的弧长)/(圆区段P10–P8的弧长)=c1/c2。
图4d示出了在方框214中(步骤D)第一组合分量OP9和第二组合分量OP10如何结合成结果频谱值OP11。这通过使组合分量OP9和OP10的复数值相加实现。
因此,依次重复地或彼此并行地实施上面根据图4描述步骤,并且与图2相结合以便获得在该方法的第二实施例中的复数的输出信息信号。
对此还额外地提出,必须始终针对新的复数频谱值对,例如针对v1(f2,t1)或v2(f2,t1)重新进行半径计算。
在上述方法中进行混合,这引起两个信息时间信号内插。这是因为c1和c2均为正且小于1。上述方法也可引起外推。在这种情况下,两个系数之一c1或c2为负,另一个系数大于1,其中仍然适用:c1+c2=1。
图5示出了用于实施如上所述的方法的混合装置的实施例。输入部502和504用于接收(N=)两个复数信息信号v1(f,t1)或v2(f,t1)。输入部502与单元508的输入部506耦合。输入部504与单元520的输入部518耦合。单元508和520一起形成用于将(N=)两个复数信息信号的频率一致的频谱值中的每一个转换成如根据图3a或图4a所述的第一和第二分量的第一单元。这意味着,在经由控制单元530的控制线路542和544的控制的影响下通过单元508或520在其输入部506或518处接收频率一致的频谱值v1(f1,t1)和v2(f1,t1)(在图3a和图4a中的OP1或OP5)且从其中通过该单元产生两个第一分量(在图3a和图4a中的OP3或OP7)和两个第二分量(在图3a和图4a中的OP4或OP8)。第一分量OP3通过单元508输送给其输出部510。第二分量OP4通过单元508输送给其输出部512。第一分量OP7通过单元520输送给其输出部522。第二分量OP8通过单元520输送给其输出部524。
为了计算在图3中的圆K的半径或在图4中的圆K’和K”的半径设置单元540。混合装置的输入部502和504与单元540的对应的输入部532或534耦合。在第二实施例的情况下单元540在经由控制单元530的控制线路546的控制的情况下如上所述从输送给输入部502和504的复数信息信号v1(f,t1)或v2(f,t1)中推导出能量EA和EB。然后单元540从能量EA和EB中推导出圆K’和K”的半径(参见图4a)且将其提供给输出部538或536。单元540的输出部538与单元508或520的输入部514和526耦合以便将圆K’的半径的值输送给单元508和520。单元540的输出部540与单元508或520的输入部516和528耦合以便将圆K”的半径的值输送给单元508和520。
在第一示例性实施例中,在单元540中,仅推导出圆K的半径的值,参见图3a,并将其输送到单元508和520。在第一实施例中仅在单元540和单元508和520之间设有一个连接线路。混合装置还包含单元548。在单元548中在经由控制单元530的控制线路558控制的情况下将由单元508和520产生的两个第一分量OP3和OP7结合以便产生第一组合分量OP9,如根据图3a和图4a所述。对此,单元508的输出部510和单元520的输出部522与单元548的对应的输入部552或554耦合。单元548也需要圆K或K’的半径值,如根据图3b和图4b所述。对此还可设置在单元540和单元548之间的耦合以便输送圆K或K’的半径值。或单元548可从输送给其的两个第一分量OP3或OP7中推到出K或K’的半径值。
为了推到出第一组合分量还需要系数c1和c2。然而,应当注意,稍后将参考图7说明要求系数小于信息信号的数量N。
通过输入部560或562提供这两个系数,并且仅通过一个输入部560或562将这一个系数提供给混合装置。这些输入部与单元548的相关输入部564或566耦合。然后在单元548的输出部556处可获得第一组合分量OP9。
混合装置还包括单元550。在单元550中在经由控制单元530的控制线路568控制的情况下将由单元508或520产生的两个第二OP4和OP8结合以便产生第二组合分量OP10,如根据图3c和图4c所述。对此,单元508的输出部512和单元520的输出部524与单元550的对应的输入部570或572耦合。单元550也需要圆K或K”的半径值,如根据图3c和图4c所述。对此还可设置在单元540和单元550之间的耦合以便输送圆K或K”的半径值。或单元550可从输送给其的两个第二分量OP4或OP8中推到出K或K”的半径值。
为了推到出第二组合分量还需要系数c1和c2。混合装置的输入部560或562对此与单元550的对应的输入部574或576耦合。此时在单元550的输出部578处提供第二组合分量OP10。
混合装置还包括单元580。在单元580中在经由控制单元530的控制线路582控制的情况下将第一和第二组合分量OP9和OP10结合以便产生结果频谱值OP11,如根据图3d和图4d所述。对此,单元548的输出部556和单元578的输出部524与单元580的对应的输入部584或586耦合。单元580的输出部588与混合装置的输出部590耦合。
控制单元530控制混合装置中的单元,使得根据参考图2所述基于生成结果频谱值的步骤重复执行两个复数信息信号的两个频率一致的频谱值,以在输出部590处获得复数的输出信息信号。或混合装置如在图5中所述,为了同时推到出结果频谱值m(f,t1)而多次实现。因此然后应实施控制单元530,以便能够并行地处理。
图6示出了用于混合三个信息信号的混合装置,三个信息信号分别首先从时域转换到频域。
在这种情况下N=3且二维线性内插或外推位置L的三个位置L1、L2和L3中给定位置之一的数学表示为
L=L1*c1+L2*c2+L3*c3
其中,c1、c2和c3是系数,且
c1+c2+c3=1。
如果对于L1、L2和L3使用输出相应麦克风信号s1(t)、s2(t)和s3(t)的麦克风的麦克风位置,则c1、c2和c3是根据本发明内插或外推的系数。
输入部602、603和604用于接收(N=)三个复数信息信号v1(f,t1)、v2(f,t1)或v3(f,t1)。输入部602与单元608的输入部606耦合。输入部603与单元617的输入部607耦合。输入部604与单元620的输入部618耦合。单元608、617和620一起形成用于将(N=)三个复数信息信号的频率一致的频谱值中的每一个转换成如根据图3a或图4a所述的第一和第二分量的第一单元。这意味着,在经由控制单元630的控制线路642、643和644的控制的影响下通过单元608或617或620在其输入部606或607或618处接收频率一致的频谱值v1(f1,t1)(在图3a和图4a中的OP1)、v2(f1,t1)(在图3a和图4a中的OP5)和v3(f1,t1)且从其中通过该单元产生三个第一分量(OP3、OP7、OP12)和三个第二分量(OP4、OP8、OP13)。第一分量OP3通过单元608输送给其输出部610。第二分量OP4通过单元608输送给其输出部612。第一分量OP7通过单元617输送给其输出部611。第二分量OP8通过单元617输送给其输出部613。第一分量OP12通过单元620输送给其输出部622。第二分量OP13通过单元620输送给其输出部624。
为了计算在图3中的圆K的半径或在图4中的圆K’和K”的半径设置单元640。混合装置的输入部602、603和604与单元640的对应的输入部632或633或634耦合。在第二实施例的情况下单元640在经由控制单元630的控制线路646的控制的情况下如下所述从输送给输入部602、603和604的复数信息信号v1(f,t1)、v2(f,t1)和v3(f,t1)中推导出能量EA和EB。
第一能量值E1(f1,t1)等于:ABS(v1(f1,t1))2
第二能量值E2(f1,t1)等于:ABS(v2(f1,t1))2
第三能量值E3(f1,t1)等于:ABS(v3(f1,t1))2
圆K的半径R此时等于:SQRT{(E1+E2+E3)/3}。
对于K’和K”的推导适用以下内容:
在这种情况下,单元640从能量值EA和EB中如下地推导出圆K’和K”的半径(参见图4a)并且将其提供给输出部638或636。
EA=(E1+E2+E3)/3+Ed
EB=(E1+E2+E3)/3-Ed
Ed应大于零。另一方面,Ed不允许过大,因为此时就不可能再将三个频谱值之一划分为分量。
圆K′的半径R′此时等于SQRT(EA)。
圆K″的半径R″此时等于SQRT(EB)。
单元640的输出部638与单元608或617或620的输入部614、615和626耦合以便将圆K’的半径的值输送给单元608、617和620。单元640的输出部636与单元608或617或620的输入部616、619和628耦合以便将圆K”的半径的值输送给单元608、617和620。
在第一示例性实施例中,在单元640中,仅推导出圆K的半径的值,参见图3a,并将其输送到单元608、617和620。在第一实施例中仅在单元640和单元608、617和620之间设有一个连接线路。
混合装置还包含单元648。在单元648中在经由控制单元630的控制线路658控制的情况下将由单元608或617或520产生的三个第一分量OP3、OP7和OP12结合以便产生第一组合分量OP19。这将根据图7进一步阐述。图7在复平面中示出了三个分量OP3、OP7和OP12以及第一组合分量OP19。分量OP3与复平面的轴线、例如水平轴线的角度等于α1。分量OP7与水平轴线的角度等于α2。分量OP12与水平轴线的角度等于α3。以及组合分量OP19与水平轴线的角度等于α4。此时在角α1、α2、α3和α4之间有以下关系:
α4=c1*α1+c2*α2+c3*α3公式(1)或
α4′=c2*α2′+c3*α3′公式(2)
其中,α4′是OP3和OP19之间的角,α2′是OP3和OP7之间的角,α3′是OP3和OP12之间的角。
在利用公式(2)来推导OP19时,c1=0,因此c2+c3=1。
对此,单元608的输出部610、单元617的输出部611和单元620的输出部622与单元648的对应的输入部652或654或655耦联。单元648也需要圆K或K’的半径值,如根据图3a和图4a所述。对此还可设置在单元640和单元648之间的耦合以便输送圆K或K’的半径值。或单元648可从输送给其的三个第一分量OP3、OP7或OP12中推到出K或K’的半径值。
如图7中所述,在单元648中从OP3、OP7或OP12中推导出OP19,其中,如公式(1)或(2)所给出地,利用三个系数c1、c2和c3或两个系数c2和c3。
三个或两个系数经由输入部660、662、663或输入部662、663输送给混合装置。该输入部与单元648的对应的输入部664或666或667耦合。然后在单元648的输出部656处提供第一组合分量OP19。
混合装置还包括单元650。在单元650中在经由控制单元630的控制线路668控制的情况下将由单元608或617或620产生的三个第二OP4、OP8和OP13结合以便以与根据图7所述相同的方式产生第二组合分量OP20。对此,单元608、613的输出部612和单元620的输出部624与单元650的对应的输入部670或672或673耦合。单元650也需要圆K或K”的半径值,如根据图3c和图4c所述。对此还可设置在单元640和单元650之间的耦合以便输送圆K或K”的半径值。或单元650可从输送给其的三个第二分量OP4、OP8或OP13中推到出K或K”的半径值。
为了推到出第二组合分量OP20还需要系数c1、c2和c3。混合装置的输入部660、662或663对此与单元650的对应的输入部674或676或677耦合。此时在单元650的输出部678处提供第二组合分量OP20。
混合装置还包括单元680。在单元680中在经由控制单元630的控制线路682控制的情况下将第一和第二组合分量OP19和OP20结合以便产生结果频谱值OP21,如根据图3d和图4d所述。对此,单元648或650的输出部656和678与单元680的对应的输入部684或686耦合。单元680的输出部688与混合装置的输出部690耦合。
控制单元630控制混合装置中的单元,使得根据参考图2所述基于生成结果频谱值的步骤重复执行三个复数信息信号的三个频率一致的频谱值,以在输出部690处获得复数的输出信息信号。或混合装置如在图6中所述,为了同时推到出结果频谱值m(f,t1)而多次实现。
因此针对N大于3,该装置可扩展成混合N大于3个的N个复数信息信号。
因此装置针对N=4获得:
-除了图6中的输入部602、603和604以外的第四输入部,以便用于接收第四复数信息信号v4(f,t1),
-用于将第四复数信息信号v4(f,t1)输送给单元640的额外的输入部的额外的线路,
-除了图6中的输入部602、603和604以外的额外的单元,
-用于通过图6中的控制单元630控制额外的单元的额外的控制线路,
-单元640中的额外的线路,以便将半径值输送给额外的单元,
-从额外的单元至图6中的单元648或650的每个额外的输入部的两个额外的输出线路,以及
-除了图6中的输入部660、662、663之外的第四输入部以用于接收第四系数c4。
类似于上述针对N=2和N=3,三维线性内插或外推位置L的三个位置L1、L2、L3和L4中给定位置之一的数学表示为
L=L1*c1+L2*c2+L3*c3+L4*c4
其中,c1、c2、c3和c4是系数,且
c1+c2+c3+c4=1。
如果对于L1、L2、L3和L4使用输出相应麦克风信号s1(t)、s2(t)、s3(t)和s4(t)的麦克风的麦克风位置,则c1、c2、c3和c4是根据本发明内插或外推的系数。下面进行总结。
将频率一致的频率值分成第一和第二分量并且结合第一或第二分量是指,假设声场由两个直接的声场的重叠构成,其中每一个分量与其中一个假设的直接声场对应。通过假设可为该分量采用混合(内插或外推),其模拟了直接声场的声场变量与空间位置的物理关系。使用该假设使得,只要声场是由来自最多两个声源的声波产生的,则混合(内插或外推的)信号可以很好地估计内插或外推的位置处的声场测量值。
由于所有第一分量的振幅相等且所有第二分量的振幅相等,因此物理关系的模拟可能非常简单,即仅限于具有平面波前的直接声场。
内插或外推分量的平均能量与所有麦克风信号的平均能量相等,意味着使用了一个附加假设,即,空间中声场变量的平均能量是恒定的。作为此附加假设的结果,如果最多两个直接声音分量的假设与实际情况不符,则内插或外推信号仍然是对内插或外推位置处声场测量变量的值的有用估计。
所有第一分量的能量相等意味着,不必对第一分量的能量进行内插或外推,但是可以简单地将第一内插分量或外推分量的能量等同于它们。后者就是这样的。结果,第一内插或外推减小为第一分量的相位的内插或外推。
对于第二分量,第二内插或外推分量,第二内插或外推以及第二分量的相位,以相同的方式同样适用。
Claims (15)
1.一种用于混合N个信息时间信号的方法,所述信息时间信号首先从时域转换成频域的N个复数信息信号之一,其中,N是大于1的整数,其中,所述方法包括以下步骤:
(a)在N个复数信息信号的频率一致的频谱值中分别转换成第一和第二分量,
(b)将在频率一致的N个频谱值中的N个第一分量结合成第一组合分量,
(c)将在频率一致的N个频谱值中的N个第二分量结合成第二组合分量,
(d)所述第一组合分量和第二组合分量结合成结果频谱值,
(e)所述步骤(a)至(d)也针对N个复数信息信号的频率一致的其他频谱值进行以便产生其他的结果频谱值,
(f)其中,由此获得的结果频谱值形成复数的输出信息信号。
2.根据权利要求1所述的方法,其特征在于,为了在步骤(b)中推导出第一组合分量,在步骤(a)中推导出的第一分量具有彼此基本相等的幅度。
3.根据权利要求1所述的方法,其特征在于,为了在步骤(c)中推导出第二组合分量,在步骤(a)中推导出的第二分量具有彼此基本相等的幅度。
4.根据权利要求1所述的方法,其特征在于,为了在步骤(b)或步骤(c)中推导出第一和第二组合分量,在步骤(a)中推导出的第一和第二分量具有彼此基本相等的幅度。
5.根据前述权利要求中任一项所述的方法,其特征在于,在步骤(a)中从复数信息信号的频谱值转换成第一和第二分量如此实现,使得所述第一分量和所述第二分量的复数值相加得到频谱值。
6.根据权利要求1所述的方法,其特征在于,在步骤(d)中所述第一和第二组合分量结合以便获得结果频谱值如此实现,使得所述第一组合分量和所述第二组合分量的复数值相加得到结果频谱值。
7.根据权利要求1所述的方法,其特征在于,将N个第一分量作为从复平面的原点出发的矢量在所述复平面中示出,并且所述矢量的端点位于所述复平面原点周围的一个圆上,使得N个信息信号的混合为c1与c2与c3与...cN的比,其中c1+c2+c3+...+cN=1,且在步骤(b)中所述N个第一分量结合以便获得所述第一组合分量如此实现,使得所述第一组合分量作为从原点出发的矢量在所述复平面中示出并且所述第一组合分量的端点位于一个圆上,其中,所述第一组合分量与所述复平面的轴线的夹角相对于N个第一分量与所述轴线的夹角的比例如下:
αC=c1*α1+c2*α2+c3*α3+...+cN*αN,
其中,αC所述第一组合分量与所述轴线的夹角且α1至αN是N个第一分量与轴线的夹角。
8.根据权利要求1所述的方法,其特征在于,将N个第二分量作为从复平面的原点出发的矢量在所述复平面中示出,并且所述矢量的端点位于复平面原点周围的一个圆上,使得N个信息信号的混合为c1与c2与c3与...cN的比,其中c1+c2+c3+...+cN=1,且在步骤(c)中所述N个第二分量结合以便获得所述第二组合分量如此实现,使得所述第二组合分量作为从原点出发的矢量在所述复平面中示出并且所述第二组合分量的端点位于一个圆上,其中,所述第二组合分量与所述复平面的轴线的夹角相对于N个第二分量与所述轴线的夹角的比例如下:
αC=c1*α1+c2*α2+c3*α3+...+cN*αN,
其中,αC是所述第二组合分量与所述轴线的夹角且α1至αN是N个第二分量与轴线的夹角。
9.根据权利要求1或7所述的方法,其特征在于,N=2,所述两个第一分量作为从复平面的原点出发的矢量在所述复平面中示出,并且所述矢量的端点位于所述复平面原点周围的一个圆上,使得二个信息信号的混合为c1/c2的比,其中c1+c2=1,且在步骤(b)中所述二个第一分量结合以便获得所述第一组合分量如此实现,使得所述第一组合分量作为从原点出发的矢量在所述复平面中示出并且所述第一组合分量的端点位于一个圆上,其中,从其中一个第一分量的端点至所述第一组合分量的端点的圆区段的弧长和从另一第一分量的端点至所述第一组合分量的端点的圆区段的弧长的比例为c1/c2。
10.根据权利要求1或8所述的方法,其特征在于,N=2,所述两个第二分量作为从复平面的原点出发的矢量在所述复平面中示出,并且所述矢量的端点位于所述复平面原点周围的一个圆上,使得二个信息信号的混合为c1/c2的比,其中c1+c2=1,且在步骤(c)中所述二个第二分量结合以便获得所述第二组合分量如此实现,使得所述第二组合分量作为从原点出发的矢量在所述复平面中示出并且所述第二组合分量的端点位于一个圆上,其中,从其中一个第二分量的端点至所述第二组合分量的端点的圆区段的弧长和从另一第二分量的端点至所述第二组合分量的端点的圆区段的弧长的比例为c1/c2。
11.根据权利要求4所述的方法,其特征在于,所述第一和第二分量作为从复平面的原点出发的矢量在所述复平面中示出,并且所述矢量的端点位于所述复平面原点周围的一个圆上,,其中,所述圆的半径如下推导:
其中,|vi(f1,t1)|是N个复数信息信号的频率一致的N个频谱值的绝对值(数值)。
13.用于实施根据前述权利要求中任一项所述的方法的混合装置,设有用于接收N个复数信息信号的输入部以及用于将N个复数信息信号混合成复数的输出信息信号的混合单元,其特征在于,所述混合单元设有:
a.用于将N个复数信息信号的频率一致的频谱值中的每一个转换成第一和第二分量的第一单元(508、520;608、617、628),
b.用于将频率一致的N个频谱值的N个第一分量结合成第一组合分量的第二单元(548;648),
c.用于将频率一致的N个频谱值的N个第二分量结合成第二组合分量的第三单元(550;650),
d.用于将所述第一和第二组合分量结合成结果频谱值的第四单元(580;680),
e.用于控制第一至第四单元以针对N个复数信息信号的频率一致的频谱值中的其他频谱值重复推导出结果频谱值或用于并行地控制多个第一、第二、第三和第四单元以从频率一致的频谱值中推导出结果频谱值的控制单元(530;630),
f.用于输送作为复数的输出信息信号的由此推导出的结果频谱值的输出部(590;690)。
14.根据权利要求13所述的用于实施根据权利要求11所述的方法的混合装置,其特征在于,所述混合单元还设有用于从N个复数信息信号的频率一致的频谱值中推导出所述半径值的单元(540;640)。
15.根据权利要求13所述的用于实施根据权利要求12所述的方法的混合装置,其特征在于,所述混合单元还设有用于从N个复数信息信号的频率一致的频谱值中推导出两个半径值的单元(540;640)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IT102017000040732 | 2017-04-12 | ||
IT102017000040732A IT201700040732A1 (it) | 2017-04-12 | 2017-04-12 | Verfahren und vorrichtung zum mischen von n informationssignalen |
PCT/DE2018/100326 WO2018188697A1 (de) | 2017-04-12 | 2018-04-10 | Verfahren und vorrichtung zum mischen von n informationssignalen |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110720226A true CN110720226A (zh) | 2020-01-21 |
CN110720226B CN110720226B (zh) | 2021-12-31 |
Family
ID=59579879
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880037021.7A Active CN110720226B (zh) | 2017-04-12 | 2018-04-10 | 用于混合n个信息信号的方法和装置 |
Country Status (7)
Country | Link |
---|---|
US (1) | US10834502B2 (zh) |
EP (1) | EP3610655B1 (zh) |
JP (1) | JP7164075B2 (zh) |
CN (1) | CN110720226B (zh) |
DE (1) | DE112018001988A5 (zh) |
IT (1) | IT201700040732A1 (zh) |
WO (1) | WO2018188697A1 (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003015464A2 (en) * | 2001-08-08 | 2003-02-20 | Dspfactory Ltd | Directional audio signal processing using an oversampled filterbank |
GB2453118A (en) * | 2007-09-25 | 2009-04-01 | Motorola Inc | Generating a speech audio signal from multiple microphones with suppressed wind noise |
DE102008004674A1 (de) * | 2007-12-17 | 2009-06-18 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Signalaufnahme mit variabler Richtcharakteristik |
DE102009052992B3 (de) * | 2009-11-12 | 2011-03-17 | Institut für Rundfunktechnik GmbH | Verfahren zum Abmischen von Mikrofonsignalen einer Tonaufnahme mit mehreren Mikrofonen |
CN102016985A (zh) * | 2008-03-04 | 2011-04-13 | 弗劳恩霍夫应用研究促进协会 | 对输入数据流进行混合以及从中产生输出数据流 |
CN103355001A (zh) * | 2010-12-10 | 2013-10-16 | 弗兰霍菲尔运输应用研究公司 | 用以利用下变频混频器来分解输入信号的装置和方法 |
US20140105416A1 (en) * | 2012-10-15 | 2014-04-17 | Nokia Corporation | Methods, apparatuses and computer program products for facilitating directional audio capture with multiple microphones |
WO2014085978A1 (en) * | 2012-12-04 | 2014-06-12 | Northwestern Polytechnical University | Low noise differential microphone arrays |
US20140270241A1 (en) * | 2013-03-15 | 2014-09-18 | CSR Technology, Inc | Method, apparatus, and manufacture for two-microphone array speech enhancement for an automotive environment |
US20140286508A1 (en) * | 2011-10-05 | 2014-09-25 | Institut Fur Rundfunktechnik Gmbh | Interpolation Circuit for Interpolating a First and a Second Microphone Signal |
WO2016160821A1 (en) * | 2015-03-30 | 2016-10-06 | Bose Corporation | Adaptive mixing of sub-band signals |
US20160379664A1 (en) * | 2015-06-25 | 2016-12-29 | Htc Corporation | Sound processing device and method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2600637A1 (en) * | 2011-12-02 | 2013-06-05 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Apparatus and method for microphone positioning based on a spatial power density |
ITTO20130028A1 (it) | 2013-01-11 | 2014-07-12 | Inst Rundfunktechnik Gmbh | Mikrofonanordnung mit verbesserter richtcharakteristik |
MX358244B (es) | 2013-05-09 | 2018-08-10 | Procter & Gamble | Composicion acondicionadora para el cuidado del cabello que comprende histidina. |
DE102013105375A1 (de) * | 2013-05-24 | 2014-11-27 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Tonsignalerzeuger, Verfahren und Computerprogramm zum Bereitstellen eines Tonsignals |
DE102014108345A1 (de) * | 2014-06-13 | 2015-12-17 | Intel Ip Corp. | Verfahren, Vorrichtung und Computerprogramm zum Bestimmen eines interpolierten, komplexwertigen Abtastwerts |
-
2017
- 2017-04-12 IT IT102017000040732A patent/IT201700040732A1/it unknown
-
2018
- 2018-04-10 CN CN201880037021.7A patent/CN110720226B/zh active Active
- 2018-04-10 DE DE112018001988.8T patent/DE112018001988A5/de not_active Withdrawn
- 2018-04-10 EP EP18730965.3A patent/EP3610655B1/de active Active
- 2018-04-10 JP JP2019554960A patent/JP7164075B2/ja active Active
- 2018-04-10 WO PCT/DE2018/100326 patent/WO2018188697A1/de unknown
- 2018-04-10 US US16/604,641 patent/US10834502B2/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003015464A2 (en) * | 2001-08-08 | 2003-02-20 | Dspfactory Ltd | Directional audio signal processing using an oversampled filterbank |
US20080112574A1 (en) * | 2001-08-08 | 2008-05-15 | Ami Semiconductor, Inc. | Directional audio signal processing using an oversampled filterbank |
GB2453118A (en) * | 2007-09-25 | 2009-04-01 | Motorola Inc | Generating a speech audio signal from multiple microphones with suppressed wind noise |
DE102008004674A1 (de) * | 2007-12-17 | 2009-06-18 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Signalaufnahme mit variabler Richtcharakteristik |
CN102016985A (zh) * | 2008-03-04 | 2011-04-13 | 弗劳恩霍夫应用研究促进协会 | 对输入数据流进行混合以及从中产生输出数据流 |
DE102009052992B3 (de) * | 2009-11-12 | 2011-03-17 | Institut für Rundfunktechnik GmbH | Verfahren zum Abmischen von Mikrofonsignalen einer Tonaufnahme mit mehreren Mikrofonen |
CN103355001A (zh) * | 2010-12-10 | 2013-10-16 | 弗兰霍菲尔运输应用研究公司 | 用以利用下变频混频器来分解输入信号的装置和方法 |
US20140286508A1 (en) * | 2011-10-05 | 2014-09-25 | Institut Fur Rundfunktechnik Gmbh | Interpolation Circuit for Interpolating a First and a Second Microphone Signal |
US20140105416A1 (en) * | 2012-10-15 | 2014-04-17 | Nokia Corporation | Methods, apparatuses and computer program products for facilitating directional audio capture with multiple microphones |
WO2014085978A1 (en) * | 2012-12-04 | 2014-06-12 | Northwestern Polytechnical University | Low noise differential microphone arrays |
US20140270241A1 (en) * | 2013-03-15 | 2014-09-18 | CSR Technology, Inc | Method, apparatus, and manufacture for two-microphone array speech enhancement for an automotive environment |
WO2016160821A1 (en) * | 2015-03-30 | 2016-10-06 | Bose Corporation | Adaptive mixing of sub-band signals |
US20160379664A1 (en) * | 2015-06-25 | 2016-12-29 | Htc Corporation | Sound processing device and method |
Also Published As
Publication number | Publication date |
---|---|
IT201700040732A1 (it) | 2018-10-12 |
DE112018001988A5 (de) | 2019-12-19 |
JP2020517157A (ja) | 2020-06-11 |
CN110720226B (zh) | 2021-12-31 |
EP3610655A1 (de) | 2020-02-19 |
EP3610655B1 (de) | 2021-01-06 |
JP7164075B2 (ja) | 2022-11-01 |
US10834502B2 (en) | 2020-11-10 |
WO2018188697A1 (de) | 2018-10-18 |
US20200162814A1 (en) | 2020-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0091926B1 (en) | Improved method and apparatus for cancelling vibrations | |
JP2008519491A5 (zh) | ||
Alazard et al. | Control of water waves | |
US10164609B2 (en) | Fractional scaling digital signal processing | |
KR101021895B1 (ko) | 음향 필드 표현을 처리하기 위한 방법 및 시스템 | |
US20090134714A1 (en) | Compensation of harmonics of electrical network | |
Dong et al. | On nonlocal problems for fractional differential equations in Banach spaces | |
CN105527516B (zh) | 一种通道校准方法及电子设备 | |
CN110720226B (zh) | 用于混合n个信息信号的方法和装置 | |
JP5010148B2 (ja) | 3次元パンニング装置 | |
US20180021023A1 (en) | Ultrasound Imaging Device and Ultrasound Signal Processing Method | |
JPS5983063A (ja) | スイツチング信号を発生するスイツチング装置および電気信号をスイツチングする方法 | |
KR20100005349A (ko) | 신호원 방향 추적 장치 및 방법 | |
JP2015021824A (ja) | 類識別フィルタ生成装置及び目標類識別装置 | |
Arend et al. | Synthesis of near-field HRTFs by directional equalization of far-field datasets | |
CN104137567A (zh) | 用于内插第一和第二麦克风信号的内插电路 | |
Sönmez et al. | Effect of load increase and power system stabilizer on stability delay margin of a generator excitation control system | |
JP2016180679A5 (zh) | ||
CN112005087A (zh) | 声学分析系统 | |
CN109164437B (zh) | 用于矢量阵列流行的连续可变小数时延估计方法及装置 | |
JP2021521732A (ja) | アレイアンテナを測定するビーム合成方法及び装置 | |
CN112511208B (zh) | 一种频域宽带波束形成多滤波器合并处理方法及系统 | |
JP2012037429A (ja) | ビームフォーミング装置用フィルタの作成方法およびビームフォーミング装置用フィルタの作成装置 | |
Udupa et al. | Power system Vertical Division State Estimation (VDSE) application to super node area | |
JP2005078559A (ja) | 特性不明システムの同定装置 |
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 |