CN105159088B - 一种穿越频率的估算方法 - Google Patents
一种穿越频率的估算方法 Download PDFInfo
- Publication number
- CN105159088B CN105159088B CN201510581362.9A CN201510581362A CN105159088B CN 105159088 B CN105159088 B CN 105159088B CN 201510581362 A CN201510581362 A CN 201510581362A CN 105159088 B CN105159088 B CN 105159088B
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- cross
- over frequency
- control system
- 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
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明提供一种穿越频率的估算方法。它在开环穿越频率的斜率为‑20dB/decade的自动控制系统中,通过注入一个测试点频率f1,再通过检测在该注入频率下的实际速度,就可以估算出系统的穿越频率,如果估算出的穿越频率与测试点的频率差别较大,则通过迭代的方式重新估算来提高精度,整个计算过程只是简单的加减乘除运算,没有指数和对数运算,因此运算过程简单,估算速度快,易于在单片机或DSP等数字处理芯片中编程实现,在实际工程中具有很高的应用价值。
Description
技术领域
本发明涉及自动控制技术领域,具体讲是用于自动控制系统的一种穿越频率的估算方法。
背景技术
穿越频率为描述控制系统中频率特性的一个指标,其定义为幅频穿越0dB处的频率,在穿越频率fc处,幅频特性增益为0dB。穿越频率是自动控制系统中一个重要指标,在自动控制系统设计中,为了保证系统的稳定性,开环系统波特图中幅频曲线在穿越频率附近的斜率一般设计为-20dB/decade。
申请公布号为CN103955143A,名称为“一种交流永磁电机鲁棒控制器的参数整定方法”的中国专利申请公开了一种穿越频率的检测方法,它通过Matlab离线绘制出永磁同步电机开环系统的波特图,再从波特图上求取处控制系统的穿越频率fc,这种借助于绘图工具的求解方法需要测量多个频率下的开环增益,数据量大,且无法实现控制系统穿越频率的智能化估算,很难推广到实际工业应用中。
2005年7月的甘肃联合大学学报,作者为张柳芳的文章“自动控制系统穿越频率的通用解法”提出了自动控制系统穿越频率的通用解法,这些解法都是基于频域系统的波特图坐标系下,求解过程中涉及计算量大的对数、指数及其它运算,计算复杂,计算机编程语言实现困难,也不利于推广。
发明内容
本发明要解决的技术问题是,克服现有技术的缺陷,提供一种运算过程简单,易于实现具有很高的实际工程应用价值的穿越频率的估算方法。
为解决上述技术问题,本发明提供一种穿越频率的估算方法,它用在开环穿越频率的斜率为-20dB/decade的自动控制系统中,包括以下步骤:
步骤1:设置一个速度扰动信号vin,所述信号vin是频率为f1,幅值为A1的正弦信号;
步骤2:把信号vin与自动控制系统的设定速度vset叠加,将叠加后的信号为速度参考信号vref,vref=vset+vin,将速度参考信号vref输入到自动控制系统中;
步骤3:测量实际速度信号v,计算速度偏差信号verr,verr=vref-v;对速度偏差信号verr进行傅里叶反变换求得在f1频率下的信号幅值Ain;
步骤4:对实际速度信号v进行傅里叶反变换求得在f1频率下的信号幅值Aout;
步骤5:利用公式fc=f1Aout/Ain求得初步估算的开环穿越频率fc;
步骤6:比较速度扰动信号的频率f1与步骤5中估算的穿越频率fc,如果|f1-fc|<ε,则跳转到步骤8,否则继续,其中ε为设定的估算误差阈值;
步骤7:将步骤1中的速度扰动信号的频率f1改为步骤5中估算出的频率fc,即取f1=fc,然后跳转到步骤2;
步骤8:估算结束,系统的开环穿越频率为fc。
采用上述方法后,本发明在开环穿越频率的斜率为-20dB/decade的自动控制系统中,通过注入一个测试点频率f1,再通过检测在该注入频率下的实际速度,就可以估算出系统的穿越频率,如果估算出的穿越频率与测试点的频率差别较大,则通过迭代的方式重新估算来提高精度,整个计算过程只是简单的加减乘除运算,没有指数和对数运算,因此运算过程简单,估算速度快,易于在单片机或DSP等数字处理芯片中编程实现,在实际工程中具有很高的应用价值,例如可以直接应用在工业产品量产化过程中的质量检测,利用本方法可以快速精确地估算出系统的穿越频率,从而可以大量地对产品进行检测,筛选出系统性能差异较大的产品,提高产品质量的一致性及可靠性。
所述自动控制系统为伺服系统。在实际的伺服控制系统产品中,开环系统的穿越频率是衡量系统性能的一个重要参数,采用本方法可以快速地估算出系统的穿越频率。
所述步骤7中,设置步骤7的最大运行次数n,若步骤7运行达到最大运行次数n,则直接跳到步骤8。在应用中,如果阈值ε设置的过小,则可能会导致不能收敛,迭代会一直进行下去,因此设置一个最大迭代次数,估算出一个较为近似的穿越频率值。
附图说明
图1本发明所应用的伺服控制系统框图;
图2为斜率为-20dB/decade的开环系统在穿越频率处的波特图;
图3为本发明实施例的Simulink仿真模型图;
图4为实施例的波特图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细的说明:
如图所示,图1为伺服控制系统框图,其中vset是给定的速度指令;vin为速度扰动信号;vref=vset+vin为速度的参考信号;v为伺服电机实际速度;verr=vref-v为速度的偏差;控制器的传递函数为C(s);伺服电机的传递函数为P(s);因此,开环系统的传递函数可以表示为C(s)*P(s)。
本发明通过以下步骤来估算该系统的穿越频率,
步骤1:设置一个速度扰动信号vin,所述信号vin是频率为f1,幅值为A1的正弦信号;即vin=A1sin(2πf1t);频率f1和幅值A1可任意设置,在实际工作中如果是测试系统的穿越频率是否合格,则一般可将频率f1设一个与设定的穿越频率接近的值。
步骤2:把信号vin与自动控制系统的设定速度vset叠加,将叠加后的信号为速度参考信号vref,vref=vset+vin,将速度参考信号vref输入到自动控制系统中;
步骤3:测量实际速度信号v,计算速度偏差信号verr,verr=vref-v;对速度偏差信号verr进行傅里叶反变换求得在f1频率下的信号幅值Ain;对verr进行傅里叶反变换后得到的信号为:其中Ain为信号的幅值,为相位差,本步骤中只用到信号幅值Ain;
步骤4:对实际速度信号v进行傅里叶反变换求得在f1频率下的信号幅值Aout;对v进行傅里叶反变换后得到的信号为:其中Aout为信号的幅值,为相位差;
步骤5:利用公式fc=f1Aout/Ain求得初步估算的开环穿越频率fc;
步骤6:比较速度扰动信号的频率f1与步骤5中估算的穿越频率fc,如果|f1-fc|<ε,则跳转到步骤8,否则继续,其中ε为设定的估算误差阈值;
步骤7:将步骤1中的速度扰动信号的频率f1改为步骤5中估算出的频率fc,即取f1=fc,然后跳转到步骤2;
步骤8:估算结束,系统的开环穿越频率为fc。
在实际的控制系统设计中,为了保证系统的稳定性,开环系统的波特图幅频曲线一般设计成-20dB/decade的斜率穿越0dB,换句话说,开环系统的波特图幅频曲线在穿越频率点附近可以表示为一个线性方程,如图2所示。该线性方程的曲线的斜率为-20dB/decade,因此只需要知道该线性方程的一个点的横坐标和纵坐标,就可以计算出该曲线在纵坐标为0dB时的频率,即系统的穿越频率。由于系统的波特图幅频曲线越靠近穿越频率点就越近似于线性,因此这个点的选取离穿越频率点越近,则穿越频率的计算就越精确。因此本方法增加了迭代来提高估算精度。
所述步骤7中,设置步骤7的最大运行次数n,若步骤7运行达到最大运行次数n,则直接跳到步骤8。在应用中,如果阈值ε设置的过小,则可能会导致不能收敛,迭代会一直进行下去,因此设置一个最大迭代次数,估算出一个较为近似的穿越频率值。
下面举出一实施例来验证本方法估算穿越频率的精确度:
引用伺服电机数学模型设计PI控制器数学模型给定速度指令为vset=500rad/s,选取穿越频率的估算误差阈值为ε=1rad/s。此处可以不设置步骤7的最大运行次数。
为验证本申请提出的估算方法,我们将利用Matlab中的simulink仿真工具,建立该实施例的仿真模型如图3所示:
步骤1:设定速度扰动信号vin,幅值A1=60,频率为ω1=5rad/s,即:vin=60sin(5t);其中ω为角频率,它与频率f的关系是:ω=2πf,由于在制图软件Matlab中常用角频率,因此此处的频率采用角频率进行表示,但不管是采用角频率还是频率,均应在本申请保护范围之内。
步骤2:测量速度偏差信号verr,通过傅里叶反变换求得在ω1频率下的信号幅值为Ain=21.77;
步骤3:测量实际速度信号v,通过傅里叶反变换求得在ω1频率下的信号幅值为Aout=41.6;
步骤4:利用公式求得估算穿越频率为:ωc=ω1Aout/Ain=5*41.6/21.77=9.55rad/s;
步骤5:比较|ω1-ωc|=4.55>ε(其中ε=1rad/s为设定的估算误差阈值);
步骤6:将速度扰动信号vin的频率改为9.55rad/s,设ω2=ωc=9.55rad/s,即:vin=60sin(9.55t);
步骤7:测量速度偏差信号verr,通过傅里叶反变换求得在ω2频率下的信号幅值为Ain=49.34;
步骤8:测量实际速度信号v,通过傅里叶反变换求得在ω2频率下的信号幅值为Aout=46.35;
步骤9:利用公式求得估算穿越频率为:ωc=ω2Aout/Ain=9.55*46.35/49.34=8.97rad/s;
步骤10:比较|ω2-ωc|=0.58<ε(其中ε=1rad/s为设定的估算误差阈值);
步骤11:估算结束,穿越频率的估算值为ωc=8.97rad/s。
为验证估算穿越频率的准确性,我们通过Matlab中的Bode()和Margin()函数对开环传递函数绘制波特图,如图4所示。从图中可以得出在幅值为0dB时,频率为9.1rad/sec。
从图中可以得到开环传递函数的实际穿越频率为9.1rad/sec,与估算值8.97rad/sec很接近。
从上述实施例中可以看出本发明的计算过程仅仅是简单的加减乘除运算,经过一次迭代后就得到与实际穿越频率很接近的估算值,如果估算误差阈值设置的更小,那么估算精度会更高,这种方法计算量小,适用于工业化中的穿越频率的快速检测,应用前景广阔。
Claims (3)
1.一种穿越频率的估算方法,其特征在于:它用在开环穿越频率的斜率为-20dB/decade的自动控制系统中,包括以下步骤:
步骤1:设置一个速度扰动信号vin,所述信号vin是频率为f1,幅值为A1的正弦信号;
步骤2:把信号vin与自动控制系统的设定速度vset叠加,将叠加后的信号为速度参考信号vref,vref=vset+vin,将速度参考信号vref输入到自动控制系统中;
步骤3:测量实际速度信号v,测量速度偏差信号verr,verr=vref-v;对速度偏差信号verr进行傅里叶反变换,求得在f1频率下的信号幅值Ain;
步骤4:对实际速度信号v进行傅里叶反变换,求得在f1频率下的信号幅值Aout;
步骤5:利用公式fc=f1Aout/Ain求得初步估算的开环穿越频率fc;
步骤6:比较速度扰动信号的频率f1与步骤5中估算的穿越频率fc,如果|f1-fc|<ε,则跳转到步骤8,否则继续,其中ε为设定的估算误差阈值;
步骤7:将步骤1中的速度扰动信号的频率f1改为步骤5中估算出的频率fc,即取f1=fc,然后跳转到步骤2;
步骤8:估算结束,系统的开环穿越频率为fc。
2.根据权利要求1所述的一种穿越频率的估算方法,其特征在于:所述自动控制系统为伺服系统。
3.根据权利要求1所述的一种穿越频率的估算方法,其特征在于:所述步骤7中,设置步骤7的最大运行次数n,若步骤7运行达到最大运行次数n,则直接跳到步骤8。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510581362.9A CN105159088B (zh) | 2015-09-14 | 2015-09-14 | 一种穿越频率的估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510581362.9A CN105159088B (zh) | 2015-09-14 | 2015-09-14 | 一种穿越频率的估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105159088A CN105159088A (zh) | 2015-12-16 |
CN105159088B true CN105159088B (zh) | 2017-11-24 |
Family
ID=54799976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510581362.9A Active CN105159088B (zh) | 2015-09-14 | 2015-09-14 | 一种穿越频率的估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105159088B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955143A (zh) * | 2014-05-13 | 2014-07-30 | 长春理工大学 | 一种交流永磁电机鲁棒控制器的参数整定方法 |
CN104635496A (zh) * | 2015-02-10 | 2015-05-20 | 华中科技大学 | 一种带串联校正的pi控制器的参数整定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001259824A1 (en) * | 2000-05-05 | 2001-11-20 | Paul Reed Smith Guitars, Limited Partnership | Analysis of and noise reduction in a complex waveform |
-
2015
- 2015-09-14 CN CN201510581362.9A patent/CN105159088B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955143A (zh) * | 2014-05-13 | 2014-07-30 | 长春理工大学 | 一种交流永磁电机鲁棒控制器的参数整定方法 |
CN104635496A (zh) * | 2015-02-10 | 2015-05-20 | 华中科技大学 | 一种带串联校正的pi控制器的参数整定方法 |
Non-Patent Citations (4)
Title |
---|
Robust controllability of interval fractional order linear time invariant systems;YangQuan Chen等;《Signal Processing 》;20060310(第63期);第2794-2802页 * |
相角裕度和幅值裕度的一种估算方法;刘美侠等;《电器电子教学学报》;20120630;第34卷(第3期);第116-117页 * |
穿越频率ωc的求解;孙慧峰等;《水利电力机械》;20031031;第25卷(第5期);第55,58页 * |
自动控制系统穿越频率的求解;张柳芳;《平顶山师专学报(自然科学)》;19980831;第13卷(第4期);第20-23页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105159088A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105404296A (zh) | 一种两轮自平衡智能车姿态控制方法 | |
CN106849799A (zh) | 一种无轴承异步电机转子磁链观测与转速辨识的方法 | |
CN101360983A (zh) | 温度推定方法及装置 | |
CN108255786A (zh) | 一种称重结果的干扰补偿计算方法及系统 | |
CN103900574A (zh) | 一种基于迭代容积卡尔曼滤波姿态估计方法 | |
CN108762083A (zh) | 一种基于加速度观测器的自动控制系统 | |
CN102455660A (zh) | 基于数字h∞pid控制器的的连续时滞系统控制方法 | |
CN106405210A (zh) | 一种采样数据异常检测方法及系统 | |
CN106939840A (zh) | 用于确定内燃机中的气体质量流量的方法和装置 | |
CN105159088B (zh) | 一种穿越频率的估算方法 | |
CN104734824B (zh) | 强噪声背景下混沌同步检测微弱周期信号的方法 | |
CN104122440A (zh) | 晶振频率漂移感知方法 | |
CN104330644B (zh) | 检测电力系统中正弦波信号幅值的方法和系统 | |
TWI481183B (zh) | 馬達參數量測方法與馬達參數量測系統 | |
CN106446405B (zh) | 一种集成电路器件神经网络建模样本选择方法及装置 | |
CN109858170A (zh) | 一种强非线性Modelica系统模型的频率特性估算系统及方法 | |
CN109510543A (zh) | 一种伺服电机负载惯量的测定方法 | |
CN114614714A (zh) | 一种无速度传感器感应电机高速域的稳定控制方法 | |
CN103701393A (zh) | 一种异步电机弱磁时转矩精度的补偿方法 | |
CN104485867B (zh) | 一种基于转差频率限制的异步电机控制方法 | |
Kychkin et al. | Avtomatizirovannaya informatsionnaya sistema polunaturnogo modelirovaniya staticheskoy nagruzki elektroprivodov | |
CN105958968B (zh) | 有源电力滤波器电气参数在线辨识系统及方法 | |
Gerhard et al. | Constructive nonlinear dynamics–Foundations and application to robust nonlinear control | |
CN109684599A (zh) | 适用于电传飞控直升机的飞行实验数据快速频域辨识方法 | |
Wang et al. | Grid Line Galloping Analysis Using Euler Angle Attitude Estimation Method |
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 |