CN115598416A - 台区采样信号的处理方法、系统、存储介质及计算机设备 - Google Patents
台区采样信号的处理方法、系统、存储介质及计算机设备 Download PDFInfo
- Publication number
- CN115598416A CN115598416A CN202211130264.XA CN202211130264A CN115598416A CN 115598416 A CN115598416 A CN 115598416A CN 202211130264 A CN202211130264 A CN 202211130264A CN 115598416 A CN115598416 A CN 115598416A
- Authority
- CN
- China
- Prior art keywords
- sampling
- sequence
- sampling sequence
- signal
- peak value
- 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.)
- Granted
Links
- 238000005070 sampling Methods 0.000 title claims abstract description 224
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000012545 processing Methods 0.000 title claims abstract description 27
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 230000006870 function Effects 0.000 claims description 33
- 238000004590 computer program Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 238000003672 processing method Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 230000002999 depolarising effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 6
- 230000001360 synchronised effect Effects 0.000 abstract description 6
- 238000009826 distribution Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
Abstract
本发明实施例提供了一种台区采样信号的处理方法、系统、存储介质及计算机设备,所述方法包括:获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;对所述第一采样序列进行插值迭代,得到第二采样序列;对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。本发明实施例通过对采样序列迭代运算,解决了采样频率与信号频率不能完全同步的问题,减少了频谱泄露现象,同时,提高了基于特征电流分析台区拓扑识别的成功率。
Description
技术领域
本实施例涉及电力电子技术领域,尤其涉及台区采样信号的处理方法、系统、存储介质和计算机设备。
背景技术
随着电力电子技术的发展,各种非线性负载在工农业生产和输用电设备中得到了广泛的应用,使得电网波形畸变,即谐波干扰越来越严重。导致基于特征电流分析台区拓扑识别越来越具有挑战性,在干扰越严重的情况下其拓扑识别成功率越低。
现有技术方案中,拓扑识别的方法是先对台区电信号AD采样,经过AD采样后,采用傅里叶变换,将时域信号转换为频域信号,提取目标频点的分量,实现其拓扑识别运算。
然而这种方式由于硬件本身的特性原因,特别是AD采样器件采集信号进行模数转换时及滤波器件滤波时数据都需要一个稳定时间,同时在模拟信号转变成数字信号的过程中,会有大量的噪声(参见图1),而且采样和量化误差使得其数据的连贯性并不能做到完全的一致,即采样频率与信号频率并不能完全同步,造成周期采样信号的相位在始端和终端不连续,进而在后级数据运算过程中,即频域变换中,极其容易出现频谱泄露现象,导致信号离散傅里叶变换的结果产生测试误差,影响谐波分析的精度。
发明内容
本发明实施例提供一种台区采样信号的处理方法、系统、存储介质和计算机设备,用于解决采样频率与信号频率不同步导致的频谱泄露问题,提高了谐波分析的精度和准确度。
本发明的技术方案为:
第一方面,本申请实施例提供了一种台区采样信号的处理方法,包括:
获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;
对所述第一采样序列进行插值迭代,得到第二采样序列;
对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
进一步地,所述对所述第一采样序列进行插值迭代,得到第二采样序列,包括:
步骤1,对所述第一采样序列中的相邻数值之间进行正弦函数或余弦函数插值,得到插值后的采样序列;
步骤2,在所述插值后的采样序列中删除原采样序列,得到迭代后的采样序列;
步骤3,判断步骤1的执行次数,当所述步骤1的执行次数少于m,则执行步骤4;当所述步骤1的执行次数为m,则执行步骤5,其中,2≤m≤20;
步骤4,将所述迭代后的采样序列作为所述第一采样序列,从步骤1开始重新执行;
步骤5,将所述迭代后的采样序列作为第二采样序列。
进一步地,对所述第一采样序列进行插值迭代之前,还包括对所述第一采样序列进行去极值滤波,所述去极值滤波的方法包括:
对所述第一采样序列中符合一定预设条件的采样点进行重新赋值;
所述符合一定预设条件的采样点包括:所述第一采样序列中的各采样点的采样值高于预设峰值的采样点。
进一步地,所述对所述第一采样序列中符合一定预设条件的采样点进行重新赋值的方法包括:
判断所述第一采样序列中是否出现采样值大于第一峰值的采样点;
若所述第一采样序列中出现采样值大于第一峰值的采样点,对采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值,将赋值后的采样序列作为第一采样序列。
进一步地,所述第一峰值的获取方法包括:
其中,A表示第一峰值,Xi为所述第一采样序列中第i个采样点的采样数值,n为所述第一采样序列中采样点的数量。
进一步地,所述若所述第一采样序列中出现采样值大于第一峰值的采样点,对采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值之后还包括:判断所述第一采样序列中每1~1.5个工频周期中波峰和波谷的数量;当有且只有一个波峰和波谷,数据不作处理;当出现两个或者两个以上的波峰和波谷,则对出现的所有波峰和波谷对应的采样点,采用正弦函数或余弦函数插值法重新赋值,将赋值后的采样序列作为第一采样序列。
进一步地,所述方法还包括:依据所述频谱信息,确定所述第一采样信号是否包括台区设备发送的特征信号。
进一步地,所述依据所述频谱信息,确定所述第一采样信号是否包括台区设备发送的特征信号,具体包括:
将f1+f2作为第一特征频率,f1-f2作为第二特征频率,其中f1为所述特征信号对应的频率,f2为工频;
根据所述频谱信息中频率与模值的对应关系,确定所述第一特征频率和/或所述第一特征频率所对应的最低信号量;
在所述最低信号量高于第一预设信号量时,确定所述第一采样信号包括台区设备发送的特征信号。
第二方面,本申请实施例还提供了一种台区采样信号的处理系统,包括:
采样信号获取单元,用于获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;
采样信号处理单元,用于对所述第一采样序列进行插值迭代,得到第二采样序列;
频谱信息提取单元,用于对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
第三方面,本申请实施例还提供了一种计算机存储介质,所述计算机存储介质存储有多条指令,所述指令适于有处理器加载并执行上述方法步骤。
第四方面,本申请实施例还提供了一种计算机设备,可包括处理器和存储器,所述存储器存储有计算机程序,所述计算机程序适于由所述处理器加载并执行上述方法步骤。
本申请一些实施例提供的技术方案带来的有益效果至少包括:通过对采样序列进行插值迭代,并对迭代后的信号进行傅里叶变换,从而得到对应的频谱信息,解决了采样频率和信号频率不能完全同步的问题,避免了频谱泄露的现象,消除了测试误差,提高了谐波的分析精度。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本实施例的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为现有技术方案中台区电信号目标频点的分量提取的方法流程图;
图2为本发明实施例的方法流程图;
图3为本发明实施例的对所述第一采样序列进行插值迭代,得到第二采样序列的方法步骤流程图;
图4为本发明示例的第一采样序列的第一次插值迭代的插值位置示意图;
图5为本发明示例中的删除原序列的结构示意图;
图6为本发明实施例的投切电流的波形;
图7为图6的频域变换后的波形;
图8为本发明实施例投切电流+台区负荷电流波形;
图9为本发明实施例去极值滤波并迭代10次的投切电流+台区负荷电流波形;
图10为本发明实施例的台区采样信号的处理系统结构示意图;
图11为本发明实施例的台区采样信号的处理装置。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
以下通过特定的具体实例说明本实施例的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本实施例的其他优点与功效。显然,所描述的实施例仅仅是本实施例一部分实施例,而不是全部的实施例。本实施例还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本实施例的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。基于本实施例中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本实施例保护的范围。
下面结合图2具体介绍本实施例的方法。
如图2所示,本实施例提供的一种台区采样信号的处理方法,包括步骤:
S101,获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;
本步骤采集的第一采样信号,可以是采集到的台区线路上的原始的电压或电流信号,对于纯阻线性负载来说,电流与电压是完全同步的,都是正弦波,但是由于居民用户的普遍存在变频空调、电磁炉等含有电力电子器件的用电器。同时还有大量的非线性负载,如电机等,使得电流并不是标准的正弦波信号。电网的频率特性基本是恒定的,其相位依然存在周期性变化,因此在台区线路上采集的原始电流或电压信号可看做“类正弦波”或“类余弦波”信号。“类正弦波”与“类余弦波”由于只是相位相差90度,但是在处理本方法的原理上基本相同,可选择其一进行说明,另外,为了计算方便,本实施例以电流信号作为处理对象进行处理,对于电力信号来说,其真实信号波形更加接近正弦函数或余弦函数,特别对于用户负载是纯阻性的,其电流波形几乎完全是正弦波。
由于在采样时,按照一定的采样频率进行采样,因此第一采样信号包括离散的第一采样序列,当第一采样信号为原始的电信号时,由于大量噪声的存在,使得采样波形存在毛刺型突发干扰,当按照一定的采样频率进行采样时,采样序列中也会伴随着噪声干扰信号的存在,当不滤除出现的干扰信号时,会对后续的迭代处理产生影响,影响其提取的频谱信息,所以在一些实施例中,作为比较优选的方式,首先针对第一采样序列中符合一定预设条件的采样点进行正弦函数或余弦函数的赋值,所述符合一定预设条件的采样点包括所述第一采样序列中的各采样点的采样值高于预设峰值的采样点。
在一些实施例中,所述对所述第一采样序列中符合一定预设条件的采样点进行重新赋值的方法包括:
判断所述第一采样序列中是否出现采样值大于第一峰值的采样点;
若所述第一采样序列中出现采样值大于第一峰值的采样点,对采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值。
在另外一些可实现的实施例中,所述第一峰值的获取方法包括:
根据公式:
获取所述第一峰值,其中,A表示第一峰值,Xi为所述第一采样序列中第i个采样点的采样数值,n为所述第一采样序列中采样点的数量。
需要说明的是,所述第一峰值并不限于某一个固定的数值,考虑到误差的存在,第一峰值的取值范围可以允许由±%5的误差范围。
可以从上式中获取第一峰值A的原因在于:例如连续的电流信号的有效值可表示为:
其峰值为:
那么当采样值是离散的采样序列时,第一峰值可由式(1)获得。
在对第一采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值后,再判断所述第一采样序列中每1~1.5个工频周期中波峰和波谷的数量,当有且只有一个波峰和波谷,数据不作处理,当出现两个或者两个以上的波峰和波谷,则对出现的所有波峰和波谷对应的采样点,采用正弦函数或余弦函数插值法重新赋值,将赋值后的采样序列作为第一采样序列,第一采样序列可用{H1,H2,H3…Hn}表示,其中,n≥1。
S102,对所述第一采样序列进行插值迭代,得到第二采样序列;
如图3所示,在一些可实现的实施例中,对所述第一采样序列进行插值迭代,得到第二采样序列可通过以下方法步骤实现:
步骤1,对所述第一采样序列中的相邻数值之间进行正弦函数或余弦函数插值,得到插值后的采样序列;
步骤2,在所述插值后的采样序列中删除原采样序列,得到迭代后的采样序列;
步骤3,判断步骤1的执行次数,当所述步骤1的执行次数少于m,则执行步骤4;当所述步骤1的执行次数为m,则执行步骤5,其中,2≤m≤20;
步骤4,将所述迭代后的采样序列作为所述第一采样序列,从步骤1开始重新执行;
步骤5,将所述迭代后的采样序列作为第二采样序列。
为了更直观形象的说明上述方法步骤,下面以一示例进行具体说明:
如图4所示,设第一采样序列为{H1,H2,H3,H4,H5},此处举例采用5个采样点。
采用正弦函数插值,
得到新的采样序列:{H1,H’1,H2,H’2,H3,H’3,H4,H’4,H5},然后再删除原采样序列,如图5所示。
得到完成了一次迭代的采样序列:{H’1,H’2,H’3,H’4}。
通过此方法,实现了一次插值迭代运算,按照上述方法,针对新的序列执行插值迭代,同时,假设设定的执行次数m为3次,判断执行次数是否少于3,若小于3,则将新的序列{H’1,H’2,H’3,H’4}作为起始状态的采样序列,继续按照上述方法插值迭代,当执行次数为3,则直接将此时的采样序列作为第二采样序列。
需要说明的是,每一次插值迭代后的采样信号均呈“类正弦波”信号。
还需要说明的是,上述插值迭代的方法,每迭代一次,每两两个采样点的数学关系就越逼近正弦函数,进而使得采样频率与信号频率趋于同步。实际应用中,数据样本也并不能无限的迭代,迭代次数过多,容易出现破坏目标信号。基于台区特性电流提取的应用,综合该应用的目标信号频率、投切电流幅值以及设备的采样率等参数,经过台区现场调试,在其迭代10次左右的时候可以达到最佳效果。
S103,对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
在一些实施例中,通过离散傅里叶变换算法提取第k次频率fk次频域分量的计算公式为:
式中:
ak、bk、ck分别表示频率为fk的谐波电流的实部、虚部和模值;
N表示参与离散傅里叶变换运算的采样点数;
n表示采样点序号;
上述方法还包括依据所述频谱信息,确定所述第一采样信号是否包括台区设备发送的特征信号,具体包括:
将f1+f2作为第一特征频率,f1-f2作为第二特征频率,其中f1为所述特征信号对应的频率,f2为工频;
根据所述频谱信息中频率与模值的对应关系,确定所述第一特征频率和/或所述第一特征频率所对应的最低信号量;
在所述最低信号量高于第一预设信号量时,确定所述第一采样信号包括台区设备发送的特征信号。由此判断第一采样信号中是否包括台区设备发送的特征信号。
下面举例说明通过离散傅里叶变换后对最后的频点分量的影响。
如图6所示,为测试数据投切电流220mA,投切频率833Hz,台区负荷电流(C相)约180A,电流互感器采集二次电流互感器(600:5)的投切电流的波形,图6中,线型1,线型2和线型3分别表示C通道、B通道和A通道的电流波形,其中,横坐标表示时间,纵坐标表示电流幅值。
图7为图6的频域变换后的波形,横坐标表示频点,纵坐标也表示电流幅值,通过傅里叶变换可以解析到第一特征频率883±5Hz和第二特征频率783±5Hz两个频点的波形。图8为其投切电流+台区负荷电流波形,横坐标与纵坐标的所表示的内容与图6相同,线型1,线型2和线型3分别表示C通道、B通道和A通道的电流波形。图8中频域变换后的目标频点的分量如下表1:
表1
作为对比,图9为去极值滤波并迭代10次的投切电流+台区负荷电流波形,频域变换后的目标频点的分量如下表2:
表2
投切次数 | 783Hz | 883Hz |
1 | 286271 | 211273 |
2 | 334879 | 186789 |
3 | 398662 | 230996 |
4 | 387737 | 224032 |
5 | 364833 | 177962 |
6 | 352860 | 207164 |
对比表1和表2的数据,通过本实施例方法采样后目标信号量(783±5Hz和883±5Hz)(电流幅值)总体要大于没有采样本实施例方法的。本实施例方法前的数据频域变换后其主要表现目标频点的分量(783±5Hz和883±5Hz)并不是峰值,峰值更多是出现在目标附近的频段,极端情况下还会出现频谱泄露现象。通过本发明解决了采样频率和信号频率不能完全同步的问题,避免了频谱泄露的现象,进一步地可以消除测试误差,提高谐波的分析精度。
如图10所示,本实施例还提供了一种台区采样信号的处理系统,包括
采样信号获取单元201,用于获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;
采样信号处理单元202,用于对所述第一采样序列进行插值迭代,得到第二采样序列;
频谱信息提取单元203,用于对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
需要说明的是,上述实施例提供的系统在执行上述方法时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成上述的全部或者部分功能。另外,上述实施例提供的系统与方法实施例属于同一构思,其体现实现过程详见方法实施例,这里不再赘述。
上述本申请实施例序号仅仅为了描述,不代表实施例的优劣。
本申请实施例还提供了一种计算机存储介质,所述计算机存储介质存储有多条指令,所述指令适于有处理器加载并执行上述方法步骤。
具体执行过程可以参见图2所示实施例的具体说明,在此不进行赘述。
本申请实施例还提供了一种计算机设备,包括处理器和存储器,所述存储器存储有计算机程序,所述计算机程序适于由所述处理器加载并执行上述方法。
请参见图11,为本申请实施例提供了台区采样信号的处理装置。
处理装置300可以包括:至少一个处理器301,至少一个网络接口304,用户接口303,存储器305,至少一个通信总线302。
其中,通信总线302用于实现这些组件之间的连接通信。
其中,用户接口303可以包括触摸屏和摄像头(Camera)。
其中,网络接口304可选的可以包括标准的有线接口、无线接口(如WI-FI接口)。
其中,处理器301可以包括一个或者多个处理核心。处理器301利用各种接口和线路连接整个多租户模式数据采集系统300内的各个部分,通过运行或执行存储在存储器305内的指令、程序、代码集或指令集,以及调用存储在存储器305内的数据,执行多租户模式数据采集系统300的各种功能和处理数据。可选的,处理器301可以采用数字信号处理(Digital Signal Processing,DSP)、现场可编程门阵列(Field-Programmable GateArray,FPGA)、可编程逻辑阵列(Programmable Logic Array,PLA)中的至少一种硬件形式来实现。处理器301可集成中央处理器(Central Processing Unit,CPU)、图像处理器(Graphics Processing Unit,GPU)和调制解调器等中的一种或几种的组合。其中,CPU主要处理操作系统、用户界面和应用程序等;GPU用于负责显示屏所需要显示的内容的渲染和绘制;调制解调器用于处理无线通信。可以理解的是,上述调制解调器也可以不集成到处理器301中,单独通过一块芯片进行实现。
其中,存储器305可以包括随机存储器(Random Access Memory,RAM),也可以包括只读存储器(Read-Only Memory)。可选的,该存储器305包括非瞬时性计算机可读介质(non-transitory computer-readable storage medium)。存储器305可用于存储指令、程序、代码、代码集或指令集。存储器305可包括存储程序区和存储数据区,其中,存储程序区可存储用于实现操作系统的指令、用于至少一个功能的指令(比如触控功能、声音播放功能、图像播放功能等)、用于实现上述各个方法实施例的指令等;存储数据区可存储上面各个方法实施例中涉及到的数据等。存储器305可选的还可以是至少一个位于远离前述处理器301的存储装置。如图11所示,作为一种计算机存储介质的存储器305中可以包括操作系统、网络通信模块、用户接口模块以及应用程序。在图11所示的系统300中,用户接口303主要用于为用户提供输入的接口,获取用户输入的数据;而处理器301可以用于调用存储器305中存储的应用程序,并具体执行如图2所示的方法,具体过程可参照图2所示,此处不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本实施例的具体实施方式,但本实施例的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本实施例揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本实施例的保护范围之内。因此,本实施例的保护范围应以权利要求的保护范围为准。
Claims (10)
1.一种台区采样信号的处理方法,其特征在于,包括:
获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;
对所述第一采样序列进行插值迭代,得到第二采样序列;
对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
2.如权利要求1所述的处理方法,其特征在于,所述对所述第一采样序列进行插值迭代,得到第二采样序列,包括:
步骤1,对所述第一采样序列中的相邻数值之间进行正弦函数或余弦函数插值,得到插值后的采样序列;
步骤2,在所述插值后的采样序列中删除原采样序列,得到迭代后的采样序列;
步骤3,判断步骤1的执行次数,当所述步骤1的执行次数少于m,则执行步骤4;当所述步骤1的执行次数为m,则执行步骤5,其中,2≤m≤20;
步骤4,将所述迭代后的采样序列作为所述第一采样序列,从步骤1开始重新执行;
步骤5,将所述迭代后的采样序列作为第二采样序列。
3.如权利要求1所述的处理方法,其特征在于,对所述第一采样序列进行插值迭代之前,还包括对所述第一采样序列进行去极值滤波;
所述去极值滤波的方法包括:
对所述第一采样序列中符合一定预设条件的采样点进行重新赋值;
所述符合一定预设条件的采样点包括:所述第一采样序列中的各采样点的采样值高于预设峰值的采样点。
4.如权利要求3所述的处理方法,其特征在于:所述对所述第一采样序列中符合一定预设条件的采样点进行重新赋值的方法包括:
判断所述第一采样序列中是否出现采样值大于第一峰值的采样点;
若所述第一采样序列中出现采样值大于第一峰值的采样点,对采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值,将赋值后的采样序列作为第一采样序列。
6.如权利要求4所述的处理方法,其特征在于:所述若所述第一采样序列中出现采样值大于第一峰值的采样点,对采样序列中大于第一峰值的采样点,采用正弦函数或余弦函数插值法重新赋值之后还包括:
判断所述第一采样序列中每1~1.5个工频周期中波峰和波谷的数量;
当有且只有一个波峰和波谷,数据不作处理;当出现两个或者两个以上的波峰和波谷,则对出现的所有波峰和波谷对应的采样点,采用正弦函数或余弦函数插值法重新赋值,将赋值后的采样序列作为第一采样序列。
7.如权利要求1所述的处理方法,其特征在于:所述方法还包括:
依据所述频谱信息,确定所述第一采样信号是否包括台区设备发送的特征信号。
8.一种台区采样信号的处理系统,其特征在于,包括
采样信号获取单元,用于获取台区线路的第一采样信号,所述第一采样信号包括离散的第一采样序列;采样信号处理单元,用于对所述第一采样序列进行插值迭代,得到第二采样序列;
频谱信息提取单元,用于对所述第二采样序列进行傅里叶变换,得到所述第一采样信号对应的频谱信息。
9.一种计算机存储介质,其特征在于,所述计算机存储介质存储有多条指令,所述指令适于有处理器加载并执行如权利要求1-7中任一项所述的方法步骤。
10.一种计算机设备,其特征在于,包括处理器和存储器,所述存储器存储有计算机程序,所述计算机程序适于由所述处理器加载并执行如权利要求1-7中任一项所述的方法步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211130264.XA CN115598416B (zh) | 2022-09-16 | 2022-09-16 | 台区采样信号的处理方法、系统、存储介质及计算机设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211130264.XA CN115598416B (zh) | 2022-09-16 | 2022-09-16 | 台区采样信号的处理方法、系统、存储介质及计算机设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115598416A true CN115598416A (zh) | 2023-01-13 |
CN115598416B CN115598416B (zh) | 2024-01-30 |
Family
ID=84843880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211130264.XA Active CN115598416B (zh) | 2022-09-16 | 2022-09-16 | 台区采样信号的处理方法、系统、存储介质及计算机设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115598416B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116094637A (zh) * | 2023-04-13 | 2023-05-09 | 成都德芯数字科技股份有限公司 | 一种中波调幅广播的应急广播指令信号识别方法及其系统 |
CN116660612A (zh) * | 2023-07-31 | 2023-08-29 | 青岛鼎信通讯科技有限公司 | 一种特征电流检测方法 |
Citations (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4928251A (en) * | 1987-04-27 | 1990-05-22 | Hewlett-Packard Company | Method and apparatus for waveform reconstruction for sampled data system |
JP2000180484A (ja) * | 1998-12-18 | 2000-06-30 | Nissin Electric Co Ltd | 高調波測定装置 |
US6114859A (en) * | 1997-07-14 | 2000-09-05 | Nissin Electric Co., Ltd. | Harmonic characteristic measuring method and harmonic characteristic measuring apparatus |
JP2002022783A (ja) * | 2000-07-07 | 2002-01-23 | Chubu Electric Power Co Inc | 高調波特性測定方法及び高調波特性測定装置 |
CN101126772A (zh) * | 2007-09-11 | 2008-02-20 | 电子科技大学 | 一种示波器的高速信号重构方法 |
CN101506668A (zh) * | 2006-08-07 | 2009-08-12 | 西门子公司 | 特别用于能量传输系统的保护或控制装置的采样电流或电压变化并形成采样值的方法 |
JP2009264753A (ja) * | 2008-04-22 | 2009-11-12 | Yokogawa Electric Corp | 高調波成分測定装置 |
CN102243272A (zh) * | 2011-04-01 | 2011-11-16 | 重庆大学 | 一种高精度采样数据同步的谐波分析方法 |
CN102520245A (zh) * | 2011-09-28 | 2012-06-27 | 天津大学 | 基于三次样条插值波形重构的微网谐波及间谐波分析方法 |
CN102638265A (zh) * | 2012-04-28 | 2012-08-15 | 重庆大学 | 一种机械系统内部信号获取方法 |
JP2012154763A (ja) * | 2011-01-26 | 2012-08-16 | Hioki Ee Corp | 測定装置および測定方法 |
CN102707132A (zh) * | 2012-05-21 | 2012-10-03 | 国电南瑞科技股份有限公司 | 一种数字化保护测控装置异常采样值动态处理方法 |
CN103399204A (zh) * | 2013-07-31 | 2013-11-20 | 武汉大学 | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 |
CN103941090A (zh) * | 2014-04-22 | 2014-07-23 | 国家电网公司 | 基于谱线能量插值的谐波测量方法 |
CN104808038A (zh) * | 2015-04-03 | 2015-07-29 | 王顺江 | 一种软仿真插值提升交流采样频率的方法 |
CN105044459A (zh) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
CN106526312A (zh) * | 2016-10-10 | 2017-03-22 | 清华大学 | 基于r‑v (ⅲ)窗fft双峰插值的电能计量方法 |
CN106932642A (zh) * | 2017-03-13 | 2017-07-07 | 威胜集团有限公司 | 电力谐波分析方法 |
CN106998236A (zh) * | 2017-02-07 | 2017-08-01 | 中国人民解放军国防科学技术大学 | 一种基于插值滤波的反馈型符号定时同步装置及方法 |
CN108132399A (zh) * | 2017-12-15 | 2018-06-08 | 国家电网公司 | 一种提高数字化变电站电能质量分析精度的简化插值算法 |
CN108414001A (zh) * | 2017-11-09 | 2018-08-17 | 中国航空工业集团公司北京长城计量测试技术研究所 | 非均匀采样正弦波形失真度的确定方法 |
CN110618338A (zh) * | 2019-11-05 | 2019-12-27 | 广东电网有限责任公司 | 一种电能质量瞬态测量方法、装置和设备 |
CN110888365A (zh) * | 2019-11-30 | 2020-03-17 | 国网辽宁省电力有限公司锦州供电公司 | 电网系统非同步采样基波数据同步化方法 |
CN111984920A (zh) * | 2020-08-31 | 2020-11-24 | 广东电网有限责任公司广州供电局 | 次/超同步谐波参数识别方法、装置、设备和介质 |
CN112748284A (zh) * | 2020-12-30 | 2021-05-04 | 清华大学 | 电力系统同步波形与广谱相量测量方法及装置 |
CN112803896A (zh) * | 2021-01-07 | 2021-05-14 | 配天机器人技术有限公司 | 一种正弦信号的生成方法、装置、设备及介质 |
CN113078977A (zh) * | 2021-03-18 | 2021-07-06 | 国网江苏省电力有限公司常州供电分公司 | 数据的同步方法、同步装置和存储介质 |
CN113255579A (zh) * | 2021-06-18 | 2021-08-13 | 上海建工集团股份有限公司 | 一种施工监测异常采集数据自动识别与处理的方法 |
CN114089034A (zh) * | 2021-10-29 | 2022-02-25 | 国网浙江省电力有限公司杭州供电公司 | 一种基于波形幅值计算反馈分析的异常采样点处理方法及装置 |
CN114280366A (zh) * | 2021-12-21 | 2022-04-05 | 中国航天科工集团八五一一研究所 | 基于改进型频率插值算法的正弦信号频率估计方法 |
CN114552584A (zh) * | 2022-01-14 | 2022-05-27 | 珠海多创科技有限公司 | 一种低压电网台区拓扑识别系统及识别方法 |
-
2022
- 2022-09-16 CN CN202211130264.XA patent/CN115598416B/zh active Active
Patent Citations (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4928251A (en) * | 1987-04-27 | 1990-05-22 | Hewlett-Packard Company | Method and apparatus for waveform reconstruction for sampled data system |
US6114859A (en) * | 1997-07-14 | 2000-09-05 | Nissin Electric Co., Ltd. | Harmonic characteristic measuring method and harmonic characteristic measuring apparatus |
JP2000180484A (ja) * | 1998-12-18 | 2000-06-30 | Nissin Electric Co Ltd | 高調波測定装置 |
JP2002022783A (ja) * | 2000-07-07 | 2002-01-23 | Chubu Electric Power Co Inc | 高調波特性測定方法及び高調波特性測定装置 |
CN101506668A (zh) * | 2006-08-07 | 2009-08-12 | 西门子公司 | 特别用于能量传输系统的保护或控制装置的采样电流或电压变化并形成采样值的方法 |
CN101126772A (zh) * | 2007-09-11 | 2008-02-20 | 电子科技大学 | 一种示波器的高速信号重构方法 |
JP2009264753A (ja) * | 2008-04-22 | 2009-11-12 | Yokogawa Electric Corp | 高調波成分測定装置 |
JP2012154763A (ja) * | 2011-01-26 | 2012-08-16 | Hioki Ee Corp | 測定装置および測定方法 |
CN102243272A (zh) * | 2011-04-01 | 2011-11-16 | 重庆大学 | 一种高精度采样数据同步的谐波分析方法 |
CN102520245A (zh) * | 2011-09-28 | 2012-06-27 | 天津大学 | 基于三次样条插值波形重构的微网谐波及间谐波分析方法 |
CN102638265A (zh) * | 2012-04-28 | 2012-08-15 | 重庆大学 | 一种机械系统内部信号获取方法 |
CN102707132A (zh) * | 2012-05-21 | 2012-10-03 | 国电南瑞科技股份有限公司 | 一种数字化保护测控装置异常采样值动态处理方法 |
CN103399204A (zh) * | 2013-07-31 | 2013-11-20 | 武汉大学 | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 |
CN103941090A (zh) * | 2014-04-22 | 2014-07-23 | 国家电网公司 | 基于谱线能量插值的谐波测量方法 |
CN104808038A (zh) * | 2015-04-03 | 2015-07-29 | 王顺江 | 一种软仿真插值提升交流采样频率的方法 |
CN105044459A (zh) * | 2015-07-21 | 2015-11-11 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
CN106526312A (zh) * | 2016-10-10 | 2017-03-22 | 清华大学 | 基于r‑v (ⅲ)窗fft双峰插值的电能计量方法 |
CN106998236A (zh) * | 2017-02-07 | 2017-08-01 | 中国人民解放军国防科学技术大学 | 一种基于插值滤波的反馈型符号定时同步装置及方法 |
CN106932642A (zh) * | 2017-03-13 | 2017-07-07 | 威胜集团有限公司 | 电力谐波分析方法 |
CN108414001A (zh) * | 2017-11-09 | 2018-08-17 | 中国航空工业集团公司北京长城计量测试技术研究所 | 非均匀采样正弦波形失真度的确定方法 |
CN108132399A (zh) * | 2017-12-15 | 2018-06-08 | 国家电网公司 | 一种提高数字化变电站电能质量分析精度的简化插值算法 |
CN110618338A (zh) * | 2019-11-05 | 2019-12-27 | 广东电网有限责任公司 | 一种电能质量瞬态测量方法、装置和设备 |
CN110888365A (zh) * | 2019-11-30 | 2020-03-17 | 国网辽宁省电力有限公司锦州供电公司 | 电网系统非同步采样基波数据同步化方法 |
CN111984920A (zh) * | 2020-08-31 | 2020-11-24 | 广东电网有限责任公司广州供电局 | 次/超同步谐波参数识别方法、装置、设备和介质 |
CN112748284A (zh) * | 2020-12-30 | 2021-05-04 | 清华大学 | 电力系统同步波形与广谱相量测量方法及装置 |
CN112803896A (zh) * | 2021-01-07 | 2021-05-14 | 配天机器人技术有限公司 | 一种正弦信号的生成方法、装置、设备及介质 |
CN113078977A (zh) * | 2021-03-18 | 2021-07-06 | 国网江苏省电力有限公司常州供电分公司 | 数据的同步方法、同步装置和存储介质 |
CN113255579A (zh) * | 2021-06-18 | 2021-08-13 | 上海建工集团股份有限公司 | 一种施工监测异常采集数据自动识别与处理的方法 |
CN114089034A (zh) * | 2021-10-29 | 2022-02-25 | 国网浙江省电力有限公司杭州供电公司 | 一种基于波形幅值计算反馈分析的异常采样点处理方法及装置 |
CN114280366A (zh) * | 2021-12-21 | 2022-04-05 | 中国航天科工集团八五一一研究所 | 基于改进型频率插值算法的正弦信号频率估计方法 |
CN114552584A (zh) * | 2022-01-14 | 2022-05-27 | 珠海多创科技有限公司 | 一种低压电网台区拓扑识别系统及识别方法 |
Non-Patent Citations (6)
Title |
---|
F. ZHOU 等: "Time-domain quasi-synchronous sampling algorithm for harmonic analysis", PROCEEDINGS OF 14TH INTERNATIONAL CONFERENCE ON HARMONICS AND QUALITY OF POWER - ICHQP 2010 * |
O V S R VARAPRASAD 等: "A novel synchronous sampling algorithm for power system harmonic analysis", 2013 ANNUAL IEEE INDIA CONFERENCE(INDICON) * |
周海波, 王学伟, 白杨: "数据插值在谐波分析算法中的应用及其实现", 电测与仪表, no. 08 * |
曾泽昊;余有灵;许维胜;: "基于插值同步算法的频谱泄漏分析与仿真", 系统仿真技术, no. 03 * |
王成悦 等: "基于数据插值法的电网谐波 FFT 分析方法", 内燃机与配件, pages 134 * |
芮白林: "用插值法实现同步采样", 安徽工业大学学报(自然科学版) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116094637A (zh) * | 2023-04-13 | 2023-05-09 | 成都德芯数字科技股份有限公司 | 一种中波调幅广播的应急广播指令信号识别方法及其系统 |
CN116660612A (zh) * | 2023-07-31 | 2023-08-29 | 青岛鼎信通讯科技有限公司 | 一种特征电流检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115598416B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115598416B (zh) | 台区采样信号的处理方法、系统、存储介质及计算机设备 | |
Morsi et al. | Power quality evaluation in smart grids considering modern distortion in electric power systems | |
Chen et al. | Impact of smart metering data aggregation on distribution system state estimation | |
CN110082612B (zh) | 超/特高压换流站直流场电场仿真方法、装置及存储介质 | |
CN110045186B (zh) | 电力系统公共连接点处的系统谐波阻抗估计方法和系统 | |
EP2938034A1 (en) | Method and system for updating predistortion coefficient | |
Altintasi et al. | Power system harmonic and interharmonic estimation using Vortex Search Algorithm | |
Valverde et al. | Integration of open source tools for studying large‐scale distribution networks | |
CN113608023A (zh) | 采样频率自适应的谐波检测方法及系统 | |
Ravindran et al. | Interharmonics in PV systems: a review of analysis and estimation methods; considerations for selection of an apt method | |
CN112051446A (zh) | 一种电力系统宽频测量的混合基fft实现方法及装置 | |
Rodrigues et al. | Experimental results on a wireless wattmeter device for the integration in home energy management systems | |
Luo et al. | A new measurement algorithm for PMU in power system based on all-phase Fourier transform | |
Hsieh et al. | Slant transform applied to electric power quality detection with field programmable gate array design enhanced | |
CN112838670A (zh) | 基于宽频测量装置的电能治理系统及谐波提取方法 | |
CN112305485B (zh) | 电容式电压互感器的谐波电压测量误差修正方法及装置 | |
CN113702739A (zh) | 一种基于小波分解和神经网络的电弧检测方法及装置 | |
CN113552415A (zh) | 一种超高次谐波测量装置及测量方法 | |
CN113158447B (zh) | 一种大步长移频电磁暂态仿真方法及系统 | |
CN113030648B (zh) | 电力电缆故障点位置确定方法、装置及终端设备 | |
JP2012042296A (ja) | 不規則なパルス列の信号処理方法および信号処理装置 | |
Kamble et al. | Estimation of various PQ indices at point of common coupling using empirical wavelet transform | |
CN107944029B (zh) | 一种基于频域特征的时序数据存储的方法和装置 | |
CN109066973B (zh) | 一种智能变电站的可视化方法、系统及终端设备 | |
Lima et al. | Noninvasive monitoring of residential loads |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |