CN112668134B - 牵引网谐波阻抗确定方法、装置及终端设备 - Google Patents
牵引网谐波阻抗确定方法、装置及终端设备 Download PDFInfo
- Publication number
- CN112668134B CN112668134B CN201910974552.5A CN201910974552A CN112668134B CN 112668134 B CN112668134 B CN 112668134B CN 201910974552 A CN201910974552 A CN 201910974552A CN 112668134 B CN112668134 B CN 112668134B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- data
- group
- harmonic data
- impedance
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 17
- 230000015654 memory Effects 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000003908 quality control method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明提供了一种牵引网谐波阻抗确定方法及装置，该方法应用于电能质量控制技术领域，所述方法包括：获取多组谐波数据；基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数；根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数；若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合；根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值。本发明提供的牵引网谐波阻抗确定方法、装置及终端设备能够提高牵引网谐波阻抗计算的准确度。
Description
技术领域
本发明属于电能质量控制技术领域，更具体地说，是涉及一种牵引网谐波阻抗确定方法、装置及终端设备。
背景技术
电力列车向牵引供电系统注入谐波的问题一直受到铁路部门及研究人员的关注。牵引网谐波阻抗作为谐波治理和优化车网阻抗匹配的关键参数，其准确计算成为铁路部门及研究人员研究的热点。
目前牵引网谐波阻抗的计算方法主要有数值计算法、谐波电流注入法、主导波动量法以及线性回归法等，然而这几种方法均存在一定的缺陷：
数值计算法是根据牵引网结构和系统参数计算牵引网谐波阻抗的，其需要知道牵引网结构和系统参数才能完成计算，但全面准确地获取牵引网结构和系统参数难度较大，因此此方法在实际工程应用中有较大的局限性。
谐波电流注入法主要通过向牵引网系统注入谐波或间谐波电流，或开断某条支路来测量谐波阻抗，但此种方法可能会对牵引网系统造成影响。
主导波动量法基于牵引网受流点的谐波电压、谐波电流的主导波动量比值特征估计系统谐波阻抗，此种方法原理简单，但是要求较高的测量精度，并且主要用于较低次谐波的阻抗识别，无法对较高次谐波的阻抗进行有效计算。
线性回归法通过建立牵引网系统的等效电路模型，确定电路中公共连接点处谐波电流、谐波电压和谐波阻抗关系的方程，通过求解该方程来计算谐波阻抗，但因其缺少对变量奇异值和相关性问题分析，应用此方法计算得到的谐波阻抗的准确度也受到了影响。
因此，如何在保证谐波阻抗计算的谐波频次适用范围的同时提高谐波阻抗的计算准确性成为亟待解决的问题。
发明内容
本发明的目的在于提供一种牵引网谐波阻抗确定方法、装置及终端设备，以提高牵引网谐波阻抗的计算准确度。
本发明实施例的第一方面，提供了一种牵引网谐波阻抗确定方法，包括：
获取多组谐波数据；
基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数；
根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数；
若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合；
根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值。
本发明实施例的第二方面，提供了一种牵引网谐波阻抗确定装置，包括：
数据获取模块，用于获取多组谐波数据；
方程求解模块，用于基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数；
决定系数确定模块，用于根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数；
最优集确定模块，用于若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合；
谐波阻抗确定模块，用于根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值。
本发明实施例的第三方面，提供了一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的牵引网谐波阻抗确定方法的步骤。
本发明实施例的第四方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述的牵引网谐波阻抗确定方法的步骤。
本发明提供的牵引网谐波阻抗确定方法、装置及终端设备的有益效果在于：本发明实施例在现有技术的基础上，增加了根据谐波数据以及谐波数据对应的回归系数确定了样本决定系数的过程，并通过样本决定系数对回归系数进行进一步筛选，得到可信度最高的回归系数的集合，也即最优系数集合，最后利用最优系数集合中回归系数的值确定牵引网的谐波阻抗值。本发明实施例通过对回归系数进行进一步筛选，剔除了回归系数中可信度较低的部分，从而有效地提高了牵引网谐波阻抗的计算准确度。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一实施例提供的牵引网谐波阻抗确定方法的流程示意图；
图2的本发明另一实施例提供的牵引网谐波阻抗确定方法的流程示意图；
图3为本发明再一实施例提供的牵引网谐波阻抗确定方法的流程示意图；
图4为本发明又一实施例提供的牵引网谐波阻抗确定方法的流程示意图；
图5为本发明一实施例提供的牵引网谐波阻抗确定装置的结构框图；
图6为本发明一实施例提供的终端设备的示意框图。
具体实施方式
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参考图1，图1为本发明一实施例提供的牵引网谐波阻抗确定方法的流程示意图。该方法包括：
S101：获取多组谐波数据。
在本实施例中，获取多组谐波数据具体可包括：获取录波数据，并对录波数据进行快速傅里叶变化得到多组谐波数据，其中多组谐波数据包括多组谐波电压数据和多组谐波电流数据。
在本实施例中，预设回归方程可根据戴维南等效电路建立，其中，建立预设回归方程的步骤可以包括：
一、建立戴维南等效电路的过程，其中，戴维南等效电路可以为：
其中，
为牵引网系统两端的谐波电压相量值，Zsh为牵引网系统的等效谐波阻抗相量值，
为牵引网系统的等效电流相量值，
为牵引网系统内部电压源相量值。
二、根据戴维南等效电路建立回归方程的过程：
将
Zsh、
按照复数形式展开分别为：
将上述复数展开部分代入到戴维南等效电路中，并按照同一复数实部和虚部相等的原则，可以得到回归方程：
S102：基于每组谐波数据对回归方程进行求解，得到每组谐波数据对应的回归系数。
在本实施例中，谐波数据包括谐波电压数据和谐波电流数据，可采用偏最小二乘法，将每组谐波数据代入到回归方程中进行求解，得到回归方程的回归系数Zshx、Zshy。
S103：根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数。
在本实施例中，可根据每组谐波数据确定谐波阻抗的实测值，根据该组谐波数据对应的回归系数确定谐波阻抗的估计值，再根据谐波阻抗的实测值和估计值确定该组谐波数据对应的样本决定系数。
S104：若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合。
在本实施例中，可根据每组谐波数据对应的样本决定系数的值对该组谐波数据的回归系数进行筛选，将样本决定系数大于预设阈值的谐波数据对应的回归系数添加至最优系数集合。
S105：根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值。
在本实施例中，每一回归系数包括谐波阻抗的实部和虚部，根据每一回归系数都可得到对应的谐波阻抗，本实施例将所有谐波阻抗的平均值作为牵引网的谐波阻抗值。
从上述描述可知，本发明实施例在现有技术的基础上，增加了根据谐波数据以及谐波数据对应的回归系数确定了样本决定系数的过程，并通过样本决定系数对回归系数进行进一步筛选，得到可信度最高的回归系数的集合，也即最优系数集合，最后利用最优系数集合中回归系数的值确定牵引网的谐波阻抗值。本发明实施例通过对回归系数进行进一步筛选，剔除了回归系数中可信度较低的部分，从而有效地提高了牵引网谐波阻抗的计算准确度。
请一并参考图1及图2，图2为本申请另一实施例提供的牵引网谐波阻抗确定方法的流程示意图。在上述实施例的基础上，步骤S101可以详述为：
S201：获取牵引网受流点的录波数据。
S202：对录波数据进行快速傅里叶变换，得到牵引网的多组谐波数据。
在本实施例中，可首先对获取到的录波数据每第一预设数量(比如10)个周波进行一次快速傅里叶变换，得到谐波数据样本集。再以第二预设数量(比如100)个样本点作为滑动窗口的宽度，每滑动一个样本点即确定一组谐波数据，从而得到牵引网的多组谐波数据。例如，一共104个样本点，则第1～100个样本点为一组谐波数据，第2～101个样本点为一组谐波数据，第3～102个样本点为一组谐波数据，第4～103个样本点为一组谐波数据，第5～104个样本点为一组谐波数据。
请一并参考图1及图3，图3为本申请再一实施例提供的牵引网谐波阻抗确定方法的流程示意图。在上述实施例的基础上，在基于每组谐波数据对回归方程进行求解之前，还可以包括：
S301：确定每组谐波数据中各个谐波数据的统计量。
在本实施例中，对于某一谐波数据x，其统计量Gn为：
其中，
为该谐波数据所在组中所有谐波数据的均值，s为谐波数据所在组中所有谐波数据的标准差。
S302：若某一谐波数据的统计量大于预设检出值，则清除该谐波数据。
在本实施例中，可首先确定检出水平α，根据检出水平α查表得到预设检出值
若对于某一组谐波数据中的某一谐波数据x，有
则确定该谐波数据x为离群值，并从该组中清除该谐波数据。也即，若某组谐波数据中某一谐波数据的统计量大于预设检出值，则从该组中清除该谐波数据。
请一并参考图1及图4，图4为本申请又一实施例提供的牵引网谐波阻抗确定方法的流程示意图。在上述实施例的基础上，步骤S103可以详述为：
S401：根据每组谐波数据确定该组谐波数据对应的阻抗实测值。
在本实施例中，一个分组内包含有n个谐波数据，每个谐波数据都包含一个谐波电压数据和一个谐波电流数据，因此可根据每个谐波数据得到一个谐波阻抗值，该谐波阻抗值即该谐波数据对应的阻抗实测值，因此该组谐波数据可得到n个阻抗实测值，此n个阻抗实测值即为该组谐波数据对应的阻抗实测值。
S402：根据该组谐波数据对应的回归系数确定该组谐波数据对应的阻抗估计值。
在本实施例中，回归系数包含有谐波阻抗的实部和虚部，在上述实施例的基础上，也即Zshx和Zshy。在得到该组谐波数据对应的回归系数后，可根据谐波阻抗的实部和虚部得到该组谐波数据对应的谐波阻抗，此谐波阻抗即为该组谐波数据对应的谐波阻抗估计值。
S403：根据该组谐波数据对应的阻抗实测值和阻抗估计值确定该组谐波数据对应的样本决定系数。
在本实施例中，样本决定系数r2的确定方法可以为：
其中，SSR为该组谐波数据的回归平方和，SST为该组谐波数据的总离差平方和，
为该组谐波数据对应的谐波估计值，n为该组谐波数据中谐波数据的个数，yi为该组谐波数据对应的谐波实测值，
为该组谐波数据中n个谐波数据对应的阻抗实测值的平均值。
在本实施例中，SSR能够反映因变量也就是阻抗估计值的不确定性，而SST是由自变量也就是实测值决定的，因此样本决定系数r2能够反映因变量的波动中能用自变量解释的比例，因此可使用样本决定系数对可信度较高的回归系数进行筛选，从而提高谐波阻抗计算的准确性。
可选地，作为本发明实施例提供的牵引网谐波阻抗确定方法的一种具体实施方式，根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值，可以包括：
根据最优系数集合中每组回归系数的值确定每组回归系数对应的谐波阻抗值，并将所有谐波阻抗值的平均值确定为牵引网的谐波阻抗值。
可选地，作为本发明实施例提供的牵引网谐波阻抗确定方法的一种具体实施方式，在上述实施例的基础上，牵引网谐波阻抗确定方法还可以包括：
若最优系数集合中的数据量未达到预设数据量或最优系数集合中每组回归系数对应的谐波数据的频次未覆盖预设频次，则重新获取多组新的谐波数据，并基于重新获取的多组新的谐波数据对最优系数集合进行更新。
在本实施例中，重新获取多组新的谐波数据，并基于重新获取的多组新的谐波数据对最优系数集合进行更新可以包括：
重新获取多组新的谐波数据，并计算每组新的谐波数据对应的样本决定系数。
基于每组新的谐波数据对应的样本决定系数对最优系数集合进行更新。
其中，最优系数集合的更新方法可以为：
若某组谐波数据对应的样本决定系数不小于最优系数集合对应的样本决定系数的平均值，则将该组谐波数据对应的回归系数添加至最优系数集合。
若某组谐波数据对应的样本决定系数小于最优系数集合对应的样本决定系数的平均值，则舍弃该组谐波数据对应的回归系数。
其中，最优系数集合包含的是回归系数，每个回归系数对应一组谐波数据，每组谐波数据对应一个样本决定系数，也即每个回归系数都对应一个样本决定系数，那么最优系数集合对应的样本决定系数的平均值也即最优系数集合中所有回归系数对应的样本决定系数的平均值。
可选地，在本发明的一个实施例中，还可根据样本决定系数的值对回归系数进行分类，得到第一类集合、第二类集合以及第三类集合。
若第一类集合中的数据量未达到预设数据量或第一类集合中每组回归系数对应的谐波数据的频次未覆盖预设频次时，则重新获取新的多组谐波数据，并基于重新获取的多组新的谐波数据对第一类集合、第二类集合以及第三类集合进行更新，最后根据第一类集合中回归系数的值确定牵引网的谐波阻抗值。
其中，根据样本决定系数的值对回归系数进行分类可以包括：
若r2∈[0.9,1]，则将该样本决定系数对应的回归系数添加至第一类集合中，此集合即相当于上述实施例中的最优系数集合。
若r2∈[0.7,0.9)，则将该样本决定系数对应的回归系数添加至第二类集合中，也即可信度次高的系数集合。
若r2∈[0.3,0.7)，则将该样本决定系数对应的回归系数添加至第三类集合中，也即中等可信度的系数集合。
若r2∈(0,0.3]，则将该样本决定系数对应的回归系数作为不可信系数舍弃。
其中，对各个类别的回归系数的更新方法可以为：
若某组谐波数据对应的样本决定系数不小于第一类集合对应的样本决定系数的平均值，则将该组谐波数据对应的回归系数添加至第一类集合。
若某组谐波数据对应的样本决定系数小于第一类集合对应的样本决定系数的平均值，则检测该样本决定系数是否小于第二类集合中对应的样本决定系数的平均值。
若某组谐波数据对应的样本决定系数不小于第二类集合对应的样本决定系数的平均值，则将该组谐波数据对应的回归系数添加至第二类集合。
若某组谐波数据对应的样本决定系数小于第二类集合对应的样本决定系数的平均值，则检测该样本决定系数是否小于第三类集合中对应的样本决定系数的平均值。
若某组谐波数据对应的样本决定系数不小于第三类集合对应的样本决定系数的平均值，则将该组谐波数据对应的回归系数添加至第三类集合。
若某组谐波数据对应的样本决定系数小于第三类集合对应的样本决定系数的平均值，则舍弃该组谐波数据对应的回归系数。
在完成各个分类集合的更新后：
若更新后的第一类集合达到预设数据量，但每组回归系数对应的谐波数据的频次仍未覆盖预设频次，则将第二类集合中的部分回归系数添加至第一类集合中。其中，所添加的第二类集合中的部分回归系数对应的谐波数据的频次未在第一类集合中的回归系数对应的谐波数据的频次范围内。
同理，若需更新第二类集合，也可将第三类集合中的部分回归系数添加至第二类集合中，此处不再赘述。
对应于上文实施例的牵引网谐波阻抗确定方法，图5为本发明一实施例提供的牵引网谐波阻抗确定装置的结构框图。为了便于说明，仅示出了与本发明实施例相关的部分。参考图5，该装置包括：数据获取模块510、方程求解模块520、决定系数确定模块530、最优集确定模块540和谐波阻抗确定模块550。
其中，数据获取模块510，用于获取多组谐波数据。
方程求解模块520，用于基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数。
决定系数确定模块530，用于根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数。
最优集确定模块540，用于若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合。
谐波阻抗确定模块550，用于根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值。
参考图5，在本发明的另一个实施例中，数据获取模块510可以包括：
数据获取单元511，用于获取牵引网受流点的录波数据。
数据转换单元512，用于对录波数据进行快速傅里叶变换，得到牵引网的多组谐波数据。
参考图5，在本发明的再一个实施例中，牵引网谐波阻抗确定装置还可以包括数据清洗模块560，数据清洗模块560可以包括：
统计量计算单元561，用于确定每组谐波数据中各个谐波数据的统计量。
数据清洗单元562，用于若某一谐波数据的统计量大于预设检出值，则清除该谐波数据。
参考图5，在本发明的又一个实施例中，决定系数确定模块530可以包括：
实测值计算单元531，用于根据每组谐波数据确定该组谐波数据对应的阻抗实测值。
估计值计算单元532，用于根据该组谐波数据对应的回归系数确定该组谐波数据对应的阻抗估计值。
决定系数确定单元533，用于根据该组谐波数据对应的阻抗实测值和阻抗估计值确定该组谐波数据对应的样本决定系数。
可选地，作为本发明实施例提供的牵引网谐波阻抗确定装置的一种具体实施方式，根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值，可以包括：
根据最优系数集合中每组回归系数的值确定每组回归系数对应的谐波阻抗值，并将所有谐波阻抗值的平均值确定为牵引网的谐波阻抗值。
可选地，作为本发明实施例提供的牵引网谐波阻抗确定装置的一种具体实施方式，若最优系数集合中的数据量未达到预设数据量或最优系数集合中每组回归系数对应的谐波数据的频次未覆盖预设频次，则重新获取多组新的谐波数据，并基于重新获取的多组新的谐波数据对最优系数集合进行更新。
参见图6，图6为本发明一实施例提供的一种终端设备的示意框图。如图6所示的本实施例中的终端600可以包括：一个或多个处理器601、一个或多个输入设备602、一个或多个输出设备603及一个或多个存储器604。上述处理器601、输入设备602、则输出设备603及存储器604通过通信总线605完成相互间的通信。存储器604用于存储计算机程序，计算机程序包括程序指令。处理器601用于执行存储器604存储的程序指令。其中，处理器601被配置用于调用程序指令执行以下操作上述各装置实施例中各模块/单元的功能，例如图5所示模块510至560的功能。
应当理解，在本发明实施例中，所称处理器601可以是中央处理单元(CentralProcessing Unit，CPU)，该处理器还可以是其他通用处理器、数字信号处理器(DigitalSignal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
输入设备602可以包括触控板、指纹采传感器(用于采集用户的指纹信息和指纹的方向信息)、麦克风等，输出设备603可以包括显示器(LCD等)、扬声器等。
该存储器604可以包括只读存储器和随机存取存储器，并向处理器601提供指令和数据。存储器604的一部分还可以包括非易失性随机存取存储器。例如，存储器604还可以存储设备类型的信息。
具体实现中，本发明实施例中所描述的处理器601、输入设备602、输出设备603可执行本发明实施例提供的牵引网谐波阻抗确定方法的第一实施例和第二实施例中所描述的实现方式，也可执行本发明实施例所描述的终端的实现方式，在此不再赘述。
在本发明的另一实施例中提供一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序包括程序指令，程序指令被处理器执行时实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，计算机程序包括计算机程序代码，计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括：能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。
计算机可读存储介质可以是前述任一实施例的终端的内部存储单元，例如终端的硬盘或内存。计算机可读存储介质也可以是终端的外部存储设备，例如终端上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。进一步地，计算机可读存储介质还可以既包括终端的内部存储单元也包括外部存储设备。计算机可读存储介质用于存储计算机程序及终端所需的其他程序和数据。计算机可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的终端和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的终端和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (9)
1.一种牵引网谐波阻抗确定方法，其特征在于，包括：
获取多组谐波数据；
基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数；
根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数；
若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合；
根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值；
其中，所述根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数，包括：
根据每组谐波数据确定该组谐波数据对应的阻抗实测值；
根据该组谐波数据对应的回归系数确定该组谐波数据对应的阻抗估计值；
根据该组谐波数据对应的阻抗实测值和阻抗估计值确定该组谐波数据对应的样本决定系数；
其中，样本决定系数r2的确定方法为：
其中，SSR为该组谐波数据的回归平方和，SST为该组谐波数据的总离差平方和，
为该组谐波数据对应的阻抗估计值，n为该组谐波数据中谐波数据的个数，yi为该组谐波数据对应的阻抗实测值，
为该组谐波数据中n个谐波数据对应的阻抗实测值的平均值。
2.如权利要求1所述的牵引网谐波阻抗确定方法，其特征在于，所述获取多组谐波数据，包括：
获取牵引网受流点的录波数据；
对录波数据进行快速傅里叶变换，得到牵引网的多组谐波数据。
3.如权利要求1所述的牵引网谐波阻抗确定方法，其特征在于，在所述基于每组谐波数据对回归方程进行求解之前，还包括：
确定每组谐波数据中各个谐波数据的统计量；
若某一谐波数据的统计量大于预设检出值，则清除该谐波数据。
4.如权利要求1所述的牵引网谐波阻抗确定方法，其特征在于，所述根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值，包括：
根据最优系数集合中每组回归系数的值确定每组回归系数对应的谐波阻抗值，并将所有谐波阻抗值的平均值确定为牵引网的谐波阻抗值。
5.如权利要求1-4任一项所述的牵引网谐波阻抗确定方法，其特征在于，在所述根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值之前，还包括：
若最优系数集合中的数据量未达到预设数据量或最优系数集合中每组回归系数对应的谐波数据的频次未覆盖预设频次，则重新获取多组新的谐波数据，并基于重新获取的多组新的谐波数据对最优系数集合进行更新。
6.一种牵引网谐波阻抗确定装置，其特征在于，包括：
数据获取模块，用于获取多组谐波数据；
方程求解模块，用于基于每组谐波数据对预设回归方程进行求解，得到每组谐波数据对应的回归系数；
决定系数确定模块，用于根据每组谐波数据以及该组谐波数据对应的回归系数确定每组谐波数据对应的样本决定系数；
最优集确定模块，用于若某一组谐波数据对应的样本决定系数大于预设阈值，则将该组谐波数据对应的回归系数添加至最优系数集合；
谐波阻抗确定模块，用于根据最优系数集合中回归系数的值确定牵引网的谐波阻抗值；
其中，所述决定系数确定模块包括：
实测值计算单元，用于根据每组谐波数据确定该组谐波数据对应的阻抗实测值；
估计值计算单元，用于根据该组谐波数据对应的回归系数确定该组谐波数据对应的阻抗估计值；
决定系数确定单元，用于根据该组谐波数据对应的阻抗实测值和阻抗估计值确定该组谐波数据对应的样本决定系数；
其中，样本决定系数r2的确定方法为：
其中，SSR为该组谐波数据的回归平方和，SST为该组谐波数据的总离差平方和，
为该组谐波数据对应的阻抗估计值，n为该组谐波数据中谐波数据的个数，yi为该组谐波数据对应的阻抗实测值，
为该组谐波数据中n个谐波数据对应的阻抗实测值的平均值。
7.如权利要求6所述的牵引网谐波阻抗确定装置，其特征在于，所述数据获取模块包括：
数据获取单元，用于获取牵引网受流点的录波数据；
数据转换单元，用于对录波数据进行快速傅里叶变换，得到牵引网的多组谐波数据。
8.一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述方法的步骤。
9.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至5任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910974552.5A CN112668134B (zh) | 2019-10-14 | 2019-10-14 | 牵引网谐波阻抗确定方法、装置及终端设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910974552.5A CN112668134B (zh) | 2019-10-14 | 2019-10-14 | 牵引网谐波阻抗确定方法、装置及终端设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112668134A CN112668134A (zh) | 2021-04-16 |
| CN112668134B true CN112668134B (zh) | 2022-11-01 |
Family
ID=75399714
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910974552.5A Active CN112668134B (zh) | 2019-10-14 | 2019-10-14 | 牵引网谐波阻抗确定方法、装置及终端设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112668134B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115290984B (zh) * | 2022-08-16 | 2024-08-13 | 西南交通大学 | 一种地铁钢轨电位限制装置等效电阻测算方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107123990A (zh) * | 2017-03-29 | 2017-09-01 | 广西电网有限责任公司电力科学研究院 | 城市轨道交通负荷牵引电缆谐波放大的治理方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102749521B (zh) * | 2012-07-18 | 2015-04-08 | 华北电力大学（保定） | 一种电力系统谐波阻抗计算方法 |
| CN105588982A (zh) * | 2014-10-27 | 2016-05-18 | 国家电网公司 | 一种电力系统谐波阻抗的测量计算方法 |
| CN105548718A (zh) * | 2015-12-14 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | 一种基于混合整体最小二乘法的系统谐波阻抗计算方法 |
| CN107576853B (zh) * | 2017-07-12 | 2019-11-05 | 国网上海市电力公司 | 基于典型相关性分析的配电网谐波阻抗计算方法 |
| CN110059295B (zh) * | 2019-04-08 | 2021-04-13 | 国网江苏省电力有限公司电力科学研究院 | 一种基于三点筛选法的系统谐波阻抗估计方法 |
-
2019
- 2019-10-14 CN CN201910974552.5A patent/CN112668134B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107123990A (zh) * | 2017-03-29 | 2017-09-01 | 广西电网有限责任公司电力科学研究院 | 城市轨道交通负荷牵引电缆谐波放大的治理方法 |
Non-Patent Citations (1)
| Title |
|---|
| Harmonic Resonance Assessment to Traction Power-Supply System Considering Train Model in China High-Speed Railway;Zhengyou He;《IEEE TRANSACTIONS ON POWER DELIVERY》;20140831;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112668134A (zh) | 2021-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109034244B (zh) | 基于电量曲线特征模型的线损异常诊断方法及装置 | |
| CN101661059B (zh) | 用户端非线性负载的谐波电压发射水平的定量分析方法及装置 | |
| CN109472231B (zh) | 变电站谐波源定位方法及装置、计算机可读存储介质 | |
| CN111008510B (zh) | 一种台区内电表识别方法及判断电表所属台区的方法 | |
| CN108802535B (zh) | 筛选方法、主干扰源识别方法及装置、服务器及存储介质 | |
| CN109064058B (zh) | 一种避雷器绝缘性能分析方法及装置 | |
| CN112598537B (zh) | 电力设备故障诊断方法、装置及终端设备 | |
| CN108900622A (zh) | 基于物联网的数据融合方法、装置与计算机可读存储介质 | |
| CN107589391A (zh) | 一种检测电能计量装置整体误差的方法、装置和系统 | |
| CN112668134B (zh) | 牵引网谐波阻抗确定方法、装置及终端设备 | |
| CN110348717B (zh) | 基于栅格粒度的基站价值评分方法和装置 | |
| Zhu et al. | Injection amplitude guidance for impedance measurement in power systems | |
| US20170160321A1 (en) | Measurement apparatus | |
| CN112730654B (zh) | 六氟化硫电气设备故障检测方法、装置及终端设备 | |
| CN114167132A (zh) | 无线终端的功耗检测方法、装置、电子设备及存储介质 | |
| CN107132500B (zh) | 一种同步相量测量单元在线校准方法与装置 | |
| CN111736648A (zh) | 输出电压控制方法、装置及终端设备 | |
| CN111210140A (zh) | 税务信用评价方法、设备及可读存储介质 | |
| CN102937676B (zh) | 负荷谐波注入污染预警的实现方法及其系统 | |
| CN113255137B (zh) | 目标对象应变数据的处理方法、装置及存储介质 | |
| CN115047294A (zh) | 配电线路的故障检测方法、装置、设备及存储介质 | |
| CN115047400A (zh) | 三相电能表准确度校验方法、系统、终端设备及介质 | |
| CN116136580A (zh) | 一种基于感知设备的线损分级检测方法及装置 | |
| CN109375055B (zh) | 一种相量测量系统的检测方法、系统、装置及储存介质 | |
| CN109638845B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |