多通道回声消除方法、多通道声音传送方法、立体回声消除器、立体声音传送装置和传递函数计算装置

摘要

立体声信号直接从扬声器(SP(L)，SP(R))再生。利用立体声信号的和信号和差信号作为参考信号，并根据参考信号与传声器收集的声音信号的互频谱计算，进行计算以得到扬声器(SP(L)，SP(R))与传声器(MC(L)，MC(R))之间的四个音频传送系统的传递函数。所得到的传递函数经逆付里叶变换以得到脉冲响应，这些脉冲响应被设置到滤波器装置(40－1至40－4)中以生成回声消除信号并进行回声消除。这解决了在多通道声音信号的回声消除技术中不确定系数的问题。

说明书

技术领域

本发明涉及多通道音频信号的回声消除技术和传递函数计算技术，以及利用新技术解决不确定系数的问题。

背景技术

在电话会议系统使用的双通道立体声音传输中，通常总是有回声消除器的不确定系数问题，为解决这一问题，已提出各种技术(见Journal ofthe Institute of Electronics，Information and Communication Enginer，vol.81，No.3，pp.266-274，March 1998)。作为解决不确定系数问题的技术之一，有一种方法是降低通道间的相关性。作为这方面的具体技术，传统上已建议的有随机噪声的相加，通过滤波器排除相关性，通道间的频率偏移，使用交错梳状滤波器，非线性处理(专利公报No.H10-190848)等等。

根据上述常规的技术，因为原始立体声信号经过处理然后再生，所以就有一个问题，即在再生信号中，或多或少发生恶化。还有，当处理过程复杂时，在再生信号中就会有延迟，以致在电话会议等等中的会话发生困难。再有，当处理过程复杂且处理电路的处理能力差时，会有这样一些情况，即尽管进行回声消除处理，但难以实时更新回声消除器的系数。

本发明提供一种多通道音频信号的回声消除技术，和传递函数计算技术，解决了常规技术中的上述问题。

发明内容

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的，并且彼此之间具有相关性的多通道声音，被所述传声器收集；所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各个扬声器再生的相应单个信号或通过适当地合成所述单个信号所得的多个合成信号，分别生成回声消除信号；并且，从所述一个或多个传声器的相应单个收集的音频信号，或通过从适当地合成所述单个收集的音频信号的所得的多个合成信号中减去所述回声消除信号，从而实现回声消除；另外，在此方法中，利用一组多个低相关合成信号作为参考信号(表示涉及用于估算传递函数或合成传递函数的那些信号)，它们相应于通过适当地合成所述多通道音频信号所得的信号，并且它们彼此之间具有比所述多通道音频信号之间低的相关性(例如，通过适当地合成所述多通道音频信号，生成彼此之间比所述多通道音频信号之间有较低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，由此设置相应的滤波器特性。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，作为参考信号，分别推导各个音频传送系统的单个传递函数或适当地合成这些单个单个收集的音频信号所得的多个合成传递函数，并设置相应的滤波器特性，因而能进行回声消除。与此相应，因为多通道信号能够不受到或很少受到作用于多通道信号而导致其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用多个低相关合成信号组作为参考信号，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算：基于多个低相关合成信号与传声器的单个收集的音频信号，或通过适当地合成所述单个收集的音频信号所得的多个合成信号之间的互频谱计算，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，进一步，基于所述互频谱计算，分别推导所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，这个计算可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱；并就每一互频谱，在预定时间周期对它们作整体平均，分别推导所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，将由所述各个扬声器再生的，彼此之间具有相关性的多通道声音，被所述传声器收集；所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给由所述各个扬声器再生的相应单个信号或通过适当地合成所述单个信号所得的多个合成信号，分别生成回声消除信号；并且，从所述一个或多个传声器的相应单个收集的音频信号，或通过适当地合成所述单个收集的音频信号的所得的多个合成信号中，减去所述回声消除信号，从而实现回声消除；在此方法中，利用与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，作为参考信号，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差，以便连续地将相应的滤波器特性更新至消除所述估算误差的值，从而实现回声消除。与此相应，由于多通道信号能够不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。还有，能够实时地更新滤波器特性。利用多个低相关合成信号组作为参考信号，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述多个低相关合成信号和从所述一个或多个传声器的相应单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数中，减去回声消除信号所得的回声消除误差信号之间的互频谱计算，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差，进一步，基于所述互频谱计算，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差的计算，可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与从所述一个或多个传声器的相应单个收集的音频信号或单个收集的音频信号而得的多个合成信号中减去回声消除信号所得的回声消除误差信号之间的互频谱；并就每一互频谱在预定时间周期对它们作整体平均，分别推导所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。进一步，检测多个低相关合成信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的多通道声音传送方法是这样一种方法，涉及两个空间，其中每个空间各形成所述多个音频传送系统，执行上述任何一种多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的多通道音频信号。与此相应，经过回声消除处理的多通道音频传输能在两个地点之间进行，它能应用于例如电话会议系统或诸如此类系统。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集，两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数分别被估算，以便设置相应的滤波器特性，通过将所述设置的滤波器特性，提供给由所述各个扬声器再生的相应单个信号或通过适当地合成所述单个信号所得的多个组信号，分别生成回声消除信号；并从所述一个或两个传声器的相应单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成信号中，减去所述回声消除信号，由此进行回声消除；另外，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，由此设置相应的滤波器特性。根据本发明，因为立体声音频信号的和信号和差信号彼此之间有低的相关性，所以利用和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的传递函数或它们的合成传递函数，以便设置相应的滤波器特性，从而能进行回声消除。与此相应，由于多通道信号能够不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算：基于和信号和差信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱计算，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。进一步，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算：通过推导所述立体音音频信号的和信号和差信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。

本发明的立体回声消除方法是这样的方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集，两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性，提供给将由所述各个扬声器再生的相应单个信号或通过适当地合成所述单个信号所得的多个组信号，分别生成回声消除信号；并从所述一个或两个传声器的相应单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成信号中，减去所述回声消除信号，由此进行回声消除；另外，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，因为立体声音频信号的和信号和差信号彼此之间有低的相关性，所以利用和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的传递函数或它们的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，因而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现原始的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差的计算，可以是例如这样的计算：基于所述立体声音频信号的和信号和差信号与从所述一个或两个传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成信号中，减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差。基于所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差的计算，可以是例如这样的计算：通过推导所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间互频谱，并就每一互频谱在预定的时间周期对它们作整体平均，分别推导所述两个或四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的估算误差。进一步，检测所述立体声音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音频传输方法是这样一种方法，涉及两个空间，它们各形成所述四个音频传送系统，执行上述任何多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的立体声音频信号。与此相应，经过回声消除处理的立体声音频传输能在两个地点之间进行，它能应用于例如电话会议系统或诸如此类系统。

在本发明中，由适当地合成将由各个扬声器再生的单个信号所得到的多个合成信号，和用作参考信号的多个低相关合成信号{在本说明书中，“低相关合成信号”包括不相关信号(不相关合成信号)的含义}，可以是例如公用的信号。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的、彼此有相关性的多通道声音(例如，多通道立体声音如两通道、四通道、…)被所述传声器收集；所述多个音频传送系统的传递函数分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给由所述各个扬声器再生的相应信号，分别生成回声消除信号；并且，从相应的所述一个或多个传声器收集的音频信号中，减去所述回声消除信号，从而实现回声消除；另外，在此方法中，利用与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间比所述多通道音频信号之间有较低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导各个音频传送系统的传递函数，由此设置相应的滤波器特性。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号作为参考信号，分别推导各个音频传送系统的传递函数，设置相应的滤波器特性，因此能进行回声消除。与此相应，多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用多个低相关合成信号组作为参考信号，分别推导各个音频传送系统的传递函数的计算，可以是例如，基于多个低相关合成信号与各个传声器收集的音频信号之间的互频谱计算，分别推导各个音频传送系统的传递函数的计算。进一步，基于所述互频谱计算，分别推导所述多个音频传送系统的传递函数的计算可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号的计算，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与各个传声器收集的音频信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导在所述多个音频传送系统中合成多个适当的系统的传递函数所得的多种合成传递函数，从而基于所述多种合成传递函数，推导所述多个音频传送系统的传递函数。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的、彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给所述扬声器再生的相应信号，分别生成回声消除信号；并且，从相应的所述一个或多个传声器收集的音频信号中，减去所述回声消除信号，从而实现回声消除；另外，在此方法中，利用与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间更低相关性的一组多个低相关合成信号，作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间比所述多通道音频信号之间有较低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导各个音频传送系统的传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导各个音频传送系统的传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除所述估算误差的值，因此能进行回声消除。与此相应，由于多通道信号能够不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话，也能实时地更新滤波器特性。例如，滤波器特性能在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)内被更新。利用多个低相关合成信号组作为参考信号，分别推导各个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算，基于所述多个低相关合成信号与从所述一个或多个传声器收集的音频信号中减去相应的回声消除信号所得的回声消除误差信号之间的互频谱计算，分别推导各个音频传送系统的传递函数的估算误差的计算。进一步，基于所述互频谱计算，分别推导所述多个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算：通过相加和相减，合成所述多通道音频信号，生成彼此之间比所述多通道音频信号之间有较低相关性的多个低相关合成信号；推导所述多个低相关合成信号与从所述一个或多个传声器收集的音频信号中，减去相应的回声消除信号所得的回声消除误差信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导在所述多个音频传送系统中合成多个适当系统的传递函数的估算误差而得的多种传递函数合成的估算误差，从而基于所述多种传递函数合成的估算误差，推导所述多个音频传送系统的传递函数的估算误差。另外，检测所述多个低相关合成信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

另外，基于所述互频谱计算，分别推导所述多个音频传送系统的传递函数的计算，可以是例如这样的计算：通过对所述多通道音频信号加以主要分量分析生成多个互相不相关的合成信号；推导所述多个不相关的合成信号与各个传声器收集的音频信号之间的互频谱；并就每一互频谱，在预定的时间周期对它们作整体平均，由此基于整体平均值，推导所述多个音频传送系统的传递函数。

还有，基于所述互频谱计算，分别推导所述多个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算：通过对所述多通道音频信号加以主要分量分析生成多个互相正交不相关的合成信号的计算，推导所述多个不相关合成信号与从所述一个或多个传声器收集的音频信号中减去相应的回声消除信号所得到的回声消除误差信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，由此基于该整体平均值推导所述多个音频传送系统的传递函数的估算误差。在这种情况下，可以这样安排，即对除了由所述扬声器再生的声音以外的一些声音输入至所述传声器所致的通话重叠进行检测，当检测到通话重叠时，使所述滤波器特性的更新周期相对地加长，反之，当未检测到通话重叠时，使所述滤波器特性的更新周期相对地缩短，这样，当通话重叠存在时，能充分地收敛估算误差，而当通话重叠不存在时，能加快估算误差的收敛。

本发明的多通道声音传送方法是这样一种方法，涉及两个空间，它们各形成所述多个音频传送系统，执行上述任何多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的多通道音频信号。与此相应，经过回声消除处理的多通道音频传输能在两个地点之间进行，它能应用于例如电话会议系统或诸如此类系统。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的传递函数分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给所述各个扬声器再生的相应信号，分别生成回声消除信号；并且，从相应的所述一个或两个传声器收集的音频信号中减去回声消除信号，从而实现回声消除；另外，在此方法中，利用所述立体声音频信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数，由此设置相应的滤波器特性。根据本发明，因为立体声音频信号的和信号和差信号彼此之间有低的相关性，所以利用和信号和差信号号作为参考信号，分别推导两个或四个音频传送系统的传递函数，以便设置相应的滤波器特性，从而能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的计算，可以是例如这样的计算，基于和信号和差信号与各个传声器收集的音频信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的计算。进一步，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与各个传声器收集的音频信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导在所述两个或四个音频传送系统中由合成多个适当的系统的传递函数而得到的多种合成传递函数，从而基于所述多种合成传递函数，推导所述两个或四个音频传送系统的传递函数。另外，可以计算四个系统情况下的互频谱计算，例如，所述和信号与第一传声器收集的音频信号之间，所述和信号与第二传声器收集的音频信号之间，所述差信号与第一传声器收集的音频信号之间，以及所述差信号与第二传声器收集的音频信号之间的各个互频谱。还有，所述合成传递函数可包括例如：第一合成传递函数，其是第一扬声器与第一传声器之间的传递函数和第二扬声器与第一传声器之间的传递函数的和；第二合成传递函数，其是第一扬声器与第一传声器之间传递函数和第二扬声器与第一传声器之间的传递函数的差；第三合成传递函数，其是第一扬声器与第二传声器之间的传递函数和第二扬声器与第二传声器之间的传递函数的和；以及第四合成传递函数，其是第一扬声器与第二传声器之间的传递函数和第二扬声器与第二传声器之间的传递函数的差。推导所述四个音频传送系统的传递函数的计算，可包括例如：根据第一合成传递函数与第二合成传递函数的和，推导第一音频传送系统的传递函数的计算；根据第一合成传递函数与第二合成传递函数之间的差，推导第二音频传送系统的传递函数的计算；根据第三合成传递函数与第四合成传递函数的和，推导第三音频传送系统的传递函数的计算；以及根据第三合成传递函数与第四合成传递函数之间的差，推导第四音频传送系统的传递函数的计算。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给所述各个扬声器再生的相应信号，分别生成回声消除信号；并且，从相应的所述一个或两个传声器收集的音频信号中减去所述回声消除信号，从而实现回声消除；另外，在此方法中，利用所述立体声音频信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，因为立体声音频信号的和信号和差信号彼此之间有低的相关性，所以利用和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，因此能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新回声消除系数(滤波器特性)。滤波器特性能够在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)内被更新。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算，基于所述立体声音频信号的和信号和差信号与由从所述一个或两个传声器收集的音频信号中减去相应的回声消除所得的各个回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算。进一步，基于所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱的计算，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导在所述两个或四个音频传送系统中，由合成多个适当的系统的传递函数而得到的多种传递函数合成估算误差，从而基于所述多种传递函数合成估算误差，推导所述两个或四个音频传送系统的传递函数的估算误差。另外，可以计算四个系统情况下的互频谱计算，例如，所述和信号与第一回声消除误差信号之间，所述和信号与第二回声消除误差信号之间，所述差信号与第一回声消除误差信号之间，和所述差信号与第二回声消除误差信号之间的各个互频谱。还有，所述传递函数合成估算误差可包括例如：第一传递函数合成估算误差，其是第一扬声器与第一传声器之间的传递函数的估算误差和第二扬声器与第一传声器之间的传递函数的估算误差的和；第二传递函数合成估算误差，其是第一扬声器与第一传声器之间的传递函数的估算误差和第二扬声器与第一传声器之间的传递函数的估算误差之间的差；第三传递函数合成估算误差，其是第一扬声器与第二传声器之间的传递函数的估算误差和第二扬声器与第二传声器之间的传递函数的估算误差的和；以及第四传递函数合成估算误差，其是每一扬声器与第二传声器之间的传递函数的估算误差和第二扬声器与第二传声器之间的传递函数的估算误差之间的差。推导所述四个音频传送系统的传递函数的估算误差的计算，可以包括例如：根据第一传递函数合成估算误差与第二传递函数合成估算误差的和，推导第一音频传送系统的传递函数的估算误差的计算；根据第一传递函数合成估算误差与第二传递函数合成估算误差之间的差，推导第二音频传送系统的传递函数的估算误差的计算；根据第三传递函数合成估算误差与第四传递函数合成估算误差的和，推导第三音频传送系统的传递函数的估算误差的计算；以及根据第三传递函数合成估算误差与第四传递函数合成估算误差之间的差，推导第四音频传送系统的传递函数的估算误差的计算。另外，检测所述立体声音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的传递函数分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给将由所述各个扬声器再生的相应信号，并从相应的所述一个或两个传声器收集的音频信号中减去所述回声消除信号，从而实现回声消除；另外，在此方法中，对所述立体声音频信号加以主要分量分析，生成彼此正交的两个不相关的合成信号；并利用一组所述两个不相关的合成信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数，由此设置相应的滤波器特性。根据本发明，因为通过对立体声音频信号加以主要分量分析生成的互相正交的两个信号是彼此不相关的，所以利用这样的两个信号，分别推导所述两个或四个音频传送系统的传递函数，由此设置相应的滤波器特性，能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述两个不相关的合成信号组作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的计算，可以是例如这样的计算，基于所述两个不相关的合成信号与各个传声器收集的音频信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的计算。进一步，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的计算，可以是例如这样的计算：推导所述两个不相关合成信号与各个传声器收集的音频信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导在所述两个或四个音频传送系统中由合成多个适当的系统所得的多种合成传递函数，从而基于所述多种合成传递函数，推导所述两个或四个音频传送系统的传递函数。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的传递函数分别被估算，以便设置相应的滤波器特性；通过将所述滤波器特性提供给所述各个扬声器再生的相应信号，并从相应的所述一个或两个传声器收集的音频信号中减去所述回声消除信号，从而实现回声消除；并且，在此方法中，对所述立体声音频信号加以主要分量分析，生成彼此正交的两个不相关的合成信号，并利用一组所述两个不相关的合成信号作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，因为通过对立体声音频信号加以主要分量分析生成的互相正交的两个信号是彼此不相关的，所以利用这样的两个信号分别推导所述两个或四个音频传送系统的传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除这个估算误差的值，因而能进行回声消除。与此相应，因为立体声音信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新回声消除系数(滤波器特性)。例如，滤波器特性能在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)内被更新。利用所述两个不相关合成信号组作为参考信号，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算，基于所述两个不相关合成信号与从所述一个或两个传声器收集的音频信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算。进一步，基于所述两个不相关合成信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的传递函数的估算误差的计算，可以是例如这样的计算：推导所述两个不相关合成信号与所述回声消除误差信号之间的互频谱，并就每一互频谱，在预定的时间周期中对它们作整体平均，以推导在所述两个或四个音频传送系统中合成多个适当的系统而得的多种传递函数合成估算误差，从而基于多种传递函数合成估算误差，推导所述两个或四个音频传送系统的传递函数的估算误差。在这种情况下，可以这样安排，即对除了由所述扬声器再生的声音以外的一些声音输入至所述传声器所致的通话重叠，进行检测，当检测到通话重叠时，使所述滤波器特性的更新周期相对地加长，反之，当未检测到通话重叠时，使所述滤波器特性的更新周期相对地缩短，这样，当通话重叠存在时，能充分地收敛估算误差，而当通话重叠不存在时，能加快估算误差的收敛。

