具体实施方式

在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解，“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。

在以下的描述中，所涉及的术语“第一\第二\第三等”或模块A、模块B、模块C等，仅用于区别类似的对象，不代表针对对象的特定排序，可以理解地，在允许的情况下可以互换特定的顺序或先后次序，以使这里描述的本申请实施例能够以除了在这里图示或描述的以外的顺序实施。

在以下的描述中，所涉及的表示步骤的标号，如S110、S120……等，并不表示一定会按此步骤执行，在允许的情况下可以互换前后步骤的顺序，或同时执行。

除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述本申请实施例的目的，不是旨在限制本申请。

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行描述。

为了精确定位探测目标的位置，除了需要知道目标到雷达阵列的距离r之外，还需要知道目标在雷达坐标系中的角度θ。假设目标距离雷达坐标系原点的距离为r。如图2a所示，对于配置了线阵毫米波雷达的系统来说，雷达坐标系的原点为线阵的中心位置。角度θ为目标和任意天线之间的连线与线阵法线之间的夹角。如图2b所示，目标在雷达坐标系中的位置可以通过距离r和角度θ确定。

在一种方案(记作为方案1)中，雷达可以根据接收的回波信号估计目标在雷达坐标系中的角度θ。如果雷达的接收天线位于信号源的远场，则接收的回波信号可以假设为平面波，探测目标视为点目标。具体地，假设每个天线单元间距为d_a，那么满足远场条件的回波信号到达相邻两个天线的路程差Δd为：

Δd＝d_asin(θ) (1)

该路程差导致的相位差为：

其中，λ为雷达毫米波的波长。

对于包含M个天线的线阵，由接收的回波信号组成的阵列响应y可以表示为：

如果空间有K个目标，每个目标与雷达的角度分别为θ₁，θ₂，…,θ_k，其中θ_k为K个目标中的第k个目标与天线的角度值。那么M个接收天线的接收信号组成的阵列响应y为：

其中，为第k个目标的回波在每个天线上的复增益。阵列响应y为K个目标的回波信号在M个接收天线上复合的阵列响应。

根据阵列响应y，采用机器学习的方法进行目标角度估计，可以采用的方法包括来波方向扫描(DOA)方法、多重信号分类(MUSIC)法、旋转不变参数估计技术(ESPRIT)等方法。以来波方向扫描(DOA)方法为例，雷达需要根据探测的角度范围构造多个导向矢量。比如，探测的角度范围为:-40度-40度，按照角度间隔为1度构造导向矢量，则导向矢量为：

其中，Θ_i为探测角度范围中的第i个角度。

将每个探测角度Θ_i对应的导向矢量与阵列响应y做相关，则第i个角度的相关值g(i)为：

相关值g(i)表示回波信号在探测角度Θ_i方向上的增益。将由公式(5)获得的每个探测角度对应的导向矢量与由公式(4)获得的阵列响应y点乘后求和，当第i个探测角度与第k个目标的回波信号的角度相同即i＝k时，根据公式(6)，同相合并，点乘积为1，此时相关值g(i)最大，阵列响应y在探测角度Θ_i方向上的增益最大。在目标数量或者信源数量K已知的条件下，则相关值g(i)中K个最大值的对应的每一个探测角度Θ_i即为K个信源的角度。

上述来波方向扫描(DOA)方法的原理是通过将每个探测角度对应的导向矢量与接收信号做相关，扫描到的最大相关值对应的探测角度即为目标在雷达坐标系中的角度。但是在一个距离上的具有相同速度的目标中，通过来波方向扫描(DOA)方法得到的一个最大相关值，对于以上两个角度很接近的目标无法精确区分目标的方位。例如在自动驾驶技术中，利用雷达探测前方的车辆的位置时，如果两辆车在前方并排行驶，一个车辆在雷达坐标系中的角度为5°，另一个车辆在雷达坐标系中的角度为5.5°，采用普通的检测算法如上述来波方向扫描(DOA)方法，很容易误认为是一辆车或者一辆大车在前方行驶，导致自动驾驶安全性下降。

