CN112462179A - 一种狭窄带通滤波方法、装置、设备及存储介质 - Google Patents

一种狭窄带通滤波方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN112462179A
CN112462179A CN202011286981.2A CN202011286981A CN112462179A CN 112462179 A CN112462179 A CN 112462179A CN 202011286981 A CN202011286981 A CN 202011286981A CN 112462179 A CN112462179 A CN 112462179A
Authority
CN
China
Prior art keywords
value
transfer function
zero
sampling
target signal
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
Application number
CN202011286981.2A
Other languages
English (en)
Inventor
叶飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Shipbuilding Power Engineering Institute Co Ltd
Original Assignee
China Shipbuilding Power Engineering Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Shipbuilding Power Engineering Institute Co Ltd filed Critical China Shipbuilding Power Engineering Institute Co Ltd
Priority to CN202011286981.2A priority Critical patent/CN112462179A/zh
Publication of CN112462179A publication Critical patent/CN112462179A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明实施例公开了一种狭窄带通滤波方法、装置、设备及存储介质。其中，该方法包括：在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数；将所述新的传递函数通过Z反变换得到差分方程；通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值。本发明实施例提供的技术方案，能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。

Description

一种狭窄带通滤波方法、装置、设备及存储介质
技术领域
本发明实施例涉及滤波技术，尤其涉及一种狭窄带通滤波方法、装置、设备及存储介质。
背景技术
在电力系统发生故障后，继电保护主要是通过对故障后电流和电压信号的基频分量进行提取分析从而能够快速定位发生故障的具体位置。
常用的基频分量提取方法是采用全波傅里叶变换算法和半波傅里叶变换算法。全波傅里叶变换算法的数据窗需要一个周期，导致了响应速度比减慢，无法满足有些继电保护的速动性要求。而半波傅里叶变换算法只需要半个周期的数据窗，虽然响应速度快，但是无法有效滤除直流分量和偶次谐波分量，同时对偶次谐波具有放大影响，对低频分量抑制效果较差，容易导致继电保护中故障的误判断，造成故障扩大化。
发明内容
本发明实施例提供了一种狭窄带通滤波方法、装置、设备及存储介质，能够有效滤除直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。
第一方面，本发明实施例提供了一种狭窄带通滤波方法，该方法包括：
在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
将所述新的传递函数通过Z反变换得到差分方程；
通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
第二方面，本发明实施例提供了一种狭窄带通滤波装置，该装置包括：
传递函数确定模块，用于在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
差分方程确定模块，用于将所述新的传递函数通过Z反变换得到差分方程；
提取值确定模块，用于通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
第三方面，本发明实施例提供了一种计算机设备，该计算机设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现本发明任意实施例所述的狭窄带通滤波方法。
第四方面，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明任意实施例所述的狭窄带通滤波方法。
本发明实施例提供了一种狭窄带通滤波方法、装置、设备及存储介质，首先在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，然后将新的传递函数通过Z反变换得到差分方程，最后通过差分方程对目标信号进行滤波，得到目标信号的基频分量提取值，能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。
附图说明
图1为本发明实施例一提供的一种狭窄带通滤波方法的流程图；
图2为本发明实施例二提供的一种狭窄带通滤波方法的流程图；
图3为本发明实施例三提供的一种狭窄带通滤波装置的结构示意图；
图4为本发明实施例四提供的一种计算机设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种狭窄带通滤波方法的流程图，本实施例可适用于对继电保护过程中电力系统故障后的电流信号或者电压信号进行滤波的情况。本实施例提供的狭窄带通滤波方法可以由本发明实施例提供的狭窄带通滤波装置来执行，该装置可以通过软件和/或硬件的方式实现，并集成在执行本方法的计算机设备中。
参见图1，本实施例的方法包括但不限于如下步骤：
S110，在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数。
其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的。
考虑到要将某次谐波分量完全滤除，需要将该谐波分量设置为零点，要使某次谐波分量输出最大，需要将该谐波分量设置为极点，由此在相位角为0时，通过零点公式可以确定出传递函数中第一个零点的值，在相位角为π时通过零点公式可以确定出传递函数中第二个零点的值，然后通过极点公式确定出传递函数中两个极点的值，得到新的传递函数。
具体的，传递函数的推导过程如下所示：
(1)、一个n阶线性时不变系统的常系数线性差分方程可用下式表示：
Figure BDA0002782683940000041
其中，x表示自变量，y表示因变量，n表示大于或者等于1的整数，a0，a1…，aN为常数，b1，b2，…，bM为常数。
(2)、对上述差分方程进行Z变换可以得到下式：
Y(Z)＝a0X(Z)+a1X(Z)Z-1+…+aNX(Z)Z-N+b1Y(Z)Z-1
+b2Y(Z)Z-2+…+bMY(Z)Z-M
其中，Z表示Z变换的算子。
(3)、对上式进行合并同类项可以得到下式：
Y(Z)(1-b1Z-1-b2Z-2-…-bMZ-M)＝X(Z)(a0+a1Z-1+a2Z-2+…+aNZ-N)
(4)获取传递函数：
Figure BDA0002782683940000051
其中，Ci表示传递函数的零点，Dk表示传递函数的极点，A表示增益。
可选的，所述零点公式和极点公式为：
Ci＝ejωTs
Dk＝RejωTs
其中，Ci表示零点，i表示零点的个数，i＝1，2，Dk表示极点，k表示极点的个数，k＝1，2，ω表示角频率，T表示采样周期，s表示拉普拉斯变换算子，R为大于零且小于1的常数；
相应的，所述第一个零点的值为1，所述第二个零点的值为-1，第一个极点的值为RejωTs，第二个极点的值为Re-jωTs。
具体的，当相位角为0(即ωT＝0)时，通过零点公式可以计算得到C1＝1；当相位角为π(即ωT＝π)时，通过零点公式可以计算得到C2＝-1；通过极点点公式可以计算得到D1＝RejωTs，由于第二个极点和第一个极点是共轭关系，由此可以得到D2＝Re-jωTs。
新的传递函数为：
Figure BDA0002782683940000061
S120，将所述新的传递函数通过Z反变换得到差分方程。
在得到了新的传递函数之后，将新的传递函数通过Z反变换可以得到差分方程。
可选的，所述差分方程的表达式为：
y(n+2)＝Ax(n+2)-Ax(n)+2Rcos(ωTs)y(n+1)-R2y(n)
其中，x表示自变量，y表示因变量，n表示大于或者等于1的整数。
具体的，将
Figure BDA0002782683940000062
进行Z反变换就可以得到上面的差分方程表达式。
S130，通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值。
所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
通过S120中的差分方程对继电保护过程中电力系统故障后的电流信号或者电压信号进行滤波，能够得到目标信号对应的基频分量提取值。
进一步的，所述R的取值范围为：0.99<R<1。
具体的，当R＝0.9987时滤波效果较佳，幅值响应较大，能够较大的保留基频分量，且低频和高频幅值响应较小，能够对低频和高频分量进行良好的抑制。
本实施例提供的技术方案，首先在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，然后将新的传递函数通过Z反变换得到差分方程，最后通过差分方程对目标信号进行滤波，得到目标信号的基频分量提取值，能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。
实施例二
图2为本发明实施例二提供的一种狭窄带通滤波方法的流程图。本发明实施例是在上述实施例的基础上进行优化。可选的，本实施例对得到目标信号的基频分量提取值的过程进行详细的解释说明。
参见图2，本实施例的方法包括但不限于如下步骤：
S210，在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数。
S220，将所述新的传递函数通过Z反变换得到差分方程。
S230，对目标信号进行采样，得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值。
其中，所述预设个数为偶数。
假设采样周期T＝20ms，对目标信号进行采样，可以得到预设个数(设为N)的采样点以及与采样点个数相同的目标信号的采样值，即x(1),x(2),…,x(N)。
S240，将所述差分方程和半波傅里叶算法结合，对所述采样值进行处理，得到所述基频分量提取值。
由于初始值会影响对采样值进行处理的速度，为了能使速度加快，此时可以将差分方程和半波傅里叶算法结合，对采样值进行处理，从而得到基频分量提取值。
可选的，所述将所述差分方程和半波傅里叶算法结合，对所述采样值进行处理，得到所述基频分量提取值，可以具体包括：针对第一采样点进行半波傅里叶计算，得到第一初始值；针对第二采样点进行半波傅里叶计算，得到第二初始值，所述第一采样点和第二采样点位于所有采样点的中间，且所述第二采样点位于所述第一采样点之后；将第一初始值、第二初始值以及所述目标信号的采样值代入所述差分方程进行计算，得到所述预设个数的点数据；对所述预设个数的点数据进行半波傅里叶计算，得到所述基频分量提取值。
具体的，由于采样点个数为N，首先针对第一采样点(即第
Figure BDA0002782683940000081
点)进行半波傅里叶计算，可以得到第一初始值，记为y1；接着针对第二采样点(即第
Figure BDA0002782683940000082
+1点)进行半波傅里叶计算，可以得到第二初始值，记为y2；然后将第一初始值、第二初始值以及目标信号的采样值(x(1),x(2),…,x(N))代入差分方程的公式中进行计算，得到y3，y4，…，yN；最后对这N个数据进行半波傅里叶计算，就能够得到基频分量提取值。
本发明实施例通过将差分方程和半波傅里叶算法结合，得到基频分量提取值，不仅收敛速度快，而且能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求。
进一步的，所述针对第一采样点进行半波傅里叶计算，得到第一初始值，可以具体包括：在确定采样时间大于采样周期的一半时，针对第一采样点进行半波傅里叶计算，得到第一初始值。
具体的，由于半波傅里叶算法需要满足采样时间大于采样周期的一半，因此在确定采样时间大于采样周期的一半时，针对第一采样点进行半波傅里叶计算，得到第一初始值。
本发明实施例通过在确定采样时间大于采样周期的一半时，计算第一初始值，使得得到的第一初始值更为准确，从而加快基频分量的提取过程。
本实施例提供的技术方案，首先在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，接着将新的传递函数通过Z反变换得到差分方程，然后对目标信号进行采样，得到预设个数的采样点以及与采样点个数相同的目标信号的采样值，最后将差分方程和半波傅里叶算法结合，对采样值进行处理，得到基频分量提取值，不仅能够加快基频分量的提取过程，而且能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。
实施例三
图3为本发明实施例三提供的一种狭窄带通滤波装置的结构示意图，如图3所示，该装置可以包括：
传递函数确定模块310，用于在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
差分方程确定模块320，用于将所述新的传递函数通过Z反变换得到差分方程；
提取值确定模块330，用于通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
本实施例提供的技术方案，首先在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，然后将新的传递函数通过Z反变换得到差分方程，最后通过差分方程对目标信号进行滤波，得到目标信号的基频分量提取值，能够有效滤除目标信号的直流分量和偶次谐波分量，满足继电保护的速动性要求，从而避免故障的误判断。
进一步的，所述零点公式和极点公式为：
Ci＝ejωTs
Dk＝RejωTs
其中，Ci表示零点，i表示零点的个数，i＝1，2，Dk表示极点，k表示极点的个数，k＝1，2，ω表示角频率，T表示采样周期，s表示拉普拉斯变换算子，R为大于零且小于1的常数；
相应的，所述第一个零点的值为1，所述第二个零点的值为-1，第一个极点的值为RejωTs，第二个极点的值为Re-jωTs。
进一步的，上述狭窄带通滤波装置，还可以包括：
目标信号采样模块，用于对目标信号进行采样，得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值，所述预设个数为偶数；
分量提取值确定模块，用于将所述差分方程和半波傅里叶算法结合，对所述采样值进行处理，得到所述基频分量提取值。
进一步的，上述分量提取值确定模块，可以具体用于：
针对第一采样点进行半波傅里叶计算，得到第一初始值；
针对第二采样点进行半波傅里叶计算，得到第二初始值，所述第一采样点和第二采样点位于所有采样点的中间，且所述第二采样点位于所述第一采样点之后；
将第一初始值、第二初始值以及所述目标信号的采样值代入所述差分方程进行计算，得到所述预设个数的点数据；
对所述预设个数的点数据进行半波傅里叶计算，得到所述基频分量提取值。
进一步的，上述狭窄带通滤波装置，还可以包括：
第一初始值确定模块，用于在确定采样时间大于采样周期的一半时，针对第一采样点进行半波傅里叶计算，得到第一初始值。
进一步的，所述差分方程的表达式为：
y(n+2)＝Ax(n+2)-Ax(n)+2Rcos(ωTs)y(n+1)-R2y(n)
其中，x表示自变量，y表示因变量，n表示大于或者等于1的整数。
进一步的，所述R的取值范围为：0.99<R<1。
本实施例提供的狭窄带通滤波装置可适用于上述任意实施例提供的狭窄带通滤波方法，具备相应的功能和有益效果。
实施例四
图4为本发明实施例四提供的一种计算机设备的结构示意图，如图4所示，该计算机设备包括处理器410、存储装置420和通信装置430；计算机设备中处理器410的数量可以是一个或多个，图4中以一个处理器410为例；计算机设备中的处理器410、存储装置420和通信装置430可以通过总线或其他方式连接，图4中以通过总线连接为例。
存储装置420作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的狭窄带通滤波方法对应的模块(例如，用于狭窄带通滤波装置装置中的传递函数确定模块310、差分方程确定模块320和提取值确定模块330)。处理器410通过运行存储在存储装置420中的软件程序、指令以及模块，从而执行计算机设备的各种功能应用以及数据处理，即实现上述的狭窄带通滤波方法。
存储装置420可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储装置420可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储装置420可进一步包括相对于处理器410远程设置的存储器，这些远程存储器可以通过网络连接至计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
通信装置430，用于实现服务器之间的网络连接或者移动数据连接。
本实施例提供的一种计算机设备可用于执行上述任意实施例提供的狭窄带通滤波方法，具备相应的功能和有益效果。
实施例五
本发明实施例五还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明任意实施例中的狭窄带通滤波方法，该方法具体包括：
在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
将所述新的传递函数通过Z反变换得到差分方程；
通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的狭窄带通滤波方法中的相关操作。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
值得注意的是，上述狭窄带通滤波装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
以上所述仅为本发明的优选实施例，并不用于限制本发明，对于本领域技术人员而言，本发明可以有各种改动和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种狭窄带通滤波方法，其特征在于，包括：
在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
将所述新的传递函数通过Z反变换得到差分方程；
通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
2.根据权利要求1所述的方法，其特征在于，所述零点公式和极点公式为：
Ci＝ejωTs
Dk＝RejωTs
其中，Ci表示零点，i表示零点的个数，i＝1，2，Dk表示极点，k表示极点的个数，k＝1，2，ω表示角频率，T表示采样周期，s表示拉普拉斯变换算子，R为大于零且小于1的常数；
相应的，所述第一个零点的值为1，所述第二个零点的值为-1，第一个极点的值为Rej ωTs，第二个极点的值为Re-jωTs。
3.根据权利要求2所述的方法，其特征在于，所述通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，包括：
对目标信号进行采样，得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值，所述预设个数为偶数；
将所述差分方程和半波傅里叶算法结合，对所述采样值进行处理，得到所述基频分量提取值。
4.根据权利要求3所述的方法，其特征在于，所述将所述差分方程和半波傅里叶算法结合，对所述采样值进行处理，得到所述基频分量提取值，包括：
针对第一采样点进行半波傅里叶计算，得到第一初始值；
针对第二采样点进行半波傅里叶计算，得到第二初始值，所述第一采样点和第二采样点位于所有采样点的中间，且所述第二采样点位于所述第一采样点之后；
将第一初始值、第二初始值以及所述目标信号的采样值代入所述差分方程进行计算，得到所述预设个数的点数据；
对所述预设个数的点数据进行半波傅里叶计算，得到所述基频分量提取值。
5.根据权利要求4所述的方法，其特征在于，所述针对第一采样点进行半波傅里叶计算，得到第一初始值，包括：
在确定采样时间大于采样周期的一半时，针对第一采样点进行半波傅里叶计算，得到第一初始值。
6.根据权利要求2所述的方法，其特征在于，所述差分方程的表达式为：
y(n+2)＝Ax(n+2)-Ax(n)+2R cos(ωTs)y(n+1)-R2y(n)
其中，x表示自变量，y表示因变量，n表示大于或者等于1的整数。
7.根据权利要求6所述的方法，其特征在于，所述R的取值范围为：0.99<R<1。
8.一种狭窄带通滤波装置，其特征在于，包括：
传递函数确定模块，用于在相位角为0时，通过零点公式确定出传递函数中第一个零点的值，和在相位角为π时通过零点公式确定出传递函数中第二个零点的值，通过极点公式确定出传递函数中两个极点的值，得到新的传递函数，其中，所述两个极点的值是共轭关系，所述传递函数是通过零极点配置法得到的；
差分方程确定模块，用于将所述新的传递函数通过Z反变换得到差分方程；
提取值确定模块，用于通过所述差分方程对目标信号进行滤波，得到所述目标信号的基频分量提取值，所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
9.一种计算机设备，其特征在于，所述计算机设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-7中任一所述的狭窄带通滤波方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-7中任一所述的狭窄带通滤波方法。
CN202011286981.2A 2020-11-17 2020-11-17 一种狭窄带通滤波方法、装置、设备及存储介质 Pending CN112462179A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011286981.2A CN112462179A (zh) 2020-11-17 2020-11-17 一种狭窄带通滤波方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011286981.2A CN112462179A (zh) 2020-11-17 2020-11-17 一种狭窄带通滤波方法、装置、设备及存储介质

Publications (1)

Publication Number Publication Date
CN112462179A true CN112462179A (zh) 2021-03-09

Family

ID=74836337

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011286981.2A Pending CN112462179A (zh) 2020-11-17 2020-11-17 一种狭窄带通滤波方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN112462179A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113595527A (zh) * 2021-07-30 2021-11-02 国光电器股份有限公司 一种滤波参数确定方法、滤波方法及相关装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4896152A (en) * 1989-03-02 1990-01-23 General Electric Company Telemetry system with a sending station using recursive filter for bandwidth limiting
CN101599754A (zh) * 2009-03-26 2009-12-09 东南大学 用于增强不对称二元调制信号的冲击滤波方法
CN102623957A (zh) * 2012-04-01 2012-08-01 积成电子股份有限公司 基于窄带滤波的变数据窗基波相量求取方法
CN103162755A (zh) * 2013-01-31 2013-06-19 西安东风机电有限公司 一种基于改进自适应算法的科里奥利流量计信号跟踪方法
CN103227640A (zh) * 2013-04-16 2013-07-31 苏州汇川技术有限公司 谐波电网相位检测方法及锁相环
CN103576120A (zh) * 2013-11-07 2014-02-12 国家电网公司 三次谐波分量准同步信息传输校验与自愈算法
CN104202052A (zh) * 2014-08-29 2014-12-10 辽宁工程技术大学 一种提高信噪比的Sigma-Delta调制器自适应混合优化方法
CN104639099A (zh) * 2015-01-30 2015-05-20 北京时代民芯科技有限公司 一种产生滤波器传递函数的方法
JP2016171552A (ja) * 2015-03-15 2016-09-23 有限会社シー・アンド・エス国際研究所 デジタルバンドパスフィルタ
CN106451369A (zh) * 2015-08-13 2017-02-22 中国电力科学研究院 一种失灵保护中的截断电流互感器拖尾电流方法
CN106684513A (zh) * 2017-01-05 2017-05-17 南通大学 一种具有谐波抑制功能的平衡式带通滤波器
CN110588939A (zh) * 2019-10-29 2019-12-20 中船动力研究院有限公司 一种混合动力电池装置和混合动力系统

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4896152A (en) * 1989-03-02 1990-01-23 General Electric Company Telemetry system with a sending station using recursive filter for bandwidth limiting
CN101599754A (zh) * 2009-03-26 2009-12-09 东南大学 用于增强不对称二元调制信号的冲击滤波方法
CN102623957A (zh) * 2012-04-01 2012-08-01 积成电子股份有限公司 基于窄带滤波的变数据窗基波相量求取方法
CN103162755A (zh) * 2013-01-31 2013-06-19 西安东风机电有限公司 一种基于改进自适应算法的科里奥利流量计信号跟踪方法
CN103227640A (zh) * 2013-04-16 2013-07-31 苏州汇川技术有限公司 谐波电网相位检测方法及锁相环
CN103576120A (zh) * 2013-11-07 2014-02-12 国家电网公司 三次谐波分量准同步信息传输校验与自愈算法
CN104202052A (zh) * 2014-08-29 2014-12-10 辽宁工程技术大学 一种提高信噪比的Sigma-Delta调制器自适应混合优化方法
CN104639099A (zh) * 2015-01-30 2015-05-20 北京时代民芯科技有限公司 一种产生滤波器传递函数的方法
JP2016171552A (ja) * 2015-03-15 2016-09-23 有限会社シー・アンド・エス国際研究所 デジタルバンドパスフィルタ
CN106451369A (zh) * 2015-08-13 2017-02-22 中国电力科学研究院 一种失灵保护中的截断电流互感器拖尾电流方法
CN106684513A (zh) * 2017-01-05 2017-05-17 南通大学 一种具有谐波抑制功能的平衡式带通滤波器
CN110588939A (zh) * 2019-10-29 2019-12-20 中船动力研究院有限公司 一种混合动力电池装置和混合动力系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘序俨: "数字递归滤波器及其应用", 大地测量与地球动力学, no. 03 *
叶飞;刘观起;胡道徐;张进;纪卉;: "提高速动保护精确性滤波新算法", 华北电力大学学报(自然科学版), no. 04, pages 18 - 22 *
詹圣益;: "嵌入式系统中数字滤波器的分析和设计", 微型电脑应用, no. 08 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113595527A (zh) * 2021-07-30 2021-11-02 国光电器股份有限公司 一种滤波参数确定方法、滤波方法及相关装置
CN113595527B (zh) * 2021-07-30 2023-10-20 国光电器股份有限公司 一种滤波参数确定方法、滤波方法及相关装置

Similar Documents

Publication Publication Date Title
CN108921424B (zh) 一种电力数据异常检测方法、装置、设备及可读存储介质
CN111431758A (zh) 云网络设备的测试方法、装置、存储介质和计算机设备
CN106558314B (zh) 一种混音处理方法和装置及设备
CN111049858B (zh) 一种基于交叉验证的基线扫描漏洞去重方法、装置及设备
CN112462179A (zh) 一种狭窄带通滤波方法、装置、设备及存储介质
EP2819025A1 (en) Method for reducing noise in data-sets of harmonic signals
CN111193633A (zh) 异常网络连接的检测方法及装置
CN109116243B (zh) 一种动力电池的异常检测方法和装置
Antonov et al. Adaptive structural analysis of input signals of digital and relay protection and automation
CN111046836B (zh) 局放信号的滤波去噪及分析方法、系统、设备及存储介质
CN103784164B (zh) 超声信号的预处理方法及系统
CN108075807B (zh) 一种小区标识的检测方法及装置
EP3382978A1 (en) Distributed denial of service analysis
CN117318671B (zh) 一种基于快速傅里叶变换的自适应滤波方法
CN109684937A (zh) 一种基于fft及数学形态法的信号去噪方法及装置
CN105510621A (zh) 一种离心压缩机临界转速获取方法和装置
CN111948454A (zh) 一种抗衰减直流分量干扰的同步相量测量方法及系统
CN113055420A (zh) Https业务识别方法、装置及计算设备
Banjac et al. Some applications of λ-method for obtaining approximations in filter design
CN110874387A (zh) 移动设备识别符共现关系稀疏图的构建方法及装置
CN112067953B (zh) 一种提升变电站局部放电定向精度的方法及系统
CN113553370A (zh) 异常检测方法、装置、电子设备及可读存储介质
Tedesco et al. Discrete-time frequency-locked-loop filters for parameters estimation of sinusoidal signals
CN114062992A (zh) 一种电流互感器一次侧检测过程超声波提取方法及装置
Bank et al. Extraction and visualization of power system interarea oscillatory modes

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