CN111781415B - 一种适用于dsp控制器的交流电压有效值实时计算方法 - Google Patents

一种适用于dsp控制器的交流电压有效值实时计算方法 Download PDF

Info

Publication number
CN111781415B
CN111781415B CN202010663956.5A CN202010663956A CN111781415B CN 111781415 B CN111781415 B CN 111781415B CN 202010663956 A CN202010663956 A CN 202010663956A CN 111781415 B CN111781415 B CN 111781415B
Authority
CN
China
Prior art keywords
effective value
alternating
voltage
value
calculation
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
Application number
CN202010663956.5A
Other languages
English (en)
Other versions
CN111781415A (zh
Inventor
高翔
赵喜洋
程焱
闫新军
陈琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Aero Electric Co Ltd
Original Assignee
Shaanxi Aero Electric Co Ltd
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 Shaanxi Aero Electric Co Ltd filed Critical Shaanxi Aero Electric Co Ltd
Priority to CN202010663956.5A priority Critical patent/CN111781415B/zh
Publication of CN111781415A publication Critical patent/CN111781415A/zh
Application granted granted Critical
Publication of CN111781415B publication Critical patent/CN111781415B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/02Measuring effective values, i.e. root-mean-square values
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Measurement Of Current Or Voltage (AREA)

Abstract

本申请属于交流电压有效值计算技术领域,特别涉及一种适用于DSP控制器的交流电压有效值实时计算方法,所述方法包括:步骤S1、获取DSP控制器三相交流电压中某相的电压值;步骤S2、构建有效值计算模型,将所述DSP控制器三相交流电压中某相的电压值输入到所述有效值计算模型中,计算得到第一交流电压有效值;步骤S3、构建低通数字滤波模型,将所述第一交流电压有效值输入到所述低通数字滤波模型中,计算得到第二交流电压有效值。本申请根据均方根递推公式计算出交流电压有效值的粗略值,再对其进行二阶级联巴特沃斯数字低通滤波后输出交流电压有效值的准确值,该方法在极大程度减小DSP存储空间站用的情况下且保证了计算的高准确性以及快速性。

Description

一种适用于DSP控制器的交流电压有效值实时计算方法
技术领域
本申请属于交流电压有效值计算技术领域,特别涉及一种适用于DSP控制器的交流电压有效值实时计算方法。
背景技术
DSP控制器由于其体积小、精度高、成本低、功耗小、数据及程序存储量大等特点,广泛应用于航空电源控制器中。在航空电源控制器软件算法中常常将交流电压的有效值作为重要的参数指标之一,实现控制、转换或保护等功能。例如,汇流条功率控制器中,汇流条的三相交流电压有效值是判断接触器动作是否正确的重要条件之一,直接影响了机上电源的供电模式转换。因此如何实现交流有效值实时计算是航空电源控制器研发过程中的一项技术难点和要点。
目前,交流电压有效值的计算方大体分为两种:硬件计算方法和软件计算方法。硬件计算方法主要是通过有效值计算芯片进行计算,例如AD637、AD8436等,这种方法受限于已经具备的硬件资源。软件计算方法多对交流电压的整个周期电压值进行采样后求其均方根值,此种方法为了达到计算的准确性要保证尽可能多的采样次数,存储数据量大,占据存储空间,且需要一个周期的时间才能完全跟随电压响应,计算实时性差。
发明内容
为了解决上述技术问题,本申请提供了一种适用于DSP控制器的交流电压有效值实时计算方法,解决交流电压有效值实时计算过程中出现的实时性差的问题。
本申请提供的适用于DSP控制器的交流电压有效值实时计算方法,所述方法包括:
步骤S1、获取DSP控制器三相交流电压中某相的电压值;
步骤S2、构建有效值计算模型,将所述DSP控制器三相交流电压中某相的电压值输入到所述有效值计算模型中,计算得到第一交流电压有效值;
步骤S3、构建低通数字滤波模型,将所述第一交流电压有效值输入到所述低通数字滤波模型中,计算得到第二交流电压有效值。
优选的是,步骤S2中,所述构建有效值计算模型包括:
将第一次采样得到的单相交流电压表示为x(1),第一次计算得到的交流电压有效值表示为y(1),则第n次采样得到的单相交流电压表示为x(n),第n次计算得到的交流电压有效值表示为y(n);
根据均方根的定义可得:
y(1)2=x(1)2
Figure BDA0002579655240000021
Figure BDA0002579655240000022
Figure BDA0002579655240000023
由以上四个公式递推得到:
Figure BDA0002579655240000024
对上式差分方程两边同时进行Z变换,可得:
Figure BDA0002579655240000025
由上式可得,有效值计算模型为:
Figure BDA0002579655240000026
Z为变换算子。
优选的是,步骤S3中,所述低通数字滤波模型为二阶级联巴特沃斯低通滤波模型。
优选的是,步骤S3中,所述构建低通数字滤波模型包括:
二阶级联巴特沃斯低通数字滤波器是由两个一阶巴特沃斯低通数字滤波器组成,两个级联一阶巴特沃斯低通数字滤波器的传递函数分别表示为:
Figure BDA0002579655240000027
Figure BDA0002579655240000028
根据Z变换可得各系数与截止频率和采样频率间的关系可以表示为:
Figure BDA0002579655240000031
Figure BDA0002579655240000032
Figure BDA0002579655240000033
Figure BDA0002579655240000034
其中,fs为截止频率,fc为采样频率。
本申请提供的适用于DSP控制器的交流电压有效值实时计算方法,根据均方根递推公式计算出交流电压有效值的粗略值,再对其进行二阶级联巴特沃斯数字低通滤波后输出交流电压有效值的准确值。
本申请仅通过交流电压的当前周期采样值、上一周期采样值以及上两个周期的交流电压有效值就能够计算出交流电压有效值的准确值,对交流电压能够实时采样计算,DSP计算执行时间小于5μs且计算误差在0.5%以内,该方法实现了交流电压有效值的实时计算,在极大程度减小DSP存储空间站用的情况下且保证了计算的高准确性以及快速性。
附图说明
图1是本申请适用于DSP控制器的交流电压有效值实时计算方法涉及的有效值计算模型示意图。
图2是本申请适用于DSP控制器的交流电压有效值实时计算方法涉及的低通数字滤波模型示意图。
图3是单相交流电压供电示意图。
图4是有效值计算模型中不同n值的输出波形示意图。
图5是有效值计算结果示意图。
图6是快速响应的有效值计算结果局部图示意图。
图7是计算精度的有效值计算结果局部图示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施方式中的附图,对本申请实施方式中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施方式是本申请一部分实施方式,而不是全部的实施方式。下面通过参考附图描述的实施方式是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都属于本申请保护的范围。下面结合附图对本申请的实施方式进行详细说明。
本申请提供了一种适用于DSP控制器的交流电压有效值实时计算方法,解决交流电压有效值实时计算过程中出现的实时性差的问题。
本申请提供的适用于DSP控制器的交流电压有效值实时计算方法,包括:
步骤S1、获取DSP控制器三相交流电压中某相的电压值;
步骤S2、构建有效值计算模型,将DSP控制器三相交流电压中某相的电压值输入到有效值计算模型中,计算得到第一交流电压有效值;
步骤S3、构建低通数字滤波模型,将第一交流电压有效值输入到低通数字滤波模型中,计算得到第二交流电压有效值。
具体的,DSP控制器通过外部或者自带的ADC模块,获得三相交流电压中某相的电压值U。
将第一次采样得到的单相交流电压表示为x(1),第一次计算得到的交流电压有效值表示为y(1),则第n次采样得到的单相交流电压表示为x(n),第n次计算得到的交流电压有效值表示为y(n);
根据均方根的定义可得:
y(1)2=x(1)2
Figure BDA0002579655240000041
Figure BDA0002579655240000042
Figure BDA0002579655240000043
由以上四个公式递推得到:
Figure BDA0002579655240000051
可以看出仅需当前周期电压采样值、上一周期电压采样值以及上两周期的电压有效值即可计算出有效值大小,无需采集多周期电压值而占用大量存储空间。
对上式差分方程两边同时进行Z变换,可得:
Figure BDA0002579655240000052
由上式可得,有效值计算模型为:
Figure BDA0002579655240000053
Z为Z域的变换算子,Z变换是对离散序列进行的一种数学变换,可将时域信号(离散时间序列)变换为在复频域的表达式。
上式中系数
Figure BDA0002579655240000054
即图1中的K1,系数
Figure BDA0002579655240000055
即图1中系数K2,均与n有关且K1和K2存在数值关系K2=1-2K1
n的取值会同时影响计算精度和响应速度两个方面,为了保证在电源控制系统应用中计算的快速响应,有效值计算模型要尽可能的增加计算响应的快速性,并为在有效值计算后增加低通滤波模型来补偿有效值计算模块的精度。
本发明技术方案的适用于DSP控制器的交流电压有效值实时计算方法,由于巴特沃斯低通滤波器在通频带内外均有平稳的幅频特性,频率响应曲线最大限度平坦且没有纹波,且在阻频带逐渐下降为零,选用二阶级联巴特沃斯低通滤波器对有效值计算结果进行滤波处理得到最终电压有效值。
二阶级联巴特沃斯低通数字滤波器是由两个一阶巴特沃斯低通数字滤波器组成,两个级联一阶巴特沃斯低通数字滤波器的传递函数分别表示为:
Figure BDA0002579655240000056
Figure BDA0002579655240000061
上式中的系数K3、K4、K5、K6均与滤波器的截止频率和采样频率有关。数字滤波器的设计一般可由模拟滤波器转换为数字滤波器,即由S域变换到Z域。根据Z变换可得各系数与截止频率和采样频率间的关系可以表示为:
Figure BDA0002579655240000062
Figure BDA0002579655240000063
Figure BDA0002579655240000064
Figure BDA0002579655240000065
其中,fs为截止频率,fc为采样频率。
截止频率的选择会影响系统的响应速度和计算精度,截止频率越低响应速度越慢,计算精度越好,反之亦然。
在本申请的一个实施方式中,设定单相交流电源的电压有效值为115V,频率为400Hz,A/D采样频率为10KHz,供电时间为0.3s至0.5s共0.2s,即200ms,如图3所示。
为了平衡计算的快速性和准确性,对n进行从0到1每隔0.01的遍历。挑选具有代表性的4组数据,n的取值分别为4.65、9.01、46.51、89.93,如图4所示。从图中可以看出,n=4.65时,系统具有最佳的快速性,但准确度最差;n=89.93时,系统的计算精度最佳,但对电压变化的响应最慢,综合快速性与计算精度,选择n=9.09即K1=0.11。有效值计算模型的传递函数为:
Figure BDA0002579655240000066
根据第一部分有效值计算的效果,综合滤波精度和响应速度将第二部分低通数字滤波部分的系数分别设定为80Hz和100Hz,则:
Figure BDA0002579655240000071
即低通数字滤波模型的传递函数分别为:
Figure BDA0002579655240000072
Figure BDA0002579655240000073
根据上述传递函数计算最终的有效值结果如图5所示,选取中间一段数据局部放大如图6和图7所示,可以看出本方法有效值计算结果降至20V仅需13.4ms,计算精度在0.5%以内。
与现有交流电压有效值计算方法相比,本发明方案首先根据均方根递推公式计算出交流电压有效值的粗略值,再对其进行二阶级联巴特沃斯数字低通滤波后输出交流电压有效值的准确值。本发明的计算方法仅通过交流电压的当前周期采样值、上一周期采样值以及上两个周期的交流电压有效值就能够计算出交流电压有效值的准确值,对交流电压能够实时采样计算,DSP计算执行时间小于5μs且计算误差在0.5%以内。该方法实现了交流电压有效值的实时计算,在极大程度减小DSP存储空间站用的情况下且保证了计算的高准确性以及快速性。

Claims (3)

1.一种适用于DSP控制器的交流电压有效值实时计算方法,其特征在于,所述方法包括:
步骤S1、获取DSP控制器三相交流电压中某相的电压值;
步骤S2、构建有效值计算模型,将所述DSP控制器三相交流电压中某相的电压值输入到所述有效值计算模型中,计算得到第一交流电压有效值;
步骤S3、构建低通数字滤波模型,将所述第一交流电压有效值输入到所述低通数字滤波模型中,计算得到第二交流电压有效值;
步骤S2中,所述构建有效值计算模型包括:
将第一次采样得到的单相交流电压表示为x(1),第一次计算得到的交流电压有效值表示为y(1),则第n次采样得到的单相交流电压表示为x(n),第n次计算得到的交流电压有效值表示为y(n);
根据均方根的定义可得:
y(1)2=x(1)2
Figure FDA0003984791580000011
Figure FDA0003984791580000012
Figure FDA0003984791580000013
由均方根的定义公式递推得到:
Figure FDA0003984791580000014
对上式差分方程两边同时进行Z变换,可得:
Figure FDA0003984791580000015
由上式可得,有效值计算模型为:
Figure FDA0003984791580000016
Z为变换算子。
2.如权利要求1所述的适用于DSP控制器的交流电压有效值实时计算方法,其特征在于,步骤S3中,所述低通数字滤波模型为二阶级联巴特沃斯低通滤波模型。
3.如权利要求2所述的适用于DSP控制器的交流电压有效值实时计算方法,其特征在于,步骤S3中,所述构建低通数字滤波模型包括:
二阶级联巴特沃斯低通数字滤波器是由两个一阶巴特沃斯低通数字滤波器组成,两个级联一阶巴特沃斯低通数字滤波器的传递函数分别表示为:
Figure FDA0003984791580000021
Figure FDA0003984791580000022
根据Z变换可得各系数与截止频率和采样频率间的关系表示为:
Figure FDA0003984791580000023
Figure FDA0003984791580000024
Figure FDA0003984791580000025
Figure FDA0003984791580000026
其中,K3、K4、K5、K6分别为传递函数中对应的系数,fs为截止频率,fc为采样频率。
CN202010663956.5A 2020-07-10 2020-07-10 一种适用于dsp控制器的交流电压有效值实时计算方法 Active CN111781415B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010663956.5A CN111781415B (zh) 2020-07-10 2020-07-10 一种适用于dsp控制器的交流电压有效值实时计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010663956.5A CN111781415B (zh) 2020-07-10 2020-07-10 一种适用于dsp控制器的交流电压有效值实时计算方法

Publications (2)

Publication Number Publication Date
CN111781415A CN111781415A (zh) 2020-10-16
CN111781415B true CN111781415B (zh) 2023-02-28

Family

ID=72767378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010663956.5A Active CN111781415B (zh) 2020-07-10 2020-07-10 一种适用于dsp控制器的交流电压有效值实时计算方法

Country Status (1)

Country Link
CN (1) CN111781415B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03264867A (ja) * 1990-03-15 1991-11-26 Hitachi Ltd ディジタル信号処理方法及びその装置
CN1607392A (zh) * 2003-10-14 2005-04-20 珠海炬力集成电路设计有限公司 基于滤波器的电压有效值的测量方法
CN101718814A (zh) * 2009-11-13 2010-06-02 国网电力科学研究院 一种基于fpga实现的安全稳定控制装置电气量测算方法
CN103033675A (zh) * 2012-12-11 2013-04-10 航天科工深圳(集团)有限公司 一种交流同步采集系统及其实现方法
CN108614148A (zh) * 2018-03-30 2018-10-02 国网江西省电力有限公司电力科学研究 一种交流电流、电压同步采样系统及其采样方法
CN110441591A (zh) * 2019-09-17 2019-11-12 贵州电网有限责任公司 一种改进的电子式互感器电流采集方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03264867A (ja) * 1990-03-15 1991-11-26 Hitachi Ltd ディジタル信号処理方法及びその装置
CN1607392A (zh) * 2003-10-14 2005-04-20 珠海炬力集成电路设计有限公司 基于滤波器的电压有效值的测量方法
CN101718814A (zh) * 2009-11-13 2010-06-02 国网电力科学研究院 一种基于fpga实现的安全稳定控制装置电气量测算方法
CN103033675A (zh) * 2012-12-11 2013-04-10 航天科工深圳(集团)有限公司 一种交流同步采集系统及其实现方法
CN108614148A (zh) * 2018-03-30 2018-10-02 国网江西省电力有限公司电力科学研究 一种交流电流、电压同步采样系统及其采样方法
CN110441591A (zh) * 2019-09-17 2019-11-12 贵州电网有限责任公司 一种改进的电子式互感器电流采集方法

Also Published As

Publication number Publication date
CN111781415A (zh) 2020-10-16

Similar Documents

Publication Publication Date Title
CN103257271A (zh) 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法
CN103762595B (zh) 一种抑制电网谐波和幅值扰动的锁频控制方法及系统
CN203287435U (zh) 一种基于stm32f107vct6的微电网谐波与间谐波检测装置
CN103956919A (zh) 电网电压不平衡三电平整流器简化模型预测控制方法
CN112310991B (zh) 基于扰动观测的储能变换器改进有限集模型预测控制方法
CN106980044A (zh) 一种适应风电接入的电力系统谐波电流估计方法
Li Design and realization of FIR digital filters based on MATLAB
CN202075344U (zh) 一种电能表及电能校正系统
CN110068729B (zh) 一种信号相量计算方法
CN116008902A (zh) 一种电能脉冲产生方法及系统
CN111781415B (zh) 一种适用于dsp控制器的交流电压有效值实时计算方法
CN104111373A (zh) 一种用于智能变电站数字化电能的计量方法
CN102928660A (zh) 基于fir数字滤波器的无功功率测量方法
CN112526208A (zh) 一种基于高耦合度迭代模型的电能质量测量系统及方法
CN110649801B (zh) 一种对母线电压的采样方法及pfc控制电路、电源转换电路
CN209787027U (zh) 一种能有效压制纹波电压的前馈采样电路及其开关电源
CN107154802A (zh) 一种校正测量数据的方法及装置
CN115639400A (zh) 一种高精度交流电能计量方法
CN101650385A (zh) 一种电力网络仪表修改采样频率的方法
CN111277163B (zh) 基于fpga和改进电压重构实现有源阻尼的故障辨识控制方法
CN113295917A (zh) 一种基于罗氏线圈与lpct组合的保护装置采样单元
CN103166224B (zh) 一种三相四线制并联有源电力滤波器输出电感优化方法
Wei Research on FCS-MPC for LC filter inverter [J]
CN105391326A (zh) 一种光伏逆变的直流分量的控制方法及装置
CN115001016B (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