本发明的立体声音频的传输的方法是这样一种方法，涉及两个空间，它们各形成所述四个音频传送系统，执行上述任何多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的立体声音频信号。与此相应，经过回声消除处理的立体声音频传输，能在两个地点之间进行，它能应用于例如电话会议系统或诸如此类的系统。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；提供给第一扬声器的音频信号分别经与第一和第二传声器相对应设置的第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号分别经与第一和第二传声器相对应设置的第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从第一传声器收集的音频信号中减去所述第一和第三回声消除信号而进行回声消除；并利用第二减法装置，从第二传声器收集的音频信号中减去所述第二和第四回声消除信号而进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于所述各个扬声器再生的立体声音频信号的和信号和差信号与所述各个传声器收集的音频信号之间的互频谱计算，分别推导与所述四个音频传送系统的传递函数相应的滤波器特性，从而将所推导的滤波器特性，分别设置至各个相应的所述第一至第第四滤波器装置。例如可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生从所述输入装置输入的立体声音频信号的和信号和差信号；和主信号传输系统，其用于不通过所述和/差信号产生装置，将所述输入装置输入的立体声音频信号传输至所述各个扬声器，其中所述传递函数计算装置，基于所述和/差信号产生装置产生的和信号和差信号与各个传声器收集的音频信号之间的互频谱计算，推导与所述四个音频传送系统的传递函数相应的滤波器特性，并将推导的滤波器特性设置至各个相应的所述第一至第四滤波器装置。另一方面，可以这样安排，立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生从所述输入装置输入的立体声音频信号的和信号和差信号；立体声音频信号解调装置，用于计算由所述和/差信号产生装置产生的和信号和差信号之间的和以及差，以便恢复原始立体声音频信号，其中由立体声音频信号解调装置恢复的立体声音频信号被传输至所述各个扬声器，并且，所述传递函数计算装置基于由所述和/差信号产生装置产生的和信号和差信号与个各个传声器收集的音频信号之间的互频谱计算，推导与所述四个音频传送系统相应的滤波器特性，并将所推导的滤波器特性分别设置至各个相应的所述第一至第四滤波器装置。另外，也可以这样安排，立体回声消除器包括：输入装置，其用于输入所述立体声音频信号的和信号和差信号；和立体声音频信号解调装置，其用于计算所述输入的和信号和差信号之间的和以及差，以便恢复原始的立体声音频信号，其中，由所述立体声音频信号解调装置恢复的立体声音频信号被传送至所述各个扬声器，并且，所述传递函数计算装置基于所述输入的和信号和差信号与各个传声器收集的音频信号之间的互频谱计算，推导与所述四个音频传送系统的传递函数相应的滤波器特性，并将推导的滤波器特性分别设置至各个相应的所述第一至第四滤波器装置。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；提供给第一扬声器的音频信号分别经与第一和第二传声器相对应设置的第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号分别经与第一和第二传声器相对应设置的第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从第一传声器收集的音频信号中减去所述第一和第三回声消除信号而进行回声消除；并利用第二减法装置，从第二传声器收集的音频信号中减去所述第二和第四回声消除信号而进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于将由所述各个扬声器再生的立体声音频信号的和信号和差信号与从所述传声器收集的音频信号中减去相应的回声消除信号而得的各个回声消除误差信号之间的互频谱计算，分别推导所述四个音频传送系统的传递函数的估算误差，由此分别将所述第一至第四滤波器装置的滤波器特性更新至消除所述估算误差的值。例如可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生从所述输入装置输入的立体音频信号的和信号和差信号；和主信号传输系统，其用于不经过所述和/差信号产生装置，将所述输入装置输入的立体声音频信号传输至所述各个扬声器，其中，所述传递函数计算装置基于由所述和/差信号产生装置产生的和信号和差信号与各个回声消除误差信号之间的互频谱计算，推导所述四个音频传送系统的传递函数的估算误差，并将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。另一种方法，可以这样安排，立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生由所述输入装置输入的立体声音频信号的和信号和差信号；和立体声音频信号解调装置，其用于计算由所述和/差信号产生装置产生的和信号和差信号之间的和以及差，以便恢复原始的立体声音频信号，其中，由所述立体声音频信号解调装置恢复的立体声音频信号传输至各个扬声器，并且，所述传递函数计算装置基于和/差信号产生装置产生的和信号和差信号与各个回声消除误差信号之间的互频谱计算，推导所述四个音频传送系统的传递函数的估算误差，并将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。或者，也可以这样安排，立体回声消除器包括：输入装置，其用于输入所述立体声音频信号的和信号和差信号；立体声音频信号解调装置，其用于计算所述输入的和信号和差信号之间的和以及差，以便恢复原始的立体声音频信号，其中，由立体声音频信号解调装置恢复的立体声音频信号被传输至所述各个扬声器，并且，传递函数计算装置基于由和/差信号产生装置产生的所述输入的和信号和差信号与各自回声消除误差信号之间的互频谱计算，推导所述四个音频传送系统的传递函数的估算误差，并将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；提供给第一扬声器的音频信号分别经与第一和第二传声器相对应设置的第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号分别经与第一和第二传声器相对应设置的第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从第一传声器收集的音频信号中减去所述第一和第三回声消除信号而进行回声消除；并利用第二减法装置，从第二传声器收集的音频信号中减去所述第二和第四回声消除信号而进行回声消除，所述立体回声消除器包括：传递函数计算装置，其用于基于对所述各个扬声器再生的立体声音频信号加以主要分量分析生成的互相正交的两个不相关的合成信号，与各个传声器收集的音频信号之间的互频谱计算，分别推导与所述四个音频传送系统的传递函数相应的滤波器特性，由此将所述推导的滤波器特性，分别设置至各个相应的第一至第四滤波器装置。例如可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；正交化装置，其用于对所述输入装置输入的立体声音频信号加以主要分量分析，生成互相正交的两个不相关的合成信号；和主信号传输系统，其用于不通过所述正交化装置，将输入装置输入的立体声音频信号传输至所述各个扬声器，其中，所述传递函数计算装置基于由所述正交化装置生成的两个不相关的合成信号与各自传声器收集的音频信号之间的互频谱计算，推导与所述四个音频传送系统的传递函数相应的滤波器特性，并将推导的滤波器特性分别设置至各个相应的所述第一至第四滤波器装置。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；提供给第一扬声器的音频信号，分别经与第一和第二传声器相对应设置的第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号，分别经与第一和第二传声器相对应设置的第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从第一传声器收集的音频信号中减去所述第一和第三回声消除信号而进行回声消除；并利用第二减法装置，从第二传声器收集的音频信号中减去所述第二和第四回声消除信号而进行回声消除，所述立体回声消除器包括：传递函数计算装置，其用于基于对所述各个扬声器再生的立体声音频信号加以主要分量分折生成的互相正交的两个不相关的合成信号，与由从所述两个传声器收集的音频信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述四个音频传送系统的传递函数的估算误差，由此将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。例如可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；正交化装置，其用于对所述输入装置输入的立体声音频信号加以主要分量分折，生成互相正交的两个不相关的合成信号；和主信号传输系统，用于不经过正交化装置，将输入装置输入的立体声音频信号传输至所述各个扬声器，其中所述传递函数计算装置基于由所述正交化装置生成的两个不相关的合成信号与各个回声消除误差信号之间的互频谱计算，推导所述四个音频传送系统的传递函数的估算误差，并将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。在这种情况下，可以这样安排，即设置通话重叠检测装置，其用于检测除了由所述扬声器再生的声音以外的一些声音输入至所述传声器造成的通话重叠，并且，当检测到通话重叠时，所述传递函数计算装置使所述滤波器特性的更新周期相对地加长，反之，当未检测到通话重叠时，使所述滤波器特性的更新周期相对地缩短，这样，当通话重叠存在时，能充分地收敛估算误差，而当通话重叠不存在时，能加快估算误差的收敛。

本发明的立体回声消除器，可进一步设置有相关检测装置，用于检测所述立体声音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音传送装置是这样的一种装置，涉及两个空间，它们各形成所述四个音频传送系统，在每一空间安排有前面所述立体回声消除器中的任何一种，以致在所述两个空间之间传输已用立体回声消除器消除回声的立体声音频信号。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的、彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给所述各个扬声器再生的单个信号的合成信号，分别生成回声消除信号；并且，从所述一个或多个传声器的单个收集的音频信号中减去所述回声消除信号，从而实现回声消除；另外，在方法中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间比所述多通道音频信号之间有较低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的多个合成传递函数，由此设置相应的滤波器特性。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数，设置相应的滤波器特性，因而能进行回声消除。与此相应，多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用多个低相关合成信号组作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的计算，可以是例如这样的计算，基于多个低相关合成信号与各个传声器的单个收集的音频信号之间的互频谱计算，分别推导所述多个音频传送系统的多个合成传递函数的计算。进一步，基于所述互频谱计算，分别推导所述多个音频传送系统的多个合成传递函数的计算，可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱；以及就每一互频谱，在预定的时间周期对它们作整体平均，分别推导所述多个音频传送系统的合成传递函数。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的、彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给所述各个扬声器再生的单个信号的合成信号，分别生成回声消除信号；并且，从所述一个或多个传声器的单个收集的音频信号中减去所述回声消除信号，从而实现回声消除；另外，在方法中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相应、且彼此之间比所述多通道音频信号之间有较低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的多个合成传递函数的估算误差，由此，将相应的滤波器特性更新至消除估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除所述估算误差的值，因此能进行回声消除。与此相应，多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新滤波器特性。利用多个低相关合成信号组作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号中减去所述回声消除信号所得到的回声消除误差信号之间的互频谱计算，分别推导多个音频传送系统的多个合成传递函数的估算误差所进行的计算。另外，基于所述互频谱计算，分别推导所述多个音频传送系统的合成传递函数的估算误差所进行的计算，可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间更低相关性的多个低相关合成信号；推导所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号中减去相应的回声消除信号所得的回声消除误差信号之间的互频谱；并就每个互频谱，在预定时间周期对它们作整个平均，以推导所述多个音频传送系统的合成传递函数的估算误差。另外，检测所述多个低相关合成信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的多通道声音传送方法是这样一种方法，涉及两个空间，它们各形成所述多个音频传送系统，在其中执行上述任何多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的多通道音频信号。因此，能在两个地点之间进行经过回声消除处理的多通道音频传输，它能应用于例如电话会议系统或诸如此类系统。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的多个合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给由各自扬声器要再生的各个信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号中减去所述回声消除信号，由此实现回声消除；并且，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数，由此设置相应的滤波器特性。根据本发明，因为立体声音频信号的和信号和差信号彼此之间具有低的相关性，所以用和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的多个合成传递函数，以便设置相应的滤波器特性，从而能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算，基于和信号和差信号与各个传声器的单个收集的音频信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算。进一步，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与各个传声器的单个收集的音频信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的多个合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由各自扬声器再生的各个信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号中，减去所述回声消除信号，由此实现回声消除，并且，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，因为立体声音频信号的和信号和差信号彼此之间具有低的相关性，所以用和信号和差信号作为参考信号，分别推导两个或四个音频传送系统的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，从而能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能够获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现自然的会话。利用所述立体声音频信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述立体声音频信号的和信号和差信号与从所述一个或两个传声器的单个收集的音频信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差。此外，基于所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差所进行的计算，可以是例如这样的计算：分别推导所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数的估算误差。另外，检测所述立体声音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音频的传输的方法是这样一种方法，涉及两个空间，它们各形成所述四个音频传送系统，执行上述任何多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的立体声音频信号。与此相应，能在两个地点之间进行经过回声消除处理的立体声音频传输，它能应用于例如会议系统或诸如此类的系统。

本发明的立体回声消除器是一种立体回声消除器，其中，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；将由所述各个扬声器再生的立体声音频信号的和信号分别经第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；将由所述各个扬声器再生的立体声音频信号的差信号分别经第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，通过由第一传声器收集的音频信号中减去第一和第三回声消除信号，进行回声消除；利用第二减去装置，通过从第二传声器收集的音频信号中减去所述第二和第四回声消除信号，进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于将由所述各个扬声器再生的立体音频信号的和信号和差信号与各个传声器收集的音频信号之间的互频谱计算，分别推导与所述四个音频传送系统的合成传递函数相对应的滤波器特性，从而将所述推导的滤波器特性，分别设置至各个相应的所述第一至第四滤波器装置。例如，可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生由所述输入装置输入的立体声音频信号的和信号和差信号；和主信号传输系统，其用于不经过所述和/差信号产生装置，将所述输入装置输入的立体声音频信号传输至所述各个扬声器，其中，所述传递函数计算装置基于由所述和/差信号产生装置产生的和信号和差信号与各个传声器收集的音频信号之间的互频谱计算，推导与所述四个音频传送系统的合成传递函数相对应的滤波器特性，并将推导的滤波器特性，分别设置至各个相应的所述第一至第四滤波器装置。

本发明的立体回声消除器是一种立体回声消除器，其中，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；将由所述各个扬声器再生的立体声音频信号的和信号分别经第一和第二滤波器装置作卷积计算，以便生成第一和第二回声消除信号；将由所述各个扬声器再生的立体声音频信号的差信号分别经第三和第四滤波器装置作卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，通过由第一传声器收集的音频信号中减去第一和第三回声消除信号，进行回声消除；利用第二减去装置，通过由第二传声器收集的音频信号中减去所述第二和第四回声消除信号，进行回声消除；所述立体回声消除器包括：传递函数计算装置，其基于将由所述各个扬声器再生的立体声音频信号的和信号和差信号与从所述两个传声器收集的音频信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述四个音频传送系统的合成传递函数的估算误差，从而，将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。例如，可以这样安排，即本发明的立体回声消除器包括：输入装置，其用于输入所述立体声音频信号；和/差信号产生装置，其用于产生由输入装置输入的立体声音频信号的和信号和差信号；和主信号传输系统，其用于不经过所述和/差信号产生装置，将从所述输入装置输入的立体声音频信号传输至所述各个扬声器，其中，所述传递函数计算装置基于由所述和/差信号产生装置产生的和信号和差信号与各个回声消除信号之间的互频谱计算，推导所述四个音频传送系统的合成传递函数的估算误差，并将所述第一至第四滤波器装置的滤波器特性分别更新至消除估算误差的值。

本发明的立体回声消除器可进一步被提供以相关检测装置，相关检测装置用于检测所述立体声音频信号的和信号和差信号，当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音传送装置是这样一种装置，涉及两个空间，它们各形成所述四个音频传送系统，在每个空间配置前面所述立体回声消除器中的任何一种，以致在所述两个空间之间传输已用所述立体回声消除器消除回声的立体音频信号。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或两个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的并且彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给所述各个扬声器再生的单个信号，分别生成回声消除信号；并且，从所述一个或多个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，从而实现回声消除；另外，在这个方法中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组所述多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的合成传递函数，由此设置相应的滤波器特性。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的合成传递函数，设置相应的滤波器特性，从而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述一组多个低相关合成信号作为参考信号分别推导所述多个音频传送系统的合成传递函数的计算，可以是例如这样的计算，基于多个低相关合成信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱计算，分别推导多个音频传送系统的合成传递函数的计算。此外，基于所述互频谱计算，分别推导所述多个音频传送系统的合成传递函数的计算，可以是例如这样的计算：通过相加和相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱；并就每一互频谱，在预定的时间周期对它们做整体平均，以推导所述多个音频传送系统的合成传递函数。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或两个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的并且彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各个扬声器再生的单个信号，分别生成回声消除信号；并且，从所述一个或多个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，从而实现回声消除；另外，在这个方法中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间有较低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的合成传递函数的估算误差，从而将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，从而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新滤波器特性。利用多个低相关合成信号组作为参考信号，分别推导所述多个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号所得的回声消除误差信号之间的互频谱计算，分别推导所述多个音频传送系统的合成传递函数的估算误差的计算。另外，基于所述互频谱计算，分别推导所述多个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号以生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号所得的回声消除误差信号之间的互频谱，并就每个互频谱，在预定的时间周期对它们作整体平均，以推导所述多个音频传送系统的合成传递函数的估算误差。另外，检测所述多个低相关合成信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的多通道声音传送方法是这样的一种方法，涉及两个空间，它们各形成所述多个音频传送系统，执行上述任何一种多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的多通道音频信号。与此相应，能在两个地点之间进行经过回声消除处理的多通道音频信号传输，它能应用于例如电话会议系统或诸如此类的系统。

本发明的立体回声消除方法是这样的一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的多个合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各个扬声器再生的单个信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，由此实现回声消除；并且，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数，由此设置相应的滤波器特性。根据本发明，因为立体声音频信号的和信号和差信号彼此之间具有低的相关性，所以用和信号和差信号信号作为参考信号，分别推导两个或四个音频传送系统的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除所述估算误差的值，从而能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能够获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算，基于和信号和差信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算。此外，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数。

本发明的立体回声消除方法是这样的一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各个扬声器再生的单个信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，由此实现回声消除；并且，在此方法中，利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间更低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，从而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现自然的会话。也能实时地更新滤波器特性。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述立体声音频信号的和信号和差信号与从所述一个或两个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号而得到的各个回声消除误差信号之间的互频谱计算。此外，基于所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱，并就每一互频谱在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数的估算误差。另外，检测所述立体音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音传输方法是这样一种方法，涉及两个空间，它们各形成所述四个音频传送系统，执行上述任何一种多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的立体声音频信号。与此相应，能够在两个地点之间进行经过回声消除处理的立体声音频的传输，它能应用于例如电话会议系统或诸如此类的系统。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；提供给第一扬声器的音频信号经由第一和第二滤波器装置分别进行的卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号经由第三和第四滤波器装置分别进行的卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从各个传声器收集的音频信号的和信号中减去所述第一和第三回声消除信号，进行回声消除；并利用第二减法装置，从各个传声器收集的差信号中减去所述第二和第四回声消除信号进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于由所述各个扬声器再生的立体声音频信号的和信号和差信号与各个传声器收集的音频信号的和信号和差信号之间的互频谱计算，分别推导与所述四个音频传送系统的合成传递函数相对应的滤波器特性，从而分别将所述推导的滤波器特性设置至各个相应的所述第一至第四滤波器装置。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，在其中设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；提供给第一扬声器的音频信号经由第一和第二滤波器装置分别进行卷积计算，以便生成第一和第二回声消除信号；提供给第二扬声器的音频信号经分别由第三和第四滤波器装置进行的卷积计算，以便生成第三和第四回声消除信号；利用第一减法装置，从各个传声器收集的音频信号的和信号中减去所述第一和第三回声消除信号，进行回声消除；并利用第二减法装置，从各个传声器收集的差信号中减去所述第二和第四回声消除信号进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于由所述各个扬声器再生的立体音频信号的和信号和差信号与从各个传声器收集的音频信号的和信号和差信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述四个音频传送系统的合成传递函数的估算误差，从而将所述第一至第四滤波器装置的滤波器特性分别更新至消除所述估算误差的值。

本发明的立体回声消除器还设置有相关检测装置，相关检测装置用于检测所述立体声音频信号的和信号和差信号，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音频信号传送装置是这样一种装置，涉及两个空间，它们各形成所述四个音频传送系统，在每一空间配置前面所述任何一种立体回声消除器，以致在两个空间之间进行已用立体回声消除器消除回声的立体声音频信号的传输。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或两个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的并且彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各自扬声器再生的单个信号的合成信号，分别生成回声消除信号；以及，从所述一个或多个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，从而实现回声消除；在其中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的合成传递函数，从而设置相应的滤波器特性。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的合成传递函数，并设置相应的滤波器特性，从而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用多个低相关合成信号组作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的计算，可以是例如这样的计算，基于多个低相关合成信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱计算，分别推导多个音频传送系统的多个合成传递函数的计算。进一步，基于所述互频谱计算，分别推导所述多个音频传送系统的合成传递函数的计算，可以是例如这样的计算：通过相加和相减，合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱；并就每一互频谱，在预定的时间周期对它们做整体平均，以推导所述多个音频传送系统的合成传递函数。

本发明的多通道回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有多个扬声器和一个或两个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的并且彼此有相关性的多通道声音被所述传声器收集；所述多个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给所述各自扬声器再生的单个信号的合成信号，分别生成回声消除信号；并且，从所述一个或多个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，从而实现回声消除；在其中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，分别推导所述多个音频传送系统的多个合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，从而能进行回声消除。与此相应，多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新滤波器特性。利用所述一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的多个合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号所得到的回声消除误差信号之间的互频谱计算，分别推导所述多个音频传送系统的合成传递函数的估算误差的计算。另外，基于所述互频谱计算，分别推导所述多个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算：通过相加或相减，合成所述多通道音频信号，以生成彼此之间比所述多通道音频信号之间有较低相关性的多个低相关合成信号；推导所述多个低相关合成信号与从所述一个或多个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号而得的回声消除误差信号之间的互频谱；以及就每个互频谱，在预定的时间周期对它们作整体平均，以推导所述多个音频传送系统的合成传递函数的估算误差。另外，检测所述多个低相关合成信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的多通道音频信号传送方法是这样一种方法，涉及两个空间，它们各形成所述多个音频传送系统，执行上述任何一种多通道回声消除方法，以致在所述两个空间之间传输已用所述方法作过回声消除的多通道音频信号。与此相应，能够在两个地点之间进行经过回声消除处理的多通道音频的传输，它能应用于例如会议系统或诸如此类系统。

本发明的立体回声消除方法是这样一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给由所述各个扬声器再生的单个信号的合成信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，从而实现回声消除；其中利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数，由此设置相应的滤波器特性。根据本发明，因为立体声音频信号的和信号和差信号彼此之间具有低的相关性，所以用和信号和差信号信号作为参考信号，分别推导两个或四个音频传送系统的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除所述估算误差的值，从而能进行回声消除。与此相应，因为立体声信号能不经过或很少经过作用于立体声信号而引起其恶化的处理，而从扬声器再生，所以能够获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算，基于和信号和差信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算。此外，基于所述互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与各个传声器的单个收集的音频信号的合成信号之间的互频谱，以及就每一互频谱，在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数。

