基于数字信号处理机的时域宽带信号频域闭环测向方法
【专利摘要】本发明公开了一种时域宽带信号频域闭环测向方法，其主要思路为：获取幅相误差校正后第m个通道中的数字宽带信号并进行离散傅立叶变换，得到幅相误差校正后第m个通道中的数字宽带信号在第j个频点处的频域信号Xm(ωj)，进而分别计算使用频域最小均方算法后第j个频点ωj处的频域窄带信号最优权矢量Wj和第j个频点ωj处的频域信号最优波束形成方向图yopt(ωj)，直到得到用频域最小均方算法后J个频点各自的频域窄带信号最优权矢量和J个频点各自的最优波束形成方向图，并进行几何平均，计算得到N个阵元接收到的时域宽带信号对应的最终波束形成方向图F，进而依次计算使用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数PLMS和PLMS的谱峰位置，最后得到时域宽带信号的角度估计值。
【专利说明】
基于数字信号处理机的时域宽带信号频域闭环测向方法
技术领域
[0001] 本发明属于通信技术领域，特别涉及一种基于数字信号处理机的时域宽带信号频 域闭环测向方法，适用于宽带信号的方向估计。
【背景技术】
[0002] 由于宽带信号入射到阵列的各个传感器上时，不但在相位上会产生差异，而且信 号的幅度(复包络)也会发生变化，使得许多窄带信号的测向算法不能应用到宽带信号的测 向中。目前，宽带信号的波达方向（Direction of Arrival，D0A)估计算法主要分为最大似 然类方法和基于信号子空间的方法，最大似然类方法的估计性能在高斯白噪声条件下是最 优的；Clark等提出的基于迭代二次型最大似然算法是一种宽带信号D0A最大似然估计方 法，该种宽带信号D0A最大似然估计方法将宽带信号近似为若干窄带正弦信号的叠加，并将 窄带ARMA模型应用到宽带信号的测向中，然后在若干窄带正弦信号的各个频点上通过对多 项式求根的方法得到相应波达方向估计值，最后通过最小二乘估计方法求得宽带信号的方 向最终估计值;针对该种宽带信号D0A最大似然估计方法过程复杂，并且计算量大的缺点， Agrawal等提出在宽带信号D0A最大似然估计方法中采用三维优化，降低了计算量，但是却 并没有在实际工程中得到应用。
[0003] 最大似然类方法运算过于复杂，并且估计结果容易在局部极值点收敛;相比于最 大似然类方法，基于信号子空间的方法虽然无法得到最优的估计结果，但却具有较低的运 算复杂度和较高的估计性能，因此成为主要的宽带信号高分辨处理技术。目前，基于空间信 号子空间的算法主要有三种:非相干信号子空间方法（ISSM)，相干信号子空间方法(CSSM) 和宽带直接处理法；ISSM利用经典的窄带子空间类高分辨算法(MUSIC)分别对各个频点下 的数据进行谱估计，该方法计算量大，且每个频点都要进行一次特征值分解，而且无法处理 相干源;基于ISSM的不足之处，Wang等提出了相干子空间方法;相干子空间方法使用一个聚 焦矩阵，将不同频率的协方差矩阵映射为中心频率的协方差矩阵，然后进行频域平滑得到 统一的协方差矩阵，最后应用窄带处理的方法估计宽带信号的波达方向；CSSM方法相比于 ISSM方法，具有运算量小、估计精度高、可实现对相干信号的处理等优点，其中聚焦矩阵的 构造是CSSM方法的关键，也是难点之一，CSSM方法的另一难点在于选择最佳的聚焦频率来 减小估计偏差，并且要进行波达方向预估计；宽带直接处理法中的典型算法包含宽带信号 子空间谱算法和基于频率模型算法。宽带直接处理方法有效避免了信号波达方向预估计和 聚焦矩阵的构造，且阵列的协方差矩阵维数的增大也更有效的描述了子空间的特性，因此 宽带直接处理法得到宽带信号的波达方向估计性能更好，但是宽带直接处理法计算量很 大，而且阵列误差对其影响明显，使得必须对阵列进行校准。
[0004] 以上的宽带信号波达方向估计方法都是基于开环算法，工程可应用性都不强。
【发明内容】
[0005] 针对以上现有技术存在的不足，本发明的目的在于提出一种基于数字信号处理机 的时域宽带信号频域闭环测向方法，该种基于数字信号处理机的时域宽带信号频域闭环测 向方法解决现有算法需要估计协方差矩阵，且运算复杂度高，以及并行射频通道数多从而 使得硬件设备复杂问题，并在保证准确估计宽带信号方位的基础上，减小运算复杂度，降低 工程上硬件复杂度。
[0006]为实现上述目的，本发明采用如下技术方案予以实现。
[0007]技术方案一：
[0008] 一种数字信号处理机，包括波束形成器、N选一开关、模拟接收机、模数转换器、数 据采集模块、数字正交插值单元、数据预处理模块、宽带信号闭环测向模块和N个耦合器，N 为自然数；
[0009] 所述数据预处理模块包含数据缓存单元和幅相误差校正单元，所述宽带信号闭环 测向模块包含宽带信号处理单元、最优权值生成单元和测向空间谱计算单元；
[0010]所述N选一开关包含N个输入端和一个输出端；
[0011] 所述模拟接收机包含第一模拟信号输入端、第二模拟信号输入端和模拟信号输出 端;所述波束形成器包含第1传输信号输入端至第(N+1)传输信号输入端，以及第一传输信 号输出端和第二传输信号输出端；
[0012] 所述模数转换器包含第三模拟信号输入端和第一数字信号输出端；
[0013] 所述数据采集模块包含第一数字信号输入端和第二数字信号输出端；
[0014] 所述数字正交插值单元包含第二数字信号输入端和第三数字信号输出端；
[0015] 所述数据缓存单元包含数字传输信号输入端和数字传输信号输出端；
[0016] 所述幅相误差校正单元包含数字传输信号幅相误差输入端和数字传输信号幅相 误差输出端；
[0017] 所述宽带信号处理单元包含宽带数字信号输入端和宽带数字信号输出端；
[0018] 所述最优权值生成单元包含窄带信号输入端和权矢量输出端；
[0019] 所述测向空间谱计算单元包含宽带信号角度估计值输出端；
[0020] 所述N个耦合器中第1个耦合器包含第1时域信号输出端和第1耦合信号输出端，所 述N个耦合器包含第1时域信号输出端至第N时域信号输出端，以及第1耦合信号输出端至第 N親合信号输出端，其中所述第1时域信号输出端至第N时域信号输出端分别电连接N选一开 关的N个输入端，所述第1耦合信号输出端至第N耦合信号输出端分别电连接波束形成器的 第1传输信号输入端至第N传输信号输入端，1 = 1，2，…，N;
[0021] N选一开关的输出端电连接模拟接收机的第一模拟信号输入端，模拟接收机的第 二模拟信号输入端电连接波束形成器的第一传输信号输出端，波束形成器的第二传输信号 输出端电连接测向空间谱计算单元的宽带信号角度估计值输出端;模拟接收机的模拟信号 输出端电连接模数转换器的第三模拟信号输入端，模数转换器的第一数字信号输出端电连 接数据采集模块的第一数字信号输入端，数据采集模块的第二数字信号输出端电连接数字 正交插值单元的第二数字信号输入端，数字正交插值单元的第三数字信号输出端电连接数 据缓存单元的数字传输信号输入端，数据缓存单元的数字传输信号输出端电连接幅相误差 校正单元的数字传输信号幅相误差输入端，幅相误差校正单元的数字传输信号幅相误差输 出端电连接宽带信号处理单元的宽带数字信号输入端，宽带信号处理单元的宽带数字信号 输出端电连接最优权值生成单元的窄带信号输入端，最优权值生成单元的权矢量输出端电 连接波束形成器的第(N+1)传输信号输入端；
[0022] 所述N个耦合器对应N个通路，分别用于获取对应阵元上的时域宽带信号，并对获 取的时域宽带信号分别进行耦合处理，得到N个通道耦合时域宽带信号，然后分别发送至波 束形成器和N选一开关；
[0023] 所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并进 行波束合成，得到合路时域信号，并将所述合路时域信号发送至模拟接收机；
[0024] 所述N选一开关，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并选 通其中一个通道耦合时域宽带信号，并发送至模拟接收机；
[0025] 所述模拟接收机用于接收所述一个通道耦合时域宽带信号和波束形成器发送过 来的合路时域信号，并计算所述一个通道耦合时域宽带信号和波束形成器发送过来的合路 时域信号之间的互相关，得到一个通道耦合时域宽带信号和合路时域信号之间的互相关， 然后发送至模数转换器；
[0026] 所述模数转换器用于接收一个通道耦合时域宽带信号和合路时域信号之间的互 相关，并将所述一个通道耦合时域宽带信号和合路时域信号之间的互相关转换为数字宽带 信号，发送至数据采集模块；
[0027] 所述数据采集模块用于接收模拟接收机发送过来的数字宽带信号，并对所述数字 宽带信号进行中频信号采集，得到中频数字宽带信号，并将所述中频数字宽带信号发送至 数字正交插值单元；
[0028] 所述数字正交插值单元用于接收数据采集模块发送过来的中频数字宽带信号，并 对所述中频数字宽带信号进行正交插值处理，得到正交插值处理后的数字宽带信号，并发 送至数据缓存模块；
[0029]所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理 后的数字宽带信号，同时将所述正交插值处理后的数字宽带信号发送至幅相误差校正单 元；
[0030] 所述幅相误差校正单元用于接收数据缓存模块缓存的所述正交插值处理后的数 字宽带信号，并对所述正交插值处理后的数字宽带信号进行幅相误差校正，得到幅相误差 校正后的数字宽带信号，然后发送至宽带信号处理单元；
[0031] 所述宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正后 的数字宽带信号，并对所述幅相误差校正后的数字宽带信号进行离散傅里叶变换，得到幅 相误差校正后的数字宽带信号在单个频点处的频域信号，然后发送至最优权值生成单元；
[0032] 所述最优权值生成单元用于接收宽带信号处理单元发送过来的幅相误差校正后 的数字宽带信号在单个频点处的频域信号，并计算使用频域最小均方算法后单个频点处的 频域窄带信号最优权矢量，并将所述单个频点处的频域窄带信号最优权矢量发送至波束形 成器计算单个频点处的最优波束形成方向图，进而计算数字宽带信号在所有频点处的各自 最优波束形成方向图，然后将所述数字宽带信号在所有频点处的各自最优波束形成方向图 发送至测向空间谱计算单元；
[0033] 所述测向空间谱计算单元用于接收波束形成器发送过来的数字宽带信号在所有 频点处的各自最优波束形成方向图，并计算接收到的时域宽带信号对应的最终波束形成方 向图，获取时域宽带信号闭环测向空间谱函数的谱峰位置，进而得到时域宽带信号的角度 估计值。
[0034] 技术方案二：
[0035] 一种时域宽带信号频域闭环测向方法，基于一种数字信号处理机，所述数字信号 处理机，包括N个耦合器、波束形成器、N选一开关、模拟接收机、模数转换器、数据采集模块、 数字正交插值单元、数据缓存单元、幅相误差校正单元、宽带信号处理单元、最优权值生成 权值单元和测向空间谱计算单元，所述时域宽带信号频域闭环测向方法，包括以下步骤：
[0036] 步骤1，获取空间N个阵元组成的天线阵列，并将每个阵元作为一个单通道，所述N 个阵元分别获取对应阵元上的时域宽带信号，并对获取的时域宽带信号分别进行耦合处 理，得到N个通道耦合时域宽带信号，然后分别发送至波束形成器和N选一开关;所述N选一 开关选取第m个通道耦合时域宽带信号并将接收到的第m个通道耦合时域宽带信号发送至 模拟接收机;所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号 并进行波束合成，得到合路时域信号，并将所述合路时域信号发送至模拟接收机;其中，m e {1，2,…，N}，N表示空间天线阵列包含的通道个数，N为自然数；
[0037] 所述模拟接收机用于接收第m个通道耦合时域宽带信号和所述合路时域信号，并 计算第m个通道耦合时域宽带信号和所述合路时域信号之间的互相关，然后将所述第m个通 道耦合时域宽带信号和所述合路时域信号之间的互相关发送至模数转换器;所述模数转换 器用于接收第m个通道耦合时域宽带信号和所述合路时域信号之间的互相关，并将第m个通 道耦合时域宽带信号和所述合路时域信号之间的互相关转换为第m个通道数字宽带信号， 并发送至数据采集模块；
[0038] 步骤2,所述数据采集模块用于接收模拟接收机发送过来的第m个通道数字宽带信 号，并对所述第m个通道数字宽带信号进行中频信号采集，得到第m个通道中频数字宽带信 号，并将第m个通道中频数字宽带信号发送至数字正交插值单元；
[0039]所述数字正交插值单元用于接收数据采集模块发送过来的第m个通道中频数字宽 带信号，并对所述第m个通道中频数字宽带信号进行正交插值处理，得到正交插值处理后第 m个通道的数字宽带信号，并发送至数据缓存模块；
[0040] 所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理 后第m个通道的数字宽带信号，同时将所述正交插值处理后第m个通道的数字宽带信号发送 至幅相误差校正单元；
[0041] 所述幅相误差校正单元用于接收数据缓存模块缓存的正交插值处理后第m个通道 的数字宽带信号，并对所述正交插值处理后第m个通道的数字宽带信号进行幅相误差校正， 得到幅相误差校正后第m个通道中的数字宽带信号，然后发送至宽带信号处理单元；
[0042]步骤3,宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正 后第m个通道中的数字宽带信号，并对所述幅相误差校正后第m个通道中的数字宽带信号进 行离散傅立叶变换，得到幅相误差校正后第m个通道中的数字宽带信号在第j个频点处的频 域信号X m(〇^)，并发送至最优权值生成单元;其中，mG{l，2,…，N}，N表示空间天线阵列包 含的通道个数，N为自然数;jG{l，2,…，J}，J表示幅相误差校正后第m个通道内的数字宽带 信号在宽带范围[wi wh]内等效为频域窄带信号叠加后包含的频点个数，们表示幅相误差校 正后第m个通道中的数字宽带信号的频率下界，w h表示幅相误差校正后第m个通道中的数字 宽带信号的频率上界；
[0043]步骤4,初始化：j表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围 [wi wh]内等效为频域窄带信号叠加后的第j个频点，且j的初始值为1，j G {1，2，…，J}，J表 示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 wh]内等效为频域窄带信号 叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信号的频率下界， wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界；
[0044]步骤5,最优权值生成单元接收宽带信号处理单元发送过来的幅相误差校正后第m 个通道中的数字宽带信号在第j个频点处的频域信号Xm(〇^)，并计算使用频域最小均方算 法后第j个频点处的频域窄带信号最优权矢量^，然后将使用频域最小均方算法后第j个 频点处的频域窄带信号最优权矢量W发送至波束形成器;波束形成器接收使用频域最小 均方算法后第j个频点处的频域窄带信号最优权矢量％并进行波束合成，计算得到N个通 道中的数字宽带信号在第j个频点处的最优波束形成方向图y〇 Pt( ;
[0045] 步骤6,令j加 1，返回步骤5,直到得到使用频域最小均方算法后第J个频点《 j处的 频域窄带信号最优权矢量Wj和N个通道中的数字宽带信号在第J个频点coj处的最优波束形 成方向图7_(〇^)，并将此时得到的1^个通道中的数字宽带信号在第1个频点《1的最优波束 形成方向图y〇 Pt( 到N个通道中的数字宽带信号在第J个频点coj处的最优波束形成方向 图3^*(〇^)，8_个通道中的数字宽带信号在】个频点各自的最优波束形成方向图，并分别 发送至测向空间谱计算单元。
[0046] 步骤7,测向空间谱计算单元接收波束形成器发送过来的N个通道中的数字宽带信 号在J个频点各自的最优波束形成方向图，并进行几何平均，计算得到N个阵元接收到的时 域宽带信号对应的最终波束形成方向图F;
[0047] 步骤8,根据N个阵元接收到的时域宽带信号对应的最终波束形成方向图F，计算得 到使用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数Plms，然后根据所述 时域宽带信号闭环测向空间谱函数Plms，获取时域宽带信号闭环测向空间谱函数Plms的谱峰 位置，进而得到时域宽带信号的角度估计值。
[0048] 本发明与现有技术相比具有以下优点：
[0049] 第一，本发明方法能够准确估计出一个信号源时宽带信号的角度，和非相干信号 子空间（ISSM)方法相比，本发明方法的测向空间谱更尖锐，测向效果更好；
[0050] 第二，本发明方法随着信噪比的增大，角度估计均方根误差越来越小，测角性能越 来越好相比已有的ISSM方法，本发明方法受信噪比的影响更小；
[0051] 第三，本发明方法进行角度估计时不需提前知道信号源的数目，没有阵型要求，对 于任意阵型均适用；
[0052] 第四，本发明方法不需要同时获取多通道数据，并能够分时通道切换获取数据样 本，实现原理简单，工程可应用性强，降低了运算量。
【附图说明】
[0053]下面结合附图和【具体实施方式】对本发明作进一步详细说明。
[0054]图1是本发明的一种时域宽带信号频域闭环测向的优化方法的实现原理图；
[0055] 图2是本发明的宽带信号入射示意图；
[0056] 图3是本发明中LMS频域自适应波束形成原理框图；
[0057] 图4是本发明中各频点运用LMS方法及几何平均后的仿真图；
[0058] 图5是本发明中宽带信号测向仿真图；
[0059] 图6是本发明中角度估计误差与信噪比关系仿真图；
[0060] 图7是本发明角度估计误差与闭环迭代次数关系仿真图。
【具体实施方式】
[0061] 参照图1，为本发明的一种时域宽带信号频域闭环测向的优化方法的实现原理图； 本发明的一种数字信号处理机，包括波束形成器、N选一开关、模拟接收机、模数转换器、数 据采集模块、数字正交插值单元、数据预处理模块、宽带信号闭环测向模块和N个耦合器，N 为自然数。
[0062]所述数据预处理模块包含数据缓存单元和幅相误差校正单元，所述宽带信号闭环 测向模块包含宽带信号处理单元、最优权值生成单元和测向空间谱计算单元。
[0063]所述N选一开关包含N个输入端和一个输出端。
[0064]所述模拟接收机包含第一模拟信号输入端、第二模拟信号输入端和模拟信号输出 端;所述波束形成器包含第1传输信号输入端至第(N+1)传输信号输入端，以及第一传输信 号输出端和第二传输信号输出端。
[0065]所述模数转换器包含第三模拟信号输入端和第一数字信号输出端。
[0066]所述数据采集模块包含第一数字信号输入端和第二数字信号输出端。
[0067]所述数字正交插值单元包含第二数字信号输入端和第三数字信号输出端。
[0068]所述数据缓存单元包含数字传输信号输入端和数字传输信号输出端。
[0069]所述幅相误差校正单元包含数字传输信号幅相误差输入端和数字传输信号幅相 误差输出端。
[0070] 所述宽带信号处理单元包含宽带数字信号输入端和宽带数字信号输出端。
[0071] 所述最优权值生成单元包含窄带信号输入端和权矢量输出端。
[0072] 所述测向空间谱计算单元包含宽带信号角度估计值输出端。
[0073] 所述N个耦合器中第1个耦合器包含第1时域信号输出端和第1耦合信号输出端，所 述N个耦合器包含第1时域信号输出端至第N时域信号输出端，以及第1耦合信号输出端至第 N親合信号输出端，其中所述第1时域信号输出端至第N时域信号输出端分别电连接N选一开 关的N个输入端，所述第1耦合信号输出端至第N耦合信号输出端分别电连接波束形成器的 第1传输信号输入端至第N传输信号输入端，1 = 1，2，…，N。
[0074] N选一开关的输出端电连接模拟接收机的第一模拟信号输入端，模拟接收机的第 二模拟信号输入端电连接波束形成器的第一传输信号输出端，波束形成器的第二传输信号 输出端电连接测向空间谱计算单元的宽带信号角度估计值输出端;模拟接收机的模拟信号 输出端电连接模数转换器的第三模拟信号输入端，模数转换器的第一数字信号输出端电连 接数据采集模块的第一数字信号输入端，数据采集模块的第二数字信号输出端电连接数字 正交插值单元的第二数字信号输入端，数字正交插值单元的第三数字信号输出端电连接数 据缓存单元的数字传输信号输入端，数据缓存单元的数字传输信号输出端电连接幅相误差 校正单元的数字传输信号幅相误差输入端，幅相误差校正单元的数字传输信号幅相误差输 出端通过数据总线电连接宽带信号处理单元的宽带数字信号输入端，宽带信号处理单元的 宽带数字信号输出端电连接最优权值生成单元的窄带信号输入端，最优权值生成单元的权 矢量输出端电连接波束形成器的第(N+1)传输信号输入端。
[0075] 所述N个耦合器对应N个通路，分别用于获取对应阵元上的时域宽带信号，并对获 取的时域宽带信号分别进行耦合处理，得到N个通道耦合时域宽带信号，然后分别发送至波 束形成器和N选一开关。
[0076] 所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并进 行波束合成，得到合路时域信号，并将所述合路时域信号发送至模拟接收机。
[0077] 所述N选一开关，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并选 通其中一个通道耦合时域宽带信号，并发送至模拟接收机。
[0078] 所述模拟接收机用于接收所述一个通道耦合时域宽带信号和波束形成器发送过 来的合路时域信号，并计算所述一个通道耦合时域宽带信号和波束形成器发送过来的合路 时域信号之间的互相关，得到一个通道耦合时域宽带信号和合路时域信号之间的互相关， 然后发送至模数转换器。
[0079] 所述模数转换器用于接收一个通道耦合时域宽带信号和合路时域信号之间的互 相关，并将所述一个通道耦合时域宽带信号和合路时域信号之间的互相关转换为数字宽带 信号，发送至数据采集模块。
[0080] 所述数据采集模块用于接收模拟接收机发送过来的数字宽带信号，并对所述数字 宽带信号进行中频信号采集，得到中频数字宽带信号，并将所述中频数字宽带信号发送至 数字正交插值单元。
[0081] 所述数字正交插值单元用于接收数据采集模块发送过来的中频数字宽带信号，并 对所述中频数字宽带信号进行正交插值处理，得到正交插值处理后的数字宽带信号，并发 送至数据缓存模块。
[0082]所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理 后的数字宽带信号，同时将所述正交插值处理后的数字宽带信号发送至幅相误差校正单 J L 〇
[0083]所述幅相误差校正单元用于接收数据缓存模块缓存的所述正交插值处理后的数 字宽带信号，并对所述正交插值处理后的数字宽带信号进行幅相误差校正，得到幅相误差 校正后的数字宽带信号，然后发送至宽带信号处理单元。
[0084]所述宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正后 的数字宽带信号，并对所述幅相误差校正后的数字宽带信号进行离散傅里叶变换，得到幅 相误差校正后的数字宽带信号在单个频点处的频域信号，然后发送至最优权值生成单元。 [0085]所述最优权值生成单元用于接收宽带信号处理单元发送过来的幅相误差校正后 的数字宽带信号在单个频点处的频域信号，并计算使用频域最小均方算法后单个频点处的 频域窄带信号最优权矢量，并将所述单个频点处的频域窄带信号最优权矢量发送至波束形 成器计算单个频点处的最优波束形成方向图，进而计算数字宽带信号在所有频点处的各自 最优波束形成方向图，然后将所述数字宽带信号在所有频点处的各自最优波束形成方向图 发送至测向空间谱计算单元。
[0086]所述测向空间谱计算单元用于接收波束形成器发送过来的数字宽带信号在所有 频点处的各自最优波束形成方向图，并计算接收到的时域宽带信号对应的最终波束形成方 向图，获取时域宽带信号闭环测向空间谱函数的谱峰位置，进而得到时域宽带信号的角度 估计值。
[0087] 一种时域宽带信号频域闭环测向方法，基于一种数字信号处理机，所述数字信号 处理机，包括N个耦合器、波束形成器、N选一开关、模拟接收机、模数转换器、数据采集模块、 数字正交插值单元、数据预处理模块和宽带信号闭环测向模块;所述数据预处理模块包含 数据缓存单元和幅相误差校正单元，所述宽带信号闭环测向模块包含宽带信号处理单元、 最优权值生成权值单元和测向空间谱计算单元，所述宽带信号频域闭环测向的优化方法， 包括以下步骤；
[0088] 步骤1，获取空间N个阵元组成的天线阵列，并将每个阵元作为一个单通道，所述N 个阵元分别获取对应阵元上的时域宽带信号，并对获取的时域宽带信号分别进行耦合处 理，得到N个通道耦合时域宽带信号，然后分别发送至波束形成器和N选一开关;所述N选一 开关选取第m个通道耦合时域宽带信号并将接收到的第m个通道耦合时域宽带信号发送至 模拟接收机;所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号 并进行波束合成，得到合路时域宽带信号，并将所述合路时域宽带信号发送至模拟接收机； 其中，mG{l，2,…，N}，N表示空间天线阵列包含的通道个数，N为自然数。
[0089] 所述模拟接收机用于接收第m个通道耦合时域宽带信号和所述合路时域宽带信 号，并计算第m个通道耦合时域宽带信号和所述合路时域宽带信号之间的互相关，然后将所 述第m个通道耦合时域宽带信号和所述合路时域宽带信号之间的互相关发送至模数转换 器;所述模数转换器用于接收第m个通道耦合时域宽带信号和所述合路时域宽带信号之间 的互相关，并将第m个通道耦合时域宽带信号和所述合路时域宽带信号之间的互相关转换 为第m个通道数字宽带信号，并发送至数据采集模块。
[0090] 具体地，获取空间N个阵元组成的天线阵列，参照图2,为本发明的宽带信号入射示 意图；在图2中，1、2、N表示天线阵列包含的阵元序号，d表示阵元间距，0 P表示第p个信号源 在每个阵元处的入射角;本发明使用N个阵元接收空间中入射的时域宽带信号，每个阵元作 为一个单通道，每个阵元包含作为一个传感器，用于检测所接收到的时域宽带信号。由于时 域宽带信号入射到阵列天线的多个传感器上时不同阵元上时，不但在相位上产生差异，而 且数字宽带信号的幅度或复包络也会发生变化，因此不能直接对数字宽带信号进行测向， 也不能直接利用数字宽带信号的测向方法对数字宽带信号进行测向；此处N=5。
[0091]步骤2,所述数据采集模块用于接收模拟接收机发送过来的第m个通道数字宽带信 号，并对所述第m个通道数字宽带信号进行中频信号采集，得到第m个通道中频数字宽带信 号，并将第m个通道中频数字宽带信号发送至数字正交插值单元。
[0092]所述数字正交插值单元用于接收数据采集模块发送过来的第m个通道中频数字宽 带信号，并对所述第m个通道中频数字宽带信号进行正交插值处理，得到正交插值处理后第 m个通道的数字宽带信号，并发送至数据缓存模块。
[0093]所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理 后第m个通道的数字宽带信号，同时将所述正交插值处理后第m个通道的数字宽带信号发送 至幅相误差校正单元；
[0094]所述幅相误差校正单元用于接收数据缓存模块缓存的正交插值处理后第m个通道 的数字宽带信号，并对所述正交插值处理后第m个通道的数字宽带信号进行幅相误差校正， 得到幅相误差校正后第m个通道中的数字宽带信号，然后发送至宽带信号处理单元。
[0095]步骤3,宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正 后第m个通道中的数字宽带信号，并对所述幅相误差校正后第m个通道中的数字宽带信号进 行离散傅立叶变换(DFT)，得到幅相误差校正后第m个通道中的数字宽带信号在第j个频点 处的频域信号X m(〇^)，并发送至最优权值生成单元;其中，mG{l，2,…，N}，N表示空间天线 阵列包含的通道个数，N为自然数;j G {1，2，…，J}，J表示幅相误差校正后第m个通道内的数 字宽带信号在宽带范围[wi wh]内等效为频域窄带信号叠加后包含的频点个数，们表示幅相 误差校正后第m个通道中的数字宽带信号的频率下界，wh表示幅相误差校正后第m个通道中 的数字宽带信号的频率上界。
[0096]具体地，宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正 后第m个通道中的数字宽带信号，并对所述幅相误差校正单元发送过来的幅相误差校正后 第m个通道中的数字宽带信号进行M'点离散傅立叶变换(DFT)，将所述幅相误差校正单元发 送过来的幅相误差校正后第m个通道中的数字宽带信号在宽带范围[ W1 wh]内划分为J个频 点各自对应频域窄带信号的叠加，们表示幅相误差校正后第m个通道中的数字宽带信号的 频率下界， Wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界，然后计算得到 幅相误差校正后第m个通道中的数字宽带信号在第j个频点处的频域信号X m(〇^)，其表达式 为：
[0098]其中，pG{l，2,…，P}，P表示第m个通道接收到的时域宽带信号包含的信号源数 目，j G {1，2，…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 Wh] 内等效为频域窄带信号叠加后包含的频点个数，SP(〇^)表示第p个信号源在第j个频点处的 频域信号，N m( coj表示第m个通道在第j个频点处的频域噪声，表示第j个频点，Tpm表示第 P个信号源到达第m个通道时相对于参考通道的时延，参考通道为N个通道中第一个通道， exp( ?)表示?的指数函数，mG{l，2,…，N}，N表示空间天线阵列包含的通道个数，M'表示 对幅相误差校正后第m个通道接收的时域宽带信号进行离散傅立叶变换时包含的点数，M' >J，且M'、J、N分别为自然数。
[0099]步骤4,初始化：j表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围 [wi wh]内等效为频域窄带信号叠加后的第j个频点，且j的初始值为1，j G {1，2，…，J}，J表 示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 wh]内等效为频域窄带信号 叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信号的频率下界， wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界。
[0100]步骤5,最优权值生成单元接收宽带信号处理单元发送过来的幅相误差校正后第m 个通道中的数字宽带信号在第j个频点处的频域信号Xm(〇^)，并使用频域最小均方算法 (LMS)，计算使用频域最小均方算法后第j个频点处的频域窄带信号最优权矢量然后 将使用频域最小均方算法后第j个频点处的频域窄带信号最优权矢量^发送至波束形成 器;波束形成器接收使用频域最小均方算法后第j个频点处的频域窄带信号最优权矢量 Wj并进行波束合成，计算得到N个通道中的数字宽带信号在第j个频点处的最优波束形成 方向图y〇pt ( 〇 j) d
[0101] 具体地，所述使用频域最小均方算法后第j个频点处的频域窄带信号最优权矢 量^和所述N个通道中的数字宽带信号在第j个频点处的最优波束形成方向图 ycipt( coj)， 其得到过程为：
[0102] 5.1初始化:n为迭代次数，且n初值为1，设定第j个频点《 j处的频域最小均方算法 调零门限为为空间天线阵列中第m个通道，且m初值为1，!11￡{1，2，一，《4表示空间天 线阵列包含的通道个数，表示第1次迭代后第1个通道中的数字宽带信号在第j个频点处 的初始权值Wj( 1，1 )，即Wj( 1，1 ) = Wq，并将初始权值Wj( 1，1 )作为静态权值Wq，Xl( ? j)表示第 1个通道中的数字宽带信号在第j个频点处的频域信号。
[0103] 5.2计算得到第n次迭代后第m个通道中的数字宽带信号在第j个频点0^处的权值 wj(m，n)，其表达式为：
[0104] wj(m,n+l)=wj(m,n)-2iiXm( ?j)e( ?j)
[0105] e(c〇j)=y(c〇j)-d(c〇j)
[0106] 其中，d(c^)表示第j个频点的频域期望信号，y(A)表示第j个频点的波束 形成方向图。
[0107] 5.3令m加1，重复子步骤5.2,直到得到第n次迭代后第N个通道中的数字宽带信号 在第j个频点处的权值巧(N，n)，此时计算第n次迭代后N个通道中的数字宽带信号在第j 个频点处的权值％(n)，并将所述第n次迭代后N个通道中的数字宽带信号在第j个频点 ?」处的权值W(n)发送至波束形成器，计算得到第n次迭代后N个通道中的数字宽带信号在 第j个频点处的波束形成方向图h(n)，然后计算得到第n次迭代后第j个频点处的频 域窄带信号自相关尺（"），其表达式分别为：
[0108] ffj(n) = [wj( 1 ,n), ??? ,wj(m,n), ??? ,wj(N,n)]
[0109] yj(n) =ffj(n)HX( ?j)
[0110] Ryj(n)=E[ffj(n)HX(?j)X(?j)Hffj(n)]
[0111] 其中，wj(m,n)表示第n次迭代后第m个通道中的数字宽带信号在第j个频点《 j处的 权值，X(coJ表示N个通道中的数字宽带信号在第j个频点处的频域信号，X(c^) = [Xl (A)，…，Xm( %)，…，XN( c^) ]T，Xm( %)表示第m个通道中的数字宽带信号在第j个频点处 的频域信号。
[0112] 5.4如果第n次迭代后第j个频点％处的频域窄带信号自相关％(?)大于第j个频 点处的频域最小均方算法调零门限\，则令n加1，返回子步骤5.2;
[0113] 如果第n次迭代后第j个频点％处的频域窄带信号自相关'(?)小于第j个频点《 j 处的频域最小均方算法调零门限则将此时得到的第n次迭代后N个通道中的数字宽带信 号在第j个频点《』处的权值W』(n)，作为使用频域最小均方算法后得到的第j个频点处的 频域窄带信号最优权矢量将第n次迭代后N个通道中的数字宽带信号在第j个频点coj处 的波束形成方向图n( n)，作为N个通道中的数字宽带信号在第j个频点处的最优波束形 成方向图y〇Pt( w」），Wj = [Wlj，…，Wmj，"_，WNj]，(你/:)二灰产尤⑷,.），X( ?j)表示N个通道中的 数字宽带信号在第j个频点处的频域信号，表示第m个通道中的数字宽带信号在第j个频 点处的权值。
[0114] 然后，参照图3,为本发明中LMS频域自适应波束形成原理框图；对幅相误差校正后 第m个通道内的数字宽带信号在宽带范围[W1 Wh]内等效为频域窄带信号叠加后包含的J个 频点，任意选取第j个频点，使用多通道选择开关，利用频域LMS方法进行第j个频点co j处 的波束形成，计算得到第j个频点处的频域信号的波束形成方向图y(?J，其表达式为:y (c^)=W/X(〇^)，jG{l，2,…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在 宽带范围[W1 Wh]内等效为频域窄带信号叠加后包含的频点个数，W」表示使用频域最小均方 算法后第j个频点处的频域窄带信号最优权值，N表示空间天线阵列包含的通道个数，X (?』)表示N个通道中的数字宽带信号在第j个频点处的频域信号。
[0115]假设第j个频点％的频域期望信号为d( %)，且第j个频点％的频域期望信号d (?j)第j个频点的波束形成方向图yOj)之间的误差为e( Wj)，则e( ?j)=y( ?j)-d (wj)，误差e( w j)的均方误差为E|e( 〇j) |2,其表达式为：
[0117]其中，上标H表示共辄转置，E[ ?]表示求?的均值，Re[ ?]表示取?的实部，Wj表 示使用频域最小均方算法后得到的第j个频点处的频域窄带信号最优权值，％表示第j 个频点％处的频域信号的自相关矩阵，且心^表示第j个频点《 j 处的频域信号与第j个频点的频域期望信号之间的互相关信号，d(c^)表示第j个频点 ?j的频域期望信号，X( coj表示第j个频点处的频域信号。
[0118]利用梯度下降法，对误差e ( )的均方误差为E | e ( ) |2关于％求导：
[0120] 用瞬时值2X(c^)，（a)代替稳态值2￡[乂（0^)，（0^)]，计算得到第」个频点〇^ 处的偏导值免^，其表达式为：
[0121] Vff：i =2X(0)j)e/
[0122] 其中，q表示第j个频点的频域期望信号d(c^)与第j个频点的波束形成方 向图y( %)之间的误差e( %)，X( %)表示第j个频点处的频点％处的频域信号，3表示求 导运算符，上标h表示共辄转置，（? r表示求?的共辄。
[0123] 步骤6,令j加1，返回步骤5,直到得到使用频域最小均方算法后第J个频点coj处的 频域窄带信号最优权矢量Wj和N个通道中的数字宽带信号在第J个频点coj处的最优波束形 成方向图7_(〇^)，并将此时得到的1^个通道中的数字宽带信号在第1个频点《1的最优波束 形成方向图y〇 Pt( 到N个通道中的数字宽带信号在第J个频点coj处的最优波束形成方向 图ycipt(c〇j)，即N个通道中的数字宽带信号在J个频点各自的最优波束形成方向图发送至测 向空间谱计算单元。
[0124] 步骤7,测向空间谱计算单元波束形成器发送过来的N个通道中的数字宽带信号在 J个频点各自的最优波束形成方向图，并进行几何平均，计算得到N个阵元接收到的时域宽 带信号对应的最终波束形成方向图F。
[0125]具体地，所述N个阵元接收到的时域宽带信号对应的最终波束形成方向图F，其表 达式为：
[0127] 其中，II表示联乘运算，yc>Pt(c^)表示N个通道中的数字宽带信号在第j个频点 处的最优波束形成方向图，W」表示使用频域最小均方算法后第j个频点处的频域窄带信 号最优权值，X」表示N个通道中的数字宽带信号在第j个频点处的频域信号X(〇^)，je{l， 2，…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[ W1 Wh]内等效为 频域窄带信号叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信 号的频率下界， Wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界。
[0128] 步骤8,根据N个阵元接收到的时域宽带信号对应的最终波束形成方向图F，计算得 到使用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数Plms，然后根据所述 时域宽带信号闭环测向空间谱函数Plms，获取时域宽带信号闭环测向空间谱函数Plms的谱峰 位置，进而得到时域宽带信号的角度估计值。
[0129] 具体地，分别假定第j个频点处的频域窄带信号的导向矢量为^(0)，第j个频点 处的信号子空间为，第j个频点处的噪声子空间为，则此时频域窄带信号最优 权值Wj与第j个频点co」处的频域窄带信号的导向矢量aj (0)的向量积Pj近似为〇，Pj = WjHaj (9)，则使用频域最小均方算法后得到的第j个频点处的频域窄带信号最优权值^垂直于 第j个频点处的信号子空间，频域窄带信号最优权值％作为第j个频点处的噪声子 空间的一个矢量，将频域窄带信号最优权值W代替第j个频点％处的噪声子空间从%， 计算得到第j个频点处的谱峰搜索空间谱函数匕(9)，其表达式为：
[0131]其中，W」表示使用频域最小均方算法后第j个频点处的频域窄带信号最优权值 Wj，aj (0)表示第j个频点《 j处的阵列导向矢量，0表示谱峰搜索空间谱函数Pj(0)进行谱峰 搜索的扫描角度，9 e [-90° 90°]。
[0132]将J个频点分别进行频域最小均方算法后得到的最终波束形成方向图数据F代替 第j个频点处的谱峰搜索空间谱函数匕(9)中的a/(0)W，/^(0)，计算得到使用频域最小 均方算法后得到的时域宽带信号闭环测向空间谱函数Plms，其表达式为：
[0134]其中，II表示联乘运算，yc>Pt(c^)表示N个通道中的数字宽带信号在第j个频点 处的最优波束形成方向图，W」表示使用频域最小均方算法后第j个频点处的频域窄带信 号最优权值，X(coJ表示N个通道中的数字宽带信号在第j个频点处的频域信号，je{l， 2，…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 Wh]内等效为 频域窄带信号叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信 号的频率下界，Wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界。
[0135] 最后根据使用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数 Plms，获取时域宽带信号闭环测向空间谱函数Plms的谱峰位置，进而得到时域宽带信号的角 度估计值。
[0136] 本发明的效果可通过以下仿真进一步说明：
[0137] ( - )仿真参数：
[0138] 选取5个阵元全向均匀线阵，信号方向为-10度，信噪比为20dB，中心频率为 100MHz，带宽为40MHz，利用DFT算法得到信号带宽内的35个频率点的窄带成分，阵元间距为 中心频率对应波长的一半，
；静态权值wq=[l 1 1 1 1]T，迭代次数为50次，单通道 接收的时域宽带信号进行离散傅立叶变换时包含的点数M'为256,快拍数256次。
[0139] (二)仿真内容及结果：
[0140]仿真1，采用本发明对DFT后J个频点分别应用频域LMS方法进行波束形成，并对J个 频点的频域信号进行波束形成后结果求几何平均得到最终给的方向图。如图4是各个频点 进行波束形成后的方向图和几何平均后的方向图。
[0141] 参照图4,为本发明中各频点运用LMS方法及几何平均后的仿真图；从图4可看出， 对于单个宽带信号，对各频点的波束形成方向图进行几何平均后在信号位置出现一个明显 凹陷，且凹陷位置深度在40dB以上。
[0142] 仿真2,采用本发明算法和已提出的ISSM方法分别对宽带信号进行测向，测向结果 如图5所示，图5是本发明中宽带信号测向仿真图；由图5可看出，利用本发明算法和已提出 的ISSM方法均可以对宽带信号实现测向，相比较而言，本发明测向空间谱更尖锐，测向效果 更好。
[0143] 仿真3,将信号的信噪比从OdB到30dB以间隔5dB变化，每次变化做100次蒙特卡洛 实验，其他仿真条件不变;分别计算两种方法在不同信噪比下的角度估计均方根误差，求角 度估计均方根误差公式为
其中CMn表示蒙特卡洛实验 的次数;本发明的测角性能受信噪比影响仿真图为图6所示，图6为本发明中角度估计误差 与信噪比关系仿真图。
[0144] 由图6可以看出本发明中，随着信噪比的增大，角度估计均方根误差越来越小，测 角性能越来越好，并且与已提出的ISSM方法相比较，本发明算法受信噪比的影响更小。
[0145] 仿真4,将闭环迭代次数从10到50以间隔5变化，每次变化做100次蒙特卡洛实验， 测角性能受LMS闭环迭代次数的影响仿真图如图7所示，图7为本发明角度估计误差与闭环 迭代次数关系仿真图。
[0146] 由图7可以看出，本发明中角度误差随LMS迭代次数的增加而减小，当迭代次数为 25次时，信号的角度估计误差为0,可以得到信号角度的准确估计值。
[0147] 仿真5，本发明算法与已提出的ISSM方法运算复杂度比较
[0148] 由图6可知，本发明算法在迭代25次时即可收敛，使用本发明方法进行测向到收敛 时所需的复乘与复加次数，与使用ISSM算法进行测向到收敛时所需的复乘与复加次数的比 较，如表1所示。
[0149] 表1
[0151 ]由表1可看出本发明方法中的复乘和复加次数，分别明显小于使用ISSM算法的复 乘和复加次数，相比于ISSM算法，本发明方法的运算复杂度大大减小。
[0152]综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。
[0153]显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精 神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围 之内，则本发明也意图包含这些改动和变型在内。
【主权项】
1. 一种数字信号处理机，其特征在于，包括波束形成器、N选一开关、模拟接收机、模数 转换器、数据采集模块、数字正交插值单元、数据预处理模块、宽带信号闭环测向模块和N个 耦合器，N为自然数； 所述数据预处理模块包含数据缓存单元和幅相误差校正单元，所述宽带信号闭环测向 模块包含宽带信号处理单元、最优权值生成单元和测向空间谱计算单元； 所述N选一开关包含N个输入端和一个输出端； 所述模拟接收机包含第一模拟信号输入端、第二模拟信号输入端和模拟信号输出端； 所述波束形成器包含第1传输信号输入端至第(N+1)传输信号输入端，以及第一传输信号输 出端和第二传输信号输出端； 所述模数转换器包含第三模拟信号输入端和第一数字信号输出端； 所述数据采集模块包含第一数字信号输入端和第二数字信号输出端； 所述数字正交插值单元包含第二数字信号输入端和第三数字信号输出端； 所述数据缓存单元包含数字传输信号输入端和数字传输信号输出端； 所述幅相误差校正单元包含数字传输信号幅相误差输入端和数字传输信号幅相误差 输出端； 所述宽带信号处理单元包含宽带数字信号输入端和宽带数字信号输出端； 所述最优权值生成单元包含窄带信号输入端和权矢量输出端； 所述测向空间谱计算单元包含宽带信号角度估计值输出端； 所述N个耦合器中第1个耦合器包含第1时域信号输出端和第1耦合信号输出端，所述N 个耦合器包含第1时域信号输出端至第N时域信号输出端，以及第1耦合信号输出端至第N耦 合信号输出端，其中所述第1时域信号输出端至第N时域信号输出端分别电连接N选一开关 的N个输入端，所述第1耦合信号输出端至第N親合信号输出端分别电连接波束形成器的第1 传输信号输入端至第N传输信号输入端，1 = 1，2，…，N; N选一开关的输出端电连接模拟接收机的第一模拟信号输入端，模拟接收机的第二模 拟信号输入端电连接波束形成器的第一传输信号输出端，波束形成器的第二传输信号输出 端电连接测向空间谱计算单元的宽带信号角度估计值输出端;模拟接收机的模拟信号输出 端电连接模数转换器的第三模拟信号输入端，模数转换器的第一数字信号输出端电连接数 据采集模块的第一数字信号输入端，数据采集模块的第二数字信号输出端电连接数字正交 插值单元的第二数字信号输入端，数字正交插值单元的第三数字信号输出端电连接数据缓 存单元的数字传输信号输入端，数据缓存单元的数字传输信号输出端电连接幅相误差校正 单元的数字传输信号幅相误差输入端，幅相误差校正单元的数字传输信号幅相误差输出端 电连接宽带信号处理单元的宽带数字信号输入端，宽带信号处理单元的宽带数字信号输出 端电连接最优权值生成单元的窄带信号输入端，最优权值生成单元的权矢量输出端电连接 波束形成器的第(N+1)传输信号输入端。2. 如权利要求1所述的一种数字信号处理机，其特征在于，所述N个耦合器对应N个通 路，分别用于获取对应阵元上的时域宽带信号，并对获取的时域宽带信号分别进行耦合处 理，得到N个通道耦合时域宽带信号，然后分别发送至波束形成器和N选一开关； 所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并进行波 束合成，得到合路时域信号，并将所述合路时域信号发送至模拟接收机； 所述N选一开关，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并选通其 中一个通道耦合时域宽带信号，并发送至模拟接收机； 所述模拟接收机用于接收所述一个通道耦合时域宽带信号和波束形成器发送过来的 合路时域信号，并计算所述一个通道耦合时域宽带信号和波束形成器发送过来的合路时域 信号之间的互相关，得到一个通道耦合时域宽带信号和合路时域信号之间的互相关，然后 发送至模数转换器； 所述模数转换器用于接收一个通道耦合时域宽带信号和合路时域信号之间的互相关， 并将所述一个通道耦合时域宽带信号和合路时域信号之间的互相关转换为数字宽带信号， 发送至数据采集模块； 所述数据采集模块用于接收模拟接收机发送过来的数字宽带信号，并对所述数字宽带 信号进行中频信号采集，得到中频数字宽带信号，并将所述中频数字宽带信号发送至数字 正交插值单元。3. 如权利要求2所述的一种数字信号处理机，其特征在于，所述数字正交插值单元用于 接收数据采集模块发送过来的中频数字宽带信号，并对所述中频数字宽带信号进行正交插 值处理，得到正交插值处理后的数字宽带信号，并发送至数据缓存模块； 所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理后的 数字宽带信号，同时将所述正交插值处理后的数字宽带信号发送至幅相误差校正单元； 所述幅相误差校正单元用于接收数据缓存模块缓存的所述正交插值处理后的数字宽 带信号，并对所述正交插值处理后的数字宽带信号进行幅相误差校正，得到幅相误差校正 后的数字宽带信号，然后发送至宽带信号处理单元； 所述宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正后的数 字宽带信号，并对所述幅相误差校正后的数字宽带信号进行离散傅里叶变换，得到幅相误 差校正后的数字宽带信号在单个频点处的频域信号，然后发送至最优权值生成单元； 所述最优权值生成单元用于接收宽带信号处理单元发送过来的幅相误差校正后的数 字宽带信号在单个频点处的频域信号，并计算使用频域最小均方算法后单个频点处的频域 窄带信号最优权矢量，并将所述单个频点处的频域窄带信号最优权矢量发送至波束形成器 计算单个频点处的最优波束形成方向图，进而计算数字宽带信号在所有频点处的各自最优 波束形成方向图，然后将所述数字宽带信号在所有频点处的各自最优波束形成方向图发送 至测向空间谱计算单元； 所述测向空间谱计算单元用于接收波束形成器发送过来的数字宽带信号在所有频点 处的各自最优波束形成方向图，并计算接收到的时域宽带信号对应的最终波束形成方向 图，获取时域宽带信号闭环测向空间谱函数的谱峰位置，进而得到时域宽带信号的角度估 计值。4. 一种时域宽带信号频域闭环测向方法，所述方法应用于如权利要求1-3中任一项所 述的一种数字信号处理机，所述数字信号处理机，包括N个耦合器、波束形成器、N选一开关、 模拟接收机、模数转换器、数据采集模块、数字正交插值单元、数据缓存单元、幅相误差校正 单元、宽带信号处理单元、最优权值生成权值单元和测向空间谱计算单元，所述时域宽带信 号频域闭环测向方法，其特征在于，包括以下步骤： 步骤1，获取空间N个阵元组成的天线阵列，并将每个阵元作为一个单通道，所述N个阵 元分别获取对应阵元上的时域宽带信号，并对获取的时域宽带信号分别进行耦合处理，得 到N个通道耦合时域宽带信号，然后分别发送至波束形成器和N选一开关;所述N选一开关选 取第m个通道耦合时域宽带信号并将接收到的第m个通道耦合时域宽带信号发送至模拟接 收机;所述波束形成器，用于接收N个耦合器发送过来的N个通道耦合时域宽带信号并进行 波束合成，得到合路时域信号，并将所述合路时域信号发送至模拟接收机；其中，me {1， 2,…，N}，N表示空间天线阵列包含的通道个数，N为自然数； 所述模拟接收机用于接收第m个通道耦合时域宽带信号和所述合路时域信号，并计算 第m个通道耦合时域宽带信号和所述合路时域信号之间的互相关，然后将所述第m个通道耦 合时域宽带信号和所述合路时域信号之间的互相关发送至模数转换器;所述模数转换器用 于接收第m个通道耦合时域宽带信号和所述合路时域信号之间的互相关，并将第m个通道耦 合时域宽带信号和所述合路时域信号之间的互相关转换为第m个通道数字宽带信号，并发 送至数据采集模块； 步骤2,所述数据采集模块用于接收模拟接收机发送过来的第m个通道数字宽带信号， 并对所述第m个通道数字宽带信号进行中频信号采集，得到第m个通道中频数字宽带信号， 并将第m个通道中频数字宽带信号发送至数字正交插值单元； 所述数字正交插值单元用于接收数据采集模块发送过来的第m个通道中频数字宽带信 号，并对所述第m个通道中频数字宽带信号进行正交插值处理，得到正交插值处理后第m个 通道的数字宽带信号，并发送至数据缓存模块； 所述数据缓存模块用于接收并缓存数字正交插值单元发送过来的正交插值处理后第m 个通道的数字宽带信号，同时将所述正交插值处理后第m个通道的数字宽带信号发送至幅 相误差校正单元； 所述幅相误差校正单元用于接收数据缓存模块缓存的正交插值处理后第m个通道的数 字宽带信号，并对所述正交插值处理后第m个通道的数字宽带信号进行幅相误差校正，得到 幅相误差校正后第m个通道中的数字宽带信号，然后发送至宽带信号处理单元； 步骤3,宽带信号处理单元用于接收幅相误差校正单元发送过来的幅相误差校正后第m 个通道中的数字宽带信号，并对所述幅相误差校正后第m个通道中的数字宽带信号进行离 散傅立叶变换，得到幅相误差校正后第m个通道中的数字宽带信号在第j个频点处的频域信 号X m( )，并发送至最优权值生成单元；j G {1，2，…，J}，J表示幅相误差校正后第m个通道 内的数字宽带信号在宽带范围[W1 wh]内等效为频域窄带信号叠加后包含的频点个数，们表 示幅相误差校正后第m个通道中的数字宽带信号的频率下界，wh表示幅相误差校正后第m个 通道中的数字宽带信号的频率上界； 步骤4,初始化:j表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 wh] 内等效为频域窄带信号叠加后的第j个频点，且j的初始值为l，jG{l，2,…，J}，J表示幅相 误差校正后第m个通道内的数字宽带信号在宽带范围[ W1 wh]内等效为频域窄带信号叠加后 包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信号的频率下界，wh表示 幅相误差校正后第m个通道中的数字宽带信号的频率上界； 步骤5,最优权值生成单元接收宽带信号处理单元发送过来的幅相误差校正后第m个通 道中的数字宽带信号在第j个频点处的频域信号Xm(〇^)，并计算使用频域最小均方算法后 第j个频点《」处的频域窄带信号最优权矢量^，然后将使用频域最小均方算法后第j个频点 处的频域窄带信号最优权矢量％发送至波束形成器;波束形成器接收使用频域最小均方 算法后第j个频点处的频域窄带信号最优权矢量^并进行波束合成，计算得到N个通道中 的数字宽带信号在第j个频点处的最优波束形成方向图y〇 Pt( ; 步骤6,令j加1，返回步骤5,直到得到使用频域最小均方算法后第J个频点《 j处的频域 窄带信号最优权矢量Wj和N个通道中的数字宽带信号在第J个频点coj处的最优波束形成方 向图7_(〇^)，并将此时得到的1^个通道中的数字宽带信号在第1个频点《 1的最优波束形成 方向图ycipt( co 〇到N个通道中的数字宽带信号在第J个频点co j处的最优波束形成方向图y_ (co j)，g卩N个通道中的数字宽带信号在J个频点各自的最优波束形成方向发送至测向空间 谱计算单元。 步骤7,测向空间谱计算单元接收波束形成器发送过来的N个通道中的数字宽带信号在 J个频点各自的最优波束形成方向图，并进行几何平均，计算得到N个阵元接收到的时域宽 带信号对应的最终波束形成方向图F; 步骤8,根据N个阵元接收到的时域宽带信号对应的最终波束形成方向图F，计算得到使 用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数Plms，然后根据所述时域 宽带信号闭环测向空间谱函数Plms，获取时域宽带信号闭环测向空间谱函数Plms的谱峰位 置，进而得到时域宽带信号的角度估计值。5. 如权利要求4所述的一种时域宽带信号频域闭环测向方法，其特征在于，在步骤3中， 所述幅相误差校正后第m个通道中的数字宽带信号在第j个频点处的频域信号X m(〇^)，其表 达式为：其中，pG{l，2,…，P}，P表示第m个通道接收到的时域宽带信号包含的信号源数目，j G {1，2，…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 wh]内等效 为频域窄带信号叠加后包含的频点个数，SP(〇^)表示第p个信号源在第j个频点处的频域信 号具（^)表示第m个通道在第j个频点处的频域噪声，^表示第j个频点， Tpm表示第p个信 号源到达第m个通道时相对于参考通道的时延，参考通道为N个通道中第一个通道，exp( ?) 表示?的指数函数，mG{l，2,…，N}，N表示空间天线阵列包含的通道个数，且J、N分别为自 然数。6. 如权利要求4所述的一种时域宽带信号频域闭环测向方法，其特征在于，在步骤5中， 所述使用频域最小均方算法后第j个频点处的频域窄带信号最优权矢量％和所述N个通 道中的数字宽带信号在第j个频点处的最优波束形成方向图y〇 Pt(c^)，其得到过程为： 5.1初始化:n为迭代次数，且n初值为1，设定第j个频点处的频域最小均方算法调零 门限为Sj;m为空间天线阵列中第m个通道，且m初值为1，!11￡{1，2，"_，《4表示空间天线阵 列包含的通道个数，w q表示第1次迭代后第1个通道中的数字宽带信号在第j个频点处的 初始权值Wj(l，l)，即^(1，1)=?<^ 1(?」)表示第1个通道中的数字宽带信号在第」个频点处 的频域信号； 5.2计算得到第n次迭代后第m个通道中的数字宽带信号在第j个频点处的权值Wj(m， n)，其表达式为： wj(m,n+l) =wj(m,n)-2iiXm( ?j)e( to j) e( W j)=y( ? j)-d( 〇 j) 其中，d( CO」)表示第j个频点的频域期望信号，y( coj表示第j个频点的波束形成 方向图； 5.3令m加1，重复子步骤5.2，直到得到第n次迭代后第N个通道中的数字宽带信号在第j 个频点处的权值Wj(N，n)，此时计算第n次迭代后N个通道中的数字宽带信号在第j个频点 处的权值~(n)，并将所述第n次迭代后N个通道中的数字宽带信号在第j个频点处的 权值W(n)发送至波束形成器，计算得到第n次迭代后N个通道中的数字宽带信号在第j个频 点处的波束形成方向图h(n)，然后计算得到第n次迭代后第j个频点处的频域窄带信 号自相关' (/7): 5.4如果第n次迭代后第j个频点处的频域窄带信号自相关尺,(》)大于第j个频点 处的频域最小均方算法调零门限\，则令n加1，返回子步骤5.2; 如果第n次迭代后第j个频点处的频域窄带信号自相关小于第j个频点处的 频域最小均方算法调零门限则将此时得到的第n次迭代后N个通道中的数字宽带信号在 第j个频点处的权值W(n)，作为使用频域最小均方算法后得到的第j个频点处的频域 窄带信号最优权矢量将第n次迭代后N个通道中的数字宽带信号在第j个频点处的波 束形成方向图yXn)，作为N个通道中的数字宽带信号在第j个频点处的最优波束形成方 向图y〇 Pt( w」），Wj = [W1」，…，Wmj，…，WNj ]，)，_(仿,)=广)，X( ?j)表示N个通道中的数字 宽带信号在第j个频点处的频域信号，ww表示第m个通道中的数字宽带信号在第j个频点 处的权值。7. 如权利要求6所述的一种时域宽带信号频域闭环测向方法，其特征在于，所述第n次 迭代后N个通道中的数字宽带信号在第j个频点处的权值％( n)、第n次迭代后N个通道中 的数字宽带信号在第j个频点处的波束形成方向图h(n)和所述第n次迭代后第j个频点 %处的频域窄带信号自相关％(〃)，其表达式分别为： ffj(n) = [wj(l,n) ,---,wj(m,n) ,??? ,wj(N,n)] yj(n)=ffj(n)HX( 〇 j) /?, (") = ￡「沙)(")"义(似,),丫(％)" %(")] 其中，wXm，n)表示第n次迭代后第m个通道中的数字宽带信号在第j个频点处的权 值，X(c^)表示N个通道中的数字宽带信号在第j个频点处的频域信号， (%)，…，Xm( %)，…，XN( %) ]T，Xm( %)表示第m个通道中的数字宽带信号在第j个频点处 的频域信号。8. 如权利要求4所述的一种时域宽带信号频域闭环测向方法，其特征在于，在步骤7中， 所述幅相误差校正后第m个通道内的数字宽带信号对应的最终波束形成方向图F，其表达式其中，II表示联乘运算，y。^%)表示N个通道中的数字宽带信号在第j个频点^处的 最优波束形成方向图，W」表示使用频域最小均方算法后第j个频点处的频域窄带信号最 优权矢量，X」表示N个通道中的数字宽带信号在第j个频点处的频域信号x(c^)，je{l， 2，…，J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 Wh]内等效为 频域窄带信号叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信 号的频率下界， Wh表示幅相误差校正后第m个通道中的数字宽带信号的频率上界。9.如权利要求4所述的一种时域宽带信号频域闭环测向方法，其特征在于，在步骤7中， 所述使用频域最小均方算法后得到的时域宽带信号闭环测向空间谱函数Plms，其表达式为：其中，II表示联乘运算，y。^%)表示N个通道中的数字宽带信号在第j个频点^处的 最优波束形成方向图，W」表示使用频域最小均方算法后第j个频点处的频域窄带信号最 优权矢量，X(c〇J表示N个通道中的数字宽带信号在第j个频点处的频域信号，je{l，2，一， J}，J表示幅相误差校正后第m个通道内的数字宽带信号在宽带范围[W1 wh]内等效为频域窄 带信号叠加后包含的频点个数，wi表示幅相误差校正后第m个通道中的数字宽带信号的频 率下界，w h表示幅相误差校正后第m个通道中的数字宽带信号的频率上界。
【文档编号】G01S3/14GK106054122SQ201610370030
【公开日】2016年10月26日
【申请日】2016年5月30日
【发明人】陶海红, 张倩, 周忠锦, 刘宝瑞, 张丹丹, 梁克强, 廖桂生
【申请人】西安电子科技大学