CN102298097B - 一种估计云闪雷电脉冲信号到达时间差的方法 - Google Patents
一种估计云闪雷电脉冲信号到达时间差的方法 Download PDFInfo
- Publication number
- CN102298097B CN102298097B CN 201110199389 CN201110199389A CN102298097B CN 102298097 B CN102298097 B CN 102298097B CN 201110199389 CN201110199389 CN 201110199389 CN 201110199389 A CN201110199389 A CN 201110199389A CN 102298097 B CN102298097 B CN 102298097B
- Authority
- CN
- China
- Prior art keywords
- pulse
- station
- lightning
- homologous
- pulses
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000001514 detection method Methods 0.000 claims abstract description 66
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims description 37
- 238000005314 correlation function Methods 0.000 claims description 26
- 230000000875 corresponding effect Effects 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 claims description 5
- 230000002596 correlated effect Effects 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 238000003491 array Methods 0.000 abstract 4
- 230000005855 radiation Effects 0.000 description 7
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 6
- 239000000523 sample Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 241000167686 Reichardia Species 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种估计云闪雷电脉冲信号到达时间差的方法,步骤为:雷电脉冲数据块拼接、波形还原、形成雷电脉冲簇;输入时间区间,提取各站的雷电脉冲簇;确定参考探测站;统计A站雷电脉冲簇中的脉冲个数;雷电脉冲关联与TDOA计算;形成用户指定时间区间内的所有同源脉冲对队列和TDOA队列;改变时间区间,计算并形成新的同源脉冲对队列和TDOA队列。本发明专门针对云闪雷电脉冲信号的特点进行设计,给出了波形关联与同源性判定、到达时间差估计的完整过程,可胜任白噪声和色噪声环境下的TDOA估计要求。本发明方法无需遍历两个探测站接收到的所有雷电脉冲,鲁棒性强,简单可靠,计算结果精度高,易于工程实现。
Description
技术领域
本发明属于雷电监测预警技术研究领域,涉及一种在由多个探测站和一个中心站组成的云闪雷电探测定位***中,在中心站精确计算雷电脉冲信号到达任意两个探测站的时间差(TDOA)的方法。
背景技术
基于信号来波方向(DOA)对雷电辐射源进行定位,其定位精度与空间分辨力将随着雷电辐射源与探测站之间距离的增大而降低。相比之下,基于信号到达时间差(TimeDifference Of Arrival,TDOA)的雷电辐射源定位技术,其定位精度与空间分辨力则将随着辐射源与探测站间距的增加而得到提高(当定位双曲面相互正交时为最佳),但当辐射源非常靠近某个探测站时,TDOA定位的空间分辨力也将急剧下降。
获取辐射源发射信号到达两个探测站的时间差有两种方法:(1)将两个接收站得到的信号到达时间(TOA)测量值作差;(2)利用互相关技术,计算一个探测站的接收信号与另一个探测站的接收信号之间的互相关函数,然后进行峰值搜索。前一种方法不仅要求各探测站之间保持严格的时钟同步,而且还须知道发射信号波形的先验知识,因而多用于规则脉冲信号情形。采用互相关技术直接计算信号到达时间差,可以完全摆脱对发射信号波形的限制,可用于对雷电脉冲到达时间差的估计,其应用的前提条件依然是各探测站之间必须保持高精度的时钟同步。
目前常用的时延估计方法有:广义相关时延估计法,如以P.R.Roth命名的ROTH加权广义相关法,以J.C.Hassab和R.E.Boucher命名的HB加权法;用三阶累积量或三阶矩的傅立叶变换的广义双谱时延估计法;利用功率谱与相关函数互为傅立叶变换的广义相位谱时延估计法;以自适应滤波或自适应迭代为特征的自适应时延估计法;基于子空间分解的超分辨时延估计法。在广义相关时延估计法中,在对两个探测站的信号进行相关处理前,先对其进行预滤波,相当于在频域进行加权,有利于增强接收信号的谱分量,提高信噪比,从而获得较高的时延估计精度。这些方法若直接用来求取雷电脉冲的TDOA,其缺点是需要遍历各探测站接收到的所有雷电脉冲信号,运算量非常大,且由于相关峰判决门限不易确定,导致出现许多虚假估计值。
目前,尚未发现对云闪雷电辐射源进行定位时,具备高精度、鲁棒性雷电脉冲信号到达时间差估计的完整方法和实现流程。
发明内容
本发明的目的在于提供一种估计云闪雷电脉冲信号到达时间差的方法,该方法无需遍历两个探测站接收到的所有雷电脉冲,其计算量和虚假估计结果明显减少,估计精度得到提高,且同时也获得了配对的TDOA信息。
本发明提供一种估计云闪雷电脉冲信号到达时间差的方法,该方法包括下述步骤:
第1步对上传到中心站的雷电波形数据块进行拼接和波形还原,得到一个个完整的雷电脉冲,形成雷电脉冲簇;
第2步中心站根据设定的时间区间提取位于该时间区间内的各探测站雷电脉冲簇;
第3步针对用户输入时间区间内、各探测站的雷电脉冲簇,中心站统计分属于每个探测站的雷电脉冲簇中的脉冲个数,由此计算出这些脉冲所对应的原始雷电波形数据块的总个数,并将该总个数最多的那个探测站作为该时间区间内进行到达时间差计算的参考探测站,命名为A站,其它探测站称之为一般探测站,命名为B站、C站、D站、......,设一般探测站的数量为S;
第4步对所述时间区间内参考探测站雷电脉冲簇中的脉冲个数Q进行统计;
第5步令q=1,s=1;
第6步取参考探测站雷电脉冲簇中的第q个脉冲,与第s个一般探测站的脉冲簇中的各个脉冲进行粗关联、细关联,得到同源脉冲配对结果,再计算同源脉冲之间的到达时间差值,最后形成完整的同源脉冲配对及其对应的到达时间差队列;
第7步令s=s+1,判断s>S是否成立,如果成立,转入第8步,否则进入第6步;
第8步令q=q+1,判断q>Q是否成立,如果成立,转入第9步,否则令s=1,转入第6步;
第9步改变时间区间,重复前述步骤(2)~(8),得到位于该时间区间内的同源脉冲对队列和到达时间差队列。
本发明方法鲁棒地估计云闪雷电脉冲信号到达任意两个探测站的时间差,它通过若干预处理和判断过程缩小可能关联的雷电脉冲信号范围直至完成雷电脉冲信号的一一配对(得到同源脉冲对),再利用广义互相关函数或脉冲的包络能量重心来最终估计出TDOA,这样可以大幅度减小TDOA估计的计算量及可能出现的虚假结果,并依据云闪雷电脉冲的持续时间长度选择不同的TDOA具体计算方法。具体而言,本发明具有以下特点:
(1)采用两级关联进行脉冲配对,方法简单,运算量低,遗漏或虚假配对的情况大幅度减小,方法鲁棒性强;
(2)依据云闪雷电脉冲信号的持续时间长短来选择采用不同的TDOA估计方法,兼顾了计算精度、计算量和结果可信度;
(3)通过在互相关函数峰值点附近进行插值,在运算量增加较少的前提下,使TDOA估计的精度突破了整数倍波形采样间隔的限制,实现了连续时差估计;
(4)本发明方法可应用于低信噪比、高信噪比环境,可应用于白噪声、色噪声环境,易于工程实现。
本发明专门针对云闪雷电脉冲信号的特点进行设计,给出了波形关联与同源性判定、到达时间差估计的完整过程,可胜任白噪声和色噪声环境下的TDOA估计要求。本发明所指方法,鲁棒性强,简单可靠,计算结果精度高,易于工程实现。
附图说明
图1是一种估计云闪雷电脉冲到达时间差的计算流程;
图2是雷电发生源与两个探测站的空间几何关系;
图3是雷电脉冲信号的包络能量重心示意图;
图4是当两个同源脉冲的长度相同时,利用包络能量重心计算TDOA的原理示意图(情况一);
图5是当两个同源脉冲的长度相同时,利用包络能量重心计算TDOA的原理示意图(情况二);
图6是当两个同源脉冲的长度不相同时,利用包络能量重心计算TDOA的原理示意图。
具体实施方式
以下结合附图对本发明的内容及具体实施方式作进一步说明。
图1给出了一种在由多个探测站和一个中心站组成的云闪雷电定位***中,形成同源脉冲队列并计算雷电信号到达任意两个探测站的时间差的方法流程,图2给出了在直角坐标系中雷电发生源与两探测站的空间几何关系,图3给出了雷电信号包络能量重心的定义,图4、图5和图6分别给出了依据两个同源脉冲的包络能量重心计算信号到达时间差(TDOA)的原理示意图。
如图1所示,本发明提供一种鲁棒的云闪雷电脉冲到达时间差的估计方法,其步骤包括:
(1)雷电脉冲数据块拼接、波形还原、形成雷电脉冲簇
为了利用各个探测站上传的雷电波形数据块进行TDOA估计,首先需对上传到中心站的雷电波形数据块进行拼接和波形还原,得到一个个完整的雷电脉冲,形成雷电脉冲簇。
以同一个探测站上传到中心站的各个雷电波形数据块为例进行说明。首先,查询该探测站上传到中心站的雷电波形数据块的起始时间、终止时间标记,若前后两个数据块的终止时间、起始时间是连续的,则将这两个数据块拼接为一个较长的数据块。依次类推,完成该探测站上传的所有数据块的起止时间查询和连续数据块的拼接。若前后两个数据块的终止时间、起始时间是不连续的,则停止本个连续数据块的拼接,重新开始下一个新数据块的拼接过程。
对各个探测站上传的所有雷电波形数据块,按上述方式完成各自全部数据块的拼接处理。
对上述拼接得到的每个数据块内的数据,按下式进行波形还原,形成雷电脉冲信号数据块:
式中,n表示脉冲波形采样点的采样时刻,f0表示接收信号的中心频率,fs表示采样率,I(n)和Q(n)分别表示雷电信号的同相分量和正交分量,x(n)表示还原出的雷电波形在时刻n时的采样值。
在每个还原出的雷电脉冲信号数据块前部标记以下信息:所属探测站的序号、探测站接收到该脉冲的起止时刻、探测站开始向中心站上传该脉冲数据的时刻、中心站接收完该脉冲数据的时刻。
对应于同一个探测站的全部雷电脉冲,形成该站的雷电脉冲簇。
(2)输入时间区间,提取各站的雷电脉冲簇
中心站以用户输入的时间区间为依据,取出各探测站雷电脉冲簇中脉冲起始时刻位于该时间区间内的雷电脉冲簇数据。
(3)确定参考探测站
针对用户输入时间区间内、各探测站的雷电脉冲簇,中心站统计分属于每个探测站的雷电脉冲簇中的脉冲个数,由此计算出这些脉冲所对应的原始雷电波形数据块的总个数,并将该总个数最多的那个探测站作为该时间区间内进行TDOA计算的参考探测站,将其命名为A站,而其它探测站命名为B站、C站、D站、......,以下计算过程均以B站为例。
注意到一个完整的云闪雷电过程的持续时间从数十微秒至数百毫秒不等,当相邻探测站的间距不超过200公里时,用户输入时间区间的基本原则是:起止时刻差不超过一整秒。
(4)统计A站雷电脉冲簇中的脉冲个数
确定了参考探测站即A站后,对该时间区间内A站雷电脉冲簇中的脉冲个数Q进行统计。
(5)雷电脉冲关联与TDOA计算
依次从A站雷电脉冲簇中取出第一个、第二个、......、第Q个脉冲,与B站、C站、D站、......等各个站的脉冲簇中的各个脉冲进行粗关联、细关联,得到同源脉冲配对结果,再计算同源脉冲之间的TDOA值,最后形成完整的同源脉冲配对及其对应的TDOA队列。以下以A站第q个脉冲与B站的脉冲簇进行同源性配对和TDOA计算过程为例,说明具体实施步骤。
(5.1)雷电脉冲簇的同源性粗关联
以A站的第q个脉冲为基准,将其与B站的脉冲簇进行粗关联的模块。此处所谓的粗关联,是指对两个雷电脉冲的波形相似性进行初步判断。
如图2,假设τS时刻在空间S点处发生雷电,其辐射的雷电脉冲分别于τA和τB时刻抵达A、B两站。则必有:
|(τA-τS)·Vc-(τB-τS)·Vc|≤DAB;
即:|(τA-τB)·Vc|≤DAB
或:
其中,DAB为A站和B站之间的直线距离,Vc为雷电波传播速度。
据此可以对探测站A、B接收到的雷电脉冲簇进行同源性粗关联:以A站脉冲簇中的第q个脉冲起始时标τA为基点,将B站脉冲簇中脉冲起始时标位于时间区间内的所有脉冲粗关联为A站该脉冲的潜在同源脉冲,形成B站潜在同源脉冲簇。
为了减少上述粗关联判决引起的遗漏情况、将B站尽可能多的脉冲纳入到粗关联潜在同源脉冲簇中,可将上述粗关联判决公式修正为:
|(τA-τS)·Vc-(τB-τS)·Vc|≤DAB+Vc·max{TA,TB}
即:
即:以A站脉冲簇中的第q个脉冲起始时标τA为基点,将B站脉冲簇中脉冲起始时标位于时间区间内的所有脉冲粗关联为A站该脉冲的潜在同源脉冲,形成B站潜在同源脉冲簇。式中,TA、TB分别为进行粗关联的A站、B站雷电脉冲的持续时间。
(5.2)雷电脉冲的同源性细关联
A站第q个脉冲与B站的潜在同源脉冲簇进行细关。此处所谓的细关联,是指对两个雷电脉冲的波形相似性进行最终判断。
在完成前述的雷电脉冲粗关联得到潜在同源脉冲簇后,再进行同源性细关联,实现A站与其它探测站雷电脉冲的一一配对(形成同源脉冲对)。以与B站潜在同源脉冲簇的配对为例,说明本发明所指的三种细关联方法:
①从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲之间的包络互相关系数ρi,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预先设定的互相关系数判决门限进行比较或直接按照取大原则,找出与A站该脉冲真正同源的B站脉冲,形成A、B站同源脉冲对。
其中,雷电脉冲信号的包络为:
A站与B站雷电脉冲信号的包络互相关系数按下式计算:
式中, I为B站潜在同源脉冲簇队列中的脉冲个数,yA(n)和yBi(n)分别为A、B站的雷电脉冲包络,其采样点数分别为MA和MBi。当MA≠MBi时(即yA(n)和yBi(n)两个脉冲的长度不相等时),须对较短的脉冲波形补零直至与另一脉冲波形的长度相等。既可在脉冲前端补零,也可在脉冲后端补零。因此,N=max{MA,MBi}。
为了在不同的信噪比条件下,兼顾细关联的虚警概率和检测概率,可对计算出的全部互相关系数ρ1,ρ2,…,ρI按大小排序,并按高于判决门限原则或取大原则筛选出同源脉冲。其中,互相关系数判决门限设定为0.5~0.9之间的一个数值。
②从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲的波形包络归一化均方偏差Ji,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预先设定的均方偏差判决门限进行比较或直接按照取小原则,找出与A站该脉冲真正同源的B站脉冲,形成同源脉冲对。
其中,雷电脉冲的波形包络归一化均方偏差为:
式中,GABi为B站脉冲包络相对于A站脉冲包络峰值的归一化系数,其值等于A站脉冲包络的最大值与B站脉冲包络的最大值之比,I为B站潜在同源脉冲簇队列中的脉冲个数。当MA≠MBi时(即yA(n)和yBi(n)两个脉冲的长度不相等时),须对较短的脉冲波形补零直至与另一脉冲的长度相等,补零的基本原则为:在较短脉冲的前、后两端同时补零。因此,N=max{MA,MBi}。
为了在不同的信噪比条件下,兼顾细关联的虚警概率和检测概率,可对计算出的全部均方偏差J1,J2,…,JI按大小排序,并按低于判决门限原则或取小原则筛选出同源脉冲。其中,均方偏差判决门限设定为0~1.0之间的一个数值。
③从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲的波形包络偏度偏差Si或(和)波形包络峰度偏差Ki,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预先设定的偏度偏差判决门限或(和)峰度偏差判决门限进行比较,找出与A站该脉冲真正同源的B站脉冲,形成同源脉冲对。其中,偏度偏差判决门限和峰度偏差判决门限均可设定为0~1.0之间的一个数值。
其中,雷电脉冲的波形包络偏度偏差和波形包络峰度偏差分别为:
式中,E(·)为计算均值。I为B站潜在同源脉冲簇队列中的脉冲个数。当MA≠MBi时(即yA(n)和yBi(n)两个脉冲的长度不相等时),须对较短的脉冲波形补零直至与另一脉冲的长度相等。既可在脉冲前端补零,也可在脉冲后端补零。因此,N=max{MA,MBi}。
为了在不同的信噪比条件下,兼顾细关联的虚警概率和检测概率,可对计算出的全部偏度偏差和峰度偏差S1,S2,…,SI和K1,K2,…,KI按大小排序,并按取小原则筛选出同源脉冲。
(5.3)对同源脉冲对进行TDOA估计
由于不同时间和空间点发生的雷电脉冲,其持续时间可能相差很大,因此本发明给出了两类计算TDOA的方法。以下以A、B站的同源脉冲对为例具体说明这两类计算方法。
①短脉冲情况下的TDOA计算
当同源脉冲对中的脉冲持续时间(脉冲长度)较短时(例如,N=200或400),按下式分别计算同源脉冲对中的A站、B站脉冲的包络能量重心:
图3说明了本发明所指的脉冲的包络能量重心的含义。将脉冲包络的时间中心点定为参考零点,其前、后段设置为相同采样点数,则能量重心为相对于该参考零点的偏移量。图4、图5和图6表明了根据两个脉冲的起始时刻和重心偏差计算脉冲到达A、B站的时间差的原理。
所以,上述同源脉冲对中的两个脉冲的能量重心实际时刻为:
因此,该雷电脉冲到达A站、B站的时间差为:
式中,MA=2NA为同源脉冲对中的A站脉冲长度(采样点数),MBi=2NBi为同源脉冲对中的B站脉冲长度(采样点数),T0为采样点间隔(T0=1/fs,fs为采样频率),tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻。
②长脉冲情况下的TDOA计算
当同源脉冲对中的脉冲持续时间(脉冲长度)较长时(例如,N=600或800或1000等情形),采用脉冲信号的互相关函数峰值点的位置作为雷电信号到达A站和B站的时间差(TDOA)。本发明所指的长脉冲情况下的TDOA估计方法,包括以下三种:
(a)对于同源脉冲对中的A站与B站脉冲,按下式计算它们的互相关函数:
由此得到该雷电脉冲到达A站和B站的时间差为:
式中,M为该段雷电信号的互相关函数的峰值点个数,对应于多源或多径雷电信号的个数;T0为采样点间隔(T0=1/fs,fs为采样频率);tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻;表示的绝对值取极大值时所对应的自变量τ之值。
当A站与B站脉冲长度不相等且分别为NA和NBi时,须首先将两个脉冲序列补零至长度均为NA+NBi-1,再计算互相关函数。对脉冲序列补零的原则:在原脉冲序列的后端补零。
由此计算最大似然加权函数:
其中,
进而计算广义互相关函数:
由此得到该雷电脉冲到达A站和B站的时间差为:
式中,M为该段雷电信号的互相关函数的峰值点个数,对应于多源或多径雷电脉冲的个数;T0为采样点间隔(T0=1/fs,fs为采样频率);tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻;IDFT[.]表示你傅里叶变换;
当A站与B站脉冲长度不相等且分别为NA和NBi时,须首先将两个脉冲序列补零至长度均为NA+NBi-1,再计算互相关函数。对脉冲序列补零的原则:在原脉冲序列的后端补零。
(c)为了进一步提高TDOA估计的精度,在对前述计算出的互相关函数进行峰值搜索时,对得到的峰值点位置进行次极大值点线性插值和相邻三点抛物线插值后再计算新的峰值点位置,并以此作为TDOA估计值输出。
设对互相关函数进行搜索得到的峰值点为dm,其互相关函数值为dm点的前、后相邻点为dm-1和dm+1;dm点的前、后相邻次极大值点为dm-k1和dm+k2,其互相关函数值分别为k1和k2为前、后相邻次极大值点到峰值点dm的距离。
对dm-k1点、dm点、dm+k2点进行线性插值,并将插值结果分别赋予dm-1点和dm+1点:
式中,M为该段雷电信号中包含的多源或多径雷电信号的个数,T0为采样点间隔(T0=1/fs,fs为采样频率);tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻。
(5.4)形成同源脉冲对队列、同源TDOA队列,计算雷电信号到达任意两个探测站的时间差。
按照前述各步骤,得到A站与B站同源脉冲对和TDOA后,将其加入到同源脉冲对队列和TDOA队列中;然后按照相同步骤得到A站与C站的同源脉冲对和TDOA,并将其加入到同源脉冲对队列和TDOA队列中,得到A站与D站的同源脉冲对和TDOA,并将其加入到同源脉冲对队列和TDOA队列中;如此重复,直至最后形成完整的同源脉冲对队列及对应得TDOA队列。
在上述所有站同源脉冲配对队列中,任取两站对应的雷电脉冲,按照前述步骤(5.3)给出的TDOA计算方法,得到雷电信号到达任意两个探测站的时间差。
(6)形成用户指定时间区间内的所有同源脉冲对队列和TDOA队列
重复前述步骤(5),在用户指定时间区间内,关联出A站的每一个脉冲与B站、C站、D站、......的同源脉冲对及其中任意两个脉冲的到达时间差TDOA,形成用户指定时间区间内的所有同源脉冲对队列和TDOA队列。
(7)改变时间区间,计算并形成新的同源脉冲对队列和TDOA队列
重复前述步骤(2)~(6),根据用户新输入的时间区间,重新确定参考探测站,然后进行同源脉冲关联和TDOA计算,最后得到该时间区间内的同源脉冲对队列和TDOA队列。
本发明不仅局限于上述具体实施方式,本领域一般技术人员根据本发明公开的内容,可以采用其它多种具体实施方式实施本发明,因此,凡是采用本发明的设计结构和思路,做一些简单的变化或更改的设计,都落入本发明保护的范围。
Claims (8)
1.一种估计云闪雷电脉冲信号到达时间差的方法,该方法包括下述步骤:
第1步对上传到中心站的雷电波形数据块进行拼接和波形还原,得到一个个完整的雷电脉冲,形成雷电脉冲簇;
第2步中心站根据设定的时间区间提取位于该时间区间内的各探测站雷电脉冲簇;
第3步针对用户输入时间区间内、各探测站的雷电脉冲簇,中心站统计分属于每个探测站的雷电脉冲簇中的脉冲个数,由此计算出这些脉冲所对应的原始雷电波形数据块的总个数,并将该总个数最多的那个探测站作为该时间区间内进行到达时间差计算的参考探测站,命名为A站,其它探测站称之为一般探测站,命名为B站、C站、D站、……,设一般探测站的数量为S;
第4步对所述时间区间内参考探测站雷电脉冲簇中的脉冲个数Q进行统计;
第5步令q=1,s=1;
第6步取参考探测站雷电脉冲簇中的第q个脉冲,与第s个一般探测站的脉冲簇中的各个脉冲进行粗关联、细关联,得到同源脉冲配对结果,再计算同源脉冲之间的到达时间差值,最后形成完整的同源脉冲配对及其对应的到达时间差队列;
第7步令s=s+1,判断s>S是否成立,如果成立,转入第8步,否则进入第6步;
第8步令q=q+1,判断q>Q是否成立,如果成立,转入第9步,否则令s=1,转入第6步;
第9步改变时间区间,重复前述步骤(2)~(8),得到位于该时间区间内的同源脉冲对队列和到达时间差队列。
2.根据权利要求1所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第1步中,按照下述方式对雷电波形数据块进行拼接:查询每个探测站上传到中心站的雷电波形数据块的起始时间、终止时间标记,若前后两个数据块的终止时间、起始时间是连续的,则将这两个数据块拼接为一个较长的数据块;依次类推,完成该探测站上传的所有数据块的起止时间查询和连续数据块的拼接;若前后两个数据块的终止时间、起始时间是不连续的,则停止对本次连续数据块的拼接,重新开始下一个新数据块的拼接过程。
4.根据权利要求1所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第6步中,粗关联是指对两个雷电脉冲的波形相似性进行初步判断,
假设τS时刻在空间S点处发生雷电,其辐射的雷电脉冲分别于τA和τB时刻抵达A、B两站,粗关联判决公式为:
|(τA-τS)·Vc-(τB-τS)·Vc|≤DAB+Vc·max{TA,TB}
即:
其中,DAB为A站和B站之间的直线距离,Vc为雷电信号传播速度,TA、TB分别为进行粗关联的A站、B站雷电脉冲的持续时间;
以A站脉冲簇中的第q个脉冲起始时标τA为基点,将B站脉冲簇中脉冲起始时标位于时间区间 内的所有脉冲粗关联为A站该脉冲的潜在同源脉冲,形成B站潜在同源脉冲簇。
5.根据权利要求4所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第6步中,细关联是指对两个雷电脉冲的波形相似性进行最终判断,其过程为:
从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲的包络互相关系数ρi,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预先设定的互相关系数判决门限进行比较或直接按照取大原则,找出与A站该脉冲真正同源的B站脉冲,形成A、B站同源脉冲对;
其中,雷电脉冲信号的包络为:
I(n)和Q(n)分别表示雷电信号的同相分量和正交分量;
A站与B站雷电脉冲信号的包络互相关系数按下式计算:
式中, I为B站潜在同源脉冲簇队列中的脉冲个数,yA(n)和yBi(n)为A、B站在采样时刻n时的雷电脉冲包络波形,其采样点数分别为MA和MBi,当MA≠MBi时,须对较短的脉冲波形补零直至与另一脉冲的长度相等,得到N=max{MA,MBi};对计算出的全部互相关系数ρ1 ,ρ2,…,ρI按大小排序,并按过门限原则或取大原则筛选出同源脉冲,其中,互相关系数判决门限设定为0.5~0.9之间的一个数值。
6.根据权利要求4所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第6步中,细关联是指对两个雷电脉冲的波形相似性进行最终判断,其过程为:
从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲的波形包络归一化均方偏差Ji,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预设均方偏差门限比较或取小原则,找出与A站该脉冲真正同源的B站脉冲,形成同源脉冲对;
其中,雷电脉冲的波形包络归一化均方偏差Ji为:
式中,N为雷电脉冲信号采样点的总个数,yA(n)和yBi(n)分别为A、B站在采样时刻n时的雷电脉冲包络波形,其采样点数分别为MA和MBi,GABi为B站脉冲包络相对于A站脉冲包络峰值的归一化系数,其值等于A站脉冲包络的最大值与B站脉冲包络的最大值之比,I为B站潜在同源脉冲簇队列中的脉冲个数;当MA≠MBi时,须对较短的脉冲波形补零直至与另一脉冲的长度相等,得到N=max{MA,MBi};
对计算出的全部均方偏差J1,J2,…,JI按大小排序,并按低于门限原则或取小原则筛选出同源脉冲,其中,均方偏差判决门限设定为0~1.0之间的一个数值。
7.根据权利要求4所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第6步中,细关联是指对两个雷电脉冲的波形相似性进行最终判断,其过程为:
从B站潜在同源脉冲簇中,按先后顺序逐次取出各个脉冲并计算其与A站该脉冲的波形包络偏度偏差Si或/和波形包络峰度偏差Ki,该计算遍历B站潜在同源脉冲簇队列中的所有脉冲,以此作为雷电脉冲是否真正同源的判决参数,通过与预设偏度偏差门限或/和峰度偏差门限比较,找出与A站该脉冲真正同源的B站脉冲,形成同源脉冲对,其中,偏度偏差判决门限和峰度偏差判决门限均为0~1.0之间的一个数值;
其中,雷电脉冲的波形包络偏度偏差Si和波形包络峰度偏差Ki分别为:
8.根据权利要求6、7或8所述的估计云闪雷电脉冲信号到达时间差的方法,其特征在于,第6步中,对同源脉冲对进行到达时间差估计过程为:
当同源脉冲对中的脉冲持续时间较短时,按下式分别计算同源脉冲对中的A站、B站脉冲的包络能量重心ΔA和ΔBi:
将脉冲包络的时间中心点定为参考零点,其前、后段设置为相同采样点数,则能量重心为相对于该参考零点的偏移量;
上述同源脉冲对中的两个脉冲的能量重心实际时刻tAG、tBGi分别为:
因此,该雷电脉冲到达A站、B站的时间差为:
式中,M′A=2NA为同源脉冲对中A站脉冲的长度,M′Bi=2NBi为同源脉冲对中B站脉冲的长度,T0为采样点间隔,tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻;
当同源脉冲对中的脉冲持续时间较长时,采用脉冲信号的互相关函数峰值点的位置作为雷电信号到达A站和B站的时间差,采用下述三种方式之一进行长脉冲情况下的到达时间差估计:
其中,E(·)为取平均值计算,xA(n)为A站还原出的雷电波形在时刻n时的采样值,xBi(n-τ)为B站还原出的雷电波形在时刻n-τ时的采样值;
式中,M为该段雷电信号的互相关函数的峰值点个数,对应于多源或多径雷电信号的个数;T0为采样点间隔;tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻;表示的绝对值取极大值时所对应的自变量τ之值;
此计算最大似然加权函数:
其中,
进而计算广义互相关函数:
由此得到该雷电脉冲到达A站和B站的时间差为:
式中,M为该段雷电信号的互相关函数的峰值点个数,对应于多源或多径雷电脉冲的个数;T0为采样点间隔;tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻;IDFT[.]表示逆傅里叶变换;第三种方式:在对前述计算出的互相关函数进行峰值搜索时,对得到的峰值点位置进行次极大值点线性插值和相邻三点抛物线插值后再计算新的峰值点位置,并以此作为到达时间差估计值输出;设对互相关函数进行搜索得到的峰值点为dm,其互相关函数值为dm点的前、后相邻点为dm-1和dm+1;dm点的前、后相邻次极大值点为dx-k1和dm+k2,其互相关函数值分别为 k1和k2为前、后相邻次极大值点到峰值点dm的距离;对dm-k1点、dm点、dm+k2点进行线性插值,并将插值结果分别赋予dm-1点和dm+1点:
式中,M为该段雷电信号中包含的多源或多径雷电信号的个数,T0为采样点间隔;tA和tBi分别表示同源脉冲对中的A站和B站脉冲的起始时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110199389 CN102298097B (zh) | 2011-07-15 | 2011-07-15 | 一种估计云闪雷电脉冲信号到达时间差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110199389 CN102298097B (zh) | 2011-07-15 | 2011-07-15 | 一种估计云闪雷电脉冲信号到达时间差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102298097A CN102298097A (zh) | 2011-12-28 |
CN102298097B true CN102298097B (zh) | 2013-02-13 |
Family
ID=45358652
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110199389 Expired - Fee Related CN102298097B (zh) | 2011-07-15 | 2011-07-15 | 一种估计云闪雷电脉冲信号到达时间差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102298097B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569913B (zh) * | 2015-01-30 | 2017-04-12 | 武汉大学 | 一种高精度全闪电定位方法 |
CN106383333A (zh) * | 2016-08-31 | 2017-02-08 | 北斗时空信息技术(北京)有限公司 | 一种基于互相关的改进时延估计方法 |
CN106879068B (zh) * | 2017-01-26 | 2020-05-08 | 浙江大学 | 一种强多径环境下信号的到达时间估计方法 |
CN108427099A (zh) * | 2017-02-15 | 2018-08-21 | 景斌 | 用于超声波室内定位的方法、设备、***以及移动终端 |
CN108520023B (zh) * | 2018-03-22 | 2021-07-20 | 合肥佳讯科技有限公司 | 一种基于混合聚类算法的雷暴核识别及追踪方法 |
CN109116301B (zh) * | 2018-08-14 | 2023-02-28 | 中国电子科技集团公司第三十八研究所 | 一种基于置信度估计的到达时间差测量方法 |
CN110346756B (zh) * | 2019-07-10 | 2023-06-06 | 四川中电昆辰科技有限公司 | 一种信号包络检测装置及方法及其到达时间修正方法 |
CN111487581B (zh) * | 2020-04-20 | 2022-02-11 | 中国人民解放军军事科学院国防科技创新研究院 | 基于舵矢量模式匹配的宽带闪电vhf辐射源测向方法 |
CN112881812B (zh) * | 2021-01-11 | 2022-02-22 | 中国气象科学研究院 | 基于机器学习编码的全闪实时定位方法及装置 |
CN114994467B (zh) * | 2022-05-26 | 2024-06-07 | 海南电网有限责任公司乐东供电局 | 一种基于长测试脉冲的电缆故障双端定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1210974A (zh) * | 1997-09-11 | 1999-03-17 | 武汉高压研究所新技术公司 | 一种脉冲电磁波到达时刻的探测方法 |
US5952969A (en) * | 1997-08-18 | 1999-09-14 | Telefonakiebolaget L M Ericsson (Publ) | Method and system for determining the position of mobile radio terminals |
CN101545932A (zh) * | 2009-04-28 | 2009-09-30 | 国网电力科学研究院 | 用于雷电探测的离线式大容量高速存储装置及其实现方法 |
-
2011
- 2011-07-15 CN CN 201110199389 patent/CN102298097B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5952969A (en) * | 1997-08-18 | 1999-09-14 | Telefonakiebolaget L M Ericsson (Publ) | Method and system for determining the position of mobile radio terminals |
CN1210974A (zh) * | 1997-09-11 | 1999-03-17 | 武汉高压研究所新技术公司 | 一种脉冲电磁波到达时刻的探测方法 |
CN101545932A (zh) * | 2009-04-28 | 2009-09-30 | 国网电力科学研究院 | 用于雷电探测的离线式大容量高速存储装置及其实现方法 |
Non-Patent Citations (4)
Title |
---|
115MHz频段雷电波形特征分析;杜海明等;《电网技术》;20101031;第34卷(第10期);60-64 * |
利用快电场变化脉冲定位进行云闪初始放电过程的研究;王东方等;《气象学报》;20091231;165-174 * |
杜海明等.115MHz频段雷电波形特征分析.《电网技术》.2010,第34卷(第10期), |
王东方等.利用快电场变化脉冲定位进行云闪初始放电过程的研究.《气象学报》.2009, |
Also Published As
Publication number | Publication date |
---|---|
CN102298097A (zh) | 2011-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102298097B (zh) | 一种估计云闪雷电脉冲信号到达时间差的方法 | |
CN105277923A (zh) | 一种单通道雷达信号分选方法 | |
CN103713277B (zh) | 一种基于位置信息场的辐射源定位方法 | |
CN105911521B (zh) | 联合无线电信号复包络和载波相位信息的超视距目标直接定位方法 | |
CN105929389A (zh) | 一种基于外辐射源时延和多普勒频率的直接定位方法 | |
CN101893698B (zh) | 噪声源测试分析方法及其装置 | |
CN102833016B (zh) | 一种在多径环境下针对同源信号时差的测定方法 | |
CN111157943B (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
CN105445699B (zh) | 一种非视距误差消除的测距方法及*** | |
CN110988834B (zh) | 一种基于自适应阈值双门限的脉冲到达时间测量方法 | |
CN113534144B (zh) | 基于时差约束残差的分布式信号分选方法 | |
CN109031190B (zh) | 一种高重频脉冲信号无源时差定位方法 | |
CN104219761A (zh) | 一种基于最大斜率的超宽带无线定位方法 | |
CN106291449A (zh) | 对称稳定分布噪声下波达方向角估计新方法 | |
CN105554882A (zh) | 基于能量检测的60GHz非视距识别与无线指纹定位方法 | |
CN106330342B (zh) | 一种低计算复杂度的水声通信多普勒因子估计方法 | |
CN111263321A (zh) | 一种提高toa室内超宽带定位精度的方法 | |
RU2439799C1 (ru) | Способ передачи и приема радиосигналов наземных радиомаяков | |
CN107728104A (zh) | 一种改进的卫星导航干扰测向方法 | |
CN109061556B (zh) | 一种基于弹性网络的稀疏迭代波达角估计方法 | |
CN114545342A (zh) | 利用多通道侦察接收机的雷达脉冲信号参数测量方法 | |
CN103686991A (zh) | 一种实现无线传感网节点定位的方法及装置 | |
CN102087313B (zh) | 一种卫星搜救信号的频率估计方法 | |
CN110208736B (zh) | 基于四阶累量的非圆信号均匀阵列波达方向角估计方法 | |
CN105738866A (zh) | 一种基于能量检测的60GHz非视距识别与无线指纹定位方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130213 |
|
CF01 | Termination of patent right due to non-payment of annual fee |