CN103792849B - 可调金属切削系统的鲁棒自适应动态面控制方法 - Google Patents
可调金属切削系统的鲁棒自适应动态面控制方法 Download PDFInfo
- Publication number
- CN103792849B CN103792849B CN201410083566.5A CN201410083566A CN103792849B CN 103792849 B CN103792849 B CN 103792849B CN 201410083566 A CN201410083566 A CN 201410083566A CN 103792849 B CN103792849 B CN 103792849B
- Authority
- CN
- China
- Prior art keywords
- formula
- function
- model
- unknown
- metal cutting
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000002184 metal Substances 0.000 title claims abstract description 31
- 238000013461 design Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 230000008901 benefit Effects 0.000 claims abstract description 5
- 230000003044 adaptive effect Effects 0.000 claims description 14
- 238000012546 transfer Methods 0.000 claims description 13
- 230000007423 decrease Effects 0.000 claims description 7
- 230000001965 increasing effect Effects 0.000 claims description 7
- 210000005036 nerve Anatomy 0.000 claims description 7
- 230000007935 neutral effect Effects 0.000 claims description 6
- 241000208340 Araliaceae Species 0.000 claims description 5
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 5
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 5
- 230000006978 adaptation Effects 0.000 claims description 5
- 235000008434 ginseng Nutrition 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound 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 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000009472 formulation Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 210000002569 neuron Anatomy 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 abstract description 5
- 238000013459 approach Methods 0.000 abstract description 3
- 238000013528 artificial neural network Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 abstract 2
- 238000004880 explosion Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 238000011217 control strategy Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 210000004218 nerve net Anatomy 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
Abstract
一种可调金属切削系统的鲁棒自适应动态面控制方法,其特点是,包括可调金属切削系统数学模型、磁滞模型、神经网络系统描述和自适应动态面控制器的设计等内容,通过引入的误差转换函数可任意指定控制精度,实现了对切刀切割深度的精确控制;系统模型未完全可知的情况下,通过RBF神经网络实现未知项的逼近;通过估计未知参数向量的范数来代替估计未知参数地向量,极大地减轻了系统的计算负担;采用自适应动态面技术结合误差转换函数与RBF的策略,放宽了对时滞的假设同时保证了跟踪误差和过渡过程能够在任意指定的指标内,消除了反推法的“微分爆炸”问题,保证系统信号半全局一致最终有界,具有方法科学合理,控制精度高,适用等优点。
Description
技术领域
本发明是一种可调金属切削系统的鲁棒自适应动态面控制方法,应用于带有PI磁滞模型和未知时滞的精密仪器系统的控制。
背景技术
随着超精密加工技术在工业、军事、民品中的广泛应用,加工过程的高精度、高质量、高效率、低成本越来越重要。其中精密加工设备中的磁滞现象是影响系统稳定的关键,系统中伴有磁滞现象时,整个控制系统将会表现出振荡、不稳定、准确性降低等现象。除了磁滞现象,时延问题也出现在涡轮喷气发动机引擎、飞机系统等领域,从而降低了控制性能。
目前解决具有磁滞输入的非线性系统控制问题有两种方法:一种方法是建立磁滞的逆过程,并将之级联到控制系统的输入端;另一种方法是采用鲁棒自适应的控制策略来消除或抑制磁滞现象。由于建立磁滞逆过程非常复杂,且磁滞的逆模型输出对其参数特别敏感,因此,不建立磁滞逆而采用鲁棒自适应控制策略则更具吸引力,其中最主要最有效的方面就是自适应反推控制方案。
反推控制方案,在处理某些参数不确定非线性和线性系统及改善过渡过程品质方面的潜力而受到广泛重视。反推控制方案仍然存在一些重要缺陷,如随对象相对阶增高,控制律高度非线性、高度复杂(称为“微分爆炸”)等,近年来,研究人员为克服“微分爆炸”的缺点,提出了一种新的控制方法:“动态面控制(Dynamic Surface Control)”。即通过引入滤波器,使得每一步控制律的设计与前一级的设计基本解耦,从而使控制律的复杂程度大幅下降。动态面控制实现了相对简单的控制律与良好的过渡过程品质结合。
处理非线性时滞问题,目前常用的方法是自适应神经网络。由于使用动态面或反推控制,系统的跟踪误差最终收敛到一个紧集而不是收敛到零,因此一个挑战性的问题是如何设计控制律保证系统L∞性能。
发明内容
本发明的目的是针对既有时滞问题又存在磁滞输入的可调金属切削系统,提供一种可调金属切削系统的鲁棒自适应动态面控制方法,它采用鲁棒自适应动态面控制、误差转换函数和神经网络相结合,能够保证其跟踪误差和过渡过程能够在预先任意给定的指标内,克服了反推控制方案中的“微分爆炸”问题,简化控制器结构,减少计算量,更便于实时控制,并消除了系统跟踪误差可能出现的“畸变”现象。
实现本发明目的所采用的技术方案是:
一种可调金属切削系统的鲁棒自适应动态面控制方法,其特征是,它包括以下内容:
1)可调金属切削系统数学模型
可调金属切削系统的数学模型包括切割系统装置本身模型、用来描述磁滞作动器的磁滞模型及未知延时环节模型,
建立可调金属切削系统:
其中x表示切割深度的波动部分,也称作模具厚度的偏移量;m,c,k分别代表切割器械的质量、阻尼系数和弹性硬度;ka为压电作动器的弹性当量,u为压电作动器的输出,FΔ为切割器械的切力变化
FΔ=h×(x-μx(t-τ))+f(x-μx(t-τ))
其中 h和μ为正常数,τ为未知时间延迟;
基于上述对金属切削系统的描述,令x=x1,则式(1)表述成反馈形式:
y=x1
其中对于可调金属切削系统(1),作如下假设:
假设1:未知延时项满足如下不等式:
其中为未知连续函数;
假设2:参考信号yr光滑、有界,且对于t≥0,属于某一紧集;
2)磁滞模型
采用Prandtl-Ishlinskii模型来表示磁滞非线性,Prandtl-Ishlinskii模型是Preisach模型的一个子集,由stop算子和play算子组成,该模型的优势主要是能准确刻画磁滞现象,同时便于实现控制器的设计和实时控制,
由于滞环的非对称性,play算子使得磁滞输出w随输入u增加或减小,设Cm[0,tE]为一个单调分段连续的泛函空间,对于任意输入u,定义fr:R→R,fr(u,w)=max(u-r,min(u+r,w)),其中r表示滞环阈值且满足r≥0,则play算子Fr[·]定义为如下形式:
Fr[u](0)=fr(u(0),0)
Fr[u](t)=fr(u(t),Fr[u](ti)) (4)
for ti≤t≤ti+1and 0≤i≤N-1
其中,0=t0<t1<…<tN=tE为[0,tE]的一个分割,使得输入u在每一个[ti,ti+1]内具有单调性,即递增或者递减,PI模型表述为如下形式:
其中p(r)是密度函数,满足p(r)≥0且密度函数p(r)通过实验辨识得到,由于p(r)随着r的增大而趋于零,为简便,选择D=∞作为式(5)的积分上界,λ为由p(r)确定的未知常数,λ=1.5,r∈[0,10]、输入u(t)=7sin(3t)/(1+t),t∈[0,2π]且Fr[u](0)=0,得到PI磁滞模型;
磁滞输出表示为:
u(t)=λv(t)+∫0 Dp(r)Fr[v](t)dr (6)
将(6)式代入(2)得:
y=x1
其中:
β=bλ (8)
pg(r)=bp(r) (9)
且pg(r)为未知磁滞密度函数,pg(r)和β将在接下来被在线估计,针对上述带有磁滞驱动的非线性系统,控制目的是采用反推法的改进-动态面方法使得闭环系统保证全局稳定,同时,通过调整控制参数实现系统输出x1=y能够跟踪参考信号yr;
3)神经网络系统描述
一般地,神经元网络是一个多输入单输出的系统,其数学表达式为:
Y=θTψ(ξ) (10)
其中,ξ∈Rn是神经元网络的输入;Y∈R是神经元网络的输出,θ∈RN是N维可调参数向量,这里N是神经元节点数;ψ(ξ):Rn→RN是非线性向量函数,且ψ(ξ)=[ψ1(ξ),…,ψN(ξ)]T,这里,ψk(ξ)是高斯函数,
引理1:对紧集上的任意连续函数通过选择适当的参数σk和ζk,对于足够大的正整数N,存在形如式(10)的神经网络系统,使得,|δ|≤δm
其中,θ*是权参数向量θ的最优值,定义δ(ξ)为逼近误差,
4)自适应动态面控制器的设计
定义系统跟踪误差
e:=x1-yr (14)
其中yr为理想跟踪轨迹,性能函数和误差转换函数将进行如下定义:性能函数被定义为一个光滑递减的正函数,使得对于所有t≥0
其中0<σ<1且 为系统稳定时跟踪误差所允许的最大值,为将式(15)转换成一等价函数,引入误差转换函数如下:
其中,S1是由误差转换函数进行转换后的转换误差,Φ(S1)是某一光滑、严格单调递增的函数,其反函数具有如下性质:
且
由式(18)可知,如果那么式(17)成立,同时再考虑和式(16),当e(0)>0时,成立,或者当e(0)<0时,成立,从而公式(15)成立,因此,经分析可知,若欲实现给定的性能指标,只需证明即可,由Φ(S1)的严格单调递增的性质,可得:
需要注意的是,将e(0)=0时的情况包含到e(0)>0或e(0)<0中处理,同时,σ不可取零,因为这会使得S1(0)无界,
第1步:令第1个面误差S1为式(19)所定义,考虑式(1)和式(14),有
根据Φ和的定义,可知Ψ>0,考虑如下二次型方程
因此,可得
其中x2d为待设计的虚拟控制信号,设计虚拟控制律:
令x2d经一阶低通滤波器获得的新变量z2如下:
其中,τ2为滤波器时间常数;
第2步:定义第2个面误差为:
S2=x2-z2 (26)
令考虑二次型函数:
其中,定义:γζ,γpr和为正设计参数;分别为 的估计值;
对于未知延时项根据假设1,下述不等式成立:
因此,由式(27),可得:
其中为假 设1中的未知连续函数;定义为如下形式:
由于式(29)中有未知项,因此,需用神经网络作为逼近器,在紧集上来逼近未知项:
其中则有:
其中α2为正设计参数,根据式(28)-(32)有:
设最终的控制律υ为:
其中
和的含义由式(29))的调参律设计为如下形式
其中σζ,σpr和为正设计参数。
本发明的可调金属切削系统的鲁棒自适应动态面控制方法,由于采用鲁棒自适应动态面控制、误差转换函数和神经网络相结合,通过引入的误差转换函数可任意指定控制精度,实现了对切刀切割深度的精确控制;系统模型未完全可知的情况下,通过RBF神经网络实现未知项的逼近;通过估计未知参数向量的范数来代替估计未知参数地向量,极大地减轻了系统的计算负担;采用自适应动态面技术结合误差转换函数与RBF的策略,放宽了对时滞的假设同时保证了跟踪误差和过渡过程能够在任意指定的指标内,消除了反推法的“微分爆炸”问题,保证系统信号半全局一致最终有界,简化了控制器结构,减少计算量,更便于实时控制,并消除了系统跟踪误差可能出现的“畸变”现象,具有方法科学合理,控制精度高,适用等优点。
附图说明
图1为可调金属切削系统的结构示意图;
图2为可调金属切削系统的结构简化示意图;
图3为PI磁滞模型曲线示意图;
图4为跟踪性能y=x1,yr=0.02sin(2.5t)曲线示意图;
图5为跟踪误差,性能指标函数曲线示意图;
图6为控制律υ与磁滞输出曲线示意图;
图7为面误差S1与S2曲线示意图;
图8为权值向量范数的估计值曲线示意图。
图2中,x表示切割深度的波动部分,也称作模具厚度的偏移量,FΔ为切割器械的切力变化,τ为时间延迟;图3横坐标表示磁滞输入u(t)=7sin(3t)/(1+t),t∈[0,2π],纵坐标表示磁滞输出w(t);图4中的横坐标表示时间,纵坐标表示输出y对目标函数yr的跟踪性能;图5中的横坐标表示时间,纵坐标表示跟踪误差e和性能指标函数图6中的横坐标表示时间,纵坐标表示控制律υ与磁滞输出;图7中的横坐标表示时间,纵坐标表示面误差S1与S2;图8中的横坐标表示时间,纵坐标表示权值向量范数。
具体实施方式
下面利用附图和实施例对发明作进一步说明。
本发明的一种可调金属切削系统的鲁棒自适应动态面控制方法,包括:
1)可调金属切削系统数学模型
可调金属切削系统的数学模型包括切割系统装置本身模型、用来描述磁滞作动器的磁滞模型及未知延时环节模型,
建立可调金属切削系统:
其中x表示切割深度的波动部分,也称作模具厚度的偏移量;m,c,k分别代表切割器械的质量、阻尼系数和弹性硬度;ka为压电作动器的弹性当量,u为压电作动器的输出,FΔ为切割器械的切力变化
FΔ=h×(x-μx(t-τ))+f(x-μx(t-τ))
其中 h和μ为正常数,τ为未知时间延迟;
基于上述对金属切削系统的描述,令x=x1,则式(1)表述成反馈形式:
y=x1
其中对于可调金属切削系统(1),作如下假设:
假设1:未知延时项满足如下不等式:
其中为未知连续函数;
假设2:参考信号yr光滑、有界,且对于t≥0,属于某一紧集;
2)磁滞模型
采用Prandtl-Ishlinskii模型来表示磁滞非线性,Prandtl-Ishlinskii模型是Preisach模型的一个子集,由stop算子和play算子组成,该模型的优势主要是能准确刻画磁滞现象,同时便于实现控制器的设计和实时控制,
由于滞环的非对称性,play算子使得磁滞输出w随输入u增加或减小,设Cm[0,tE]为一个单调分段连续的泛函空间,对于任意输入u,定义fr:R→R,fr(u,w)=max(u-r,min(u+r,w)),其中r表示滞环阈值且满足r≥0,则play算子Fr[·]定义为如下形式:
Fr[u](0)=fr(u(0),0)
Fr[u](t)=fr(u(t),Fr[u](ti)) (4)
for ti≤t≤ti+1and 0≤i≤N-1
其中,0=t0<t1<…<tN=tE为[0,tE]的一个分割,使得输入u在每一个[ti,ti+1]内具有单调性,即递增或者递减,PI模型表述为如下形式:
w(t)=λu(t)+∫0 Dp(r)Fr[u](t)dr (5)
其中p(r)是密度函数,满足p(r)≥0且密度函数p(r)通过实验辨识得到,由于p(r)随着r的增大而趋于零,为简便,选择D=∞作为式(5)的积分上界,λ为由p(r)确定的未知常数,λ=1.5,r∈[0,10]、输入u(t)=7sin(3t)/(1+t),t∈[0,2π]且Fr[u](0)=0,得到PI磁滞模型;
磁滞输出表示为:
u(t)=λυ(t)+∫0 Dp(r)Fr[υ](t)dr (6)
将(6)式代入(2)得:
y=x1
其中:
β=bλ (8)
pg(r)=bp(r) (9)
且pg(r)为未知磁滞密度函数,pg(r)和β将在接下来被在线估计,针对上述带有磁滞驱动的非线性系统,控制目的是采用反推法的改进-动态面方法使得闭环系统保证全局稳定,同时,通过调整控制参数实现系统输出x1=y能够跟踪参考信号yr;
3)神经网络系统描述
一般地,神经元网络是一个多输入单输出的系统,其数学表达式为:
Y=θTψ(ξ) (10)其中,ξ∈Rn是神经元网络的输入;Y∈R是神经元网络的输出,θ∈RN是N维可调参数向量,这里N是神经元节点数;ψ(ξ):Rn→RN是非线性向量函数,且ψ(ξ)=[ψ1(ξ),…,ψN(ξ)]T,这里,ψk(ξ)是高斯函数,
引理1:对紧集上的任意连续函数通过选择适当的参数σk和ζk,对于足够大的正整数N,存在形如式(10)的神经网络系统,使得,|δ|≤δm
其中,θ*是权参数向量θ的最优值,定义δ(ξ)为逼近误差,
4)自适应动态面控制器的设计
定义系统跟踪误差
e:=x1-yr (14)其中yr为理想跟踪轨迹,性能函数和误差转换函数将进行如下定义:性能函数被定义为一个光滑递减的正函数,使得对于所有t≥0
其中0<σ<1且 为系统稳定时跟踪误差所允许的最大值,为将式(15)转换成一等价函数,引入误差转换函数如下:
其中,S1是由误差转换函数进行转换后的转换误差,Φ(S1)是某一光滑、严格单调递增的函数,其反函数具有如下性质:
且
由式(18)可知,如果那么式(17)成立,同时再考虑和式(16),当e(0)>0时,成立,或者当e(0)<0时,成立,从而公式(15)成立,因此,经分析可知,若欲实现给定的性能指标,只需证明即可,由Φ(S1)的严格单调递增的性质,可得:
需要注意的是,将e(0)=0时的情况包含到e(0)>0或e(0)<0中处理,同时,σ不可取零,因为这会使得S1(0)无界,
第1步:令第1个面误差S1为式(19)所定义,考虑式(1)和式(14),有
根据Φ和的定义,可知Ψ>0,考虑如下二次型方程
因此,可得
其中x2d为待设计的虚拟控制信号,设计虚拟控制律:
令x2d经一阶低通滤波器获得的新变量z2如下:
其中,τ2为滤波器时间常数;
第2步:定义第2个面误差为:
S2=x2-z2 (26)
令考虑二次型函数:
其中,定义:γζ,γpr和为正设计参数;分别为 的估计值;
对于未知延时项根据假设1,下述不等式成立:
因此,由式(27),可得:
其中为假 设1中的未知连续函数;定义为如下形式:
由于式(29)中有未知项,因此,需用神经网络作为逼近器,在紧集上来逼近未知项:
其中则有:
其中α2为正设计参数,根据式(28)-(32)有:
设最终的控制律υ为:
其中
和的含义由式(29))的调参律设计为如下形式
其中σζ,σpr和为正设计参数。
具体实施例:
具有磁滞输入及延时现象的金属切削系统的数学模型如式(1)、(2)所示,式中PI磁滞模型如式(6)所示,且λ=8。系统参数为m=26kg,c=1.2,k=6800N/m,f(x-μx(t-τ))=0.18sin(x-μx(t-1)),μ=1,h=4500,控制目的是使系统输出x1=y能够跟踪参考信号yr=0.02sin(2.5t)。神经网络的输入向量为神经网络的基函数ψ(ξ1)选择31个节点,中心点σj,j=1,…,31,均匀分布在[-3,3]×[-3,3]×[-3,3]]且宽度为ηj=1。为基函数ψ(ξ2)选择61个节点,中心点ζj,j=1,…,61,均匀分布在[-3,3]×[-3,3]×[-3,3],且宽度为σj=1。
系统中各个参数的调参律的初始值为此外,设计参数k1=7.5,k2=2,α2=2.5,γζ=1.5,γpr=3;滤波器的时间常数τ2=0.15,小增益参数σpr=0.00001;系统各个状态的初始值为x1(0)=-0.1,x2(0)=0。仿真中,给定性能指标函数
为避免使控制信号避免产生“抖震”现象,使用:
来代替符号函数sgn(S2/0.01)。
仿真结果如图4到图8所示。图4中良好的跟踪性能表明,本发明的可调金属切削系统的鲁棒自适应动态面控制方法可很好克服磁滞现象、时滞现象以及外部扰动对控制系统的不良影响。图5展示了系统的跟踪误差不超过0.05按照指定性能指标收敛,保持在预先指定的性能函数的约束范围内。图6给出了系统的控制信号υ(t)和在此信号下磁滞输出u(t)。图7给出了面误差信号S1和S2,S1的有界性恰好说明了跟踪误差将按指定的性能函数来收敛。图8为权值向量范数的估计值图4到图8说明本发明的可调金属切削系统的鲁棒自适应动态面控制方法可以很好的处理输出磁滞现象与系统状态的延时,证明了所提控制策略的可行性。
具体实施方式仅是对本发明的说明,并不构成对权利要求保护范围的限制,本领域技术人员不经过创造性劳动的等同替代,均在本发明保护范围内。
Claims (1)
1.可调金属切削系统的鲁棒自适应动态面控制方法,其特征是,它包括以下内容:
1)可调金属切削系统数学模型
可调金属切削系统的数学模型包括切割系统装置本身模型、用来描述磁滞作动器的磁滞模型及未知延时环节模型,
建立可调金属切削系统:
其中x表示切割深度的波动部分,也称作模具厚度的偏移量;m,c,k分别代表切割器械的质量、阻尼系数和弹性硬度;ka为压电作动器的弹性当量,u为压电作动器的输出,FΔ为切割器械的切力变化
FΔ=h×(x-μx(t-τ))+f(x-μx(t-τ))
其中 h和μ为正常数,τ为未知时间延迟;
基于上述对金属切削系统的描述,令x=x1,则式(1)表述成反馈形式:
y=x1
其中对于可调金属切削系统(1),作如下假设:
假设1:未知延时项i=1,…,n,满足如下不等式:
其中为未知连续函数;
假设2:参考信号yr光滑、有界,且对于t≥0,属于某一紧集;
2)磁滞模型
采用Prandtl-Ishlinskii模型来表示磁滞非线性,Prandtl-Ishlinskii模型是Preisach模型的一个子集,由stop算子和play算子组成,该模型的优势主要是能准确刻画磁滞现象,同时便于实现控制器的设计和实时控制,
由于滞环的非对称性,play算子使得磁滞输出w随输入u增加或减小,设Cm[0,tE]为一个单调分段连续的泛函空间,对于任意输入u,定义fr:R→R, fr(u,w)=max(u-r,min(u+r,w)),其中r表示滞环阈值且满足r≥0,则play算子Fr[·]定义为如下形式:
Fr[u](0)=fr(u(0),0)
Fr[u](t)=fr(u(t),Fr[u](ti)) (4)
for ti≤t≤ti+1and 0≤i≤N-1
其中,0=t0<t1<…<tN=tE为[0,tE]的一个分割,使得输入u在每一个[ti,ti+1]内具有单调性,即递增或者递减,PI模型表述为如下形式:
其中p(r)是密度函数,满足p(r)≥0且密度函数p(r)通过实验辨识得到,由于p(r)随着r的增大而趋于零,为简便,选择D=∞作为式(5)的积分上界,λ为由p(r)确定的未知常数,λ=1.5,r∈[0,10]、输入u(t)=7sin(3t)/(1+t),t∈[0,2π]且Fr[u](0)=0,得到PI磁滞模型;
磁滞输出表示为:
将(6)式代入(2)得:
y=x1
其中:
β=bλ (8)
pg(r)=bp(r) (9)
且pg(r)为未知磁滞密度函数,pg(r)和β将在接下来被在线估计,针对上述带有磁滞驱动的非线性系统,控制目的是采用反推法的改进-动态面方法使得闭环系统保证全局稳定,同时,通过调整控制参数实现系统输出x1=y能够跟踪参考信号yr;
3)神经网络系统描述
一般地,神经元网络是一个多输入单输出的系统,其数学表达式为:
Y=θTψ(ξ) (10)
其中,ξ∈Rn是神经元网络的输入;Y∈R是神经元网络的输出,θ∈RN是N维可调参数 向量,这里N是神经元节点数;ψ(ξ):Rn→RN是非线性向量函数,且ψ(ξ)=[ψ1(ξ),…,ψN(ξ)]T,这里,ψk(ξ)是高斯函数,
引理1:对紧集上的任意连续函数通过选择适当的参数σk和ζk,对于足够大的正整数N,存在形如式(10)的神经网络系统,使得,
其中,θ*是权参数向量θ的最优值,定义δ(ξ)为逼近误差,
4)自适应动态面控制器的设计
定义系统跟踪误差
e:=x1-yr (14)
其中yr为理想跟踪轨迹,性能函数和误差转换函数将进行如下定义:性能函数被定义为一个光滑递减的正函数,使得对于所有t≥0
其中0<σ<1且 为系统稳定时跟踪误差所允许的最大值,为将式(15)转换成一等价函数,引入误差转换函数如下:
其中,S1是由误差转换函数进行转换后的转换误差,Φ(S1)是某一光滑、严格单调递增的函数,其反函数具有如下性质:
且
由式(18)可知,如果那么式(17)成立,同时再考虑和式(16),当e(0)>0时,成立,或者当e(0)<0时,成立,从而公式(15)成立,因此,经分析可知,若欲实现给定的性能指标,只需证明即可,由Φ(S1)的严格单调递增的性质,可得:
需要注意的是,将e(0)=0时的情况包含到e(0)>0或e(0)<0中处理,同时,σ不可取零,因为这会使得S1(0)无界,
第1步:令第1个面误差S1为式(19)所定义,考虑式(1)和式(14),有
根据Φ和的定义,可知Ψ>0,考虑如下二次型方程
因此,可得
其中x2d为待设计的虚拟控制信号,设计虚拟控制律:
令x2d经一阶低通滤波器获得的新变量z2如下:
其中,τ2为滤波器时间常数;
第2步:定义第2个面误差为:
S2=x2-z2 (26)
令考虑二次型函数:
其中,定义:γζ,γpr和为正设计参 数;分别为 的估计值;
对于未知延时项根据假设1,下述不等式成立:
因此,由式(27),可得:
其中为假设1中的未知连续函数;定义为如下形式:
由于式(29)中有未知项,因此,需用神经网络作为逼近器,在紧集上来逼近未知项:
其中则有:
其中α2为正设计参数,根据式(28)-(32)有:
设最终的控制律v为:
其中
和的含义由式(29))的调参律设计为如下形式
其中σζ,σpr和为正设计参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410083566.5A CN103792849B (zh) | 2014-03-09 | 2014-03-09 | 可调金属切削系统的鲁棒自适应动态面控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410083566.5A CN103792849B (zh) | 2014-03-09 | 2014-03-09 | 可调金属切削系统的鲁棒自适应动态面控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103792849A CN103792849A (zh) | 2014-05-14 |
| CN103792849B true CN103792849B (zh) | 2017-02-08 |
Family
ID=50668643
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410083566.5A Expired - Fee Related CN103792849B (zh) | 2014-03-09 | 2014-03-09 | 可调金属切削系统的鲁棒自适应动态面控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103792849B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104749952B (zh) * | 2015-04-16 | 2017-10-17 | 大连海事大学 | 一种自适应动态面控制器结构及设计方法 |
| CN106054670B (zh) * | 2016-07-18 | 2019-05-14 | 华南理工大学 | 一种基于时滞的超磁致驱动器回滞建模方法 |
| CN106773694B (zh) * | 2016-12-26 | 2019-10-22 | 东北电力大学 | 压电精密位置平台自适应输出反馈逆控制方法 |
| CN111766775B (zh) * | 2020-05-06 | 2023-03-28 | 东北电力大学 | 具有未知饱和pi回滞的非线性系统动态面隐逆控制器 |
| CN114384800B (zh) * | 2021-12-09 | 2023-09-12 | 上海工程技术大学 | 一种具有输入信号延时的未知非线性系统反推控制方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012006431A2 (en) * | 2010-07-09 | 2012-01-12 | Temple University - Of The Commonwealth System Of Higher Education | Apparatus and method for surface and subsurface tactile sensation imaging |
| CN102522945A (zh) * | 2012-01-10 | 2012-06-27 | 江苏大学 | 基于多神经网络逆模型的多相电机容错控制方法及系统 |
| CN102540882A (zh) * | 2012-03-01 | 2012-07-04 | 北京航空航天大学 | 一种基于最小参数学习法的飞行器航迹倾角控制方法 |
| CN102566417A (zh) * | 2012-02-17 | 2012-07-11 | 南京电力设备质量性能检验中心 | 一种柔性关节机械臂的动态面控制方法 |
| CN102591203A (zh) * | 2012-02-17 | 2012-07-18 | 南京国电环保设备有限公司 | 一种伺服电机的基于微分器的直接神经网络控制方法 |
| CN102749843A (zh) * | 2012-07-24 | 2012-10-24 | 大连海事大学 | 一种自适应反馈保护动态面控制器结构及设计方法 |
-
2014
- 2014-03-09 CN CN201410083566.5A patent/CN103792849B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012006431A2 (en) * | 2010-07-09 | 2012-01-12 | Temple University - Of The Commonwealth System Of Higher Education | Apparatus and method for surface and subsurface tactile sensation imaging |
| CN102522945A (zh) * | 2012-01-10 | 2012-06-27 | 江苏大学 | 基于多神经网络逆模型的多相电机容错控制方法及系统 |
| CN102566417A (zh) * | 2012-02-17 | 2012-07-11 | 南京电力设备质量性能检验中心 | 一种柔性关节机械臂的动态面控制方法 |
| CN102591203A (zh) * | 2012-02-17 | 2012-07-18 | 南京国电环保设备有限公司 | 一种伺服电机的基于微分器的直接神经网络控制方法 |
| CN102540882A (zh) * | 2012-03-01 | 2012-07-04 | 北京航空航天大学 | 一种基于最小参数学习法的飞行器航迹倾角控制方法 |
| CN102749843A (zh) * | 2012-07-24 | 2012-10-24 | 大连海事大学 | 一种自适应反馈保护动态面控制器结构及设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103792849A (zh) | 2014-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103792849B (zh) | 可调金属切削系统的鲁棒自适应动态面控制方法 | |
| Dulău et al. | Fractional order controllers versus integer order controllers | |
| CN110161841A (zh) | 一种适用于暂冲式跨声速风洞的前馈-模糊pid控制方法 | |
| CN103217902B (zh) | 一种基于干扰观测器的指令滤波反步控制方法 | |
| CN103207568B (zh) | 一种抗舵机饱和的船舶航向自适应控制方法 | |
| CN106647271A (zh) | 基于神经网络理论的非线性系统自适应比例积分控制方法 | |
| CN106444368A (zh) | 具有输入非线性的近空间飞行器预设性能姿态跟踪控制方法 | |
| van Duist et al. | FLOreS-Fractional order loop shaping MATLAB toolbox | |
| Adnan et al. | Self-tuning fuzzy PID controller for electro-hydraulic cylinder | |
| Zhmud et al. | Calculation of PID-regulator for MISO system with the method of numerical optimization | |
| Durmaz et al. | Sliding mode control for non-linear systems with adaptive sliding surfaces | |
| Kim et al. | PID controller design with time response specifications | |
| CN103809446B (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
| CN110609568A (zh) | 一种大型无人机uav的强自耦pi协同控制方法 | |
| CN109828586A (zh) | 基于非线性不确定性的无人船舶鲁棒h∞航向控制方法 | |
| Di Teodoro et al. | Sliding-mode controller based on fractional order calculus for chemical processes | |
| Var et al. | An internal model control based design method for single input fuzzy PID controllers | |
| Lin et al. | Robust fuzzy-neural sliding-mode controller design via network structure adaptation | |
| Yuasa et al. | Optimal control designs for systems with input saturations and rate limiters | |
| CN117724340A (zh) | 一种基于分数阶自适应步长固定时间二阶优化算法的无人机pid控制律参数整定方法 | |
| Rabiei et al. | Constrained multivariable dynamic matrix control for a class of fractional-order system | |
| CN115857320A (zh) | 单输入单输出非线性未知动力系统的优化跟踪控制方法 | |
| Deniz et al. | Experimental verification of lead-lag compensators on a twin rotor system | |
| Huzmezan | Theory and aerospace applications of constrained model based predictive control | |
| CN106647270A (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 | ||
| 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: 20170208 |