CN101355367A - 干扰消除的方法和装置 - Google Patents

Abstract

本发明的实施例公开了一种干扰消除方法和装置，其通过设置滤波器系数向量的边界约束条件、或滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件为非零常数；并根据该边界约束条件和代价函数获取滤波器系数向量和滤波后信道冲击响应系数向量。通过使用本发明的实施例，设置滤波器系数向量的边界约束条件、或滤波后信道冲击响应系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，可以在较低的复杂度下得到较好的性能。

Description

干扰消除的方法和装置

技术领域

本发明涉及通讯技术领域，尤其涉及一种干扰消除的方法和装置。

背景技术

为了增加GSM(Global System for Mobile Communications，全球移动通信系统)系统网络的性能，现有技术中提出了DARP(Downlink Advanced ReceiverPerformance，下行链路高级接收机性能)解决方案，通过使用较低的频率复用因子有效的提高网络容量，但同时必须使用抑制邻小区的干扰消除算法。

现有技术中的干扰消除算法基于MMSE(Min Mean Square Error，最小均方误差)准则实现，将接收信号经过滤波器p滤波后，令下列(1)式的代价函数最小：

$J=E\left\{{\left(\underset{n=0}{\overset{\mathrm{Qp}-1}{\mathrm{\Σ}}}p\right[n\left]r\right[k-n]-\underset{n=0}{\overset{\mathrm{Qh}-1}{\mathrm{\Σ}}}h[n\left]d(k-n-k_0)\right)}^2\right\}---\left(1\right)$

其中：r[k]为时刻k的接收信号，d[k]为发射数据符号，p[n]为滤波器系数，h[n]为滤波后信道冲击响应系数，Q_p为滤波器抽头数，Q_h为信道冲击响应抽头数，k₀为固定延时。

将上述代价函数(1)写为矩阵的形式，即：

J＝||Rp-Dh||²                       (2)

其中，R为接收信号矩阵，D为数据符号矩阵，p为滤波器系数向量，h为滤波后信道冲击响应系数向量。

其中 $\text{D=}\left(\begin{array}{ccc}d(k-k_0)& ...& d(k-k_0-Q_d)\\ ...& ...& ...\\ d(k-k_0+K)& ...& d(k-k_0+K-Q_d)\end{array}\right)$ $h=\left(\begin{array}{c}h\left(0\right)\\ h\left(1\right)\\ .\\ .\\ .\\ h\left(Q_d\right)\end{array}\right)---\left(3\right)$

现有算法都是在这个准则的基础上，建立方程进而求解滤波器系数p。

发明人在实现本发明的过程中，发现现有技术至少存在以下问题：

要在J＝||Rp-Dh||²代价函数下求得p，h，需要对p或h进行约束，以避免零解，下面以约束h为例进行说明，当然也可以对p进行约束。

约束||h||²＝1，则：

p＝(R^TR)^-1R^TDh                          (4)

h＝eigvec_min(D^TD-D^TR(R^TR)^-1R^TD)         (5)

其中eigvec_min(A)表示矩阵A最小特征值对应的特征向量。

发明人在实现本发明的过程中，发现现有技术中使用反幂法以及迭代求解最小特征值对应的特征向量时，需要较大的复杂度及实现时间。

发明内容

本发明的实施例提供一种干扰消除的方法和装置，用于在较低的复杂度下得到较好的干扰消除性能。

本发明的实施例提供一种干扰消除的方法，包括：

设置滤波后信道冲击响应系数向量的边界约束条件，所述边界条件中滤波后信道冲击响应系数向量h中的一个元素h(i)的值为非零常数；

根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系；

根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

本发明的实施例还提供一种干扰消除方法，包括：

设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数；

根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系；

根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

本发明的实施例还提供一种干扰消除装置，包括：

第一设置单元，用于设置滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件为非零常数；

第一变换单元，用于根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系；

第一获取单元，用于根据所述第一变换单元获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第一干扰消除单元，用于根据所述第一获取单元获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

本发明的实施例还提供一种干扰消除装置，包括：