为提高自动驾驶的安全性能，在雷达检测中，雷达的信号处理系统可以在进行准确的信源数量估计的基础上，采用某些特殊算法，比如角度超分辨算法，来检测处两辆车的位置。因此准确的信源数量估计往往是目标方位检测的基础。

在另一个方案中(记作方案2)，对信源数量估计采用信息论准则(AIC)或者最小描述长度(MDL)算法。在这些方法中需要用到阵列响应y的相关矩阵。

相关矩阵也叫相关系数矩阵，是由矩阵各列间的相关系数构成，相关矩阵第i行第j列的元素是原矩阵第i行第j列的相关系数。相关系数是反映变量之间相关关系密切程度的统计指标。

例如，如同方案1的实施场景中，雷达的接收天线位于信号源的远场，接收的回波信号为平面波，探测目标视为点目标，记为信源。对于包含M个天线的线阵，由接收的回波信号组成的阵列响应为y。对阵列响应y做N次快拍，所谓“快拍”是指对阵列响应在不同时刻进行采样。根据公式(4)获取M个接收天线的接收信号获得阵列响应y，则第j个快拍时刻采样获得的阵列响应为y^j：

其中，为第j个快拍时刻第i个天线接收的回波信号。阵列响应y^j为第j个快拍时刻M个天线接收的回波信号组成的向量。第j个快拍时刻的阵列响应y^j与其共轭转置(y^j)^H相乘获得第j个快拍时刻的M×M的自相关矩阵，将N个快拍的自相关矩阵相加后求平均值，获得相关矩阵R为：

在公式(8)中，R为M×M的相关矩阵，N为快拍的次数，(y^j)^H为第j个快拍时刻阵列响应y^j的共轭转置向量。

然后，对相关矩阵R做特征值分解:

R＝U∑U^H (9)

其中，U为相关相关矩阵R的分解矩阵，U^H为分解矩阵U的共轭转置矩阵；对角线矩阵Σ为：

对角线矩阵Σ的维度为M×M，对角线元素l_i为相关矩阵R的第i个特征值；其反映了理想的无噪声的回波信号数量的特征值。