本发明的立体回声消除方法是这样的一种方法，其中，涉及一个空间，在该空间内设置有两个扬声器和一个或两个传声器，形成两个或四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述传声器收集；所述两个或四个音频传送系统的合成传递函数分别被估算，以便设置相应的滤波器特性；通过将所述设置的滤波器特性提供给将由所述各个扬声器再生的单个信号的合成信号，分别生成回声消除信号；并从所述一个或两个传声器的单个收集的音频信号的合成信号中减去所述回声消除信号，由此实现回声消除；并且，其中利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的多个合成传递函数的估算误差，由此将相应的滤波器特性更新至消除所述估算误差的值。根据本发明，利用与适当地合成彼此有相关性的多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号，分别推导所述多个音频传送系统的合成传递函数的估算误差，以便逐次地将相应的滤波器特性更新至消除估算误差的值，从而能进行回声消除。与此相应，由于多通道信号能不经过或很少经过作用于多通道信号而引起其恶化的处理，而从扬声器再生，所以能获得优良的再生音调质量。另外，在再生信号中没有或只有很少延迟。因此，当应用于电话会议系统或诸如此类系统时，能实现很自然的会话。也能实时地更新滤波器特性。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算，基于所述立体声音频信号的和信号和差信号与从所述一个或两个传声器的单个收集的音频信号的合成信号中减去相应的回声消除信号所得到的各个回声消除误差信号之间的互频谱计算的计算。此外，基于所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱计算，分别推导所述两个或四个音频传送系统的合成传递函数的估算误差的计算，可以是例如这样的计算：推导所述立体声音频信号的和信号和差信号与所述回声消除误差信号之间的互频谱，并就每一互频谱在预定的时间周期对它们作整体平均，以推导所述两个或四个音频传送系统的合成传递函数的估算误差。另外，检测所述立体音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音传输方法是这样一种方法，涉及两个空间，它们各形成所述四个音频传送系统，执行上述任何一种多通道回声消除方法，以致在所述两个空间之间传输已用所述方法消除回声的立体声音频信号。与此相应，能够在两个地点之间进行经过回声消除处理的立体声音频信号的传输，它能应用于例如电话会议系统或诸如此类的系统。

本发明的立体回声消除器是这样一种立体回声消除器，其中，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；由所述各个扬声器再生的立体声音频信号的和信号分别经第一和第二滤波器装置作卷积计算，结果生成第一和第二回声消除信号；由所述各个扬声器再生的立体音频信号的差信号分别经第三和第四滤波器装置作卷积计算，结果生成第三和第四回声消除信号；利用第一减法装置，通过从各个传声器收集的音频信号的和信号中减去所述第一和第三回声消除信号，进行回声消除；利用第二减法装置，通过从各个传声器收集的音频信号差信号中减去所述第二和第四回声消除信号，进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于将由所述各自扬声器再生的立体音频信号的和信号和差信号与各个传声器收集的音频信号的和信号和差信号之间的互频谱计算，分别推导与所述四个音频传送系统的合成传递函数相应的滤波器特性，由此，将所述推导的滤波器特性分别设置至各个相应的所述第一至第四滤波器装置。

本发明的立体回声消除器是一种立体回声消除器，其中，涉及一个空间，该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由各个扬声器再生的立体声音被各个传声器收集；由各个扬声器再生的立体声音频信号的和信号分别经第一和第二滤波器装置作卷积计算，结果生成第一和第二回声消除信号；由所述各个扬声器再生的立体音频信号的差信号分别经第三和第四滤波器装置作卷积计算，结果生成第三和第四回声消除信号；利用第一减法装置，通过从各个传声器收集的音频信号中减去所述第一和第三回声消除信号，进行回声消除；同时，利用第二减法装置，从各个传声器收集的音频信号差信号中减去所述第二和第四回声消除信号，进行回声消除；所述立体回声消除器包括：传递函数计算装置，其用于基于将由所述各个扬声器再生的立体音频信号的和信号和差信号与从所述各个传声器收集的音频信号的和信号和差信号中减去相应的回声消除信号所得的各个回声消除误差信号之间的互频谱计算，分别推导所述四个音频传送系统的合成传递函数的估算误差，从而，将所述第一至第四滤波器装置的滤波器特性更新至消除所述估算误差的值。

本发明的立体回声消除器，可进一步设置有相关检测装置，相关检测装置用于检测所述立体声音频信号的和信号和差信号之间的相关性，并当所述相关值不小于预定值时，停止更新所述滤波器特性，以防止回声消除误差信号意外地增加。

本发明的立体声音频信号传送装置是这样一种装置，涉及两个空间，它们各形成所述四个音频传送系统，在每一空间配置上述任何一种所述立体回声消除器，以致在两个空间之间进行已用立体回声消除器消除回声的立体声音频信号的传输。

本发明的传递函数计算装置是这样一种传递函数计算装置，涉及一个空间，在该空间内设置有多个扬声器和一个或多个传声器，形成多个音频传送系统，在音频传送系统中，由所述各个扬声器再生的并且彼此有相关性的多通道声音被所述传声器收集；估算所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，其中，利用与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号作为参考信号(例如，通过适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号，并用一组所述多个低相关合成信号作为参考信号，或者直接输入与通过适当地合成所述多通道音频信号所得的信号相对应、且彼此之间具有比所述多通道音频信号之间低的相关性的一组多个低相关合成信号，并用这组所述多个低相关合成信号作为参考信号，或诸如此类)，估算各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。利用多个低相关合成信号组作为参考信号，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算，基于多个低相关合成信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成信号之间的互频谱计算，分别推导各个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算。此外，基于所述互频谱计算，分别推导所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算：分别推导所述多个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，通过相加或相减，适当地合成所述多通道音频信号，生成彼此之间具有比所述多通道音频信号之间低的相关性的多个低相关合成信号；推导所述多个低相关合成信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱，并就每一互频谱，在预定的时间周期对它们作整体平均。基于所述互频谱计算，分别推导所述多个音频传送系统的单个传递函数的计算，也可以是这样的计算：通过对所述多通道音频信号加以主要分量分析，生成多个互相正交的不相关的合成信号，分别推导所述多个音频传送系统的单个传递函数；推导所述多个不相关的合成信号与传声器的单个收集的音频信号之间的互频谱；以及就每一互频谱，在预定的时间周期对它们作整体平均。

本发明的传递函数计算装置是这样一种传递函数计算装置，涉及一个空间，该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；估算所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数，其中，利用所述立体声音频信号的和信号和差信号作为参考信号，估算所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得到的多个合成传递函数。利用所述立体声音频信号的和信号和差信号作为参考信号，分别推导所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是例如这样的计算，基于所述和信号和所述差信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得到的多个合成信号之间的互频谱计算，分别推导所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得到的多个合成传递函数的计算。进一步，基于所述互频谱计算，分别推导所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数的计算，可以是这样的计算：推导所述立体声音频信号的和信号和差信号与传声器的单个收集的音频信号或通过适当地合成所述单个收集的音频信号所得的多个合成传递函数之间的互频谱，以及就每一互频谱，在预定的时间周期对它们作整体平均，由此分别推导所述四个音频传送系统的单个传递函数或通过适当地合成所述单个传递函数所得的多个合成传递函数。

本发明的传递函数计算装置是这样一种传递函数计算装置，涉及一个空间，在该空间内设置有两个扬声器和两个传声器，形成四个音频传送系统，在音频传送系统中，由所述各个扬声器再生的立体声音被所述各个传声器收集；估算所述四个音频传送系统的单个传递函数，其中，通过对所述立体声音频信号加以主要分量分析，生成互相正交的两个不相关的合成信号；并利用一组所述两个不相关的合成信号作为参考信号，估算所述四个音频传送系统的单个传递函数。利用所述两个不相关的合成信号作为参考信号，分别推导所述四个音频传送系统的单个传递函数的计算，可以是例如这样的计算，基于所述两个不相关的合成信号与传声器的单个收集的音频信号之间的互频谱计算。进一步，基于所述互频谱计算，分别推导所述四个音频传送系统的单个传递函数的计算，可以是例如这样的计算：推导所述两个不相关的合成信号与传声器的单个收集的音频信号之间的互频谱；以及就每一互频谱，在预定的时间周期对它们作整体平均，从而分别推导所述四个音频传送系统的单个传递函数。在这种情况下，当检测到通话重叠时，即除了由所述扬声器再生的声音以外的一些声音也输入至所述传声器时，将进行整体平均的时间周期相对地加长，而当未检测到通话重叠时，将整体平均的时间周期相对缩短，这样，当通话重叠存在时，能充分地收敛估算误差，而当通话重叠不存在时，加快估算误差的收敛。

在本发明中，“音频信号”不限于人的声音，而是复盖音频频带中的全部音频信号。

附图说明

图1示出图2的立体回声消除器16、24中一种结构示例的方块图。

图2示出本发明的立体声音传送装置实施例的方块图。

图3示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图4示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图5示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图6示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图7示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图8示出图1的立体回声消除器16、24的回声消除性能的模拟测量结果曲线图。

图9示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图10示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图11示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图12示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图13示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图14示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图15示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图16示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图17示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图18示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图19示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图20示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图21示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图22示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图23示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图24示出图1、9至10的结构修改的方块图。

图25示出图12至15的结构修改的方块图。

图26示出图16至19的结构修改的方块图。

图27示出图20和21的结构修改的方块图。

图28示出图2的立体回声消除器16、24中的另一结构示例的方块图。

图29示出在图28立体回声消除器中，正交处理和推导脉冲响应或其估算误差用的单位时间间隔示例的时序图。

图30示出图28的正交化滤波器500的功能块示例图。

图31示出图28的传递函数计算502的示例图。

图32示出图31的传递函数计算502的修改功能块图。

图33示出当没有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图34示出当没有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图35示出当没有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图36示出当没有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图37示出当有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图38示出当有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图39示出当有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图40示出当有通话重叠时，对图28的立体回声消除器16、24的回声消除量时域变化的模拟测量结果曲线图。

图41示出当没有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图42示出当没有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图43示出当没有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图44示出当没有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图45示出当有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图46示出当有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图47示出当有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图48示出当有通话重叠时，对图28的立体回声消除器16、24的传递函数估算误差时域变化的模拟测量结果曲线图。

图49示出根据重迭处理的整体平均处理的示例图。

图50示出图28的立体回声消除器16、24的修改方块图。

图51示出图1中的传声器数目修改为一个时的结构示例方块图。

具体实施方式

下面，将描述本发明的实施例。图2示出根据本发明的两个方向立体声音传送装置的整体结构。这种装置执行地点A和地点B之间的两个方向立体传输，应用于例如电话会议系统。在地点A，两个扬声器SP-A(L)和SP-A(R)以及两个传声器MC-A(L)和MC-A(R)安排在一个空间内。传声器MC-A(L)和MC-A(R)收集的音频信号，在A/D转换器12和14中分别转换为数字信号，并在立体回声消除器16中加以回声消除处理，然后在CODEC(编码器和译码器)18中调制后，通过有线或无线电传输线路20传送至地点B。在地点B，两个扬声器SP-B(L)和SP-B(R)以及两个传声器MC-B(L)和MC-B(R)安排在一个空间内。顺便说一下，假定一种情况是：参与者在地点A和B以面对面的方式，通过电话会议系统通话，扬声器和传声器的安排是：由一个地点的参与者的左边的传声器收集的声音，从另一地点的参与者右边的扬声器再生，相反，由一个地点的参与者右边的传声器收集的声音，从另一地点的参与者左边的扬声器再生。具体地说，当扬声器SP-A(L)和传声器MC-A(L)安排在地点A的参与者左边，而扬声器SP-A(R)和传声器MC-A(R)安排在地点A的参与者右边时，则扬声器SP-B(L)和传声器MC-B(L)安排在地点B参与者右边，而扬声器SP-B(R)和传声器MC-B(R)安排在地点B的参与者左边。

从地点A发送的信号输入至CODEC 22，这里，传声器MC-A(L)和MC-A(R)收集的音频信号被解调。这些被解调的传声器MC-A(L)和MC-A(R)收集的信号，经过立体回声消除器24，在D/A转换器26和28中分别转换为模拟信号，并在扬声器SP-B(L)和SP-B(R)分别再生。在地点B传声器MC-B(L)和MC-B(R)收集的音频信号，分别在A/D转换器30和32中转换为数字信号，在立体回声消除器24作回声消除处理，然后在CODEC 22中被调制，并通过传输线路20发送至地点A。发送至地点A的信号输入至CODEC 18，在那里，传声器MC-B(L)和MC-B(R)收集的音频信号被解调。这些被解调的传声器MC-B(L)和MC-B(R)收集的音频信号，经过立体回声消除器16，在D/A转换器34和36分别转换为模拟信号，并在扬声器SP-A(L)和SP-A(R)分别再生。

图1示出立体回声消除器16、24的结构示例。从配对侧发送来的、输入至线路输入端LI(L)和LI(R)的左/右两通道立体信号x_L和x_R，照原样(即不通过和/差信号产生装置52)，并分别在扬声器SP(L){表示SP-A(L)或SP-B(L)}和SP(R){表示SP-A(R)或SP-B(R)}再生。

滤波器装置40-1至40-4由例如FIR滤波器形成。它们当中，滤波器装置40-1设置的脉冲响应，与扬声器SP(L)和传声器MC(L){表示MC-A(L)或MC-B(L)}之间的传递函数相对应，用这个脉冲响应，进行将从声音输出端SO(L)输出的信号的卷积计算，由此产生回声消除信号EC1，其与这样得到的一个信号相对应，这个信号从声音输出端SO(L)输出，在扬声器SP(L)再生，由传声器MC(L)收集并输入至声音输入端SI(L)。滤波器装置40-2设置的脉冲响应与扬声器SP(L)和传声器MC(R){表示MC-A(R)或MC-B(R)}之间的传递函数相对应。用这个脉冲响应，进行从声音输出端SO(L)输出的信号的卷积计算，由此产生回声消除信号EC2，其与这样得到的一个信号相对应，这个信号从声音输出端SO(L)输出，在扬声器SP(L)再生，由传声器MC(L)收集并输入至声音输入端SI(R)。滤波器装置40-3设置的脉冲响应，与扬声器SP(R)和传声器MC(L)之间的传递函数相对应。用这个脉冲响应，进行从声音输出端SO(R)输出的信号的卷积计算，由此产生回声消除信号EC3，其与这样得到的一个信号相对应，这个信号从声音输出端SO(R)输出，在扬声器SP(R)再生，由传声器MC(L)收集并输入至声音输入端SI(L)。滤波器装置40-4设置的脉冲响应，与扬声器SP(R)和传声器MC(R)之间的传递函数相对应。用这个脉冲响应，进行从声音输出端SO(R)输出的信号的卷积计算，由此产生回声消除信号EC4，其与这样得到的一个信号相对应，这个信号从声音输出端SO(R)输出，在扬声器SP(R)再生，由传声器MC(R)收集并输入至声音输入端SI(R)。

加法器44进行EC1+EC3的计算。加法器46进行EC2+EC4的计算。减法器48从由声音输入端SI(L)输入的传声器MC(L)收集的音频信号中减去回声消除信号EC1+EC3，由此进行回声消除。减法器50从由声音输出端SI(R)输入的传声器MC(R)收集的音频信号中减去回声消除信号EC2+EC4，由此进行回声消除。各个左和右通道的这些回声消除信号，分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点发送。

和/差信号产生装置52利用加法器54，对输入至线路输入端LI(L)和LI(R)的左/右两个通道立体信号x_L和x_R进行相加，以产生和信号x_L+x_R，利用减法器56对它们进行相减，以产生差信号x_L-x_R，和信号x_L+x_R，和差信号x_L-x_R它们之间的相关一般是低的，经常是大致不相关。传递函数计算装置58执行由52产生的和信号x_L和x_R和差信号x_L-x_R与由减法器48和50输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，设置滤波器装置40-1至40-4的滤波器特性(脉冲响应)。具体地说，刚启动系统时，滤波器装置40-1至40-4的滤波器特性未被设置，即系数全部设置为零，以致回声消除信号EC1至EC4是零，因此，传声器MC(L)和MC(R)它们本身的收集的音频信号从减法器48和50输出。所以在这时，传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_L+x_R和差信号x_L-x_R，与由减法器48和50输出的传声器MC(L)和MC(R)收集的音频信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，分别推导扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数，并执行将滤波器装置40-1至40-4的滤波器特性初始设置为与这些传递函数相对应的值。在初始设置之后，因为回声消除信号是由滤波器装置40-1至40-4产生的，所以，回声消除误差信号e_L和e_R，其相应于传声器MC(L)和MC(R)收集的音频信号与回声消除信号EC1至EC4之间的差信号从减法器48和50输出。因此在这时，传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_L+x_R和差信号x_L-x_R与由减法器48和50输出的回声消除误差信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，分别推导扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数，并将滤波器装置40-1至40-4的滤波器特性分别更新为消除估算误差的值。通过每一规定时间周期的重复这种更新操作，回声消除误差能收敛到最小值，另外，即使传递函数由于传声器位置之类的移动而变化，通过不断地更新有关于它的滤波器装置40-1至40-4的滤波器特性，回声消除误差也能收敛到最小值。

相关检测装置60基于相关值计算或诸如此类计算，检测和信号x_L+x_R和与差信号x_L-x_R之间的相关，当相关值不小于规定值时，停止更新上述滤波器特性。当相关值变得低于规定值时，则重新启动上述滤波器特性的更新。顺便说一下，关于用于推导和信号x_L+x_R与差信号x_L-x_R之间的相关的具体技术，任何检测两个信号相关的已知技术都可使用。

这里，将描述通过传递函数计算装置58为滤波器装置40-1至40-4设置滤波器特性(脉冲响应)。在下面的描述中，传递函数和滤波器特性用下述符号表示。

·Hxx：传递函数(频率轴表示)。

·hxx：对应于Hxx的脉冲响应(时间轴表示)。

·H^xx：估算的传递函数(给滤波器设置的传递函数)。

·h^xx：对应于H^xx的脉冲响应。

·ΔHxx：传递函数估算误差。

·Δhxx：对应于ΔHxx的脉冲响应。

(注意：适合的符号分派给xx。)

和信号和差信号分别定义如下。

和信号：x_M(＝x_L+x_R)

差信号：x_S(＝x_L-x_R)

扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间四个音频传送系统的传递函数分别定义如下。

H_LL：从扬声器SP(L)至传声器MC(L)的系统传递函数。

H_LR：从扬声器SP(L)至传声器MC(R)的系统传递函数。

H_RL：从扬声器SP(R)至传声器MC(L)的系统传递函数。

H_RR：从扬声器SP(R)至传声器MC(R)的系统传递函数。

立体回声消除器16、24的线路输入端LI(L)和LI(R)中的输入信号x_L和x_R，用下式代替

x_L＝(x_M+x_S)/2

x_R＝(x_M-x_S)/2

然后，在传递函数计算装置58中进行下面所示的计算。从减法器48和50输出的信号e_L和e_R{在滤波器装置40-1至40-4初始设置之前，他们其实是MC(L)和MC(R)收集的音频信号，而在初始设置之后，它们是回声消除误差信号}，假定信号x_M、x_S、e_L和e_R的频率轴表示，分别给定为X_M、X_S、E_L和E_R，成为

E_L＝{(X_M+X_S)HLL/2}+{(X_M-X_S)H_RL/2}

E_R＝{(X_M+X_S)H_LR/2}+{(X_M-X_S)H_RR/2}

因此

2E_L＝(X_M+X_S)H_LL+(X_M-X_S)H_RL

   ＝X_M(H_LL+H_RL)+X_S(H_LL-H_RL)…(1)

2E_R＝(X_M+X_S)H_LR+(X_M-X_S)H_RR

   ＝X_M(H_LR+H_RR)+X_S(H_LR-H_RR)…(2)

当等式(1)两边乘以X_M和X_S的复数共轭X_M^*和X_S^*(即推导互频谱计算)以及整体平均时，分别得到

SX_M^*·2E_L＝SX_M^*·X_M(H_LL+H_RL)+SX_M^*·X_S(H_LL-H_RL)…(3)

SX_S^*·2E_L＝SX_S^*·X_M(H_LL+H_RL)+SX_S^*·X_S(H_LL-H_RL)…(4)类似地，当等式(2)两边乘以X_M和X_S的复数共轭X_M^*和X_S^*时，分别得到

SX_M^*·2E_R＝SX_M^*·X_M(H_LR+H_RR)+SX_M^*·X_S(H_LR-H_RR)…(5)

SX_S*·2E_R＝SX_S^*·X_M(H_LR+H_RR)+SX_S^*·X_S(H_LR-H_RR)…(6)

在等式(3)至(6)中，因为X_M和X_S相互大致不相关，所以具有X_M^*·X_S或X_S^*·X_M的项在整体平均时变成近似于零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(3)至(6)分别变成

SX_M^*·2E_L＝S|X_M|²(H_LL+H_RL)…(3′)

SX_S^*·2E_L＝S|X_S|²(H_LL-H_RL)…(4′)

SX_M^*·2E_R＝S|X_M|²(H_LR+H_RR)…(5′)

SX_S^*·2E_R＝S|X_S|²(H_LR-H_RR)…(6′)

通过变换等式(3′)至(6′)，分别导出两个传递函数合成形式中的下述每一合成传递函数

H_LL+H_RL＝SX_M^*·2E_L/S|X_M|²…(3″)

H_LL-H_RL＝SX_S^*·2E_L/S|X_S|²…(4″)