第二设置单元，用于设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数；

第二变换单元，用于根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系；

第二获取单元，用于根据所述第二变换单元获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第二干扰消除单元，用于根据所述第二获取单元获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

与现有技术相比，本发明的实施例具有以下优点：

通过设置滤波器系数向量的边界约束条件、或滤波后信道冲击响应系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，在较低的复杂度下得到较好的性能。

附图说明

图1是本发明实施例中干扰消除方法的流程图；

图2是本发明实施例中干扰消除方法的另一流程图；

图3是本发明实施例中干扰消除装置的结构示意图；

图4是本发明实施例中干扰消除装置的另一结构示意图。

具体实施方式

本发明的实施例中提供一种干扰消除方法，通过设置滤波后信道冲击响应系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，可以在较低的复杂度下得到较好的性能。如图1所示，该干扰消除方法包括：

步骤s101、设置滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件中滤波后信道冲击响应系数向量h中的一个元素h(i)的值为非零常数。

步骤s102、根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系。

步骤s103、根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量。

步骤s104、根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

以下结合一个具体的应用场景，描述本发明的实施例中干扰消除方法的具体实现。以设置滤波后信道冲击响应系数向量的边界约束条件为例，该约束条件为：对于滤波后信道冲击响应系数向量h，设置边界约束条件为滤波后信道冲击响应系数向量h中的一个元素的值为1。这里只是以边界约束条件为1进行说明，在实际应用中设置为其它非零常数(包括实数和复数)均可。

以约束h(0)＝1为例，则代价函数Rp＝Dh+n可变形为：

$d=\mathrm{Rp}-\tilde{D}\tilde{h}+n---\left(6\right)$

其中n为干扰噪声向量，比D少第一列，比h少了第一个元素，d为D的第一列：

$\tilde{D}=\left(\begin{array}{ccc}d(k-k_0-1)& ...& d(k-k_0-Q_d)\\ ...& ...& ...\\ d(k-k_0+K-1)& ...& d(k-k_0+K-Q_d)\end{array}\right)$ $\tilde{h}=\left(\begin{array}{c}h\left(1\right)\\ h\left(2\right)\\ .\\ .\\ .\\ h\left(Q_d\right)\end{array}\right)---\left(7\right)$

可将p，h看作一个向量，上式写为：

$d=(R,-\tilde{D}){(p^T,{\tilde{h}}^T)}^T+n---\left(8\right)$

然后由LS(Linear-least Square，最小二乘)算法，可以得：

${(p^T,{\tilde{h}}^T)}^T={\left({(R,-\tilde{D})}^T(R,-\tilde{D})\right)}^{-1}{(R,-\tilde{D})}^{T}d---\left(9\right)$

上述式(9)中，左侧矩阵

的两部分分别与p，h相关，而式(9)的右侧表达式中，d为发射数据符号，R为接收信号矩阵，D为数据符号矩阵，即式(9)中右侧的各变量均为已知。因此根据式(9)的右侧表达式的计算结果，可以对应的得到滤波器系数向量p、滤波后信道冲击响应系数向量h的值。根据该滤波器系数向量p、滤波后信道冲击响应系数向量h的值进行干扰消除。

上述实施例以设置的约束条件为h(0)＝1为例，说明本发明实施例中的干扰消除方法。除了该h(0)＝1的约束条件外，还可以设置其他约束条件，以设置h(i)＝1(i为大于0的正整数)的约束条件为例，则代价函数Rp＝Dh+n可变形为：

$d=\mathrm{Rp}-\tilde{D}\tilde{h}+n---\left(10\right)$

其中n为干扰噪声向量，比D少第i列，

比h少了第i个元素，d为D的第i列：

$\tilde{D}=\left(\begin{array}{cccccc}d(k-k_0)& ...& d(k-k_0-i+1)& d(k-k_0-i-1)& ...& d(k-k_0-Q_d)\\ ...& ...& ...& ...& ...& ...\\ d(k-k_0+K)& ...& d(k-k_0+K-i+1)& d(k-k_0+K-i-1)& ...& d(k-k_0+K-Q_d)\end{array}\right)$

$\tilde{h}=\left(\begin{array}{c}h\left(0\right)\\ .\\ .\\ .\\ h(i-1)\\ h(i+1)\\ .\\ .\\ .\\ h\left(Q_d\right)\end{array}\right)---\left(11\right)$

可将p，h看作一个向量，上式写为：

$d=(R,-\tilde{D}){(p^T,{\tilde{h}}^T)}^T+n---\left(12\right)$

然后由LS算法，可以得：

${(p^T,{\tilde{h}}^T)}^T={\left({(R,-\tilde{D})}^T(R,-\tilde{D})\right)}^{-1}{(R,-\tilde{D})}^{T}d---\left(13\right)$

根据该式(13)同样可以对应的得到滤波器系数向量p、滤波后信道冲击响应系数向量h的值。

本发明的另一实施例中还提供一种干扰消除方法，通过设置滤波器系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，可以在较低的复杂度下得到较好的性能。如图2所示，该干扰消除方法包括：

步骤s201、设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数。

步骤s202、根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系。

步骤s203、根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量。

步骤s204、根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

以下结合一个具体的应用场景，描述本发明的实施例中干扰消除方法的具体实现。以对于滤波器系数向量p设置边界约束条件为例，设置边界约束条件为滤波器系数向量p中的一个元素的值为1。这里只是以边界约束条件为1进行说明，在实际应用中设置为其它非零常数(包括实数和复数)均可。

例如设置p(0)＝1，则代价函数Rp＝Dh+n可变形为：

$r=\tilde{R}\tilde{p}-\mathrm{Dh}+n---\left(14\right)$

使用上述实施例中相似的步骤，同样可以得到：

${({\tilde{p}}^T,h^T)}^T={\left({(\tilde{R},-D)}^T(\tilde{R},-D)\right)}^{-1}{(\tilde{R},-D)}^{T}d---\left(15\right)$

其中n为干扰噪声向量，

比R少第一列，比r少了第一个元素，r为R的第一列。

上述式(15)中，左侧矩阵

的两部分分别与p，h相关，而式(15)的右侧表达式中，d为发射数据符号，R为接收信号矩阵，D为数据符号矩阵，即式(15)中右侧的各变量均为已知。因此根据式(15)的右侧表达式的计算结果，可以对应的得到滤波器系数向量p、滤波后信道冲击响应系数向量h的值。根据该滤波器系数向量p、滤波后信道冲击响应系数向量h的值进行干扰消除。

上述实施例以设置的约束条件为p(0)＝1为例，说明本发明实施例中的干扰消除方法。除了该p(0)＝1的约束条件外，还可以设置其他约束条件，以设置p(j)＝1(j为大于0的正整数)的约束条件为例，则代价函数Rp＝Dh+n可变形为：

$r=\tilde{R}\tilde{p}-\mathrm{Dh}+n---\left(16\right)$

使用上述实施例中相似的步骤，同样可以得到：

${({\tilde{p}}^T,h^T)}^T={\left({(\tilde{R},-D)}^T(\tilde{R},-D)\right)}^{-1}{(\tilde{R},-D)}^{T}d---\left(17\right)$

其中n为干扰噪声向量，

比R少第j列，

比r少了第j个元素，r为R的第j列。

根据该式(17)同样可以对应的得到滤波器系数向量p、滤波后信道冲击响应系数向量h的值。

本发明的实施例提供的干扰消除方法中，通过设置滤波器系数向量的边界约束条件、或滤波后信道冲击响应系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，在较低的复杂度下得到较好的性能。以本发明中的LS算法为例，设矩阵阶数为N，则LS算法的求解复杂度与N³相当。而对于现有技术中的反幂法和迭代法，反幂法的求解复杂度与N³相当，一次迭代的求解复杂度与N²相当。与现有技术中一次求解需要一个反幂法运算以及至少数十次迭代相比，本发明实施例提供的方法显著降低了求解的复杂度。

