一种波达方向估计方法和装置
阅读：643发布：2020-05-11
专利汇可以提供一种波达方向估计方法和装置专利检索，专利查询，专利分析的服务。并且本 发明 实施例 提供了一种 波达方向 估计方法和装置。 电子 设备分别获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的 信号 源发出的信号， 叠加 整合获取的各个信号，得到两个L型天线阵列接收的整体信号，确定该整体信号的协方差矩阵，对该协方差矩阵进行特征值分解，得到整体噪声子空间；提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵，对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间。根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。通过本发明实施例提供的技术方案，充分利用了x轴接收信号的协方差信息，提高了DOA估计的精准度，更满足实际应用场景的需要。，下面是一种波达方向估计方法和装置专利的具体信息内容。
1.一种波达方向DOA估计方法，其特征在于，应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素；所述方法包括：
获取所述第一x轴天线子阵列、所述第一y轴天线子阵列、所述第二x轴天线子阵列和所述第二y轴天线子阵列分别接收的信号源发出的信号；
叠加整合获取的各个信号，得到所述两个L型天线阵列接收的整体信号；
确定所述两个L型天线阵列接收的整体信号的协方差矩阵；
对所述协方差矩阵进行特征值分解，得到整体噪声子空间；
提取所述协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵；
对所述x轴协方差矩阵进行特征值分解，得到x轴噪声子空间；
根据所述x轴噪声子空间和所述整体噪声子空间，估计接收信号的DOA。
2.根据权利要求1所述的方法，其特征在于，所述叠加整合获取的各个信号，得到所述两个L型天线阵列接收的整体信号的步骤，包括：
叠加整合所述第一x轴天线子阵列接收的信号和所述第二x轴天线子阵列接收的信号，得到x轴接收信号；
叠加整合所述第一y轴天线子阵列接收的信号和所述第二y轴天线子阵列接收的信号，得到y轴接收信号；
叠加整合所述x轴接收信号和所述y轴接收信号，得到所述两个L型天线阵列接收的整体信号。
3.根据权利要求1所述的方法，其特征在于，所述确定所述两个L型天线阵列接收的整体信号的协方差矩阵，包括：
利用以下公式，确定所述两个L型天线阵列接收的整体信号的协方差矩阵
其中，F为对所述整体信号的采样值个数，z(t)为所述两个L型天线阵列接收的整体信号的第t个采样值。
4.根据权利要求1所述的方法，其特征在于，所述根据所述x轴噪声子空间和所述整体噪声子空间，估计接收信号的DOA的步骤，包括：
根据所述x轴噪声子空间，构造第一谱函数；
对所述第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
根据所述每一信号的第一变量和所述整体噪声子空间，构造第二谱函数；
对所述第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
根据所述第一变量和所述第二变量，估计接收信号的DOA。
5.根据权利要求4所述的方法，其特征在于，所述根据所述x轴噪声子空间，构造第一谱函数的步骤，包括：
利用以下公式，构造第一谱函数f(μ)：
其中， 为所述x轴噪声子空间，ax(μk)是x轴接收信号的导向矢量，μk为第k个信号的第一变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
6.根据权利要求4所述的方法，其特征在于，所述根据所述每一信号的第一变量和所述整体噪声子空间，构造第二谱函数的步骤，包括：
利用以下公式，构造第二谱函数f(ν)：
其中， 为所述整体噪声子空间， 为第k个信号的第一变量μk的估计值 的导向矢量，ay(νk)是y轴接收信号的导向矢量，νk为第k个信号的第二变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
7.根据权利要求4所述的方法，其特征在于，所述根据所述第一变量和所述第二变量，估计接收信号的DOA的步骤，包括：
利用以下公式，估计接收到的第k个信号与z轴的夹角θk和第k个信号与xOy平面的夹角其中， 为所述第k个信号的第一变量μk的估计值， 为所述第k个信号的第二变量νk的估计值， 为所述第k个信号与z轴的夹角θk的估计值， 为所述第k个信号与xOy平面的夹角 的估计值，1≤k≤K，K表示接收到所述信号源发出的信号的总个数；
将所述估计值 与 作为接收信号的DOA。
8.一种波达方向DOA估计装置，其特征在于，应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素；所述装置包括：
信号接收单元，用于获取所述第一x轴天线子阵列、所述第一y轴天线子阵列、所述第二x轴天线子阵列和所述第二y轴天线子阵列分别接收的信号源发出的信号；
信号叠加单元，用于叠加整合获取的各个信号，得到所述两个L型天线阵列接收的整体信号；
矩阵计算单元，用于确定所述两个L型天线阵列接收的整体信号的协方差矩阵；
第一特征值分解单元，用于对所述协方差矩阵进行特征值分解，得到整体噪声子空间；
矩阵提取单元，用于提取所述协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵；
第二特征值分解单元，用于对所述x轴协方差矩阵进行特征值分解，得到x轴噪声子空间；
DOA估计单元，用于根据所述x轴噪声子空间和所述整体噪声子空间，估计接收信号的DOA。
9.根据权利要求8所述的装置，其特征在于，所述DOA估计单元，包括：
第一谱函数构造子单元，用于根据所述x轴噪声子空间，构造第一谱函数；
第一峰值搜索子单元，用于对所述第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
第二谱函数构造子单元，用于根据所述每一信号第一变量和所述整体噪声子空间，构造第二谱函数；
第二峰值搜索子单元，用于对所述第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
DOA估算子单元，用于根据所述第一变量和所述第二变量，估计接收信号的DOA。
10.一种电子设备，其特征在于，包括处理器和机器可读存储介质，所述机器可读存储介质存储有能够被所述处理器执行的机器可执行指令，所述处理器被所述机器可执行指令促使：实现权利要求1-7任一所述的方法步骤。
一种波达方向估计方法和装置
技术领域
背景技术
[0002] DOA(Direction of Arrival，波达方向)定位技术，是指通过处理接收到的信号，获取信号的到达方向，进而估计信号源的距离信息和方位信息。目前，DOA估计算法主要包括：针对均匀线性天线阵列的多重信号分类(MUSIC)算法和针对均匀矩形天线阵列的传统二维MUSIC算法。针对位于x轴、y轴的传统互素L型天线阵列，可以沿x轴、y轴进行扩展，得到展开的互素L型天线阵列，使用降低复杂度的MUSIC(RC-MUSIC)算法代替传统二维MUSIC算法进行DOA估计，大大降低运算复杂度。
[0003] 然而，RC-MUSIC算法在第一步使用MUSIC算法进行DOA估计时，需要使用天线阵列在x轴上分布的阵元所接收的信号(以下简称x轴接收信号)的噪声子空间。其使用的方法是选取阵列接收信号总的噪声子空间的一半作为x轴噪声子空间(以下简称x轴噪声子空间)。但这样并未完全利用阵列在x轴接收信号的协方差信息，故获取的噪声子空间不够准确，所以用该噪声子空间进行DOA估计时会带来较大误差。因此，目前RC-MUSIC算法的估计精度不能满足于实际应用，有待提高。
发明内容
[0004] 本发明实施例的目的在于提供一种波达方向估计方法和装置，以实现充分利用x轴接收信号的协方差信息，提高DOA估计的精准度，更满足实际应用场景的需要。具体技术方案如下：
[0005] 第一方面，本发明实施例提供了一种波达方向估计方法，应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素；包括：
[0006] 获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号；
[0007] 叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号；
[0008] 确定两个L型天线阵列接收的整体信号的协方差矩阵；
[0009] 对该协方差矩阵进行特征值分解，得到整体噪声子空间；
[0010] 提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵；
[0011] 对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间；
[0012] 根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0013] 可选的，叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号的步骤，包括：
[0014] 叠加整合第一x轴天线子阵列接收的信号和第二x轴天线子阵列接收的信号，得到x轴接收信号；
[0015] 叠加整合第一y轴天线子阵列接收的信号和第二y轴天线子阵列接收的信号，得到y轴接收信号；
[0016] 叠加整合x轴接收信号和y轴接收信号，得到两个L型天线阵列接收的整体信号。
[0017] 可选的，确定两个L型天线阵列接收的整体信号的协方差矩阵，包括：
[0018] 利用以下公式，确定两个L型天线阵列接收的整体信号的协方差矩阵[0019]
[0020] 其中，F为对所述整体信号的采样值个数，z(t)为所述两个L型天线阵列接收的整体信号的第t个采样值。
[0021] 可选的，根据x轴噪声子空间和所述整体噪声子空间，估计接收信号的DOA的步骤，包括：
[0022] 根据x轴噪声子空间，构造第一谱函数；
[0023] 对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0024] 根据每一信号的第一变量和整体噪声子空间，构造第二谱函数；
[0025] 对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0026] 根据第一变量和第二变量，估计接收信号的DOA。
[0027] 可选的，根据x轴噪声子空间，构造第一谱函数的步骤，包括：
[0028] 利用以下公式，构造第一谱函数f(μ)：
[0029]
[0030] 其中， 为x轴噪声子空间，ax(μk)是x轴接收信号的导向矢量，μk为第k个信号的第一变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0031] 可选的，根据每一信号的第一变量和整体噪声子空间，构造第二谱函数的步骤，包括：
[0032] 利用以下公式，构造第二谱函数f(ν)：
[0033]
[0034] 其中， 为整体噪声子空间， 为第k个信号的第一变量μk的估计值 的导向矢量，ay(νk)是y轴接收信号的导向矢量，νk为第k个信号的第二变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0035] 可选的，根据第一变量和第二变量，估计接收信号的DOA的步骤，包括：
[0036] 利用以下公式，估计接收到的第k个信号与z轴的夹角θk和该信号与xOy平面的夹角
[0037]
[0038]
[0039] 其中， 为第k个信号的第一变量μk的估计值， 为第k个信号的第二变量νk的估计值， 为第k个信号与z轴的夹角θk的估计值， 为第k个信号与xOy平面的夹角 的估计值，1≤k≤K，K表示接收到信号源发出的信号的总个数；
[0040] 将该估计值 与 作为接收信号的DOA。
[0041] 第二方面，本发明实施例提供了一种波达方向估计装置，应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素；包括：
[0042] 信号接收单元，用于获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号；
[0043] 信号叠加单元，用于叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号；
[0044] 矩阵计算单元，用于确定两个L型天线阵列接收的整体信号的协方差矩阵；
[0045] 第一特征值分解单元，用于对该协方差矩阵进行特征值分解，得到整体噪声子空间；
[0046] 矩阵提取单元，用于提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵；
[0047] 第二特征值分解单元，用于对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间；
[0048] DOA估计单元，用于根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0049] 可选的，信号叠加单元具体用于：
[0050] 叠加整合第一x轴天线子阵列接收的信号和第二x轴天线子阵列接收的信号，得到x轴接收信号；
[0051] 叠加整合第一y轴天线子阵列接收的信号和第二y轴天线子阵列接收的信号，得到y轴接收信号；
[0052] 叠加整合x轴接收信号和所述y轴接收信号，得到两个L型天线阵列接收的整体信号。
[0053] 可选的，矩阵计算单元具体用于：
[0054] 利用以下公式，确定两个L型天线阵列接收的整体信号的协方差矩阵[0055]
[0056] 其中，t＝1,2,…,F，F为对所述整体信号的采样值个数，z(t)为两个L型天线阵列接收的整体信号，zH(t)为z(t)的共轭转置。
[0057] 可选的，DOA估计单元具体用于：
[0058] 根据x轴噪声子空间，构造第一谱函数；
[0059] 对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0060] 根据每一信号的第一变量和整体噪声子空间，构造第二谱函数；
[0061] 对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0062] 根据第一变量和第二变量，估计接收信号的DOA。
[0063] 可选的，DOA估计单元，包括：
[0064] 第一谱函数构造子单元，用于根据x轴噪声子空间，构造第一谱函数；
[0065] 第一峰值搜索子单元，用于对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0066] 第二谱函数构造子单元，用于根据每一信号第一变量和所述整体噪声子空间，构造第二谱函数；
[0067] 第二峰值搜索子单元，用于对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0068] DOA估算子单元，用于根据第一变量和第二变量，估计接收信号的DOA。
[0069] 可选的，第一谱函数构造子单元具体用于：
[0070] 利用以下公式，构造第一谱函数f(μ)：
[0071]
[0072] 其中， 为所述x轴噪声子空间，ax(μk)是x轴接收信号的导向矢量，μk为第k个信号的第一变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0073] 可选的，DOA估算子单元具体用于：
[0074] 利用以下公式，估计接收到的第k个信号与z轴的夹角θk和与xOy平面的夹角[0075]
[0076]
[0077] 其中， 为第k个信号的第一变量μk的估计值， 为第k个信号的第二变量νk的估计值， 为第k个信号与z轴的夹角θk的估计值， 为第k个信号与xOy平面的夹角 的估计值，1≤k≤K，K表示接收到信号源发出的信号的总个数；
[0078] 将该估计值 与 作为接收信号的DOA。
[0079] 第三方面，本发明实施例提供了一种电子设备，包括处理器和机器可读存储介质，机器可读存储介质存储有能够被所述处理器执行的机器可执行指令，处理器被机器可执行指令促使：实现上述第一方面提供的波达方向估计方法的任一步骤。
[0080] 第四方面，本发明实施例提供了一种机器可读存储介质，存储有能够被处理器执行的机器可执行指令，处理器被所述机器可执行指令促使：实现上述第一方面提供的波达方向估计方法的任一步骤。
[0081] 由上述技术方案可见，本发明实施例提供了一种DOA估计方法和装置，应用于与两个L型天线阵列连接的电子设备，两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。电子设备可以获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号，叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号，确定该整体信号的协方差矩阵，对该协方差矩阵进行特征值分解，得到整体噪声子空间；提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵，对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间。根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0082] 在本发明实施例中，两个L型天线阵列在估计接收信号的波达方向DOA时，不再选取接收信号整个噪声子空间的一半作为x轴噪声子空间，而是提取x轴接收信号的协方差矩阵进行特征值分解，得到x轴噪声子空间，而这样就充分利用了x轴接收信号的协方差信息，提高了DOA估计的精准度，更满足实际应用场景的需要。
[0084] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0085] 图1为本发明实施例提供的DOA估计方法的一种流程示意图；
[0086] 图2为两个L型天线阵列的结构示意图；
[0087] 图3为本发明实施例提供的DOA估计方法的另一种流程示意图；
[0088] 图4为本发明实施例提供的DOA估计装置的一种结构示意图；
[0089] 图5为本发明实施例提供的DOA估计装置的另一种结构示意图；
[0090] 图6为本发明实施例提供的电子设备的一种结构示意图。
具体实施方式
[0091] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0092] 下面通过具体实施例，对本发明进行详细说明。
[0093] 参见图1，图1为本发明实施例提供的DOA估计方法的一种流程图。该方法应用于与两个L型天线阵列连接的电子设备，两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。该方法包括如下步骤。
[0094] 步骤101，获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号。
[0095] 在一个实施例中，第一x轴天线子阵列和第一y轴天线子阵列是互相垂直的均匀直线天线阵列，且阵元彼此间的距离为N2λ/2；第二x轴天线子阵列和第二y轴天线子阵列也是互相垂直的均匀直线天线阵列，且阵元彼此间的距离为N1λ/2。其中N1和N2是两个互素的整数，λ为载波的波长。如图2所示，图2为两个L型天线阵列的结构示意图。其中，第一x轴天线子阵列位于x轴正半轴，第一y轴天线子阵列位于y轴正半轴，第二x轴天线子阵列位于x轴负半轴，第二y轴天线子阵列位于y轴负半轴。两个L型天线阵列接收到信号源发出的信号为K个不相关的、波长为λ的窄带远场信号sk(t)(k＝1,2,...,K)，分别用θk和 表示第k个入射信号与z轴的夹角(俯仰角)和与xOy平面的夹角(方位角)，它们的范围是θk∈(0,90°)以及[0096] 在另一个实施例中，获取的第一x轴天线子阵列、第二x轴天线子阵列、第一y轴天线子阵列和第二y轴天线子阵列接收的信号分别表示为
[0097]
[0098]
[0099]
[0100]
[0101] 其中，t＝1,2,…,F，F为对每个接收信号的采样值个数；s(t)为信号矢量，s(t)＝[s1(t),s2(t),...,sK(t)]T；nx1(t)、nx2(t)、ny1(t)和ny2(t)都是均值为0且方差为 的加性高斯白噪声向量； 和 分别表示第一x轴天线子阵列、第二x轴天线子阵列、第一y轴天线子阵列和第二y轴天线子阵列的流型矩阵； 表示复数域。另外，上述流型矩阵也可以表示为：
[0102]
[0103]
[0104]
[0105]
[0106] 其中 和 分别为第一x轴天线子阵列、第二x轴天线子阵列、第一y轴天线子阵列和第二y轴天线子阵列的导向矢量，该导向矢量中包含接收信号的波达方向θk和 上述导向矢量又可以分别表示为：
[0107]
[0108]
[0109]
[0110]
[0111] 其中di表示天线阵列中阵元之间的距离。
[0112] 步骤102，叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号。
[0113] 在本步骤中，电子设备在获取到上述四个天线子阵列接收的信号后，对获取的信号进行叠加整合，得到该天线阵列接收的整体信号。
[0114] 在一个实施例中，电子设备把x轴、y轴接收到的信号分别进行叠加整合，得到：
[0115]
[0116]
[0117] 其中，Ax、Ay分别表示该天线阵列在x轴、y轴的流型矩阵，记为：
[0118]
[0119]
[0120] 其中， 和 分别为该天线阵列在x轴、y轴的导向矢量，记为：
[0121]
[0122]
[0123] 电子设备再把x轴接收信号与y轴接收信号进行叠加整合，得到两个L型天线阵列接收的整体信号：
[0124]
[0125] 其中，整体信号的流型矩阵为导向矢量为
[0126] 步骤103，确定两个L型天线阵列接收的整体信号的协方差矩阵
[0127] 在本步骤中，电子设备对前述叠加整合得到的整体信号z(t)求协方差矩阵。
[0128] 一个可选的实施例中，整体信号的协方差矩阵可以表示为：
[0129]
[0130] 其中，Rs＝E[s(t)sH(t)]是信源信号的协方差矩阵，I为单位矩阵。
[0131] 在另一个可选的实施例中，可以用整体信号的F个采样值的时间平均来近似代替整体信号的协方差矩阵。具体的，整体信号的协方差矩阵可以近似表示为：
[0132]
[0133] 其中，F为对所述整体信号的采样值个数，z(t)为所述两个L型天线阵列接收的整体信号的第t个采样值，zH(t)为z(t)的共轭转置。
[0134] 步骤104，对该协方差矩阵进行特征值分解，得到整体噪声子空间。
[0135] 在本步骤中，对该协方差矩阵 进行特征值分解，得到整体噪声子空间 具体的，使用以下公式对 进行特征值分解：
[0136]
[0137] 其中， 表示整体噪声子空间； 表示整体信号子空间； 是主对角线为信号功率的对角矩阵；
是主对角线为噪声功率的对角矩阵。
是主对角线为噪声功率的对角矩阵。
[0138] 步骤105，提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵。
[0139] 在本步骤中，提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵，记为
[0140] 在本实施例中，对步骤105与步骤104的执行顺序不作具体限定。步骤105可以在步骤104之后执行，也可以在步骤104之前执行，还可以与步骤104同时执行。
[0141] 步骤106，对x轴协方差矩阵 进行特征值分解，得到x轴噪声子空间。
[0142] 在本步骤中，对x轴接收信号的协方差矩阵 进行特征值分解：
[0143]
[0144] 其中， 为x轴噪声子空间， 为x轴信号子空间， 为该x轴协方差矩阵的特征值， 为该x轴协方差矩阵的信号特征值， 为该x轴协方差矩阵的噪声特征值。
[0145] 步骤107，根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0146] 参见图3，图3为本发明实施例提供的DOA估计方法的另一种流程示意图。该方法应用于与两个L型天线阵列连接的电子设备，两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。该方法包括如下步骤。
[0147] 步骤301，获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号。步骤301与步骤101相同。
[0148] 步骤302，叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号。步骤302与步骤102相同。
[0149] 步骤303，确定两个L型天线阵列接收的整体信号的协方差矩阵。步骤303与步骤103相同。
[0150] 步骤304，对该协方差矩阵进行特征值分解，得到整体噪声子空间。步骤304与步骤104相同。
[0151] 步骤305，提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵。步骤305与步骤105相同。
[0152] 在本实施例中，对步骤305与步骤304的执行顺序不作具体限定。步骤305可以在步骤304之后执行，也可以在步骤304之前执行，还可以与步骤304同时执行。
[0153] 步骤306，对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间。步骤306与步骤106相同。
[0154] 步骤307，根据x轴噪声子空间，构造第一谱函数。
[0155] 在一个实施例中，定义第一转换方程 结合θk和 的取值范围θk∈(0,90°)和 可得uk∈(-1,1)。将μk记为第k个信号的第一变量，则x轴接收信号的导向矢量可表示为ax(μk)，记为：
[0156]
[0157] 其中， 由于x轴接收信号的导向矢量ax(μk)与信号子空间相同，又由于信号子空间与噪声子空间正交，所以当第k个信号的第一变量μk满足波达方向时，x轴接收信号的导向矢量ax(μk)与噪声子空间 正交。因此，利用以下公式，构造第一变量的第一谱函数f(μ)：
[0158]
[0159] 其中， 为x轴噪声子空间，ax(μk)是x轴接收信号的导向矢量，μk为第k个信号的第一变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0160] 步骤308，对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量。
[0161] 在本步骤中，对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；其中，由步骤307可知，每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量。
[0162] 在一个实施例中，电子设备可以按照从小到大的顺序以固定步进搜索第k个信号的第一变量μk的估计值。具体的，电子设备按照从小到大的顺序以固定步进搜索第k个信号的第一变量μk的可能取值，当第k个信号的第一变量μk的一个可能取值满足波达方向时，x轴接收信号的导向矢量ax(μk)会与x轴噪声子空间 正交，即 此时第一谱函数的分母为0，谱函数输出峰值，电子设备记录下该μk的值，再继续往后进行搜索。依此类推，共能得到K个估计值，记为
[0163] 在另一个实施例中，电子设备可以按照从大到小的顺序以固定步进搜索第k个信号的第一变量μk的可能取值，来确定第k个信号的第一变量的估计值。本发明实施例对此不作具体限定。
[0164] 步骤309，根据每一信号的第一变量和整体噪声子空间，构造第二谱函数。
[0165] 在一个实施例中，定义第二转换方程 结合θk和 的取值范围θk∈(0,90°)和 可得vk∈(0,1)。将vk记为第k个信号的第二变量，则y轴接收信号的导向矢量可表示为ay(νk)，记为：
[0166]
[0167] 其中， 由于两个L型天线阵列接收的整体信号的导向矢量与信号子空间相同，又由于信号子空间与噪声子空间正交且第k个信号的第一变量在步骤308中得出其估计值 所以当第k个信号的第二变量νk满足波达方向时，整体信号的导向矢量 与整体噪声子空间 正交。因此，利用以下公式，构造第二变量的第二谱函数f(ν)：
[0168]
[0169] 其中， 为整体噪声子空间， 为第k个信号的第一变量μk的估计值 的导向矢量，ay(νk)是y轴接收信号的导向矢量，νk为第k个信号的第二变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0170] 步骤310，对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量。
[0171] 在本步骤中，对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；其中，由步骤309可知，每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量。
[0172] 在一个实施例中，电子设备可以按照从小到大的顺序以固定步进搜索第k个信号的第二变量νk的可能取值，来确定第k个信号的第二变量的估计值。具体的，电子设备按照从小到大的顺序以固定步进搜索第k个信号的第二变量νk的可能取值，当第k个信号的第二变量νk的一个可能取值满足波达方向时，整体信号的导向矢量 会与整体噪声子空间 正交，即 此时第二谱函数的分母为0，谱函数输出峰值，电子设备记录下该νk的值，再继续往后进行搜索。依此类推，共能得到K个估计值，记为[0173] 在另一个实施例中，电子设备可以按照从大到小的顺序以固定步进搜索第k个信号的第二变量νk的可能取值，来确定第k个信号的第二变量的估计值。本发明实施例对此不作具体限定。
[0174] 步骤311，根据第一变量和第二变量，估计接收信号的DOA。
[0175] 在本步骤中，电子设备根据第k个信号的第一变量的估计值 和第k个信号的第二变量的估计值 估计出接收到的第k个信号与z轴的夹角 和该信号与xOy平面的夹角将估计值 和 即作为接收信号的DOA。
[0176] 在一个实施例中，因为在步骤310中，有估计值 作为角度搜索的约束条件，所以对上述第二谱函数进行峰值搜索时， 会和 自动配对。所以信号z轴的夹角(俯仰角)的估计值 和与xOy平面的夹角(方位角)的估计值 也能够自动配对。利用以下公式，估计接收到的第k个信号与z轴的夹角θk和该信号与xOy平面的夹角
[0177]
[0178]
[0179] 其中， 为所述第k个信号的第一变量μk的估计值， 为所述第k个信号的第二变量νk的估计值， 为所述第k个信号与z轴的夹角θk的估计值， 为所述第k个信号与xOy平面的夹角 的估计值，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0180] 通过本发明实施例提供的DOA估计方法，互素L型天线阵列在估计接收信号的波达方向DOA时，选取x轴接收信号的协方差矩阵进行特征值分解，得到x轴噪声子空间，而不再选取接收信号整个噪声子空间的一半来作为x轴噪声子空间，这样就充分利用了x轴接收信号的协方差信息，提高了算法估计的精准度，使该算法更满足实际应用场景的需要。
[0181] 基于同一种发明构思，根据上述本发明实施例提供的一种DOA估计方法，本发明实施例还提供了一种DOA估计装置。如图4所示，图4为本发明实施例提供的DOA估计装置的一种结构示意图。该装置应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。该装置包含以下单元。
[0182] 信号接收单元401，用于获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号。
[0183] 信号叠加单元402，用于叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号。
[0184] 矩阵计算单元403，用于确定两个L型天线阵列接收的整体信号的协方差矩阵。
[0185] 第一特征值分解单元404，用于对该协方差矩阵进行特征值分解，得到整体噪声子空间。
[0186] 矩阵提取单元405，用于提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵。
[0187] 第二特征值分解单元406，用于对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间。
[0188] DOA估计单元407，用于根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0189] 一个实施例中，信号叠加单元402具体可以用于：
[0190] 叠加整合第一x轴天线子阵列接收的信号和第二x轴天线子阵列接收的信号，得到x轴接收信号；
[0191] 叠加整合第一y轴天线子阵列接收的信号和第二y轴天线子阵列接收的信号，得到y轴接收信号；
[0192] 叠加整合x轴接收信号和所述y轴接收信号，得到两个L型天线阵列接收的整体信号。
[0193] 一个实施例中，矩阵计算单元403，具体可以用于：
[0194] 利用以下公式，确定两个L型天线阵列接收的整体信号的协方差矩阵[0195]
[0196] 其中，F为对所述整体信号的采样值个数，z(t)为所述两个L型天线阵列接收的整体信号的第t个采样值，zH(t)为z(t)的共轭转置。
[0197] 一个实施例中，DOA估计单元407，具体可以用于：
[0198] 根据x轴噪声子空间，构造第一谱函数；
[0199] 对第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0200] 根据每一信号的第一变量和整体噪声子空间，构造第二谱函数；
[0201] 对第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0202] 根据第一变量和第二变量，估计接收信号的DOA。
[0203] 通过本发明实施例提供的DOA估计装置，互素L型天线阵列在估计接收信号的波达方向DOA时，选取x轴接收信号的协方差矩阵进行特征值分解，得到x轴噪声子空间，而不再选取接收信号整个噪声子空间的一半来作为x轴噪声子空间，这样就充分利用了x轴接收信号的协方差信息，提高了算法估计的精准度，使该算法更满足实际应用场景的需要。
[0204] 参见图5，图5为本发明实施例提供的DOA估计装置的另一种结构示意图。该装置应用于与两个L型天线阵列连接的电子设备，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。该装置包含以下单元。
[0205] 信号接收单元501，用于获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号。单元501与单元401相同。
[0206] 信号叠加单元502，用于叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号。单元502与单元402相同。
[0207] 矩阵计算单元503，用于确定两个L型天线阵列接收的整体信号的协方差矩阵。单元503与单元403相同。
[0208] 第一特征值分解单元504，用于对该协方差矩阵进行特征值分解，得到整体噪声子空间。单元504与单元404相同。
[0209] 矩阵提取单元505，用于提取该协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵。单元505与单元405相同。
[0210] 第二特征值分解单元506，用于对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间。单元506与单元406相同。
[0211] DOA估计单元507，用于根据x轴噪声子空间，估计接收信号的DOA。单元507与单元407相同。
[0212] 第一谱函数构造子单元508，用于根据所述x轴噪声子空间，构造第一谱函数。
[0213] 第一峰值搜索子单元509，用于对所述第一谱函数进行峰值搜索，确定每一信号的第一变量的估计值；每一信号的第一变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0214] 第二谱函数构造子单元510，用于根据所述每一信号的第一变量和所述整体噪声子空间，构造第二谱函数；
[0215] 第二峰值搜索子单元511，用于对所述第二谱函数进行峰值搜索，确定每一信号的第二变量的估计值；每一信号的第二变量是根据该信号与z轴的夹角和该信号与xOy平面的夹角确定的变量；
[0216] DOA估算子单元512，用于根据所述第一变量和所述第二变量，估计接收信号的DOA。
[0217] 一个实施例中，第一谱函数构造子单元507，具体可以用于：
[0218] 利用以下公式，构造第一谱函数f(μ)：
[0219]
[0220] 其中， 为x轴噪声子空间，ax(μk)是x轴接收信号的导向矢量，μk为第k个信号的第一变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0221] 一个实施例中，第二谱函数构造子单元509，具体可以用于：
[0222] 利用以下公式，构造第二谱函数f(ν)：
[0223]
[0224] 其中， 为整体噪声子空间， 为第k个信号的第一变量μk的估计值 的导向矢量，ay(νk)是y轴接收信号的导向矢量，νk为第k个信号的第二变量，1≤k≤K，K表示接收到所述信号源发出的信号的总个数。
[0225] 一个实施例中，DOA估算子单元511，具体可以用于：
[0226] 利用以下公式，估计接收到的第k个信号与z轴的夹角θk和该信号与xOy平面的夹角
[0227]
[0228]
[0229] 其中， 为第k个信号的第一变量μk的估计值， 为第k个信号的第二变量νk的估计值， 为第k个信号与z轴的夹角θk的估计值， 为第k个信号与xOy平面的夹角 的估计值，1≤k≤K，K表示接收到信号源发出的信号的总个数；
[0230] 将该估计值 与 作为接收信号的DOA。
[0231] 通过本发明实施例提供的DOA估计装置，互素L型天线阵列在估计接收信号的波达方向DOA时，选取x轴接收信号的协方差矩阵进行特征值分解，得到x轴噪声子空间，而不再选取接收信号整个噪声子空间的一半来作为x轴噪声子空间，这样就充分利用了x轴接收信号的协方差信息，提高了算法估计的精准度，使该算法更满足实际应用场景的需要。
[0232] 基于同一种发明构思，根据上述本发明实施例提供的一种DOA估计方法，本发明实施例还提供了一种电子设备，如图6所示。该电子设备与两个L型天线阵列601相连接，所述两个L型天线阵列的原点相同，其中，一个L型天线阵列包括第一x轴天线子阵列和第一y轴天线子阵列，另一个L型天线阵列包括第二x轴天线子阵列和第二y轴天线子阵列，所述第一x轴天线子阵列和第一y轴天线子阵列分别包括N1个阵元，所述第二x轴天线子阵列和第二y轴天线子阵列分别包括N2个阵元，N1和N2互素。该电子设备包括转换装置602、处理器603、通信接口604、存储器605和通信总线606，其中，转换装置602、处理器603、通信接口604、存储器605和通信总线606完成相互间的通信。
[0233] 天线阵列601，用于接收信号源发出的模拟信号；
[0234] 转换装置602，用于将上述天线阵列接收的模拟信号转换成数字信号；
[0235] 存储器605，用于存放计算机程序；
[0236] 处理器603，用于执行存储器605上所存放的程序时，实现如下步骤：
[0237] 获取第一x轴天线子阵列、第一y轴天线子阵列、第二x轴天线子阵列和第二y轴天线子阵列分别接收的信号源发出的信号；
[0238] 叠加整合获取的各个信号，得到两个L型天线阵列接收的整体信号；
[0239] 确定两个L型天线阵列接收的整体信号的协方差矩阵；
[0240] 对所述协方差矩阵进行特征值分解，得到整体噪声子空间；
[0241] 提取协方差矩阵的前N1+N2行和前N1+N2列的元素，作为x轴协方差矩阵；
[0242] 对x轴协方差矩阵进行特征值分解，得到x轴噪声子空间；
[0243] 根据上述x轴噪声子空间和整体噪声子空间，估计接收信号的DOA。
[0244] 通过本发明实施例提供的电子设备，互素L型天线阵列在估计接收信号的波达方向DOA时，选取x轴接收信号的协方差矩阵进行特征值分解，得到x轴噪声子空间，而不再选取接收信号整个噪声子空间的一半来作为x轴噪声子空间，这样就充分利用了x轴接收信号的协方差信息，提高了算法估计的精准度，使该算法更满足实际应用场景的需要。
[0245] 上述电子设备提到的通信总线可以是外设部件互连标准(Peripheral Component Interconnect，PCI)总线或扩展工业标准结构(Extended Industry Standard Architecture，EISA)总线等。该通信总线可以分为地址总线、数据总线、控制总线等。为便于表示，图中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。
[0246] 通信接口用于上述电子设备与其他设备之间的通信。
[0247] 存储器可以包括随机存取存储器(Random Access Memory，RAM)，也可以包括非易失性存储器(Non-Volatile Memory，NVM)，例如至少一个磁盘存储器。可选的，存储器还可以是至少一个位于远离前述处理器的存储装置。
[0248] 上述的处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(Digital Signal Processing，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
[0249] 在本发明提供的又一实施例中，还提供了一种计算机可读存储介质，该计算机可读存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现上述任一DOA估计方法的步骤。
[0250] 在本发明提供的又一实施例中，还提供了一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行上述实施例中任一DOA估计方法。
[0251] 在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
[0252] 需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0253] 本说明书中的各个实施例均采用相关的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于DOA估计装置、电子设备、机器可读存储介质实施例而言，由于其基本相似于DOA估计方法实施例，所以描述的比较简单，相关之处参见DOA估计方法实施例的部分说明即可。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 非均匀阵列欠定波达方向估计方法 | 2020-05-12 | 911 |
| 电磁信号的二维波达方向角估计方法 | 2020-05-12 | 584 |
| 分布式子阵波达方向估计方法 | 2020-05-13 | 660 |
| 一种波达方向估计方法和装置 | 2020-05-11 | 643 |
| 一种宽带阵列波达方向估计方法 | 2020-05-11 | 816 |
| 一种基于空域自适应滤波的波达方向估计方法 | 2020-05-13 | 720 |
| 一种波达方向的检测方法及装置 | 2020-05-11 | 169 |
| 单次搜索宽带信号波达方向估计方法 | 2020-05-11 | 933 |
| 杂波背景下天线阵列的信号波达方向估计方法 | 2020-05-13 | 397 |
| 多信源波达方向估计方法 | 2020-05-12 | 935 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