H_LR+H_RR＝SX_M^*·2E_R/S|X_M|²…(5″)

H_LR-H_RR＝SX_S^*·2E_R/S|X_S|²…(6″)

当等式(3″)和(4″)的相应边相加在一起时，得到

H_LL＝(SX_M^*·E_L/S|X_M|²)+(SX_S^*·E_L/S|X_S|²)…(7)

当等式(3″)和(4″)的相应边之间相减时，得到

H_RL＝(SX_M^*·E_L/S|X_M|²)-(SX_S^*·E_L/S|X_S|²)…(8)

当等式(5″)和(6″)的相应边相加在一起时，得到

H_LR＝(SX_M^*·E_R/S|X_M|²)+(SX_S^*·E_R/S|X_S|²)…(9)

当等式(5″)和(6″)的相应边之间相减时，得到

H_RR＝(SX_M^*·E_R/S|X_M|²)-(SX_S^*·E_R/S|X_S|²)…(10)

因此，传递函数H_LL、H_RL、H_LR和H_RR分别被推导出来。对这些推导出来的传递函数加以逆付里叶变换所得的脉冲响应h_LL、h_RL、h_LR和h_RR，分别是将设置到滤波器装置40-1、40-2、40-3和40-4的滤波器特性。因此，传递函数计算装置58基于输入的信号x_M、差信号x_S和输入的减法器48和50的输出信号e_L和e_R，从等式(7)至(10)得出各个传递函数H_LL、H_RL、H_LR和H_RR，通过对这些所得传递函数加以逆付里叶变换，推导脉冲响应h_LL、h_RL、h_LR和h_RR，并分别将所得的脉冲响应设置给滤波器装置40-1、40-2、40-3和40-4，作为h^_LL、h^_RL、h^_LR和h^_RR，并进一步通过在每一适当确定的规定时间周期(即进行整体平均的时间周期)中重复这种计算，更新脉冲响应。

当上述脉冲响应更新技术利用估算误差参数来说明时，其过程变成如下。传声器MC(L)和MC(R)的输出信号(所收集的音频信号)，假定x_L、x_R、y_L和y_R分别给定为X_L、X_R、Y_L和Y_R，则变成

Y_L＝(X_M+X_S)H_LL/2+(X_M-X_S)H_RL/2

Y_R＝(X_M+X_S)H_LR/2+(X_M-X_S)H_RR/2

从减法器48输出的信号E_L变成

E_L＝Y_L-(X_L·H^_LL+X_R·H^_RL)

  ＝{(X_M+X_S)H_LL/2+(X_M-X_S)H_RL/2}-{(X_M+X_S)H^_LL/2+(X_M-X_S)H^_RL/2}

因此

2E_L＝X_M(H_LL+H_RL-H^_LL-H^_RL)+X_S(H_LL-H_RL-H^_LL+H^_RL)…(11)

从减法器50输出的信号E_R变成

E_R＝Y_R-(X_L·H^_LR+X_R·H^_RR)

  ＝{(X_M+X_S)H_LR/2+(X_M-X_S)H_RR/2}-{(X_M+X_S)H^_LR/2+(X_M-X_S)H^_RR/2}

因此

2E_R＝X_M(H_LR+H_RR-H^_LR-H^_RR)+X_S(H_LR-H_RR-H^_LR+H^_RR)…(12)

当给出如下估算误差

ΔH_LL＝H_LL-H^_LL

ΔH_RL＝H_RL-H^_RL

ΔH_LR＝H_LR-H^_LR

ΔH_RR＝H_RR-H^_RR时，等式(11)和(12)变成：

2E_L＝X_M(ΔH_LL+ΔH_RL)+X_S(ΔH_LL-ΔH_RL)…(11′)

2E_R＝X_M(ΔH_LR+ΔH_RR)+X_S(ΔH_LR-ΔH_RR)…(12′)

当等式(11′)的两边乘以X_M和X_S的复数共轭X_M^*和X_S^*(即计算互频谱计算)并整体平均时，则分别得到

SX_M^*·2E_L＝SX_M^*·X_M(ΔH_LL+ΔH_RL)+SX_M^*·X_S(ΔH_LL-ΔH_RL)…(13)

SX_S^*·2E_L＝SX_S^*·X_M(ΔH_LL+ΔH_RL)+SX_S^*·X_S(ΔH_LL-ΔH_RL)…(14)类似的，当等式(12′)两边乘以X_M和X_S的复数共轭X_M^*和X_S^*时，分别得到

SX_M^*·2E_R＝SX_M^*·X_M(ΔH_LR+ΔH_RR)+SX_M^*·X_S(ΔH_LR-ΔH_RR)…(15)

SX_S^*·2E_R＝SX_S^*·X_M(ΔH_LR+ΔH_RR)+SX_S^*·X_S(ΔH_LR-ΔH_RR)…(16)

在等式(13)至(16)中，因为X_M和X_S相互间大致不相关，所以当整体平均时，具有X_M^*·X_S或X_S^*·X_M的项近似地变成零。进一步

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(13)至(16)分别变成

SX_M^*·2E_L＝S|X_M|²(ΔH_LL+ΔH_RL)…(13′)

SX_S^*·2E_L＝S|X_S|²(ΔH_LL+ΔH_RL)…(14′)

SX_M^*·2E_R＝S|X_M|²(ΔH_LR+ΔH_RR)…(15′)

SX_S^*·2E_R＝S|X_S|²(ΔH_LR+ΔH_RR)…(16′)

从等式(13′)至(16′)分别导出：

ΔH_LL＝SX_M^*·E_L/S|X_M|²+SX_S^*·E_L/S|X_S|²…(17)

ΔH_RL＝SX_M^*·E_L/S|X_M|²-SX_S^*·E_L/S|X_S|²…(18)

ΔH_LR＝SX_M^*·E_R/S|X_M|²+SX_S^*·E_R/S|X_S|²…(19)

ΔH_RR＝SX_M^*·E_R/S|X_M|²-SX_S^*·E_R/S|X_S|²…(20)

使用由等式(17)至(20)导出的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR，在每一适当确定的规定时间周期(例如进行整体平均的时间周期)中更新滤波器装置40-1、40-2、40-3和40-4的滤波器特性。例如，假定利用与导出的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR相应的脉冲响应Δh_LL、Δh_RL、Δh_LR和Δh_RR，在第k次更新之后的脉冲响应h_LL、h_RL、h_LR和h_RR给定为h_LL(k)、h_RL(k)、h_LR(k)和h_RR(k)，则得到

h_LL(k+1)＝h_LL(k)+aΔh_LL…(21)

h_RL(k+1)＝h_RL(k)+aΔh_RL…(22)

h_LR(k+1)＝h_LR(k)+aΔh_LR…(23)

h_RR(k+1)＝h_RR(k)+aΔh_RR…(24)

这里，a是适当设置的收敛系数。

利用这些更新等式，导出第k+1个脉冲响应h_LL(k+1)、h_RL(k+1)、h_LR(k+1)和h_RR(k+1)，并分别设置到滤波器装置40-1、40-2、40-3和40-4，在每一适当确定的规定时间周期(即进行整个平均的时间周期)中，重复上述操作。

利用信号x_L和x_R或和信号和差信号x_M和x_S作为图1的立体回声消除器16、24的参考信号而实现模拟的结果，表示在各个音频传送系统的图3至8中。作为信号x_L和x_R，所用的是基于人声的立体音频信号。在图3至8中，横轴表示块数(一块表示进行整体平均的时间周期，在模拟中设置为约2.3秒)，在每一块中更新滤波器特性。第一块不设置滤波器特性，在第二块执行初始设置，然后，每块进行更新。纵轴表示传递函数的估算误差(dB)，未设置滤波器特性的初始状态定义为0dB。表1示出图3至8各个模拟的条件。

(表1)

  图号    参考信号    双向通话  传送系统变化  图3    x_L，x_R    否    否  图4    x_M，x_S    否    否  图5    x_L，x_R    是    否  图6    x_M，x_S    是    否  图7    x_L，x_R    是    是  图8    x_M，x_S    是    是

在表1中，“双向通话”表示从扬声器SP(L)和SP(R)再生的声音，以及房间中的人说话的声音，同时被传声器MC(L)和MC(R)收集的情况。因为电话会议系统一般是在发生双向通话的情况下使用的，所以要求即使在发生双向通话时，也能得到足够的回声消除性能。另一方面，在表1中，“传送系统变化”表示改变传递函数，而估算误差是收敛的，假定传声器位置移动或诸如此类。作为回声消除器的操作，要求即使估算误差由于传递函数的变化而暂时地增加，它也应再次趋于收敛。

考虑图3至8的模拟结果。比较图3和4，当信号x_L和x_R用作参考信号(图3)时，因为信号x_L和x_R之间的相关程度高，所以从操作开始，最多在第20块(约45秒)，估算误差仅下降到约-15至-25dB。比较起来，当信号x_M和x_S用作参考信号(图4)时，因为信号x_M和x_S之间的相关程度低，所以从操作开始，在第20块，误差消除下降到约-35至-45dB，因此看起来能获得足够的回声消除性能。比较表示出现双向通话的图5和6，当信号x_L和x_R用作参考信号(图5)时，从操作开始，最多在第20块，估算误差仅下降到约-10至-20dB。比较起来，当信号x_M和x_S用作参考信号(图6)时，估算误差下降到约-23至-30dB，因此可以看出，即使双向通话出现，也能获得足够的回声消除性能。这是因为，由于在房间中人说话的声音与从扬声器SP(L)和SP(R)再生的声音是不相关的，所以当传递函数计算装置58计算各个系统的传递函数时，房间中人说话的声音通过上述互频谱计算计算和整体平均计算而被消除，以致各个系统的传递函数的推导能不受双向通话影响。比较表示双向通话和传送系统的变化都出现的图7和8，当信号x_L和x_R用作参考信号(图7)时，在第11块给传送系统一个变化之后，估算误差的收敛是差的，以致估算误差在第20块最多仅下降到约-5至-15dB。比较起来，当信号x_M和x_S用作参考信号(图8)时，估算误差的收敛是很好的，甚至在第11块将变化给至传送系统之后，结果估算误差在第20块下降到约-17至-26dB，因此可以看出，即使双向通话和传送系统的变化都存在，也能得到足够的回声消除性能。

图9示出图2的立体回声消除器16、24中的另一种结构示例，其中，和/差信号产生装置配置在向扬声器的传输线路上。相同符号用于指示与图1的上述结构通用的部分。从配对侧地点传输来的并输入至线路输入端LI(L)和LI(R)的左/右两个通道立体声信号x_L和x_R被输入至和/差信号产生装置52。和/差信号产生装置52用加法器54对立体声信号x_L和x_L进行相加，以产生和信号x_M(＝x_L+x_R)，同时用减法器56对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者也可以是x_R-x_L)}。立体音频信号解调装置62，用加法器64对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器66给其结果乘以系数1/2，以恢复原始信号x_L，同时用减法器68对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器70对其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传递函数计算装置58，执行由和/差信号产生装置52产生的和信号x_M和差信号x_S与由减法器48和50输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置40-1至40-4的滤波器特性的设置和更新。它们的操作与所述图1的结构的那些操作相同。其他部件的操作也与所述图1的结构的那些相同。

图10示出图2立体回声消除器16、24的另一种结构示例，其中，在图2的地点A和B之间的传输是以和信号x_M和差信号x_S的信号形式，而不是立体声信号x_L和x_R的信号形式实现。相同的符号用于指示与图1或9的上述结构通用的部分。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_R(或者也可以是x_R-x_L)}输入至立体声信号解调装置62。立体音频信号解调装置62，用加法器64对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器66给其结果乘以系数1/2，以恢复原始信号x_L，同时用减法器68对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器70对其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_M和差信号x_S与由减法器48和50输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置40-1至40-4的滤波器特性的设置和更新。它们的操作与所述图1或9的结构的那些操作相同。和/差信号产生装置72利用加法器73，将减法器48和50输出的信号e_L和e_R相加，以产生和信号e_M(＝e_L+e_R)，同时用减法器75对它们进行相减，以产生差信号e_S{＝e_L-e_R(或者也可以是e_R-e_L)}，然后将它们发送至配对侧的地点。其它部件的操作与所述图1或9的结构的那些操作相同。

图11示出图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示图1，9或10的上述结构通用的部分。从配对侧地点传输来的并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声信号解调装置62。立体音频信号解调装置62，用加法器64对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器66给其结果乘以系数1/2，以恢复原始信号x_L，同时用减法器68对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器70对其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

和/差信号产生装置52用加法器54对由立体声信号解调装置62恢复的立体声信号进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器56对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_M和差信号x_S，与由减法器48和50输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置40-1至40-4的滤波器特性的设置和更新。它们的操作与所述图1、9或10的结构的那些操作相同。和/差信号产生装置72利用加法器73，将减法器48和50输出的信号e_L和e_R相加，以产生和信号e_M(＝e_L+e_R)，同时用减法器75对它们进行相减，以产生差信号e_S{＝e_L-e_R(或者也可以是e_R-e_L)}，然后将它们发送至配对侧的地点。其他部件的操作与所述图1、9或10的结构的那些操作相同。

图12示出立体回声消除器16、24的结构示例。从配对侧的地点传送来的并输入至线路输入端LI(L)和LI(R)的左/右两个通道立体声信号x_L和x_R，照原样从声音输出端SO(L)和SO(R)输出(即不通过和/差信号产生装置152)，并分别在扬声器SP(L)和SP(R)再生。

和/差信号产生装置152用加法器154对输入至线路输入端LI(L)和LI(R)的左/右两个通道的立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器156对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。

滤波器装置140-1至140-4由FIR滤波器构成。这些滤波器140-1至140-4各设置这样的脉冲响应，其分别与扬声器SP(L)和SP(R)以及传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数H_LL、H_LR、H_RL和H_RR中的两种合适系统传递函数的合成形式相对应，并利用这个脉冲响应进行和与差信号(低相关合成信号)的卷积计算，由此分别产生回声消除信号EC1至EC4。

加法器144进行EC1+EC3的计算。加法器146进行EC2+EC4的计算。减法器148从由声音输入端SI(L)输入的传声器MC(L)收集的音频信号中减去回声消除信号EC1+EC3，由此进行回声消除。减法器150从由声音输入端SI(R)输入的传声器MC(R)收集的音频信号中减去回声消除信号EC2+EC4，由此进行回声消除。左和右通道的这些回声消除信号，各分别从线路输出端LO(L)和LO(R)输出，向配对侧的地点发送。

传递函数计算装置158执行由和/差信号产生装置152产生的信号x_M和差信号x_L与由减法器148和150输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置140-1至140-4的滤波器特性(脉冲响应)的设置和更新。具体地说，刚启动系统时，未设置滤波器装置140-1至140-4的滤波器特性，即系数全部设置为零，致使回声消除信号EC1至EC4是零，因此，传声器MC(L)和MC(R)它们本身收集的音频信号从减法器148和150输出。所以在这时，传递函数计算装置158执行由和/差信号产生装置152产生的和信号x_M和差信号x_S，与减法器148和150输出的所收集的传声器MC(L)和MC(R)收集的音频信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，分别推导出多个合成的传递函数，每一传递函数在形式上是扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数H_LL、H_LR、H_RL和H_RR中两种适用的系统的传递函数的合成，并执行将滤波器装置140-1至140-4的滤波器特性初始设置为与这个合成传递函数相对应的值。在初始设置之后，因为回声消除信号是由滤波器装置140-1至140-4产生的，所以回声消除误差信号e_L和e_R，相应于传声器MC(L)和MC(R)收集的音频信号与回声消除EC1至EC4之间的差信号，从减法器148和150输出。因此，在这时，传递函数计算装置158执行由和/并信号产生装置152产生的和信号x_M和差信号x_S，与由减法器148和150输出的回声消除误差信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，分别推导上述合成传递函数的估算误差，并将滤波器装置140-1至140-4的滤波器特性更新为消除这个估算误差的值。通过每一规定时间周期重复这种更新操作，回声消除误差可收敛至最小值。另外，即使传递函数由于传声器位置的移动或诸如此类而变化，但通过借助于上面的连续更新与其有关的滤波器装置140-1至140-4的滤波器特性，回声消除误差也可收敛于最小值。

相关检测装置160基于相关计算或此类计算，检测和信号x_M和差信号x_S之间的相关，当相关值不小于规定的值时，停止更新上述滤波器特性。当相关值变成低于规定值时，重新启动上述滤波器特性的更新。

这里，将描述用传递函数计算装置158设置滤波器装置140-1至140-4的滤波器特性(脉冲响应)。在传递函数计算装置158中，进行下面的计算。

(在固定型操作的情况下)

信号x_L和x_R是

x_L＝(x_M+x_S)/2

x_R＝(x_M-x_S)/2

因此，传声器MC(L)和MC(R)的输出信号Y_L和Y_R变成

Y_L＝(X_M+X_S)H_LL/2+(X_M-X_S)H_RL/2

  ＝X_M{(H_LL+H_RL)/2}+X_S{(H_LL+H_RL)/2}…(25)

Y_R＝(X_M+X_S)H_LR/2+(X_M-X_S)H_RR/2

  ＝X_M{(H_LR+H_RR)/2}+X_S{(H_LR+H_RR)/2}…(26)

当合成的传递函数给定为

H_ML＝(H_LL+H_RL)/2

H_SL＝(H_LL+H_RL)/2

H_MR＝(H_LR+H_RR)/2

H_SR＝(H_LR+H_RR)/2

等式(25)和(26)分别变成

Y_L＝X_M·H_ML+X_S·H_SL…(25′)

Y_R＝X_M·H_MR+X_S·H_SR…(26′)

当等式(25′)和(26′)的两边都乘以X_M和X_S的复数共轭X_M^*和X_S^*并整体平均时，则分别得到

SX_M^*·Y_L＝SX_M^*·X_M·H_ML+SX_M^*·X_S·H_SL…(27)

SX_S^*·Y_L＝SX_S^*·X_M·H_ML+SX_S^*·X_S·H_SL…(28)

SX_M^*·Y_R＝SX_M^*·X_M·H_MR+SX_M^*·X_S·H_SR…(29)

SX_S^*·Y_R＝SX_S^*·X_M·H_MR+SX_S^*·X_S·H_SR…(30)

在等式(27)至(30)中，因为X_M和X_S彼此大致不相关，所以具有X_M^*·X_S或X_S^*·X_M的这一项在整体平均时近似变成零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(27)至(30)分别变成

SX_M^*Y_L＝S|X_M|²H_ML…(27′)

SX_S^*Y_L＝S|X_S|²H_SL…(28′)

SX_M^*Y_R＝S|X_M|²H_MR…(29′)

SX_S^*Y_R＝S|X_S|²H_SR…(30′)

从等式(27′)至(30′)，分别推导出

H_ML＝SX_M^*·Y_L/S|X_M|²…(31)

H_SL＝SX_S^*·Y_L/S|X_S|²…(32)

H_MR＝SX_M^*·Y_R/S|X_M|²…(33)

H_SR＝SX_S^*·Y_R/S|X_S|²…(34)

通过对这些推导出来的合成传递函数H_ML、H_SL、H_MR和H_SR加以逆付里叶变换而得到的脉冲响应h_ML、h_SL、和h_SR，分别是将设置到滤波器装置140-1、140-2、140-3和140-4的滤波器特性。因此，传递函数计算装置158基于输入的和信号x_M、差信号x_S和传声器MC(L)和MC(R)的输出信号y_L和y_R，通过对这些导出的合成传递函数加以逆付里叶变换，从等式(31)到(34)中推导脉冲响应h_ML、h_SL、h_MR和h_SR，并将导出的脉冲响应分别设置到滤波器装置140-1、140-2、140-3和140-4，进一步，通过在每一适当确定的规定时间周期(例如进行整体平均的时间周期)内重复这一计算，来更新脉冲响应。

(在自适应型操作的情况下)

假定，为滤波器装置140-1、140-2、140-3和140-4设置的滤波器特性给定为H^_ML、H^_SL、H^_MR和H^_SR(当用脉冲响应的术语表示时为h^_ML、h^_SL、h^_MR和h^_SR)，从减法器148和150输出的信号e_L和e_R变成

