CN106416168B - 一种信号处理方法及装置 - Google Patents

一种信号处理方法及装置 Download PDF

Info

Publication number
CN106416168B
CN106416168B CN201480073837.7A CN201480073837A CN106416168B CN 106416168 B CN106416168 B CN 106416168B CN 201480073837 A CN201480073837 A CN 201480073837A CN 106416168 B CN106416168 B CN 106416168B
Authority
CN
China
Prior art keywords
reference signal
frequency domain
formula
signal
processor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201480073837.7A
Other languages
English (en)
Other versions
CN106416168A (zh
Inventor
陈特彦
刘晟
程宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Publication of CN106416168A publication Critical patent/CN106416168A/zh
Application granted granted Critical
Publication of CN106416168B publication Critical patent/CN106416168B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00—Arrangements affording multiple use of the transmission path
    • H04L5/14—Two-way operation using the same type of signal, i.e. duplex
    • H04L5/1461—Suppression of signals in the return path, i.e. bidirectional control circuits
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38—Demodulator circuits; Receiver circuits

Abstract

本发明公开了一种信号处理方法及装置，涉及通信领域，能在通信系统存在IQ不平衡时实现自干扰抵消。具体方案为：信号处理装置获取数字基带参考信号、自干扰参考信号和频域基带信号；参照数字基带参考信号获得互为镜像的基本参考信号和镜像参考信号；参照基本参考信号、镜像参考信号和自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；参照第一综合响应的估算值、第二综合响应的估算值、基本参考信号和镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。本发明用于信号处理。

Description

一种信号处理方法及装置
技术领域
本发明涉及通信领域，尤其涉及一种信号处理方法及装置。
背景技术
在移动蜂窝通信系统、无线局域网(Wireless Local Area Network，WLAN)、固定无线接入(Fixed Wireless Access，FWA)等无线通信系统中，基站(Base Station，BS)或接入点(Access Point，AP)、中继站(Relay Station，RS)以及用户设备(User Equipment，UE)等通信节点通常具有发射自身信号和接收其它通信节点信号的能力。由于无线信号在无线信道中的衰减非常大，与自身的发射信号相比，来自通信对端的信号到达接收端时信号已非常微弱，例如，移动蜂窝通信系统中一个通信节点的收发信号功率差达到80dB-140dB甚至更大，因此，为了避免同一信号发射机和信号接收机的发射信号对接收信号的自干扰，无线信号的发送和接收通常采用不同的频段或时间段加以区分。例如，在频分双工(Frequency Division Duplex，FDD)中，发送和接收使用相隔一定保护频带的不同频段进行通信，在时分双工(Time Division Duplex，FDD)中，发送和接收则使用相隔一定保护时间间隔的不同时间段进行通信，其中，FDD系统中的保护频带和FDD系统中的保护时间间隔都是为了保证接收和发送之间充分地隔离，避免发送对接收造成干扰。
无线全双工技术不同于现有的FDD或TDD技术，可以在相同无线信道上同时进行接收与发送操作，这样理论上无线全双工技术的频谱效率是FDD或TDD技术的两倍。实现无线全双工的前提在于尽可能地避免、降低与消除同一信号发射机和信号接收机的发射信号对接收信号的强干扰(称为自干扰，Self-interference)，使之不对有用信号的正确接收造成影响。
现有无线全双工系统中，发射通道的DAC(Digital to Analog Converter数字模拟转换器)、上变频及功率放大，以及接收通道的低噪声放大器(Low Noise Amplifier，LNA)、下变频及ADC(Analog to Digital Converter，模拟数字转换器)等是现有信号发射机和信号接收机的中射频单元的功能模块。对发射信号造成的自干扰抵消是通过空间干扰抑制、射频前端模拟干扰抵消、数字干扰抵消等单元来完成的。在现有无线全双工系统中，模拟自干扰抵消主要是针对经过自干扰主径的干扰信号进行消除。数字干扰抵消是模拟干扰抵消的补充，对数字接收信号中残余的自干扰信号在基带进行消除。
现有的数字干扰抵消技术主要基于自干扰信号模型进行干扰重构和抵消。然而，由于受到信号发射机和信号接收机的非理想特性的影响，接收到的自干扰信号并不能完全通过理想信号模型表示进而重构和抵消。因此，现有技术的干扰抵消性能会受限于信号发射机和信号接收机的非理想特性。以存在IQ(in-phase quadrature，同相正交)不平衡的情况为例，在一般的通信系统中，由IQ不平衡而引起的镜像干扰信号功率会比主信号功率低20～30dB，因此在不消除IQ不平衡的情况下也能正常通信。然而，对于现有的无线全双工系统而言，由于接收到的自干扰信号功率要远大于通信对端发送过来的有用信号功率，在经过射频干扰抵消后，其功率差依然大于30dB，因此数字干扰抵消模块必须完成大于40dB的自干扰抵消，才能确保与通信对端的正常通信。然而现有数字干扰抵消技术并不能有效地抵消由IQ不平衡性所引起的镜像干扰信号，使其无法实现大于40dB的自干扰抵消，也就不能实现在无线全双工系统中的自干扰抵消。
现有数字干扰抵消技术无法在通信系统存在IQ不平衡时有效的实现自干扰抵消，数字干扰抵消性能较低。
发明内容
本发明的实施例提供一种信号处理方法及装置，该信号处理方法及装置，能在通信系统存在IQ不平衡时实现自干扰抵消，提高数字干扰抵消性能。
为达到上述目的，本发明的实施例采用如下技术方案：
第一方面，提供一种信号处理装置，包括：
获取单元，用于获取数字基带参考信号、自干扰参考信号和频域基带信号；
转化单元，用于参照所述获取单元获取的所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
估算单元，用于参照所述转化单元获得的所述基本参考信号、所述转化单元获得的所述镜像参考信号和所述获取单元获取的所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；
重构单元，用于参照所述估算单元估算的所述第一综合响应的估算值、所述估算单元估算的所述第二综合响应的估算值、所述转化单元获得的所述基本参考信号和所述转化单元获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述获取单元获取的所述频域基带信号进行自干扰抵消。
结合第一方面，在第一种可能的实现方式中，所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述估算单元具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述转化单元获得的所述基本频域参考信号，S*(-f)为所述转化单元获得的所述镜像频域参考信号，Z(f)为所述获取单元获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值；
结合第一方面，在第二种可能的实现方式中，所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述估算单元具体用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述转化单元获得的所述基本时域参考信号，s*(-t)为所述转化单元获得的所述镜像时域参考信号，Z(t)为所述获取单元获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；Z(f)为Z(t)的傅里叶变换值。
结合第一方面的第一种可能的实现方式，在第三种可能的实现方式中，所述估算单元包括：
第一估算子单元，用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述转化单元获取的第i帧的S(f)、S*(-f)和所述获取单元获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)；
第二估算子单元，用于将所述转化单元从第i帧起获取的N帧的S(f)、S*(-f)和所述获取单元从第i帧起获取的N帧的Z(f)依次代入所述第一估算子单元获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
第三估算子单元，用于依据所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
第四估算子单元，用于结合所述第三估算子单元计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第一方面的第二种可能的实现方式，在第四种可能的实现方式中，所述估算单元包括：
第一估算子单元，用于依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述转化单元获取的第i帧的S(f)、S*(-f)和所述获取单元获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)；
第二估算子单元，用于将所述转化单元从第i帧起获取的N帧的S(f)、S*(-f)和所述获取单元从第i帧起获取的N帧的Z(f)依次代入所述第一估算子单元获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
第三估算子单元，用于依据所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
第四估算子单元，用于结合所述第三估算子单元计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第一方面的第三种或第四种可能的实现方式，在第五种可能的实现方式中，所述第三估算子单元具体用于将所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元，具体用于结合所述第三估算子单元计算出的所述第四公式与计算出和
结合第一方面的第三种或第四种可能的实现方式，在第六种可能的实现方式中，所述第三估算子单元具体用于将所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元，具体用于结合所述第三估算子单元计算出的所述第四代替式与计算出和
结合第一方面或第一种至第六种可能的实现方式中的任一实现方式，在第七种可能的实现方式中，所述重构单元具体用于：依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为所述自干扰信号，为所述估算单元估算的所述第一综合相应的估算值，为所述估算单元估算的所述第二综合响应的估算值，S(f)为所述转化单元获得的所述基本频域参考信号，S*(-f)为所述转化单元获得的所述镜像频域参考信号。
第二方面，提供一种信号处理装置，包括：接收器、处理器、存储器及总线，其中，所述接收器、所述处理器及所述存储器通过所述总线相互连接，
其中，所述接收器，用于获取数字基带参考信号、自干扰参考信号和频域基带信号；
所述处理器，用于参照所述接收器获取的所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
所述处理器还用于参照所述处理器获得的所述基本参考信号、所述处理器获得的所述镜像参考信号和所述接收器获取的所述自干扰参考信号估算第一综合响应的估算值和第二综合响应的估算值；
所述处理器还用于参照所述处理器估算的所述第一综合响应的估算值、所述处理器估算的所述第二综合响应的估算值、所述处理器获得的所述基本参考信号和所述处理器获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述接收器获取的所述频域基带信号进行自干扰抵消。
结合第二方面，在第一种可能的实现方式中，所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算出H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述处理器获得的所述基本频域参考信号，S*(-f)为所述处理器获得的所述镜像频域参考信号，Z(f)为所述接收器获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值；
结合第二方面，在第二种可能的实现方式中，所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述处理器具体用于依据第一代替式估算出H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述处理器获得的所述基本时域参考信号，s*(-t)为所述处理器获得的所述镜像时域参考信号，Z(t)为所述接收器获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；Z(f)为Z(t)的傅里叶变换值。
结合第二方面的第一种可能的实现方式，在第三种可能的实现方式中，所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值包括：
所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述处理器获取的第i帧的S(f)、S*(-f)和所述接收器获取的第i帧Z(f)之间的关系表达式，S′(f)＝S*(-f)；
所述处理器还用于将所述处理器从第i帧起获取的N帧的S(f)、S*(-f)和所述接收器从第i帧起获取的N帧的Z(f)代入所述处理器获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述处理器还用于依据所述处理器计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述处理器还用于结合所述处理器计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第二方面的第二种可能的实现方式，在第四种可能的实现方式中，所述处理器具体用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值包括：
所述处理器具体用于依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述处理器获取的第i帧的S(f)、S*(-f)和所述接收器获取的第i帧Z(f)之间的关系表达式，S′(f)＝S*(-f)；
所述处理器还用于将所述处理器从第i帧起获取的N帧的S(f)、S*(-f)和所述接收器从第i帧起获取的N帧的Z(f)代入所述处理器获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述处理器还用于依据所述处理器计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述处理器还用于结合所述处理器计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第二方面的第三种或第四种可能的实现方式，在第五种可能的实现方式中，所述处理器具体用于将所述处理器计算出的矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述处理器还具体用于结合所述处理器计算出的第四公式与计算出和
结合第二方面的第三种或第四种可能的实现方式，在第六种可能的实现方式中，所述处理器具体用于将所述处理器计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述处理器还具体用于结合所述处理器计算出的所述第四代替式与计算出和
结合第二方面或第一种至第六种可能的实现方式中的任一实现方式，在第七种可能的实现方式中，所述处理器还具体用于：依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为频域的所述自干扰信号，为所述处理器估算的频域的所述第一综合相应的估算值，为所述处理器估算的频域的所述第二综合响应的估算值，S(f)为所述处理器获得的所述基本频域参考信号，S*(-f)为所述处理器获得的所述镜像频域参考信号。
第三方面，提供一种信号处理方法，包括：
信号处理装置获取数字基带参考信号、自干扰参考信号和频域基带信号；
所述信号处理装置参照所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
所述信号处理装置参照所述基本参考信号、所述镜像参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；
所述信号处理装置参照所述第一综合响应的估算值、所述第二综合响应的估算值、所述基本参考信号和所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
结合第三方面，在第一种可能的实现方式中，所述信号处理装置参照所述基本频域参考信号、所述镜像频域参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值，包括：
所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号，Z(f)为频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
结合第三方面，在第二种可能的实现方式中，所述信号处理装置参照所述基本频域参考信号、所述镜像频域参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值，包括：
所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述信号处理装置依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述基本时域参考信号，s*(-t)为所述镜像时域参考信号，Z(t)为所述获取单元获取的时域的所述自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；Z(f)为Z(t)的傅里叶变换值。
结合第三方面的第一种可能的实现方式，在第三种可能的实现方式中，所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值包括：
所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述信号处理装置获取的第i帧的S(f)、S*(-f)和Z(f)之间的关系表达式，其中S′(f)＝S*(-f)；
将所述信号处理装置从第i帧起获取的N帧的S(f)、S*(-f)和Z(f)依次代入第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第三方面的第二种可能的实现方式，在第四种可能的实现方式中，所述信号处理装置依据第一代替式估算H1(f)的估算值和H2(f)的估算值包括：
所述信号处理装置依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述信号处理装置获取的第i帧的S(f)、S*(-f)和Z(f)之间的关系表达式，其中S′(f)＝S*(-f)；
将所述信号处理装置从第i帧起获取的N帧的S(f)、S*(-f)和Z(f)依次代入第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和
结合第三方面的第三种或第四种可能的实现方式，在第五种可能的实现方式中，所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，具体包括：
所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和具体包括：
所述信号处理装置结合第四公式与计算出和
结合第三方面的第三种或第四种可能的实现方式，在第六种可能的实现方式中，所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，具体包括：
所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和具体包括：
所述信号处理装置结合所述第四代替式与计算出和
结合第三方面或第一种至第六种可能的实现方式中的任一实现方式，在第七种可能的实现方式中，所述信号处理装置参照所述第一综合响应、所述第二综合响应、所述基本频域参考信号和所述镜像频域参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消，包括：
所述信号处理装置依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为频域的所述自干扰信号，为频域的所述第一综合相应的估算值，为频域的所述第二综合响应的估算值，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号。
本发明实施例提供的信号处理方法及装置，能计算出受IQ不平衡影响的自干扰信号用于自干扰抵消，从而在通信系统存在IQ不平衡时实现自干扰抵消，提高了数字干扰抵消性能。
附图说明
图1为本发明的实施例提供的一种信号处理装置的结构示意图；
图2为本发明的实施例提供的另一种信号处理装置的结构示意图；
图3为本发明的另一实施例提供的一种信号处理装置的结构示意图；
图4为本发明的实施例提供的一种信号处理方法的流程示意图；
图5为本发明的实施例提供的一种信号处理方法的运行示意图。
具体实施方式
本发明实施例的信号处理装置可以设置于或本身即为采用无线全双工技术的接入终端。接入终端也可以称为系统、用户单元、用户站、移动站、移动台、远方站、远程终端、移动设备、用户终端、终端、无线通信设备、用户代理、用户装置或用户设备(UE，UserEquipment)。接入终端可以是蜂窝电话、无绳电话、SIP(Session Initiation Protocol，会话启动协议)电话、WLL(Wireless Local Loop，无线本地环路)站、PDA(Personal DigitalAssistant，个人数字处理)、具有无线通信功能的手持设备、车载设备、可穿戴设备、计算设备或连接到无线调制解调器的其它处理设备。
此外，本发明实施例的信号处理装置还可以设置于或本身即为采用无线全双工技术的基站。基站可用于与移动设备通信，基站可以是WiFi的AP(Access Point，无线接入点)，或者是GSM(Global System of Mobile communication，全球移动通讯)或CDMA(CodeDivision Multiple Access，码分多址)中的BTS(Base Transceiver Station，基站)，也可以是WCDMA(Wideband Code Division Multiple Access，宽带码分多址)中的NB(NodeB，基站)，还可以是LTE(Long Term Evolution，长期演进)中的eNB或eNodeB(Evolutional NodeB，演进型基站)，或者中继站或接入点，或者未来5G网络中的基站设备等。
本发明的实施例提供一种信号处理装置，参照图1所示，用于实现通信领域的自干扰抵消，包括以下结构：
获取单元101，用于获取数字基带参考信号、自干扰参考信号和频域基带信号。
其中，所述自干扰参考信号包括频域的自干扰参考信号或者时域的自干扰参考信号。
转化单元102，用于参照所述获取单元101获取的所述数字基带参考信号获得基本参考信号和镜像参考信号。
其中，所述镜像参考信号为所述基本参考信号的镜像信号。所述基本参考信号包括基本频域参考信号和基本时域参考信号中的至少一种，所述镜像参考信号包括镜像频域参考信号和镜像时域参考信号中的至少一种。
估算单元103，用于参照所述转化单元102获得的所述基本参考信号、所述转化单元102获得的所述镜像参考信号和所述获取单元101获取的所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值。
其中，所述第一综合响应包括频域的第一综合响应和时域的第一综合响应中的至少一种，所述第二综合响应包括频域的第二综合响应和时域的第二综合响应中的至少一种，所述自干扰参考信号包括频域的自干扰参考信号和时域的自干扰参考信号中的至少一种。
可选的，所述估算单元103可以用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述转化单元102获得的所述基本频域参考信号，S*(-f)为所述转化单元102获得的所述镜像频域参考信号，Z(f)为所述获取单元101获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
具体的，当所述转化单元102获得所述基本频域参考信号S(f)、所述转化单元102获得所述镜像频域参考信号S*(-f)和所述获取单元101获取频域的所述自干扰参考信号Z(f)时，所述估算单元103可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算出H1(f)的估算值和H2(f)的估算值
可选的，所述估算单元103还可以用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述转化单元102获得的所述基本时域参考信号，s*(-t)为所述转化单元102获得的所述镜像时域参考信号，Z(t)为所述获取单元101获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)可以为H1(t)的傅里叶变换值，H2(f)可以为H2(t)的傅里叶变换值，表示卷积；H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
具体的，当所述转化单元102获得所述基本时域参考信号s(t)、所述转化单元102获得所述镜像时域参考信号s*(-t)和所述获取单元101获取时域的所述自干扰参考信号Z(t)时，所述估算单元103可以依据第一代替式估算出H1(t)的估算值和H2(t)的估算值
因为H1(t)可以为H1(f)对应的时域表达形式，H2(t)可以为H2(f)对应的时域表达形式，相应的可以为对应的时域表达形式，可以为对应的时域表达形式，所以可以对进行傅里叶变换或者以傅里叶级数的方式得到也可以对进行傅里叶变换或者以傅里叶级数的方式得到
其中，当以上信号为周期信号时，以傅立叶级数的形式分别将和转化为和当以上信号为非周期信号时，以傅立叶变换的形式分别将和转化为和此时，所述估算单元103还可以将所述第一代替式转化为第一公式，并将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，结合所述第一公式得到和其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
这里需要说明的是，所述估算单元103的工作内容不限于上述内容，其他可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者第一代替式估算出H1(f)的估算值和H2(f)的估算值的工作内容，均可以由所述估算单元103完成，这些工作内容可以为本发明所涵盖。
重构单元104，用于参照所述估算单元103估算的所述第一综合响应的估算值、所述估算单元103估算的所述第二综合响应的估算值、所述转化单元102获得的所述基本参考信号和所述转化单元102获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述获取单元101获取的所述频域基带信号进行自干扰抵消。
具体的，所述重构单元104可以用于依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
其中，为频域的所述自干扰信号，为所述估算单元103估算的频域的所述第一综合相应的估算值，为所述估算单元103估算的频域的所述第二综合响应的估算值，S(f)为所述转化单元102获得的所述基本频域参考信号，S*(-f)为所述转化单元102获得的所述镜像频域参考信号。
根据对所述估算单元103工作内容的上述描述，可选的，参照图2所示，所述估算单元103可以包括：
第一估算子单元1031，用于依据所述第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
其中所述第二公式为所述转化单元102获取的第i帧的S(f)、S*(-f)和所述获取单元101获取的第i帧Z(f)之间的关系表达式，S′(f)＝S*(-f)。
具体的，当所述转化单元102获得所述基本频域参考信号S(f)、所述转化单元102获得所述镜像频域参考信号S*(-f)和所述获取单元101获取频域的所述自干扰参考信号Z(f)时，所述第一估算子单元1031可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
更进一步，可选的，所述第一估算子单元1031，还可以用于依据所述第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
具体的，当所述转化单元102获得所述基本时域参考信号s(t)、所述转化单元102获得所述镜像时域参考信号s*(-t)和所述获取单元101获取时域的所述自干扰参考信号Z(t)时，所述第一估算子单元1031也可以依据第一代替式进行频域转化得到第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，并将将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，然后依据由第一代替式转换得来的第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
需要说明的是，所述第一估算子单元1031的工作内容不限于以上所述，凡由依据所述第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者所述第一代替式获取所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i的工作内容均可为所述第一估算子单元1031的工作内容，可以为本发明所涵盖。
第二估算子单元1032，用于将所述转化单元102从第i帧起获取的N帧内的S(f)、S*(-f)和所述获取单元101从第i帧起获取的N帧的Z(f)一次代入所述第一估算子单元1031获取的所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组，写为矩阵形式：表示为Z(f)＝S(f)H(f)，其中，
第三估算子单元1033，用于依据所述第二估算子单元1032计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
第四估算子单元1034，用于结合所述第三估算子单元1033计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
更进一步，可选的，所述第三估算子单元1033，具体用于将所述第二估算子单元1032计算出的矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元1034，具体用于结合所述第三估算子单元1033计算出的所述第四公式与计算出和
可选的，所述第三估算子单元1033还可以具体用于将所述第二估算子单元1032计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1为RLS(Recursive least-mean-square，递归最小均方差)算法中的遗忘因子，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元1034，同样还可以具体用于结合所述第三估算子单元1033计算出的所述第四代替式与计算出和
这里需要说明的是，因为所述第二估算子单元1032得到了N阶的线性方程组Z(f)＝S(f)H(f)，而通过这一方程组求解H(f)得到H(f)的估算值矩阵可以有很多算法，由不同求解方法得出的所述表示Z(f)、S(f)和之间关系的表达式也就有很多种。所述第三公式是根据LS(Least-mean-square，最小方差)算法得出的表示Z(f)、S(f)和之间关系表达式，所述第三代替式是根据RLS(Recursive least-mean-square，递归最小均方差)得到的表示Z(f)、S(f)和之间关系表达式，其他能表示Z(f)、S(f)和之间关系的表达式，也同样可以由所述第三估算子单元1033计算得出，以便所述第四估算子单元1034结合所述第三估算子单元1033计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和这些工作内容也同样可以被本发明所涵盖。
本发明实施例提供的信号处理装置，能计算出受IQ不平衡影响的自干扰信号用于自干扰抵消，从而在通信系统存在IQ不平衡时实现自干扰抵消，提高了数字干扰抵消性能。
本发明的实施例提供一种信号处理装置，参照图3所示，该信号处理装置可以嵌入或本身就是微处理计算机，比如：通用计算机、客户定制机、手机终端或平板机等便携设备，该信号处理装置3001包括：至少一个接收器3011、处理器3012、存储器3013和总线3014，该至少一个接收器3011、处理器3012和存储器313通过总线3014连接并完成相互间的通信。
该总线3014可以是ISA(Industry Standard Architecture，工业标准体系结构)总线、PCI(Peripheral Component，外部设备互连)总线或EISA(Extended IndustryStandard Architecture，扩展工业标准体系结构)总线等。该总线3014可以分为地址总线、数据总线、控制总线等。为便于表示，图3中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。其中：
存储器3013用于存储可执行程序代码，该程序代码包括计算机操作指令。存储器3013可能包括高速RAM存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
处理器3012可能是一个中央处理器(Central Processing Unit，简称为CPU)，或者是特定集成电路(Application Specific Integrated Circuit，简称为ASIC)，或者是被配置成实施本发明实施例的一个或多个集成电路。
其中，所述接收器3011，用于获取数字基带参考信号、自干扰参考信号和频域基带信号。
其中，所述自干扰参考信号包括频域的自干扰参考信号或者时域的自干扰参考信号。
所述处理器3012，用于参照所述接收器3011获取的所述数字基带参考信号获得基本参考信号和镜像参考信号。
其中，所述镜像参考信号为所述基本参考信号的镜像信号。所述基本参考信号包括基本频域参考信号和基本时域参考信号中的至少一种，所述镜像参考信号包括镜像频域参考信号和镜像时域参考信号中的至少一种。
所述处理器3012还用于参照所述处理器3012获得的所述基本参考信号、所述处理器3012获得的所述镜像参考信号和所述接收器3011获取的所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值。
其中，所述第一综合响应包括频域的第一综合响应和时域的第一综合响应中的至少一种，所述第二综合响应包括频域的第二综合响应和时域的第二综合响应中的至少一种，所述自干扰参考信号包括频域的自干扰参考信号和时域的自干扰参考信号中的至少一种。
可选的，所述处理器3012可以用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述处理器3012获得的所述基本频域参考信号，S*(-f)为所述处理器3012获得的所述镜像频域参考信号，Z(f)为所述接收器3011获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
具体的，当所述处理器3012获得所述基本频域参考信号S(f)、所述处理器3012获得所述镜像频域参考信号S*(-f)和所述接收器3011获取频域的所述自干扰参考信号Z(f)时，所述处理器3012可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算出H1(f)的估算值和H2(f)的估算值
可选的，所述处理器3012可以用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述处理器3012获得的所述基本时域参考信号，s*(-t)为所述处理器3012获得的所述镜像时域参考信号，Z(t)为所述接收器3011获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)可以为H1(t)的傅里叶变换值，H2(f)可以为H2(t)的傅里叶变换值，表示卷积；H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
具体的，当所述处理器3012获得所述基本时域参考信号s(t)、所述处理器3012获得所述镜像时域参考信号s*(-t)和所述接收器3011获取时域的所述自干扰参考信号Z(t)时，所述处理器3012可以依据第一代替式估算出H1(t)的估算值和H2(t)的估算值
因为H1(t)可以为H1(f)对应的时域表达形式，H2(t)可以为H2(f)对应的时域表达形式，相应的可以为对应的时域表达形式，可以为对应的时域表达形式，所以可以对进行傅里叶变换或者以傅里叶级数的方式得到也可以对进行傅里叶变换或者以傅里叶级数的方式得到
其中，当以上信号为周期信号时，以傅立叶级数的形式分别将和转化为和当以上信号为非周期信号时，以傅立叶变换的形式分别将和转化为和此时，所述处理器3012还可以将所述第一代替式转化为第一公式，并将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，结合所述第一公式得到和其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
这里需要说明的是，所述处理器3012的工作内容不限于以上两个方面，其他可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者第一代替式估算出H1(f)的估算值和H2(f)的估算值的工作内容，均可以由所述处理器3012完成，这些工作内容可以为本发明所涵盖。
所述处理器3012还用于参照所述处理器3012估算的所述第一综合响应的估算值、所述处理器3012估算的所述第二综合响应的估算值、所述处理器3012获得的所述基本参考信号和所述处理器3012获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置3001对所述接收器3011获取的所述频域基带信号进行自干扰抵消。
具体的，所述处理器3012还可以用于：依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
其中，为频域的所述自干扰信号，为所述处理器3012估算的频域的所述第一综合相应的估算值，为所述处理器3012估算的频域的所述第二综合响应的估算值，S(f)为所述处理器3012获得的所述基本频域参考信号，S*(-f)为所述处理器3012获得的所述镜像频域参考信号。
可选的，所述处理器3012具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值可以包括以下内容：
所述处理器3012具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
其中所述第二公式为所述处理器3012获取的第i帧的S(f)、S*(-f)和所述接收器3011获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)。
具体的，当所述处理器3012获得所述基本频域参考信号S(f)、所述处理器3012获得所述镜像频域参考信号S*(-f)和所述接收器3011获取频域的所述自干扰参考信号Z(f)时，所述处理器3012可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S(f)}i。
更进一步，可选的，所述处理器3012还可以具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
具体的，当所述处理器3012获得所述基本时域参考信号s(t)、所述处理器3012获得所述镜像时域参考信号s*(-t)和所述接收器3011获取时域的所述自干扰参考信号Z(t)时，所述处理器3012也可以依据第一代替式进行频域转化得到第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，并将将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，然后依据由第一代替式转换得来的第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
需要说明的是，所述处理器3012的工作内容不限于以上所述，凡由依据所述第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者所述第一代替式获取所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i的工作内容均可为所述处理器3012的工作内容，可以为本发明所涵盖。
所述处理器3012还用于将所述处理器3012从第i帧起获取的N帧的S(f)、S*(-f)和所述接收器3011从第i帧起获取的N帧的Z(f)代入所述处理器3012获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组，写为矩阵形式：
表示为Z(f)＝S(f)H(f)，其中，
所述处理器3012，还用于依据所述处理器3012计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述处理器3012，还用于结合所述处理器3012计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
更进一步，可选的，所述处理器3012还具体用于将上述矩阵Z(f)、S(f)和H(f)代入所述第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述处理器3012还具体用于结合所述处理器3012计算出的所述第四公式与计算出和
可选的，所述处理器3012还可以具体用于将所述处理器3012计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1为RLS(递归最小均方差)算法中的遗忘因子， SH(f)为S(f)的转置共轭矩阵，
所述处理器3012，同样还可以具体用于结合所述处理器3012计算出的所述第四代替式与计算出和
这里需要说明的是，因为所述处理器3012得到了N阶的线性方程组Z(f)＝S(f)H(f)，而通过这一方程组求解H(f)得到H(f)的估算值矩阵可以有很多算法，由不同求解方法得出的所述表示Z(f)、S(f)和之间关系的表达式也就有很多种。所述第三公式是根据LS(最小方差)算法得出的表示Z(f)、S(f)和之间关系表达式，所述第三代替式是根据RLS(递归最小均方差)得到的表示Z(f)、S(f)和之间关系表达式，其他能表示Z(f)、S(f)和之间关系的表达式，也同样可以由所述处理器3012计算得出，以便所述处理器3012结合所述处理器3012计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和这些工作内容也同样可以被本发明所涵盖。
本发明实施例提供的信号处理装置，能计算出受IQ不平衡影响的自干扰信号用于自干扰抵消，从而在通信系统存在IQ不平衡时实现自干扰抵消，提高了数字干扰抵消性能。
通过以上的实施方式的描述，所属领域的技术人员可以清楚地了解到本发明可以用硬件实现，或固件实现，或它们的组合方式来实现。当使用软件实现时，可以将上述功能存储在计算机可读介质中或作为计算机可读介质上的一个或多个指令或代码进行传输。计算机可读介质包括计算机存储介质和通信介质，其中通信介质包括便于从一个地方向另一个地方传送计算机程序的任何介质。存储介质可以是计算机能够存取的任何可用介质。以此为例但不限于：计算机可读介质可以包括RAM(Random Access Memory，随机存储器)、ROM(Read Only Memory，只读内存)、EEPROM(Electrically Erasable Programmable ReadOnly Memory，电可擦可编程只读存储器)、CD-ROM(Compact Disc Read Only Memory，即只读光盘)或其他光盘存储、磁盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质。此外，任何连接可以适当的成为计算机可读介质。例如，如果软件是使用同轴电缆、光纤光缆、双绞线、DSL(Digital Subscriber Line，数字用户专线)或者诸如红外线、无线电和微波之类的无线技术从网站、服务器或者其他远程源传输的，那么同轴电缆、光纤光缆、双绞线、DSL或者诸如红外线、无线和微波之类的无线技术包括在所属介质的定影中。如本发明所使用的，盘和碟包括CD(Compact Disc，压缩光碟)、激光碟、光碟、DVD碟(Digital Versatile Disc，数字通用光)、软盘和蓝光光碟，其中盘通常磁性的复制数据，而碟则用激光来光学的复制数据。上面的组合也应当包括在计算机可读介质的保护范围之内。
本发明的实施例提供一种信号处理方法，参照图4所示，用于实现通信领域的自干扰抵消，可以用于本发明的上述实施例提供的信号处理装置，为方便理解，可以同时参照图5所示，包括以下步骤：
401、信号处理装置获取数字基带参考信号、自干扰参考信号和频域基带信号。
其中，所述自干扰参考信号包括频域的自干扰参考信号或者时域的自干扰参考信号。
402、所述信号处理装置参照所述数字基带参考信号获得基本参考信号和镜像参考信号。
其中，所述镜像参考信号为所述基本参考信号的镜像信号。所述基本参考信号包括基本频域参考信号和者基本时域参考信号的至少一种，所述镜像参考信号包括镜像频域参考信号和镜像时域参考信号的至少一种。
具体的，因为所述镜像参考信号为所述基本参考信号的镜像信号，且参照所述数字基带参考信号获得所述基本参考信号的方法为现有技术，所以获得所述基本参考信号和所述镜像参考信号的具体方法在此不再赘述。
403、所述信号处理装置参照所述基本参考信号、所述镜像参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值。
其中，所述第一综合响应包括频域的第一综合响应和时域的第一综合响应中的至少一种，所述第二综合响应包括频域的第二综合响应和时域的第二综合响应中的至少一种，所述自干扰参考信号包括频域的自干扰参考信号和时域的自干扰参考信号中的至少一种。
可选的，所述信号处理装置可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值其中，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号，Z(f)为频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
具体的，当所述信号处理装置获得所述基本频域参考信号S(f)、所述镜像频域参考信号S*(-f)和频域的所述自干扰参考信号Z(f)时，可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算出H1(f)的估算值和H2(f)的估算值
可选的，所述信号处理装置可以依据第一代替式估算H1(f)的估算值和H2(f)的估算值其中，s(t)为所述基本时域参考信号，s*(-t)为所述镜像时域参考信号，Z(t)为所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)可以为H1(t)的傅里叶变换值，H2(f)可以为H2(t)的傅里叶变换值表示卷积。
具体的，当所述信号处理装置获得所述基本时域参考信号s(t)、所述镜像时域参考信号s*(-t)和时域的所述自干扰参考信号Z(t)时，所述信号处理装置可以依据第一代替式估算出H1(t)的估算值和H2(t)的估算值
因为H1(t)可以为H1(f)对应的时域表达形式，H2(t)可以为H2(f)对应的时域表达形式，相应的可以为对应的时域表达形式，可以为对应的时域表达形式，所以可以对进行傅里叶变换或者以傅里叶级数的方式得到也可以对进行傅里叶变换或者以傅里叶级数的方式得到
其中，当以上信号为周期信号时，以傅立叶级数的形式分别将和转化为和当以上信号为非周期信号时，以傅立叶变换的形式分别将和转化为和此时，所述信号处理装置还可以将所述第一代替式转化为第一公式，并将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，结合所述第一公式得到和其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
这里需要说明的是，所述信号处理方法不限于以上所述，其他可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者第一代替式估算出H1(f)的估算值和H2(f)的估算值的方法，均可以为本发明所涵盖。
进一步的，所述步骤403可以由以下具体步骤实现：
4031、所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
其中，所述第二公式为所述信号处理装置获取的第i帧的S(f)、S*(-f)和Z(f)之间的关系表达式，S′(f)＝S*(-f)。
具体的，当所述信号处理装置获得所述基本频域参考信号S(f)、所述镜像频域参考信号S*(-f)和频域的所述自干扰参考信号Z(f)时，可以依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
可选的，依据所述步骤403的可选方案，所述步骤4031还可以由下述步骤4031a代替：
4031a、所述信号处理装置依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
具体的，当所述信号处理装置获得所述基本时域参考信号s(t)、所述镜像时域参考信号s*(-t)和时域的所述自干扰参考信号Z(t)时，也可以依据第一代替式进行频域转化得到第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，并将将s(t)、s*(-t)和Z(t)分别进行频域转化，得到S(f)、S*(-f)和Z(f)，然后依据由第一代替式转换得来的第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，得出所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i。
需要说明的是，所述信号处理方法不限于以上所述，凡由依据所述第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)或者所述第一代替式获取所述第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i的方法均可以为本发明所涵盖。
4032、将所述信号处理装置从第i帧起获取的N帧的基本参考信号、镜像参考信号和自干扰参考信号依次代入第二公式，计算出N阶线性方程组。
其中步骤4032具体为将所述信号处理装置从第i帧起获取的N帧的S(f)、S*(-f)和Z(f)依次代入第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组，写为矩阵形式：表示为Z(f)＝S(f)H(f)，其中，
其中上述N阶线性方程组为矩阵Z(f)、S(f)和H(f)之间的关系。
4033、所述信号处理装置依据所述N阶线性方程组计算出表示自干扰参考信号矩阵、基本参考信号和镜像参考信号矩阵、综合响应估算值矩阵矩阵之间关系的第一表达式。
其中步骤4033具体为所述信号处理装置依据所述N阶线性方程组中的所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的第一表达式，其中，是H(f)的估算值矩阵，
4034、所述信号处理装置结合所述第一表达式和综合响应估算值矩阵计算出第一综合响应的估算值和第二综合相应的估算值。
其中，所述步骤4034具体为所述信号处理装置结合所述第一表达式和综合响应估算值矩阵计算出和
可选的，上述步骤4033和步骤4034可以由如下步骤4033a和步骤4034a具体实现，步骤4033a和步骤4034a之间的关系等同于步骤4033和步骤4034之间的关系：
4033a、所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
4034a、所述信号处理装置结合第四公式与计算出和
可选的，上述步骤4033和步骤4034也可以由如下步骤4033b和步骤4034b具体实现，步骤4033b和步骤4034b之间的关系等同于步骤4033和步骤4034之间的关系：
4033b、所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1为RLS(递归最小均方差)算法中的遗忘因子， SH(f)为S(f)的转置共轭矩阵，
4034b、所述信号处理装置结合所述第四代替式与计算出和
这里需要说明的是，因为所述步骤4032得到了N阶的线性方程组Z(f)＝S(f)H(f)，而通过这一方程组求解H(f)得到H(f)的估算值矩阵可以有很多算法，由不同求解方法得出的所述表示Z(f)、S(f)和之间关系的表达式也就有很多种。所述第三公式是根据LS(最小方差)算法得出的表示Z(f)、S(f)和之间关系表达式，所述第三代替式是根据RLS(递归最小均方差)得到的表示Z(f)、S(f)和之间关系表达式，其他能表示Z(f)、S(f)和之间关系的表达式，也同样可以由所述步骤4033计算得出，以便所述步骤4034结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和这些工作内容也同样可以被本发明所涵盖。
上述步骤中，为所述第一综合相应H1(f)的估算值，为所述第二综合响应H2(f)的估算值。
404、所述信号处理装置参照所述第一综合响应的估算值、所述第二综合响应的估算值、所述基本参考信号和所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
具体的，所述信号处理装置可以依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
其中，为频域的所述自干扰信号，为频域的所述第一综合相应的估算值，为频域的所述第二综合响应的估算值，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号。
需要说明的是，作为所述步骤403可选方案的实现基础，所述第一公式可以由如下过程推导：
所述信号发射机存在IQ不平衡时，发送的存在IQ不平衡的频域基带信号可以用第一参考公式X(f)＝G1(f)S(f)+G2(f)S*(-f)表示，其中，S(f)为所述基本频域参考信号，S*(-f)为镜像频域参考信号，G1(f)和G2(f)为包含所述信号发射机IQ通道幅度和相位不平衡的频率响应。
所述信号接收机不存在IQ不平衡时，接收的所述频域基带信号中自干扰信号Z(f)可用第二参考公式Z(f)＝G1(f)H(f)S(f)+G2(f)H(f)S*(-f)表示，其中，H(f)为自干扰信道的频域响应。
当所述信号接收机也存在IQ不平衡的情况时，参照第二参考公式，可知，所述信号接收机接收的频域基带信号中的自干扰信号Z(f)可以表示为第三参考公式Z(f)＝K1(f)H(f)X(f)+K2(f)H*(-f)X*(-f)，其中，K1(f)和K2(f)为包含信号接收机IQ通道幅度和相位不平衡的频率响应，X(f)为信号发射机发送的存在IQ不平衡的频域基带信号。
将第一参考公式X(f)＝G1(f)S(f)+G2(f)S*(-f)代入第三参考公式Z(f)＝K1(f)H(f)X(f)+K2(f)H*(-f)X*(-f)，可以得到第四参考公式：
令第四参考公式式中：
H1(f)＝K1(f)H(f)G1(f)+K2(f)H*(-f)G2 *(-f)，
H2(f)＝K1(f)H(f)G2(f)+K2(f)H*(-f)G1 *(-f)，
得到第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，其中，Z(f)为所述信号接收机接收的频域基带信号中的自干扰信号，等于频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，H2(f)为频域的所述第二综合响应，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号。
本发明实施例提供的信号处理方法，能通过分析自干扰参考信号和数字基带信号，计算出受IQ不平衡影响的自干扰信号用于自干扰抵消，从而在通信系统存在IQ不平衡时实现自干扰抵消，提高数字干扰抵消性能。
根据上述实施例可以知道的是本发明的实施例可以直接用于无线全双工系统中相互通信的信号发射机和信号接收机均存在IQ不平衡时的情景，特征是信号处理装置获取的是不存在IQ不平衡的数字基带参考信号、存在IQ不平衡的自干扰参考信号和存在IQ不平衡的频域基带信号。最后，计算存在IQ不平衡的自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
根据上述实施例还可以知道的是，本发明的实施例一种变通的实施例也可以具体用于在MIMO(Multiple-Input Multiple-Output，多输入多输出)场景下的无线全双工系统中相互通信的信号发射机和信号接收机均存在IQ不平衡时的情景，以2x2MIMO(2x2Multiple-Input Multiple-Output，双发射天线结合双接收天线的多输入多输出)场景下的无线全双工系统中的一个接收天线为例，包括以下内容：
参照所述步骤401，所述信号处理装置获取不存在IQ不平衡的数字基带参考信号、存在IQ不平衡的自干扰参考信号Z(f)和存在IQ不平衡的频域基带信号。其中，所述自干扰参考信号包括频域的自干扰参考信号Z(f)和时域的自干扰参考信号Z(t)中的至少一种。
参照所述步骤402，所述信号处理装置参照所述数字基带参考信号获得第一基本频域参考信号S(f)、第一镜像频域参考信号S*(-f)、第二基本频域参考信号Y(f)和第二镜像频域参考信号Y*(-f)。
其中，第一基本频域参考信号S(f)和第一镜像频域参考信号S*(-f)可以是在2x2MIMO场景中从第一信号接收链获取的，第二基本频域参考信号Y(f)和第二镜像频域参考信号Y*(-f)可以是2x2MIMO场景中从第二信号接收链获取的信号。所述第二基本频域参考信号Y(f)与所述第二镜像频域参考信号Y*(-f)之间的关系和S(f)与S*(-f)之间的关系类似，Y(f)的获得与S(f)类似，Y*(-f)的获得与S*(-f)类似。如果获取的是所述基本时域参考信号和所述镜像时域参考信号，则获取方法与此类同。以下主要以获得频域的各种信号为例进行说明。
参照所述步骤403及其可选方案，所述信号处理装置参照所述基本频域参考信号S(f)、所述镜像频域参考信号S*(-f)、所述第二基本频域参考信号Y(f)、所述第二镜像频域参考信号Y*(-f)和所述自干扰参考信号Z(f)，估算第一综合响应H1(f)的估算值第二综合响应H2(f)的估算值第三综合响应H3(f)的估算值和第四综合响应H4(f)的估算值
可选的，当所述信号处理装置获得第一基本频域参考信号S(f)、第一镜像频域参考信号S*(-f)、第二基本频域参考信号Y(f)、第二镜像频域参考信号Y*(-f)和频域的所述自干扰参考信号Z(f)时，可以依据所述第六公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)+H3(f)Y(f)+H4(f)Y*(-f)，估算出H1(f)的估算值H2(f)的估算值H3(f)的估算值和H4(f)的估算值
所述第六公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)+H3(f)Y(f)+H4(f)Y*(-f)的推导过程类似所述第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)的推导过程。
可选的，当所述信号处理装置获得第一基本时域参考信号、第一镜像时域参考信号、第二基本时域参考信号、第二镜像时域参考信号和时域的所述自干扰参考信号Z(t)时，第六公式时域转换后的其他表达式，如第六代替式估算出H1(t)的估算值H2(t)的估算值H3(t)的估算值和H4(t)的估算值
其中，第六公式的用法类似第一公式，和之间关系类似和之间关系，和之间关系类似和之间关系，Y(t)和Y(f)之间关系类似s(t)和S(f)之间关系，Y*(-t)和Y*(-f)之间关系类似s*(-t)和S*(-f)之间关系，Z(t)类似Z(f)，其中，Z(f)可以为Z(t)的傅里叶变换值，Z(f)为所述频域的所述自干扰参考信号。
参照步骤403及具体细化步骤4032，所述信号处理装置可以依据第六公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)+H3(f)Y(f)+H4(f)Y*(-f)，令S′(f)＝S*(-f)，Y′(f)＝Y*(-f)，计算出N阶线性方程组，写为矩阵形式：
表示为Z(f)＝S(f)H(f)，其中，
参照步骤步骤4033a，将本实施例中的上述矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，
参照步骤步骤4034a，结合第四公式与计算出和
可选的，上述参照步骤4033a和步骤4034a完成的工作，也可以参照步骤4033b和步骤4034b具体实现：
参照步骤4033b，所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1为RLS(递归最小均方差)算法中的遗忘因子，
参照步骤4034b，所述信号处理装置结合所述第四代替式与计算出和
这里需要说明的是，其他能可以通过参照所述步骤4033计算得出所述第一表达式，以便所述信号处理装置参照所述步骤4034结合所述第一表达式和计算出和的工作内容，也同样可以被本发明所涵盖。
最后，参照所述步骤404及实施例的上述内容，所述信号处理装置参照S(f)、S*(-f)、Y(f)和Y*(-f)结合第六公式及变形，计算存在IQ不平衡的频域的自干扰信号以便于所述信号处理装置对所述频域基带信号进行自干扰抵消。
由上可知，本发明实施例提供的信号处理方法，能计算出多输入多输出场景下的无线全双工系统中受IQ不平衡影响的自干扰信号用于自干扰抵消，从而在无线全双工系统中相互通信的信号发射机和信号接收机中的至少一个存在IQ不平衡时实现自干扰抵消，提高了无线全双工系统的数字干扰抵消性能。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (21)

1.一种信号处理装置，其特征在于，包括：
获取单元，用于获取数字基带参考信号、自干扰参考信号和频域基带信号；
转化单元，用于参照所述获取单元获取的所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
估算单元，用于参照所述转化单元获得的所述基本参考信号、所述转化单元获得的所述镜像参考信号和所述获取单元获取的所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；
重构单元，用于参照所述估算单元估算的所述第一综合响应的估算值、所述估算单元估算的所述第二综合响应的估算值、所述转化单元获得的所述基本参考信号和所述转化单元获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述获取单元获取的所述频域基带信号进行自干扰抵消；
其中，所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述估算单元具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述转化单元获得的所述基本频域参考信号，S*(-f)为所述转化单元获得的所述镜像频域参考信号，Z(f)为所述获取单元获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
2.根据权利要求1所述的装置，其特征在于，所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述估算单元具体用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述转化单元获得的所述基本时域参考信号，s*(-t)为所述转化单元获得的所述镜像时域参考信号，Z(t)为所述获取单元获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
3.根据权利要求1所述的装置，其特征在于，所述估算单元包括：
第一估算子单元，用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述转化单元获取的第i帧的S(f)、S*(-f)和所述获取单元获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)；
第二估算子单元，用于将所述转化单元从第i帧起获取的N帧的S(f)、S*(-f)和所述获取单元从第i帧起获取的N帧的Z(f)依次代入所述第一估算子单元获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
第三估算子单元，用于依据所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
第四估算子单元，用于结合所述第三估算子单元计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
4.根据权利要求2所述的装置，其特征在于，所述估算单元包括：
第一估算子单元，用于依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述转化单元获取的第i帧的S(f)、S*(-f)和所述获取单元获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)，Z(f)为Z(t)的傅里叶变换值，Z(f)为频域的所述自干扰参考信号；
第二估算子单元，用于将所述转化单元从第i帧起获取的N帧的S(f)、S*(-f)和所述获取单元从第i帧起获取的N帧的Z(f)依次代入所述第一估算子单元获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
第三估算子单元，用于依据所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
第四估算子单元，用于结合所述第三估算子单元计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
5.根据权利要求3或4所述的装置，其特征在于，包括：
所述第三估算子单元具体用于将所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元，具体用于结合所述第三估算子单元计算出的所述第四公式与计算出和
6.根据权利要求3或4所述的装置，其特征在于，包括：
所述第三估算子单元具体用于将所述第二估算子单元计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述第四估算子单元，具体用于结合所述第三估算子单元计算出的所述第四代替式与计算出和
7.根据权利要求1-4任一项所述的装置，其特征在于，所述重构单元具体用于：依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为频域的所述自干扰信号，为所述估算单元估算的频域的所述第一综合相应的估算值，为所述估算单元估算的频域的所述第二综合响应的估算值，S(f)为所述转化单元获得的所述基本频域参考信号，S*(-f)为所述转化单元获得的所述镜像频域参考信号。
8.一种信号处理装置，其特征在于，包括：接收器、处理器、存储器及总线，其中，所述接收器、所述处理器及所述存储器通过所述总线相互连接，
其中，所述接收器，用于获取数字基带参考信号、自干扰参考信号和频域基带信号；
所述处理器，用于参照所述接收器获取的所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
所述处理器还用于参照所述处理器获得的所述基本参考信号、所述处理器获得的所述镜像参考信号和所述接收器获取的所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；
所述处理器还用于参照所述处理器估算的所述第一综合响应的估算值、所述处理器估算的所述第二综合响应的估算值、所述处理器获得的所述基本参考信号和所述处理器获得的所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述接收器获取的所述频域基带信号进行自干扰抵消；
其中，所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述处理器获得的所述基本频域参考信号，S*(-f)为所述处理器获得的所述镜像频域参考信号，Z(f)为所述接收器获取的频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
9.根据权利要求8所述的装置，其特征在于，所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述处理器具体用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述处理器获得的所述基本时域参考信号，s*(-t)为所述处理器获得的所述镜像时域参考信号，Z(t)为所述接收器获取的所述时域的自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
10.根据权利要求8所述的装置，其特征在于，所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值包括：
所述处理器具体用于依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述处理器获取的第i帧的S(f)、S*(-f)和所述接收器获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)；
所述处理器还用于将所述处理器从第i帧起获取的N帧的S(f)、S*(-f)和所述接收器从第i帧起获取的N帧的Z(f)代入所述处理器获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述处理器还用于依据所述处理器计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述处理器还用于结合所述处理器计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
11.根据权利要求9所述的装置，其特征在于，所述处理器具体用于依据第一代替式估算H1(f)的估算值和H2(f)的估算值包括：
所述处理器具体用于依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述处理器获取的第i帧的S(f)、S*(-f)和所述接收器获取的第i帧的Z(f)之间的关系表达式，S′(f)＝S*(-f)，Z(f)为Z(t)的傅里叶变换值，Z(f)为频域的所述自干扰参考信号；
所述处理器还用于将所述处理器从第i帧起获取的N帧的S(f)、S*(-f)和所述接收器从第i帧起获取的N帧的Z(f)代入所述处理器获取的第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述处理器还用于依据所述处理器计算出的矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述处理器还用于结合所述处理器计算出的所述表示Z(f)、S(f)和之间关系的表达式和计算出和
12.根据权利要求10或11所述的装置，其特征在于，包括：
所述处理器具体用于将所述处理器计算出的矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述处理器还具体用于结合所述处理器计算出的第四公式与计算出和
13.根据权利要求10或11所述的装置，其特征在于，包括：
所述处理器具体用于将所述处理器计算出的矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述处理器还具体用于结合所述处理器计算出的所述第四代替式与计算出和
14.根据权利要求8-11任一项所述的装置，其特征在于，所述处理器还具体用于：依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为频域的所述自干扰信号，为所述处理器估算的频域的所述第一综合相应的估算值，为所述处理器估算的频域的所述第二综合响应的估算值，S(f)为所述处理器获得的所述基本频域参考信号，S*(-f)为所述处理器获得的所述镜像频域参考信号。
15.一种信号处理方法，其特征在于，包括：
信号处理装置获取数字基带参考信号、自干扰参考信号和频域基带信号；
所述信号处理装置参照所述数字基带参考信号获得基本参考信号和镜像参考信号，所述镜像参考信号为所述基本参考信号的镜像信号；
所述信号处理装置参照所述基本参考信号、所述镜像参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值；
所述信号处理装置参照所述第一综合响应的估算值、所述第二综合响应的估算值、所述基本参考信号和所述镜像参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，所述基本参考信号包括基本频域参考信号；所述镜像参考信号包括镜像频域参考信号；所述第一综合响应包括频域的第一综合响应；所述第二综合响应包括频域的第二综合响应；所述自干扰参考信号包括频域的自干扰参考信号；
所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值
其中，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号，Z(f)为频域的所述自干扰参考信号，H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
16.根据权利要求15所述的方法，其特征在于，所述信号处理装置参照所述基本频域参考信号、所述镜像频域参考信号和所述自干扰参考信号，估算第一综合响应的估算值和第二综合响应的估算值，包括：
所述基本参考信号包括基本时域参考信号；所述镜像参考信号包括镜像时域参考信号；所述第一综合响应包括时域的第一综合响应；所述第二综合响应包括时域的第二综合响应；所述自干扰参考信号包括时域的自干扰参考信号；
所述信号处理装置依据第一代替式估算H1(f)的估算值和H2(f)的估算值
其中，s(t)为所述基本时域参考信号，s*(-t)为所述镜像时域参考信号，Z(t)为所述获取单元获取的时域的所述自干扰参考信号，H1(t)为时域的所述第一综合响应，H2(t)为时域的所述第二综合响应，其中H1(f)为H1(t)的傅里叶变换值，H2(f)为H2(t)的傅里叶变换值；H1(f)为频域的所述第一综合响应，为频域的所述第一综合相应H1(f)的估算值，H2(f)为频域的所述第二综合响应，为频域的所述第二综合响应H2(f)的估算值。
17.根据权利要求15所述的方法，其特征在于，所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)，估算H1(f)的估算值和H2(f)的估算值包括：
所述信号处理装置依据第一公式Z(f)＝H1(f)S(f)+H2(f)S*(-f)获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述信号处理装置获取的第i帧的S(f)、S*(-f)和Z(f)之间的关系表达式，S′(f)＝S*(-f)；
将所述信号处理装置从第i帧起获取的N帧的S(f)、S*(-f)和Z(f)依次代入第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和
18.根据权利要求16所述的方法，其特征在于，所述信号处理装置依据第一代替式估算H1(f)的估算值和H2(f)的估算值包括：
所述信号处理装置依据第一代替式获取第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，其中所述第二公式为所述信号处理装置获取的第i帧的S(f)、S*(-f)和Z(f)之间的关系表达式，S′(f)＝S*(-f)，Z(f)为Z(t)的傅里叶变换值，Z(f)为频域的所述自干扰参考信号；
将所述信号处理装置从第i帧起获取的N帧的S(f)、S*(-f)和Z(f)依次代入第二公式{Z(f)}i＝H1(f){S(f)}i+H2(f){S′(f)}i，计算出N阶线性方程组Z(f)＝S(f)H(f)，其中，
所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和
19.根据权利要求17或18所述的方法，其特征在于，所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，具体包括：
所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三公式计算出第四公式其中，SH(f)为S(f)的转置共轭矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和具体包括：
所述信号处理装置结合第四公式与计算出和
20.根据权利要求17或18所述的方法，其特征在于，所述信号处理装置依据所述矩阵Z(f)、S(f)和H(f)计算出表示Z(f)、S(f)和之间关系的表达式，其中，是H(f)的估算值矩阵，具体包括：
所述信号处理装置将矩阵Z(f)、S(f)和H(f)代入第三代替式计算出第四代替式其中，0＜λ＜1，SH(f)为S(f)的转置共轭矩阵，
所述信号处理装置结合所述表示Z(f)、S(f)和之间关系的表达式和计算出和具体包括：
所述信号处理装置结合所述第四代替式与计算出和
21.根据权利要求15-18任一项所述的方法，其特征在于，所述信号处理装置参照所述第一综合响应、所述第二综合响应、所述基本频域参考信号和所述镜像频域参考信号来计算自干扰信号，以便于所述信号处理装置对所述频域基带信号进行自干扰抵消，包括：
所述信号处理装置依据第五公式计算以便于所述信号处理装置对所述频域基带信号进行自干扰抵消；
其中，为频域的所述自干扰信号，为频域的所述第一综合相应的估算值，为频域的所述第二综合响应的估算值，S(f)为所述基本频域参考信号，S*(-f)为所述镜像频域参考信号。
CN201480073837.7A 2014-05-09 2014-05-09 一种信号处理方法及装置 Active CN106416168B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2014/077169 WO2015168941A1 (zh) 2014-05-09 2014-05-09 一种信号处理方法及装置