本发明的实施例中还提供一种干扰消除装置，如图3所示，包括：

第一设置单元11，用于设置滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件中滤波后信道冲击响应系数向量h中的一个元素h(i)的值为非零常数；

第一变换单元12，用于根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系；

第一获取单元13，用于根据所述第一变换单元12获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第一干扰消除单元14，用于根据所述第一获取单元13获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

具体的，本发明的实施例中：

第一设置单元11可以具体用于：对于滤波后信道冲击响应系数向量h，设置边界约束条件为滤波后信道冲击响应系数向量h中的一个元素h(i)的值为1。

第一变换单元12可以具体用于：对代价函数进行变换，获取到所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系具体为：

$d=\mathrm{Rp}-\tilde{D}\tilde{h}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵；所述约束条件下的滤波后信道冲击响应系数向量

比h少第i个元素，所述约束条件下的数据符号矩阵

比D少第i列，d为D的第i列。

第一获取单元13可以具体用于，获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${(p^T,{\tilde{h}}^T)}^T={\left({(R,-\tilde{D})}^T(R,-\tilde{D})\right)}^{-1}{(R,-\tilde{D})}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

本发明的实施例还提供一种干扰消除装置，如图4所示，包括：

第二设置单元21，用于设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数；

第二变换单元22，用于根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系；

第二获取单元23，用于根据第二变换单元22获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第二干扰消除单元24，用于根据第二获取单元23获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

具体的，本发明的实施例中：

第二设置单元21可以具体用于：对于滤波器系数向量p，设置边界约束条件为滤波器系数向量p中的一个元素p(j)的值为1。

第二变换单元22可以具体用于，对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系具体为：

$r=\tilde{R}\tilde{p}-\mathrm{Dh}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵；所述约束条件下的滤波器系数向量比p少第j个元素，所述约束条件下的接收信号矩阵

比R少第j列，r为R的第j列。

第二获取单元23可以具体用于，获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${({\tilde{p}}^T,h^T)}^T={\left({(\tilde{R},-D)}^T(\tilde{R},-D)\right)}^{-1}{(\tilde{R},-D)}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

本发明的实施例提供的干扰消除装置中，通过设置滤波器系数向量的边界约束条件、或滤波后信道冲击响应系数向量的边界约束条件，可以快速简便的获取滤波器系数向量和滤波后信道冲击响应系数向量并进行干扰消除，在较低的复杂度下得到较好的性能。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本发明可以通过硬件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，本发明的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM，U盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。

以上公开的仅为本发明的几个具体实施例，但是，本发明并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。

Claims (16)

1、一种干扰消除方法，其特征在于，包括：

设置滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件中滤波后信道冲击响应系数向量h中的一个元素h(i)的值为非零常数；

根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系；

根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

2、如权利要求1所述的方法，其特征在于，所述滤波后信道冲击响应系数向量的边界约束条件具体为：

对于滤波后信道冲击响应系数向量h，设置边界约束条件为滤波后信道冲击响应系数向量h中的一个元素h(i)的值为1。

3、如权利要求2所述的方法，其特征在于，所述获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系具体为：

$d=\mathrm{Rp}-\tilde{D}\tilde{h}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵，n为干扰噪声向量；所述约束条件下的滤波后信道冲击响应系数向量

比h少第i个元素，所述约束条件下的数据符号矩阵比D少第i列，d为D的第i列。

4、如权利要求3所述的方法，其特征在于，所述根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${(p^T,{\tilde{h}}^T)}^T={\left({(R,-\tilde{D})}^T(R,-\tilde{D})\right)}^{-1}{(R,-\tilde{D})}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

5、一种干扰消除方法，其特征在于，包括：

设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数；

根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系；

根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

根据所述获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

6、如权利要求5所述的方法，其特征在于，所述滤波器系数向量的边界约束条件具体为：

对于滤波器系数向量p，设置边界约束条件为滤波器系数向量p中的一个元素p(j)的值为1。

7、如权利要求5所述的方法，其特征在于，所述获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系具体为：

$r=\tilde{R}\tilde{p}-\mathrm{Dh}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵，n为干扰噪声向量；所述约束条件下的滤波器系数向量

比p少第j个元素，所述约束条件下的接收信号矩阵

比R少第j列，r为R的第j列。

8、如权利要求7所述的方法，其特征在于，根据所述关系获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${({\tilde{p}}^T,h^T)}^T={\left({(\tilde{R},-D)}^T(\tilde{R},-D)\right)}^{-1}{(\tilde{R},-D)}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

9、一种干扰消除装置，其特征在于，包括：

第一设置单元，用于设置滤波后信道冲击响应系数向量的边界约束条件，所述边界约束条件中滤波后信道冲击响应系数向量h中的一个元素h(i)的值为非零常数；

第一变换单元，用于根据所述约束条件对代价函数进行变换，获取所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系；

第一获取单元，用于根据所述第一变换单元获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第一干扰消除单元，用于根据所述第一获取单元获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

10、如权利要求9所述的干扰消除装置，其特征在于，所述第一设置单元具体用于：对于滤波后信道冲击响应系数向量h，设置边界约束条件为滤波后信道冲击响应系数向量h中的一个元素h(i)的值为1。

11、如权利要求10所述的干扰消除装置，其特征在于，所述第一变换单元具体用于，对代价函数进行变换，获取到所述约束条件下的滤波后信道冲击响应系数向量、所述约束条件下的数据符号矩阵、滤波器系数向量以及接收信号矩阵的关系具体为：

$d=\mathrm{Rp}-\tilde{D}\tilde{h}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵，n为干扰噪声向量；所述约束条件下的滤波后信道冲击响应系数向量比h少第i个元素，所述约束条件下的数据符号矩阵

比D少第i列，d为D的第i列。

12、如权利要求11所述的干扰消除装置，其特征在于，其特征在于，所述第一获取单元具体用于，获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${(p^T,{\tilde{h}}^T)}^T={\left({(R,-\tilde{D})}^T(R,-\tilde{D})\right)}^{-1}{(R,-\tilde{D})}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

13、一种干扰消除装置，其特征在于，包括：

第二设置单元，用于设置滤波器系数向量的边界约束条件，所述边界约束条件中滤波器系数向量p中的一个元素p(j)的值为为非零常数；

第二变换单元，用于根据所述约束条件对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系；

第二获取单元，用于根据所述第二变换单元获取的关系获取滤波器系数向量和滤波后信道冲击响应系数向量；

第二干扰消除单元，用于根据所述第二获取单元获取的滤波器系数向量和滤波后信道冲击响应系数向量进行干扰消除。

14、如权利要求13所述的干扰消除装置，其特征在于，所述第二设置单元具体用于：对于滤波器系数向量p，设置边界约束条件为滤波器系数向量p中的一个元素p(j)的值为1。

15、如权利要求14所述的干扰消除装置，其特征在于，所述第二变换单元具体用于，对代价函数进行变换，获取滤波后信道冲击响应系数向量、数据符号矩阵、所述约束条件下的滤波器系数向量以及所述约束条件下的接收信号矩阵的关系具体为：

$r=\tilde{R}\tilde{p}-\mathrm{Dh}+n$

其中h为滤波后信道冲击响应系数向量，D为数据符号矩阵，p为滤波器系数向量，R为接收信号矩阵，n为干扰噪声向量；所述约束条件下的滤波器系数向量

比p少第j个元素，所述约束条件下的接收信号矩阵比R少第j列，r为R的第j列。

16、如权利要求15所述的干扰消除装置，其特征在于，所述第二获取单元具体用于，获取滤波器系数向量和滤波后信道冲击响应系数向量具体为：

${({\tilde{p}}^T,h^T)}^T={\left({(\tilde{R},-D)}^T(\tilde{R},-D)\right)}^{-1}{(\tilde{R},-D)}^{T}d$

根据所述

的值获取滤波器系数向量p和滤波后信道冲击响应系数向量h。

Images