E_L＝Y_L-(X_M·H^_ML+X_S·H^_SL)

  ＝{(X_M+X_S)H_LL/2+(X_M-X_S)H_RL/2}-{(X_M·H^_ML+X_S·H^_RL)

  ＝X_M{(H_LL+H_RL)/2-H^_ML}+X_S{(H_LL-H_RL)/2-H^_SL}…(35)

E_R＝Y_R-(X_M·H^_MR+X_S·H^_SR)

  ＝{(X_M+X_S)H_LR/2+(X_M-X_S)H_RR/2}-{(X_M·H^_MR+X_S·H^_SR)

  ＝X_M{(H_LR+H_RR)/2-H^_MR}+X_S{(H_LR-H_RR)/2-H^_SR}…(36)

当合成传递函数给定为

H_ML＝(H_LL+H_RL)/2

H_SL＝(H_LL-H_RL)/2

H_MR＝(H_LR+H_RR)/2

H_SR＝(H_LR-H_RR)/2时，等式(35)和(36)分别变成

E_L＝X_M(H_ML-H^_ML)+X_S(H_SL-H^_SL)…(35′)

E_R＝X_M(H_MR-H^_MR)+X_S(H_SR-H^_SR)…(36′)

当合成传递函数的估算误差给定为

ΔH_ML＝H_ML-H^_ML

ΔH_SL＝H_SL-H^_SL

ΔH_MR＝H_MR-H^_MR

ΔH_SR＝H_SR-H^_SR时，等式(35′)和(36′)分别变成

E_L＝X_M·ΔH_ML+X_S·ΔH_SL…(35″)

E_R＝X_M·ΔH_MR+X_S·ΔH_SR…(36″)

当等式(35″)和(36″)乘以X_M和X_S的复数共轭X_M^*和X_S^*并整体平均时，分别得到

SX_M^*·E_L＝SX_M^*·X_M·ΔH_ML+SX_M^*·X_S·ΔH_SL…(37)

SX_S^*·E_L＝SX_S^*·X_M·ΔH_ML+SX_S^*·X_S·ΔH_SL…(38)

SX_M^*·E_R＝SX_M^*·X_M·ΔH_MR+SX_M^*·X_S·ΔH_SR…(39)

SX_S^*·E_R＝SX_S^*·X_M·ΔH_MR+SX_S^*·X_S·ΔH_SR…(40)

在等式(37)至(40)中，因为X_M和X_S彼此大致不相关，所以当整体平均时，具有X_M^*·X_S或X_S^*·X_M这一项近似为零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(37)至(40)分别变成

SX_M^*·E_L＝S|X_M|²ΔH_ML…(37′)

SX_S^*·E_L＝S|X_S|²ΔH_SL…(38′)

SX_M^*·E_R＝S|X_M|²ΔH_MR…(39′)

SX_S^*·E_R＝S|X_S|²ΔH_SR…(40′)

从等式(37′)至(40′)，分别推导出

ΔH_ML＝SX_M^*·E_L/S|X_M|²…(41)

ΔH_SL＝SX_S^*·E_L/S|X_S|²…(42)

ΔH_MR＝SX_M^*·E_R/S|X_M|²…(43)

ΔH_SR＝SX_S^*·E_R/S|X_S|²…(44)

用从等式(41)至(44)得到的估算误差ΔH_ML、ΔH_SL、ΔH_MR和ΔH_SR，在每一适当确定的规定时间周期(例如进行整体平均的时间周期)进行滤波器装置140-1、140-2、140-3和140-4的滤波器特性的更新。例如，假定利用与得到的估算误差ΔH_ML、ΔH_SL、ΔH_MR和ΔH_SR相应的脉冲响应Δh_ML、Δh_SL、Δh_MR和Δh_SR，在第k次更新后脉冲响应h_ML、h_SL、h_MR和h_SR给定为h_ML(k)、h_SL(k)、h_MR(k)和h_SR(k)，则有

h_ML(k+1)＝h_ML(k)+aΔh_ML…(45)

h_SL(k+1)＝h_SL(k)+aΔh_SL…(46)

h_MR(k+1)＝h_MR(k)+aΔh_MR…(47)

h_SR(k+1)＝h_SR(k)+aΔh_SR…(48)

用这些更新的等式，导出第(k+1)次脉冲响应h_ML(k+1)、h_SL(k+1)、h_MR(k+1)和h_SR(k+1)，并分别设置到滤波器装置140-1、140-2、140-3和140-4，这种操作在每适当确定的规定时间周期(例如进行整体平均的时间周期)重复进行。

图13示出图2的立体回声消除器16、24中另一种结构示例，其中和/差信号产生装置配置在去扬声器的传输线路上。有关图12的上述结构的那些共同位置上，使用相同的符号。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号输入至和/差信号产生装置152。和/差信号产生装置152利用加法器154对立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器156对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。立体声音频信号解调装置162用加法器164对和信号和差信号x_M和x_S进行相加，并进一步利用系统乘法器166给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器168对和和差信号x_M和x_S进行相减，并进一步利用系数乘法器170给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传递函数计算装置158执行由和/差信号产生装置152产生的和信号x_M和差信号x_S与由减法器148和150输出的信号e_L和e_R之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置140-1至140-4的滤波器特性的设置和更新。它们的操作与所述图12的结构的那些操作相同，其他部件的操作，也与所述图12的结构的那些操作相同。

图14示出图2的立体回声消除器16、24中的另一种结构示例，其中，在图2的地点A和B之间的传输，是以和信号x_M和差信号x_S的信号形式，而不是立体声信号x_L和x_R的信号形式实现的。相同的符号用于指示与图12或13的上述结构的通用部分。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声音频信号解调装置162。立体声音频信号解调装置162，用加法器164对和和差信号x_M和x_S进行相加，并进一步利用系数乘法器166给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器168对和和差信号x_M和x_S进行相减，并进一步利用系数乘法器170给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)再生。

传递函数计算装置158，执行从线路输入端LI(L)和LI(R)输入的和信号x_M和差信号x_S与从减不器148和150输出的信号e_L和e_R的互频谱计算，并基于这个互频谱计算，对滤波器装置140-1至140-4的滤波器特性进行设置和更新。它们的操作与所述图12或13的结构的那些操作相同。和/差信号产生装置172利用加法器173对从减法器148和150输出的信号e_L和e_R进行相加，以产生和信号e_M(＝e_L+e_R)，同时，用减法器175对它们进行相减，以产生差信号e_S{＝e_L-e_R(或者它也可以是e_R-e_L)}，然后将它们发送至配对侧的地点。其他部件的操作与所述图12和13的结构的那些操作相同。

图15示出图2立体回声消除器16、24的结构示例。相同的符号用于指示图12、13或14的上述结构的通用部分。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声音频信号解调装置162。立体声音频信号解调装置162利用加法器164对和信号和差信号x_M和x_S进行相加，并进一步利用系数乘法器166给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器168对和信号和差信号x_M和x_S进行相减，并进一步利用系数乘法器170给其结果乘以1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)再生。

和/差信号产生装置152利用加法器154对由立体声音频信号解调装置162恢复的立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，由减法器156对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。传递函数计算装置158，执行由和/差信号产生装置152产生的和信号x_M和差信号x_S与从减法器148和150输出的信号e_L和e_R之间的互频谱计算，基于这个互频谱计算，进行滤波器装置140-1至140-4的滤波器特性的设置和更新。它们的操作与所述图12、13和14的结构的那些操作相同。和/差信号产生装置172利用加法器173对从148和150输出的信号e_L和e_R进行相加，以产生和信号e_M(＝e_L+e_R)，同时，用减法器175对它们进行相减，以产生差信号e_S{＝e_L-e_R(或者它也可以是e_R-e_L)}，然后将它们发送至配对侧的地点。其他部件的操作与所述图12、13和14的结构的那些操作相同。

图16表示立体回声消除器16、24的另一种结构示例。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声音频信号解调装置262。立体声音频信号解调装置262利用加法器264对和信号和差信号x_M和x_S进行相加，并进一步利用系数乘法器266给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器268对和信号和差信号x_M和x_S进行相减，并进一步利用系数乘法器270给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)再生。

传声器MC(L)和MC(R)收集的音频信号输入至和/差信号产生装置272。和/差信号产生装置272利用加法器273对传声器收集的音频信号y_L和y_R进行相加，以产生和信号y_M，同时，用减法器275对它们进行相减，以产生差信号y_S。

滤波器装置240-1至240-4例如是由FIR滤波器形成。这些滤波器240-1至240-4各设置这样的脉冲响应，其分别与扬声器SP(L)和SP(R)和传声器MC(L)和MC(R)之间的四个音频传送系统H_LL、H_LR、H_RL和H_RR中两个合适系统的合成传递函数的合成形式相对应，并利用这个脉冲响应进行左/右两通道立体声信号x_L和x_R的卷积计算，由此分别产生回声消除信号EC1至EC4。

加法器244进行EC1+EC3的计算。加法器246进行EC2+EC4的计算。减法器248从和信号y_M减去回声消除信号EC1+EC3，从而进行回声消除。减法器250从差信号y_S减去EC2+EC4，从而进行回声消除。从减法器248和250输出的信号e_M和e_S分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。

传递函数计算装置258，执行由线路输入端LI(L)和LI(R)输出的信号x_M和差信号x_S与由减法器248和250输出的信号e_M和e_S之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置240-1至240-4的滤波器特性(脉冲响应)的设置和更新。具体地说，刚启动系统时，未设置滤波器装置240-1至240-4的滤波器特性，即系数全部设置为零，以致回声消除信号EC1至EC4为零，因此，从和/差信号产生装置272输出的信号y_M和差信号y_S，照原样从减法器248和250输出。所以，这时传递函数计算装置258，执行由线路输入端LI(L)和LI(R)输出的和信号x_M和差信号x_S与由减法器248和250输出的y_M和y_S之间的互频谱计算，基于这个互频谱计算，分别推导多个合成传递函数，每个都是扬声器SP(L)和SP(R)与MC(L)和MC(R)之间的四个音频传送系统的H_LL、H_LR、H_RL和H_RR中的两个合适系统的合成传递函数的形式，并将滤波器装置240-1至240-4的滤波器特性，初始设置为这个合成传递函数相应的值。在初始设置之后，因为回声消除信号是由滤波器装置240-1至240-4产生的，所以，和/差信号产生装置272输出的和信号y_M和差信号y_S与回声消除信号EC1至EC4之间的差信号相应的回声消除误差信号e_M和e_S，从减法器248和250输出。因此，这时传递函数计算装置258执行由线路输入端LI(L)和LI(R)输入的和信号x_M和差信号x_S与由减法器248和250输出的回声消除误差消除差信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，分别推导上述合成传递函数，并分别更新滤波器装置240-1至240-4。通过在每一规定时间周期重复这种更新的操作，回声消除误差可收敛至最小值。另外，即使传递函数由于传声器位置移动或诸如此类而变化，回声消除误差也可通过逐次地更新与其有关的滤波器装置240-1至240-4的滤波器特性，收敛至最小值。

相关检测装置260基于相关值计算或类似的计算，检测和信号x_M和差信号x_S之间的相关性，当相关值不小于规定的值时，停止更新上述的滤波器特性。当相关值变成低于规定值时，则重新启动上述滤波器特性的更新。

这里，将描述通过传递函数计算装置258设置滤波器装置240-1至240-4的滤波器特性(脉冲响应)。在传递函数计算装置258中，进行如下计算。

(在固定型操作的情况下)

假定y_M和y_S的频率轴表示分别给定为Y_M和Y_S，传声器MC(L)和MC(R)的输出信号y_L和y_R的和信号和差信号y_M和y_S变成

Y_M＝Y_L+Y_R

  ＝(X_L·H_LL+X_R·H_RL)+(X_L·H_LR+X_R·H_RR)

  ＝X_L(H_LL+H_LR)+X_R(H_RL+H_RR)…(49)

Y_S＝Y_L-Y_R

  ＝(X_L·H_LL+X_R·H_RL)-(X_L·H_LR+X_R·H_RR)

  ＝X_L(H_LL-H_LR)+X_R(H_RL-H_RR)…(50)

当合成传递函数给定为

H_LM＝H_LL+H_LR

H_RM＝H_RL+H_RR

H_LS＝H_LL-H_LR

H_RS＝H_RL-H_RR

时，则

X_L＝(X_M+X_S)/2

X_R＝(X_M-X_S)/2

因此，等式(49)和(50)分别变成

Y_M＝(X_M+X_S)H_LM/2+(X_M-X_S)H_RM/2

  ＝X_M(H_LM+H_RM)/2+X_S(H_LM-H_RM)/2

Y_S＝(X_M+X_S)H_LS/2+(X_M-X_S)H_RS/2

  ＝X_M(H_LS+H_RS)/2+X_S(H_LS-H_RS)/2

所以，

2Y_M＝X_M(H_LM+H_RM)+X_S(H_LM-H_RM)…(49′)

2Y_S＝X_M(H_LS+H_RS)+X_S(H_LS-H_RS)…(50′)

当等式(49′)和(50′)的两边都乘以X_M和X_S的复数共轭的X_M^*和X_S^*并整体平均时，分别得到

SX_M^*·Y_M＝SX_M^*·X_M(H_LM+H_RM)+SX_M^*·X_S(H_LM-H_RM)…(51)

SX_S^*·Y_M＝SX_S^*·X_M(H_LM+H_RM)+SX_S^*·X_S(H_LM-H_RM)…(52)

SX_M^*·Y_S＝SX_M^*·X_M(H_LS+H_RS)+SX_M^*·X_S(H_LS-H_RS)…(53)

SX_S^*·Y_S＝SX_S^*·X_M(H_LS+H_RS)+SX_S^*·X_S(H_LS-H_RS)…(54)

在等式(51)至(54)中，因为X_M和X_S彼此大致不相关，所以当整体平均时，具有X_M^*·X_S或X_S^*·X_M的项近似变成零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(51)至(54)分别变成

SX_M^*·Y_M＝S|X_M|²(H_LM+H_RM)…(51′)

SX_S^*·Y_M＝S|X_S|²(H_LM-H_RM)…(52′)

SX_M^*·Y_S＝S|X_M|²(H_LS+H_RS)…(53′)

SX_S^*·Y_S＝S|X_S|²(H_LS-H_RS)…(54′)

通过变换等式(51′)至(54′)，分别推导出下列合成传递函数

H_LM+H_RM＝SX_M^*·Y_M/S|X_M|²…(51″)

H_LM-H_RM＝SX_S^*·Y_M/S|X_S|²…(52″)

H_LS+H_RS＝SX_M^*·Y_S/S|X_M|²…(53″)

H_LS-H_RS＝SX_S^*·Y_S/S|X_S|²…(54″)

从等式(51″)至(54″)分别推导出

H_LM＝SX_M^*·Y_M/S|X_M|²+SX_S^*·Y_M/S|X_S|²…(55)

H_RM＝SX_M^*·Y_M/S|X_M|²-SX_S*·Y_M/S|X_S|²…(56)

H_LS＝SX_M^*·Y_S/S|X_M|²+SX_S*·Y_S/S|X_S|²…(57)

H_RS＝SX_M*·Y_S/S|X_M|²-SX_S^*·Y_S/S|X_S|²…(58)

对这些导出的合成传递函数H_LM、H_RM、H_LS和H_RS加以逆付里叶变换得到的脉冲响应h_LM、h_RM、h_LS和h_RS，将被分别设置到滤波器装置240-1、240-2、240-3和240-4的滤波器特性。传递函数计算设置258，基于输入至线路输入端LI(L)和LI(R)的和信号x_M和差信号x_S，以及从和/差信号产生装置272输出的和信号y_M和差信号y_S，由等式(55)至(58)推导出各个合成传递函数H_LM、H_RM、H_LS和H_RS，然后对这些推导出的合成传递函数加以逆付里叶变换，推导出脉冲响应h_LM、h_RM、h_LS和h_RS，将所推导出的脉冲响应分别设置到滤波器装置240-1、240-2、240-3和240-4，并进一步，通过在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复这一操作，来更新脉冲响应。

(在自适应型操作情况下)

假定e_M和e_S的频率轴表示分别给定为E_M和E_S，并且设置到滤波器装置240-1、240-2、240-3和240-4的滤波器特性给定为H^_LM、H^_RM、H^_LS和H^_RS(当用脉冲响应的术语表示时为h^_LM、h^_RM、h^_LS和h^_RS)，从减法器248和250输出的信号e_M和e_S则变成

E_M＝{L(H_LL+H_LR)+X_R(H_RL+H_RR)}-(X_L·H^_LM+X_R·H^_RM)

  ＝L{(H_LL+H_LR)-H^_LM}+X_R{(H_RL+H_RR)-H^_RM}…(59)

E_S＝{L(H_LL-H_LR)+X_R(H_RL-H_RR)}-(X_L·H^_LS+X_R·H^_RS)

   ＝L{(H_LL-H_LR)-H^_LS}+X_R{(H_RL-H_RR)-H^_RS}…(60)

当合成传递函数给定为

H_LM＝H_LL+H_LR

H_RM＝H_RL+H_RR

H_LS＝H_LL-H_LR

H_RS＝H_RL-H_RR

则

X_L＝(X_M+X_S)/2

X_R＝(X_M-X_S)/2

因此，等式(59)和(60)分别变成

E_M＝(X_M+X_S)·(H_LM-H^_LM)/2+(X_M-X_S)·(H_RM-H^_RM)/2…(61)

E_S＝(X_M+X_S)·(H_LS-H^_LS)/2+(X_M-X_S)·(H_RS-H^_RS)/2…(62)

当合成传递函数的估算误差给定为

ΔH_LM＝H_LM-H^_LM

ΔH_RM＝H_RM-H^_RM

ΔH_LS＝H_LS-H^_LS

ΔH_RS＝H_RS-H^_RS

等式(61)和(62)分别变成

E_M＝X_M(ΔH_LM+ΔH_RM)/2+X_S(ΔH_LM-ΔH_RM)/2

E_S＝X_M(ΔH_LS+ΔH_RS)/2+X_S(ΔH_LS-ΔH_RS)/2

因此

2E_M＝X_M(ΔH_LM+ΔH_RM)+X_S(ΔH_LM-ΔH_RM)…(61′)

2E_S＝X_M(ΔH_LS+ΔH_RS)+X_S(ΔH_LS-ΔH_RS)…(62′)

当等式(61′)和(62′)的两边都乘以X_M和X_S的复数共轭的X_M^*和X_S^*并整体平均时，则分别得到

SX_M^*·2E_M＝SX_M^*·X_M(ΔH_LM+ΔH_RM)+SX_M^*·X_S(ΔH_LM-ΔH_RM)…(63)

SX_S^*·2E_M＝SX_S^*·X_M(ΔH_LM+ΔH_RM)+SX_S^*·X_S(ΔH_LM-ΔH_RM)…(64)

SX_M^*·2E_S＝SX_M^*·X_M(ΔH_LS+ΔH_RS)+SX_M^*·X_S(ΔH_LS-ΔH_RS)…(65)

SX_S^*·2E_S＝SX_S^*·X_M(ΔH_LS+ΔH_RS)+SX_S^*·X_S(ΔH_LS-ΔH_RS)…(64)

在等式(63)至(66)中，因为X_M和X_S彼此大致不相关，所以当整体平均时，具有X_M^*·X_S或X_S^*·X_M的这一项近似变成零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(63)至(66)分别变成

SX_M^*·2E_M＝S|X_M|²(ΔH_LM+ΔH_RM)…(63′)

SX_S^*·2E_M＝S|X_S|²(ΔH_LM-ΔH_RM)…(64′)

SX_M^*·2E_S＝S|X_M|²(ΔH_LS+ΔH_RS)…(65′)

SX_S^*·2E_S＝S|X_S|²(ΔH_LS-ΔH_RS)…(66′)

通过变换等式(63′)至(66′)，分别推导出下列合成传递函数

ΔH_LM+ΔH_RM＝SX_M^*·2E_M/S|X_M|²…(63″)

ΔH_LM-ΔH_RM＝SX_S^*·2E_M/S|X_S|²…(64″)

ΔH_LS+ΔH_RS＝SX_M^*·2E_S/S|X_M|²…(65″)

ΔH_LS-ΔH_RS＝SX_S^*·2E_S/S|X_S|²…(64″)

从等式(63″)至(66″)，分别推导出

ΔH_LM＝SX_M^*·2E_M/S|X_M|²+SX_S^*·2E_M/S|X_S|²…(67)

ΔH_RM＝SX_M^*·2E_M/S|X_M|²-SX_S*·2E_M/S|X_S|²…(68)

ΔH_LS＝SX_M^*·2E_S/S|X_M|²+SX_S*·2E_S/S|X_S|²…(69)

ΔH_RS＝SX_M^*·2E_S/S|X_M|²-SX_S^*·2E_S/S|X_S|²…(70)

利用由等式(67)至(70)推导出的估算误差ΔH_LM、ΔH_RM、ΔH_LS和ΔH_RS，在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)更新滤波器装置240-1、240-2、240-3和240-4的滤波器特性。例如，假定用与推导出的估算误差ΔH_LM、ΔH_RM、ΔH_LS和ΔH_RS相应的脉冲响应Δh_LM、Δh_RM、Δh_LS和Δh_RS，在第k次更新之后的脉冲响应h_LM、h_RM、h_LS和h_RS给定为h_LM(k)、h_RM(k)、h_LS(k)和h_RS(k)，则有

h_LM(k+1)＝h_LM(k)+aΔh_LM…(71)

h_RM(k+1)＝h_RM(k)+aΔh_RM…(72)

h_LS(k+1)＝h_LS(k)+aΔh_LS…(73)

h_RS(k+1)＝h_RS(k)+aΔh_RS…(74).

利用这些更新的等式，分别推导出第k+1次脉冲响应h_LM(k+1)、h_RM(k+1)、h_LS和h_RS，并设置到滤波器装置240-1、240-2、240-3和240-4，这种操作在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复进行。

图17示出图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示与图16的上述结构的通用部分。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声x_L和x_R，分别照原样从声音输出端SO(L)和SO(R)输出(即不经过和/差信号产生装置252)，并在扬声器SP(L)和SP(R)中再生。和/差信号产生装置252用加法器254对这个立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器256对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。

传递函数计算装置258，执行由和/差信号产生装置252产生的信号x_M和差信号x_S与由减法器248和250输出的信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置240-1至240-4的滤波器特性的设置和更新。它们的操作与所述图16的结构的那些操作相同。从减法器248和250输出的信号e_M和e_S输入至立体声音频信号解调装置282。立体声音频信号解调装置282用加法器284对信号e_M和e_S进行相加，并进一步用系数乘法器286给其结果乘以系数1/2，以恢复左通道信号e_L，同时，用减法器288对信号e_M和e_S进行相减，并进一步用系数乘法器290给其结果乘以系数1/2，以恢复右通道信号e_R。所恢复的信号e_L和e_R分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图16的结构的那些操作相同。

图18示出图2的立体回声消除器16、24的另一种结构示例。相同的符号用于指示与图16或17的上述结构的通用部分。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声音频信号解调装置262。立体声音频信号解调装置262用加法器264对和信号和差信号x_M和x_S进行相加，并进一步利用系数乘法器266给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器268对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器270给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。和/差信号产生装置252用加法器254对这个立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器256对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。

传递函数计算装置258执行由和/差信号产生装置252产生的和信号x_M和差信号x_S与由减法器248和250输出的信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置240-1至240-4的滤波器特性的设置和更新。它们的操作与所述图16或17的结构的那些操作相同。由减法器输出的信号e_M和e_S分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图16或17的结构的那些操作相同。

图19示出图2立体回声消除器16、14的另一种结构示例。相同符号用于指示与图16、17和18的上述结构的通用部分。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R输入至和/差信号产生装置252。和/差信号产生装置252用加法器254对立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器256对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。所产生的和信号x_M和差信号x_S输入至立体声音频信号解调装置262。立体声音频信号解调装置262用加法器264对和信号和差信号x_M和x_S进行相加，进一步用系数乘法器266给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器268对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器270给其结果乘以1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传递函数计算装置258执行由和/差信号产生装置252产生的和信号x_M和差信号x_S与由减法器248和250输出的信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置240-1至240-4的滤波器特性的设置和更新。它们的操作与所述图16、17或18的结构的那些操作相同。从减法器248和250输出的信号e_M和e_S输入至立体声音频信号解调装置282。立体声音频信号解调装置282用加法器284对信号e_M和e_S进行相加，并进一步利用系数乘法器286给其结果乘以1/2，以恢复左通道信号e_L，同时，用减法器288对信号e_M和e_S进行相减，并进一步用系数乘法器290给其结果乘以系数1/2，以恢复右通道信号e_R。所恢复的信号e_L和e_R分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图16、17或18的结构的那些操作相同。

图20示出立体回声消除器16、24的另一种结构示例。从配对侧地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_S(或者也可以是x_R-x_L)}输入至立体声音频信号解调装置362。立体声音频信号解调装置362用加法器364对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器366给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器368对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器370给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传声器MC(L)和MC(R)收集的信号y_L和y_R输入至和/差信号产生装置372。和/差信号产生装置372用加法器373对传声器收集的音频信号y_L和y_R进行相加，以产生和信号y_M，同时，用减法器375对它们进行相减，以产生差信号y_S。

滤波器装置340-1至340-4由FIR滤波器形成。这些滤波器装置340-1至340-4各设置这样的脉冲响应，这些脉冲响应分别与形式为扬声器SP(L)和SP(R)和传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数H_LL、H_LR、H_RL和H_RR合成的合成传递函数相对应，并用这些脉冲响应分别进行由线路输入端LI(L)和LI(R)输入的和信号x_M和差信号x_S的卷积计算，由此产生回声消除信号EC1至EC4。

加法器344进行EC1+EC3的计算。加法器346进行EC2+EC4的计算。减法器348从和信号y_M中减去回声消除信号EC1+EC3，从而实现回声消除。减法器350从差信号y_S中减去回声消除信号EC2+EC4，从而实现回声消除。由减法器348和350输出的信号e_M和e_S分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。

传递函数计算装置358执行由线路输入端LI(L)和LI(R)输出的和信号x_M和差信号x_S与由减法器348和350输出的信号e_M和e_S之间的互频谱计算，并基于这个互频谱计算，进行滤波器装置340-1至340-4的滤波器特性(脉冲响应)的设置和更新。具体地说，刚启动系统时，未设置滤波器装置340-1至340-4的滤波器特性，即系数全部设置为零，以致回声消除信号EC1至EC4是零，因此，从和/差信号产生装置372输出的和信号y_M和差信号y_S照原样从减法器348和350输出。所以，这时传递函数计算装置358执行由线路输入端LI(L)和LI(R)输入的和信号x_M和差信号x_S与由减法器348和350输出的和信号e_M和差信号e_S之间的互频谱计算，基于这个互频谱计算，分别推导出多个合成传递函数，每个在形式上是扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数H_LL、H_LR、H_RL和H_RR的合成，并将滤波器装置340-1至340-4的滤波器特性，初始设置为与这个合成传递函数相对应的值。在初始设置之后，因为回声消除信号是由滤波器装置340-1至340-4产生的，所以，与从和/差信号产生装置372输出的和信号y_M和差信号y_S与回声消除信号EC1至EC4之间的差信号相对应的回声消除误差信号e_M和e_S，从减法器348和350输出。因此，这时传递函数计算装置358执行从线路输入端LI(L)和LI(R)输入的和信号x_M和差信号x_S与从减法器348和350输出的回声消除误差信号e_M和e_S之间的互频谱计算，并基于这个互频谱计算，分别推导出上述合成传递函数的估算误差，并分别将滤波器装置340-1至340-4的滤波器特性更新到消除这个估算误差的值。通过每一规定时间周期重复这一更新操作，回声消除误差可收敛至最小值。另外，即使传递函数由于传声器位置的移动或诸如此类而变化，但回声消除误差，通过连续更新与其有关的滤波器装置340-1至340-4的滤波器特性，也可收敛至最小值。

相关检测装置360基于相关值计算或诸如此类的计算，检测和信号x_M与差信号x_S之间的相关性，当相关值不小于规定的值时，停止更新上述滤波器特性。当相关值变成低于规定值时，则重新启动上述滤波器特性的更新。

这里，将描述由传递函数计算装置358为滤波器装置340-1至340-4设置的滤波器特性(脉冲响应)。在传递函数计算装置358中，进行如下的计算。

(在固定型操作的情况下)

传声器MC(L)和MC(R)的输出信号y_L和y_R的和信号和差信号y_M和y_S变成

Y_M＝Y_L+Y_R

  ＝(X_L·H_LL+X_R·H_RL)+(X_L·H_LR+X_R·H_RR)

  ＝X_L(H_LL+H_LR)+X_R(H_RL+H_RR)…(71)

Y_S＝Y_L-Y_R

  ＝(X_L·H_LL+X_R·H_RL)-(X_L·H_LR+X_R·H_RR)

  ＝X_L(H_LL-H_LR)+X_R(H_RL-H_RR)…(72)

X_L＝(X_M+X_S)/2

X_R＝(X_M-X_S)/2

因此，等式(71)和(72)分别变成

Y_M＝(X_M+X_S)·(H_LL+H_LR)/2+(X_M-X_S)·(H_RL+H_RR)/2

  ＝X_M(H_LL+H_LR+H_RL+H_RR)/2+X_S(H_LL+H_LR-H_RL-H_RR)/2…(71′)

Y_S＝(X_M+X_S)·(H_LL+H_LR)/2+(X_M-X_S)·(H_RL-H_RR)/2

  ＝X_M(H_LL-H_LR+H_RL-H_RR)/2+X_S(H_LL-H_LR-H_RL+H_RR)/2…(72′)

当合成传递函数给定为

H_MM＝(H_LL+H_LR+H_RL+H_RR)/2

H_SM＝(H_LL+H_LR-H_RL-H_RR)/2

H_MS＝(H_LL-H_LR+H_RL-H_RR)/2

H_SS＝(H_LL-H_LR-H_RL+H_RR)/2时，等式(71′)和(72′)分别变成

Y_M＝X_M·H_MM+X_S·H_SM…(71″)

Y_S＝X_M·H_MS+X_S·H_SS…(72″)

当等式(71″)和(72″)两边同时乘以x_M和x_S的复数共轭的X_M^*和X_S^*并整体平均时，则分别得到

SX_M^*·Y_M＝SX_M^*·X_M·H_MM+SX_M^*·X_S·H_SM…(73)

SX_S^*·Y_M＝SX_S^*·X_M·H_MM+SX_S^*·X_S·H_SM…(74)

SX_M^*·Y_S＝SX_M^*·X_M·HM_S+SX_M^*·X_S·H_SS…(75)

SX_S^*·Y_S＝SX_S^*·X_M·H_MS+SX_S^*·X_S·H_SS…(76)

在等式(73)至(76)中，因为X_M和X_S彼此大致不相关，所以，当整体平均时具有X_M^*·X_S或X_S^*·X_M这一项近似变成零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(73)至(76)分别变成

SX_M^*·Y_M＝S|X_M|²·H_MM…(73′)

SX_S^*·Y_M＝S|X_S|²·H_SM…(74′)

SX_M^*·Y_S＝S|X_M|²·H_MS…(75′)

SX_S^*·Y_S＝S|X_S|²·H_SS…(74′)

从等式(73′)至(76′)，分别推导出

H_MM＝SX_M^*·Y_M/S|X_M |²…(77)

H_SM＝SX_S^*·Y_M/S|X_S|²…(78)

H_MS＝SX_M^*·Y_S/S|X_M|²…(79)

H_SS＝SX_S^*·Y_S/S|X_S|²…(80)

对这些推导出的合成传递函数H_MM、H_SM、H_MS和H_SS加以逆付里叶变换得到的脉冲响应h_MM、h_SM、h_MS和h_SS，分别是为滤波器装置340-1、340-2、340-3和340-4设置的滤波器特性。因此，传递函数计算装置358基于输入至线路输入端LI(L)和LI(R)的和信号x_M和差信号x_S以及从和/差信号产生装置372输出的和信号y_M和差信号y_S，由等式(77)至(80)推导各自的合成传递函数H_MM、H_SM、H_MS和H_SS对这些推导出的合成传递函数加以付里叶变换，推导脉冲响应h_MM、h_SM、h_MS和h_SS；将推导出的脉冲响应分别设置到滤波器装置340-1、340-2、340-3和340-4；进一步，通过在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复这一计算，更新脉冲响应。

(在自适应型操作的情况下)

假定滤波器装置340-1、340-2、340-3和340-4设置的滤波器特性给定为H^_MM、H^_SM、H^_MS和H^_SS(当以脉冲响应术语表示时为h^_MM、h^_SM、h^_MS和h^_SS)，从减法器348和350输出的信号e_M和e_S则变成

E_M＝{(X_M+X_S)·(H_LL+H_LR)/2+(X_M-X_S)·(H_RL+H_RR)/2}

    -(X_M·H^_MM+X_S·H^_SM)

  ＝X_M[{H_LL+H_LR+H_RL+H_RR)/2}-H^_MM]

    +X_S[{(H_LL+H_LR-H_RL-H_RR)/2}-H^_SM]…(81)

E_S＝{(X_M+X_S)·(H_LL-H_LR)/2+(X_M-X_S)·(H_RL-H_RR)/2}

    -(X_M·H^_MS+X_S·H^_SS)

＝X_M[{(H_LL-H_LR+H_RL-H_RR)/2}-H^_MS]

    +X_S[{(H_LL-H_LR-H_RL+H_RR)/2}-H^_SS]…(82)

当合成传递函数给定为

H_MM＝(H_LL+H_LR+H_RL+H_RR)/2

H_SM＝(H_LL+H_LR-H_RL-H_RR)/2

H_MS＝(H_LL-H_LR+H_RL-H_RR)/2

H_SS＝(H_LL-H_LR-H_RL+H_RR)/2时，等式(81)和(82)分别变成

E_M＝X_M(H_MM-H^_MM)+X_S(H_SM-H^_SM)…(81′)

E_S＝X_M(H_MS-H^_MS)+X_S(H_SS-H^_SS)…(82′)

当给定

ΔH_MM＝H_MM-H^_MM

ΔH_SM＝H_SM-H^_SM

ΔH_MS＝H_MS-H^_MS

ΔH_SS＝H_SS-H^_SS时，等式(81′)和(82′)分别变成

E_M＝X_M·ΔH_MM+X_S·ΔH_SM…(81″)

E_S＝X_M·ΔH_MS+X_S·ΔH_SS…(82″)

当等式(81″)和(82″)两边同时乘以X_M和X_S的复数共轭X_M^*和X_S^*并整体平均时，分别得到

SX_M^*·E_M＝SX_M^*·X_M·ΔH_MM+SX_M^*·X_S·ΔH_SM…(83)

SX_S^*·E_M＝SX_S^*·X_M·ΔH_MM+SX_S^*·X_S·ΔH_SM…(84)

SX_M*·E_S＝SX_M^*·X_M·ΔH_MS+SX_M^*·X_S·ΔH_SS…(85)

SX_S^*·E_S＝SX_S^*·X_M·ΔH_MS+SX_S^*·X_S·ΔH_SS…(86)

在等式(83)至(86)中，因为X_M和X_S彼此大致不相关，所以当整体平均时，具有X_M^*·X_S或X_S^*·H_M这一项近似变成零。另外，

X_M^*·X_M＝|X_M|²

X_S^*·X_S＝|X_S|²

因此，等式(83)至(86)分别变成

SX_M^*·E_M＝S|X_M|²·ΔH_MM…(83′)

SX_S^*·E_M＝S|X_S|²·ΔH_SM…(84′)

SX_M^*·E_S＝S|X_M|²·ΔH_MS…(85′)

SX_S^*·E_S＝S|X_S|²·ΔH_SS…(86′)

从等式(83′)至(86′)，分别推导出

ΔH_MM＝SX_M^*·E_M/S|X_M|²…(87)

ΔH_SM＝SX_S^*·E_M/S|X_S|²…(88)

ΔH_MS＝SX_M^*·E_S/S|X_M|²…(89)

ΔH_SS＝SX_S*·E_S/S|X_S|²…(90)

用从等式(87)至(90)推导出的估算误差ΔH_MM、ΔH_SM、ΔH_MS和ΔH_SS，在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)更新滤波器装置340-1、340-2、340-3和340-4的滤波器特性。例如，假定利用与所推导的估算误差ΔH_MM、ΔH_SM、ΔH_MS和ΔH_SS相对应的脉冲响应Δh^_MM、Δh^_SM、Δh^_MS和Δh^_SS，在第k次更新之后的脉冲响应h_MM、h_SM、h_MS和h_SS给定为h_MM(k)、h_SM(k)、h_MS(k)和h_SS(k)，则得到

h_MM(k+1)＝h_MM(k)+aΔH_MM…(91)

h_SM(k+1)＝h_SM(k)+aΔh_SM…(92)

h_MS(k+1)＝h_MS(k)+aΔh_MS…(93)

h_SS(k+1)＝h_SS(k)+aΔh_SS…(94)

利用这些更新的等式，推导出第(k+1)次脉冲响应h_MM(k+1)、h_SM(k+1)、h_MS(k+1)和h_SS(k+1)，并分别设置到滤波器装置340-1、340-2、340-3和340-4，在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复进行这一操作。

图21表示图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示图20上述结构的通用部分。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R，照原样从声音输出端SO(L)和SO(R)输出(即不经过和/差信号产生装置352)，并分别在扬声器SP(L)和SP(R)再生。和/差信号产生装置352用加法器354对这些立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时用减法器356对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。

传递函数计算装置358执行由和/差信号产生装置352产生的和信号x_M和差信号x_S与由减法器348和350输出的信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置340-1至340-4的滤波器特性的设置和更新。它们的操作与所述图20的结构的那些操作相同。从减法器348和350输出的信号e_M和e_S输入至立体声音频信号解调装置382。立体声音频信号解调装置382用加法器384对信号e_M和e_S进行相加，并进一步用系数乘法器386给其结果乘以系数1/2，以恢复左通道信号e_L，同时，用减法器388对信号e_M和e_S进行相减，并进一步利用系数乘法器390给其结果乘以系数1/2，以恢复右通道信号e_R。所恢复的信号e_L和e_R分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图20的结构的那些操作相同。

图22示出图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示图20或21的上述结构的通用部分。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的和信号x_M(＝x_L+x_R)和差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}，输入至立体声音频信号解调装置362。立体声音频信号解调装置362用加法器364对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器366给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器368对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器370给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。和/差信号产生装置352用加法器354对这样的立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时用减法器356对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。

传递函数计算装置358执行由和/差信号产生装置352产生的和信号x_M和差信号x_S与由减法器348和350输出的e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置340-1至340-4的滤波器特性的设置和更新。它们的操作与所述图20或21的结构的那些操作相同。从减法器348和350输出的信号e_M和e，分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图20或21的结构的那些操作相同。

图23示出图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示图20、21或22上述结构的通用部件。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R输入至和/差信号产生装置352。和/差信号产生装置352用加法器354对立体声信号x_L和x_R进行相加，以产生和信号x_M(＝x_L+x_R)，同时，用减法器356对它们进行相减，以产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。所产生的和信号x_M和差信号x_S输入至立体声音频信号解调装置362。立体声音频信号解调装置362用加法器364对和信号和差信号x_M和x_S进行相加，并进一步用系数乘法器366给其结果乘以系数1/2，以恢复原始信号x_L，同时，用减法器368对和信号和差信号x_M和x_S进行相减，并进一步用系数乘法器370给其结果乘以系数1/2，以恢复原始信号x_R。所恢复的信号x_L和x_R分别从声音输出端SO(L)和SO(R)输出，并在扬声器SP(L)和SP(R)中再生。

传递函数计算装置358执行由和/差信号产生装置352产生的和信号x_M和差信号x_S与由减法器348和350输出的信号e_M和e_S之间的互频谱计算，基于这个互频谱计算，进行滤波器装置340-1至340-4的滤波器特性的设置和更新。它们的操作与所述图20、21或22的结构的那些操作相同。从减法器348和350输出的信号e_M和e_S输入至立体声音频信号解调装置382。立体声音频信号解调装置382用加法器364对信号e_M和e_S进行相加，并进一步用系数乘法器386给其结果乘以系数1/2，以恢复左通道信号e_L，同时，用减法器388对信号e_M和e_S进行相减，并进一步利用系数乘法器390给其结果乘以系数1/2，以恢复右通道信号e_R。所恢复的信号e_L和e_R，分别从线路输出端LO(L)和LO(R)输出，并向配对侧的地点传输。其他部件的操作与所述图20、21或22的结构的那些操作相同。

本发明除上述实例所示的结构以外，还可以有各种不同的结构。例如，图1、9至11的结构的部件可改变为像图24那样。具体地说，和/差信号产生装置402安排在加法器44和46与减法器48和50之间，以便用加法器404对回声消除信号EC1+EC3和EC2+EC4进行相加，以产生(EC1+EC3)+(EC2+EC4)，同时，用减法器406对它们进行相减，以产生(EC1+EC3)-(EC2+EC4)。另外，和/差信号产生装置408安排在声音输入端SI(L)和SI(R)与减法器48和50之间，以便用加法器410对传声器MC(L)和MC(R)收集的音频信号y_L和y_R进行相加，以产生和信号y_M，同时，用减法器412对它们进行相减，以产生差信号y_S。减法器48从和信号y_M中减去回声消除信号(EC1+EC3)+(EC2+EC4)，以实现回声消除。减法器50从差信号y_S中减去回声消除信号(EC1+EC3)-(EC2+EC4)，以实现回声消除。从减法器48和50输出的信号e_M和e_S输入至立体声音频信号解调装置414。立体声音频信号解调装置414用加法器416对信号e_M和e_S进行相加，并进一步用系数乘法器418给其结果乘以系数1/2，以恢复左通道信号e_L，同时，用减法器420对信号e_M和e_S进行相减，并进一步用系数乘法器给其结果乘以系数1/2，以恢复右通道信号e_R。其他部件的操作与所述图1、9至10的那些操作相同。