Publications (2)

Publication Number Publication Date
CN106416168A CN106416168A (zh) 2017-02-15
CN106416168B true CN106416168B (zh) 2019-10-18

Family

ID=54392022

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480073837.7A Active CN106416168B (zh) 2014-05-09 2014-05-09 一种信号处理方法及装置

Country Status (3)

Country Link
US (1) US10193684B2 (zh)
CN (1) CN106416168B (zh)
WO (1) WO2015168941A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017022961A1 (ko) * 2015-07-31 2017-02-09 엘지전자 주식회사 Fdr 방식을 이용하는 통신 장치가 비선형 자기간섭 신호의 채널 추정을 위한 참조신호를 전송하는 방법
EP3343857B1 (en) * 2015-08-26 2020-03-18 LG Electronics Inc. Method for changing rs mode in environment where operation is done by fdr scheme, and device therefor
CN112787962A (zh) * 2020-12-28 2021-05-11 上海复旦通讯股份有限公司 基于导频参考信号的信道时域特征提取方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101785190A (zh) * 2007-08-02 2010-07-21 Lm爱立信电话有限公司 Iq失衡图像抑制
CN102057644A (zh) * 2008-04-11 2011-05-11 爱立信电话股份有限公司 用于估计iq不平衡的方法和布置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7248625B2 (en) * 2002-09-05 2007-07-24 Silicon Storage Technology, Inc. Compensation of I-Q imbalance in digital transceivers
JP4261578B2 (ja) * 2006-12-27 2009-04-30 株式会社東芝 無線通信装置及び受信方法
US8565352B2 (en) * 2010-05-03 2013-10-22 Telefonaktiebolaget L M Ericsson (Publ) Digital IQ imbalance compensation for dual-carrier double conversion receiver
US9184902B2 (en) * 2012-04-25 2015-11-10 Nec Laboratories America, Inc. Interference cancellation for full-duplex communications

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101785190A (zh) * 2007-08-02 2010-07-21 Lm爱立信电话有限公司 Iq失衡图像抑制
CN102057644A (zh) * 2008-04-11 2011-05-11 爱立信电话股份有限公司 用于估计iq不平衡的方法和布置

