CN113037293A - 相位突变检测方法、装置、设备及存储介质 - Google Patents
相位突变检测方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN113037293A CN113037293A CN202110186387.4A CN202110186387A CN113037293A CN 113037293 A CN113037293 A CN 113037293A CN 202110186387 A CN202110186387 A CN 202110186387A CN 113037293 A CN113037293 A CN 113037293A
- Authority
- CN
- China
- Prior art keywords
- sampling
- data file
- processed
- phase jump
- sum
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 55
- 238000005070 sampling Methods 0.000 claims abstract description 283
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 38
- 230000035772 mutation Effects 0.000 claims abstract description 26
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 230000011218 segmentation Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Measuring Phase Differences (AREA)
Abstract
本发明公开了一种相位突变检测方法、装置、设备及存储介质。本发明通过读取待处理数据文件对应的采样点序列，根据预设频率的余弦信号对采样点序列进行混频处理，然后进行滤波处理，获得单频信号，根据单频信号中每个采样点的采样值确定相邻两个采样点的采样值平方和，在采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。本发明通过对待处理数据文件中的采样点序列进行混频以及滤波处理，得到单频信号，再获取单频信号中相邻两个采样点的采样值平方和，根据采样值平方和所在的区域范围将相位突变自动检测出来，相较于现有的通过人工逐段放大查看搜索相位突变点的方式，本发明上述方式能够更加准确地检测相位突变。
Description
技术领域
本发明涉及设备性能检验技术领域，尤其涉及一种相位突变检测方法、装置、设备及存储介质。
背景技术
目前，在采集设备性能检验时，通常会关注高速大容量数据存储落盘的正确性，即采集样点是否有丢失，测试信号通常为正弦信号，当丢失一个样点时也会发生相位突变，当分辨率足够精确时，会在频谱图上出现明显的光点。但是，对于大容量的落盘数据，通过人工逐段查看寻找相位不连续点过于耗时耗精力，劳动强度大，自动化程度低，并且容易出现漏警。在现有技术中，对于设备高速大容量采集存储数据正确性的测试，并不具备通用的性能稳定的测试方法。
上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。
发明内容
本发明的主要目的在于提供了一种相位突变检测方法、装置、设备及存储介质，旨在解决现有技术对高速大容量采集存储数据的相位突变检测不准确的技术问题。
为实现上述目的，本发明提供了一种相位突变方法，所述相位突变检测方法包括：
读取待处理数据文件对应的采样点序列；
根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
对所述双频信号进行滤波处理，以获得单频信号；
获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
优选地，所述读取待处理数据文件对应的采样点序列的步骤，具体包括：
对待处理数据文件进行分段处理，获得分段后的数据文件；
从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列。
优选地，所述在所述平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变的步骤之后，还包括：
确定所述相位突变对应的采样值平方和，并返回所述从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列的步骤，直至得到所述待处理数据文件中所有发生相位突变的采样值平方和集合；
获取所述采样值平方和集合中的局部极值；
根据所述局部极值确定所述采样值平方和对应的采样点在所述待处理数据文件中的索引。
优选地，所述获取所述采样值平方和集合中的局部极值的步骤，具体包括：
获取所述采样值平方和集合中的所有目标采样值平方和；
判断所述目标采样值平方和是否小于预设位置的采样值平方和；
若是，则所述目标采样值平方和为所述采样值平方和集合中的局部极值。
优选地，所述对所述待处理数据文件进行分段处理，获得分段后的数据文件的步骤，具体包括：
根据所述待处理数据文件确定采样点总数，根据预设处理能力确定每段检测采样点数；
根据所述采样点总数和所述每段检测采样点数确定所述待处理数据文件的段数；
根据所述段数对所述待处理数据文件进行分段处理，获得分段后的数据文件。
优选地，所述根据所述局部极值确定所述值平方和对应的采样点在所述待处理数据文件中的索引的步骤之后，还包括：
根据所述索引从所述待处理数据文件中截取预设长度的数据文件；
基于所述预设长度的数据文件验证所述索引对应的采样点是否出现相位突变。
优选地，所述在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变的步骤之前，还包括：
根据所述采样值平方和确定均值和均方差；
根据所述均值和均方差确定目标阈值范围。
此外，为实现上述目的，本发明还提出一种相位突变检测装置，所述相位突变检测装置包括：
数据读取模块，用于读取待处理数据文件对应的采样点序列；
频率搬移模块，用于根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
带通滤波模块，用于对所述获得双频信号进行滤波处理，以获得单频信号；
数据计算模块，用于获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
突变判断模块，用于在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
此外，为实现上述目的，本发明还提出一种相位突变检测设备，所述相位突变检测设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的相位突变检测程序，所述相位突变检测程序配置为实现如上文所述的相位突变检测方法的步骤。
此外，为实现上述目的，本发明还提出一种存储介质，所述存储介质上存储有相位突变检测程序，所述相位突变检测程序被处理器执行时实现如上文所述的相位突变检测方法的步骤。
本发明通过读取待处理数据文件对应的采样点序列，根据预设频率的余弦信号对采样点序列进行混频处理，获得双频信号，对双频信号进行滤波处理，获得单频信号，获取单频信号中每个采样点的采样值，根据每个采样点的采样值确定相邻两个采样点的采样值平方和，在采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。本发明通过对待处理数据文件中的采样点序列进行混频以及滤波处理，得到单频信号，再获取单频信号中相邻两个采样点的采样值平方和，根据采样值平方和所在的区域范围将相位突变自动检测出来，相较于现有的通过人工逐段放大查看搜索相位突变点的方式，本发明上述方式能够更加准确地检测相位突变。
附图说明
图1是本发明实施例方案涉及的硬件运行环境的相位突变检测设备的结构示意图；
图2为本发明相位突变检测方法第一实施例的流程示意图；
图3为本发明相位突变检测方法第一实施例采样点序列的频谱示意图；
图4为本发明相位突变检测方法第一实施例双频信号的频谱示意图；
图5为本发明相位突变检测方法第一实施例带通滤波器幅频和相频特性示意图；
图6为本发明相位突变检测方法第一实施例单频信号的频谱示意图；
图7为本发明相位突变检测方法第二实施例的流程示意图；
图8为本发明相位突变检测方法第三实施例的流程示意图；
图9为本发明相位突变检测方法第三实施例相邻两个采样点的采样值平方和示意图；
图10为本发明相位突变检测方法第三实施例回溯验证示意图；
图11为本发明相位突变检测装置第一实施例的结构框图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
参照图1，图1为本发明实施例方案涉及的硬件运行环境的相位突变检测设备结构示意图。
如图1所示，该相位突变检测设备可以包括：处理器1001，例如中央处理器(Central Processing Unit，CPU)，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(WIreless-FIdelity，WI-FI)接口)。存储器1005可以是高速的随机存取存储器(RandomAccess Memory，RAM)存储器，也可以是稳定的非易失性存储器(Non-Volatile Memory，NVM)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的结构并不构成对相位突变检测设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种存储介质的存储器1005中可以包括操作系统、数据存储模块、网络通信模块、用户接口模块以及相位突变检测程序。
在图1所示的相位突变检测设备中，网络接口1004主要用于与网络服务器进行数据通信；用户接口1003主要用于与用户进行数据交互；本发明相位突变检测设备中的处理器1001、存储器1005可以设置在相位突变检测设备中，所述相位突变检测设备通过处理器1001调用存储器1005中存储的相位突变检测程序，并执行本发明实施例提供的相位突变检测方法。
本发明实施例提供了一种相位突变检测方法，参照图2，图2为本发明相位突变检测方法第一实施例的流程示意图。
本实施例中，所述相位突变检测方法包括以下步骤：
步骤S10：读取待处理数据文件对应的采样点序列；
需要说明的是，本实施例方法的执行主体可以是高速大容量数据采集设备，以下简称采集设备，上述采集设备可以应用于军事、航空航天、通信、探测等需要采集高速大容量数据的领域，本实施例对此不做具体限制。
可理解的是，所述待处理数据文件可以是由上述采集设备采集的数据文件，所述采样点序列指的是待处理数据文件中每个采样点对应的序列。
在具体实现中，参照图3，图3为本发明相位突变检测方法第一实施例采样点序列的频谱示意图。如图3所示，本实施例中，采集设备可以通过COOLEDIT软件查看采样点序列对应的频谱图，便于用户更方便地观察上述采样点序列，例如：在图3中，采集到的采样点序列对应的正弦信号的频率为70MHz，采样率为200MHz。
步骤S20：根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
需要说明的是，所述混频处理是指将采样点序列与预设频率的余弦信号进行相乘，获得新生成的双频信号。
可理解的是，所述预设频率的余弦信号可以对采样点序列进行混频处理，获得双频信号，以使双频信号的频率变为采样率的四分之一，并且获得双频信号的和频和差频，和频是指两个信号的频率之和，差频是指两个信号的频率之差，例如：两个信号的频率分别为30Hz和50Hz，和频为80Hz，差频为20Hz。
在具体实现中，预设频率是根据采样率来设置的，例如：上述采样点序列对应的采样率为200MHz，采样率的四分之一即50MHz，要获得50MHz的频率，就要更改采样点序列对应的频率，即进行混频处理，如果是和频为50MHz，余弦信号的频率就需要设置为-20MHz，不符合实际场景，所以需要设置差频为50MHz，在本实施例中，设置预设频率为20MHz，即余弦信号的频率为20MHz，如图4所示，图4为本发明相位突变检测方法第一实施例双频信号的频谱示意图，其中，双频信号的和频为90MHz(70MHz+20MHz)，差频为50MHz(70MHz-20MHz)。
可理解的是，上述内容只列举了将差频设置为采样率的四分之一，在实际操作中，也可以根据应用场景将和频设置为采样率的四分之一，本实施例对此不做具体限制。
步骤S30：对所述双频信号进行滤波处理，以获得单频信号；
应理解的是，因为所述双频信号有两个频率，不便于具体分析，因此需要从双频信号中滤出只有一个频率的单频信号，参考图6，图6为本发明相位突变检测方法第一实施例单频信号的频谱示意图。
需要说明的是，可以通过带通滤波器进行滤波处理，也可以通过高通滤波器或者低通滤波器，本实施例对此不做具体限制。
在具体实现中，采集设备需要滤出和频或差频为采样率的四分之一的信号，例如：上述双频信号的差频为50MHz，由此可以设置带通滤波器的通带频率范围为[0.4π，0.6π]，只需要使50MHz的频率在带通滤波器的通带频率范围里，就可以实现滤波，也就是说，还可以设置带通滤波器的通带频率范围为[0.3π，0.7π]，带通滤波器的阶数可以设置为1000，也可以设置为1500，本实施例对此不做具体限制，参考图5，图5为本发明相位突变检测方法第一实施例带通滤波器幅频和相频特性示意图。
步骤S40：获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
需要说明的是，可以通过COOLEDIT软件获取所述单频信号中每个采样点的采样值，也可以通过其他具有类似功能的软件获取采样值，本实施例对此不做具体限制。
可理解的是，相邻两个采样点的采样值平方和是根据每个采样点的采样值确定的，例如：采样到了5个采样点，就需要获取4个采样值平方和，假设这5个采样点的采样值分别为0,0.1,0.3,0.5,0.6，本实施例需要获取所有相邻的两个采样点的采样值平方和，即0.01(02+0.12)，0.10(0.12+0.32)，0.34(0.32+0.52),0.61(0.52+0.62)。
步骤S50：在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
应理解的是，所述采样值平方和的数值大部分都在一定范围内，因为通过COOLEDIT软件可以看到所述采样点序列可以拟合为一个正弦信号，如果采样值平方和没有在大部分采样值平方和所在的范围内，即可判定存在相位突变。
本实施例采集设备通过读取待处理数据对应的采样点序列，根据预设频率的余弦信号对采样点序列进行混频处理，获得双频信号，对双频信号进行滤波处理，获得单频信号，获取单频信号中每个采样点的采样值，根据每个采样点的采样值确定相邻两个采样点的采样值平方和，在采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。本发明通过对待处理数据中的采样点序列进行混频以及滤波处理，得到单频信号，再获取单频信号中相邻两个采样点的采样值平方和，根据采样值平方和所在的区域范围将相位突变自动检测出来，相较于现有的通过人工逐段放大查看搜索相位突变点的方式，本发明上述方式能够更加准确地检测相位突变。
进一步地，为了确定所述目标阈值范围，本实施例中，上述步骤S50之前，还包括：
步骤S1：根据所述采样值平方和确定均值和均方差；
需要说明的是，由于待处理数据对应的采样点序列中出现相位突变的点很少，可以只选取起始的一段序列，序列长度可以设置为1M，也可以设置为2M，本实施例对此不做具体限制。
在具体实现中，所述均值为每相邻两个采样点的采样值平方和的平均值，所述均方差为每相邻两个采样点的采样值平方和与均值之差的平方值的平均数的算术平方根，例如：5个采样点的采样值分别为0,0.1,0.3,0.5,0.6，每相邻两个采样点的采样值平方和分别为0.01,0.10,0.34,0.61，均值为(0.01+0.10+0.34+0.61)/4＝0.265，方差为((0.01-0.265)2+(0.10-0.265)2+(0.34-0.265)2+(0.61-0.265)2/4＝0.054，均方差为
在实际处理过程中，所要处理的序列不可能只有5个采样值，对于多个采样值也是进行相同的操作，上述举例仅供理解。
S2：根据所述均值和均方差确定目标阈值范围。
需要说明的是，所述目标阈值范围的上限为均值加上均方差的8倍，所述目标阈值范围的下限为均值减去均方差的8倍。例如：延续上面所列举的例子，均值为0.265，均方差为0.232，目标阈值范围的上限为0.265+8*0.232＝2.121，目标阈值范围的下限为0.265-8*0.232＝-1.591，即目标阈值范围为[-1.591，2.121]。
本实施例通过确定每相邻两个采样点的采样值平方和的均值和均方差，然后根据均值和均方差确定目标阈值范围，而不是仅仅通过均值来确定目标阈值范围，本实施例能够提高相位突变检测的准确性，保证了相位突变检测的可靠性。
参考图7，图7为本发明相位突变检测方法第二实施例的流程示意图。
基于上述第一实施例，在本实施例中，所述步骤S10包括：
步骤S101：对待处理数据文件进行分段处理，获得分段后的数据文件；
需要说明的是，所述分段处理是指把待处理数据文件分段成为若干个小数据文件，分段方式可以是根据待处理数据文件的大小确定的，例如，待处理数据文件大小为32M，就可以分段成为16个小数据文件，每个小数据文件大小为2M；分段方式还可以是根据待处理数据文件中的标识确定的，例如，待处理数据文件里有10个标识，就可以根据标识把待处理数据文件分成10个小数据文件，本实施例对此不做具体限制。
可理解的是，采集设备通常采集到的数据文件里面的数据是极大容量的数据，本实施例通过对待处理数据文件进行分段处理，获得若干个小文件，以便更加快速地处理。
可理解的是，为了便于对待处理数据文件进行分段处理，本实施例中，上述步骤S101，可具体包括：
步骤S1011：根据所述待处理数据文件确定采样点总数，根据预设处理能力确定每段检测采样点数；
需要说明的是，采样点总数是根据待处理文件的大小来决定的，在具体实现中，在把待处理文件写入COOLEDIT软件时，会设置待处理文件的存储格式，根据存储格式可以确定采样点总数，所述预设处理能力指的是上述采集设备的存储能力和计算能力，通过采集设备的存储能力和计算能力可以确定每段检测采样点数。
步骤S1012：根据所述采样点总数和所述每段检测采样点数确定所述待处理数据文件的段数；
在具体实现中，段数可由采样点总数除以每段检测采样点数获得，例如：采样点总数为32M，每段检测采样点数为1M，段数则为32段，如果遇到不能除尽的情况，则向上取整，例如：采样点总数为32M，每段检测采样点数为1.5M，段数则为22段。
步骤S1013：根据所述段数对所述待处理数据文件进行分段处理，获得分段后的数据文件。
可理解的是，根据段数就可以对待处理数据文件进行分段处理，例如：段数为32段，就会把待处理文件分成32段小文件。
步骤S102：从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列。
本实施例通过采样点总数和采集设备每段检测采样点数确定段数，对待处理数据文件进行分段，并分段读取待处理数据文件中的采样点序列，把大容量数据文件分段成若干个小数据文件，使得采集设备能够根据自身的处理能力进行分段，并且使采集设备工作达到饱和状态，能够更加快速地检测相位突变点，提高了相位突变检测的准确性。
参考图8，图8为本发明相位突变检测方法第三实施例的流程示意图。
基于上述各实施例，在本实施例中，所述步骤S50之后，所述方法还包括：
步骤S60：确定所述相位突变对应的采样值平方和，并返回所述从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列的步骤，直至得到所述待处理数据文件中所有发生相位突变的采样值平方和集合；
需要说明的是，所述采样值平方和集合里面包含待处理数据文件中所有发生相位突变的采样值平方和。
可理解的是，在检测到出现相位突变后，需要记录相位突变对应的采样值平方和，以实现精准定位相位突变点。
在具体实现中，检测完第一段数据文件是否发生相位突变后，需要对后面没有检测的数据文件进行处理，以获得完整的待处理文件的检测结果，所以需返回步骤S102。
可理解的是，在读取选取的数据文件对应的采样点序列时，需要判断该段数据文件是否为最后一段数据文件，如果是，则只需要读取剩余的采样点序列，因为可能会出现最后一段数据文件与前面的数据文件长度不一致的情况，如果不是最后一段数据文件，则依旧读取预先设置的采样点数，例如：除了最后一段数据文件采样点数为1M，其余数据文件采样点数都为1.5M，在读取的过程中，判断是否为最后一段数据文件，若是，则读取1M的采样点数，若不是，则读取1.5M的采样点数，然后继续读取下一段数据文件。
应理解的是，在实际应用场景中，由于与带通滤波器进行线性卷积会导致相位不连续的情况，即相位连续性被破坏的点数是滤波器点数的两倍，所以分段时需要采取段间重叠的策略，段间重叠的采样点数为滤波器点数的两倍，这样就可以检测到所有采样点，例如：设置段间重叠点数为2000，则第一段读取1～10001之间的采样点，第二段读取8001～18001之间的采样点，第二段读取16001～26001之间的采样点。
步骤S70：获取所述采样值平方和集合中的局部极值；
需要说明的是，所述局部极值可以是局部极低值，也可以是局部极高值，根据具体情况具体分析，本实施例对此不做具体限制，在本实施例中，优选为局部极低值。
如图9所示，图9为本发明相位突变检测方法第三实施例相邻两个采样点的采样值平方和示意图，可以看出相邻两个采样点的采样值平方和存在局部极值。
可理解的是，所述局部极值可以是通过比较当前区域的所有采样点的采样值来获得，也可以是通过比较局部区域的采样点的采样值来获得，本实施例对此不做限制。
进一步地，为了获取采样值的局部极值，本实施例中，上述步骤S70，可具体包括：
步骤S701：获取所述采样值平方和集合中的所有目标采样值平方和；
步骤S702：判断所述目标采样值平方和是否小于预设位置的采样值平方和；
需要说明的是，所述预设位置为目标采样点的前面三个以及后面三个采样点，当目标采样点的采样值同时小于前面三个采样点的采样值以及后面单个采样点的采样值时，才会判定当前目标采样点的采样值为局部极值。
步骤S703：若是，则所述目标采样值平方和为所述采样值平方和集合中的局部极值。
本实施例通过获取所述采样值平方和集合中的所有目标采样值平方和，然后判断所述目标采样值平方和是否小于预设位置的采样值平方和，若是，则所述目标采样值平方和为所述采样值平方和集合中的局部极值，这种方式能够更加快速、准确地确定局部极值，并且提高了相位突变检测的准确性。
步骤S80：根据所述局部极值确定所述采样值平方和对应的采样点在所述待处理数据文件中的索引。
需要说明的是，所述采样值平方和对应的采样点由于是通过相邻两个采样点的采样值平方和确定的，所以采样点可能会出现两个，对应的索引任取一个即可。
在具体实现中，用户可以获取局部最低值对应的采样点在当前段的数据文件中的索引，再加上前面所有段的数据文件中的采样点就可以得到所述局部极值确定所述采样值平方和对应的采样点在所述待处理数据文件中的索引，例如：局部极值对应的采样点在当前段的数据文件中的索引为200，在此之前有两段数据文件，每段数据文件有1000个采样点，局部极值对应的采样点在待处理数据文件中的索引即为1000*2+200＝2200。
为了实现回溯验证，本实施例中，上述步骤S80之后，还包括：
步骤S90：根据所述索引从所述待处理文件中截取预设长度的数据文件；
需要说明的是，所述预设长度的数据文件为索引前后的1000个采样点对应的数据文件，例如，索引为2000，就需要截取1500～2500之间的采样点。
步骤S100：基于所述预设长度的数据文件验证所述索引对应的采样点是否出现相位突变。
在具体实现中，可以将上述预设长度的数据文件写入COOLEDIT软件中，并进行查看，如图10所示，图10为本发明相位突变检测方法第三实施例回溯验证示意图，图中有光点的位置即产生相位突变的位置，可以验证之前检测到的相位突变点的位置是否正确。
本实施例通过确定相位突变对应的采样点，并返回从所述分段后的数据文件中选取一段采样点序列，并读取选取的数据文件对应的采样点序列的步骤，直至得到待处理数据文件中所有发生相位突变的采样点集合，并确定采样点集合里面的所有采样点对应的索引，能够实现精准定位相位突变点所在的位置，并确定相位突变点在待处理文件中所有采样点的索引，方便用户获取到待处理数据文件中相位突变点的具体位置，还能进一步验证相位突变点检测是否正确。
此外，本发明实施例还提出一种存储介质，所述存储介质上存储有相位突变检测程序，所述相位突变检测程序被处理器执行时实现如上文所述的相位突变检测方法的步骤。
参照图11，图11为本发明相位突变检测装置第一实施例的结构框图。
如图11所示，本发明实施例提出的相位突变检测装置包括：
数据读取模块1101，用于读取待处理数据文件对应的采样点序列；
频率搬移模块1102，用于根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
带通滤波模块1103，用于对所述获得双频信号进行滤波处理，以获得单频信号；
数据计算模块1104，用于获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
突变判断模块1105，用于在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
本实施例通过读取待处理数据文件对应的采样点序列，根据预设频率的余弦信号对采样点序列进行混频处理，获得双频信号，对双频信号进行滤波处理，获得单频信号，获取单频信号中每个采样点的采样值，根据每个采样点的采样值确定相邻两个采样点的采样值平方和，在采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。本发明通过对待处理数据文件中的采样点序列进行混频以及滤波处理，得到单频信号，再获取单频信号中相邻两个采样点的采样值平方和，根据采样值平方和所在的区域范围将相位突变自动检测出来，相较于现有的通过人工逐段放大查看搜索相位突变点的方式，本发明上述方式能够更加准确地检测相位突变。
本发明相位突变检测装置的其他实施例或具体实现方式可参照上述各方法实施例，此处不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如只读存储器/随机存取存储器、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种相位突变检测方法，其特征在于，所述相位突变检测方法包括：
读取待处理数据文件对应的采样点序列；
根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
对所述双频信号进行滤波处理，以获得单频信号；
获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
2.如权利要求1所述的相位突变检测方法，其特征在于，所述读取待处理数据文件对应的采样点序列的步骤，具体包括：
对待处理数据文件进行分段处理，获得分段后的数据文件；
从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列。
3.如权利要求2所述的相位突变检测方法，其特征在于，所述在所述平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变的步骤之后，还包括：
确定所述相位突变对应的采样值平方和，并返回所述从所述分段后的数据文件中选取一段数据文件，并读取选取的数据文件对应的采样点序列的步骤，直至得到所述待处理数据文件中所有发生相位突变的采样值平方和集合；
获取所述采样值平方和集合中的局部极值；
根据所述局部极值确定所述采样值平方和对应的采样点在所述待处理数据文件中的索引。
4.如权利要求3所述的相位突变检测方法，其特征在于，所述获取所述采样值平方和集合中的局部极值的步骤，具体包括：
获取所述采样值平方和集合中的所有目标采样值平方和；
判断所述目标采样值平方和是否小于预设位置的采样值平方和；
若是，则所述目标采样值平方和为所述采样值平方和集合中的局部极值。
5.如权利要求2所述的相位突变检测方法，其特征在于，所述对所述待处理数据文件进行分段处理，获得分段后的数据文件的步骤，具体包括：
根据所述待处理数据文件确定采样点总数，根据预设处理能力确定每段检测采样点数；
根据所述采样点总数和所述每段检测采样点数确定所述待处理数据文件的段数；
根据所述段数对所述待处理数据文件进行分段处理，获得分段后的数据文件。
6.如权利要求3所述的相位突变检测方法，其特征在于，所述根据所述局部极值确定所述采样值平方和对应的采样点在所述待处理数据文件中的索引的步骤之后，还包括：
根据所述索引从所述待处理数据文件中截取预设长度的数据文件；
基于所述预设长度的数据文件验证所述索引对应的采样点是否出现相位突变。
7.如权利要求1至6任一项所述的相位突变检测方法，其特征在于，所述在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变的步骤之前，还包括：
根据所述采样值平方和确定均值和均方差；
根据所述均值和均方差确定目标阈值范围。
8.一种相位突变检测装置，其特征在于，所述相位突变检测装置包括：
数据读取模块，用于读取待处理数据文件对应的采样点序列；
频率搬移模块，用于根据预设频率的余弦信号对所述采样点序列进行混频处理，以获得双频信号；
带通滤波模块，用于对所述获得双频信号进行滤波处理，以获得单频信号；
数据计算模块，用于获取所述单频信号中每个采样点的采样值，根据所述采样值确定相邻两个采样点的采样值平方和；
突变判断模块，用于在所述采样值平方和不属于目标阈值范围时，判定所述采样点序列存在相位突变。
9.一种相位突变检测设备，其特征在于，所述相位突变检测设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的相位突变检测程序，所述相位突变检测程序配置为实现如权利要求1至7中任一项所述的相位突变检测方法的步骤。
10.一种存储介质，其特征在于，所述存储介质上存储有相位突变检测程序，所述相位突变检测程序被处理器执行时实现如权利要求1至7任一项所述的相位突变检测方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110186387.4A CN113037293A (zh) | 2021-02-07 | 2021-02-07 | 相位突变检测方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110186387.4A CN113037293A (zh) | 2021-02-07 | 2021-02-07 | 相位突变检测方法、装置、设备及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113037293A true CN113037293A (zh) | 2021-06-25 |
Family
ID=76461238
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110186387.4A Pending CN113037293A (zh) | 2021-02-07 | 2021-02-07 | 相位突变检测方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113037293A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114109798A (zh) * | 2021-11-26 | 2022-03-01 | 广东美的暖通设备有限公司 | 频率确定方法、频率确定装置、压缩机系统和存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103856366A (zh) * | 2012-12-06 | 2014-06-11 | 腾讯科技（深圳）有限公司 | 一种平台数据监控方法及系统 |
| US20160019587A1 (en) * | 2012-12-30 | 2016-01-21 | Certona Corporation | Extracting predictive segments from sampled data |
| CN106405285A (zh) * | 2016-08-30 | 2017-02-15 | 华北电力大学 | 一种电力系统故障录波数据突变时刻检测方法及系统 |
| CN106708694A (zh) * | 2015-11-18 | 2017-05-24 | 腾讯科技（深圳）有限公司 | 一种定位毛刺异常点的方法、装置及计算设备 |
| CN109981517A (zh) * | 2019-01-22 | 2019-07-05 | 西安电子科技大学 | 一种基于fpga的qpsk神经网络解调器及其控制方法 |
| CN110579682A (zh) * | 2019-08-22 | 2019-12-17 | 国网安徽省电力有限公司 | 一种故障录波数据暂态同源对比方法及装置 |
-
2021
- 2021-02-07 CN CN202110186387.4A patent/CN113037293A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103856366A (zh) * | 2012-12-06 | 2014-06-11 | 腾讯科技（深圳）有限公司 | 一种平台数据监控方法及系统 |
| US20160019587A1 (en) * | 2012-12-30 | 2016-01-21 | Certona Corporation | Extracting predictive segments from sampled data |
| CN106708694A (zh) * | 2015-11-18 | 2017-05-24 | 腾讯科技（深圳）有限公司 | 一种定位毛刺异常点的方法、装置及计算设备 |
| CN106405285A (zh) * | 2016-08-30 | 2017-02-15 | 华北电力大学 | 一种电力系统故障录波数据突变时刻检测方法及系统 |
| CN109981517A (zh) * | 2019-01-22 | 2019-07-05 | 西安电子科技大学 | 一种基于fpga的qpsk神经网络解调器及其控制方法 |
| CN110579682A (zh) * | 2019-08-22 | 2019-12-17 | 国网安徽省电力有限公司 | 一种故障录波数据暂态同源对比方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 安涛: "一种基于FPGA的检测单频信号相位突变的实现方法", 《第三十二届中国（天津）2018’IT、网络、信息技术、电子、仪器仪表创新学术会议论文集》 * |
| 张初础: "《银行业计算机辅助审计的技术与方法》", 31 July 2009 * |
| 曹锋: "《虚拟仪器及应用》", 31 March 2018 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114109798A (zh) * | 2021-11-26 | 2022-03-01 | 广东美的暖通设备有限公司 | 频率确定方法、频率确定装置、压缩机系统和存储介质 |
| CN114109798B (zh) * | 2021-11-26 | 2023-08-29 | 广东美的暖通设备有限公司 | 频率确定方法、频率确定装置、压缩机系统和存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109522895B (zh) | 兴趣点位置校验方法、装置、服务器和计算机可读介质 | |
| CN109996269A (zh) | 一种无线通信网络异常原因确定方法、装置、设备及介质 | |
| US7548201B2 (en) | Method and system for automatically analyzing and modifying cable television signal leak information | |
| CN113037293A (zh) | 相位突变检测方法、装置、设备及存储介质 | |
| US6801162B1 (en) | Doppler-based automated direction finding system and method for locating cable television signal leaks | |
| CN111367958A (zh) | 一种数据分析方法、装置、电子设备及存储介质 | |
| CN113140109B (zh) | 一种路测数据处理方法、装置、计算机设备及存储介质 | |
| CN109271453B (zh) | 一种数据库容量的确定方法和装置 | |
| CN114332012A (zh) | 缺陷检测方法、装置、设备及计算机可读存储介质 | |
| CN110968657A (zh) | 感兴趣区域查询方法、装置、设备及计算机可读介质 | |
| CN117252837A (zh) | 一种晶圆测试的数据处理方法、装置、介质和电子设备 | |
| CN116660176A (zh) | 傅里叶光谱自动基线校正方法、装置和存储介质 | |
| CN109308235B (zh) | 点检方法及相关装置 | |
| CN115902706A (zh) | 避雷器的缺陷分析方法、装置、设备及存储介质 | |
| CN114125864B (zh) | 基于mdt的邻区检测方法及装置 | |
| CN113055239B (zh) | 数据显示方法、装置、设备及介质 | |
| CN114113933A (zh) | 局部放电信号滤波方法、装置、计算机设备和存储介质 | |
| CN114109794B (zh) | 转向确定方法、装置、可读存储介质和空气调节设备 | |
| CN111710347A (zh) | 音频数据分析方法、电子设备及存储介质 | |
| CN111159172A (zh) | 一种数据处理方法、装置及电子设备 | |
| JP5391635B2 (ja) | 解析装置、データ保存方法およびデータ保存プログラム | |
| CN117371385B (zh) | 邻近图形的确认方法、电子设备及存储介质 | |
| KR102544621B1 (ko) | 시계열 데이터의 에러 값 보정을 위한 필터링 및 성능 비교 시스템 및 방법 | |
| WO2022143713A1 (zh) | 一种v-slam地图校验方法、装置及设备 | |
| CN108430045B (zh) | 一种字符串修剪方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210625 |