另一方面，图12至15的结构部件可改变为像图25那样。具体地说，和/差信号产生装置432安排在加法器144和146与减法器148和150之间，以便用加法器434对回声消除信号EC1+EC3和EC2+EC4进行相加，以产生回声消除信号(EC1+EC3)+(EC2+EC4)，同时，用减法器436对它们进行相减，以产生(EC1+EC3)-(EC2+EC4)。另外，和/差信号产生装置438安排在声音输入端SI(L)和SI(R)与减法器148和150之间，以便用加法器440对传声器MC(L)和MC(R)收集的音频信号y_L和y_R进行相加，以产生和信号y_M，同时，用减法器442对它们进行相减，以产生差信号y_S。减法器148从和信号y_M中减去回声消除信号(EC1+EC3)+(EC2+EC4)以实现回声消除。减法器150从差信号y_S中减去回声消除信号(EC1+EC3)-(EC2+EC4)，以实现回声消除。从减法器148和150输出的信号e_M和e_S输入至立体声音频信号解调装置444。立体声音频信号解调装置444用加法器446对信号e_M和e_S进行相加，并进一步用系数乘法器448给其结果乘以系数1/2，以恢复左通道信号e_L，同时，用减法器450对信号e_M和e_S进行相减，并进一步用系数乘法器452给其结果乘以系数1/2，以恢复右通道信号e_R。其他部件的操作与所述图12至15的结构的那些操作相同。

另一方面，图16至19的部分结构可改变为像图26那样。具体地说，和/差信号产生装置462安排在加法器244和246与减法器248和250之间，以便用加法器464对回声消除信号EC1+EC3和EC2+EC4进行相加，以产生(EC1+EC3)+(EC2+EC4)，同时，用减法器466对它们进行相减，以产生(EC1+EC3)-(EC2+EC4)。图16至19的和/差信号产生装置272不需要。减法器248从传声器MC(L)收集的音频信号y_L中减去回声消除信号(EC1+EC3)+(EC2+EC4)以实现回声消除。减法器250从传声器MC(R)收集的音频信号y_R中减去回声消除信号(EC1+EC3)-(EC2+EC4)，以实现回声消除。和/差信号产生装置468安排在减法器248和250的输出侧，以便用加法器470对减法器248和250的输出信号e_L和e_R进行相加，以产生和信号e_M，同时，用减法器472对它们进行相减，以产生差信号e_S。其他部件的操作与所述图16至19的那些操作相同。

另一方面，图20至21的部分结构可改变为像图27那样。具体地说，和/差信号产生装置482安排在加法器344和346与减法器348和350之间，以便用加法器484对回声消除信号EC1+EC3和EC2+EC4进行相加，以产生回声消除信号(EC1+EC3)+(EC2+EC4)，同时，用减法器486对它们进行相减，以产生(EC1+EC3)-(EC2+EC4)。图20至23的和/差信号产生装置372不需要。减法器348从传声器MC(L)收集的音频信号y_L中减去回声消除信号(EC1+EC3)+(EC2+EC4)以实现回声消除。减法器350从传声器MC(R)收集的音频信号y_R中减去回声消除信号(EC1+EC3)-(EC2+EC4)，以实现回声消除。和/差信号产生装置488安排在减法器348和350的输出侧，以便用加法器490对减法器348和350的输出信号e_L和e_R进行相加，以产生和信号e_M，同时，用减法器492对它们进行相减，以产生差信号e_S。其他部件的操作与所述图20至23的那些操作相同。

图28示出图2立体回声消除器16、24的另一种结构示例。相同的符号用于指示与图1的上述结构的通用部分。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R，照原样从声音输出端SO(L)和SO(R)输出(不经过正交化滤波器500)，并分别在扬声器SP(L)和SP(R)再生。

滤波器装置40-1设置有与扬声器SP(L)和传声器MC(L)之间的传递函数相应的脉冲响应，滤波器装置40-2设置有与扬声器SP(L)和传声器MC(R)之间的传递函数相应的脉冲响应，滤波器装置40-3设置有与扬声器SP(R)和传声器MC(L)之间的传递函数相应的脉冲响应，滤波器装置40-4设置有与扬声器SP(R)和传声器MC(R)之间的传递函数相应的脉冲响应。

正交化滤波器500在每一规定的时间周期，对立体声信号x_L和x_R进行主要分量的分析，并将这个输入的立体声x_L和x_R转换为彼此正交的两个信号。传递函数计算装置502执行由正交化滤波器500产生的互相正交的两个信号与从减法器48和50输出的信号e_L和e_R之间的互频谱计算，基于这个互频谱计算，设置滤波器装置40-1至40-4的滤波器特性(脉冲响应)。具体地说，刚启动系统时，未设置滤波器装置40-1至40-4的滤波器特性，即系数全部设置为零，以致回声消除信号EC1至EC4为零，因此，传声器MC(L)和MC(R)收集的音频信号它们本身从减法器48和50输出。所以，这时传递函数计算装置502，执行正交化滤波器500产生的互相正交的两个信号与由减法器48和50输出的传声器MC(L)和MC(R)收集的音频信号e_L和e_R之间的互频谱计算，基于这个互频谱计算，分别推导扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统的传递函数，并将滤波器装置40-1至40-4的滤波器特性，初始设置为与这个传递函数相对应的值。在初始设置之后，因为回声消除信号是由滤波器装置40-1至40-4产生的，所以，传声器MC(L)和MC(R)收集的音频信号与回声消除信号EC1至EC4之间的差信号所对应的回声消除误差信号e_L和e_R，从减法器48和50输出。所以，这时传递函数计算装置502，执行由正交化滤波器500产生的回声消除误差信号e_L和e_R之间的互频谱计算，基于这个互频谱计算，分别推导扬声器SP(L)和SP(R)与传声器MC(L)和MC(R)之间的四个音频传送系统传递函数的估算误差，并分别将滤波器装置40-1至40-4的滤波器特性，更新为消除估算误差的值。通过每一规定时间周期重复这一更新操作，回声消除误差可收敛至最小值。另外，即使传递函数由于传声器位置移动或诸如此类而变化，回声消除误差也可通过逐次地更新与其有关的滤波器装置40-1至40-4的滤波器特性，而收敛至最小值。

根据基于输入至线路输入端LI(L)和LI(R)的左/右双通道立体声信号x_L和x_R与由线路输出端LO(L)和LO(R)输出的信号e_L和e_R之间比较的已知技术，通话重叠检测装置504检测那些除扬声器SP(L)和SP(R)再生的声音以外的声音输入至传声器MC(L)和MC(R)的通话重叠。传递函数计算装置502在检测到通话重叠时，使滤波器装置40-1至40-4的滤波器特性的更新周期相对较长，而在未检测到通话重叠时，使滤波器特性的更新周期相对较短。这在通话重叠存在时，能充分收敛估算误差，并在通话重叠不存在时，加快估算误差的收敛。

下面，将描述正交化滤波器500的正交化处理。正交化处理在输入立体声信号的每一规定时间周期中执行。这里，正交化处理就图29所示的每一帧(例如512个样本)来进行(但是，重叠处理可按图49所示的方式进行，将在后面描述)。给出两个矢量如[等式1]所示，即具有输入至正交化滤波器500的左通道输入信号x_L的一帧样本组作为其元素的矢量，和具有右通道输入信号x_R的一帧样本组作为其元素的矢量：

[等式1]

$>>>>x>\to >>L>>=>>(>>x>>L>1>>>,>>x>>L>2>>>,>>x>>L>3>>>,>·>·>·>,>>x>Ln>>)>>>s>$

$>>>>x>\to >>R>>=>>(>>x>>R>1>>>,>>x>>R>2>>>,>>x>>R>3>>>,>·>·>·>,>>x>Rn>>)>>>s>$

因为x_L和x_R是立体声信号，所以它们互相之间有相关性。正交处理的实现，是每帧用x_L和x_R作为两个变量，对由这两个变量合成而成的样本组进行主要分量分析，以推导互相正交的第一主要分量和第二主要分量的本征矢量，并将这两个变量合成而成的各个样本，分别投影至所推导的第一主要分量和第二主要分量的本征矢量。

下面将描述正交处理计算的具体内容，现在，假定观察矩阵B给定为

[等式2]

$>>B>=>>$ >\begin{array}{}>>>>x>>L>1>>>,>>x>>L>2>>>,>>x>>L>3>>>,>·>·>·>,>>x>Ln>>>>>>>x>>R>1>>>,>>x>>R>2>>>,>>x>>R>3>>>,>·>·>·>,>>x>Rn>>>>>>s>\end{array}$$

则B的协方差矩阵S变成

[等式3]

$>>S>=>>1>>n>->1>>>>BB>T>>>s>(BT是B的转置矩阵)$

$>>=>>1>>n>->1>>>>$ >>\begin{array}{}>>>>x>>L>1>>>,>>x>>L>2>>>,>>x>>L>3>>>,>·>·>·>,>>x>Ln>>>>>>>x>>R>1>>>,>>x>>R>2>>>,>>x>>R>3>>>,>·>·>·>,>>x>Rn>>>>>>\begin{array}{}>>>>x>>L>1>>>>>>x>>R>1>>>>>>>>x>>L>2>>>>>>x>>R>2>>>>>>>·>>>·>>>>>·>>>·>>>>>>x>Ln>>>>>x>Rn>>>>>>>s>\end{array}\end{array}$$

$>>=>>1>>n>->1>>>\begin{array}{}>>>>\Sigma >>k>=>1>>n>>>>x>Lk>>2>>>>>\Sigma >>k>=>1>>n>>>x>Lk>>>x>Rk>>>>>>>\Sigma >>k>=>1>>n>>>x>Lk>>>x>Rk>>>>>\Sigma >>k>=>1>>n>>>>x>Lk>>2>>>>>>>s>\end{array}$

$>>=>\begin{array}{}>>>>S>11>>>>>S>12>>>>>>>S>21>>>>>S>22>>>>>>>s>\end{array}$

(S₁₁是x_L的方差，S₂₂是x_R的方差，和S₁₂(＝S₂₁)是x_L和x_R的协方差)

因此，从[等式4]

$>>\begin{array}{}>>>>S>11>>->\lambda >>>>S>12>>>>>>>S>21>>>>>S>22>>->\lambda >>>>>=>0>>s>\end{array}$

可得：

[等式5]

(S₁₁-λ)(S₂₂-λ)-S₁₂S₂₁＝0

因此[等式6]

$>>\lambda >=>>>>S>11>>+>>S>22>>\pm >>>>(>>S>11>>+>>S>22>>)>>2>>->4>>(>>S>11>>>S>22>>->>>S>12>>2>>)>\; >>2>>>s>$

导出本征值λ的两个解。

假定具有较大方差的两个本征值之一(第一主要分量的本征值)是λ₁，与本征值λ₁相应的本征矢量U_max则是

[等式8]

$>>\begin{array}{}>>>>u>11>>>>>>>u>12>>>>>>,>>s>\end{array}$

其建立

[等式7]

$>>\begin{array}{}>>>>S>11>>>>>S>12>>>>>>>S>21>>>>>S>22>>>>>>\begin{array}{}>>>>u>11>>>>>>>u>12>>>>>>=>>\lambda >1>>\begin{array}{}>>>>u>11>>>>>>>u>12>>>>>>>s>\end{array}\end{array}\end{array}$

这里u₁₁²+u₁₂²＝1

通过解u₁₁和u₁₂，推导出

[等式9]

$>>>u>11>>=>\pm >>>S>12>>>>>S>12>>2>>->>>(>>\lambda >1>>->>S>11>>)>>2>\; >>>s>$

$>>>u>12>>=>\pm >>>>\lambda >1>>->>S>11>>>>>>S>12>>2>>->>>(>>\lambda >1>>->>S>11>>)>>2>\; >>>s>$

(同一阶双重符号)

不管u₁₁和u₁₂为正或负的符号，由第一主要分量表示的轴是相同的。

另一方面，假定具有较小方差两个本征值之一(第二主要分量)是λ₂，与本征值λ₂相应的本征矢量U_min是

[等式11]

$>>\begin{array}{}>>>>u>21>>>>>>>u>22>>>>>>,>>s>\end{array}$

其建立

[等式10]

$>>\begin{array}{}>>>>S>11>>>>>S>12>>>>>>>S>21>>>>>S>22>>>>>>\begin{array}{}>>>>u>21>>>>>>>u>22>>>>>>=>>\lambda >2>>\begin{array}{}>>>>u>21>>>>>>>u>22>>>>>>>s>\end{array}\end{array}\end{array}$

这里u₂₁²+u₂₂²＝1

通过解u₂₁和u₂₂，推导出

[等式12]

$>>>u>21>>=>\pm >>>S>12>>>>>S>12>>2>>->>>(>>\lambda >2>>->>S>11>>)>>2>\; >>>s>$

$>>>u>22>>=>\pm >>>>\lambda >2>>->>S>11>>>>>>S>12>>2>>->>>(>>\lambda >2>>->>S>11>>)>>2>\; >>>s>$

(同一阶双重符号)

不管u₂₁和u₂₂为正或负的符号，由第二主要分量表示的轴是相同的。

根据上述，第一主要分量和第二主要分量推导如下：

[等式13]

$>>>>U>\to >>max>>=>\begin{array}{}>>>>u>11>>>>>>>u>12>>>>>>,>>>U>\to >>min>>=>\begin{array}{}>>>>u>21>>>>>>>u>22>>>>>>>s>\end{array}\end{array}$

至于从协方差矩阵得到的本征矢量U_max和U_min，不能预测本征矢量本来该出现在哪个象限。如果每一帧本征矢量出现的象限有变化，则在以的每块对U_max′作整块平均时，U_max′它们本身互相抵消。如果本征矢量出现象限有变化，则以后在每块对U_min作整体平均时，U_min′它们本身互相抵消。因此，为了使本征矢量U_max和U_min出现的象限固定，需进行转换操作。例如，当将本征矢量U_max固定至第一象限，并将本征矢量U_min固定至第四象限时，可通过如下转换操作实现。

·本征矢量U_max和U_min之一存在于第一象限(正，正)或第三象限(负，负)，给定为U_min′。和

·本征矢量U_max和U_min之一存在于第二象限(负，正)或第四象限(正，负)，给定为U_min′。然后，转换按这样进行：

·当U_max′存在于第一象限时，U_max＝U_max′

·当U_max′存在于第三象限时，U_max＝-U_max′

·当U_min′存在于第二象限时，U_min＝-U_min′

·当U_min′存在于第四象限时，U_min＝U_min′。

通过这样的转换操作，能固定本征矢量U_max和U_min出现的象限。由[等式14]

$>>>b>\to >>=>\begin{array}{}>>>>x>Ln>>>>>>>x>Rn>>>>>>>s>\end{array}$

给出的观察矩阵B的列矢量，投影到上述所得的第一主要分量的本征矢量U_max和第二主要分量的本征矢量U_min上。观察矢量B投影到本征矢量U_max上得到的输出信号x_max的值，推导为

[等式15]

$>>>x>max>>=>>b>\to >>·>>>U>\to >>max>>>s>(·表示内积)$

另外，观察矩阵B投影到本征矢量U_min上得到的输出信号x_min的值，推导为

[等式16]

$>>>x>min>>=>>b>\to >>·>>>U>\to >>min>>>s>(·表示内积)$

下面，将描述传递函数计算装置502的传递函数计算处理。(在固定型操作情况下)

正交化滤波器500的滤波器特性给定为U。滤波器特性U是将输入信号x_L和x_R投影到基于主要分量分析的互不相关信号的x_max和x_min上的特性，可建立下列关系：

[等式17]

滤波器特性U的逆特性给定为V。逆滤波器特性V是将信号x_max和x_min恢复到原始信号x_L和x_R的特性，可建立下列关系：

[等式18]

传声器MC(L)和MC(R)的输出信号y_L和y_R表示为

Y_L＝H_LL·X_L+H_RL·X_R…(97)

Y_R＝H_LR·X_L+H_RR·X_R…(98)

从等式(95)和(96)，有

[等式19]

$>>\begin{array}{}>>>>x>L>>>>>>>x>R>>>>>>=>\begin{array}{}>>>>v>11>>>>>v>12>>>>>>>v>21>>>>>v>22>>>>>>\begin{array}{}>>>>x>max>>>>>>>x>min>>>>>>·>·>·>>(>99>)>>>s>\end{array}\end{array}\end{array}$

因此，假定x_max和x_min的频率轴表示分别给定为X_max和X_min，则等式(97)和(98)分别变成

Y_L＝H_LL(v₁₁·X_max+v₁₂·X_min)+H_RL(v₂₁·X_max+v₂₂·X_min)…(97′)

Y_R＝H_LR(v₁₁·X_max+v₁₂·X_min)+H_RR(v₂₁·X_max+v₂₂·X_min)…(98′)

当等式(97′)的两边同时乘以X_max和X_min的复数共轭的X^*_max和X^*_min(即推导互谱)并整体平均时，由于X_max和X_min是互相正交的，所以X^*_max·X_min和X^*_min·X_max的期望值分别变成零，因而得到下列两个等式：

E[X^*_max·Y_L]＝E[X^*_max·H_LL·v₁₁·X_max

              +X^*_max·H_RL·v₂₁·X_max]…(100)

E[X^*_min·Y_L]＝E[X^*_min·H_LL·v₁₂·X_min

              +X^*_min·H_RL·v₂₂·X_min]…(101)

(注意：E[]表示整体平均)。

类似地，当等式(98′)的两边同时乘以X_max和X_min的复数共轭的X^*_max和X^*_min并整体平均时，则得到下列的两个等式：

E[X^*_max·Y_R]＝E[X^*_max·H_LR·v₁₁·X_max

              +X^*_max·H_RR·v₂₁·X_max]…(102)

