CN106226829A - 一种地质结构面产状测定装置 - Google Patents

Abstract

本发明公开了一种地质结构面产状测定装置，所述底座通过铰链与上盖活动连接，所述上盖的内部中间位置设有反光镜，所述反光镜上设有透视孔，所述底座上设有水平刻度盘，所述水平刻度盘的内部设有磁针，所述磁针的一侧设有底盘水准器，所述磁针的另一侧设有水准器和刻度指示器，所述水准器的右上方设有激光测距器，所述激光测距器关于底座的中心对称位置设有磁针螺旋，所述底座的一端设有瞄准觇板。该发明设计科学合理，对地质结构面产状测试的可靠性稳定，测试的误差较少，提高了地质结构面产状测定的精准度，可多次重复操作，携带方便易于保存。

Description

一种地质结构面产状测定装置

技术领域

本发明属于地质测定技术领域，具体涉及一种地质结构面产状测定装置。

背景技术

随着我国经济的发展，自然资源的利用率越来越高，由于各地地质条件的差异，必须先进行地质岩石勘测才能进一步开发新的工程建设，为方便资源的开发和利用需要进驻一定的工具和手段进行科学探究来收集数据和资料。

地质勘查对于国民经济的发展具有重要的意义，属于一种基础性服务行业，对于我国经济较好较快发展具有重要的推动作用。改革开放至今，我国地质勘查行业取得了长足发展，很多勘查单位在项目、队伍建设等方面取得了骄人的成绩，但地质勘查业在发展过程中也遇到了很多瓶颈，这阻碍了地质勘查业的进一步发展，我国的地质勘查取得了重大成果。从地质勘查矿产资源的概况、勘查现状及其未来的发展趋势来看，当下，针对我国的地质结构面产状测定的测定设备，相对较少，而且有的设备携带不方便，有的设备测量误差较大，极大地影响着我的地质工程的顺利进展，因此，发明一种地质结构面产状测定装置，很有应用前景。

发明内容

本发明的目的在于提供一种地质结构面产状测定装置，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：一种地质结构面产状测定装置，包括底座，所述底座通过铰链与上盖活动连接，所述上盖的内部中间位置设有反光镜，所述反光镜上设有透视孔，所述底座上设有水平刻度盘，所述水平刻度盘的内部设有磁针，所述磁针的一侧设有底盘水准器，所述磁针的另一侧设有水准器和刻度指示器，所述水准器的右上方设有激光测距器，所述激光测距器关于底座的中心对称位置设有磁针螺旋，所述底座的一端设有瞄准觇板。

优选的，所述水准器包括垂直水准器和水平水准器。

优选的，所述透视孔设置为椭圆形。

优选的，所述瞄准觇板的中间位置设有瞄准孔。

优选的，所述激光测距器上设置有信号接收器，该信号接收器包括：

一接收器主体；

一波导管，该波导管包括一开口端以及一连接端，该连接端连接该接收器主体；

一传输单元，该传输单元设于该波导管之中，并靠近该连接端；以及

一介质片，该介质片设于该波导管之中，其中，该介质片包括一第一段部以及一第二段部，该第一段部靠近该传输单元，该第二段部靠近该开口端，该第一段部具有一第一介电常数，该第二段部具有一第二介电常数，该第一介电常数低于该第二介电常数。

优选的，该传输单元包括一传输针，该介质片与该传输针之间形成有一夹角，该夹角为55度；

该第一段部包括聚丙烯，该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物；

该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物以及聚丙烯；

该第一段部与该第二段部的长度比介于1∶0.8～1∶1之间；

该介质片包括一第一端以及一第二端，该第一端朝向该传输单元，该第二端朝向该波导管的该开口端，该第一端与该传输单元间形成有一第一间距，该第二端与该波导管的该开口端之间形成有一第二间距；

该第一间距为6.95～7.15mm，该第二间距为9.3～9.4mm；

该第一端上形成有一第一缺口，该第二端上形成有一第二缺口；

该波导管包括一导槽，该介质片沿该导槽滑动而插入于该波导管之中。

优选的，激光测距器上设置有同步正交跳频信号盲源分离模块，所述同步正交跳频信号盲源分离模块的同步正交跳频信号盲源分离方法包括以下步骤：

步骤一，利用含有M个阵元的阵列天线接收来自多个同步正交跳频电台的跳频信号，对每一路接收信号进行采样，得到采样后的M路离散时域混合信号m＝1,2,…,M；

