CN102426416A - 控制对象频域幅值特性辩识方法和装置 - Google Patents

控制对象频域幅值特性辩识方法和装置 Download PDF

Info

Publication number
CN102426416A
CN102426416A CN2011102601555A CN201110260155A CN102426416A CN 102426416 A CN102426416 A CN 102426416A CN 2011102601555 A CN2011102601555 A CN 2011102601555A CN 201110260155 A CN201110260155 A CN 201110260155A CN 102426416 A CN102426416 A CN 102426416A
Authority
CN
China
Prior art keywords
signal
frequency
amplitude
frequency signal
controlling object
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.)
Granted
Application number
CN2011102601555A
Other languages
English (en)
Other versions
CN102426416B (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.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid 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 Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201110260155.5A priority Critical patent/CN102426416B/zh
Publication of CN102426416A publication Critical patent/CN102426416A/zh
Application granted granted Critical
Publication of CN102426416B publication Critical patent/CN102426416B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了控制对象频域幅值特性辩识方法和装置。该方法包括:在由调节器和控制对象组成的调节回路的输入信号中加入频谱辨识信号;对所述调节器的输出信号进行滤波,获取参考频率信号;对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。采用本发明,可以节约频谱资源地辩识控制对象的频域幅值特性,不需要使用到热工对象全部连续的频谱资源,只需使用到有限数量非连续的单频率点频谱,并通过使用滤波器滤出这些单频率点频谱信号,排除绝大部分的干扰信号。

Description

控制对象频域幅值特性辩识方法和装置
技术领域
本发明涉及热工控制技术领域,尤其涉及控制对象频域幅值特性辩识方法和装置。
背景技术
火电厂热工控制对象具有大惯性、纯延时、以及时变性等特点,特别是控制对象的特性难以深入准确了解等因素,使得系统难以控制。如何有效克服各种不利影响,从而提高热工控制回路的调节品质,一直是热工控制领域不断追求的目标之一。为提高热工对象的控制品质,对控制对象模型辩识是非常重要的。
现有的对象模型辩识技术主要是基于时域方面对过程对象模型的特性进行辩识,主要存在以下三大问题:其一,时域辨识的实现算法复杂,因此,难以在现有DCS(Distributed Control System,分布式控制系统)中应用。其二,时域辩识的精度受过程干扰影响很大;热工对象是需要使用频谱资源的,其频谱的基本特点是:以零频率为极点,从有限高的高频到无限低的低频连续的频谱。时域模型辩识技术最根本的问题在于:需要使用到热工对象全部连续的频谱资源,在实验室的仿真过程中,可以做到这些频谱信号无任何噪声污染,能够获得良好的辨识效果。但在实际的应用过程中,工业过程信号并不干净,对时域模型辩识影响很大,甚至使辨识的效果很差。其三是激励信号问题。为了获得较高的辩识精度,需要有较高幅度的充分激励信号,这对过程工艺正常生产影响较大。综合这些问题,使得现有的模型辩识技术难以在实践中得到广泛应用。
发明内容
为了解决上述技术问题,本发明提供了控制对象频域幅值特性辩识方法和装置,能够非常节约频谱资源地实现控制对象的频域幅值特性的辨识。
本发明提供了控制对象频域幅值特性辩识方法,包括:
在由调节器和控制对象组成的调节回路的输入信号中加入频谱辨识信号;
对所述调节器的输出信号进行滤波,获取参考频率信号;
对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
相应地,本发明还提供了控制对象频域特性辩识装置,包括:
与调节器和控制对象组成的调节回路输入端相连的辨识信号发生模块,用于在调节回路的输入信号中加入频谱辨识信号;
与所述调节器输出端相连的参考频率滤波模块,用于对所述调节器的输出信号进行滤波,获取参考频率信号;
与控制对象输出端相连的过程频率滤波模块,用于对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
与所述参考频率滤波模块、所述过程频率滤波模块分别相连的幅值比较模块,用于比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
实施本发明,具有如下有益效果:
本发明通过在由PID调节器和控制对象构成的闭环控制系统中,向过程给定信号加入方波频谱辩识信号。这样,在PID调节器的输出信号和过程输出信号中、均包含与频谱辩识信号相同的频谱成分,仅是频谱的幅值特性相对频谱辩识信号发生了变化,这种变化与PID调节器特性和控制对象特性有着必然的内在联系。通过对PID调节器输出信号中频谱成分与过程输出信号中频谱成分的幅值比较,从而辩识控制对象的频域幅值特性。其重要的用途之一是用于对象模型传递函数的计算。
与时域辩识方法不同,频域辩识方法非常节约频谱资源,不需要使用到热工对象全部连续的频谱资源,只使用到有限数量非连续的单频率点频谱,数量上仅占到热工对象全部连续频谱资源的极小部分。可以使用带通滤波器滤出这些单频率点频谱信号,将绝大部分的干扰信号排除在外。
附图说明
图1是本发明控制对象频域幅值特性辩识方法的流程图;
图2是本发明获取参考频率信号的滤波方法的流程图;
图3是本发明获取过程频率信号的滤波方法的流程图;
图4是本发明控制对象频域幅值特性辩识装置的频率滤波模块的实施效果示意图;
图5是本发明获取控制对象频域传递幅值系数的幅值比较方法的流程图;
图6是本发明控制对象频域幅值特性辩识装置的示意图;
图7是本发明控制对象频域幅值特性辩识装置的参考频率滤波模块的示意图;
图8是本发明控制对象频域幅值特性辩识装置的过程频率滤波模块的示意图;
图9是本发明控制对象频域幅值特性辩识装置的幅值比较模块的结构示意图;
图10是本发明控制对象频域幅值特性辩识装置的实施效果之一示意图;
图11是本发明控制对象频域幅值特性辩识装置的实施效果之二示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
图1是本发明控制对象频域幅值特性辩识方法的流程图,包括:
S101:在由调节器和控制对象组成的调节回路的输入信号中加入频谱辨识信号;
S102:对所述调节器的输出信号进行滤波,获取参考频率信号,所述参考频率信号包含与所述频谱辨识信号相同的频谱成分;
S103:对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号,所述过程频率信号包含与所述频谱辨识信号相同的频谱成分;
S104:比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
在由调节器和控制对象组成的调节回路中,调节回路的输入端接入的是过程给定信号,向所述过程给定信号中加入频谱辩识信号。所述调节器的输出端与所述控制对象的输入端相连,故此,在调节器的输出信号和控制对象的输出信号当中,均包含与所述频谱辩识信号相同的频谱成分,仅是频谱的幅值特性相对频谱辩识信号发生了变化,这种变化与调节器特性和控制对象特性有着必然的内在联系。
通过将所述调节器的输出信号和所述控制对象的输出信号进行滤波,可以将绝大部分的干扰信号排除在外。与现有技术的时域辨识方法相比,本发明的抗干扰能力更强,能够更准确地对控制对象进行辨识。
经过上述滤波,排除绝大部分干扰之后,将所述参考频率信号与所述过程频率信号进行幅值比较,获取控制对象频域传递幅值系数,从而辨识出控制对象的频域特性。与现有技术的时域辩识方法不同,本发明采用的频域辩识方法非常节约频谱资源,不需要使用到热工对象全部连续的频谱资源,只使用到有限数量非连续的单频率点频谱,数量上仅占到热工对象全部连续频谱资源的极小部分。
另外,本发明还具有实现算法简单,容易在现有DCS中应用,对准确认识了解控制对象的频域幅值特性具有重要意义。另外,其重要的用途还包括用于对象模型传递函数的计算。
所述频谱辨识信号包括由基波成分和预设数次谐波成分组成的方波辨识信号。
典型的频谱辨识信号包括方波信号、脉冲激发信号、正弦波信号、余弦波信号等等。优选地,本发明采用方波辨识信号,其目的在于方波频谱辩识信号除了上升沿和下降沿有变化外、信号的中间过程无变化,另外所加方波频谱辩识信号的幅值较小,对过程工艺影响较小。
理论上,由于方波的偶次频率为低电平,而奇次频率为高电平,所以,所述方波辩识信号成分包括基波成分和无限次的奇次谐波成分。但是,本发明不排除可以通过电路转换技术,将方波的偶次频率转换为高电平,而将奇次频率转为为低电平,一切基于本发明思想所作的类似修改或等同变换,应包含在本发明的保护范围之内。
在基于方波的奇次频率为高电平的前提下,本发明所述的预设数次谐波成分,具体为预设数奇次谐波成分。由于高次数谐波的幅值会迅速下降,高阶对象对高次数谐波会迅速衰减,又由于A/D转换器位数限制和量化噪声限制,辨识高次数获取的频率信号的信噪比并不能得到保证。故此,为了满足工程实际需要,优选地,本发明采取由基波成分和第3次至第21次的奇次谐波成分组成的方波辨识信号。
图2是本发明获取参考频率信号的滤波方法的流程图,与图1相比,图2是步骤S102的具体细化步骤,包括:
S201:将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号。
将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号,具体包括,步骤S202至步骤S211。
S202:将所述调节器的输出信号中的3次谐波成分进行带通滤波,获取3次谐波的参考频率信号。
S203:将所述调节器的输出信号中的5次谐波成分进行带通滤波,获取5次谐波的参考频率信号。
S204:将所述调节器的输出信号中的7次谐波成分进行带通滤波,获取7次谐波的参考频率信号。
S205:将所述调节器的输出信号中的9次谐波成分进行带通滤波,获取9次谐波的参考频率信号。
S206:将所述调节器的输出信号中的11次谐波成分进行带通滤波,获取11次谐波的参考频率信号。
S207:将所述调节器的输出信号中的13次谐波成分进行带通滤波,获取13次谐波的参考频率信号。
S208:将所述调节器的输出信号中的15次谐波成分进行带通滤波,获取15次谐波的参考频率信号。
S209:将所述调节器的输出信号中的17次谐波成分进行带通滤波,获取17次谐波的参考频率信号。
S210:将所述调节器的输出信号中的19次谐波成分进行带通滤波,获取19次谐波的参考频率信号。
S211:将所述调节器的输出信号中的21次谐波成分进行带通滤波,获取21次谐波的参考频率信号。
图3是本发明获取过程频率信号的滤波方法的流程图,与图1相比,图3是步骤S103的具体细化步骤,包括:
S301:将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
将所述控制对象的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的过程频率信号,具体包括,步骤S302至步骤S311。
S302:将所述控制对象的输出信号中的3次谐波成分进行带通滤波,获取3次谐波的过程频率信号。
S303:将所述控制对象的输出信号中的5次谐波成分进行带通滤波,获取5次谐波的过程频率信号。
S304:将所述控制对象的输出信号中的7次谐波成分进行带通滤波,获取7次谐波的过程频率信号。
S305:将所述控制对象的输出信号中的9次谐波成分进行带通滤波,获取9次谐波的过程频率信号。
S306:将所述控制对象的输出信号中的11次谐波成分进行带通滤波,获取11次谐波的过程频率信号。
S307:将所述控制对象的输出信号中的13次谐波成分进行带通滤波,获取13次谐波的过程频率信号。
S308:将所述控制对象的输出信号中的15次谐波成分进行带通滤波,获取15次谐波的过程频率信号。
S309:将所述控制对象的输出信号中的17次谐波成分进行带通滤波,获取17次谐波的过程频率信号。
S310:将所述控制对象的输出信号中的19次谐波成分进行带通滤波,获取19次谐波的过程频率信号。
S311:将所述控制对象的输出信号中的21次谐波成分进行带通滤波,获取21次谐波的过程频率信号。
需要说明的是,如图2、图3所述的滤波方法采用的是带通滤波,优选地,可采用LCR滤波器。对上述参考频率信号或过程频率信号进行滤波。由于干扰信号也是由各种频率成分构成的,而“带通滤波”最大特点是可将大部分不需要的频率成分排出在外。频域辨识最大优势也在于强抗干扰性,只取有用信号,将不需要的信号排除在外。不排除干扰信号频率与辨识频率相同或非常接近问题,这主要取决于信噪比的大小,如果信噪比足够高,也完全可以获得良好的性能,还因为辨识点数较多,计算传递函数3到4个点已足够了。另外,设计滤波器对相邻谐波或基波频率的衰减率大于40分贝,带宽0.003rad/s,这样获得的频率信号比较干净,有较好的抗干扰性和响应时间,从加入方波开始,需要2-3个方波周期,滤波器输出幅值才会稳定。但频域辨识取的是稳态结果,响应时间不是最重要的。也可设计滤波器对相邻谐波或基波频率的衰减率大于80分贝,带宽0.0005rad/s,可获得更好的性能,但需要6-7个方波周期,滤波器输出幅值才会稳定。故此,优选地,可采用LCR滤波器进行四级滤波,可以达到上述性能指标。
图4是本发明控制对象频域幅值特性辩识装置的频率滤波模块的实施效果示意图,选择参数下表所示:
Figure 685791DEST_PATH_IMAGE001
 图5是本发明获取控制对象频域传递幅值系数的幅值比较方法的流程图,与图1相比,图5是步骤S104的具体细化步骤,包括:
S401:将参考频率信号、过程频率信号分别进行延时运算,获取延时信号。
对参考频率信号进行Tn/4纯延时运算,作用是使参考频率信号相位滞后π/2,在后面流程中产生一个与参考频率信号正峰值时刻相对应的脉冲信号,得到参考频率信号的Tn/4纯延时信号,表达式如下:
设参考频率信号为
Figure 606473DEST_PATH_IMAGE002
其中,
Figure 1683DEST_PATH_IMAGE003
为参考频率信号的幅值,
Figure 939025DEST_PATH_IMAGE004
为频率信号的角频率,频率信号周期为:
Figure 27067DEST_PATH_IMAGE005
则 
Figure 107150DEST_PATH_IMAGE006
对过程频率信号进行Tn/4纯延时运算,作用是使过程频率信号相位滞后π/2,在后面流程中产生一个与过程频率信号正峰值时刻相对应的脉冲信号,得到过程频率信号的Tn/4纯延时信号,表达式如下:
设:过程频率信号为  
其中,
Figure 83251DEST_PATH_IMAGE008
为控制对象在所述频率的传递幅值系数,
Figure 404511DEST_PATH_IMAGE009
为控制对象在所述频率的滞后角弧度值,
Figure 706310DEST_PATH_IMAGE003
为参考频率信号的幅值,频率信号周期为:
Figure 708901DEST_PATH_IMAGE010
则:
Figure 294603DEST_PATH_IMAGE011
S402:比较延时信号的幅值与零值,分别获取参考频率信号的二进制逻辑信号、过程频率信号的二进制逻辑信号。
将上述参考频率信号的Tn/4纯延时信号幅值进行
Figure 773383DEST_PATH_IMAGE012
0值比较,得到参考频率信号的Tn/4纯延时信号的二进制逻辑信号,表达式如下:
Figure 77325DEST_PATH_IMAGE013
其中,A2为1的时间段是
Figure 431077DEST_PATH_IMAGE014
,(i为非负整数)。
将上述过程频率信号的Tn/4纯延时信号幅值进行
Figure 808969DEST_PATH_IMAGE012
0值比较,得到过程频率信号的Tn/4纯延时信号的二进制逻辑信号,表达式如下:
Figure 472031DEST_PATH_IMAGE015
其中,B2为1的时间段是:
Figure 246958DEST_PATH_IMAGE016
,(i为非负整数)。
S403:将二进制逻辑信号的上升沿转换为脉冲输出信号,分别获取参考频率信号的参考脉冲输出信号、过程频率信号的过程脉冲输出信号。
将上述参考频率信号的二进制逻辑信号的上升沿转换为参考脉冲输出信号,脉冲信号宽度较小,将脉冲信号为1看成为一个时刻,该时刻与参考频率信号正峰值出现时刻相对应,表达式如下:
Figure 653669DEST_PATH_IMAGE017
其中,
Figure 964695DEST_PATH_IMAGE018
为单位脉冲函数,i为非负整数。
将上述过程频率信号的二进制逻辑信号的上升沿转换为过程脉冲输出信号,脉冲信号宽度较小,将脉冲信号为1看成为一个时刻,该时刻与过程频率信号正峰值出现时刻相对应,表达式如下:
Figure 470763DEST_PATH_IMAGE019
上式中
Figure 749298DEST_PATH_IMAGE018
为单位脉冲函数,i为非负整数。
S404:以参考脉冲输出信号对该参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号。
用上述参考脉冲输出信号,对参考频率信号正峰值进行采样保持,具体地,在参考脉冲输出信号为1时刻、对参考频率信号正峰值进行采样,表达式如下:
参考脉冲输出信号为1时刻为:
Figure 946317DEST_PATH_IMAGE020
,将其代入下式为:
Figure 298801DEST_PATH_IMAGE021
则采样时刻:
Figure 38087DEST_PATH_IMAGE022
在参考脉冲输出信号1消失后、对采样信号进行保持,得到参考频率信号的幅值信号。表达式如下:
S405:以过程脉冲输出信号对该过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号。
用上述过程脉冲输出信号,对过频率信号正峰值进行采样保持,具体地,在过程脉冲输出信号为1时刻、对过程频率信号正峰值进行采样,表达式如下:
过程脉冲输出信号为1时刻为: 
Figure 975267DEST_PATH_IMAGE024
,将其代入下式为:
Figure 244574DEST_PATH_IMAGE025
则采样时刻:
Figure 669608DEST_PATH_IMAGE026
在过程脉冲输出信号1消失后、对采样信号进行保持,得到过程频率信号的幅值信号,表达式如下:
Figure 657156DEST_PATH_IMAGE027
S406:将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
将上述过程频率信号的幅值信号除以参考频率信号的幅值信号,获得控制对象在所述频率传递幅值系数值,表达式如下:
Figure 756830DEST_PATH_IMAGE028
这里采用的是一种同步峰值检波技术,直接得到输入频率信号的峰值输出,如果输入峰值稳定,输出也稳定,即输出信号中没有波动成分。现有的一些检波技术得到的是输入信号峰值的平均值或接近峰值输出,而且输入峰值稳定,输出信号中也存在波动成分,即不能完全消除波动成分。
需要说明的是,对上述基波成分的参考频率信号、各奇次谐波成分的参考频率信号、基波成分的过程频率信号、各奇次谐波成分的过程频率信号分别进行上述步骤S401至步骤S405,获得一组基波成分的参考/过程频率信号的幅值、一组3次谐波成分的参考/过程频率信号的幅值…一组21次谐波成分的参考/过程频率信号的幅值,最后,将各组幅值进行上述步骤S406,得到一组控制对象的传递幅值系数,从而获得基波频率到21次奇次谐波频率范围内控制对象的频域幅值特性。对准确认识了解控制对象的特性具有重要意义,更重要的用途是用于对象模型传递函数的计算。
图6是本发明控制对象频域幅值特性辩识装置的示意图,包括:
典型的,由调节器和控制对象组成的调节回路;
与调节器和控制对象组成的调节回路输入端相连的辨识信号发生模块,用于在调节回路的输入信号中加入频谱辨识信号;
与所述调节器输出端相连的参考频率滤波模块,用于对所述调节器的输出信号进行滤波,获取参考频率信号;
与控制对象输出端相连的过程频率滤波模块,用于对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
与所述参考频率滤波模块、所述过程频率滤波模块分别相连的幅值比较模块,用于比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
如图6所示的辨识信号发生模块包括:
方波函数单元,用于输出方波辨识信号,所述方波辨识信号包括基波成分和各个谐波成分。
图7是本发明控制对象频域幅值特性辩识装置的参考频率滤波模块的示意图,所述参考频率滤波模块包括:
参考基波滤波单元,用于将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号;
至少一个参考谐波滤波单元,用于将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号。具体地,包括,参考3次谐波滤波单元,用于对3次谐波成分的参考频率信号进行滤波;参考5次谐波滤波单元,用于对5次谐波成分的参考频率信号进行滤波;参考7次谐波滤波单元,用于对7次谐波成分的参考频率信号进行滤波;参考9次谐波滤波单元,用于对9次谐波成分的参考频率信号进行滤波;参考11次谐波滤波单元,用于对11次谐波成分的参考频率信号进行滤波;参考13次谐波滤波单元,用于对13次谐波成分的参考频率信号进行滤波;参考15次谐波滤波单元,用于对15次谐波成分的参考频率信号进行滤波;参考17次谐波滤波单元,用于对17次谐波成分的参考频率信号进行滤波;参考19次谐波滤波单元,用于对19次谐波成分的参考频率信号进行滤波;参考21次谐波滤波单元,用于对21次谐波成分的参考频率信号进行滤波。
图8是本发明控制对象频域幅值特性辩识装置的过程频率滤波模块的示意图,所述过程频率滤波模块包括:
过程基波滤波单元,用于将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
至少一个过程谐波滤波单元,用于将所述控制对象的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的过程频率信号。具体地,包括,过程3次谐波滤波单元,用于对3次谐波成分的过程频率信号进行滤波;过程5次谐波滤波单元,用于对5次谐波成分的过程频率信号进行滤波;过程7次谐波滤波单元,用于对7次谐波成分的过程频率信号进行滤波;过程9次谐波滤波单元,用于对9次谐波成分的过程频率信号进行滤波;过程11次谐波滤波单元,用于对11次谐波成分的过程频率信号进行滤波;过程13次谐波滤波单元,用于对13次谐波成分的过程频率信号进行滤波;过程15次谐波滤波单元,用于对15次谐波成分的过程频率信号进行滤波;过程17次谐波滤波单元,用于对17次谐波成分的过程频率信号进行滤波;过程19次谐波滤波单元,用于对19次谐波成分的过程频率信号进行滤波;过程21次谐波滤波单元,用于对21次谐波成分的过程频率信号进行滤波。
图9是本发明控制对象频域幅值特性辩识装置的幅值比较模块的结构示意图。如图9所示,所述幅值比较模块包括过程幅值检测模块,参考幅值检测模块,与所述过程幅值检测模块、所述参考幅值检测模块分别相连的幅值比较单元C1;
所述参考幅值检测模块包括参考延时单元B1,用于将参考频率信号进行延时运算,获取参考频率信号的延时信号;与所述参考延时单元B1相连的参考比较单元B2,用于比较参考频率信号的延时信号的幅值与零值,获取参考频率信号的二进制逻辑信号;与所述参考比较单元B2相连的参考上升沿脉冲单元B3,用于通过参考频率信号的二进制逻辑信号获取参考频率信号的参考脉冲输出信号;与所述参考上升沿脉冲单元B3相连的参考数据采集保持单元B4,用于以参考频率信号的参考脉冲输出信号对该参考数据采集保持单元接入的参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号;
所述过程幅值检测模块包括过程延时单元A1,用于将过程频率信号进行延时运算,获取过程频率信号的延时信号;与所述过程延时单元A1相连的过程比较单元A2,用于比较过程频率信号的延时信号的幅值与零值,获取过程频率信号的二进制逻辑信号;与所述过程比较单元A2相连的过程上升沿脉冲单元A3,用于通过过程频率信号的二进制逻辑信号获取过程频率信号的过程脉冲输出信号;与所述过程上升沿脉冲单元A3相连的过程数据采集保持单元A4,用于以过程频率信号的过程脉冲输出信号对该过程数据采集保持单元接入的过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号;
所述幅值比较单元C1分别与所述过程幅值检测模块的过程数据采集保持单元A4、所述参考幅值检测模块的参考数据采集保持单元B4相连,用于将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。优选地,所述比较单元包括除法器,用于完成上述除法运算。
需要说明的是,各个滤波和幅值比较需要互相对应。例如,经过图7的参考基波滤波单元滤波后获得的参考基波频率信号,与经过图8的参考过程基波滤波单元滤波后获得的过程基波频率信号相对应,将基波频率的过程频率信号的幅值信号除以与所述基波频率相对应的参考频率信号的幅值信号,获得控制对象在所述基波频率的传递幅值系数。其余,各奇次谐波频率的参考频率信号和各奇次谐波频率的过程频率信号如此类推,最后,获得一组控制对象频域传递幅值系数,从而获得基波频率到21次奇次谐波频率范围内控制对象的频域幅值特性。对准确认识了解控制对象的特性具有重要意义,更重要的用途是用于对象模型传递函数的计算。
图10是本发明控制对象频域幅值特性辩识装置的实施效果之一示意图。
如图10所示,为对传递函数为
Figure 398421DEST_PATH_IMAGE029
的二阶惯性惯性环节幅值频率特性仿真辩识结果。图中实线为理论幅值频率特性曲线,粗点为仿真辩识结果,可见辩识结果与理论特性完全符合,具有很高的辩识精度。
图11是本发明控制对象频域幅值特性辩识装置的实施效果之二示意图。
如图11所示,为对传递函数为的二阶惯性惯性环节在有干扰环境下的辩识精度仿真结果,具体在参考滤波模块组输入中加入带宽为100Hz,平均幅值为方波辨识信号幅值100%的高频高斯干扰信号。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的硬件平台的方式来实现,当然也可以全部通过硬件来实施。基于这样的理解,本发明的技术方案对背景技术做出贡献的全部或者部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
以上所述的本发明实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。

Claims (10)

1.一种控制对象频域幅值特性辩识方法,其特征在于,包括:
在由调节器和控制对象组成的调节回路的输入信号中加入频谱辨识信号;
对所述调节器的输出信号进行滤波,获取参考频率信号;
对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
2.根据权利要求1所述的控制对象频域幅值特性辩识方法,其特征在于:
所述频谱辨识信号包括由基波成分和预设数次谐波成分组成的方波辨识信号。
3.根据权利要求2所述的控制对象频域幅值特性辩识方法,其特征在于,对所述调节器的输出信号进行滤波,获取参考频率信号的步骤,包括:
将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号;
将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号。
4.根据权利要求3所述的控制对象频域幅值特性辩识方法,其特征在于,对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号的步骤,包括:
将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
将所述控制对象的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的过程频率信号。
5.根据权利要求4所述的控制对象频域幅值特性辩识方法,其特征在于,比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数的步骤,包括:
将参考频率信号、过程频率信号分别进行延时运算,获取延时信号;
比较延时信号的幅值与零值,分别获取参考频率信号的二进制逻辑信号、过程频率信号的二进制逻辑信号;
将二进制逻辑信号的上升沿转换为脉冲输出信号,分别获取参考频率信号的参考脉冲输出信号、过程频率信号的过程脉冲输出信号;
以参考脉冲输出信号对该参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号;
以过程脉冲输出信号对该过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号;
将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
6.一种控制对象频域幅值特性辩识装置,其特征在于,包括:
与调节器和控制对象组成的调节回路输入端相连的辨识信号发生模块,用于在调节回路的输入信号中加入频谱辨识信号;
与所述调节器输出端相连的参考频率滤波模块,用于对所述调节器的输出信号进行滤波,获取参考频率信号;
与控制对象输出端相连的过程频率滤波模块,用于对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
与所述参考频率滤波模块、所述过程频率滤波模块分别相连的幅值比较模块,用于比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数。
7.根据权利要求6所述的控制对象频域幅值特性辩识装置,其特征在于,所述辨识信号发生模块包括:
方波函数单元,用于输出方波辨识信号,所述方波辨识信号包括基波成分和各个谐波成分。
8.根据权利要求7所述的控制对象频域幅值特性辩识装置,其特征在于,所述参考频率滤波模块包括:
参考基波滤波单元,用于将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号;
至少一个参考谐波滤波单元,用于将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号。
9.根据权利要求8所述的控制对象频域幅值特性辩识装置,其特征在于,所述过程频率滤波模块包括:
过程基波滤波模块,用于将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
至少一个过程谐波滤波单元,用于将所述控制对象的输出信号中的谐波成分进行带通滤波,获取各个谐波的过程频率信号。
10.根据权利要求9所述的控制对象频域幅值特性辩识装置,其特征在于:所述幅值比较模块包括过程幅值检测模块,参考幅值检测模块,与所述过程幅值检测模块、所述参考幅值检测模块分别相连的幅值比较单元;
所述参考幅值检测模块包括参考延时单元,用于将参考频率信号进行延时运算,获取参考频率信号的延时信号;与所述参考延时单元相连的参考比较单元,用于比较参考频率信号的延时信号的幅值与零值,获取参考频率信号的二进制逻辑信号;与所述参考比较单元相连的参考上升沿脉冲单元,用于通过参考频率信号的二进制逻辑信号获取参考频率信号的参考脉冲输出信号;与所述参考上升沿脉冲单元相连的参考数据采集保持单元,用于以参考频率信号的参考脉冲输出信号对该参考数据采集保持单元接入的参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号;
所述过程幅值检测模块包括过程延时单元,用于将过程频率信号进行延时运算,获取过程频率信号的延时信号;与所述过程延时单元相连的过程比较单元,用于比较过程频率信号的延时信号的幅值与零值,获取过程频率信号的二进制逻辑信号;与所述过程比较单元相连的过程上升沿脉冲单元,用于通过过程频率信号的二进制逻辑信号获取过程频率信号的过程脉冲输出信号;与所述过程上升沿脉冲单元相连的过程数据采集保持单元,用于以过程频率信号的过程脉冲输出信号对该过程数据采集保持单元接入的过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号;
所述幅值比较单元分别与所述过程幅值检测模块的过程数据采集保持单元、所述参考幅值检测模块的参考数据采集保持单元相连,用于将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
CN201110260155.5A 2011-09-05 2011-09-05 控制对象频域幅值特性辩识方法和装置 Active CN102426416B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110260155.5A CN102426416B (zh) 2011-09-05 2011-09-05 控制对象频域幅值特性辩识方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110260155.5A CN102426416B (zh) 2011-09-05 2011-09-05 控制对象频域幅值特性辩识方法和装置

Publications (2)

Publication Number Publication Date
CN102426416A true CN102426416A (zh) 2012-04-25
CN102426416B CN102426416B (zh) 2014-04-09

Family

ID=45960416

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110260155.5A Active CN102426416B (zh) 2011-09-05 2011-09-05 控制对象频域幅值特性辩识方法和装置

Country Status (1)

Country Link
CN (1) CN102426416B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707142A (zh) * 2012-05-23 2012-10-03 广东电网公司电力科学研究院 控制对象的频谱分析方法和装置
CN103217291A (zh) * 2013-01-06 2013-07-24 国电联合动力技术有限公司 一种风力发电机组故障诊断方法及系统
CN103713558A (zh) * 2013-12-31 2014-04-09 西安龙腾新能源科技发展有限公司 一种用于数字控制电路的环路增益扫描方法
CN104267259A (zh) * 2014-09-30 2015-01-07 广东电网有限责任公司电力科学研究院 电力谐波测量方法和系统
CN105182800A (zh) * 2015-10-09 2015-12-23 宁波罗杰克智能科技有限公司 一种幅频特性曲线的绘制方法
CN104267257B (zh) * 2014-09-30 2017-02-01 广东电网有限责任公司电力科学研究院 基于点频滤波器的信号幅频谱检测方法和系统
CN104502706B (zh) * 2014-12-29 2017-05-03 广东电网有限责任公司电力科学研究院 电力信号的谐波幅值测量方法和系统
CN109489747A (zh) * 2018-12-24 2019-03-19 合肥工业大学 一种基于谐波分析的电磁流量计信号处理方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1271476A (zh) * 1997-09-24 2000-10-25 诺基亚网络有限公司 集成振荡器的自动调谐
CN1694345A (zh) * 2005-06-03 2005-11-09 孙力 步进电动机正弦波驱动器
US20060240786A1 (en) * 2002-10-31 2006-10-26 Xiaowei Liu Method and system for broadband predistortion linearization
CN101124721A (zh) * 2005-03-25 2008-02-13 飞思卡尔半导体公司 数字时间常数跟踪技术和装置
CN101329383A (zh) * 2006-06-19 2008-12-24 安科特纳有限责任公司 采用频域反射计的网络设备检测
CN101697084A (zh) * 2009-10-19 2010-04-21 大连海事大学 一种基于rls滤波器的电液伺服系统随机振动控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1271476A (zh) * 1997-09-24 2000-10-25 诺基亚网络有限公司 集成振荡器的自动调谐
US20060240786A1 (en) * 2002-10-31 2006-10-26 Xiaowei Liu Method and system for broadband predistortion linearization
CN101124721A (zh) * 2005-03-25 2008-02-13 飞思卡尔半导体公司 数字时间常数跟踪技术和装置
CN1694345A (zh) * 2005-06-03 2005-11-09 孙力 步进电动机正弦波驱动器
CN101329383A (zh) * 2006-06-19 2008-12-24 安科特纳有限责任公司 采用频域反射计的网络设备检测
CN101697084A (zh) * 2009-10-19 2010-04-21 大连海事大学 一种基于rls滤波器的电液伺服系统随机振动控制方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707142A (zh) * 2012-05-23 2012-10-03 广东电网公司电力科学研究院 控制对象的频谱分析方法和装置
CN102707142B (zh) * 2012-05-23 2015-02-11 广东电网公司电力科学研究院 控制对象的频谱分析方法和装置
CN103217291B (zh) * 2013-01-06 2016-05-04 国电联合动力技术有限公司 一种风力发电机组故障诊断方法及系统
CN103217291A (zh) * 2013-01-06 2013-07-24 国电联合动力技术有限公司 一种风力发电机组故障诊断方法及系统
CN103713558A (zh) * 2013-12-31 2014-04-09 西安龙腾新能源科技发展有限公司 一种用于数字控制电路的环路增益扫描方法
CN103713558B (zh) * 2013-12-31 2016-06-15 西安龙腾新能源科技发展有限公司 一种用于数字控制电路的环路增益扫描方法
CN104267259A (zh) * 2014-09-30 2015-01-07 广东电网有限责任公司电力科学研究院 电力谐波测量方法和系统
CN104267259B (zh) * 2014-09-30 2017-01-18 广东电网有限责任公司电力科学研究院 电力谐波测量方法和系统
CN104267257B (zh) * 2014-09-30 2017-02-01 广东电网有限责任公司电力科学研究院 基于点频滤波器的信号幅频谱检测方法和系统
CN104502706B (zh) * 2014-12-29 2017-05-03 广东电网有限责任公司电力科学研究院 电力信号的谐波幅值测量方法和系统
CN105182800A (zh) * 2015-10-09 2015-12-23 宁波罗杰克智能科技有限公司 一种幅频特性曲线的绘制方法
CN109489747A (zh) * 2018-12-24 2019-03-19 合肥工业大学 一种基于谐波分析的电磁流量计信号处理方法
CN109489747B (zh) * 2018-12-24 2020-03-20 合肥工业大学 一种基于谐波分析的电磁流量计信号处理方法

Also Published As

Publication number Publication date
CN102426416B (zh) 2014-04-09

Similar Documents

Publication Publication Date Title
CN102426416A (zh) 控制对象频域幅值特性辩识方法和装置
CN104979845B (zh) 基于多级线性最优理论的多频段高压直流输电附加阻尼控制方法
CN105207236A (zh) 基于svg的抑制次同步振荡自适应控制系统
CN105137757A (zh) 一种具有频率自适应能力的重复控制器及控制方法
CN107786201A (zh) 一种基于锁频环的二阶广义积分器结构及锁相环同步方法
CN109768565B (zh) 适用于柔性直流的无源阻抗适配器参数设计方法及装置
CN101887238B (zh) 一种特定次重复控制器及控制方法
CN105515021B (zh) 多模式附加次/超同步振荡控制方法和控制系统
CN108493961B (zh) 可再生能源制氢系统抑制次/超同步振荡的方法及系统
CN102411301B (zh) 控制对象频域相位特性辩识方法和装置
CN106300386B (zh) 基于svg动态抑制电网次同步振荡的频率闭环控制方法
CN110365018B (zh) 一种自适应宽频带正负序分离方法
CN102420559A (zh) 基于系统辨识和遗传算法的发电机广域阻尼控制方法
CN107607784A (zh) 一种开环同步方法
CN103199533A (zh) 基于两点式工况的无源滤波器优化设计方法
CN104111373A (zh) 一种用于智能变电站数字化电能的计量方法
CN206788245U (zh) 一种基于专用数字滤波器的谐波电能计量模块
CN109088617A (zh) 比率可变数字重采样滤波器
CN112054549A (zh) 一种直流输电系统稳定性提升的方法
CN103401447A (zh) 一种级联型pwm整流器的调制方法
CN103956736A (zh) 主动式电力滤波装置多环谐波控制方法及系统
CN100446377C (zh) 一种晶闸管控制电抗器的控制器
CN110336279A (zh) 基于并网变换器的电力系统振荡自适应抑制方法、系统及介质
CN210867493U (zh) 一种舰载用电设备的多级谐波滤波器
CN109599867A (zh) 基于双环结构的三相并网变流器电流环干扰抑制控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Patentee after: Electric Power Research Institute of Guangdong Power Grid Co.,Ltd.

Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Patentee before: ELECTRIC POWER RESEARCH INSTITUTE OF GUANGDONG POWER GRID Corp.

CP01 Change in the name or title of a patent holder