E[X^*_min·Y_R]＝E[X^*_min·H_LR·v₁₂·X_min

              +X^*_min·H_RR·v₂₂·X_min…(103)

这里，如果在进行整体平均的时间周期中本征值λ的变化小，则可建立

[等式20]

因而将等式(100)至(103)重写为下列等式(100′)至(103′)：

E[X^*_max·Y_L]～E[|X_max|²]·H_LL·E[v₁₁]

              +E[|X_max|²]·H_RL·E[v₂₁]…(100′)

E[X^*_min·Y_L]～E[|X_min|²]·H_LL·E[v₁₂]

              +E[|X_min|²]·H_RL·E[v₂₂]…(101′)

E[X^*_max·Y_R]～E[|X_max|²]·H_LR·E[v₁₁]

              +E[|X_max|²]·H_RR·E[v₂₁]…(102′)

E[X^*_min·Y_R]～E[|X_min|²]·H_LR·E[v₁₂]

              +E[|X_min|²]·H_RR·E[v₂₂]…(103′)

当分别将等式(100′)和(102′)两边同时除以E[|X_max|²]，以及等式(101′)和(103′)同时除以E[|X_min|²]时，则得

E[X^*_max·Y_L]/E[|X_max|²]～H_LL·E[v₁₁]+H_RL·E[v₂₁]

E[X^*_min·Y_L]/E[|X_min|²]～H_LL·E[v₁₂]+H_RL·E[v₂₂]

E[X^*_max·Y_R]/E[|X_max|²]～H_LR·E[v₁₁]+H_RR·E[v₂₁]

E[X^*_min·Y_R]/E[|X_min|²]～H_LR·E[v₁₂]+H_RR·E[v₂₂]

因此，得到

[等式21]

$>>\begin{array}{}>>>>>E>[sup>>X>max>*sup>>·>>Y>L>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>Y>L>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>>>E>[sup>>X>max>*sup>>·>>Y>R>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>Y>R>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>\approx >\begin{array}{}>>>>H>LL>>>>>H>RL>>>>>>>H>LR>>>>>H>RR>>>>>>\begin{array}{}>>>E>[>>v>11>>]>>>E>[>>v>12>>]>>>>>E>[>>v>21>>]>>>E>[>>v>22>>]>>>>>·>·>·>>(>105>)>>>s>\end{array}\end{array}\end{array}$

从E[U]·E[V]～I，等式(105)重写为

[等式22]

$>>\begin{array}{}>>>>H>LL>>>>>H>RL>>>>>>>H>LR>>>>>H>RR>>>>>>\approx >\begin{array}{}>>>>>E>[sup>>X>max>*sup>>·>>Y>L>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>Y>L>>]>>>E>[>>X>min>>>|>2>>]>>>>>>>>>E>[sup>>X>max>*sup>>·>>Y>R>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>Y>R>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>\begin{array}{}>>>E>[>>u>11>>]>>>E>[>>u>12>>]>>>>>E>[>>u>21>>]>>>E>[>>u>22>>]>>>>>·>·>·>>(>106>)>>>s>\end{array}\end{array}\end{array}$

对由等式(106)推导出的传递函数H_LL、H_LR、H_RL和H_RR加以逆付里叶变换得到的脉冲响应h_LL、h_RL、h_LR和h_RR，分别是要设置到滤波器装置40-1、40-2、40-3和40-4的滤波器特性。因此，传递函数计算装置502，基于从正交化滤波器500输出的信号x_max和x_min、正交化滤波器500的滤波器特性U和传声器MC(L)和MC(R)的输出信号y_L和y_R，推导各个传递函数H_LL、H_RL、H_LR和H_RR；通过对这些推导出的传递函数加以逆付里叶变换，推导脉冲响应h_LL、h_RL、h_LR和h_RR；分别将推导出的脉冲响应设置到滤波器装置40-1、40-2、40-3和40-4；并进一步，通过在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复这一计算，更新脉冲响应。

(在自适应型操作的情况下)

假定设置到滤波器装置40-1、40-2、40-3和40-4的滤波器特性给定为H^_LL、H^_RL、H^_LR和H^_RR(当用脉冲响应的术语表示时为h^_LL、h^_RL、h^_LR和h^_RR)，则从图28的减法器48和50输出的信号e_L和e_R表示为

E_L＝(H_LL·X_L+H_RL·X_R)-(H^_LL·X_L+H^_RL·X_R)…(107)

E_R＝(H_LR·X_L+H_RR·X_R)-(H^_LR·X_L+H^_RR·X_R)…(108)

从上述等式(99)，等式(107)和(108)分别变为

E_L＝(H_LL-H^_LL)(v₁₁·X_max+v₁₂·X_min)

    +(H_RL-H^_RL)(v₂₁·X_max+v₂₂·X_min)…(107′)

E_R＝(H_LR-H^_LR)(v₁₁·X_max+v₁₂·X_min)

    +(H_RR-H^_RR)(v₂₁·X_max+v₂₂·X_min)…(108′)

当传递函数的估算误差给定为

ΔH_LL＝H_LL-H^_LL

ΔH_RL＝H_RL-H^_RL

ΔH_LR＝H_LR-H^_LR

ΔH_RR＝H_RR-H^_RR时，等式(107′)和(108′)分别变为

E_L＝ΔH_LL(v₁₁·X_max+v₁₂·X_min)

    +ΔH_RL(v₂₁·X_max+v₂₂·X_min)…(107″)

E_R＝ΔH_LR(v₁₁·X_max+v₁₂·X_min)

    +ΔH_RR(v₂₁·X_max+v₂₂·X_min)…(108″)

当等式(107″)的两边同时乘以X_max和X_min的复数共轭的X^*_max和X^*_min(即推导互谱)并整体平均时，由于X_max和X_min互相正交，所以X^*_max·X_min和X^*_min·X_max的期望值分别变为零，因而得到下列两个等式

(注意：E[]表示整体平均)：

E[X^*_max·E_L]＝E[X^*_max·ΔH_LL·v₁₁·X_max

               +X^*_max·ΔH_RL·v₂₁·X_max]…(109)

E[X^*_min·E_L]＝E[X^*_min·ΔH_LL·v₁₂·X_min

               +X^*_min·ΔH_RL·v₂₂·X_min]…(110)

类似地，当等式(108″)两边同时乘以X_max和X_min的复数共轭的X^*_max和X^*_min并整体平均时，得到下列两个等式：

E[X^*_max·E_R]＝E[X^*_max·ΔH_LR·v₁₁·X_max

               +X^*_max·ΔH_RR·v₂₁·X_max]…(111)

E[X^*_min·E_R]＝E[X^*_min·ΔH_LR·v₁₂·X_min

               +X^*_min·ΔH_RR·v₂₂·X_min]…(112)

这里，如果在进行整体平均时间周期中，本征值λ的变化小，则可建立上述等式(104)，因而等式(109)至(112)重写为下列等式(109′)至(112′)：

E[X^*_max·E_L]～E[|X_max|²]·ΔH_LL·E[v₁₁]

               +E[|X_max|²]·ΔH_RL·E[v₂₁]…(109′)

E[X^*_min·E_L]～E[|X_min|²]·ΔH_LL·E[v₁₂]

               +E[|X_min|²]·ΔH_RL·E[v₂₂]…(110′)

E[X^*_max·E_R]～E[|X_max|²]·ΔH_LR·E[v₁₁]

               +E[|X_max|²]·ΔH_RR·E[v₂₁]…(111′)

E[X^*_min·E_R]～E[|X_min|²]·ΔH_LR·E[v₁₂]

               +E[|X_min|²·ΔH_RR·E[v₂₂]…(112′)

当等式(109′)和(111′)的两边同时除以E[|X_max|²]和等式(110′)和(112′)两边同时除以E[|X_min|²]时，则有

E[X^*_max·E_L]/E[|X_max|²]～ΔH_LL·E[v₁₁]+ΔH_RL·E[v₂₁]

E[X^*_min·E_L]/E[|X_min|²]～ΔH_LL·E[v₁₂]+ΔH_RL·E[v₂₂]

E[X^*_max·E_R]/E[|X_max|²]～ΔH_LR·E[v₁₁]+ΔH_RR·E[v₂₁]

E[X^*_min·E_R]/E[|X_min|²]～ΔH_LR·E[v₁₂]+ΔH_RR·E[v₂₂]

因此，得到[等式23]

$>>\begin{array}{}>>>>>E>[sup>>X>max>*sup>>·>>E>L>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>E>L>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>>>E>[sup>>X>max>*sup>>·>>E>R>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>E>R>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>\approx >\begin{array}{}>>>\Delta >>H>LL>>>>>\Delta H>RL>>>>>>\Delta >>H>LR>>>>\Delta >>H>RR>>>>>>\begin{array}{}>>>E>[>>v>11>>]>>>E>[>>v>12>>]>>>>>E>[>>v>21>>]>>>E>[>>v>22>>]>>>>>·>·>·>>(>113>)>>>s>\end{array}\end{array}\end{array}$

从E[U]·E[V]～I，等式(113)重写为[等式24]

$>>\begin{array}{}>>>>\Delta H>LL>>>>\Delta >>H>RL>>>>>>\Delta >>H>LR>>>>\Delta >>H>RR>>>>>>\approx >\begin{array}{}>>>>>E>[sup>>X>max>*sup>>·>>E>L>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>E>L>>]>>>E>[>>X>min>>>|>2>>]>>>>>>>>>E>[sup>>X>max>*sup>>·>>E>R>>]>>>E>[>|>>X>max>>>|>2>>]>>>>>>>E>[sup>>X>min>*sup>>·>>E>R>>]>>>E>[>|>>X>min>>>|>2>>]>>>>>>>\begin{array}{}>>>E>[>>u>11>>]>>>E>[>>u>12>>]>>>>>E>[>>u>21>>]>>>E>[>>u>22>>]>>>>>·>·>·>>(>114>)>>>s>\end{array}\end{array}\end{array}$

利用从等式(114)推导的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR，在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)更新滤波器装置40-1、40-2、40-3和40-4的滤波器特性。例如，假定在第k次更新之后的脉冲响应h_LL、h_RL、h_LR和h_RR给定为h_LL(k)、h_RL(k)、h_LR(k)和h_RR(k)，利用与推导的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR相应的脉冲响应Δh_LL、Δh_RL、Δh_LR和Δh_RR，则有

h_LL(k+1)＝h_LL(k)+aΔh_LL

h_RL(k+1)＝h_RL(k)+aΔh_RL

h_LR(k+1)＝h_LR(k)+aΔh_LR

h_RR(k+1)＝h_RR(k)+aΔh_RR

利用这些更新的等式，推导出第(k+1)次脉冲响应h_LL(k+1)、h_RL(k+1)、h_LR(k+1)和h_RR(k+1)，并分别设置到滤波器装置40-1、40-2、40-3和40-4，这种操作在每一适当确定的规定时间周期(例如，进行整体平均的时间周期)重复。

图30示出图28的正交化滤波器500的功能方块图。立体声信号x_L和x_R分别从输入端506和508输入。协方差矩阵计算装置510推导每一帧的输入立体声信号x_L和x_R的协方差矩阵S。基于推导的协方差矩阵S，本征矢量计算装置512推导每一帧的第一和第二主要分量的本征矢量U_max和U_min的系数u₁₁、u₁₂、u₂₁和u₂₂。推导出的系数u₁₁、u₁₂、u₂₁和u₂₂分别设置到系数乘法器514、516、518和520。系数乘法器514和516对输入信号x_L分别乘以系数u₁₁和u₂₁，得到x_L·u₁₁和x_L·u₂₁。系数乘法器518和520对输入信号x_R分别乘以系数u₁₂和u₂₂，得到x_R·u₁₂和x_R·u₂₂。加法器522计算x_L·u₁₁+X_R·u₁₂，作为将输入立体声信x_L和x_R投影在本征矢量U_max上所得到的信号x_max。加法器524计算x_L·u₂₁+X_R·u₂₂，作为将立体声信号x_L和x_R投影在本征矢量U_min上所得的信号X_min。信号x_max和x_min分别从输出端526和528输出，同时本征矢量的系数u₁₁、u₂₁、u₁₂和u₂₂从输出端530输出。

图31示出图28传递函数计算装置502的功能方块图。信号x_max和x_min分别从输入端532和534输入，同时本征矢量系数u₁₁、u₂₁、u₁₂和u₂₂从输入端536输入。FFT装置538对信号x_max加以快速付里叶变换。复数共轭的计算装置540计算X_max复数共轭的X^*_max。功率谱计算装置542计算X_max·X^*_max＝|X_max|²。整体平均装置544计算E[|X_max|²]。FFT装置546对信号x_min加以快速付里叶变换。复数共轭的计算装置548计算X_min复数共轭的X^*_min。功率谱计算装置550计算X_min·X*_min＝|X_min|²。整体平均装置552计算E[|X_min|²]。

输入端554输入传声器MC(L)的输出信号y_L(或减法器48的输出信号E_L)。FFT装置558对信号y_L(或信号EL)加以快速付里叶变换。互频谱计算装置560推导X^*_max·Y_L(或X^*_max·E_L)，整体平均装置562推导E[X^*_max·Y_L](或E[X^*_max·E_L])。互频谱计算装置564推导X^*_min·Y_L(或X^*_min·E_L)，整体平均装置566推导E[X^*_min·Y_L](或E[X^*_min·E_L])。

输入端568输入传声器MC(R)的输出信号y_R(或减法器50的输出信号E_R)。FFT装置570对信号y_R(或信号E_R)加以快速付里叶变换。互频谱计算装置572推导X^*_max·Y_R(或X^*_max·E_R)，整体平均装置574推导E[X^*_max·Y_R](或E[X^*_max·E_R])。互频谱计算装置576推导X^*_min·Y_R(或X^*_min·E_R)，整体平均装置578推导E[X^*_min·Y_R](或E[X^*_min·E_R])。

合成传递函数计算装置580，基于上述推导出的E[|X_max|²]、E[|X_min|²]、E[X^*_max·Y_L](或E[X^*_max·E_L])、E[X^*_min·Y_L](或E[X^*_min·E_L])、E[X^*_max·Y_R](或E[X^*_max·E_R])和E[X^*_min·Y_R](或E[X^*_min·E_R])，推导等式(106){或等式(114)}右边的第一项。平均装置582各个地平均本征矢量的系数u₁₁、u₂₁、u₁₂和u₂₂，以推导E[v₁₁]、E[v₁₂]、E[v₂₁]和E[v₂₂]。

基于合成传递函数计算装置580和平均装置582的输出，传递函数计算装置584对等式(106){或等式(114)}右边进行计算。以推导各个传递函数H_LL、H_RL、H_LR和H_RR(或它们的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR)。逆FFT装置586对推导的传递函数H_LL、H_RL、H_LR和H_RR(或它们的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR)加以逆快速付里叶变换，以推导相应的脉冲响应h_LL、h_RL、h_LR和h_RR(或它们的估算误差ΔH_LL、ΔH_RL、ΔH_LR和ΔH_RR)，并分别将它们从输出端588、590、592和594输出。如图32所示，传递函数计算装置584和逆FFT装置586它们之间的位置可以交换。

对于如此构成的图28的立体回声消除器16、25，进行自适应型操作模拟的结果，表示在与一个音频传送系统相联系的图33至48中。图33至40表示回声消除量的时域变化，而图41至48表示传递函数估算误差的时域变化。这里，模拟是在如下条件下实施的。

·取样频率：11.025kHz

·一帧中的样本数：4096个样本

·一块中的帧数：可变(2帧、4帧、8帧、16帧)

·滤波器特性的更新周期：每一块(在一块中帧数为两帧的情况下约0.75秒，在4帧情况下约为1.5秒，在8帧的情况下约为3秒，在16帧情况下约为6秒)

·一块中整体平均的平均次数：一块中的帧数为两帧的情况下为31次，在4帧的情况下为63次，在8帧的情况下为127次，在16帧的情况下为255次{为了增加平均的次数，如图49所示，一帧划分成16个间隔，在连续地重叠数据提取时，与一帧相应的数据从每一划分间隔的开头提取，由此推导有待整体平均的各个参数值，并在一块中对推导出的各个参数值分别进行整体平均。所以，整体平均的平均次数N变成这样一个值，即N＝(16×一块中的帧数-1)。}

在任何情况下，滤波器特性不设置在第一块中，初始设置在第二块中执行，此后，每块进行更新。纵坐标(dB)轴的定义以0dB表示未设置滤波器特性的初始状态。涉及图33至48的模拟条件的差别如下：

[表2]

  图号  通话重叠存在/不存在  一块中的帧数  图33、图41  不存在  2  图34、图42  4  图35、图43  8  图36、图44  16  图37、图45  存在  2  图38、图46  4  图39、图47  8  图40、图48  16

下面讨论图33至48的模拟结果：

(1)在通话重叠不存在的情况下

因为回声消除量的上升速度是图33＞图34＞图35＞图36，所以回声消除量的上升速度随着一块(更新周期)变短(每块帧数变小)而加快。当一块的帧数为2、4或8时，所得到回声消除量约为25dB(图33、34或35)。当一块的帧数增加即为16帧时，回声消除量的减少需要一个长的时间(图36)。因为估算误差收敛速度是图41＞图42＞图43＞图44，所以估算误差的收敛速度随着块的变短而加快。

(2)在通话重叠存在的情况下

当一块的帧数是2或4帧时，回声消除量不增加(图37或38)，并且估算误差不收敛(图45或46)，因此不能估算。当一块的帧数是8时，所得到回声消除的数量约为15dB(图39)，并且估算误差收敛至约-6dB(图47)。当一块的帧数为16时，所得到回声消除量约为17dB(图40)，并且估算误差收敛至约10dB，因此能够获得十分稳定的估算。

从上述模拟结果可以看出如下几点：

(a)当未检测到通话重叠时，估算误差的收敛可通过相对地缩短滤波器特性的更新周期而加快。

(b)当检测到通话重叠时，估算误差可通过相对地延长滤波器特性的更新周期而充分地收敛。

因此，如上所述，当检测到通话重叠时，图28的传递函数计算装置502使滤波器装置40-1至40-4的滤波器特性的更新周期相对地延长，相反，当未检测到通话重叠时，使滤波器特性更新周期相对地缩短。这能使当通话重叠存在时充分收敛估算误差，而当未存在通话重叠时，加快估算误差的收敛。

图50示出图28立体回声消除器16、24的修改例。相同的符号用于指示与图28的通用部分。在这个修改例中，正交滤波器安排在扬声器SP(L)和SP(R)的信号线路上。逆滤波器596具有正交化滤波器500的逆特性{前面的等式(96)的逆滤波器特性V}，由此将正交化滤波器500的输出信号x_max和x_min恢复为原始信号x_L和x_R，并将它们馈送至扬声器SP(L)和SP(R)。

在上述实施例中，扬声器的数目是两个，传声器的数目是两个。但是，也可以构成传声器的数目是一个，扬声器的数目是两个。图51示出作为用这种方式修改图1的结果的结构示例。相同的符号用于指示与图1的通用部分。从配对侧的地点传输来的，并输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R，照原样从声音输出端SO(L)和SO(R)输出(即不经过和/差信号产生装置52)，并分别在扬声器SP(L)和SP(R)再生。

滤波器装置40-1设置有与扬声器SP(L)和传声器MC之间的传递函数H_L相应的脉冲响应，并用这个脉冲响应，进行从声音输出端SO(L)输出的信号x_L的卷积计算，从而产生与信号y_L相应的回声消除信号EC1，信号y_L是由声音输出端SO(L)输出的信号x_L在扬声器SP(L)再生，由传声器MC收集，并输入至声音输入端SI所得到的信号。滤波器装置40-3设置为具有与扬声器SP(R)和传声器MC之间的传递函数H_R相应的脉冲响应，并用这个脉冲响应，进行从声音输出端SO(R)输出的信号x_R的卷积计算，从而产生与信号y_R相应的回声消除信号EC3，信号y_R是由声音输出端SO(R)输出的信号x_R在扬声器SP(R)再生，由传声器MC收集，并输入至声音输入端SI所得到的信号。加法器44进行EC1+EC3的计算。减法器48从由输入声音输入端SI的传声器MC收集的音频信号y(＝y_L+y_R)中减去回声消除信号EC1+EC3，从而实现回声消除。已消除回声的信号e(＝e_L+e_R)从线路输出端LO输出，并向配对侧的地点传输。

和/差信号产生装置52用加法器54对输入至线路输入端LI(L)和LI(R)的左/右两通道立体声信号x_L和x_R进行相加，以便产生和信号x_M(＝x_L+x_R)，同时用减法器56对它们进行相减，以便产生差信号x_S{＝x_L-x_R(或者它也可以是x_R-x_L)}。传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_M和差信号x_S与从减法器48输出的信号e之间的互频谱计算，并基于这个互频谱计算，设置滤波器装置40-1至40-3的滤波器特性(脉冲响应)。具体地说，刚启动系统时，未设置滤波器装置40-1和40-3的滤波器特性，即系数全部设置为零，以致回声消除信号EC1和EC3为零，因此，传声器MC收集的音频信号本身从减法器48输出。所以，在这时，传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_M和差信号x_S与从减法器输出的传声器MC收集的音频信号e之间的互频谱计算，并基于这个互频谱计算，分别推导扬声器SP(L)和SP(R)与传声器MC之间的两个音频传送系统的传递函数，并将滤波器装置40-1和40-3的滤波器特性初始设置为与这些传递函数相应的值。在初始设置之后，因为回声消除信号是由滤波器装置40-1和40-3产生的，所以从减法器48输出与传声器MC收集的音频信号和回声消除信号EC1+EC3之间的差信号相应的回声消除误差信号e。因此，这时传递函数计算装置58执行由和/差信号产生装置52产生的和信号x_M和差信号x_S与从减法器48输出的回声消除误差信号e之间的互频谱计算，并基于这个互频谱计算，推导扬声器SP(L)和SP(R)与传声器MC之间的两个音频传送系统的传递函数的估算误差，并分别将滤波器装置40-1和40-3的滤波器特性更新至消除估算误差的值。通过每一规定时间周期重复这种更新操作，回声消除误差可收敛至最小值。另外，即使传递函数由于传声器位置移动或诸如此类而变化，通过连续地更新与其有关的滤波器装置40-1和40-3的滤波器特性，回声消除误差也可收敛至最小值。

相关检测装置60基于相关值计算或诸如此类，检测和信号x_M与差信号x_S之间的相关性，当相关值不小于规定值时，停止更新上述滤波器特性。当相关值变成低于规定值时，重新启动上述滤波器特性的更新。也是在除图1以外的一些实施例中，可构成扬声器数目为两个，而传声器数目为一个。

表示在上述各实施例中的立体回声消除器16、24可由专用硬件形成，或者也可经过通用计算机中的软件处理来实现。例如，图1等所示的各个方块的功能，可由组成计算机的CPU(中央处理单元)和存储装置如RAM或ROM来完成。也就是说，CPU根据存储在存储装置如ROM或RAM中的程序，可以起回声消除器的作用。

在上述实施例中，所描述的是关于对两通道立体声信号进行处理的情况。但是，对于那些彼此有相关性的三通道或更多通道的立体声信号，也能应用本发明的技术进行回声消除。