Also Published As

Publication number Publication date
WO2015168941A1 (zh) 2015-11-12
US10193684B2 (en) 2019-01-29
US20170063518A1 (en) 2017-03-02
CN106416168A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN106664273B (zh) Mimo同信道全双工收发器中的干扰抵消
Hakkarainen et al. Widely-linear beamforming and RF impairment suppression in massive antenna arrays
US10038471B2 (en) Method and apparatus for canceling self-interference
CN113055119B (zh) 全双工多小区网络中的干扰消除方法及装置
CN105704339B (zh) 用于消除非线性失真的变换核心
JP2017529715A (ja) 干渉除去装置および方法
KR101517424B1 (ko) 전송 다이버시티 방법, 및 관련 장치 및 시스템
WO2016155464A1 (en) Joint radio-frequency/baseband self-interference cancellation methods
EP3065361B1 (en) Digital interference cancellation apparatus and method for wireless full duplex system and transceiver
CN108270541A (zh) 传送宽带宽数据帧的装置、方法和系统
CN109428624B (zh) 信号处理方法、网络设备、终端及计算机可读存储介质
CN106416168B (zh) 一种信号处理方法及装置
CN113924734A (zh) 使用单个tci状态的多波束操作
WO2021169819A1 (en) Method and device for uplink transmission
CN107852381A (zh) 用于在无线通信系统中估计非线性自干扰信道的方法及其装置
CN111480298B (zh) 通过电荷共享的混叠抑制
CN105393458B (zh) 一种干扰抵消方法和装置
WO2021081847A1 (zh) 一种信道测量方法和通信装置
CN106256143A (zh) 用于干扰消除的方法和终端
CN107395222A (zh) 消除失真的方法及用于用户设备的电路
US11909455B2 (en) Global equalizer self interference cancelation for hybrid MIMO systems
US11050495B2 (en) Electronic device including transceiver for calibrating I/Q imbalance in millimeter wave communication system and method of operating same
CN111130682B (zh) 一种信号处理方法及装置
CN115118358B (zh) 干扰源角度的估计方法、装置、设备和存储介质
WO2021062758A1 (zh) 用于多输入多输出mimo通信的装置及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant