高通滤波器的电磁暂态建模方法、装置及电子设备与流程

1.本发明涉及滤波器的建模技术领域，尤其涉及一种高通滤波器的电磁暂态建模方法、装置及电子设备。


背景技术：

2.通过对电力系统相关元件进行准确的建模和详细、快速的电磁暂态仿真分析，可以准确地掌握电力系统的动态特征，为电力系统的规划设计、建设发展和运行实践提供重要的帮助。
3.高通滤波器是由电容和电阻组成的滤波电路，可以使得超过截止频率的低频信号正常通过，而小于截止频率的高频信号被阻隔、减弱。高通滤波器的电磁暂态模型是电力系统电磁暂态仿真中不可或缺的模型组件。
4.目前，与滤波器建模相关技术中，没有能够对高通滤波器进行快速、准确进行电磁暂态建模的方法。


技术实现要素：

5.基于此，有必要针对上述问题，本发明提出了一种高通滤波器的电磁暂态建模方法、装置及电子设备，可以快速、便捷、准确的对高通滤波器进行电磁暂态建模，有利于电力系统电磁暂态仿真。
6.第一方面，本发明提供了一种高通滤波器的电磁暂态建模方法包括：首先，根据高通滤波器确定高通滤波器的频域模型，其次，根据拉普拉斯反变换和高通滤波器的频域模型，得到高通滤波器的分时段时域积分模型；最后，根据梯形法则和高通滤波器的分时段时域积分模型，得到高通滤波器的分段时域电磁暂态模型。
7.在第一方面的一种可能的实现方式中，该电磁暂态建模方法还包括：
8.构建高通滤波器的分段时域电磁暂态模型补偿分量；
9.利用高通滤波器的分段时域电磁暂态模型补偿分量，对高通滤波器的分段时域积分模型进行补偿，得到高通滤波器的分段时域电磁暂态修正模型。
10.在第一方面的一种可能的实现方式中，上述根据拉普拉斯反变换和高通滤波器的频域模型，得到高通滤波器的分时段时域积分模型包括：
11.利用拉普拉斯反变换对高通滤波器的频域模型进行求解得到高通滤波器的基础时域积分模型；
12.根据冲激函数的特性，将高通滤波器的时域积分模型按照时间变量的取值范围划分为两部分，得到高通滤波器的分段时域积分模型。
13.在第一方面的一种可能的实现方式中，上述根据梯形法则和高通滤波器的分时段时域积分模型，得到高通滤波器的分段时域电磁暂态模型包括：
14.利用梯形法则对所述高通滤波器的分时段时域积分模型进行求解，得到高通滤波器的时域电磁暂态基础模型；
15.根据冲激函数和求极限方法，将高通滤波器的时域电磁暂态基础模型按照时间变量的取值范围划分为三部分，得到高通滤波器的分段时域电磁暂态模型。
16.第二方面，本发明提供了一种建模装置，用于对高通滤波器进行电磁暂态建模，该建模装置包括：频域模型确定模块、时域模型计算模块和电磁暂态模型计算模块，其中，频域模型确定模块，用于根据高通滤波器确定高通滤波器的频域模型；时域模型计算模块，用于根据拉普拉斯反变换和高通滤波器的频域模型，得到高通滤波器的分时段时域积分模型；电磁暂态模型计算模块，用于根据梯形法则和高通滤波器的分时段时域积分模型，得到高通滤波器的分段时域电磁暂态模型。
17.在第二方面的一种可能的实现方式中，该建模装置还包括：电磁暂态模型补偿模块，用于执行以下操作：
18.构建高通滤波器的分段时域电磁暂态模型补偿分量；
19.利用高通滤波器的分段时域电磁暂态模型补偿分量，对高通滤波器的分段时域积分模型进行补偿，得到高通滤波器的分段时域电磁暂态修正模型。
20.第三方面，本发明还提供了一种电子设备，包括处理器和存储器，该存储器上存储有可执行代码，当可执行代码被处理器执行时，使处理器执行如上述第一方面及第一方面中任意一种实现方式中所述的方法。
21.第四方面提供了一种计算机可读存储介质，其上存储有可执行代码，当可执行代码被电子设备的处理器执行时，使处理器执行如上述第一方面及第一方面中任意一种实现方式中所述的方法。
22.本发明提供的技术方案具有以下有益效果：
23.通过拉普拉斯反变换和梯形法则求解高通滤波器的时域电磁暂态基础模型，以实现对高通滤波器的电磁暂态建模，既能快速、便捷的建立高通滤波器的分段时域电磁暂态模型，又从数学角度提供了严谨的求解公式，数学逻辑的严谨性可以确保建立分段时域电磁暂态模型的准确性，有利于电力系统电磁暂态仿真及其实际应用。
24.进一步的，本发明的技术方案，构建并通过高通滤波器的分段时域电磁暂态模型补偿分量，可以对高通滤波器的分段时域积分模型的模型精度进行修正，可以有效地提高修正后得到的高通滤波器的分段时域积分修正模型的实用性。
25.应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。
附图说明
26.通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
27.图1为本发明实施例中提供的一个一阶高通滤波器电路结构示意图；
28.图2为本发明实施例中提供的高通滤波器的电磁暂态建模方法的一个实施例示意图；
29.图3为本发明实施例中提供的高通滤波器的电磁暂态建模方法的另一个实施例示意图；
30.图4为本发明实施例中提供的建模装置的一个结构示意图；
31.图5为本发明实施例中提供的建模装置的另一个结构示意图；
32.图6为本发明实施例中提供的电子设备的一个组成结构示意图。
具体实施方式
33.下面将参照附图更详细地描述本发明的实施方式。虽然附图中显示了本发明的实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
34.在本发明使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
35.应当理解，尽管在本发明可能采用术语“第一”、“第二”、“第三”等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
36.滤波器是一种能使有用频率信号通过而同时抑制无用频率信号的电子装置，常用于信号处理、数据传输和干扰抑制等方面。高通滤波器是容许低于截止频率的信号通过，但高于截止频率的信号不能通过的电子滤波装置。高通滤波器根据谐波被过滤的次数可分为一阶高通滤波器、二阶高通滤波器等，其中几阶代表谐波被过滤几次，例如一阶高通滤波器是将谐波进行一次过滤的高通滤波器。
37.为了便于理解本发明中高通滤波器的的电磁暂态建模方法，下面先以一阶高通滤波器为例，对本发明中的高通滤波器进行说明。
38.图1为本发明实施例中提供的一个一阶高通滤波器电路结构示意图。
39.如图1所示，本发明实施例中提供的一阶高通滤波器电路包括：
40.滤波电阻r，滤波电容c，反馈电阻r1，平衡电阻r2和集成运算放大器(简称集成运放)，其中ui为输入电压，u0为输出电压。
41.其中滤波电阻r一端接地，另一端与集成运放的正极输入端连接；滤波电容c连接在集成运放的正极输入端；反馈电阻r1连接在集成运放的负极输入端和输出端之间；平衡电阻r2一段接地，另一端与集成运放的负极输入端连接。
42.对高通滤波器的电力结构进行基本介绍之后，下面对本发明中高通滤波器的的电磁暂态建模方法进行说明。
43.图2为本发明实施例中提供的高通滤波器的电磁暂态建模方法的一个实施例示意图。
44.如图2所示，本发明实施例中高通滤波器的电磁暂态建模方法包括：
45.201、根据高通滤波器确定高通滤波器的频域模型。
46.图1所示一阶高通滤波器，该高通滤波器的频域模型(即频域中的高通滤波器模型)为：
[0047][0048]
其中，y(s)为高通滤波器的频域输出变量，x(s)为高通滤波器的频域输入变量；s为拉普拉斯算子，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数。
[0049]
202、根据拉普拉斯反变换和高通滤波器的频域模型，得到高通滤波器的分时段时域积分模型。
[0050]
具体可选的，首先利用拉普拉斯反变换对高通滤波器的频域模型进行求解得到高通滤波器的基础时域积分模型；
[0051]
例如，利用拉普拉斯反变换对步骤201中高通滤波器的频域模型进行求解可以得到高通滤波器的基础时域积分模型如下：
[0052][0053]
其次根据冲激函数的特性，将高通滤波器的时域积分模型按照时间变量的取值范围划分为两部分，得到高通滤波器的分段时域积分模型。
[0054]
例如，根据冲激函数的特性，对上述公式所示的基础时域积分模型按照时间变量t的取值范围划分为两部分，可以得到高通滤波器的分段时域积分模型如下：
[0055][0056]
其中，t为时间变量，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数。
[0057]
203、根据梯形法则和高通滤波器的分时段时域积分模型，得到高通滤波器的分段时域电磁暂态模型。
[0058]
具体可选的，首先利用梯形法则对高通滤波器的分时段时域积分模型进行求解，得到高通滤波器的时域电磁暂态基础模型。
[0059]
例如，利用梯形法则对步骤202中求解得到的分时段时域积分模型进行求解得到如下时域电磁暂态基础模型：
[0060][0061]
其次，根据冲激函数和求极限方法，将高通滤波器的时域电磁暂态基础模型按照时间变量的取值范围划分为三部分，得到高通滤波器的分段时域电磁暂态模型。
[0062]
例如，根据冲激函数和求极限方法，将上述公式中所示的时域电磁暂态基础模型按照时间变量t的取值范围划分为三部分，可以得到如下分段时域电磁暂态模型：
[0063][0064]
其中，t为时间变量，δt为时间步长，δt

为时间步长的右极限，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，y(t-δt)为高通滤波器t-δt时刻的时域输出变量，x(t-δt)为高通滤波器t-δt时刻的时域输入变量，x(0)为高通滤波器0时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，δ(t-δt)为t-δt时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数，b为常数。
[0065]
综上所述，本发明实施例中的高通滤波器的电磁暂态建模方法，可以通过拉普拉斯反变换和梯形法求解高通滤波器的分段时域电磁暂态模型，以实现对高通滤波器的电磁暂态建模，既能快速、便捷的建立高通滤波器的分段时域电磁暂态模型，又从数学角度提供了严谨的求解公式，数学逻辑的严谨性可以确保建立分段时域电磁暂态模型的准确性，有利于电力系统电磁暂态仿真及其实际应用。
[0066]
进一步的，在对电磁暂态模型精度要求更高的应用场景中，还可以对上述高通滤波器的分段时域电磁暂态模型的模型精度进行修正，以提高其精度。
[0067]
图3为本发明实施例中提供的高通滤波器的电磁暂态建模方法的另一个实施例示意图。
[0068]
如图3所示，本发明实施例中提供的高通滤波器的电磁暂态建模方法包括：
[0069]
301、获取高通滤波器的分段时域电磁暂态模型。
[0070]
获取通过图2对应实施例中的技术方案得到的高通滤波器的分段时域电磁暂态模型：
[0071][0072]
302、构建高通滤波器的分段时域电磁暂态模型补偿分量。
[0073]
具体的，构建高通滤波器的分段时域电磁暂态模型补偿分量为：
[0074][0075]
其中，t为时间变量，δt为时间步长，δt

为时间步长的右极限，y’(t)为高通滤波器t时刻的分段时域输出变量补偿分量，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数。
[0076]
303、利用高通滤波器的分段时域电磁暂态模型补偿分量，对高通滤波器的分段时域积分模型进行补偿，得到高通滤波器的分段时域电磁暂态修正模型。
[0077]
通过步骤302中构建的分段时域电磁暂态模型补偿分量，对步骤301中高通滤波器的分段时域积分模型进行补偿，得到的高通滤波器的分段时域电磁暂态修正模型如下：
[0078][0079]
其中，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，y(t-δt)为高通滤波器t-δt时刻的时域输出变量，x(t-δt)为高通滤波器t-δt时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，δ(t-δt)为t-δt时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数，b为常数。
[0080]
综上所述，本发明实施例中的高通滤波器的电磁暂态建模方法中通过构建高通滤波器的分段时域电磁暂态模型补偿分量，可以对高通滤波器的分段时域积分模型的模型精度进行修正，可以有效地提高修正后得到的高通滤波器的分段时域积分修正模型的实用性。
[0081]
与前述应用功能实现方法实施例相对应，本发明还提供了一种用于高通滤波器的电磁暂态建模的建模装置、电子设备及相应的实施例。
[0082]
图4为本发明实施例中提供的建模装置的一个实施例结构示意图。
[0083]
如图4所示，本发明实施例中的建模装置40包括：
[0084]
频域模型确定模块401、时域模型计算模块402和电磁暂态模型计算模块403；
[0085]
其中，频域模型确定模块401，用于根据高通滤波器确定高通滤波器的频域模型；
[0086]
时域模型计算模块402，用于根据拉普拉斯反变换和高通滤波器的频域模型，得到高通滤波器的分时段时域积分模型；
[0087]
电磁暂态模型计算模块403，用于根据梯形法则和高通滤波器的分时段时域积分模型，得到高通滤波器的分段时域电磁暂态模型。
[0088]
可选的，在本发明实施例的一种实施方式中，时域模型计算模块402具体用于执行
以下操作：利用拉普拉斯反变换对高通滤波器的频域模型进行求解得到高通滤波器的基础时域积分模型；根据冲激函数的特性，将高通滤波器的时域积分模型按照时间变量的取值范围划分为两部分，得到高通滤波器的分段时域积分模型。
[0089]
可选的，在本发明实施例的一种实施方式中，电磁暂态模型计算模块403具体用于执行以下操作：利用梯形法则对高通滤波器的分时段时域积分模型进行求解，得到高通滤波器的时域电磁暂态基础模型；根据冲激函数和求极限方法，将高通滤波器的时域电磁暂态基础模型按照时间变量的取值范围划分为三部分，得到高通滤波器的分段时域电磁暂态模型。
[0090]
可选的，在本发明实施例的一种实施方式中，频域模型确定模块401获取的高通滤波器的频域模型为：
[0091][0092]
其中，y(s)为高通滤波器的频域输出变量，x(s)为高通滤波器的频域输入变量；s为拉普拉斯算子，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数。
[0093]
进一步可选的，在本发明实施例的一种实施方式中，时域模型计算模块402计算得到的高通滤波器的基础时域积分模型为：
[0094][0095]
时域模型计算模块402进一步对上述基础时域积分模型进行求解得到高通滤波器的分时段时域积分模型为：
[0096][0097]
其中，t为时间变量，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数，b为常数。
[0098]
进一步可选的，在本发明实施例的一种实施方式中，电磁暂态模型计算模块403计算得到的高通滤波器的时域电磁暂态基础模型为：
[0099][0100]
电磁暂态模型计算模块403进一步对上述时域电磁暂态基础模型进行求解得到高通滤波器的分段时域电磁暂态模型为：
[0101][0102]
其中，t为时间变量，δt为时间步长，δt

为时间步长的右极限，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，y(t-δt)为高通滤波器t-δt时刻的时域输出变量，x(t-δt)为高通滤波器t-δt时刻的时域输入变量，x(0)为高通滤波器0时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，δ(t-δt)为t-δt时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数，b为常数。
[0103]
综上所述，本发明实施例中的建模装置及各个模块，可以通过拉普拉斯反变换和梯形法求解高通滤波器的分段时域电磁暂态模型，以实现对高通滤波器的电磁暂态建模，既能快速、便捷的建立高通滤波器的分段时域电磁暂态模型，又从数学角度提供了严谨的求解公式，数学逻辑的严谨性可以确保建立分段时域电磁暂态模型的准确性，有利于电力系统电磁暂态仿真及其实际应用
[0104]
图5为本发明实施例中提供的建模装置的另一个实施例结构示意图。
[0105]
如图5所示，本发明实施例中的建模装置40包括：频域模型确定模块401、时域模型计算模块402、电磁暂态模型计算模块403和电磁暂态模型补偿模块404；
[0106]
其中频域模型确定模块401、时域模型计算模块402、电磁暂态模型计算模块403执行的操作、功能与上述图4中的操作一样，此处不再赘述。
[0107]
电磁暂态模型补偿模块404用于执行以下操作：构建高通滤波器的分段时域电磁暂态模型补偿分量；利用高通滤波器的分段时域电磁暂态模型补偿分量，对高通滤波器的分段时域积分模型进行补偿，得到高通滤波器的分段时域电磁暂态修正模型。
[0108]
可选的，在本发明实施例的一种实施方式中，电磁暂态模型补偿模块404构建的高通滤波器的分段时域电磁暂态模型补偿分量为：
[0109][0110]
进一步，电磁暂态模型补偿模块404基于上述分段时域电磁暂态模型补偿分量对分段时域积分模型进行补偿后得到高通滤波器的分段时域电磁暂态修正模型为：
[0111][0112]
其中，y(t)为高通滤波器t时刻的时域输出变量，x(t)为高通滤波器t时刻的时域输入变量，y(t-δt)为高通滤波器t-δt时刻的时域输出变量，x(t-δt)为高通滤波器t-δt时刻的时域输入变量，δ(t)为t时刻的时域冲激函数值，δ(t-δt)为t-δt时刻的时域冲激函数值，t为高通滤波器时间常数，a为高通滤波器截止角频率系数，g为高通滤波器增益倍数，e为自然常数，b为常数。
[0113]
综上所述，本发明实施例中建模装置中的电磁暂态模型补偿模块404通过构建高通滤波器的分段时域电磁暂态模型补偿分量，可以对高通滤波器的分段时域积分模型的模型精度进行修正，可以有效地提高修正后得到的高通滤波器的分段时域积分修正模型的实用性。
[0114]
图6为本发明实施例中提供的电子设备的一个实施例结构示意图。
[0115]
如图6所示，本发明实施例中电子设备60包括存储器601和处理器602。存储器上存储有可执行代码，当可执行代码被处理器执行时，使处理器执行上述任一实施例中的方法。
[0116]
处理器602可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0117]
存储器601可以包括各种类型的存储单元，例如系统内存、只读存储器(rom)，和永久存储装置。其中，rom可以存储处理器602或者计算机的其他模块需要的静态数据或者指令。永久存储装置可以是可读写的存储装置。永久存储装置可以是即使计算机断电后也不会失去存储的指令和数据的非易失性存储设备。在一些实施方式中，永久性存储装置采用大容量存储装置(例如磁或光盘、闪存)作为永久存储装置。另外一些实施方式中，永久性存储装置可以是可移除的存储设备(例如软盘、光驱)。系统内存可以是可读写存储设备或者易失性可读写存储设备，例如动态随机访问内存。系统内存可以存储一些或者所有处理器在运行时需要的指令和数据。此外，存储器601可以包括任意计算机可读存储媒介的组合，包括各种类型的半导体存储芯片(dram，sram，sdram，闪存，可编程只读存储器)，磁盘和/或光盘也可以采用。在一些实施方式中，存储器601可以包括可读和/或写的可移除的存储设备，例如激光唱片(cd)、只读数字多功能光盘(例如dvd-rom，双层dvd-rom)、只读蓝光光盘、超密度光盘、闪存卡(例如sd卡、min sd卡、micro-sd卡等等)、磁性软盘等等。计算机可读存储媒介不包含载波和通过无线或有线传输的瞬间电子信号。
[0118]
存储器601上存储有可执行代码，当可执行代码被处理器602处理时，可以使处理
器602执行上文述及的方法中的部分或全部。
[0119]
此外，根据本发明的方法还可以实现为一种计算机程序或计算机程序产品，该计算机程序或计算机程序产品包括用于执行本发明的上述方法中部分或全部步骤的计算机程序代码指令。
[0120]
或者，本发明还可以实施为一种计算机可读存储介质(或机器可读存储介质)，其上存储有可执行代码(或计算机程序、或计算机指令代码)，当可执行代码(或计算机程序、或计算机指令代码)被电子设备(或电子设备、服务器等)的处理器执行时，使处理器执行根据本发明的上述方法的各个步骤的部分或全部。
[0121]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。