步骤二，对M路离散时域混合信号进行重叠加窗短时傅里叶变换，得到M个混合信号的时频域矩阵

p＝0,1,…,P-1,q＝0,1,…,N_fft-1，其中P表示总的窗数，N_fft表示FFT变换长度；

步骤三，对步骤二中得到的跳频混合信号时频域矩阵进行预处理；对进行去低能量预处理，即在每一采样时刻p，将幅值小于门限ε的值置0，得到门限ε的设定可根据接收信号的平均能量来确定；

找出p时刻(p＝0,1,2,…P-1)非零的时频域数据，用表示，其中表示p时刻时频响应非0时对应的频率索引，对这些非零数据归一化预处理，得到预处理后的向量b(p,q)＝[b₁(p,q),b₂(p,q),…,b_M(p,q)]^T，其中

步骤四，利用聚类算法估计每一跳的跳变时刻以及各跳对应的归一化的混合矩阵列向量、跳频频率；

步骤五，根据步骤四估计得到的归一化混合矩阵列向量估计时频域跳频源信号；对所有采样时刻索引p判断该时刻索引属于哪一跳，具体方法为：如果则表示时刻p属于第l跳；如果则表示时刻p属于第1跳；对第l(l＝1,2,…)跳的所有时刻p_l，估计该跳各跳频源信号的时频域数据，计算公式如下：

$\left\{\begin{array}{c}{\tilde{S}}_j(p_l,q)=\frac{1}{\left|\right|{\hat{a}}_j\left(l\right)|{|}^2}\·{\hat{a}}_j^H\left(l\right)\×\left[\begin{array}{c}{\tilde{X}}_1(p_l,q)\\ {\tilde{X}}_2(p_l,q)\\ .\\ .\\ .\\ {\tilde{X}}_M(p_l,q)\end{array}\right]\\ j=\underset{j_0=1,2,\mathrm{...},\hat{N}}{\arg \max }\left(\right|{\[{\tilde{X}}_1\left(p_l,q\right),{\tilde{X}}_2\left(p_l,q\right),\mathrm{...},{\tilde{X}}_M\left(p_l,q\right)\]}^H\×{\hat{a}}_{j_0}\left(l\right)\left|\right)\\ {\tilde{S}}_m(p_l,q)=0,m=1,2,\mathrm{...},M,m\≠j\\ q=0,1,2,\mathrm{...},N_{fft}-1\end{array}\right.$

步骤六，对不同跳频点之间的时频域跳频源信号进行拼接；估计第l跳对应的个入射角度，用表示第l跳第n个源信号对应的入射角度，的计算公式如下：

${\widehat{\mathrm{\θ}}}_n\left(l\right)=\frac{1}{M-1}\underset{m=2}{\overset{M}{\Σ}}{\sin }^{-1}\[\frac{angle\left({\hat{a}}_{n,m}\right(l)/{\hat{a}}_{n,m-1}\left(l\right))*c}{2\mathrm{\π}{\hat{f}}_{c,n}\left(l\right)d}\],n=1,2,\mathrm{...},\hat{N}$

表示第l跳估计得到的第n个混合矩阵列向量的第m个元素，c表示光速，即v_c＝3×10⁸米/秒；判断第l(l＝2，3，…)跳估计的源信号与第一跳估计的源信号之间的对应关系，判断公式如下：

${m_n}^{\left(l\right)}=\underset{m}{\arg \min }|{\widehat{\mathrm{\θ}}}_m^{\left(l\right)}-{\widehat{\mathrm{\θ}}}_n^{\left(1\right)}|,n=1,2,\mathrm{...},\hat{N}$

其中m_n ^(l)表示第l跳估计的第m_n ^(l)个信号与第一跳估计的第n个信号属于同一个源信号；将不同跳频点估计到的属于同一个源信号的信号拼接在一起，作为最终的时频域源信号估计，用Y_n(p,q)表示第n个源信号在时频点(p,q)上的时频域估计值，p＝0,1,2,....,P，q＝0,1,2,...,N_fft-1，即

步骤七，根据源信号时频域估计值，恢复时域跳频源信号；对每一采样时刻p(p＝0，1，2，…)的频域数据Y_n(p,q),q＝0,1,2,…,N_fft-1做N_fft点的IFFT变换，得到p采样时刻对应的时域跳频源信号，用y_n(p,q_t)(q_t＝0,1,2,…,N_fft-1)表示；对上述所有时刻得到的时域跳频源信号y_n(p,q_t)进行合并处理，得到最终的时域跳频源信号估计，具体公式如下：

$s_n\&lsqb;kC:(k+1)C-1\&rsqb;=\left\{\begin{array}{cc}\underset{m=0}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k<K_c\\ \underset{m=k-K_c+1}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k\&GreaterEqual;K_c\end{array}\right.,k=0,1,2,\mathrm{...}$

这里K_c＝N_fft/C，C为短时傅里叶变换加窗间隔的采样点数，N_fft为FFT变换的长度。

本发明的技术效果和优点：该地质结构面产状测定装置，磁针的一侧的底盘水准器，对地质结构面产状测试的可靠性稳定，底座上的水平刻度盘，使得测试的误差较少，提高了地质结构面产状测定的精准度，该发明设计科学合理，对地质结构面产状测试的可靠性稳定，测试的误差较少，提高了地质结构面产状测定的精准度，可多次重复操作，携带方便易于保存。

附图说明

图1为本发明的结构示意图。

图中：1底座、2水准器、3激光测距器、4反光镜、5上盖、6透视孔、7水平刻度盘、8底盘水准器、9磁针螺旋、10瞄准觇板、11磁针、12刻度指示器。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提供了如图1所示的一种地质结构面产状测定装置，包括底座1，所述底座1通过铰链与上盖5活动连接，所述上盖5的内部中间位置设有反光镜4，所述反光镜4上设有透视孔6，所述底座1上设有水平刻度盘7，所述水平刻度盘7的内部设有磁针11，所述磁针11的一侧设有底盘水准器8，所述磁针11的另一侧设有水准器2和刻度指示器12，所述水准器2的右上方设有激光测距器3，所述激光测距器3关于底座1的中心对称位置设有磁针螺旋9，所述底座1的一端设有瞄准觇板10。

进一步地，所述水准器2包括垂直水准器和水平水准器。

进一步地，所述透视孔6设置为椭圆形。

进一步地，所述瞄准觇板10的中间位置设有瞄准孔。

优选的，所述激光测距器上设置有信号接收器，该信号接收器包括：

一接收器主体；

一波导管，该波导管包括一开口端以及一连接端，该连接端连接该接收器主体；

一传输单元，该传输单元设于该波导管之中，并靠近该连接端；以及

一介质片，该介质片设于该波导管之中，其中，该介质片包括一第一段部以及一第二段部，该第一段部靠近该传输单元，该第二段部靠近该开口端，该第一段部具有一第一介电常数，该第二段部具有一第二介电常数，该第一介电常数低于该第二介电常数。

优选的，该传输单元包括一传输针，该介质片与该传输针之间形成有一夹角，该夹角为55度；

该第一段部包括聚丙烯，该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物；

该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物以及聚丙烯；

该第一段部与该第二段部的长度比介于1∶0.8～1∶1之间；

该介质片包括一第一端以及一第二端，该第一端朝向该传输单元，该第二端朝向该波导管的该开口端，该第一端与该传输单元间形成有一第一间距，该第二端与该波导管的该开口端之间形成有一第二间距；

该第一间距为6.95～7.15mm，该第二间距为9.3～9.4mm；

该第一端上形成有一第一缺口，该第二端上形成有一第二缺口；

该波导管包括一导槽，该介质片沿该导槽滑动而插入于该波导管之中。

优选的，激光测距器上设置有同步正交跳频信号盲源分离模块，所述同步正交跳频信号盲源分离模块的同步正交跳频信号盲源分离方法包括以下步骤：

步骤一，利用含有M个阵元的阵列天线接收来自多个同步正交跳频电台的跳频信号，对每一路接收信号进行采样，得到采样后的M路离散时域混合信号m＝1,2,…,M；

步骤二，对M路离散时域混合信号进行重叠加窗短时傅里叶变换，得到M个混合信号的时频域矩阵

p＝0,1,…,P-1,q＝0,1,…,N_fft-1，其中P表示总的窗数，N_fft表示FFT变换长度；

步骤三，对步骤二中得到的跳频混合信号时频域矩阵进行预处理；对进行去低能量预处理，即在每一采样时刻p，将幅值小于门限ε的值置0，得到门限ε的设定可根据接收信号的平均能量来确定；

找出p时刻(p＝0,1,2,…P-1)非零的时频域数据，用表示，其中表示p时刻时频响应非0时对应的频率索引，对这些非零数据归一化预处理，得到预处理后的向量b(p,q)＝[b₁(p,q),b₂(p,q),…,b_M(p,q)]^T，其中

步骤四，利用聚类算法估计每一跳的跳变时刻以及各跳对应的归一化的混合矩阵列向量、跳频频率；

步骤五，根据步骤四估计得到的归一化混合矩阵列向量估计时频域跳频源信号；对所有采样时刻索引p判断该时刻索引属于哪一跳，具体方法为：如果则表示时刻p属于第l跳；如果则表示时刻p属于第1跳；对第l(l＝1,2,…)跳的所有时刻p_l，估计该跳各跳频源信号的时频域数据，计算公式如下：

$\left\{\begin{array}{c}{\tilde{S}}_j(p_l,q)=\frac{1}{\left|\right|{\hat{a}}_j\left(l\right)|{|}^2}\&CenterDot;{\hat{a}}_j^H\left(l\right)\&times;\left[\begin{array}{c}{\tilde{X}}_1(p_l,q)\\ {\tilde{X}}_2(p_l,q)\\ .\\ .\\ .\\ {\tilde{X}}_M(p_l,q)\end{array}\right]\\ j=\underset{j_0=1,2,\mathrm{...},\hat{N}}{\arg \max }\left(\right|{\&lsqb;{\tilde{X}}_1\left(p_l,q\right),{\tilde{X}}_2\left(p_l,q\right),\mathrm{...},{\tilde{X}}_M\left(p_l,q\right)\&rsqb;}^H\&times;{\hat{a}}_{j_0}\left(l\right)\left|\right)\\ {\tilde{S}}_m(p_l,q)=0,m=1,2,\mathrm{...},M,m\&NotEqual;j\\ q=0,1,2,\mathrm{...},N_{fft}-1\end{array}\right.$

步骤六，对不同跳频点之间的时频域跳频源信号进行拼接；估计第l跳对应的个入射角度，用表示第l跳第n个源信号对应的入射角度，的计算公式如下：

${\widehat{\mathrm{\&theta;}}}_n\left(l\right)=\frac{1}{M-1}\underset{m=2}{\overset{M}{\&Sigma;}}{\sin }^{-1}\&lsqb;\frac{angle\left({\hat{a}}_{n,m}\right(l)/{\hat{a}}_{n,m-1}\left(l\right))*c}{2\mathrm{\&pi;}{\hat{f}}_{c,n}\left(l\right)d}\&rsqb;,n=1,2,\mathrm{...},\hat{N}$

表示第l跳估计得到的第n个混合矩阵列向量的第m个元素，c表示光速，即v_c＝3×10⁸米/秒；判断第l(l＝2，3，…)跳估计的源信号与第一跳估计的源信号之间的对应关系，判断公式如下：

${m_n}^{\left(l\right)}=\underset{m}{\arg \min }|{\mathrm{\&theta;}}_m^{\left(l\right)}-{\widehat{\mathrm{\&theta;}}}_n^{\left(1\right)}|,n=1,2,\mathrm{...},\hat{N}$

其中m_n ^(l)表示第l跳估计的第m_n ^(l)个信号与第一跳估计的第n个信号属于同一个源信号；将不同跳频点估计到的属于同一个源信号的信号拼接在一起，作为最终的时频域源信号估计，用Y_n(p,q)表示第n个源信号在时频点(p,q)上的时频域估计值，p＝0,1,2,....,P，q＝0,1,2,...,N_fft-1，即

步骤七，根据源信号时频域估计值，恢复时域跳频源信号；对每一采样时刻p(p＝0，1，2，…)的频域数据Y_n(p,q),q＝0,1,2,…,N_fft-1做N_fft点的IFFT变换，得到p采样时刻对应的时域跳频源信号，用y_n(p,q_t)(q_t＝0,1,2,…,N_fft-1)表示；对上述所有时刻得到的时域跳频源信号y_n(p,q_t)进行合并处理，得到最终的时域跳频源信号估计，具体公式如下：

$s_n\&lsqb;kC:(k+1)C-1\&rsqb;=\left\{\begin{array}{cc}\underset{m=0}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k<K_c\\ \underset{m=k-K_c+1}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k\&GreaterEqual;K_c\end{array}\right.,k=0,1,2,\mathrm{...}$

这里K_c＝N_fft/C，C为短时傅里叶变换加窗间隔的采样点数，N_fft为FFT变换的长度。

该地质结构面产状测定装置，磁针11的一侧的底盘水准器8，对地质结构面产状测试的可靠性稳定，底座1上的水平刻度盘7，使得测试的误差较少，提高了地质结构面产状测定的精准度，该发明设计科学合理，对地质结构面产状测试的可靠性稳定，测试的误差较少，提高了地质结构面产状测定的精准度，可多次重复操作，携带方便易于保存。

最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种地质结构面产状测定装置，包括底座(1)，其特征在于：所述底座(1)通过铰链与上盖(5)活动连接，所述上盖(5)的内部中间位置设有反光镜(4)，所述反光镜(4)上设有透视孔(6)，所述底座(1)上设有水平刻度盘(7)，所述水平刻度盘(7)的内部设有磁针(11)，所述磁针(11)的一侧设有底盘水准器(8)，所述磁针(11)的另一侧设有水准器(2)和刻度指示器(12)，所述水准器(2)的右上方设有激光测距器(3)，所述激光测距器(3)关于底座(1)的中心对称位置设有磁针螺旋(9)，所述底座(1)的一端设有瞄准觇板(10)；

所述水准器(2)包括垂直水准器和水平水准器；

所述透视孔(6)设置为椭圆形；

所述瞄准觇板(10)的中间位置设有瞄准孔。

2.如权利要求1所述的地质结构面产状测定装置，其特征在于，所述激光测距器上设置有信号接收器，该信号接收器包括：

一接收器主体；

一波导管，该波导管包括一开口端以及一连接端，该连接端连接该接收器主体；

一传输单元，该传输单元设于该波导管之中，并靠近该连接端；以及

一介质片，该介质片设于该波导管之中，其中，该介质片包括一第一段部以及一第二段部，该第一段部靠近该传输单元，该第二段部靠近该开口端，该第一段部具有一第一介电常数，该第二段部具有一第二介电常数，该第一介电常数低于该第二介电常数。

3.如权利要求2所述的地质结构面产状测定装置，其特征在于，该传输单元包括一传输针，该介质片与该传输针之间形成有一夹角，该夹角为55度；

该第一段部包括聚丙烯，该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物；

该第二段部包括丙烯腈-丁二烯-苯乙烯共聚物以及聚丙烯；

该第一段部与该第二段部的长度比介于1∶0.8～1∶1之间；

该介质片包括一第一端以及一第二端，该第一端朝向该传输单元，该第二端朝向该波导管的该开口端，该第一端与该传输单元间形成有一第一间距，该第二端与该波导管的该开口端之间形成有一第二间距；

该第一间距为6.95～7.15mm，该第二间距为9.3～9.4mm；

该第一端上形成有一第一缺口，该第二端上形成有一第二缺口；

该波导管包括一导槽，该介质片沿该导槽滑动而插入于该波导管之中。

4.如权利要求1所述的地质结构面产状测定装置，其特征在于，激光测距器上设置有同步正交跳频信号盲源分离模块，所述同步正交跳频信号盲源分离模块的同步正交跳频信号盲源分离方法包括以下步骤：

步骤一，利用含有M个阵元的阵列天线接收来自多个同步正交跳频电台的跳频信号，对每一路接收信号进行采样，得到采样后的M路离散时域混合信号m＝1,2,…,M；

步骤二，对M路离散时域混合信号进行重叠加窗短时傅里叶变换，得到M个混合信号的时频域矩阵p＝0,1,…,P-1,q＝0,1,…,N_fft-1，其中P表示总的窗数，N_fft表示FFT变换长度；

步骤三，对步骤二中得到的跳频混合信号时频域矩阵进行预处理；对进行去低能量预处理，即在每一采样时刻p，将幅值小于门限ε的值置0，得到门限ε的设定可根据接收信号的平均能量来确定；

找出p时刻(p＝0,1,2,…P-1)非零的时频域数据，用表示，其中表示p时刻时频响应非0时对应的频率索引，对这些非零数据归一化预处理，得到预处理后的向量b(p,q)＝[b₁(p,q),b₂(p,q),…,b_M(p,q)]^T，其中

步骤四，利用聚类算法估计每一跳的跳变时刻以及各跳对应的归一化的混合矩阵列向量、跳频频率；

步骤五，根据步骤四估计得到的归一化混合矩阵列向量估计时频域跳频源信号；对所有采样时刻索引p判断该时刻索引属于哪一跳，具体方法为：如果则表示时刻p属于第l跳；如果则表示时刻p属于第1跳；对第l(l＝1,2,…)跳的所有时刻p_l，估计该跳各跳频源信号的时频域数据，计算公式如下：

$\left\{\begin{array}{c}{\tilde{S}}_j(p_l,q)=\frac{1}{\left|\right|{\hat{a}}_j\left(l\right)|{|}^2}\&CenterDot;{\hat{a}}_j^H\left(l\right)\&times;\left[\begin{array}{c}{\tilde{X}}_1(p_l,q)\\ {\tilde{X}}_2(p_l,q)\\ .\\ .\\ .\\ {\tilde{X}}_M(p_l,q)\end{array}\right]\\ j=\underset{j_0=1,2,\mathrm{...},\hat{N}}{arg\max }\left(\right|{\&lsqb;{\tilde{X}}_1\left(p_l,q\right),{\tilde{X}}_2\left(p_l,q\right),\mathrm{...},{\tilde{X}}_M\left(p_l,q\right)\&rsqb;}^H\&times;{\hat{a}}_{j_0}\left(l\right)\left|\right)\\ {\tilde{S}}_m(p_l,q)=0,m=1,2,\mathrm{...},M,m\&NotEqual;j\\ q=0,1,2,\mathrm{...},N_{fft}-1\end{array}\right.$

步骤六，对不同跳频点之间的时频域跳频源信号进行拼接；估计第l跳对应的个入射角度，用表示第l跳第n个源信号对应的入射角度，的计算公式如下：

${\widehat{\mathrm{\&theta;}}}_n\left(l\right)=\frac{1}{M-1}\underset{m=2}{\overset{M}{\&Sigma;}}{\sin }^{-1}\&lsqb;\frac{angle\left({\hat{a}}_{n,m}\right(l)/{\hat{a}}_{n,m-1}\left(l\right))*c}{2\mathrm{\&pi;}{\hat{f}}_{c,n}\left(l\right)d}\&rsqb;,n=1,2,\mathrm{...},\hat{N}$

表示第l跳估计得到的第n个混合矩阵列向量的第m个元素，c表示光速，即v_c＝3×10⁸米/秒；判断第l(l＝2，3，…)跳估计的源信号与第一跳估计的源信号之间的对应关系，判断公式如下：

${m_n}^{\left(l\right)}=\underset{m}{\arg \min }|{\widehat{\mathrm{\&theta;}}}_m^{\left(l\right)}-{\widehat{\mathrm{\&theta;}}}_n^{\left(1\right)}|,n=1,2,...,\hat{N}$

其中m_n ^(l)表示第l跳估计的第m_n ^(l)个信号与第一跳估计的第n个信号属于同一个源信号；将不同跳频点估计到的属于同一个源信号的信号拼接在一起，作为最终的时频域源信号估计，用Y_n(p,q)表示第n个源信号在时频点(p,q)上的时频域估计值，p＝0,1,2,....,P，q＝0,1,2,...,N_fft-1，即

步骤七，根据源信号时频域估计值，恢复时域跳频源信号；对每一采样时刻p(p＝0，1，2，…)的频域数据Y_n(p,q),q＝0,1,2,…,N_fft-1做N_fft点的IFFT变换，得到p采样时刻对应的时域跳频源信号，用y_n(p,q_t)(q_t＝0,1,2,…,N_fft-1)表示；对上述所有时刻得到的时域跳频源信号y_n(p,q_t)进行合并处理，得到最终的时域跳频源信号估计，具体公式如下：

$s_n\&lsqb;kC:(k+1)C-1\&rsqb;=\left\{\begin{array}{cc}\underset{m=0}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k<K_c\\ \underset{m=k-K_c+1}{\overset{k}{\&Sigma;}}{y}_n\&lsqb;m,(k-m)C:(k-m+1)C-1\&rsqb;& k\&GreaterEqual;K_c\end{array}\right.,k=0,1,2,\mathrm{...}$

这里K_c＝N_fft/C，C为短时傅里叶变换加窗间隔的采样点数，N_fft为FFT变换的长度。