提取对角线矩阵Σ的对角线元素l_i，将提取出的对角线元素按照从大到小的顺序排列，得到特征向量l＝{l₁,…,l_i,…,l_M'}：其中l₁≥l₂≥…≥l_M'。

根据相关矩阵R的特征向量l，可以以信息论准则(AIC)算法计算估计的信源数量Ns为：

其中，N_s＝0,1,2,…,M-1。则信源数量为：

信源数量为使AIC(N_s)达到最小值时的变量N_s的取值。

在该方案2中，估计信源数量算法的性能置信度取决于快拍的次数N的选取。当N取值较大的时候，该算法在高信噪比区间性能较好，置信度较高；但是低信噪比区间性能较差，置信度较低。当快拍的次数N值取较小值时候，该算法在低信噪比区间性能较好，置信度较高；但是在高信噪比区间存在估计偏差，置信度较低。即估计准确率不能逼近1。

由于方案1不能精确地分辨角度很接近的多个目标，而方案2进行信源数量估计时低信噪比区间和高信噪比区间性能存在偏差，定位精度均不理想，因此本申请提出采用多分类的支持向量机(SVM)估计信源的数量，下面介绍本申请的构思。

本申请实施例提供一种信源数量估计的方法对天线阵列的阵列响应进行处理得到相关矩阵，根据相关矩阵获得用于检测信源数量估计的特征向量，将特征向量输入支持向量机(SVM)，从而合理地估计出一个距离、速度单元上的多个信源的数量。

本申请实施例提供的信源数量估计的方法具体应用于雷达系统信号处理中的角度谱估计模块，根据支持向量机(SVM)获得信源数量，在后续的角度谱估计中可以采用某些特殊算法，比如角度超分辨算法，获得具体目标方位。后续的角度谱估计中如何根据信源数量进行角度和距离的计算处理，最终获得具体目标方位的算法本申请不做阐述。

下面介绍支持向量机(SVM)的原理。

支持向量机(SVM)是一类按照监督学习的方式对数据进行分类的分类器。假设分类器的输入向量为X＝{X₁,…,X_i,…,X_N}，所对应的分类为z＝{z₁,…,z_i,…,z_N}，其中第i个输入向量为X_i＝[x₁,x₂,…,x_n]。在二分类问题中，z_i∈{-1,1},z_i值为-1表示负类，z_i值为1表示正类。

在SVM的分类问题中，需要找到参数ω和b，使得：

z_i(ω^TX_i+b)≥1 (12)

其中ω为超平面的法向量，ω^T为超平面的法向量ω的转置向量；b为截距。

在某些环境中，X_i是线性不可分的，这时候往往需要将X_i映射到高维空间中，此时：

z_i(ω^Tφ(X_i)+b)≥1 (13)

其中，ω^Tφ(·)为核函数。

SVM算法最初是为二值分类问题设计的，当处理多类问题时，需要训练合适的多类分类器，常见的方法有一对多法(one-versus-rest,简称OVR SVMs)和一对一法(one-versus-one,简称OVO SVMs或者pairwise)。

在一对多法中，训练时依次把某个类别的样本归为一类，其他剩余的样本归为另一类，这样m个类别的样本就构造出了m个SVM。分类时将未知样本分类为具有最大分类函数值的那类。

例如假设有四个类别要划分(也就是4个标签)，分别是A、B、C和D。在抽取训练集的时候，分别抽取以下训练集：

(1)A所对应的向量作为正集，B，C，D所对应的向量作为负集；

(2)B所对应的向量作为正集，A，C，D所对应的向量作为负集；

(3)C所对应的向量作为正集，A，B，D所对应的向量作为负集；

(4)D所对应的向量作为正集，A，B，C所对应的向量作为负集；

使用这四个训练集分别进行训练，然后得到四个训练结果文件作为四个训练好的支持向量机。

在测试的时候，将测试向量分别利用这四个训练好的支持向量机进行测试，每个支持向量机输出一个评估值：f1(x)，f2(x)，f3(x)和f4(x)。该评估值用来评估输入的测试向量与相应的训练集的相似度，最终的分类结果是这四个评估值中最大的一个对应的类别。即Max(f1(x)，f2(x)，f3(x)和f4(x))。

下面结合图3-图11对本申请实施例提供的一种信源数量估计的方法和装置进行具体阐述。

图3为毫米波雷达接收系统的应用场景示意图。如图3所示，在毫米波雷达接收系统中，射频前端模块接收到回波信号Y，回波信号Y进入信号处理单元处理后获得估计的信源的数量，再经过数据处理单元进行处理得到角度和距离，最终获得具体信源位置。信号处理单元包括距离谱估计、速度谱估计、角度谱估计等模块，本申请的信源数量估计的方法和装置主要应用于角度谱估计模块。

要合理地估计信源数量，首先需要获得具有相同距离和相同速度的分辨目标上的回波信号Y在雷达天线阵列的阵列响应y；然后计算阵列响应y的相关矩阵R后，对相关矩阵R做特征值分解，从而提取相关矩阵R的特征向量；最后利用提取出的特征向量，通过支持向量机(SVM)估计的信源数量。

图4为本申请实施例提供的一种信源数量估计的方法的流程图。该方法的执行主体可以为具有计算处理能力的设备，服务器或者电子装置，如同方案1的实施场景，雷达阵列的接收天线位于信号源的远场，接收的回波信号为平面波，探测目标视为点目标，记为信源。下面针对每个步骤对本申请进行具体阐述。

如图4所示，首先执行步骤S401,获取雷达天线阵列的阵列响应。

具体地，该雷达阵列为M个等间距排列的天线单元组成的线阵，天线索引为1，…,M，以雷达的天线单元的法线方向为0度方向。该雷达阵列接收回波信号的示意图可以参考图5所示。假设在雷达阵列的探测前方位于同一个距离的区域内有N个具有相同速度的点目标，记为信源1、信源2、…,信源N，信源1的方位角为θ₁，信源2的方位角为θ₂，…,信源N的方位角为θ_n。由具有M个天线单元的雷达阵列接收到的阵列响应y为：

其中，y₁,y₂,…,y_M分别为M个天线单元的接收的回波信号，y_M为第M个天线的接收信号,M ≥1；阵列响应y是由M个天线单元接收的回波信号组成的行/列向量；s_n为信源n的复信号响应系数，其取值跟信源的材质有关；N_n为噪声向量，其取值由器件本身导致的热噪声决定；A(θ_n)为导向矢量，在M个天线单元上方位角为θ_n的回波信号的导向矢量A(θ_n)为：

当雷达阵列为面阵时，如果该面阵的排列为等间距的W行M列，W、M为任意自然数，则该雷达面阵包括W×M个天线单元，天线单元的索引为1，2，…,W×M。将该雷达面阵等效成W个行向量的拼接，每个行向量由M个天线单元的回波信号组成，由公式(14)获得每个线阵单元的阵列响应y，将W个阵列响应y拼接得到该雷达面阵阵列响应其中y_W×M为第W×M天线单元的接收信号。

当雷达阵列的排列为圆阵时，如果该圆阵包括M个沿圆周均匀排列的天线单元，M为任意自然数，可将该沿圆周均匀排列的M个天线单元拉伸为等间距的直线排列，由公式(14)获得则M个天线单元的雷达圆阵接收到的阵列响应y为

获取阵列响应y后，执行步骤S402，对阵列响应y计算，获得相关矩阵R。

在实际中考虑到接收信号是有限长度的，可以根据阵列响应y的N次快拍，采用时间-空间的评估方法计算相关矩阵R。具体地，对由公式(14)获取的阵列响应y做J次快拍采样，第j个快拍时刻的阵列响应y^j与其共轭转置(y^j)^H相乘获得第j个快拍时刻的M×M的自相关矩阵，将J个快拍的自相关矩阵相加后求平均值，参考公式(8)获得相关矩阵R为：

其中，y^j为第j个快拍时刻的阵列响应， 为第j个快拍时刻第i个天线接收的回波信号,为第j个快拍中第M个天线的响应，(y^j)^H为第j个快拍时刻的阵列响应y^j的共轭转置向量。

在上述步骤S402中，采用N次拍照计算出来的相关矩阵R维度较大，虽然提取的特征向量准确性高，但是需要对阵列响应y采样N次，增加了采样操作的工作量和数据量，计算较为复杂。

作为一种可选的实施方式，也可以对步骤S401采集到的阵列响应y做单次快拍，采用平滑分组的方式获取空间相关矩阵R。具体地，在步骤S402中，可以对阵列响应y做一次快拍采样；具体做法为，将阵列响应y的接收信号平滑分组，每M'个天线的阵列响应为一组，则M个阵列响应分为M-M'+1组向量,将采集到的阵列响应转化为向量y′：

例如，在一个雷达阵面为16个天线的线阵，取8个天线为一组，则将将采集到的阵列响应平滑分为由9个行向量组成的y′:

则阵列响应y的相关矩阵R为：

其中，s为索引值，M'为每组天线的数量，M'<M。

采用该方式计算阵列响应y的相关矩阵R，只需采样一次，获得的相关矩阵R的维度低，减少了计算量，通过平滑分组的方式评估空间相关矩阵R,可以精确地定向。该方式适用于根据雷达线阵和面阵的阵列响应y的相关矩阵R，但不适用于雷达圆阵。

得到相关矩阵R后，执行步骤S403，对相关矩阵R做特征值分解，提取相关矩阵R的特征向量。具体流程如图6所示，可以通过执行以下步骤S4031-4035实现。

S4031，对相关矩阵R做特征值分解，如公式(9):

R＝U∑U^H

其中，U为相关R的分解矩阵，U^H为U的共轭转置矩阵；Σ为理想的无噪声的阵列响应的特征值方阵，如公式(10)：

在特征值方阵Σ中对角线元素l₁ … l_i … l_M'为反映理想的无噪声的接收信号数量的特征值。

S4032，提取特征矩阵Σ的对角线元素，得到相关矩阵R的特征向量l＝[l₁ … l_i… l_M']^T，其中M'＝M。

为了减少支持向量机(SVM)的计算量，可以对特征向量l＝[l₁ … l_i … l_M']^T通过函数处理进行值域映射，执行步骤S4033-S4035。

S4033，以特征向量l＝[l₁ … l_i … l_M']^T为第一特征向量，对第一特征向量的每一个元素通过函数处理进行映射，得到特征函数x＝f(l)，其中，f(.)表示一种特征函数，特征函数可以为对数函数。

S4034，计算特征函数x＝f(l)，获得每一个特征值l_i对应的特征函数值x_i。

在一个可行的示例中，可以对第一特征向量的每一个元素通过对数函数处理进行映射，以每一个特征值l_i为对数函数的自变量，计算所述对数函数的值x_i，从而获得第二特征向量x＝[x₁ … x_n]^T。具体地，对数特征函数为x_i＝C*log(l_i)或者x_i＝C*log10(l_i)，其中C为常数。

采用该对数函数可以将特征参数从线性域转化为log域，减少最大值和最小值之间的差距，使特征函数的值域集中，离散性小，减少后续支持向量机(SVM)的计算量。

在另一个可行的示例中，可以先比较所述特征向量中的每一个特征值，获得最大特征值，然后计算每一个特征值与所述最大特征值的比值，以比值为对数函数的自变量，计算对数函数的值，获得每一个特征值对应的特征参数值x_i。

具体地，特征函数为x_i＝C*log(l_i/max(l))或者x_i＝C*log10(l_i/max(l))，其中C为常数。

通过这种对特征向量的归一化处理和对数函数映射，使得特征参数的分布域集中，离散性小，减少后续支持向量机(SVM)的计算量。

在第三个可行的示例中，可以对特征向量中的每一个特征值进行排序，以排序后的特征向量中的每一个特征值为对数函数的自变量，计算对数函数的值，获得每一个特征值对应的特征参数值x_i。

具体地，对特征向量l_i进行从大到小排序，获得特征向量l'＝sort(l₁ … l_n),则第i个特征参数可以为x_i＝C*log(l_i')；或者，x_i＝C*log10(l_i'),或者x_i＝C*log(l'_i/max(l))或者，x_i＝C*log10(l'_i/max(l))。

排序后采用该特征函数可以将特征参数从线性域转化为log域，通过对特征向量的归一化处理，可以进一步地减少最大值和最小值之间的差距，使特征函数的分布域集中，离散性小，利于计算。

S4035，以每一个特征函数的值x_i为元素组成特征向量x＝[x₁ … x_i … x_n]^T。其中x_i为第i个特征值l_i对应的特征参数值。

获得特征向量后，执行步骤S404，将特征向量输入训练好的与每一个分类项对应的支持向量机，输出对应信源数量的各个分类项的评估值，将评估值最大的分类项对应信源数量作为估计的信源数量。

特征向量可以为步骤S4032提取出的相关矩阵R的特征向量l＝[l₁… l_i…l_M']^T，优选为通过步骤S4035对数函数映射计算后获得的特征向量x＝[x₁ … x_i … x_n]^T。

最后，将特征向量输入支持向量机(SVM)，确定信源数量。具体如图7所示，包括以下步骤S4041-S4043。

S4041，预估计信源数量，对估计的信源数量进行分类，获得与信源数量对应数量的分类项。

具体地，预估计信源的数量为N_s，当预先估计信源数量N_s＝4时,在数量4以内，对信源数量可能出现的情况进行分类，可能是0个，1个、2个、3个或4个，获得的分类项s分别为s＝0,1,2,3,4。每一个分类项对应一个支持向量机。

S4042，将特征向量分别输入训练好的与每一个分类项对应的支持向量机，输出所述每一个分类项的估计值φ₀,φ₁,…,φ_n-1。

训练好的与每一个分类项对应的支持向量机可以以评估值φ₀,φ₁,…,φ_m对应的算法来表示。例如，信源数量为m的分类项对应的支持向量机输出的评估值φ_m的算法为：

φ_m＝H_m×K(D_m,x)+b_m (18)

其中，H_m为采用支持向量机训练出来的参数向量，维度为1×L。K(D_m,x)为支持向量机中的核函数，维度为L×1，L为支持向量机中支持向量的数量；D_m为采用支持向量机训练出来的支持向量；x为输入的特征向量。b_m为采用支持向量机训练出来的数值变量，维度为1×1。

则信源数量为m的分类项对应的支持向量机可以表示为：H_m×K(D_m,x)+b_m。

图8为利用训练好的支持向量机进行信源数量估计的示意图。如图8所示，将特征向量分别输入训练好的与每一个分类项对应的支持向量机，分别输出对应的信源数量的评估值φ₀,φ₁,…,φ_n-1。

支持向量机的核函数可以为线性核函数、多项式核函数、高斯核函数或sigmiod核函数之一。

其中线性核函数为：该线性核函数的函数值为以特征向量x为自变量，支持向量D_m与自变量的乘积。

多项式核函数为：该线性核函数的函数值为以第二特征向量x为自变量，支持向量D_m与自变量的乘积与1求和后的x次幂。

高斯核函数为：其中γ为配置参数，diag(.)为取对角线的操作。

sigmiod核函数为： 其中β和θ为配置参数，s为变量。

需要特别指出的是：对于线性核函数，信源数值的估值φ_m可以进一步表示为P_m为线性支持向量。在具体的实施过程中，D_m、H_m和b_m为通过支持向量机训练出来的参数。其下标m＝0,1,2,…,n-1。

最后执行步骤S4043，根据每一个分类项的估计值确定信源数量的值。

具体地，比较每一个分类项的估计值φ₀,φ₁,…,φ_n-1，获得最大的估计值，信源数量的取值N_s为φ₀,φ₁,…,φ_n-1中最大值所对应的下标，即

N_s＝argmax_m{φ₀,φ₁,…,φ_n-1} (19)。

本申请实施例提供的信源数量估计的方法，还包括训练与每一个分类项对应的支持向量机的步骤。在本申请的实施例中，估计信源数量属于多分类问题，因此在支持向量机的训练中通过一对多方式实现多分类问题。

具体地，假设有Q组训练向量，训练集为∪＝{d⁽¹⁾,d⁽²⁾,…,d^(Q)},其中d^(q)为M'×1的训练向量。Q组训练向量对应的信源数量分别为y⁽¹⁾,y⁽²⁾,…,y^(Q)，信源数量集为Y＝{y⁽¹⁾,y⁽²⁾,…,y^(Q)},其中y^(q)∈{0,1,…,M'-1}。

在本申请实施例中，要训练信源数量为m的参数D_m、H_m和b_m，可以将信源数量为m，y^(q)＝m的向量选择出来，并将其向量对应的真值r^(q)置位1，将y^(q)＝m以外的训练向量对应的真值r^(q)置位-1，即：

设定每一个支持向量机的核函数维度为L×1，L为支持向量机中支持向量的数量；将训练集∪＝{d⁽¹⁾,d⁽²⁾,…,d^(Q)}和真值集R＝{r⁽¹⁾,r⁽²⁾,…,r^(Q)}送入支持向量机进行训练，使得r^(q)(H_q×K(D_q,d^(q))+b_q)≥1。

通过训练获得L个支持向量其中支持向量矩阵为参数向量矩阵为其中为对应的真值，a₀,a₁,…a_L为训练的变量，b_q为训练出来的参数。

作为一种可选的实施方式，用于训练或者估计的输入向量可以为或者y的函数f(y)d＝y或者d＝f(y)，特别的，其余过程与上述实施用例类似。

在仿真环境下，分别采用本申请提供的信源数量估计的方法和改进的MDL方法进行信源数量估计的置信度对比。在相同条件下，雷达阵列为16个天线的线阵，以波长为天线间距，信源数量分别为：1、2和3个，其中2个信源以3度角间隔，3个7度角间隔排列。

图9a为使用本申请提供的信源数量估计的方法进行信源数量估计结果的置信度曲线图。如图9a所示，图中横坐标为信噪比(SNR)，纵坐标为估计结果的准确率，坐标系内分别示出了一个信源的估计结果置信度曲线、两个信源的估计结果置信度曲线和三个信源的估计结果置信度曲线。对置信度曲线分析可得，在信噪比为10的环境下，一个信源的估计结果准确率为99％，两个信源的估计结果的准确率为95％，三个信源的估计结果的准确率为45％。随着信噪比增加，信源数量的估计结果的准确率增加，在信噪比为20的环境下，一个和两个信源的估计结果准确率为100％，三个信源的估计结果的准确率接近为90％。

图9b为使用改进的MDL进行信源数量估计的置信度曲线图。如图9b所示，在信噪比为10的环境下，两个信源的估计结果的准确率为80％，三个信源的估计结果的准确率只有28％；在信噪比为20的环境下，一个和两个信源的估计结果准确率为100％，三个信源的估计结果的准确率接近为80％。

通过图9a和图9b对比分析可得，使用本申请提供的一种信源数量估计的方法进行信源数量估计结果具有较高的置信度，能够解决传统信源数量估计方法中存在的低信噪比区间性能和高信噪比区间性能存在偏差较大的问题。

本申请实施例提供的信源数量估计的方法，通过大规模数据训练，无论是在高信噪比还是低信噪比下都具有良好性能。在支持向量机采用线性核函数时，本方法的计算复杂度很低。

本申请实施例提供的信源数量估计的方法不但可以根据同一个距离、速度单元上天线阵列响应估计信源数量，还可以估计其他维度的信源数量，比如同一个距离、方位角度方向上速度维度的信源数量，或者同一个距离单元上，速度维度和角度维度上总共的信源数量等等。

本申请实施例提供一种信源数量估计的装置，如图10所示，所述装置包括：数据采集模块1001，用于获取雷达天线阵列的阵列响应，所述雷达天线阵列包括至少一个天线单元；相关矩阵计算模块1002，用于对所述阵列响应进行自相关计算得到相关矩阵；特征提取模块1003，用于对所述相关矩阵做特征值分解，提取特征向量；评估值计算模块1004，用于设定预估计的信源数量，对所述预估计的信源数量分类，获得至少一个分类项；和信源数量确定模块1005，用于将所述特征向量输入与所述至少一个分类项对应的支持向量机，输出对应所述至少一个分类项的评估值；比较所述至少一个分类项的评估值，获得其中最大的评估值；根据所述最大的评估值对应的分类项确定信源数量。

本申请实施例提供的信源数量估计的装置，还包括训练支持向量机模块1006，所述模块根据不同的信源数量将训练集中的训练向量分组；所述训练向量为带有信源数量标签的特征向量；根据所述分组选择所述每一个分类项对应组别的训练向量以所述对应组别的训练向量为正；以训练集中的所述对应组别以外的训练向量为负，形成真值集；将所述训练集和真值集输入与每一个分类项对应的所述支持向量机，训练所述支持向量机的支持向量、参数向量和变量,获得训练好的与每一个分类项对应的支持向量机。

本申请实施例提供一种电子装置1100，包括处理器1101和存储器1102；所述处理器1101用于执行所述存储器1102所存储的计算机执行指令，所述处理器1101运行所述计算机执行指令执行上述任意实施例所述的信源数量估计的方法。

本申请实施例提供一种存储介质1103，包括可读存储介质和存储在所述可读存储介质中的计算机程序，所述计算机程序用于实现上述任意一实施例所述的信源数量估计的方法。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请实施例的范围。

此外，本申请实施例的各个方面或特征可以实现成方法、装置或使用标准编程和/或工程技术的制品。本申请中使用的术语“制品”涵盖可从任何计算机可读器件、载体或介质访问的计算机程序。例如，计算机可读介质可以包括，但不限于:磁存储器件(例如，硬盘、软盘或磁带等)，光盘(例如，压缩盘(compact disc，CD)、数字通用盘(digital versatiledisc，DVD)等)，智能卡和闪存器件(例如，可擦写可编程只读存储器(erasableprogrammable read-only memory，EPROM)、卡、棒或钥匙驱动器等)。另外，本文描述的各种存储介质可代表用于存储信息的一个或多个设备和/或其它机器可读介质。术语“机器可读介质”可包括但不限于，无线信道和能够存储、包含和/或承载指令和/或数据的各种其它介质。应当理解的是，在本申请实施例的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者接入网设备等)执行本申请实施例各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅为本申请实施例的具体实施方式，但本申请实施例的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请实施例揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请实施例的保护范围之内。