CN104049541A - 一种直流电机鲁棒控制器的参数整定方法 - Google Patents
一种直流电机鲁棒控制器的参数整定方法 Download PDFInfo
- Publication number
- CN104049541A CN104049541A CN201410201321.8A CN201410201321A CN104049541A CN 104049541 A CN104049541 A CN 104049541A CN 201410201321 A CN201410201321 A CN 201410201321A CN 104049541 A CN104049541 A CN 104049541A
- Authority
- CN
- China
- Prior art keywords
- omega
- direct current
- lambda
- coefficient
- matlab
- 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
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明涉及一种直流电机鲁棒控制器的参数整定方法,属于分数阶自动控制技术领域。应用于位置随动***中直流电机,其传递函数为其待整定FO[PI]控制器传递函数形式利用MATLAB画出被控对象P(s)的伯德图,求得在频率ωc处的模值m和相角n,利用C(s)P(s)=G(s),利用MATLAB求解关于积分阶次λ的方程,利用稳定性条件,求得比例系数和微分系数利用所得λ以及得到Kp和Ki代入本发明有益效果是减少了控制器参数整定的计算量,简化了分数阶FO[PI]鲁棒控制器参数整定过程。
Description
技术领域
本发明属于分数阶自动控制技术领域,主要涉及一种基于MATLAB分数阶FO[PI]鲁棒控制器参数整定方法。
背景技术
工业现代化发展水平是衡量一个国家综合国力水平的重要因素,电机是这些工业设备的动力来源,是设备正常运行的保障,这就使得对电机控制的研究就显得尤为迫切,开发具有高位置精度、响应速度快、高可靠性的伺服控制器已成为研究热点。
MATLAB是矩阵实验室的简称,用于算法开发、数据可视化、数据分析以及控制仿真等,尤其是近年来,MATLAB在控制***仿真、分析和设计方面得到了广泛应用。用MATLAB语言编程效率高,程序调试十分方便,可大大缩减软件开发周期。
随着分数阶控制理论的发展,证实了分数阶控制器具有比传统整数阶控制器更好的响应能力和抗干扰能力,可以使控制***获得更好的动态性能和鲁棒性。由于直流电机伺服控制***具有控制简单、稳定性高、调速范围宽等优点,在实际***中得到了广泛应用。近年来,许多学者将分数阶FO[PI]控制器的做为直流电动机的伺服控制器,以获得更好的动态性能和鲁棒性。
直流电机伺服***的性能不仅和所选用的控制器的结构有关,而且还取决于伺服控制器的参数。然而,因为分数阶FO[PI]控制器多了一个可调参数λ,使得FO[PI]控制器参数整定过程复杂,运算量大。对不同的被控对象,参数整定方程需要重新推导和计算,使得控制器设计变得繁琐而又费时。
直流电机以其优良的性能被广泛地应用于位置随动***中。其传递函数为 若忽略电枢电感La及粘性阻尼系数B,传递函数可以近似为其中Tm表示电动机机电时间常数,Ke分别表示反电动势系数,Kt表示电磁力矩系数,J是电动机的转动惯量,La是电枢电感,B是粘性摩擦系数,s为拉普拉斯算子;运行实例中,选定位置随动***中常用直流电 动机作为控制对象,不失一般性,其数学模型传递函数可表示为其中,T为时间常数。
发明内容
本发明提供一种直流电机鲁棒控制器的参数整定方法,以解决控制器设计存在的繁琐和费时的问题。
应用于位置随动***中直流电机,其传递函数为 若忽略电枢电感La及粘性阻尼系数B,传递函数可以近似为其中Tm表示电动机机电时间常数,Ke分别表示反电动势系数,Kt表示电磁力矩系数,J是电动机的转动惯量,La是电枢电感,B是粘性摩擦系数,s为拉普拉斯算子;运行实例中,选定位置随动***中常用直流电动机作为控制对象,不失一般性,其数学模型传递函数可表示为其中,T为时间常数;
所述直流电机鲁棒控制器参数整定方法,包括以下步骤:
(1)对于直流电机被控对象的数学模型传递函数P(s),其待整定FO[PI]控制器传递函数形式待整定参数为比例系数Kp,积分系数Ki和积分阶次λ,并给定需校正穿越频率ωc和需保持稳定的相位裕度φm;
(2)利用MATLAB画出被控对象P(s)的伯德图,求得在频率ωc处的模值m和相角n,同时可以求得被控对象频率ωc在相位变化率
(3)利用C(s)P(s)=G(s),鲁棒性条件:
得到:
注意到其中
利用MATLAB求解关于积分阶次λ的方程(3);
(4)利用稳定性条件:在开环***穿越频率ωc处相位裕度为φm;
C(jωc)P(jωc)=1∠φm-180° (4)
根据步骤(2)所求得的被控对象在频率ωc处的模值m和相角n,得到
由步骤(3)中
令A=10-m/20,从而得到
可以求得比例系数 和微分系数
(5)利用步骤(3)中MATLAB所得λ以及步骤(4)得到Kp和Ki代入 即完成了分数阶FO[PI]鲁棒控制器参数整定。
本发明的有益效果:减少了控制器参数整定的计算量,简化了分数阶FO[PI]鲁棒控制器参数整定过程。由于本发明中参数方程仅与A、θ和有关,针对不同的电机被控对象,本发明参数方程不需要重新推导和计算。其中A、θ和可以由步骤2输入被控对象传递函数利用MATLAB函数指令求得,从而本发明可以快速进行 分数阶FO[PI]鲁棒控制器参数整定。除此之外,本发明采用的基于MATLAB分数阶FO[PI]鲁棒控制器参数整定方法求得参数唯一且有效。
附图说明
图1是基于MATLAB分数阶FO[PI]鲁棒控制器参数整定方法的流程图;
图2是具体实施例1被控对象伯德图在频率10rad/s的模值和相角;
图3是具体实施例1基于MATLAB求解的FO[PI]控制器参数λ;
图4是具体实施例1中的所设计的开环***伯德图;
图5是具体实施例1中整个闭环控制***的阶跃响应图;其中,三条曲线是开环***增益分别为0.9、1和1.1是的阶跃响应曲线。
图6是具体实施例2被控对象伯德图在频率30rad/s的模值和相角;
图7是具体实施例2基于MATLAB求解的控制器参数λ;
图8是具体实施例2中的所设计的开环***伯德图;
图9是具体实施例2中整个闭环控制***的阶跃响应图;其中,三条曲线是开环***增益分别为0.9、1和1.1是的阶跃响应曲线。
图10是具体实施例3被控对象伯德图在频率10rad/s的模值和相角;
图11是具体实施例3基于MATLAB求解的控制器参数λ;
图12是具体实施例3中的所设计的开环***伯德图;
图13是具体实施例3中整个闭环控制***的阶跃响应图;其中,三条曲线是开环***增益分别为0.9、1和1.1是的阶跃响应曲线。
具体实施方式
下面结合具体实施例以及附图对本发明做进一步详细说明。在此,本发明的具体实施例及其说明用于解释本发明,但并不作为对本发明的限定。
实施例1
1.假设直流电机被控对象***的数学模型传递函数其中T=0.4。并给定穿越频率ωc=10rad/s和需保持稳定的相位裕度φm=70°。
2.求得被控对象在频率10rad/s处的模值为-12.3dB和相角-76°,同时可以求得被控对象频率10rad/s在相位变化率
3.可以求得θ=φm-n-180°=-34°得到方程:
用MATLAB图解法求得分数阶积分阶次λ=0.5564。
4.可以求得A=10-m/20=4.1210,从而求得比例系数和积分系数
5.所求分数阶FO[PI]鲁棒控制器为
图4为所设计的开环***的伯德图;其中,从图中可以看出***在穿越频率ωc附近的相角裕度保持恒定。
图5为所设计的控制***的阶跃响应图其中,三条曲线在开环***增益分别为0.9、1和1.1的情况下***也能保持稳定的输出超调量,即利用本发明所列方法整定出FO[PI]结构的分数阶控制器具有非常好的鲁棒特性。
通过实施例1,可知基于MATLAB分数阶FO[PD]鲁棒控制器参数整定方法,针对不同的被控对象可以快速进行参数整定。
实施例2
1.直流电机被控对象***的数学模型传递函数其中时间常数变成为T=0.1。并给定穿越频率ωc=30rad/s和需保持稳定的相位裕度φm=70°。
2.求得被控对象在频率30rad/s处的模值为-10dB和相角-71.5°,同时可以求得被控对象频率30rad/s在相位变化率
3.可以求得θ=φm-n-180°=-38.5°得到方程:
用MATLAB图解法求得分数阶积分阶次λ=0.6657。
4.可以求得A=10-m/20=3.1623,从而求得比例系数和积分系数
5.所求分数阶FO[PI]鲁棒控制器为
图8为所设计的开环***的伯德图;其中,从图中可以看出***在穿越频率ωc附近的相角裕度保持恒定。
图9为所设计的控制***的阶跃响应图其中,三条曲线在开环***增益分别为0.9、1和1.1的情况下***也能保持稳定的输出超调量,即利用本发明所列方法整定出FO[PI]结构的分数阶控制器具有非常好的鲁棒特性。
通过实施例2,可知基于MATLAB分数阶FO[PD]鲁棒控制器参数整定方法,针对不同的被控对象可以快速进行参数整定。
实施例3
1.在实际***中往往存在延迟,假设直流电机被控对象的数学模型传递函数 其中T=0.4,L=0.01。并给定穿越频率ωc=10rad/s和需保持稳定的相位裕度φm=50°。
2.求得被控对象在频率10rad/s处的模值为-12.3dB和相角-81.8°,同时可以求得被控对象频率10rad/s在相位变化率
3.可以求得θ=φm-n-180°=-48.2°得到方程:
用MATLAB图解法求得分数阶积分阶次λ=0.7904。
4.可以求得A=10-m/20=4.1210,从而求得比例系数和积分系数
5.分数阶FO[PI]鲁棒控制器为
图12为所设计的开环***的伯德图;其中,从图中可以看出***在穿越频率ωc附近的相角裕度保持恒定。
图13为所设计的控制***的阶跃响应图其中,三条曲线在开环***增益分别为0.9、1和1.1的情况下***也能保持稳定的输出超调量,即利用本发明所列方法整定出FO[PI]结构的分数阶控制器具有非常好的鲁棒特性。
通过实施例3,可知基于MATLAB分数阶FO[PD]鲁棒控制器参数整定方法,针对不同的被控对象可以快速进行参数整定。
Claims (1)
1.一种直流电机鲁棒控制器的参数整定方法,
应用于位置随动***中直流电机,其传递函数为 若忽略电枢电感La及粘性阻尼系数B,传递函数可以近似为其中Tm表示电动机机电时间常数,Ke分别表示反电动势系数,Kt表示电磁力矩系数,J是电动机的转动惯量,La是电枢电感,B是粘性摩擦系数,s为拉普拉斯算子;运行实例中,选定位置随动***中常用直流电动机作为控制对象,不失一般性,其数学模型传递函数可表示为 其中,T为时间常数;
其特征在于所述直流电机鲁棒控制器参数整定方法,包括以下步骤:
(1)对于直流电机被控对象的数学模型传递函数P(s),其待整定FO[PI]控制器传递函数形式待整定参数为比例系数Kp,积分系数Ki和积分阶次λ,并给定需校正穿越频率ωc和需保持稳定的相位裕度φm;
(2)利用MATLAB画出被控对象P(s)的伯德图,求得在频率ωc处的模值m和相角n,同时可以求得被控对象频率ωc在相位变化率
(3)利用C(s)P(s)=G(s),鲁棒性条件:
得到:
注意到其中
利用MATLAB求解关于积分阶次λ的方程(3);
(4)利用稳定性条件:在开环***穿越频率ωc处相位裕度为φm;
C(jωc)P(jωc)=1∠φm-180° (4)
根据步骤(2)所求得的被控对象在频率ωc处的模值m和相角n,得到
由步骤(3)中
令A=10-m/20,从而得到
可以求得比例系数 和微分系数
(5)利用步骤(3)中MATLAB所得λ以及步骤(4)得到Kp和Ki代入即完成了分数阶FO[PI]鲁棒控制器参数整定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410201321.8A CN104049541B (zh) | 2014-05-13 | 2014-05-13 | 一种直流电机鲁棒控制器的参数整定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410201321.8A CN104049541B (zh) | 2014-05-13 | 2014-05-13 | 一种直流电机鲁棒控制器的参数整定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104049541A true CN104049541A (zh) | 2014-09-17 |
CN104049541B CN104049541B (zh) | 2016-12-07 |
Family
ID=51502569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410201321.8A Active CN104049541B (zh) | 2014-05-13 | 2014-05-13 | 一种直流电机鲁棒控制器的参数整定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104049541B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109946978A (zh) * | 2019-04-25 | 2019-06-28 | 广东省智能机器人研究院 | 一种考虑延时环节的伺服***分数阶模型辨识方法 |
CN110361964A (zh) * | 2019-07-15 | 2019-10-22 | 深圳大学 | 一种伺服驱动控制器参数整定方法 |
CN111473049A (zh) * | 2020-04-17 | 2020-07-31 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10011607A1 (de) * | 2000-03-10 | 2001-09-20 | Siemens Ag | Verfahren und Vorrichtung sowie Computerprogrammprodukt zum Betrieb einer technischen Anlage |
WO2007035559A2 (en) * | 2005-09-19 | 2007-03-29 | Cleveland State University | Controllers, observers, and applications thereof |
US20110084638A1 (en) * | 2009-10-12 | 2011-04-14 | Gm Global Technology Operations, Inc. | Methods, systems and apparatus for dynamically controlling an electric motor that drives an oil pump |
CN102944996A (zh) * | 2012-07-20 | 2013-02-27 | 长春理工大学 | 伺服***控制器及控制方法 |
CN102944995A (zh) * | 2012-07-20 | 2013-02-27 | 长春理工大学 | 伺服***控制器及控制方法 |
CN103092069A (zh) * | 2013-01-28 | 2013-05-08 | 上海交通大学 | 基于参数稳定域的PIλDμ控制器参数整定方法 |
CN103558755A (zh) * | 2013-11-05 | 2014-02-05 | 四川理工学院 | 分数阶积分pid控制器整定和自整定方法 |
-
2014
- 2014-05-13 CN CN201410201321.8A patent/CN104049541B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10011607A1 (de) * | 2000-03-10 | 2001-09-20 | Siemens Ag | Verfahren und Vorrichtung sowie Computerprogrammprodukt zum Betrieb einer technischen Anlage |
WO2007035559A2 (en) * | 2005-09-19 | 2007-03-29 | Cleveland State University | Controllers, observers, and applications thereof |
WO2007035559A3 (en) * | 2005-09-19 | 2009-05-22 | Univ State Cleveland | Controllers, observers, and applications thereof |
US20110084638A1 (en) * | 2009-10-12 | 2011-04-14 | Gm Global Technology Operations, Inc. | Methods, systems and apparatus for dynamically controlling an electric motor that drives an oil pump |
CN102944996A (zh) * | 2012-07-20 | 2013-02-27 | 长春理工大学 | 伺服***控制器及控制方法 |
CN102944995A (zh) * | 2012-07-20 | 2013-02-27 | 长春理工大学 | 伺服***控制器及控制方法 |
CN103092069A (zh) * | 2013-01-28 | 2013-05-08 | 上海交通大学 | 基于参数稳定域的PIλDμ控制器参数整定方法 |
CN103558755A (zh) * | 2013-11-05 | 2014-02-05 | 四川理工学院 | 分数阶积分pid控制器整定和自整定方法 |
Non-Patent Citations (1)
Title |
---|
L.H.KEEL 等: "Controller Synthesis Free of Analytical Models: Three Term Controllers", 《IEEE TRANSACTIONS ON AUTOMATIC CONTROL》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109946978A (zh) * | 2019-04-25 | 2019-06-28 | 广东省智能机器人研究院 | 一种考虑延时环节的伺服***分数阶模型辨识方法 |
CN109946978B (zh) * | 2019-04-25 | 2022-02-08 | 广东省智能机器人研究院 | 一种考虑延时环节的伺服***分数阶模型辨识方法 |
CN110361964A (zh) * | 2019-07-15 | 2019-10-22 | 深圳大学 | 一种伺服驱动控制器参数整定方法 |
CN110361964B (zh) * | 2019-07-15 | 2022-02-22 | 深圳大学 | 一种伺服驱动控制器参数整定方法 |
CN111473049A (zh) * | 2020-04-17 | 2020-07-31 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
CN111473049B (zh) * | 2020-04-17 | 2021-08-20 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104049541B (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104777746B (zh) | 一种增强型的增益鲁棒分数阶pid控制器参数整定方法 | |
CN106911274B (zh) | 一种原动机调速***附加阻尼器控制方法 | |
CN104503228B (zh) | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 | |
CN103338003A (zh) | 一种电机负载转矩及惯量同时在线辨识的方法 | |
CN114944663A (zh) | 一种基于构网型变流器的储能***控制方法及*** | |
CN103401256A (zh) | 火力发电厂大型汽轮机功频电液控制***参数的优化方法 | |
CN108599649A (zh) | Pmsm位置伺服***高阶对象控制器设计及参数确定方法 | |
CN101915203A (zh) | 改善水轮发电机组功角振荡的阻尼注入控制方法 | |
CN105511262A (zh) | 分数阶鲁棒控制器的参数快速整定方法及*** | |
CN106505591B (zh) | 一种多机电力***超低频振荡分析方法 | |
CN104049541A (zh) | 一种直流电机鲁棒控制器的参数整定方法 | |
CN105512367A (zh) | 并网火电机组一次调频转速不等率的临界稳定值确定方法 | |
CN103955143A (zh) | 一种交流永磁电机鲁棒控制器的参数整定方法 | |
CN104135205A (zh) | 一种感应电机最大转矩电流比控制方法 | |
CN111092440B (zh) | 一种用于协调控制水电机组一次调频与agc的方法及*** | |
CN103117692B (zh) | 多种外部干扰下的带有机械弹性储能的永磁电动机组控制方法 | |
CN106877769A (zh) | 一种伺服电机速度控制器增益参数自整定的方法 | |
CN103825516B (zh) | 一种同步电机的复合控制器 | |
CN107181437A (zh) | 一种同步发电机励磁控制方法 | |
CN107276474A (zh) | 一种基于模糊滑模的同步发电机励磁控制方法 | |
CN103605290B (zh) | 基于向量方法的鲁棒分数阶PDμ控制器的参数整定方法 | |
Van den Broeck et al. | Model predictive control for time-optimal point-to-point motion control | |
CN103631135B (zh) | 基于向量方法的分数阶PIλ控制器的参数整定方法 | |
CN101813053A (zh) | 改善水轮发电机组功角振荡的互联结构修正控制方法 | |
CN108983618A (zh) | 基于凸多胞形的pmsm鲁棒h∞输出反馈控制设计方法 |
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 |