CN104407197B - 一种基于三角函数迭代的信号相量测量的方法 - Google Patents
一种基于三角函数迭代的信号相量测量的方法 Download PDFInfo
- Publication number
- CN104407197B CN104407197B CN201410706261.5A CN201410706261A CN104407197B CN 104407197 B CN104407197 B CN 104407197B CN 201410706261 A CN201410706261 A CN 201410706261A CN 104407197 B CN104407197 B CN 104407197B
- Authority
- CN
- China
- Prior art keywords
- component
- iteration
- signal
- electric power
- power 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.)
- Expired - Fee Related
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000005070 sampling Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000012804 iterative process Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241001274660 Modulus Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种基于三角函数迭代的信号相量测量的方法,在时域对信号进行采样得到N点离散序列,采用准同步采样算法估计采样信号的基波频率,根据三角函数和角公式对时域信号进行三角基函数分解,构建离散序列信号矩阵模型,利用基波频率估计值初始化矩阵模型参数,选取分量的一阶导数的最大值作为迭代下降因子,建立迭代调整方程,根据并行迭代算法求出信号矩阵模型系数,从而得到信号相量的测量结果。基于三角函数迭代的信号相量测量的方法可以有效避免传统傅里叶变换的频谱泄漏,选择一阶导数最大值构建迭代调整方程,算法计算简便,能确保收敛过程快速稳定。
Description
技术领域
本发明涉及信号相量测量领域,具体是一种基于三角函数迭代的信号相量测量的方法。
背景技术
改革开放以来,我国的经济得到高速发展,与之相关的电力工业也得到了长足的进步。随着电力电子技术的发展,非线性的电力电子装置、半导体器件在电力系统中得到广泛应用,但是非线性设备的使用带来的谐波问题,对电力系统的稳定和电能质量构成的威胁愈发严重,而电力信号相量实时和精确测量的实现,能有效解决谐波问题带来的威胁,对电力系统的监视、分析、控制以及保护等也能起到很大的改善。
当前主要的信号相量测量算法可以分为三大类:(1)傅里叶变换法,包括离散傅里叶变换法、递归傅里叶变换法和差分傅里叶变换法,其特点算法稳定性好,准确度较高,但在系统频率变化较快时易产生较大的误差;(2)过零检测法,其特点是原理比较简单,易于在DSP(数字信号处理器)和嵌入式等硬件设计上实现,但实时性较差,易受谐波、噪声分量的影响;(3)参数优化估计法,包括其特点是检测精度高,收敛速度快,收敛过程稳定,但是计算量大,对硬件计算设备来说,处理的数据还是比较多。
傅里叶变换算法、过零检测法和参数估计法这三种方法各有自己独特的优缺点,在实际测量工作中使用广泛。但是由于DSP和嵌入式等数据处理硬件水平的提高,研究一种计算简便、精确度高、易于实现的信号相量测量方法对电力系统的监控和保护具有重要的意义。
发明内容
本发明的目的是为了克服上述信号相量测量方法的不足,提出一种基于三角函数迭代的信号相量测量的方法,可以有效避免传统傅里叶变换的频谱泄漏,选择一阶导数最大值构建迭代调整方程,算法计算简便,能确保收敛过程快速稳定,从而实现信号相量的准确测量。
为解决上述技术问题,本发明提出的解决方案为:在时域对信号进行采样得到N点离散序列,采用准同步采样算法估计采样信号的基波频率,根据三角函数和角公式对时域信号进行三角基函数分解,构建离散序列信号矩阵模型,利用基波频率估计值初始化矩阵模型参数,选取分量的一阶导数的最大值作为迭代下降因子,建立迭代调整方程,根据并行迭代算法求出信号矩阵模型系数,从而得到信号相量的测量结果。
本发明提出的基于三角函数迭代的信号相量测量的方法包括以下步骤:
步骤一:选择合适的采样频率fs和采样长度N对信号进行采样,采样频率fs应不低于信号最高频率的2倍,得到信号离散采样序列U(n),n=0,1,2,…,N-1;
步骤二:采用准同步采样算法对采样序列U(n)进行基波频率估计,得到基波频率的估计值f0,准同步采样算法单次迭代的采样点数H一般不少于20,迭代次数L一般不超过9次;
步骤三:利用三角函数和角公式对时域离散序列U(n)进行三角基函数分解,构建矩阵测量模型U=ZD+ES+GC,其中U为模型输出向量,D、S、C分别为直流分量、正弦分量和余弦分量,Z、E、G分别为各分量系数,利用基波频率估计值f0初始化矩阵测量模型的各项参数,并作为迭代的初始值;
步骤四:选取分量的一阶导数的最大值作为迭代下降因子,对基波频率f及矩阵测量模型参数直流分量D、正弦分量S、余弦分量C分别构造迭代调整方程,并行迭代求出参数值,并更新矩阵测量模型U的各项参数,迭代调整方程为:
其中△U为迭代方程的估计准则,η为下降系数;
步骤五:每次更新参数后,计算估计准则△U<ε是否满足条件,满足条件退出迭代过程,否则根据迭代方程循环迭代,其中ε是本方法预设的迭代允许误差,根据IEC61000系列标准中的信号相量测量精度要求,ε的取值可以设置为不超过10-20;
步骤六:当△U<ε成立时,迭代过程结束,由测量模型系数矩阵中的正弦分量S、余弦分量C,计算出信号的幅值Au和相位从而获得信号相量的测量结果。
所述的方法,步骤三中初始化直流分量D=0,直流分量系数Z=[1,1,...,1...]N,正弦分量S=[0,1,0,1,...,0,1...]M,余弦分量C=[1,0,1,0,...,1,0...]M,其中M为信号所含谐波最高次数,下标代表数据的长度。令正弦分量和余弦分量系数E、G可分别表示为
所述的方法,步骤四中定义误差函数ζ=U(n)-U,则估计准则△U表示为下降系数η应根据实际迭代参数确定,选取分量的一阶导数的最大值作为迭代下降因子,能确保收敛过程快速稳定。
所述的方法,步骤六中信号的幅值Au和相位可由和求出,S[0]、C[0]表示一维矩阵的第一个元素,信号相量可表示为
本发明基于三角函数迭代的信号相量测量的方法,可以有效避免传统傅里叶变换的频谱泄漏,选择一阶导数最大值构建迭代调整方程,算法计算简便,能确保收敛过程快速稳定,在高精度要求的条件下,算法计算简便,且在实际的测量工作中,易于在嵌入式和DSP中编程实现。
附图说明
图1是本发明处理流程的原理框图;
图2是本发明中实现基于三角函数迭代的信号相量测量方法的程序流程图;
图3是本发明中准同步采样算法迭代示意图。
具体实施方式
本发明提出了一种基于三角函数迭代的信号相量测量的方法。以下结合附图作详细说明:
本实施例的处理流程的原理框图如图1所示,信号通过调理电路后经过A/D(模数转换)后转化为数字量,并送入DSP处理器进行处理,DSP处理器完成三角函数迭代算法的运算,得到各个测量参数,从而完成信号相量的测量。
如图2,本实施例的一种基于三角函数迭代的信号相量测量的方法的具体步骤如下:
1.以采样频率fs=1000Hz对信号进行采样,一个包含有多项整数次谐波的时域信号可表示为:
其中,U0为信号的直流分量,M为信号所含有谐波的最高次数;m为谐波次数;f为基波频率;Um和分别为第m次谐波的幅值和初相角。本实施例中信号基波频率为50.5Hz,基波幅值为220v,初相角为0.05°;
通过三角函数和角公式分解式(1),令采样可得到N=500个样本的信号离散信号U(n),n=0,1,2,…,N-1:
2.对获得的采样序列U(n)采用准同步采样算法,迭代过程如图3所示,单次迭代点数H为20点,迭代次数L为5次,得到基波频率的估计值f0=50.499999999999993;
3.利用三角函数和角公式对时域信号进行三角基函数分解,对信号离散信号U(n)构建矩阵模型,可得:
U=ZD+ES+GC (3)
式中,其中U为模型输出向量;D、S、C分别为直流分量、正弦分量和余弦分量;Z、E、G分别为各分量系数;
4.根据构建的矩阵模型,初始化直流分量D=0,直流分量系数Z=[1,1,...,1...]N,余弦分量S=[0,1,0,1,...,0,1...]M,正弦分量C=[1,0,1,0,...,1,0...]M,其中下标代表数据的长度。令根据式(4)、(5)可以计算得出正弦分量系数E和余弦分量系数G的初始参数值,从而得出模型输出相量U:
定义误差函数ζ=U(n)-U,估计准则△U可由式(6)计算得出:
5.为使收敛过程快速稳定,选取分量的一阶导数的最大值作为迭代下降因子,对基波频率f及矩阵测量模型参数直流分量D、正弦分量S、余弦分量C分别构造迭代调整方程:
由迭代调整方程可迭代更新f、D、S、C的值,并由新的估计值由式(3)、(4)、(5)、(6)重新计算矩阵模型中E、G、U和估计准则△U的参数值;
6.令迭代允许误差ε=10-20,每次更新完参数后,判断估计准则是否满足△U<ε的精度要求,满足条件退出迭代过程,否则退回到步骤5,本实施例估计准则精度达到9.055470139575152e-21;
7.当△U<ε成立,迭代过程结束,由正弦分量和余弦分量的第一个元素根据式(11)和(12)得出信号的幅值和相位:
本实施例中,Au=2.199999999999995e+02,则信号相量可表示为
本实施例迭代次数为197次,迭代用时0.032925s,可见运算量小且快速收敛。
至此,完成了信号相量的测量。
Claims (4)
1.一种基于三角函数迭代的信号相量测量的方法,电力信号通过调理电路后经过模数转换后转化为数字量,并送入DSP处理器进行处理,DSP处理器完成三角函数迭代算法的运算,得到各个测量参数,从而完成电力信号相量的测量,其特征在于,DSP处理器完成三角函数迭代算法的运算包括以下步骤:
步骤一:选择合适的采样频率fs和采样长度N对电力信号进行采样,采样频率fs应不低于电力信号最高频率的2倍,得到电力信号离散采样序列U(n),n=0,1,2,···,N-1;
步骤二:采用准同步采样算法对采样序列U(n)进行基波频率估计,得到基波频率的估计值f0,准同步采样算法单次迭代的采样点数H不少于20,迭代次数L不超过9次;
步骤三:利用三角函数和角公式对时域离散序列U(n)进行三角基函数分解,构建矩阵测量模型U=ZD+ES+GC,其中U为模型输出向量,D、S、C分别为直流分量、正弦分量和余弦分量,Z、E、G分别为各分量系数,利用基波频率估计值f0初始化矩阵测量模型的各项参数,并作为迭代的初始值;
步骤四:选取直流分量、正弦分量和余弦分量的一阶导数的最大值作为迭代下降因子,对基波频率f及矩阵测量模型参数直流分量D、正弦分量S、余弦分量C分别构造迭代调整方程,并行迭代求出参数值,并更新矩阵测量模型U的各项参数,迭代调整方程为: 其中ΔU为迭代方程的估计准则,h为下降系数;
步骤五:每次更新参数后,计算估计准则ΔU<ε是否满足条件,满足条件退出迭代过程,否则根据迭代方程循环迭代,其中ε是本方法预设的迭代允许误差,根据IEC61000系列标准中的电力信号相量测量精度要求,ε的取值设置为不超过10-20;
步骤六:当ΔU<ε成立时,迭代过程结束,由测量模型系数矩阵中的正弦分量S、余弦分量C,计算出电力信号的幅值Au和相位,从而获得电力信号相量的测量结果。
2.根据权利要求1所述的方法,其特征在于,步骤三中初始化直流分量D=0,直流分量系数Z=[1,1,...,1...]N,正弦分量S=[0,1,0,1,...,0,1...]M,余弦分量C=[1,0,1,0,...,1,0...]M,其中M为电力信号所含谐波最高次数,下标代表数据的长度。令正弦分量和余弦分量系数E、G可分别表示为
3.根据权利要求1所述的方法,其特征在于,步骤四中定义误差函数ζ=U(n)-U,则估计准则ΔU表示为下降系数η应根据实际迭代参数确定,选取直流分量、正弦分量和余弦分量的一阶导数的最大值作为迭代下降因子,能确保收敛过程快速稳定。
4.根据权利要求1所述的方法,其特征在于,步骤六中电力信号的幅值Au和相位可由和求出,S[0]、C[0]表示一维矩阵的第一个元素,电力信号相量可表示为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410706261.5A CN104407197B (zh) | 2014-11-27 | 2014-11-27 | 一种基于三角函数迭代的信号相量测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410706261.5A CN104407197B (zh) | 2014-11-27 | 2014-11-27 | 一种基于三角函数迭代的信号相量测量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104407197A CN104407197A (zh) | 2015-03-11 |
CN104407197B true CN104407197B (zh) | 2017-06-27 |
Family
ID=52644842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410706261.5A Expired - Fee Related CN104407197B (zh) | 2014-11-27 | 2014-11-27 | 一种基于三角函数迭代的信号相量测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104407197B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106406830B (zh) * | 2015-07-29 | 2021-08-03 | 腾讯科技(深圳)有限公司 | 一种准周期信号的处理方法、装置和计算机可读存储介质 |
CN107329732B (zh) * | 2016-04-29 | 2021-07-16 | 中科寒武纪科技股份有限公司 | 一种用于执行多种超越函数运算的装置和方法 |
CN107766293B (zh) * | 2017-09-28 | 2021-04-23 | 中国计量科学研究院 | 部分采样数据规则性缺失时的信号频谱分析方法及系统 |
CN114965367B (zh) * | 2022-06-14 | 2024-06-25 | 广东工业大学 | 一种用于光学层析测量的混叠正弦波信号分离方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1332313A1 (ru) * | 1986-01-22 | 1987-08-23 | Кировский Политехнический Институт | Устройство дл вычислени тригонометрических функций |
CN102004186A (zh) * | 2010-11-08 | 2011-04-06 | 北京航空航天大学 | 一种消除频谱泄漏的高准确度正弦信号测量方法与装置 |
CN102095933A (zh) * | 2010-12-07 | 2011-06-15 | 天津大学 | 一种测相装置及其控制方法 |
CN102135570A (zh) * | 2011-02-25 | 2011-07-27 | 上海思源弘瑞自动化有限公司 | 用于智能变电站的同步相量测量方法及其装置 |
CN103869162A (zh) * | 2014-03-05 | 2014-06-18 | 湖南大学 | 一种基于时域准同步的动态信号相量测量方法 |
-
2014
- 2014-11-27 CN CN201410706261.5A patent/CN104407197B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1332313A1 (ru) * | 1986-01-22 | 1987-08-23 | Кировский Политехнический Институт | Устройство дл вычислени тригонометрических функций |
CN102004186A (zh) * | 2010-11-08 | 2011-04-06 | 北京航空航天大学 | 一种消除频谱泄漏的高准确度正弦信号测量方法与装置 |
CN102095933A (zh) * | 2010-12-07 | 2011-06-15 | 天津大学 | 一种测相装置及其控制方法 |
CN102135570A (zh) * | 2011-02-25 | 2011-07-27 | 上海思源弘瑞自动化有限公司 | 用于智能变电站的同步相量测量方法及其装置 |
CN103869162A (zh) * | 2014-03-05 | 2014-06-18 | 湖南大学 | 一种基于时域准同步的动态信号相量测量方法 |
Non-Patent Citations (1)
Title |
---|
一种高精度的电力系统相量测量方法;王超 等;《电力系统保护与控制》;20091001;第37卷(第19期);第36-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104407197A (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN101806832B (zh) | 一种低频率信号的频率测量方法 | |
CN104535836B (zh) | 电力信号的基波频率测量方法和系统 | |
CN104407197B (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN103399213B (zh) | 一种基于基波频率修正的高压设备介质损耗因数计算方法 | |
Tu et al. | CMF signal processing method based on feedback corrected ANF and Hilbert transformation | |
CN109633262A (zh) | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 | |
CN105629060B (zh) | 基于最优基带滤波的电网频率测量方法和装置 | |
CN102955068B (zh) | 一种基于压缩采样正交匹配追踪的谐波检测方法 | |
CN104316768B (zh) | 一种三相不平衡扰动源定位的负序阻抗参数估算方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN106018956B (zh) | 一种加窗谱线插值的电力系统频率计算方法 | |
Radonjic et al. | Stochastic measurement of power grid frequency using a two-bit A/D converter | |
CN104483539B (zh) | 一种基于泰勒展开式的有功功率快速测量方法 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN105486921A (zh) | 凯撒三阶互卷积窗三谱线插值的谐波与间谐波检测方法 | |
CN104181391A (zh) | 数字功率计谐波检测的方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN104569581B (zh) | 一种电网频率测量的多水平集单周期估计方法 | |
CN105067882B (zh) | 一种电量幅值的测定方法 | |
CN110068729A (zh) | 一种信号相量计算方法 | |
CN109239463B (zh) | 一种基于线性修正算法的介质损耗测量方法 | |
CN107167658B (zh) | 一种高精度抗干扰的电力系统基波频率及相位测量方法 | |
CN103293379B (zh) | 基于有效值的apf谐波检测方法及其直流侧电压的控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170627 |