CN101403774B - 一种基于非同步采样的谐波分析方法 - Google Patents

一种基于非同步采样的谐波分析方法 Download PDF

Info

Publication number
CN101403774B
CN101403774B CN200810194747XA CN200810194747A CN101403774B CN 101403774 B CN101403774 B CN 101403774B CN 200810194747X A CN200810194747X A CN 200810194747XA CN 200810194747 A CN200810194747 A CN 200810194747A CN 101403774 B CN101403774 B CN 101403774B
Authority
CN
China
Prior art keywords
tau
cos
sin
frequency
signal
Prior art date
Application number
CN200810194747XA
Other languages
English (en)
Other versions
CN101403774A (zh
Inventor
王柏林
Original Assignee
江苏中凌高科技有限公司
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 江苏中凌高科技有限公司 filed Critical 江苏中凌高科技有限公司
Priority to CN200810194747XA priority Critical patent/CN101403774B/zh
Publication of CN101403774A publication Critical patent/CN101403774A/zh
Application granted granted Critical
Publication of CN101403774B publication Critical patent/CN101403774B/zh

Links

Abstract

本发明涉及信号谐波检测技术领域内的一种基于非同步采样的谐波分析方法,包括如下步骤:1)将预设的若干与信号频率一一对应的权系数依次存入存储器中;2)对谐波用恒定不变的采样周期进行采样,采样信号经模数转换后输出给数字信号处理器;3)数字信号处理器中将从多路模数转换器来的数字信号进行处理,根据其频率,确定其从存储器中所需要调用的权系数的位置;4)在数字信号处理器中,将从存储器中调用的相应的权系数与输入的数字信号相乘后进行FFT处理;5)将FFT处理的结果进行显示或打印。该方法克服了谐波分析时的同步偏差问题,获得了实时、高精度谐波分析信号。可应用于机械、电力、信息、地质、水文等领域的信号进行谐波分析。

Description

一种基于非同步采样的谐波分析方法
技术领域
[0001] 本发明涉及一种信号谐波分析方法,它适用于频率变化的‘近似周期信号’(比如 电力信号)的谐波分析。
背景技术
[0002] 机械、电力、信息、地质、水文等领域常常要对信号(或波形)进行谐波分析,能完 成谐波分析的仪器称为谐波分析仪(或谐波分析仪)。谐波分析方法的主要工作包括:对 连续时间信号采样和ADC(多路模数转换器)得到采样序列、对采样序列进行FFT(快速富 里叶变换)得到谐波分析结果(各次谐波的幅值和相位信息)、谐波分析结果的显示和传送 (通信)等。
[0003] 谐波分析仪的精度主要取决于两方面:一是ADC的精度,二是采样的‘同步’程 度——最好是采样周期乘采样点数刚好等于整数个信号周期——即‘同步采样’。实际上, 同步采样是很难实现的,因为信号(比如电力信号)的频率始终是变化的,并且数字信号处 理器(即DSP)的定时器总是有最小分辨率的。
[0004] 为了在信号频率变化时维持采样的同步或近似同步,大多数谐波分析方法都设置 了 ‘同步环节’,同步环节的作用就是根据信号频率的变化实时调整采样周期的长短。同步 环节只能在一定程度上减小同步偏差,不可能实现真正的同步采样,因为:
[0005] (1)同步环节基于频率跟踪,频率跟踪总是有时延的;
[0006] (2) DSP的定时器有最小分辨率,调整采样周期的精度受最小分辨率的限制;
[0007] (3)高档的谐波分析仪采样频率很高(几KHz到几十KHz)——采样周期很短(微 秒级),
[0008] 由DSP来精确调整采样周期是很困难的;
[0009] (4)实际工程中信号的频率无时无刻不在变化。
发明内容
[0010] 本发明的目的是提供一种克服同步偏差的新方法——一种基于非同步采样的谐 波分析方法,以提高谐波分析的精度。
[0011] 为了实现上述目的,本发明提供的一种基于非同步采样的谐波分析方法,包括如 下步骤:
[0012] 1)将预设的若干与信号频率一一对应的权系数依次存入存储器中;所述权系数 是指在一定的范围内将信号频率细分为K个离散值f2,…,fk,针对每一个离散值fk代 入方程FWG = I2M+1中,并离线解出对应的权系数Wk(k = 1,2,…,K) ;W是加权矩阵,W = (1土叫[奴1)奴2)"1洲]^个实数《(11) (η = 1,2,…,N)称为一组权系数;
[0013] 2)对谐波用恒定不变的采样周期进行采样,采样后得到的模拟信号送入多路模数 转换器(ADC)中进行模数转换,再将模数转换后得到的采样序列(数字信号)输出给数字 信号处理器(DSP);[0014] 3)数字信号处理器中将从多路模数转换器来的数字信号进行处理,根据其频率, 确定其从存储器中所需要调用的权系数的位置;
[0015] 4)在数字信号处理器中,将从存储器中调用的相应的权系数Wk与输入的数字信号 相乘后进行FFT处理;
[0016] 5)将FFT处理的结果通过数字信号处理器的通讯接口输出给显示设备或打印设 备进行显示或打印。
[0017] 式中
"1 1 1 … 1 : 2 2 cos r0 2cos(2r0) ... 2cos[(7V-1)t0]
0 2 sin τ n 2sin(2r。)…2 sin[(jV - 1)τ0 ]
[0018]<formula>formula see original document page 4</formula>[0019] <formula>formula see original document page 4</formula>
[0020] τ 0 = 2S π /N, τ = τ 0f/f0
[0021] 其中,F e R(2M+1)xn是一个常数矩阵、G e Rnx(2M+1)是一个随信号频率的变化而变化 的矩阵,I2M+1是(2M+1)维单位阵;S为正整数,M是谐波的最高次数,f是信号的基波频率; fo定义为‘同步频率’——同步频率fo满足:
[0022] 1) f0十分接近信号的额定频率或平均频率,2)f = f0时刚好达到同步采样,f0是常 数。
[0023] 本发明中,权系数是预先设定并存入存储器中的;根据连续等间隔采集的N个采 样数据y (1),y (2),…,y (N),计算出当前信号频率f,当前信号频率满足fV^f彡fk,就从 W表中取出wk,再对这N个采样数据分别‘加权’,对加权后的N个数据进行标准FFT处理, 从FFT处理的结果求出各次谐波参数&和九。该方法是一种全新的、实时、高精度谐波分析 方法,其结构上的显著特点是没有‘同步环节’,DSP控制ADC用恒定不变的采样周期τ。进 行采样;该方法可应用于机械、电力、信息、地质、水文等领域的信号进行谐波分析,例如对 电网的电压、电流进行谐波分析,提高了谐波分析精度,减小了同步误差;能够实现快速、实 时处理。
附图说明
[0024] 图1为本发明应用于三相电网谐波分析的工作原理图。 具体实施方式
[0025] 下面结合附图具体介绍一种基于非同步采样的三相电网的谐波分析方法,其步骤如下:
[0026] 1)将预设的若干与信号频率一一对应的权系数依次存入存储器中;所述权系数 是指在一定的范围内将信号频率细分为K个离散值f2,…,fk,针对每一个离散值fk代 入方程FWG = I2M+1中,并离线解出对应的权系数Wk (k = 1,2, -,K);
[0027] 2)对谐波用恒定不变的采样周期进行采样,对三相电压Ua、Ub、Uc以相应的电压 互感器Pta、PTb、PTc进行采样,对三相电流la、lb、Ic以电流互感器Cta、CTb、CTc进行采 样,采样后得到的信号送入多路模数转换器中进行模数转换,再将模数转换后得到的采样 序列输出给数字信号处理器;
[0028] 3)数字信号处理器中将从多路模数转换器来的数字信号进行处理,根据其频率, 确定其从存储器中所需要调用的权系数的位置;
[0029] 4)在数字信号处理器中,将从存储器中调用的相应的权系数Wk与输入的数字信号 相乘后进行FFT处理;
[0030] 5)将FFT处理的结果通过数字信号处理器的通讯接口输出给显示设备或打印设 备进行显示或打印。
[0031] 下面对上述步骤进行进一步阐述:
[0032] 本发明提出一种基于非同步采样的谐波分析方法,
[0033] y = [y(l)y(2)-y(N)]T
[0034] 如果直接对采样序列y进行FFT (快速富里叶变换)处理,就会产生众所周知的 ‘同步误差’,因为现在采用的是非同步采样方式。
[0035] 为了消除同步误差,我们不是直接对采样序列进行FFT,而是首先对采样序列进行 预处理——加权:
[0036] 加权时,预先构造若干权系数,预先构造权系数的办法是:在一定的范围内将信号 频率细分为K个离散值&,f1;…,fK,针对每一个离散值fk(k= 1,2,…,K)离线算出Wk(k =1,2,…,K),将K组权系数按一定规律存在DSP的flash存储器中。
[0037] 可以证明,只要权系数满足:
[0038] FffG = I2M+1 (1)
[0039] 从y的FFT结果求出的信号各次谐波的幅值和相角的估计值就是无偏的——与同 步采样时得到的结果一样,式中
Ί ι ι … ι “
2 2 cos r0 2cos(2r0) ··· 2cos[(jV-l)r0]
0 2 sin r0 2sin(2T0) ... 2sin[(JV —1)τ0]
[0040] F-. . . . .
2 2Cos(Mr0) 2cos(2Mr0) ··· 2cos[(iV-l)Mr0] 0 2 Sin(Mr0) 2sin(2M*r0) ···[0041] <formula>formula see original document page 6</formula>
[0043] 其中,F e R(2M+1)xn是一个常数矩阵、G e Rnx(2M+1)是一个随信号频率的变化而变化 的矩阵,12„+1是(2M+1)维单位阵。S为正整数,M是谐波的最高次数,f是信号的基波频率。 f0定义为‘同步频率’——同步频率fo满足=Dftl十分接近信号的额定频率或平均频率,2) f = fo时刚好达到同步采样,fo是常数。
[0044] 可以证明,只要N彡(2M+1)2,方程(1)中W就一定有解。满足方程(1)的权系数 不是常数,因为信号频率f变化时矩阵G会随之变化,而F和I2M+1都是常数矩阵,显然W也 随信号频率的变化而变化。
[0045] 然后从存储器中调用的相应的权系数Wk与输入的数字信号相乘,公式为:
[0046] y = Wy (2)
[0047]式中,W 是加权矩阵,W = diag[w(l) w⑵…w(N)],N 个实数 w(η) (η = 1,2,…, N)称为一组权系数。公式(2)的过程就是从存储器中调用权系数的过程。
[0048] 然后,对加权序列y进行FFT处理:
[0049] X = fft (y) (3)
[0050] 进而,从N维复向量X就可以求出信号各次谐波的幅值和相角的估计值。然后,进 行显示或打印。
[0051] 这种基于非同步采样的谐波分析方法运行时,DSP对采样控制和信号处理的步骤 是:
[0052] (1)控制ADC用恒定不变的采样周期进行采样,一次连续采样N点;
[0053] (2) 一次采样结束后测出当时的信号频率f ;
[0054] (3)查权系数表读出当前的权系数,如果当前信号频率满足fk—^f ^ fk,按预定规 则从权系数表中读出相应的wk(η) (η = 1,2,…,N);
[0055] (4)对 N 个采样数据加权,即 y (n) = wk(n) Xy(n) (η = 1,2,…,N);
[0056] (5)对加权后的N个数据y (1)、y (2) · · · y (N)执行FFT。
[0057] (6)从FFT结果求出各次谐波幅值和相位。
[0058] 步骤(1)〜(6)反复循环执行。显然,权系数的频率细分越密、测频精度越高,谐 波分析精度越高。仿真表明:如果频率细分间隔为0. 5Hz、测频精度为0. 1%,谐波分析精度 将优于0.1%。
[0059] 总结起来,本发明提出的基于非同步采样的谐波分析方法,有如下特点:
[0060] 1)方法简单、易实现;
[0061] 2)采样控制既简单又精确——可以不受DSP定时器的最小分辨率的影响并且采 样周期恒定不变;
[0062] 3)谐波分析的精度很高——几乎完全消除了同步误差;
[0063] 与标准FFT相比,DSP增加的数据处理工作量不大——与传统的‘加窗’差不多一能够实现快速、实时处理。
[0064] 本发明并不局限于上述实施例,该谐波分析方法还可以应用到其他 领域对相应的 其他谐波进行实时分析。

Claims (1)

  1. 一种基于非同步采样的谐波分析方法,包括如下步骤:1)将预设的若干与信号频率一一对应的权系数依次存入存储器中;所述权系数是指在一定的范围内将信号频率细分为K个离散值f1,f2,…,fk,针对每一个离散值fk代入方程FWG=I2M+1中,并离线解出对应的权系数Wk(k=1,2,…,K);W是加权矩阵,W=diag[w(1)w(2)…w(N)],N个实数w(n)(n=1,2,…,N)称为一组权系数;2)对谐波用恒定不变的采样周期进行采样,一次连续采样N点,采样后得到的模拟信号送入多路模数转换器中进行模数转换,再将模数转换后得到的采样序列输出给数字信号处理器;3)数字信号处理器中将从多路模数转换器来的数字信号进行处理,根据其频率,确定其从存储器中所需要调用的权系数的位置;4)在数字信号处理器中,将从存储器中调用的相应的权系数Wk与输入的数字信号相乘后进行FFT处理;5)将FFT处理的结果通过数字信号处理器的通讯接口输出给显示设备或打印设备进行显示或打印;式中 <mrow> <mi>F</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>2</mn> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>2</mn> <mi>s</mi> <msub> <mi>in&tau;</mi> <mn>0</mn> </msub> </mtd> <mtd> <mn>2</mn> <mi>sin</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>2</mn> <mi>sin</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>2</mn> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <mrow> <mo>(</mo> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>2</mn> <mi>cos</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>2</mn> <mi>sin</mi> <mrow> <mo>(</mo> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>2</mn> <mi>sin</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>2</mn> <mi>sin</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>M</mi> <msub> <mi>&tau;</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> </mtable> </mfenced> </mrow> <mrow> <mi>G</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mi>M&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mi>M&tau;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mi>cos</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>M&tau;</mi> <mo>)</mo> </mrow> </mtd> <mtd> <mi>sin</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>M&tau;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>cos</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&tau;</mi> <mo>]</mo> </mtd> <mtd> <mi>sin</mi> <mo>[</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&tau;</mi> <mo>]</mo> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mi>cos</mi> <mo>[</mo> <mi>M</mi> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&tau;</mi> <mo>]</mo> </mtd> <mtd> <mi>sin</mi> <mo>[</mo> <mi>M</mi> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&tau;</mi> <mo>]</mo> </mtd> </mtr> </mtable> </mfenced> </mrow>τ0=2Sπ/N,τ=τ0f/f0其中,F∈R(2M+1)×N是一个常数矩阵、G∈RN×(2M+1)是一个随信号频率的变化而变化的矩阵,I2M+1是(2M+1)维单位阵;S为正整数,M是谐波的最高次数,f是信号的基波频率;f0定义为‘同步频率’——同步频率f0满足:1)f0十分接近信号的额定频率或平均频率,2)f=f0时刚好达到同步采样,f0是常数。
CN200810194747XA 2008-11-07 2008-11-07 一种基于非同步采样的谐波分析方法 CN101403774B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810194747XA CN101403774B (zh) 2008-11-07 2008-11-07 一种基于非同步采样的谐波分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810194747XA CN101403774B (zh) 2008-11-07 2008-11-07 一种基于非同步采样的谐波分析方法

Publications (2)

Publication Number Publication Date
CN101403774A CN101403774A (zh) 2009-04-08
CN101403774B true CN101403774B (zh) 2010-08-25

Family

ID=40537847

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810194747XA CN101403774B (zh) 2008-11-07 2008-11-07 一种基于非同步采样的谐波分析方法

Country Status (1)

Country Link
CN (1) CN101403774B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915874B (zh) * 2010-07-20 2013-01-09 北海市深蓝科技发展有限责任公司 一种基于傅立叶变换的谐波检测方法
CN102035554B (zh) * 2010-11-19 2014-04-23 钜泉光电科技(上海)股份有限公司 电能计量和分析系统、方法和模数转换电路
CN102353839B (zh) * 2011-07-18 2013-05-29 华北电力大学(保定) 基于多层前馈神经网络的电力系统谐波分析方法
CN103454494B (zh) * 2011-08-24 2016-01-20 常州顺创电气科技有限公司 一种高精度的谐波分析方法
CN103547328B (zh) * 2012-05-22 2015-12-23 深圳市英威腾电气股份有限公司 谐波检测方法及相关装置
KR101812915B1 (ko) * 2014-04-15 2017-12-27 미쓰비시덴키 가부시키가이샤 회로 차단기
CN103969508B (zh) * 2014-05-07 2016-06-15 江苏天浩达科技有限公司 一种实时高精密的电力谐波分析方法及装置
CN105137183A (zh) * 2015-09-15 2015-12-09 湖北工业大学 一种电力系统谐波分析方法及系统
CN105872380A (zh) * 2016-04-26 2016-08-17 乐视控股(北京)有限公司 一种视频图像采集方法及装置
CN105973272B (zh) * 2016-07-28 2018-01-19 深圳市康必达控制技术有限公司 一种提高电能计量精度的方法
CN106597099B (zh) * 2016-12-15 2019-02-01 国网新疆电力公司电力科学研究院 一种基于参数估计的谐波分析系统
CN108680787A (zh) * 2018-05-23 2018-10-19 成都玖锦科技有限公司 基于fpga的实时频谱分析方法

Also Published As

Publication number Publication date
CN101403774A (zh) 2009-04-08

Similar Documents

Publication Publication Date Title
Romano et al. Enhanced interpolated-DFT for synchrophasor estimation in FPGAs: Theory, implementation, and validation of a PMU prototype
Phadke et al. Synchronized phasor and frequency measurement under transient conditions
Asiminoaei et al. A digital controlled PV-inverter with grid impedance estimation for ENS detection
Testa et al. Interharmonics: Theory and modeling
Zhu Exact harmonics/interharmonics calculation using adaptive window width
US6700515B2 (en) Digitizer apparatus and semiconductor testing apparatus
Roscoe Exploring the relative performance of frequency-tracking and fixed-filter phasor measurement unit algorithms under C37. 118 test procedures, the effects of interharmonics, and initial attempts at merging P-class response with M-class filtering
CN105223418B (zh) 次同步和超同步谐波相量的测量方法及测量装置
Ren et al. Real-time power system frequency and phasors estimation using recursive wavelet transform
Lavopa et al. Real-time estimation of fundamental frequency and harmonics for active shunt power filters in aircraft electrical systems
Lin Inter-harmonic identification using group-harmonic weighting approach based on the FFT
Stenbakken A Wideband Sampling Wattmeter1
CN103454497B (zh) 基于改进加窗离散傅立叶变换的相位差测量方法
Zygarlicki et al. A reduced Prony's method in power-quality analysis—parameters selection
Zeng et al. Harmonic phasor analysis based on improved FFT algorithm
CN100561162C (zh) 一种虚拟振动台检测信号处理方法及其设备
Ma et al. Identification and tracking of harmonic sources in a power system using a Kalman filter
Wen et al. Simple interpolated FFT algorithm based on minimize sidelobe windows for power-harmonic analysis
JP4947173B2 (ja) 交流インピーダンス測定装置
CN103543335B (zh) 一种同步相量测量方法
CN101806832B (zh) 一种低频率信号的频率测量方法
CN102435844B (zh) 一种频率无关的正弦信号相量计算方法
CN104007316B (zh) 一种欠采样速率下的高精度频率测量方法及其测量仪
CN100353169C (zh) 测试电子元件的方法及其仪器
Zhan et al. A Clarke transformation-based DFT phasor and frequency algorithm for wide frequency range

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20090814

Address after: Postal code 168, Shu Dong Road, Weiyang District, Jiangsu, Yangzhou Province, China: 225008

Applicant after: Jiangsu Zhongling Hi-Tech Development Co., Ltd.

Address before: Postcode 168, Shu Dong Road, Yangzhou, Jiangsu Province, China: 225008

Applicant before: Yangzhou Zhongling High-tech Development Co., Ltd.

ASS Succession or assignment of patent right

Owner name: JIANGSU ZHONGLING HIGH SCIENCE AND TECHNOLOGY CO.,

Free format text: FORMER OWNER: YANGZHOU ZHONGLING HIGH TECHNOLOGY DEVELOPMENT CO., LTD.

Effective date: 20090814

C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee

Owner name: JIANGSU ZHONGLING HIGH-TECH CO., LTD.

Free format text: FORMER NAME: JIANGSU ZHONGLING HI-TECH DEVELOPMENT CO., LTD.

CP01 Change in the name or title of a patent holder

Address after: 225008, No. 168 East Shu Shu Road, Weiyang District, Jiangsu, Yangzhou

Patentee after: Jiangsu ATrend high-tech incorporated company

Address before: 225008, No. 168 East Shu Shu Road, Weiyang District, Jiangsu, Yangzhou

Patentee before: Jiangsu Zhongling Hi-Tech Development Co., Ltd.