CN114859120B - 谐波分析方法 - Google Patents
谐波分析方法 Download PDFInfo
- Publication number
- CN114859120B CN114859120B CN202210683887.3A CN202210683887A CN114859120B CN 114859120 B CN114859120 B CN 114859120B CN 202210683887 A CN202210683887 A CN 202210683887A CN 114859120 B CN114859120 B CN 114859120B
- Authority
- CN
- China
- Prior art keywords
- sampling
- vector
- harmonic
- zero
- 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.)
- Active
Links
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
- 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
- G01R23/165—Spectrum analysis; Fourier analysis using filters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Abstract
本发明公开了一种谐波分析方法，以固定采样频率进行采样，还包括以下步骤：步骤1、根据采样信号计算信号的工作频率估计值；步骤2、根据工作频率的估计值动态调整采样窗口长度；步骤3、根据工作频率估计值、采样频率、采样窗口长度和谐波分析范围计算初始角度向量和增量角度向量；步骤4、缓存长度一段采样信号作为计算输入信号，计算输入信号的与初始角度向量的长度相等；根据计算输入信号、初始角度向量和增量角度向量计算出谐波分析范围内的谐波含量。本发明突破了常规FFT方法数据点数的限制，解决了常规方法在频偏情况下结果不准确的问题，同时还具有硬件成本低、不影响电能计量、运算量低、计算速度快、易于工程实现等优点。
Description
技术领域
本发明涉及一种对电网中的谐波进行分析的方法。
背景技术
电网的谐波信息是衡量电网供电质量的重要数据。常规的谐波分析方法是以固定采样间隔，采集若干个整数周波的数据后进行快速傅里叶变换(FFT)，根据变换后的各条谱线信息计算谐波含量。常规的FFT方法具有以下缺陷：
1、每次进行FFT的数据点数必须为2的整数次幂，比如512、1024等。IEC61000-4-30中推荐累计10周波数据进行一次谐波分析，如果累计点数为1024，那么每周波采样点数为102.4点，这对采样器件提出了较高的要求。
2、在频率偏移(频偏)的情况下，比如电网供电质量较差、实际频率为49Hz，如果以固定采样间隔，采集固定点数进行FFT，由于存在栅栏效应和频谱泄露，即使采用加窗和插值拟合的方法，也无法得到很高的精度。
针对常规FFT方法存在的缺陷，现有技术提出了另一种谐波分析的方法，它是基于硬件频率跟踪，动态调整采样频率，即无论电网频率如何变化，通过调整采样间隔，每周波都采样固定的点数，如512点。这种方法克服了上述常规FFT方法存在的两个缺陷，但是也有其局限性：
1、硬件频率跟踪常采用FPGA等可编程器件实现，硬件成本较高。
2、带有谐波分析功能的设备一般都会有电能计量功能，而电能的累计是通过瞬时功率的积分计算得到的。如果修改了采样间隔，那么在频率变化的这段时间内，积分间隔是错误的，会导致计量错误，这对于计量设备来说是不可接受的。
发明内容
本发明提出了一种谐波分析方法，其目的是：(1)突破常规FFT方法数据点数必须为2的整数次幂的限制；(2)解决常规方法在频偏情况下，计算得到的谐波数据不准确的问题；(3)克服使用硬件频率跟踪元件所带来的硬件成本过高的问题；(4)避免采样间隔变化导致的电能计量错误；(5)降低运算量，提高计算效率，易于工程实现。
本发明技术方案如下：
一种谐波分析方法，以固定采样频率进行采样，还包括以下步骤：
步骤1、根据采样信号计算信号的工作频率估计值；
步骤2、根据工作频率的估计值动态调整采样窗口长度；
步骤3、根据工作频率估计值、采样频率、采样窗口长度和谐波分析范围计算初始角度向量和增量角度向量；
步骤4、缓存长度一段采样信号作为计算输入信号，计算输入信号的与初始角度向量的长度相等；根据计算输入信号、初始角度向量和增量角度向量计算出谐波分析范围内的谐波含量；
执行步骤2至4的过程中，同步执行步骤1，当工作频率估计值发生变化时，中止步骤2至4的执行过程，跳转至步骤2重新计算。
作为所述的谐波分析方法的进一步改进：步骤1中，采用过零法计算工作频率估计值：
步骤1-1、将采样信号通过低通滤波器，获得仅包含基波的信号；
步骤1-2、初始化过零点计数Cnt-zero为0，经历时间值T-zero为0；
步骤1-3、遍历各采样点，直至过零点计数Cnt-zero达到要求点数：若当前采样点的值大于0并且上一个采样点的值小于0，则当前点视为过零点，Cnt-zero加1，并且累加从上一个过零点到当前过零点所经历的时间至T-zero；
步骤1-4、将Cnt-zero除以T-zero，将得到的值作为工作频率估计值fm。
作为所述的谐波分析方法的进一步改进：步骤2中，设工作频率估计值为fm，采样频率为fs，则每周波的采样点数Nm为：
式中运算符号round(·，n)表示对·四舍五入并保留n位小数；
然后，按以下原则选择采样窗口覆盖的周波数Nw：
式中min(Nw)表示选择满足限定条件的Nw的最小值，表示对·向下取整操作，Threshold表示预设的门限值；
采样窗口长度即采样窗口中包含的采样点的数量NL＝round(Nm×Nw，0)。
作为所述的谐波分析方法的进一步改进：Threshold＝0.3。
作为所述的谐波分析方法的进一步改进：步骤3中，设工作频率估计值为fm，采样频率为fs，采样窗口长度为NL，谐波分析范围中的谐波次数以等差数列的形式表示为：
[Indexstart：Indexstop：interval]
式中Indexstart表示首个谐波次数，Indexstop表示最后一个谐波次数，interval表示谐波次数的间隔值；
则初始角度向量中的各元素分别为：
增量角度向量中的各元素分别为：
上式中Amp为角度扩大倍数；
对于初始角度向量M和增量角度向量P中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp。
作为所述的谐波分析方法的进一步改进：在执行步骤4之前，先将计算所需的三角函数值存储到查询表tcos中，具体步骤为：
设角度缩小倍数Amp′，Amp′＜Amp，以Amp′/Amp为分辨率取0度至450度的所有角度值；对于各角度α，以α*Amp/Amp′作为查询表的序号值，以cos(α)为该序号值n＝α*Amp/Amp′对应的三角函数值，即tcos[n]＝cos(α)。
作为所述的谐波分析方法的进一步改进：步骤4的具体步骤为：
步骤4-1、分配用于存放向量的长度为NL的空间，对向量/>的各元素进行初始化：
Di＝Mi i＝0，1，2，...，NL-1
步骤4-2、分配两个长度均为NL的空间，分别用于存放实部转换向量和虚部转换向量/>对实部转换向量/>和虚部转换向量/>进行初始化：
TRi＝tcos[Di/Amp′]i＝0，1，2，...，NL-1
步骤4-3、设k＝Indexstart，以Rek、Imk分别表示第k次谐波含量的实部和虚部，则第k次谐波含量的实部和虚部的计算方式如下：
式中为用于存放计算输入信号的输入信号向量，该向量长度为NL；
步骤4-4、将向量与向量/>相加并保存至向量/>
谐波次数k递增：k＝k+inteFval；
如果k＞Indexstop，则退出，否则继续执行步骤4-5；
步骤4-5、遍历当前向量中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp；
步骤4-6、根据当前的刷新实部转换向量/>和虚部转换向量/>
TRi＝tcos[Di/Amp′]i＝0，1，2，...，NL-1
步骤4-7、计算第k次谐波含量的实部和虚部：
步骤4-8、跳转至步骤4-4。
相对于现有技术，本发明具有以下有益效果：(1)采用本发明方法，数据点数可为任意整数，突破了常规FFT方法中必须为2的整数次幂的限制；(2)本方法根据工作频率估计值动态调整采样窗口长度，避免了栅栏效应，并且通过合理设置Nw，将频偏时频谱泄露造成的误差控制在可接受的范围内，保证了谐波数据的准确；(3)本方法可以指定某个感兴趣的谐波范围进行分析，其余谐波含量不计算，降低了运算量，提高了运算速度；(4)根据运算数据的周期特性，通过将查表法、倍数扩大和向量增量计算等手段相结合，将计算过程中的浮点运算和三角函数运算转换为整数加法运算和查表运算，既易于工程实现，也大大降低了运算量，并且减少了RAM开销；(5)使用固定采样间隔，无需增加硬件，降低了硬件成本，并且不会因为采样间隔变化导致计量错误。
附图说明
图1为本方法的流程示意图。
具体实施方式
下面结合附图详细说明本发明的技术方案：
如图1，一种谐波分析方法，以固定采样频率fs对电网信号进行采样，还包括以下步骤：
步骤1、根据采样信号计算信号的工作频率估计值。
本实施例中，采用过零法计算工作频率估计值：
步骤1-1、将采样信号通过低通滤波器，获得仅包含基波的信号；
步骤1-2、初始化过零点计数Cnt-zero为0，经历时间值T-zero为0；
步骤1-3、遍历各采样点，直至过零点计数Cnt-zero达到要求点数：若当前采样点的值大于0并且上一个采样点的值小于0，则当前点视为过零点，Cnt-zero加1，并且累加从上一个过零点到当前过零点所经历的时间至T-zero；
步骤1-4、将Cnt-zero除以T-zero，将得到的值作为工作频率估计值fm。
步骤2、根据工作频率的估计值动态调整采样窗口长度。
具体的，每周波的采样点数Nm为：
式中运算符号round(·，n)表示对·四舍五入并保留n位小数；
然后，按以下原则选择采样窗口覆盖的周波数Nw：
式中min(Nw)表示选择满足限定条件的Nw的最小值，表示对·向下取整操作，Threshold表示预设的门限值。优选的，Threshold＝0.3。Nw一般取值为10或者20。
采样窗口长度即采样窗口中包含的采样点的数量NL＝round(Nm×Nw，0)。
步骤3、根据工作频率估计值、采样频率、采样窗口长度和谐波分析范围计算初始角度向量和增量角度向量。
具体的，谐波分析范围中的谐波次数以等差数列的形式表示为：
[Indexstart：Indexstop：interval]
式中Indexstart表示首个谐波次数，Indexstop表示最后一个谐波次数，interval表示谐波次数的间隔值。
举例说明：
[1∶50∶1]表示从1次谐波，即基波开始，依次分析2、3、4直至50次谐波。
则初始角度向量中的各元素分别为：
增量角度向量中的各元素分别为：
上式中Amp为角度扩大倍数，本实施例中，Amp＝10000。
对于初始角度向量M和增量角度向量P中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp。
在执行步骤4之前，先将计算所需的三角函数值存储到查询表tcos中，从而在步骤4中可以直接通过查表的方式获得三角函数值，降低实时计算量。设置查询表的具体步骤为：
设角度缩小倍数Amp′，Amp′＜Amp，以Amp′/Amp为分辨率取0度至450度的所有角度值。本实施例中，Amp′＝100，因此角度的分辨率为0.01度。对于各角度α，以α*100作为查询表的序号值(为整数值)，以cos(α)为该序号值n＝α*100对应的三角函数值，即tcos[n]＝cos(α)。
说明：此处余弦查询表的建立，相对于前面的步骤1至3，没有严格的先后顺序要求。只要在步骤4开始之前完成查询表的建立即可。
步骤4、缓存长度一段采样信号作为计算输入信号，计算输入信号的与初始角度向量的长度相等；根据计算输入信号、初始角度向量和增量角度向量计算出谐波分析范围内的谐波含量。
具体步骤为：
步骤4-1、分配用于存放向量的长度为NL的空间，对向量/>的各元素进行初始化：
Di＝Mi i＝0，1，2，...，NL-1
步骤4-2、分配两个长度均为NL的空间，分别用于存放实部转换向量和虚部转换向量/>对实部转换向量/>和虚部转换向量/>进行初始化：
TRi＝tcos[Di/100]i＝0，1，2，...，NL-1
步骤4-3、设k＝Indexstart，以Rek、Imk分别表示第k次谐波含量的实部和虚部，则第k次谐波含量的实部和虚部的计算方式如下：
式中为用于存放计算输入信号的输入信号向量，该向量长度为NL；
步骤4-4、将向量与向量/>相加并保存至向量/>
谐波次数k递增：k＝k+inteFVal；
如果k＞Indexstop，则退出，否则继续执行步骤4-5；
步骤4-5、遍历当前向量中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp；
步骤4-6、根据当前的刷新实部转换向量/>和虚部转换向量/>
TRi＝tcos[Di/100]i＝0，1，2，...，NL-1
步骤4-7、计算第k次谐波含量的实部和虚部：
步骤4-8、跳转至步骤4-4。
执行步骤2至4的过程中，同步执行步骤1，当工作频率估计值发生变化时，中止步骤2至4的执行过程，跳转至步骤2重新计算。
本发明的步骤4中的谐波运算只需要两次整数加法、一次整数减法和一次整数除法，以及两次查表操作，大大提升了算法的计算速度.并且从空间占用上看，步骤4-1和4-2中的向量空间可以重复使用，无需重复开辟，节省了大量的内存空间。
实际工程应用中，在STM32H743单片机上使用本发明中的方法进行1280点的1至50次谐波分析耗时约2.4ms，并且本发明计算得到的谐波分析数值相对于期望值，误差在0.05％以内，突破了传统FFT对数据点数必须为2的整数次幂的限制，解决了频偏时误差过大的缺陷，其运算精度和运算速度已具备在嵌入式系统中应用的可行性。即使后续需要增加谐波的分析范围，比如分析1-60次谐波，其运算时间的增加也是线性的，约为1-50次运算时间的1.2倍。
Claims (6)
1.一种谐波分析方法，以固定采样频率进行采样，其特征在于还包括以下步骤：
步骤1、根据采样信号计算信号的工作频率估计值；
步骤2、根据工作频率的估计值动态调整采样窗口长度；
步骤3、根据工作频率估计值、采样频率、采样窗口长度和谐波分析范围计算初始角度向量和增量角度向量；
步骤3中，设工作频率估计值为fm，采样频率为fs，采样窗口长度为NL，谐波分析范围中的谐波次数以等差数列的形式表示为：
[Indexstart:Indexstop:interval]
式中Indexstart表示首个谐波次数，Indexstop表示最后一个谐波次数，interval表示谐波次数的间隔值；
则初始角度向量中的各元素分别为：
增量角度向量中的各元素分别为：
上式中Amp为角度扩大倍数；
对于初始角度向量M和增量角度向量P中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp；
步骤4、缓存长度一段采样信号作为计算输入信号，计算输入信号的与初始角度向量的长度相等；根据计算输入信号、初始角度向量和增量角度向量计算出谐波分析范围内的谐波含量；
执行步骤2至4的过程中，同步执行步骤1，当工作频率估计值发生变化时，中止步骤2至4的执行过程，跳转至步骤2重新计算。
2.如权利要求1所述的谐波分析方法，其特征在于：步骤1中，采用过零法计算工作频率估计值：
步骤1-1、将采样信号通过低通滤波器，获得仅包含基波的信号；
步骤1-2、初始化过零点计数Cnt-zero为0，经历时间值T-zero为0；
步骤1-3、遍历各采样点，直至过零点计数Cnt-zero达到要求点数：若当前采样点的值大于0并且上一个采样点的值小于0，则当前点视为过零点，Cnt-zero加1，并且累加从上一个过零点到当前过零点所经历的时间至T-zero；
步骤1-4、将Cnt-zero除以T-zero，将得到的值作为工作频率估计值fm。
3.如权利要求1所述的谐波分析方法，其特征在于：步骤2中，设工作频率估计值为fm，采样频率为fs，则每周波的采样点数Nm为：
式中运算符号round(·,n)表示对·四舍五入并保留n位小数；
然后，按以下原则选择采样窗口覆盖的周波数Nw：
式中min(Nw)表示选择满足限定条件的Nw的最小值，表示对·向下取整操作，Threshold表示预设的门限值；
采样窗口长度即采样窗口中包含的采样点的数量NL＝round(Nm×Nw,0)。
4.如权利要求3所述的谐波分析方法，其特征在于：Threshold＝0.3。
5.如权利要求1所述的谐波分析方法，其特征在于：在执行步骤4之前，先将计算所需的三角函数值存储到查询表tcos中，具体步骤为：设角度缩小倍数Amp′，Amp′＜Amp，以Amp′/Amp为分辨率取0度至450度的所有角度值；对于各角度α，以α*Amp/Amp′作为查询表的序号值，以cosα)为该序号值n＝α*Amp/Amp′对应的三角函数值，即tcos[n]＝cos(α)。
6.如权利要求5所述的谐波分析方法，其特征在于：步骤4的具体步骤为：
步骤4-1、分配用于存放向量的长度为NL的空间，对向量/>的各元素进行初始化：
Di＝Mi i＝0,1,2,…,NL-1
步骤4-2、分配两个长度均为NL的空间，分别用于存放实部转换向量和虚部转换向量/>对实部转换向量/>和虚部转换向量/>进行初始化：
TRi＝tcos[Di/Amp′] i＝0,1,2,…,NL-1
步骤4-3、设k＝Indexstart，以Rek、Imk分别表示第k次谐波含量的实部和虚部，则第k次谐波含量的实部和虚部的计算方式如下：
式中为用于存放计算输入信号的输入信号向量，该向量长度为NL；
步骤4-4、将向量与向量/>相加并保存至向量/>
谐波次数k递增：k＝k+interval；
如果k>Indexstop，则退出，否则继续执行步骤4-5；
步骤4-5、遍历当前向量中的各元素，若某个元素大于360×Amp，则将该元素减去360×Amp；
步骤4-6、根据当前的刷新实部转换向量/>和虚部转换向量/>
TRi＝tcos[Di/Amp′] i＝0,1,2,…,NL-1
步骤4-7、计算第k次谐波含量的实部和虚部：
步骤4-8、跳转至步骤4-4。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210683887.3A CN114859120B (zh) | 2022-06-17 | 2022-06-17 | 谐波分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210683887.3A CN114859120B (zh) | 2022-06-17 | 2022-06-17 | 谐波分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114859120A CN114859120A (zh) | 2022-08-05 |
| CN114859120B true CN114859120B (zh) | 2023-09-15 |
Family
ID=82625412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210683887.3A Active CN114859120B (zh) | 2022-06-17 | 2022-06-17 | 谐波分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114859120B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0450668A (ja) * | 1990-06-12 | 1992-02-19 | Hioki Ee Corp | 高調波電力計 |
| CN102435844A (zh) * | 2011-11-01 | 2012-05-02 | 南京磐能电力科技股份有限公司 | 一种频率无关的正弦信号相量计算方法 |
| CN102680785A (zh) * | 2012-04-26 | 2012-09-19 | 中国人民解放军海军工程大学 | 基于自适应变窗的同步相量测量方法 |
| CN103033685A (zh) * | 2012-11-27 | 2013-04-10 | 国网智能电网研究院 | 基于分频段双路滤波的宽频带谐波采集、测量系统及方法 |
| CN104037728A (zh) * | 2014-05-05 | 2014-09-10 | 珠海许继电气有限公司 | 基于软件测频与谐波分析的配电线路单相接地保护控制方法 |
| CN104181391A (zh) * | 2014-08-21 | 2014-12-03 | 艾德克斯电子（南京）有限公司 | 数字功率计谐波检测的方法 |
| CN105846697A (zh) * | 2016-05-05 | 2016-08-10 | 苏州汇川技术有限公司 | 一种低开关频率下pwm整流控制方法以及系统 |
| CN108037361A (zh) * | 2017-12-05 | 2018-05-15 | 南京福致通电气自动化有限公司 | 一种基于滑动窗dft的高精度谐波参数估计方法 |
| CN108957133A (zh) * | 2018-08-13 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | 变电站电容器采样频率动态评估和选取方法 |
| CN109374966A (zh) * | 2018-10-23 | 2019-02-22 | 国网重庆市电力公司电力科学研究院 | 一种电网频率估计方法 |
| CN114252837A (zh) * | 2021-12-17 | 2022-03-29 | 常州同惠电子股份有限公司 | 用于数字功率计谐波测量的高精度均等取样系统及方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130173189A1 (en) * | 2011-12-29 | 2013-07-04 | University of Connecticut, a public institution of higher education | Robust high resolution spectrum estimation method for accurate phasor, harmonic and interharmonic measurement in power systems |
-
2022
- 2022-06-17 CN CN202210683887.3A patent/CN114859120B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0450668A (ja) * | 1990-06-12 | 1992-02-19 | Hioki Ee Corp | 高調波電力計 |
| CN102435844A (zh) * | 2011-11-01 | 2012-05-02 | 南京磐能电力科技股份有限公司 | 一种频率无关的正弦信号相量计算方法 |
| CN102680785A (zh) * | 2012-04-26 | 2012-09-19 | 中国人民解放军海军工程大学 | 基于自适应变窗的同步相量测量方法 |
| CN103033685A (zh) * | 2012-11-27 | 2013-04-10 | 国网智能电网研究院 | 基于分频段双路滤波的宽频带谐波采集、测量系统及方法 |
| CN104037728A (zh) * | 2014-05-05 | 2014-09-10 | 珠海许继电气有限公司 | 基于软件测频与谐波分析的配电线路单相接地保护控制方法 |
| CN104181391A (zh) * | 2014-08-21 | 2014-12-03 | 艾德克斯电子（南京）有限公司 | 数字功率计谐波检测的方法 |
| CN105846697A (zh) * | 2016-05-05 | 2016-08-10 | 苏州汇川技术有限公司 | 一种低开关频率下pwm整流控制方法以及系统 |
| CN108037361A (zh) * | 2017-12-05 | 2018-05-15 | 南京福致通电气自动化有限公司 | 一种基于滑动窗dft的高精度谐波参数估计方法 |
| CN108957133A (zh) * | 2018-08-13 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | 变电站电容器采样频率动态评估和选取方法 |
| CN109374966A (zh) * | 2018-10-23 | 2019-02-22 | 国网重庆市电力公司电力科学研究院 | 一种电网频率估计方法 |
| CN114252837A (zh) * | 2021-12-17 | 2022-03-29 | 常州同惠电子股份有限公司 | 用于数字功率计谐波测量的高精度均等取样系统及方法 |
Non-Patent Citations (2)
| Title |
|---|
| 减少频谱泄漏的一种新的等角度间隔采样递推算法;陈明凯 等;电工技术学报;第20卷(第08期);94-98 * |
| 陈明凯 等.减少频谱泄漏的一种新的等角度间隔采样递推算法.电工技术学报.2005,第20卷(第08期),94-98. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114859120A (zh) | 2022-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107748734B (zh) | 一种解析-经验模态分解方法 | |
| CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及系统 | |
| CN103941087B (zh) | 欠采样速率下的高频余弦信号的频率测量方法及其装置 | |
| CN103207319A (zh) | 数字化变电站电力信号非同步采样条件下的谐波测量方法 | |
| CN104897961A (zh) | 基于互乘法窗函数的三谱线插值fft谐波分析方法及系统 | |
| Paris et al. | The sinefit spectral envelope preprocessor | |
| CN109633266B (zh) | 一种频率测量方法、系统、装置及计算机可读存储介质 | |
| CN103119453B (zh) | 基于二次型的数字频率估计 | |
| CN111984920B (zh) | 次/超同步谐波参数识别方法、装置、设备和介质 | |
| CN108776263B (zh) | 基于高阶汉宁自卷积窗及改进插值算法的谐波检测方法 | |
| CN114460527B (zh) | 关联度延拓Hilbert移相电子式互感器校验仪溯源方法和系统 | |
| CN104181391A (zh) | 数字功率计谐波检测的方法 | |
| CN115575707A (zh) | 基于改进fft算法与小波变换结合的谐波检测装置及方法 | |
| CN114859120B (zh) | 谐波分析方法 | |
| US8223057B2 (en) | Quantizing sampled inputs using fixed frequency analog to digital conversions through interpolation | |
| CN110068729A (zh) | 一种信号相量计算方法 | |
| CN114113864A (zh) | 一种频响测量用单点采样优化方法及系统 | |
| Roscoe et al. | Fast-responding measurements of power system harmonics using discrete and fast fourier transforms with low spectral leakage | |
| US7444249B2 (en) | Digital power metering system for reducing spectral leakage when determining the frequencies present in a power signal | |
| CN103543331A (zh) | 一种计算电信号谐波和间谐波的方法 | |
| CN102539913B (zh) | 高次谐波成分测定装置 | |
| CN107942139A (zh) | 一种新型电力谐波参数软件同步采样方法 | |
| CN111999635A (zh) | 一种基于4项5阶Nuttall窗的板卡故障信号分析方法及终端 | |
| CN111444613A (zh) | 一种基于改进emd电力系统谐波分析方法、装置及存储介质 | |
| CN105550516A (zh) | 基于Adams算法的频率自适应谐波间谐波分析方法 |
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 |