CN101507151B - 用于对接收信号进行分类的方法和装置 - Google Patents
用于对接收信号进行分类的方法和装置 Download PDFInfo
- Publication number
- CN101507151B CN101507151B CN2007800304753A CN200780030475A CN101507151B CN 101507151 B CN101507151 B CN 101507151B CN 2007800304753 A CN2007800304753 A CN 2007800304753A CN 200780030475 A CN200780030475 A CN 200780030475A CN 101507151 B CN101507151 B CN 101507151B
- Authority
- CN
- China
- Prior art keywords
- function
- signal
- matrix
- eigenvalue
- value
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/345—Interference values
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
- Circuits Of Receivers In General (AREA)
Abstract
本发明的实施例提供了一种用于对接收信号进行分类的方法。该方法包括:确定接收信号的信号值的协方差矩阵,并确定该协方差矩阵的特征值矩阵。所述特征值矩阵包括所述协方差矩阵的特征值。根据所述协方差矩阵的特征值中的至少一个特征值确定第一函数。根据所述协方差矩阵的特征值中的至少一个特征值确定第二函数,其中所述第二函数不同于所述第一函数。根据所述第一函数的值和所述第二函数的值之间的比较结果,接收信号被分类为包括数据的信号或噪声信号。
Description
技术领域
本发明的实施方式一般涉及对接收信号进行分类的方法,具体涉及将接收信号分类为包括数据的信号或噪声信号的方法。本发明的另外一种实施方式涉及分类器和计算机程序产品。
背景技术
认知无线电能够感测宽频带上的频谱环境并且充分利用该信息以提供最好地满足用户的通信需求的无线链路。由于无线电并非对于任何预分配的频率具有初选权,所以认知无线电操作的特殊之处在于要求无线电能够感测大频谱范围上的环境并适配于该环境。即,必须动态地检测主要用户的信号的存在。
能量(EG)检测是主要且基本的信号检测方法(感测算法)。与相干检测不同,能量检测不需要待检测信号的任何信息并且对未知的多径衰减具有鲁棒性。然而,因为能量检测依赖于对精确的噪声功率的了解,所以该方法易受噪声不确定性的影响。存在几个不确定性:(1)元件的非线性;(2)元件中的热噪声(非均匀的、时变的);(3)由于其它用户(无意的(附近的)或有意的(远处的))的发射所致的噪声。因此实际上,难以(几乎不可能)获得精确的噪声功率。
发明内容
本发明的一个实施例提供了一种对接收信号进行分类的方法。该方法包括:确定接收信号的信号值的协方差矩阵,并确定该协方差矩阵的特征值矩阵。所述特征值矩阵包括所述协方差矩阵的特征值。由所述协方差矩阵的特征值中的至少一个特征值确定第一函数。由所述协方差矩阵的特征值中的至少一个特征值确定第二函数,其中所述第二函数不同于所述第一函数。根据所述第一函数的值和所述第二函数的值之间的比较结果,接收信号被分类为包括数据的信号或噪声信号。
发明内容
本发明的一个实施例提供了一种对接收信号进行分类的方法。该方法包括:确定接收信号的信号值的协方差矩阵,并确定该协方差矩阵的特征值矩阵。所述特征值矩阵包括所述协方差矩阵的特征值。由所述协方差矩阵的特征值中的至少一个特征值确定第一函数。由所述协方差矩阵的特征值中的至少一个特征值确定第二函数,其中所述第二函数不同于所述第一函数。根据所述第一函数的值和所述第二函数的值之间的比较结果,接收信号被分类为包括数据的信号或噪声信号。
本发明的另一个实施例提供了一种对接收信号进行分类的分类器。该分类器包括用于确定接收信号的信号值的协方差矩阵的第一确定电路以及用于确定所述协方差矩阵的特征值矩阵的第二确定电路。所述特征值矩阵包括所述协方差矩阵的特征值。该分类器还包括用于由所述协方差矩阵的特征值中的至少一个特征值确定第一函数的第三确定电路,以及用于由所述协方差矩阵的特征值中的至少一个特征值确定第二函数的第四确定电路。所述第二函数确定为不同于所述第一函数。该分类器还包括用于根据所述第一函数的值和所述第二函数的值之间的比较结果将接收信号分类成包括数据的信号或噪声信号的分类电路。
附图说明
在附图中,在不同的附图中相似的附图标记通常指代相同的部分。这些附图不必按比例绘制,而重点在于用来解释本发明的原理。在以下说明中,参照附图详细说明本发明的各实施例,其中:
图1示出了表示根据本发明一个实施例的方法的流程图;
图2示出了表示根据本发明另一个实施例的方法的流程图;
图3A和图3B示出了表示根据本发明其它实施例的方法的流程图;
图4示出了根据本发明实施例的分类器;
图5示出了由无线麦克信号的仿真得出的检测概率;
图6A~图6F示出了由不同的DTV信号的仿真得出的检测概率;
图7示出了由多输入多接收机***中的仿真得出的检测概率;以及
图8示出了由单输入多接收机***中的仿真得出的检测概率。
具体实施方式
本发明的各实施例基于数据信号和噪声信号之间的统计特性的差异。例如,根据由接收信号所确定的特征值分布特性,本发明的实施例将接收信号分成包括数据的信号或噪声信号。根据本发明各实施例的方法和装置可用于各种信号检测场合而无需了解信号、信道和噪声的功率。
本发明的一个实施例涉及用于对接收信号进行分类的方法。该方法包括确定接收信号的信号值的协方差矩阵以及确定所述协方差矩阵的特征值矩阵。所述特征值矩阵包括所述协方差矩阵的特征值。第一函数由所述协方差矩阵的特征值中的至少一个特征值确定。第二函数由所述协方差矩阵的特征值中的至少一个特征值确定,其中所述第二函数不同于所述第一函数。根据第一函数的值和第二函数的值之间的比较结果,接收信号被分成包括数据的信号或噪声信号。
在一个实施例中,接收信号被采样。例如,接收信号可能是连续信号,并且被采样器采样为离散信号。采样器可以设置为以选定的采样率对接收信号进行采样。在一个实施例中,当想要处于具有预定带宽的频带中的信号时,采样率可选定为不小于所述预定带宽。在另一个实施例中,由接收机所接收的接收信号可以远高于对输入至接收机的输入信号进行采样的采样率的采样率被采样。
在一个实施例中,接收信号可以被滤波器滤波。该滤波器可以是模拟滤波器或数字滤波器。在一个实施例中,该滤波器可以是具有预定带宽的带通滤波器。该滤波器也可选择其它类型的滤波器。在另一个实施例中,已采样的接收信号可被滤波器滤波,该滤波器可以是模拟滤波器或数字滤波器。
根据本发明的实施例,接收信号的信号值x(n)的协方差矩阵R(NS)由以下等式确定:
其中NS表示协方差矩阵中信号值x(n)的个数。接收信号的信号值x(n)的协方差矩阵R(NS)也可以由不同形式的等式来确定。
其中Q表示厄密(Hermitian)矩阵。特征值矩阵 也可以根据本发明另一个实施例中的不同形式的公式来确定。在一个实施例中,厄密(Hermitian)矩阵Q可以由对接收信号进行滤波的滤波器来确定。在另一个实施例中,厄密(Hermitian)矩阵Q可以由噪声的协方差矩阵来确定。
根据实施例,第一函数是多个特征值的组合。在一个示例中,第一函数是多个特征值的线性组合,在另一个示例中,第一函数是多个特征值的非线性组合,例如为二次多项式、三次多项式或更高次的多项式。在一个实施例中,第一函数至少包括特征值矩阵的最大特征值。在另一个实施例中,第一函数包括特征值矩阵的预定个数的最大特征值,例如包括特征值矩阵的二个、三个、四个、五个…最大特征值。
根据另一个实施例,确定第一函数的值包括确定特征值矩阵的迹。因此,第一函数的值被确定为特征值矩阵的迹。
第二函数是多个特征值的组合,该组合可以是多个特征值的线性组合或者例如二次多项式、三次多项式或更高次的多项式的非线性组合。在一个实施例中,第二函数至少包括特征值矩阵的最小特征值。在另一个实施例中,第二函数包括特征值矩阵的预定个数的最小特征值,例如包括特征值矩阵的二个、三个、四个、五个…最小特征值。
根据本发明的实施例,利用比例因子在第一函数的值和第二函数的值之间进行比较。在一个实施例中,比例因子可以选定为满足误报概率的要求,例如获得不高于预定的误报概率的误报概率。根据一个实施例, 若第一函数的值与第二函数的值的比值大于选定的比例因子,则接收信号被分类为包括数据的信号。否则,接收信号被分类为噪声信号。所选定的比例因子也可称为阈值。
本发明的另一个实施例涉及对接收信号进行分类的分类器。所述分类器包括用于确定接收信号的信号值的协方差矩阵的第一确定电路以及用于确定协方差矩阵的特征值矩阵的第二确定电路。所述特征值矩阵包括协方差矩阵的特征值。所述分类器也可以包括由协方差矩阵的特征值中的至少一个特征值来确定第一函数的第三确定电路以及从协方差矩阵的特征值中的至少一个特征值来确定第二函数的第四确定电路。第二函数被确定为不同于第一函数。所述分类器还包括根据第一函数的值和第二函数的值之间的比较结果将接收信号分成包括数据的信号或噪声信号的分类电路。
在本说明书的上下文中,电路可以理解为包括例如实现为专用集成电路(ASIC)或硬件逻辑门结构的硬件电路,或者例如现场可编程门阵列(FPGA)的可编程电路,或例如实现为微处理器的任何其它的可编程逻辑电路,这些可编程电路被编程以提供相应的电路功能。
在一个实施例中,所述分类器还包括用于对接收信号进行采样的采样器。所述采样器可以设置为以选定的采样率对接收信号进行采样。在一个实施例中,当想要处于具有预定带宽的频带中的信号时,采样率可选为不小于预定带宽。在另一个实施例中,采样器可以设置为以远高于对输入至接收机的输入信号进行采样的采样率的采样率来采样接收信号。
在实施例中,所述分类器包括用于对接收信号进行滤波的滤波器。该滤波器可以是模拟滤波器或数字滤波器。在一个实施例中,该滤波器可选为具有预定带宽的带通滤波器。该滤波器也可选择其它类型的滤波器。在另一个实施例中,该滤波器用于对已采样的接收信号进行滤波。
根据一个实施例,第一确定电路设置为根据以下等式确定接收信号的信号值x(n)的协方差矩阵R(NS):
其中NS表示协方差矩阵中信号值x(n)的个数。第一确定电路也可设置为根据不同形式的等式来确定接收信号的协方差矩阵。
其中Q表示厄密(Hermitian)矩阵。在其它的示例中也可以选择不同的等式来确定特征值矩阵。在一个实施例中,厄密(Hermitian)矩阵Q可以由对接收信号进行滤波的滤波器来确定。在另一个实施例中,厄密(Hermitian)矩阵Q可以由噪声的协方差矩阵来确定。
第三确定电路可将第一函数确定为多个特征值的组合。在一个示例中,第一函数被确定为多个特征值的线性组合,在另一个示例中,第一函数被确定为多个特征值的非线性组合。在一个实施例中,第一函数至少包括特征值矩阵的最大特征值。在另一个实施例中,第一函数包括特征值矩阵的预定个数的最大特征值,例如包括特征值矩阵的二个、三个、四个、五个…最大特征值。
根据另一个实施例,第三确定电路设置为利用特征值矩阵的迹确定第一函数的值。于是,第一函数的值被确定为特征值矩阵的迹。
第四确定电路可将第二函数确定为多个特征值的组合,该组合可以是多个特征值的线性组合或者非线性组合。在一个实施例中,第二函数至少包括特征值矩阵的最小特征值。在另一个实施例中,第二函数包括特征值矩阵的预定个数的最小特征值,例如包括特征值矩阵的二个、三个、四个、五个…最小特征值。
在一个实施例中,分类电路设置为在比较第一函数的值和第二函数的值时利用比例因子。比例因子可以选为满足预定的标准。例如,可以通过选择比例因子以实现误报概率不高于预定的误报概率。根据一个实 施例,若第一函数的值与第二函数的值的比值大于选定的比例因子,则分类电路将接收信号分类为包括数据的信号。否则,分类电路将接收信号分类为噪声信号。
本发明的又一个实施例涉及包括如上所述的分类器的电子装置。在一个实施例中,电子装置还包括用于接收信号的接收机。在另一个实施例中,电子装置可包括用于接收信号的多个接收机。
在一个实施例中,电子装置设置为选自一组电子装置中的电子装置。这些电子装置的例子包括通信终端或设置为播放至少一种媒体的媒体播放装置,例如视频播放装置。
本发明的又一个实施例涉及用于对接收信号进行分类的计算机程序产品。所述计算机程序产品被设置用于当处理器执行该计算机程序产品时对接收信号进行分类。接收信号的分类包括确定接收信号的信号值的协方差矩阵以及确定协方差矩阵的特征值矩阵。所述特征值矩阵包括协方差矩阵的特征值。第一函数由协方差矩阵的特征值中的至少一个特征值确定。第二函数由协方差矩阵的特征值中的至少一个特征值确定,其中,第二函数不同于第一函数。根据第一函数的值和第二函数的值之间的比较结果,接收信号被分成包括数据的信号或噪声信号。
所述计算机程序产品可以任何编程语言来实现,从而可以通过执行以选定的编程语言写成的多条指令对接收信号进行分类。
图1示出了表示根据本发明一个实施例的方法的流程图。对于接收信号,在102中确定接收信号的信号值的协方差矩阵。然后在104中确定协方差矩阵的特征值矩阵,其中,该特征值矩阵包括协方差矩阵的特征值。在106中,由协方差矩阵的特征值中的至少一个特征值确定第一函数。在108中,由协方差矩阵的特征值中的至少一个特征值确定第二函数。在110中,根据第一函数的值和第二函数的值之间的比较结果,接收信号被分成包括数据的信号或噪声信号。
图2示出了表示根据本发明另一个实施例的方法的流程图。在202中,接收信号被采样并滤波。然后在204中确定比例因子γ。在一个实施例中,比例因子γ的确定可以基于选定的平滑因子和预定的误报概率。在206中, 确定已采样和滤波的接收信号的信号值的协方差矩阵,在208中确定协方差矩阵的特征值矩阵。特征值矩阵包括多个特征值,这些特征值例如以λ1≥λ2≥…≥λL的顺序排列。在210中,由协方差矩阵的特征值中的至少一个特征值确定第一函数F1(λ1,λ2,…,λL)和第二函数F2(λ1,λ2,…,λL)。然后在212中确定第一函数的值Γ1和第二函数的值Γ2。最后在214中,比较值Γ1和Γ2。若Γ1/Γ2>γ,则接收信号被分类为包括数据的信号。否则,接收信号被分类为噪声信号。
图3A示出了表示根据本发明实施例的方法的流程图。在302中接收信号被采样并滤波,在304中确定比例因子γ。在306中,确定已采样并滤波的接收信号的信号值的协方差矩阵,在308中确定协方差矩阵的特征值矩阵。特征值矩阵包括协方差矩阵的多个特征值。在310中,确定最大特征值λmax和最小特征值λmin。最后在312中,比较值λmax和λmin以对接收信号进行分类。若λmax/λmin>γ,则接收信号被分类为包括数据的信号。否则,接收信号被分类为噪声信号。
图3B示出了表示根据本发明又一个实施例的方法的流程图。在352中,接收信号被采样并滤波。在354中确定平滑因子L和比例因子γ,其中在一个实施例中,比例因子可以根据平滑因子和预定的误报概率来确定。在356中,确定已采样并滤波的接收信号的信号值的协方差矩阵,在358中确定协方差矩阵的特征值矩阵。特征值矩阵包括协方差矩阵的多个特征值。在360中,确定特征值矩阵的迹ρ和特征值矩阵的最小特征值λmin。最后在362中,比较迹ρ和最小特征值λmin以对接收信号进行分类。若ρ/λmin>γL,则接收信号被分类为包括数据的信号。否则,接收信号被分类为噪声信号。
图4示出了根据本发明实施例的用于对接收信号401进行分类的分类器400。分类器400包括用于确定接收信号401的信号值的协方差矩阵的第一确定电路402以及用于确定协方差矩阵的特征值矩阵的第二确定电路404。所述特征值矩阵包括协方差矩阵的特征值。分类器400还可包括用于由协方差矩阵的特征值中的至少一个特征值来确定第一函数的第三确定电路406以及用于由协方差矩阵的特征值中的至少一个特征值来确定第二函数的第四确定电路408。第二函数被确定为不同于第一函数。分类 器400还包括分类电路410,该分类电路410根据第一函数的值和第二函数的值之间的比较结果将接收信号分成包括数据的信号或噪声信号。
分类器400以及各电路402、404、406、408和410可以是硬件电路或硬件逻辑门结构,或者可编程电路,或任何其它的可编程逻辑电路,这些可编程电路被编程以提供如上所述的相应的电路功能。分类器400可被包括在电子装置中,所述电子装置用于对该电子装置所接收的信号进行分类。
在另一个实施例中,分类器可以是包括指令的计算机程序产品,该计算机程序产品被设置用于当处理器执行该计算机程序产品时根据本发明的上述实施例对接收信号进行分类。
信号和噪声的统计学协方差矩阵通常是不同的。噪声的统计学协方差矩阵是通过接收滤波器来确定的。因而其结构对于接收机是已知的。根据该结构,接收信号的协方差矩阵可以被变换成特征值矩阵。当接收信号是噪声信号时,特征值矩阵的最大和最小特征值相等。然而,当接收信号包括数据时,特征值矩阵的最大、平均和最小特征值是不同的。根据所谓的随机矩阵理论(RMT),特征值分布特性用于信号检测。于是,根据本发明的实施例的方法可以应用于各种信号检测场合,而无需了解信号、信道和噪声的功率。
示例1
在示例性示例中,令y(t)为时间连续的接收信号。假设需要中心频率为fc且带宽为W的频带,则接收信号y(t)以采样率fs被采样,其中fs≥W。令Ts=1/fs为采样周期。则接收到的离散信号为x(n)=y(nTs)。存在两个假设:H0:不存在数据;和H1:存在数据。因此在两个假设下的接收信号采样分别如以下所示:
H0:x(n)=η(n)
H1:x(n)=s(n)+η(n),
其中s(n)为所发送的信号,例如是通过无线信道(包括衰减和多路效应)的信号,η(n)是白噪声采样。应当指出,s(n)可以是多个信号的叠加。
接收信号通常经过带通滤波器。令f(k),k=0,1,…,K为 的带通滤波器。令g(n)为滤波器的自相关,即,
定义
滤波之后,接收信号变为:
令
则
选择平滑因子L并定义
其中NS表示信号值x(n)的个数。
定义L×(L+K)矩阵为
令G=HH H。将该矩阵分解为G=QQ H,其中Q为L×L厄密(Hermitian)矩阵。
通过以下步骤说明根据本发明的一个实施例的对接收信号进行分类的方法:
步骤1:如上所述对接收信号进行采样和滤波。
步骤2:选择平滑因子L并确定比例因子γ(也称为阈值),其中比例因子γ应确定为满足误报概率的要求。
步骤3:根据以下等式确定接收信号的信号值的协方差矩阵R(NS):
步骤5:确定接收信号的特征值矩阵的特征值并将这些特征值表示为λ1≥λ2≥…≥λL。
步骤6:由至少一个特征值确定第一函数F1和第二函数F2:
Γ1(Ns)=F1(λ1,λ2,…,λL)
Γ2(Ns)=F2(λ1,λ2,…,λL)
其中F1和F2是具有L个变量的正函数,Γ1(Ns)和Γ2(Ns)分别是第一函数和第二函数的值。函数F1和F2可以是多个特征值的组合,例如线性组合。函数F1和F2也可以分别包括多个最大和最小特征值。
步骤7:根据Γ1(Ns)和Γ2(Ns)之间的比较结果对接收信号进行分类。若Γ1(Ns)/Γ2(Ns)>γ,则存在信号(即接收信号包括数据)。否则,不存在信号(即接收信号为噪声信号)。
示例2
该示例称为最大最小特征值(MME)检测,其中最大特征值与最小特征值的比较结果用于对接收信号进行分类,对其说明如下。
步骤1:如以上示例1中所述对接收信号进行采样和滤波。
步骤2:选择平滑因子L并确定比例因子γ(也称为阈值)。根据以下等式,比例因子γ确定为满足误报概率的要求:
其中F1为所谓的1阶Tracy-Wisdom分布,P0为预定的误报概率,α为正数。α应选为满足误报概率的要求。Tracy-Wisdom分布的值在表1中给出。
t | -3.90 | -3.18 | -2.78 | -1.91 | -1.27 | -0.59 | 0.45 | 0.98 | 2.02 |
F1(t) | 0.01 | 0.05 | 0.10 | 0.30 | 0.50 | 0.70 | 0.90 | 0.95 | 0.99 |
表1:1阶Tracy-Wisdom分布的数值表
在后续部分中将说明用于比例因子的等式(1)的确定。
步骤3:确定协方差矩阵:
步骤4:变换协方差矩阵以确定特征值矩阵:
步骤5:确定特征值矩阵 的最大特征值和最小特征值并分别表示为λmax和λmin。
步骤6:根据最大特征值和最小特征值的比较结果确定信号的存在。若λmax/λmin>γ,则存在信号。否则,不存在信号。
示例3
由于特征值的平均值与协方差矩阵的迹成比例,所以可将迹和最小特征值进行比较以检测接收信号中是否包括数据。本发明的这一实施例称为具有最小特征值的迹(TME)检测,以下详细说明。因为迹表示平均能量,所以TME检测也称为具有最小特征值的能量(EME)检测。
步骤1:如以上示例1中所述对接收信号进行采样和滤波。
步骤2:选择平滑因子L并确定比例因子γ(也称为阈值)。比例因子γ确定为满足误报概率的要求。
步骤3:确定协方差矩阵:
步骤4:变换协方差矩阵以确定特征值矩阵:
步骤6:根据迹和最小特征值的比较结果确定信号的存在。若ρ/λmin>γL,则接收信号被分类为包括数据的信号。否则,接收信号被分类为噪声信号。
下文给出了利用根据以上示例2和示例3中的本发明实施例的方法针对无线麦克信号和DTV信号的仿真。
针对无线麦克信号的仿真
在本发明的实施例中,FM调制的无线麦克信号可记作:
其中fc为中心频率,fΔ为频率偏差,wm(τ)为源信号。
频率偏差选为100KHz。wm(τ)生成为在(-1,1)中均匀分布的实数。我们假定信号已下转换为中心频率为fIF=5.381119MHz的IF(中频)。采样率为21.524476MHz。带宽为6MHz的带通滤波器为具有89抽头(tap)的升余弦滤波器。信号与白噪声经过相同的滤波器。感测时间为9.30毫秒(ms)。平滑因子选为L=10。根据预定的所要求的误报概率Pfa=0.1设定比例因子或阈值(基于随机矩阵理论),Pfa=0.1对于信号是固定的。比例因子与噪声功率无关。
图5示出了根据仿真得出的检测概率,其中“EG-xdB”表示具有x-dB噪声不确定性的能量检测。下面的表2示出了根据能量检测的方法和根据本发明实施例的方法得出的误报概率,其中MME是以上示例2中的最大最小检测,TME是以上示例3中的具有最小特征值的迹检测。
EG-2dB | EG-1.5dB | EG-1dB | EG-0.5dB | EG-0dB (无不确定性) | TME/EME | MME |
0.497 | 0.497 | 0.496 | 0.483 | 0.108 | 0.081 | 0.086 |
表2 误报概率(感测时间9.63ms)
从图5和表2的仿真结果中可以看出,根据本发明实施例的方法得到了比能量检测方法更高的检测概率和更低的误报概率。在图5中,MME甚至优于具有确切噪声功率(无噪声不确定性)的能量检测。实际上,噪声不确定性通常是存在的。如图5所示,若存在大于0.5dB的噪声不确定性,则能量检测的检测概率比根据本发明实施例的方法的检测概率差得多。
而且,根据本发明实施例的方法以及无噪声不确定性的能量检测方法的误报概率满足误报概率(Pfa<=0.1)的要求。但是,如表2所示,利用无噪声不确定性的能量检测方法所得的误报概率远远超出了限定值。因此,与根据本发明实施例的方法相比,在实际的存在噪声不确定性的情况下能量检测方法是不可靠的。
用于所获取的DTV信号的仿真
该仿真是基于“频谱感测仿真模型”的。
所获取的DTV信号经过升余弦滤波器(带宽为6MHz,滚降系数为1/2,89抽头)。白噪声被叠加以获得不同的SNR等级。所使用的采样次数为400000(对应于18.60ms)。平滑因子选为L=10。根据所要求的误报概率Pfa=0.1设定阈值(根据随机矩阵理论),Pfa=0.1对于信号是固定的。阈值与噪声功率无关。
表3示出了根据能量检测方法和根据本发明实施例的方法得出的误报概率。图6A~图6F示出了对于不同DTV信号的检测概率。
EG-2dB | EG-1.5dB | EG-1dB | EG-0.5dB | EG-0dB (无不确定性) | TME/EME | MME |
0.496 | 0.496 | 0.491 | 0.481 | 0.095 | 0.029 | 0.077 |
表3 误报概率(白噪声,感测时间18.60ms)
从表3以及图6A~图6F的仿真结果中可以看出,根据本发明实施例的TME/EME和MME方法的性能优于能量检测方法。
因此,在本发明的实施例中,提供了在无线电或通信***中感测或检测信号的方法,例如用于在事先不了解信号、信号所在的信道以及接收信号的噪声功率的情况下检测接收信号中数据的存在。
以下,在示例4中说明了在多输入多输出(MIMO)***中本发明的实施例。
示例4
假定存在P个信号sj(n),(j=1,2,…,P)和M个接收机。若接收信号以输入信号的采样率的M倍的采样率被采样,即接收信号被过采样,则这也可适用于单接收机的情形。以下,这两种情况予以相同对待。接收机i的接收信号为xi(n),(i=1,2,…,M),
其中hij(k)为从用户j至天线i的信道响应,Nij为信道hij(k)的阶次,ηi(n)为信道噪声。应当指出,信号、信道响应和噪声可以是实数或虚数。令 如果必要则对hij(k)补零,并定义
表达式 n=0,1,…可以改写为如下的矢量形式:
考虑L连续输出并定义
则
其中H是定义为以下形式的ML×(N+PL) 矩阵
对已发送的符号以及信道噪声的统计属性作出如下假设:
(1)噪声为白噪声,
(2)噪声和已发送的数据信号是非互相关的。
令R(Ns)为如下形式的接收信号的样本协方差矩阵:
其中Ns为所收集的样本的个数。若Ns较大,则根据以上的假设(1)和(2),可以确定为
其中R s为输入信号的统计协方差矩阵, ση 2为噪声的方差,I ML为ML阶的单位矩阵。
令R和HR s H H的特征值分别为λ1≥λ2≥…≥λML和ρ1≥ρ2≥…≥ρML。显然,
若不存在信号,即 则 如存在信号,则ρ1>ρML,从而λ1>λML。因此,通过判断λ1/λML>γ(γ是阈值)是否成立,便可以检测信号的存在。
通过以下步骤说明根据本发明实施例的最大最小特征值(MME)检测对接收信号进行分类的方法:
步骤1:根据以下等式确定接收信号的信号值的协方差矩阵R(Ns):
步骤2:确定协方差矩阵R(Ns)的最大和最小特征值,即λmax和λmin。
步骤3:根据最大特征值和最小特征值的比较结果对接收信号进行分类。若λmax/λmin>γ,则存在信号(即接收信号包括数据)。否则,不存在信号(即接收信号为噪声信号)。这里γ>1,并且在下文中说明γ的确定。
比例因子或阈值γ的确定
实际上,通常只有有限个数的样本是可用的。因此,样本协方差矩阵R(Ns)可能是与统计协方差矩阵R相去甚远。R(Ns)的特征值分布变得很 复杂。令Pd为检测概率,即若存在信号时,将接收信号中包括数据判断为“是”的概率。令Pfa为误报的概率,即不存在信号时,将接收信号中包括数据判断为“是”的概率。对于好的检测算法而言,Pd应当大而Pfa应当小。显然地,阈值γ越大,Pfa越小,同时Pd也越小。这是一种“零和”博弈。阈值γ的选择是Pd与Pfa之间的折衷。在根据本发明实施例的方法中,设定了Pfa的限制,即Pfa≤P0。则可以确定极限阈值γ0,从而Pfa=P0。因此,满足γ≥γ0的阈值可以被确定为以上根据本发明实施例的示例1-4中所使用的阈值或比例因子。下面解释γ0的确定。
当不存在信号时,R(Ns)变为R η(Ns),噪声的样本协方差矩阵定义为
R η(Ns)为专用维希特(Wishart)随机矩阵。近年来,在数学以及通信和物理领域中,对随机矩阵的谱(特征值分布)的研究是很热门的课题。已经知道了随机矩阵R η(Ns)的有序特征值的联合概率密度函数(PDF)。但是,由于PDF的表达式很复杂,所以对于有序特征值的边缘PDF还没有找到闭式表达式。在以下定理中分别说明实数和复数矩阵的最大特征值的分布。
定理1:假设噪声为实数。令
定理2:假设噪声为复数。令
并且
应当指出,对于大的Ns,μ和μ′,ν和ν′几乎相等。这意味着,实数和复数矩阵的最大特征值的均值和方差几乎相等。但是,它们的极限分布是不同的。
在以下定理3中说明最小特征值的极限。
定理3:假设 (0<y<1),则
这里,ML视为固定的并且Ns趋于无穷大。根据定理,可得出以下结果:
令F1为1阶Tracy-Wisdom分布的累积分布函数(CDF)(有时简称为分布函数)。那么,对于实数信号,误报概率为
为了得到Pfa≤P0,阈值应当确定为满足下式
由此导出
或者,等同于,
根据μ和ν的定义,可根据以下等式确定阈值或比例因子γ:
对于复数信号,区别只在于函数F1应当由F2代替,F2为2阶Tracy-Wisdom分布的CDF。
Tracy-Wisdom分布被认为是某些随机矩阵的最大特征值的极限定律。对于分布函数没有闭式表达式。但是,根据数值计算,以上示例2的表1给出了F1在某些点的值。表1可用于计算某些点的F1 -1(y)。例如,
因此,在如上所述示例1-4的本发明的实施例中所使用的阈值或比例因子γ可以根据这一节来确定,其中γ的等式依据不同的情形而稍作修改,例如修改为示例2中所示的等式(1),其中信号是由单接收机接收的。
根据示例4的在多接收机***中的仿真
在下文中,根据以下等式,信噪比(SNR)表示平均的接收信号功率与平均的噪声功率的比:
10000个样本用于计算样本协方差矩阵。所有结果通过1000次蒙特卡罗实现得到平均值。预定的要求的误报概率Pfa≤0.1(P0=0.1)。那么对于MME方法,根据等式(2)确定阈值γ0。为了进行比较,也仿真了针对同一***的能量(EG)检测。
对于能量检测,例如根据Dyspan 2005中A.Sahai和D.Cabric的″Spectrum sensing:fundamental limits and practical challenges″所述的理论确定阈值。能量检测需要预先知道噪声功率。由于噪声不确定性,估计的噪声功率可能与实际的噪声功率不同。令估计的噪声功率为 假定α(单位为dB)均匀地分布在区间[-B,B]中,其中B称为噪声不确定性因子。
(1)考虑2输入4接收机的***(M=2,P=4)。信道阶次为N1=N2=4(5抽头)。平滑因子选为L=5。
图7示出了MME方法和能量检测(有或无噪声不确定性)的检测概率。如果确切地知道噪声方差(B=0),则能量检测方法具有良好的性能。在理想噪声功率(大约3dB)的情况下,MME方法稍劣于能量检测方法。但是,噪声不确定性通常是存在的。如图7所示,若存在噪声不确定性(B>0),则能量检测比根据本发明实施例的方法差得多。
表4示出了误报概率(Pfa)。可以看出,根据本发明实施例的MME方法的Pfa满足Pfa≤0.1的要求。无噪声不确定性的能量检测的Pfa接近于要求的Pfa,但具有噪声不确定性的能量检测的Pfa远远超出了限定值。这意味着能量检测在实际的具有噪声不确定性的情况下是很不可靠的。
方法 | EG(2dB) | EG(1.5dB) | EG(1dB) | EG(0.5dB) | EG(无不确定性) | MME |
Pfa | 0.478 | 0.479 | 0.451 | 0.480 | 0.101 | 0.095 |
表4:误报概率(M=4,P=2)
(2)考虑单输入2接收机的***(M=2,P=1)。信道具有8抽头。平滑因子选为L=9。
图8示出了MME方法和能量检测(有或无噪声不确定性)的检测概率。表5示出了误报概率。从图8和表5中可以看出,与能量检测方法相比,根据本发明实施例的MME方法具有更好的性能。
方法 | EG(2dB) | EG(1.5dB) | EG(1dB) | EG(0.5dB) | EG(无不确定性) | MME |
Pfa | 0.479 | 0.474 | 0.441 | 0.473 | 0.094 | 0.098 |
表5:误报概率(M=2,P=1)
与相干检测方法相比,根据本发明实施例的基于特征值的检测不需要任何关于信号的信息,并且不需要同步。本发明的实施例也适用于多路传播。与能量检测方法相比,如以上仿真结果所示,本发明的实施例不受噪声不确定性影响,并且可获得较好的性能。
因此,提出本发明的实施例以在没有信号、信道、噪声等级和SNR的信息的情况下感测或检测信号。本发明的实施例可以应用于所有信号(例如,DTV、无线麦克等)。而且,由于阈值与信号和噪声功率无关,所以相同的阈值或比例因子可用于所有的信号。在检测概率和误报概率方面,本发明的实施例也表现出良好的性能。
虽然已参照特定的实施例具体示出和说明了本发明,但本领域技术人员应当理解,在不脱离所附的权利要求中限定的本发明的精神和范围的情况下,可作出各种形式上和细节上的变化。因此本发明的范围由所附权利要求书表明,并且落入权利要求书的等同物的意思和范围内的所有变化都应包含在内。
Claims (37)
1.一种用于对接收信号进行分类的方法,该方法包括:
确定所述接收信号的信号值的协方差矩阵;
确定所述协方差矩阵的特征值矩阵,所述特征值矩阵包括所述协方差矩阵的特征值;
由所述协方差矩阵的特征值中的至少一个特征值确定第一函数;
由所述协方差矩阵的特征值中的至少一个特征值确定第二函数,其中所述第二函数不同于所述第一函数;以及
根据所述第一函数的值和所述第二函数的值之间的比较结果将所述接收信号分类为包括数据的信号或噪声信号。
2.如权利要求1所述的方法,还包括:对所述接收信号进行采样。
3.如权利要求1所述的方法,还包括:对所述接收信号进行滤波。
4.如权利要求2所述的方法,还包括:对已采样的所述接收信号进行滤波。
5.如权利要求1所述的方法,其中,所述接收信号的信号值x(n)的协方差矩阵R(NS)根据以下等式确定:
其中NS表示所述协方差矩阵中信号值x(n)的个数。
6.如权利要求1所述的方法,其中,所述特征值矩阵通过将所述协方差矩阵变换为特征值矩阵来确定。
8.如权利要求1所述的方法,其中,所述第一函数是多个所述特征值的组合。
9.如权利要求8所述的方法,其中,所述第一函数是多个所述特征值的线性组合。
10.如权利要求1所述的方法,其中,所述第一函数至少包括所述特征值矩阵的最大特征值。
11.如权利要求10所述的方法,其中,所述第一函数包括所述特征值矩阵的预定个数的最大特征值。
12.如权利要求1所述的方法,其中,所述的确定所述第一函数的值包括确定所述特征值矩阵的特征值的迹。
13.如权利要求1所述的方法,其中,所述第二函数是多个所述特征值的组合。
14.如权利要求13所述的方法,其中,所述第二函数是多个所述特征值的线性组合。
15.如权利要求1所述的方法,其中,所述第二函数至少包括所述特征值矩阵的最小特征值。
16.如权利要求15所述的方法,其中,所述第二函数包括所述特征值矩阵的预定个数的最小特征值。
17.如权利要求1所述的方法,其中,在所述第一函数的值和所述第二函数的值之间的比较中使用比例因子。
18.一种对接收信号进行分类的分类器,该分类器包括:
用于确定所述接收信号的信号值的协方差矩阵的第一确定电路;
用于确定所述协方差矩阵的特征值矩阵的第二确定电路,所述特征值矩阵包括所述协方差矩阵的特征值;
用于由所述协方差矩阵的特征值中的至少一个特征值确定第一函数的第三确定电路;
用于由所述协方差矩阵的特征值中的至少一个特征值确定第二函数的第四确定电路,其中所述第二函数不同于所述第一函数;以及
用于根据所述第一函数的值和所述第二函数的值之间的比较结果将所述接收信号分类成包括数据的信号或噪声信号的分类电路。
19.如权利要求18所述的分类器,还包括对所述接收信号进行采样的采样器。
20.如权利要求18所述的分类器,还包括对所述接收信号进行滤波的滤波器。
21.如权利要求19所述的分类器,还包括对已采样的所述接收信号进行滤波的滤波器。
22.如权利要求18所述的分类器,其中,所述第一确定电路设置为根据以下等式确定所述接收信号的信号值x(n)的协方差矩阵R(NS):
其中NS表示所述协方差矩阵中信号值x(n)的个数。
23.如权利要求18所述的分类器,其中,所述第二确定电路设置为通过将所述协方差矩阵变换为特征值矩阵来确定所述特征值矩阵。
25.如权利要求18所述的分类器,其中,所述第一函数是多个所述特征值的组合。
26.如权利要求25所述的分类器,其中,所述第一函数是多个所述特征值的线性组合。
27.如权利要求18所述的分类器,其中,所述第一函数至少包括所述特征值矩阵的最大特征值。
28.如权利要求27所述的分类器,其中,所述第一函数包括所述特征值矩阵的预定个数的最大特征值。
29.如权利要求18所述的分类器,其中,所述第三确定电路设置为利用所述特征值矩阵的特征值的迹确定所述第一函数的值。
30.如权利要求18所述的分类器,其中,所述第二函数是多个所述特征值的组合。
31.如权利要求30所述的分类器,其中,所述第二函数是多个所述特征值的线性组合。
32.如权利要求18所述的分类器,其中,所述第二函数至少包括所述特征值矩阵的最小特征值。
33.如权利要求32所述的分类器,其中,所述第二函数包括所述特征值矩阵的预定个数的最小特征值。
34.如权利要求18所述的分类器,其中,所述分类器电路设置为在所述第一函数的值和所述第二函数的值之间的比较中使用比例因子。
35.一种电子装置,其包括权利要求18至34中任一项的分类器。
36.如权利要求35所述的电子装置,还包括用于接收信号的接收机。
37.如权利要求35所述的电子装置,其设置为选自以下一组电子装置:
-通信终端;
-设置为播放至少一种媒体的媒体播放装置。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US80658106P | 2006-07-05 | 2006-07-05 | |
US60/806,581 | 2006-07-05 | ||
US82583706P | 2006-09-15 | 2006-09-15 | |
US60/825,837 | 2006-09-15 | ||
PCT/SG2007/000198 WO2008004987A1 (en) | 2006-07-05 | 2007-07-05 | Method, device and computer program for classifying a received signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101507151A CN101507151A (zh) | 2009-08-12 |
CN101507151B true CN101507151B (zh) | 2013-01-09 |
Family
ID=38894855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007800304753A Active CN101507151B (zh) | 2006-07-05 | 2007-07-05 | 用于对接收信号进行分类的方法和装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8315297B2 (zh) |
CN (1) | CN101507151B (zh) |
SG (1) | SG173349A1 (zh) |
WO (1) | WO2008004987A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8379706B2 (en) * | 2009-07-28 | 2013-02-19 | Qualcomm Incorporated | Signal and noise power estimation |
US10321426B2 (en) * | 2009-10-20 | 2019-06-11 | Samsung Electronics Co., Ltd. | Apparatus and method for detecting synchronization and signals using block data processing |
CN101710848B (zh) * | 2009-11-13 | 2012-08-29 | 北京邮电大学 | 基于模糊积分和最优化理论的合作式频谱感知方法 |
US8386849B2 (en) * | 2010-01-29 | 2013-02-26 | Honeywell International Inc. | Noisy monitor detection and intermittent fault isolation |
CN102158254B (zh) * | 2011-02-16 | 2014-02-26 | 华为技术有限公司 | 信号接收方法及装置 |
CN102571237B (zh) * | 2011-12-31 | 2014-08-13 | 中国科学技术大学 | 一种基于周期图的无线信号检测方法 |
AT513834B1 (de) * | 2013-03-01 | 2014-08-15 | Cellstrom Gmbh | Elastomerer Endrahmen einer Redox-Durchflussbatterie |
US9823486B2 (en) * | 2013-03-11 | 2017-11-21 | Stc. Unm | Rotating point-spread function (PSF) design for three-dimensional imaging |
AU2014235359A1 (en) * | 2013-03-15 | 2015-10-08 | Nextnav, Llc | Methods and systems for improving time of arrival determination |
US9568591B2 (en) * | 2014-11-10 | 2017-02-14 | Peter Dan Morley | Method for search radar processing using random matrix theory |
CN105827559A (zh) * | 2015-01-08 | 2016-08-03 | 中兴通讯股份有限公司 | 天线合并算法的选择方法及装置 |
EP4000191A1 (en) * | 2019-07-17 | 2022-05-25 | Telefonaktiebolaget LM Ericsson (publ) | Distortion detection with multiple antennas |
CN112260779A (zh) * | 2020-09-25 | 2021-01-22 | 广东电网有限责任公司江门供电局 | 一种小型移动主用户的信号检测方法 |
CN114019449B (zh) * | 2022-01-10 | 2022-04-19 | 南京理工大学 | 信号源的波达方向估计方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6087974A (en) * | 1998-08-03 | 2000-07-11 | Lockheed Martin Corporation | Monopulse system for target location |
CN1449204A (zh) * | 2002-03-29 | 2003-10-15 | 上海贝尔有限公司 | 一种基于阵列天线移动通信***的信噪比测量方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
USH374H (en) * | 1987-02-09 | 1987-11-03 | The United States Of America As Represented By The Secretary Of The Army | Optimum multiple target detection and resolution |
US5333191A (en) * | 1993-04-22 | 1994-07-26 | Interdigital Technology Corporation | Detection of multifrequency tone signals |
JP2000131408A (ja) * | 1998-10-23 | 2000-05-12 | Mitsubishi Electric Corp | 方向探知装置 |
JP3920483B2 (ja) * | 1998-12-28 | 2007-05-30 | 株式会社東芝 | 電波到来方向推定方法およびアンテナ装置 |
JP2001359101A (ja) * | 2000-06-09 | 2001-12-26 | Toshiba Corp | フラッシュバック処理方法とこの方法を用いたデジタル映像信号処理装置及びデジタル放送受信装置 |
US6985434B2 (en) * | 2000-09-01 | 2006-01-10 | Nortel Networks Limited | Adaptive time diversity and spatial diversity for OFDM |
US20030204378A1 (en) * | 2001-06-26 | 2003-10-30 | Cai Khiem V. | System and method for forming a beam and creating nulls with an adaptive array antenna using orthogonal eigen-weighting |
US20040086161A1 (en) * | 2002-11-05 | 2004-05-06 | Radhika Sivaramakrishna | Automated detection of lung nodules from multi-slice CT image data |
JP4287670B2 (ja) * | 2003-02-18 | 2009-07-01 | パナソニック株式会社 | 通信装置及び通信方法 |
-
2007
- 2007-07-05 SG SG2011049210A patent/SG173349A1/en unknown
- 2007-07-05 CN CN2007800304753A patent/CN101507151B/zh active Active
- 2007-07-05 WO PCT/SG2007/000198 patent/WO2008004987A1/en active Application Filing
- 2007-07-05 US US12/307,071 patent/US8315297B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6087974A (en) * | 1998-08-03 | 2000-07-11 | Lockheed Martin Corporation | Monopulse system for target location |
CN1449204A (zh) * | 2002-03-29 | 2003-10-15 | 上海贝尔有限公司 | 一种基于阵列天线移动通信***的信噪比测量方法 |
Non-Patent Citations (2)
Title |
---|
BLAIR D.CARLSON.COVARIANCE MATRIX ESTIMATION ERRORS AND DIAGONAL LOADING IN ADAPTIVE ARRAYS.《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》.1988,第24卷(第4期),397-401. * |
Franz B.Tuteur.A NEW METHOD FOR SIGNAL DETECTION AND ESTIMATION USING THE EIGENSTRUCTURE OF THE COVARIANCE DIFFERENCE.《IEEE International Conference on ICASSP"86》.1986,2811-2814. * |
Also Published As
Publication number | Publication date |
---|---|
US20100015935A1 (en) | 2010-01-21 |
CN101507151A (zh) | 2009-08-12 |
US8315297B2 (en) | 2012-11-20 |
WO2008004987A1 (en) | 2008-01-10 |
SG173349A1 (en) | 2011-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101507151B (zh) | 用于对接收信号进行分类的方法和装置 | |
Zeng et al. | Spectrum-sensing algorithms for cognitive radio based on statistical covariances | |
US8259783B2 (en) | Method of determining as to whether a received signal includes an information signal | |
CN107820255B (zh) | 一种改进的协方差绝对值协作频谱感知方法 | |
CN103346845B (zh) | 基于快速傅里叶变换的盲频谱感知方法及装置 | |
US20100014616A1 (en) | Detection of Wideband Interference | |
US20150003508A1 (en) | Cooperative spectrum sensing apparatus and cooperative spectrum sensing method | |
CN101359930A (zh) | 认知无线电***中基于最大特征值的频谱感知方法 | |
CN108322277B (zh) | 一种基于协方差矩阵反特征值的频谱感知方法 | |
CN103141067A (zh) | 频带识别方法、装置和计算机程序产品以及性能评估方法、装置和计算机程序产品 | |
US9967117B2 (en) | Cooperative spectrum sensing system using sub-nyquist sampling and method thereof | |
CN104780006A (zh) | 基于最小错误概率准则的频谱检测器软融合方法 | |
CN110932806A (zh) | 一种alpha稳定噪声衰落信道下多天线频谱感知方法 | |
US20210409146A1 (en) | System and method for monitoring wireless communication channel by using cooperative jamming and spoofing | |
US8433009B2 (en) | Method for determining as to whether a received signal includes a data signal | |
CN108111213B (zh) | 一种用于多天线的频谱感知方法 | |
CN104965959A (zh) | 一种无线通信***抗电磁干扰性能评估方法 | |
CN104363064A (zh) | 一种基于优选用户的协作频谱感知方法 | |
Kortun et al. | Exact performance analysis of blindly combined energy detection for spectrum sensing | |
Smith et al. | Distribution and characteristic functions for correlated complex Wishart matrices | |
Hamid et al. | Spectrum sensing through spectrum discriminator and maximum minimum eigenvalue detector: A comparative study | |
CN105429913B (zh) | 基于特征值的多电平检测与识别方法 | |
CN103138856B (zh) | 一种检测干扰的方法及装置 | |
Ramírez et al. | Multiantenna spectrum sensing: Detection of spatial correlation among time-series with unknown spectra | |
CN109600181A (zh) | 一种用于多天线的频谱感知方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |