CN102426416B - 控制对象频域幅值特性辩识方法和装置 - Google Patents
控制对象频域幅值特性辩识方法和装置 Download PDFInfo
- Publication number
- CN102426416B CN102426416B CN201110260155.5A CN201110260155A CN102426416B CN 102426416 B CN102426416 B CN 102426416B CN 201110260155 A CN201110260155 A CN 201110260155A CN 102426416 B CN102426416 B CN 102426416B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- amplitude
- frequency signal
- control 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 213
- 230000008569 process Effects 0.000 claims abstract description 183
- 238000001914 filtration Methods 0.000 claims abstract description 143
- 238000001228 spectrum Methods 0.000 claims abstract description 36
- 230000001105 regulatory effect Effects 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims description 19
- 230000005540 biological transmission Effects 0.000 claims description 18
- 230000000630 rising effect Effects 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 13
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 230000005428 wave function Effects 0.000 claims description 2
- 238000012850 discrimination method Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 14
- 230000003595 spectral effect Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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是本发明控制对象频域幅值特性辩识装置的频率滤波模块的实施效果示意图,选择参数下表所示:
图5是本发明获取控制对象频域传递幅值系数的幅值比较方法的流程图,与图1相比,图5是步骤S104的具体细化步骤,包括:
S401:将参考频率信号、过程频率信号分别进行延时运算,获取延时信号。
对参考频率信号进行Tn/4纯延时运算,作用是使参考频率信号相位滞后π/2,在后面流程中产生一个与参考频率信号正峰值时刻相对应的脉冲信号,得到参考频率信号的Tn/4纯延时信号,表达式如下:
则 。
对过程频率信号进行Tn/4纯延时运算,作用是使过程频率信号相位滞后π/2,在后面流程中产生一个与过程频率信号正峰值时刻相对应的脉冲信号,得到过程频率信号的Tn/4纯延时信号,表达式如下:
设:过程频率信号为
S402:比较延时信号的幅值与零值,分别获取参考频率信号的二进制逻辑信号、过程频率信号的二进制逻辑信号。
将上述参考频率信号的Tn/4纯延时信号幅值进行0值比较,得到参考频率信号的Tn/4纯延时信号的二进制逻辑信号,表达式如下:
将上述过程频率信号的Tn/4纯延时信号幅值进行0值比较,得到过程频率信号的Tn/4纯延时信号的二进制逻辑信号,表达式如下:
其中,B2为1的时间段是:
S403:将二进制逻辑信号的上升沿转换为脉冲输出信号,分别获取参考频率信号的参考脉冲输出信号、过程频率信号的过程脉冲输出信号。
将上述参考频率信号的二进制逻辑信号的上升沿转换为参考脉冲输出信号,脉冲信号宽度较小,将脉冲信号为1看成为一个时刻,该时刻与参考频率信号正峰值出现时刻相对应,表达式如下:
将上述过程频率信号的二进制逻辑信号的上升沿转换为过程脉冲输出信号,脉冲信号宽度较小,将脉冲信号为1看成为一个时刻,该时刻与过程频率信号正峰值出现时刻相对应,表达式如下:
S404:以参考脉冲输出信号对该参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号。
用上述参考脉冲输出信号,对参考频率信号正峰值进行采样保持,具体地,在参考脉冲输出信号为1时刻、对参考频率信号正峰值进行采样,表达式如下:
参考脉冲输出信号为1时刻为:,将其代入下式为:
在参考脉冲输出信号1消失后、对采样信号进行保持,得到参考频率信号的幅值信号。表达式如下:
S405:以过程脉冲输出信号对该过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号。
用上述过程脉冲输出信号,对过频率信号正峰值进行采样保持,具体地,在过程脉冲输出信号为1时刻、对过程频率信号正峰值进行采样,表达式如下:
在过程脉冲输出信号1消失后、对采样信号进行保持,得到过程频率信号的幅值信号,表达式如下:
S406:将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
将上述过程频率信号的幅值信号除以参考频率信号的幅值信号,获得控制对象在所述频率传递幅值系数值,表达式如下:
这里采用的是一种同步峰值检波技术,直接得到输入频率信号的峰值输出,如果输入峰值稳定,输出也稳定,即输出信号中没有波动成分。现有的一些检波技术得到的是输入信号峰值的平均值或接近峰值输出,而且输入峰值稳定,输出信号中也存在波动成分,即不能完全消除波动成分。
需要说明的是,对上述基波成分的参考频率信号、各奇次谐波成分的参考频率信号、基波成分的过程频率信号、各奇次谐波成分的过程频率信号分别进行上述步骤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是本发明控制对象频域幅值特性辩识装置的实施效果之一示意图。
图11是本发明控制对象频域幅值特性辩识装置的实施效果之二示意图。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的硬件平台的方式来实现,当然也可以全部通过硬件来实施。基于这样的理解,本发明的技术方案对背景技术做出贡献的全部或者部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
以上所述的本发明实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。
Claims (8)
1.一种控制对象频域幅值特性辩识方法,其特征在于,包括:
在由调节器和控制对象组成的调节回路的输入信号中加入频谱辨识信号;
对所述调节器的输出信号进行滤波,获取参考频率信号;
对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数;
比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数的步骤,包括:
将参考频率信号、过程频率信号分别进行延时运算,获取延时信号;
比较延时信号的幅值与零值,分别获取参考频率信号的二进制逻辑信号、过程频率信号的二进制逻辑信号;
将二进制逻辑信号的上升沿转换为脉冲输出信号,分别获取参考频率信号的参考脉冲输出信号、过程频率信号的过程脉冲输出信号;
以参考脉冲输出信号对该参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号;
以过程脉冲输出信号对该过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号;
将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
2.根据权利要求1所述的控制对象频域幅值特性辩识方法,其特征在于:
所述频谱辨识信号包括由基波成分和预设数次谐波成分组成的方波辨识信号。
3.根据权利要求2所述的控制对象频域幅值特性辩识方法,其特征在于,对所述调节器的输出信号进行滤波,获取参考频率信号的步骤,包括:
将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号;
将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号。
4.根据权利要求3所述的控制对象频域幅值特性辩识方法,其特征在于,对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号的步骤,包括:
将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
将所述控制对象的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的过程频率信号。
5.一种控制对象频域幅值特性辩识装置,其特征在于,包括:
与调节器和控制对象组成的调节回路输入端相连的辨识信号发生模块,用于在调节回路的输入信号中加入频谱辨识信号;
与所述调节器输出端相连的参考频率滤波模块,用于对所述调节器的输出信号进行滤波,获取参考频率信号;
与控制对象输出端相连的过程频率滤波模块,用于对与所述调节器输出端相连的控制对象的输出信号进行滤波,获取过程频率信号;
与所述参考频率滤波模块、所述过程频率滤波模块分别相连的幅值比较模块,用于比较所述参考频率信号的幅值与所述过程频率信号的幅值,获取控制对象频域传递幅值系数;
所述幅值比较模块包括过程幅值检测模块,参考幅值检测模块,与所述过程幅值检测模块、所述参考幅值检测模块分别相连的幅值比较单元;
所述参考幅值检测模块包括参考延时单元,用于将参考频率信号进行延时运算,获取参考频率信号的延时信号;与所述参考延时单元相连的参考比较单元,用于比较参考频率信号的延时信号的幅值与零值,获取参考频率信号的二进制逻辑信号;与所述参考比较单元相连的参考上升沿脉冲单元,用于通过参考频率信号的二进制逻辑信号获取参考频率信号的参考脉冲输出信号;与所述参考上升沿脉冲单元相连的参考数据采集保持单元,用于以参考频率信号的参考脉冲输出信号对该参考数据采集保持单元接入的参考频率信号的正峰值进行采样保持,获取该参考频率信号的幅值信号;
所述过程幅值检测模块包括过程延时单元,用于将过程频率信号进行延时运算,获取过程频率信号的延时信号;与所述过程延时单元相连的过程比较单元,用于比较过程频率信号的延时信号的幅值与零值,获取过程频率信号的二进制逻辑信号;与所述过程比较单元相连的过程上升沿脉冲单元,用于通过过程频率信号的二进制逻辑信号获取过程频率信号的过程脉冲输出信号;与所述过程上升沿脉冲单元相连的过程数据采集保持单元,用于以过程频率信号的过程脉冲输出信号对该过程数据采集保持单元接入的过程频率信号的正峰值进行采样保持,获取该过程频率信号的幅值信号;
所述幅值比较单元分别与所述过程幅值检测模块的过程数据采集保持单元、所述参考幅值检测模块的参考数据采集保持单元相连,用于将过程频率信号的幅值信号除以与所述频率相对应的参考频率信号的幅值信号,获得控制对象在该频率的传递幅值系数,其中,所述该频率为基波频率或各个谐波频率。
6.根据权利要求5所述的控制对象频域幅值特性辩识装置,其特征在于,所述辨识信号发生模块包括:
方波函数单元,用于输出方波辨识信号,所述方波辨识信号包括基波成分和各个谐波成分。
7.根据权利要求6所述的控制对象频域幅值特性辩识装置,其特征在于,所述参考频率滤波模块包括:
参考基波滤波单元,用于将所述调节器的输出信号中的基波成分进行带通滤波,获取基波的参考频率信号;
至少一个参考谐波滤波单元,用于将所述调节器的输出信号中的各个谐波成分进行带通滤波,获取各个谐波的参考频率信号。
8.根据权利要求7所述的控制对象频域幅值特性辩识装置,其特征在于,所述过程频率滤波模块包括:
过程基波滤波模块,用于将所述控制对象的输出信号中的基波成分进行带通滤波,获取基波的过程频率信号;
至少一个过程谐波滤波单元,用于将所述控制对象的输出信号中的谐波成分进行带通滤波,获取各个谐波的过程频率信号。
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 CN102426416A (zh) | 2012-04-25 |
CN102426416B true 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) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707142B (zh) * | 2012-05-23 | 2015-02-11 | 广东电网公司电力科学研究院 | 控制对象的频谱分析方法和装置 |
CN103217291B (zh) * | 2013-01-06 | 2016-05-04 | 国电联合动力技术有限公司 | 一种风力发电机组故障诊断方法及*** |
CN103713558B (zh) * | 2013-12-31 | 2016-06-15 | 西安龙腾新能源科技发展有限公司 | 一种用于数字控制电路的环路增益扫描方法 |
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 | 宁波罗杰克智能科技有限公司 | 一种幅频特性曲线的绘制方法 |
CN109489747B (zh) * | 2018-12-24 | 2020-03-20 | 合肥工业大学 | 一种基于谐波分析的电磁流量计信号处理方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI107093B (fi) * | 1997-09-24 | 2001-05-31 | Nokia Networks Oy | Integroidun oskillaattorin automaattinen virittäminen |
JP4091047B2 (ja) * | 2002-10-31 | 2008-05-28 | 深▲川▼市中▲興▼通▲訊▼股▲分▼有限公司 | 広帯域プリディストーション線形化の方法およびシステム |
US7190213B2 (en) * | 2005-03-25 | 2007-03-13 | Freescale Semiconductor, Inc. | Digital time constant tracking technique and apparatus |
CN100596014C (zh) * | 2005-06-03 | 2010-03-24 | 孙力 | 步进电动机正弦波驱动器 |
KR101363081B1 (ko) * | 2006-06-19 | 2014-02-21 | 제이디에스 유니페이즈 코포레이션 | 주파수 도메인 반사계측기를 사용한 네트워크 기기 검출 |
CN101697084B (zh) * | 2009-10-19 | 2011-12-28 | 大连海事大学 | 一种基于rls滤波器的电液伺服***随机振动控制方法 |
-
2011
- 2011-09-05 CN CN201110260155.5A patent/CN102426416B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102426416A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102426416B (zh) | 控制对象频域幅值特性辩识方法和装置 | |
CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN109831184B (zh) | 一种带通滤波方法、***、存储介质及终端设备 | |
CN112946363A (zh) | 一种基于神经网络的电网阻抗在线辨识方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN104217112B (zh) | 一种基于多类型信号的电力***低频振荡分析方法 | |
CN103163372A (zh) | 采用hht的电力***谐波分析方法 | |
CN109142863B (zh) | 一种电力***测频方法及*** | |
CN101858938A (zh) | 基于自适应滤波原理的瞬时频率测量方法 | |
CN109828154B (zh) | 一种分频段复合正交脉冲注入的三相电网阻抗测量方法 | |
Li | Design and realization of FIR digital filters based on MATLAB | |
CN102411301B (zh) | 控制对象频域相位特性辩识方法和装置 | |
CN117169590B (zh) | 一种基于软件变采样率的电力谐波分析的方法和装置 | |
CN112255457B (zh) | 适用于自动准同期装置的相角差测量方法 | |
CN116885718A (zh) | 电网杂质电流远程治理方法及*** | |
WO2024109031A1 (zh) | 一种基于hht变换的电压质量扰动检测方法 | |
CN106053940B (zh) | 一种基于方波傅里叶级数分解的谐波分析方法 | |
Borisov et al. | A computationally efficient RDFT-based reference signal generator for active compensators | |
Fusco et al. | A study on short-term sea profile prediction for wave energy applications | |
Hu | Design and verification of FIR filter based on Matlab and DSP | |
CN106324342A (zh) | 一种基于查表的谐波检测方法 | |
CN112362968B (zh) | 一种基于预调制cdsc与sdft的单相谐波实时提取方法 | |
CN104849548A (zh) | 一种电力***瞬时频率监测方法及*** | |
CN112485522B (zh) | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 | |
CN112986680A (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 | ||
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. |