CN105159068B - 一种改进型二阶系统的参数设计方法 - Google Patents
一种改进型二阶系统的参数设计方法 Download PDFInfo
- Publication number
- CN105159068B CN105159068B CN201510439427.6A CN201510439427A CN105159068B CN 105159068 B CN105159068 B CN 105159068B CN 201510439427 A CN201510439427 A CN 201510439427A CN 105159068 B CN105159068 B CN 105159068B
- Authority
- CN
- China
- Prior art keywords
- order system
- improved
- value
- formula
- typical
- 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
- 238000013461 design Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000013016 damping Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 8
- 230000014509 gene expression Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种改进型二阶系统的参数设计方法,包括以下步骤,在(0.447,1)范围内选取ξd的取值,在范围内选取ξ的取值,根据给定的上升时间求出ωn,根据ξd、ξ、ωn求出改进型二阶系统模型。本发明克服了传统参数设计方法计算复杂的缺点。
Description
技术领域
本发明属于自适应控制技术领域,特别涉及了一种改进型二阶系统的参数设计方法。
背景技术
近年来,作为自适应控制技术的一个分支,模型参考自适应控制技术,以其良好的控制效果得到了广泛的应用,当被控对象为二阶系统时,通常选择参考模型也为二阶系统,此时,需要设计具有一定动静态性能的二阶系统。除此之外,其他应用场合也会遇到设计二阶系统模型的需要。
目前最常见的标准二阶系统模型的传递函数如式(1)所示,这里称之为典型二阶系统。
其中,ωn为系统的固有频率,ξ为系统的阻尼比。
对于典型二阶系统,在欠阻尼(0<ξ<1)情况下,单位阶跃输入下的超调量σ、上升时间tr和调节时间ts是重要性能指标,可以作为二阶系统参数设计的依据,具体设计步骤为:
1)给定超调量指标数据,按照式(2)确定阻尼比。实际设计中,也可先给定阻尼比,再代入式(2)求解出超调量,看是否满足要求。
2)给定上升时间指标tr,并结合步骤1)中确定的ξ,代入式(3)求得系统的固有频率ωn。
3)将ξ、ωn值代入式(4)计算调节时间,检验调节时间ts是否满足要求,其中,Δ为允许的误差带宽度,可人为给定,例如选取Δ为0.01。
4)将满足条件的ξ、ωn带入式(1)即可得到典型二阶系统模型。
在典型二阶系统的基础上加入微分环节可得图1所示二阶系统模型,这里称之为改进型二阶系统,其传递函数为:
为便于表述,做如下定义:
z=1/Td (6)
其中,称ξd为等效阻尼比,Td为微分系数。与典型二阶系统类似,也可用上升时间tr,峰值时间tp和超调量σ等指标描述改进型二阶系统的性能,相关表达式如下:
其中,
可以看出,上述表达式都为超越方程,而且每个性能指标都与多个参数有关,因此,如果按照传统方法确定式(5)中的系统参数,即首先给定性能指标,然后根据式(8)计算确定二阶系统参数Td、ξ、ωn,计算量将会很大。
改进型二阶系统参数设计存在的问题总结如下:
1)大多采用试凑法,计算方法不够严谨,需要经过反复多次才能完成设计;
2)单位阶跃响应的性能指标表达式复杂,先定指标数据,再求系统参数的话,需要求解超越方程,计算难度大。
发明内容
本发明提供一种改进型二阶系统的参数设计方法,克服传统参数设计方法计算复杂的缺点。
本发明针对上述背景技术中提到的改进型二阶系统,技术方案为:
(1)在(0.447,1)范围内,选择ξd的取值;
(2)确定ξd后,在范围内,选取ξ的取值;
(3)确定ξ后,给定系统的上升时间tr,代入式(10),求得ωn,
(4)将ξd、ξ、ωn代入式(11)得到z,
(5)将ξd、ωn、z代入式(5),得到改进型二阶系统模型。
本发明的一种优选方案,在步骤(1)中,在(0.633,1)范围内选取ξd的取值。
本发明的一种优选方案,在步骤(1)中,选取ξd的取值为0.95,在步骤(2)中,选取ξ的取值为0.2。
本发明的一种优选方案,在步骤(1)中,选取ξd的取值为0.8,在步骤(2)中,选取ξ的取值为0.3。
采用上述技术方案带来的有益效果:
1)结合工程设计的要求,选取合适性能指标,确定合理的参数,保证系统稳定运行,动、静态性能良好;
2)避免求解超越方程,参数设计过程更加简单;
3)步骤明确,便于操作;
附图说明
图1是改进型二阶系统的结构框图;
图2是典型二阶系统与改进型二阶系统单位阶跃响应曲线比较图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
相较于典型二阶系统,改进型二阶系统增加了对输出误差的微分运算,因此,在ξ、ωn保持不变情况下,改进型方案将具有更好的动态性能,这非常适合用于对动态性能要求较高的场合。
首先,为发挥微分环节的优势,并考虑微分环节可能对系统稳定带来的不利影响,在工程设计中要对Td的取值进行适当限制,所以令z=kzξωn,kz∈(1,10),结合式(7)可以得到式(12)所示不等式组:
为了提高效果,在本实施例中,限制kz∈(2,5),从而得到ξd∈(0.633,1)。
式(8)显示ξd不再是超调量σ的唯一决定因素,但对σ仍有重要影响。其中,指数部分主导σ的取值,这意味着,增大ξd,并同时减小ωn,则可以取得相同的超调量。考虑到欠阻尼情况下,ωn是ξ的增函数,因此,在工程设计中,为保持相同的超调量σ,可在增大ξd取值的同时,降低ξ的取值。
根据上述分析,在给定超调量情况下,如果ξd取值较小,则ξ,或ωn的取值将会比较大,结合式(7)可知,此时z也将比较大,再结合式(6)可知微分系数Td将比较小,从而降低响应的快速性。因此,在相同超调情况下,为进一步改善响应的快速性,可尽量提高ξd的取值,从而间接提高Td的取值。当然,过高的Td会降低系统的稳定性,所以ξd可以取为0.95.
另外,分析式(12)可以发现,ξd越大,则ξ的取值范围越大,不论下限值ξmin,还是上限值ξmax都随ξd的增大而增大。
基于上述分析,兼顾设计的可操作性,总结得到发明内容中记载的改进型二阶系统的参数设计方法。
为验证前面的理论分析和参数设计方法的正确性,在matlab/simulink下进行了仿真。仿真主要包括以下两方面内容:对改进型二阶系统参数设计方案相关理论分析的验证,对改进型二阶系统参数设计方法的验证。
首先,对相关结论进行验证
这里按照上述要求对ξd和ξ进行取值,根据相应步骤构建各个不同参数下的改进型二阶系统模型。以单位阶跃信号为输入,观察各个二阶模型的动态响应效果。其中ξd和ξ的取值、各个二阶系统模型表达式以及动态性能指标仿真数据等相关信息如表1所示。
表1
通过对表1中的数据的分析,可以证明前述相关结论的正确性:
1)ξd不再是超调量的唯一决定因素,但对σ仍有重要影响。为避免超调太大,ξd取值不宜太低。
2)适当增大ξd的取值,与此同时降低ξ的取值,可获得同样的超调量;
3)在相同超调情况下,为进一步改善响应的快速性,可尽量提高ξd的取值;
4)ξd越大,则ξ的取值范围越大,不论下限值ξmin,还是上限值ξmax都随ξd的增大而增大
为了达到不同的超调效果或响应速度,可以选取不同的参数。这里选取超调量为5%作为参数选取依据。满足条件的改进型二阶系统模型相关信息如表2所示:
表2
在超调量相同的情况下,希望获得更快的响应速度,所以这里选取ξd=0.95,ξ=0.2所对应的改进型二阶系统模型作为最终的选取目标。
然后,对改进型二阶系统参数设计方法进行验证
以单位阶跃信号为输入,比较相同超调量下典型二阶系统和改进型二阶系统的动态响应效果。
其中改进型二阶系统选取的是ξd=0.95,ξ=0.2,ωn=3.6,z=2.4所对应的二阶系统模型,模型表达式:
典型二阶系统模型选取的是ξ=0.6877,ωn=3.6所对应的二阶系统模型,模型表达式:
仿真结果如图2(横坐标为时间,纵坐标为幅值)所示。从图2的仿真结果中可以看出,在超调量相同的情况下,改进型二阶系统的单位阶跃响应的上升时间和调节时间均比典型二阶系统的小,即改进型二阶系统的动态响应效果更好,尤其体现在快速性方面,证实所提出的改进型二阶系统参数设计方法的正确性。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (6)
1.一种改进型二阶系统的参数设计方法,所述改进型二阶系统是在典型二阶系统的基础上加入微分环节,典型二阶系统的闭环传递函数如式(1)所示,改进型二阶系统的闭环传递函数如式(2)所示,
式(1)、(2)中,ξ为典型二阶系统的阻尼比,ωn为典型二阶系统的固有频率,z为微分系数的倒数,ξd为等效阻尼比,
其特征在于,包括以下步骤:
(1)在(0.447,1)范围内,选取ξd的取值;
(2)在范围内,选取ξ的取值;
(3)给定系统的上升时间tr,根据式(3)求得ωn,
(4)将ξd、ξ、ωn代入式(4)得到z,
(5)将ξd、ωn、z代入式(2)得到改进型二阶系统模型。
2.根据权利要求1所述一种改进型二阶系统的参数设计方法,其特征在于:在步骤(1)中,在(0.633,1)范围内选取ξd的取值。
3.根据权利要求1所述一种改进型二阶系统的参数设计方法,其特征在于:在步骤(1)中,选取ξd的取值为0.95。
4.根据权利要求3所述一种改进型二阶系统的参数设计方法,其特征在于:在步骤(2)中,选取ξ的取值为0.2。
5.根据权利要求1所述一种改进型二阶系统的参数设计方法,其特征在于:在步骤(1)中,选取ξd的取值为0.8。
6.根据权利要求5所述一种改进型二阶系统的参数设计方法,其特征在于:在步骤(2)中,选取ξ的取值为0.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510439427.6A CN105159068B (zh) | 2015-07-23 | 2015-07-23 | 一种改进型二阶系统的参数设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510439427.6A CN105159068B (zh) | 2015-07-23 | 2015-07-23 | 一种改进型二阶系统的参数设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105159068A CN105159068A (zh) | 2015-12-16 |
CN105159068B true CN105159068B (zh) | 2017-08-11 |
Family
ID=54799956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510439427.6A Expired - Fee Related CN105159068B (zh) | 2015-07-23 | 2015-07-23 | 一种改进型二阶系统的参数设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105159068B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101075862A (zh) * | 2007-06-22 | 2007-11-21 | 清华大学 | 无线传感器网络中基于锁相环的时间同步方法 |
CN103853046A (zh) * | 2014-02-14 | 2014-06-11 | 广东工业大学 | 一种压电陶瓷驱动器的自适应学习控制方法 |
CN104242686A (zh) * | 2013-06-18 | 2014-12-24 | 无锡市幅瑞自动化科技有限公司 | 单位功率因数pwm整流器双闭环pi调节器 |
CN104408330A (zh) * | 2014-12-25 | 2015-03-11 | 武汉大学 | 水轮机调节系统高阶数学模型的降阶方法 |
-
2015
- 2015-07-23 CN CN201510439427.6A patent/CN105159068B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101075862A (zh) * | 2007-06-22 | 2007-11-21 | 清华大学 | 无线传感器网络中基于锁相环的时间同步方法 |
CN104242686A (zh) * | 2013-06-18 | 2014-12-24 | 无锡市幅瑞自动化科技有限公司 | 单位功率因数pwm整流器双闭环pi调节器 |
CN103853046A (zh) * | 2014-02-14 | 2014-06-11 | 广东工业大学 | 一种压电陶瓷驱动器的自适应学习控制方法 |
CN104408330A (zh) * | 2014-12-25 | 2015-03-11 | 武汉大学 | 水轮机调节系统高阶数学模型的降阶方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105159068A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109255160B (zh) | 基于神经网络的单元延时预测方法和单元延时灵敏度计算方法 | |
CN113219836B (zh) | 一种分数阶复值忆阻神经网络的投影同步方法及其应用 | |
CN105093931A (zh) | 一种航空发动机非线性系统控制器设计方法 | |
JP2005070940A (ja) | Pidパラメータ調整装置 | |
CN103888104B (zh) | Fir数字滤波器设计方法和系统 | |
CN110119588B (zh) | 基于扩展卡尔曼滤波状态估计值的在线优化设计方法 | |
CN108549227B (zh) | 一种基于分数阶血红细胞模型的时滞反馈控制器设计方法 | |
CN106227964B (zh) | 基于扩张状态观测器的迟滞非线性系统参数辨识方法 | |
CN109188907A (zh) | 一种应用于稳定平台控制系统的遗传退火粒子群混合算法及其稳定平台控制系统 | |
CN107291998B (zh) | 一种基于滑模思想的新型低通滤波器设计方法 | |
CN110649588B (zh) | 一种柔性直流输电控制系统攻击量化评估方法 | |
CN105929865A (zh) | 一种直线伺服系统机械谐振控制方法 | |
CN105159068B (zh) | 一种改进型二阶系统的参数设计方法 | |
CN114841106A (zh) | 一种基于规则引导的遗传算法的集成电路优化方法及系统 | |
CN103792849A (zh) | 可调金属切削系统的鲁棒自适应动态面控制方法 | |
CN108170032A (zh) | 一种提高压电叠堆式驱动器定位精度的方法 | |
Pandey et al. | Multiple models and second level adaptation for a class of nonlinear systems with nonlinear parameterization | |
CN106599527A (zh) | 一种评分变动的处理方法、装置及评分系统 | |
Adamou-Mitiche et al. | Internal projections and L 2 norm for optimal digital filters model reduction: a comparative study | |
CN112346342B (zh) | 一种非仿射动力学系统的单网络自适应评价设计方法 | |
CN115079573A (zh) | 一种非线性系统的高阶扩展强跟踪滤波器 | |
CN107256425B (zh) | 随机权网络泛化能力改进方法及装置 | |
CN109946963B (zh) | 一种判断多回路控制系统裕度的方法 | |
CN106934123A (zh) | 一种基于递归卷积的电路瞬态响应计算方法 | |
CN112750495A (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 | ||
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: 20170811 |