CN106125574A - 基于dpi模型的压电陶瓷微定位平台建模方法 - Google Patents
基于dpi模型的压电陶瓷微定位平台建模方法 Download PDFInfo
- Publication number
- CN106125574A CN106125574A CN201610583363.1A CN201610583363A CN106125574A CN 106125574 A CN106125574 A CN 106125574A CN 201610583363 A CN201610583363 A CN 201610583363A CN 106125574 A CN106125574 A CN 106125574A
- Authority
- CN
- China
- Prior art keywords
- dpi
- model
- operator
- output
- piezoelectric ceramics
- 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
- 239000000919 ceramic Substances 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000013528 artificial neural network Methods 0.000 claims description 39
- 230000006870 function Effects 0.000 claims description 24
- 210000002569 neuron Anatomy 0.000 claims description 17
- 238000012549 training Methods 0.000 claims description 11
- 238000013519 translation Methods 0.000 claims description 10
- 230000005284 excitation Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 230000003585 interneuronal effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000006073 displacement reaction Methods 0.000 description 24
- 238000010586 diagram Methods 0.000 description 19
- 238000002474 experimental method Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000003247 decreasing effect Effects 0.000 description 5
- 238000005312 nonlinear dynamic Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000004218 nerve net Anatomy 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000036316 preload Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Geometry (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
一种基于DPI模型的压电陶瓷微定位平台建模方法,属于控制工程技术领域。本发明的目的是建立DPI模型对压电陶瓷微定位平台进行动态迟滞建模的基于DPI模型的压电陶瓷微定位平台建模方法。本发明根据一阶微分方程输入输出的类迟滞特性,对算子的表达式进行一阶微分处理构建DPI算子,DPI迟滞算子,再通过一系列DPI算子与对应的权值的乘积再求和获得DPI模型。本发明DPI模型与传统PI模型的实验结果对比,可以看出,DPI模型是与输入电压频率相关的动态模型,可以更好的描述压电陶瓷微定位平台迟滞非线性的动态特性,且适用于幅值有变化的任何驱动电压信号,为压电陶瓷微定位平台在精密定位控制中应用奠定了基础。
Description
技术领域
本发明属于控制工程技术领域。
背景技术
随着纳米技术的飞速发展,压电陶瓷微定位平台逐渐成为精密制造装备中的核心器件,具有分辨率高,响应速度快,输出位移大等优点。但是固有的迟滞非线性严重影响了其定位精度。为了提高压电陶瓷微定位平台的定位精度,选择合适的方法消除迟滞非线性成为很重手段。近年来,针对迟滞非线性现象,国内外学者就迟滞非线性建模和控制等方面开展了大量的研究。Xu Q等人针对一种压电驱动显微操作平台,提出了一种新的迟滞逆补偿控制方法,首先建立平台的Dahl模型,之后提出了一种基于逆模型的前馈控制,仿真结果显示,前馈控制的最大误差为5.69μm,输入与输出之间呈现一一对应的线性关系。Lu Y等人为了消除压电陶瓷微定位平台的迟滞非线性对控制精度的影响,首先建立平台的Preisach模型,并设计了一种基于Preisach模型的前馈控制器,仿真结果表明,前馈控制最大误差为0.32μm。田艳兵等人根据超精密定位平台的特性,首先建立了平台的迟滞模型,设计了前馈控制,并在此基础上加上PID控制作为反馈控制,将平台的定位误差减小到0.72μm,平均误差为0.45μm。Li Y等人以压电驱动微操作平台为研究对象,首先建立Bouc-Wen模型描述其迟滞非线性,并用粒子群优化算法对模型进行辨识。魏强等人针对显微镜工作台的高精度要求,提出了一种基于神经网络的PID控制,用三层神经网络自学习PID代替传统的PID,并采用BP算法对神经网络进行在线训练,与常规PID相比,误差率由2.78%减小到1.39%。
PI模型是由俄国数学家Krasnosel’skii在Preisach迟滞算子的基础上提出。该模型是由基本迟滞算子通过加权叠加,从而建立被控对象的迟滞模型。意大利的ScipioneBobbio等对PI模型与Preisach模型进行比较,并将play算子与stop算子用于迟滞现象建模中。PI模型具有多值映射和非局部记忆的性质,计算量较小,能够应用到实时控制中等优点,因此在对复杂迟滞非线性特性表述方面被广泛应用。压电陶瓷微定位平台的迟滞特性依赖于输入电压信号的频率,是一种动态的迟滞特性,而传统的PI模型是一种静态模型,无法精确描述压电陶瓷微定位平台的迟滞非线性动态特性,对传统PI模型进行改进以适应实际的需求成为一项很大挑战。
发明内容
本发明的目的是建立DPI模型对压电陶瓷微定位平台进行动态迟滞建模的基于DPI模型的压电陶瓷微定位平台建模方法。
本发明DPI模型的建立:
根据一阶微分方程输入输出的类迟滞特性,对算子的表达式进行一阶微分处理构建DPI算子,DPI迟滞算子的最终表达式为:
(3)
式中,为DPI迟滞算子在t时刻的输出值,为DPI算子在t-1时刻的输出值;为DPI算子的阈值,为惯性因子;
DPI模型阈值采取非等分阈值确定表述为:
(4)
式中,为DPI算子输入,为常数,,分别代表前半升程和后半升程第个和第算子的序号;
DPI模型为:
(5)
式中,为DPI模型t时刻的输出值,为迟滞算子的个数,为DPI模型在t时刻第个迟滞算子输出,为DPI模型第个迟滞算子对应的权值。
本发明DPI模型权值参数的求取:
将已求得的个DPI迟滞算子值作为小波神经网络的输入即,其中分别代表DPI动态迟滞模型中的第1个,第2个,第m个DPI算子的输出;定义网络输入层的第个输入到隐含层第个神经元间的权值为,隐含层第个神经元到输出层的第个输出之间的权值为,,,分别代表小波神经网络的输入层,隐含层,输出层的神经元序号;t时刻小波神经网络的第个输出为:
(6)
其中,代表隐含层神经元总数,代表隐含层的激励函数,通过公式(7)—(9)求得:
(7)
(8)
(9)
其中,式(9)为Mymorlet小波函数,代表输入层神经元的总数即DPI算子的个数;此处定义为小波函数的伸缩因子,为小波函数的平移因子,则公式(8)的含义为将输入的DPI算子乘以输入层与隐含层的权值后做平移和伸缩变换,在通过公式(7)计算求出激励函数;
在辨识DPI权值参数时,小波神经网络的误差函数采取最小平方误差的形式,表示为::
(10)
式中,代表t时刻实验数据即为理想参考输出,代表着t时刻DPI模型的输出,小波神经网络权值、,节点伸缩参数、平移参数在每次训练过程中对应的修正量,,,分别按如下公式求得:
(11)
(12)
(13)
(14)
依照梯度下降法,小波神经网络t时刻的对应参数调整输出根据以下公式求得:
(15)
(16)
(17)
(18)
式中,为学习系数,为动量因子;
最终求得的DPI动态迟滞模型所有的权值参数,其中代表第1个算子对应的权值,代表第个算子对应的权值,代表第个算子对应的权值。
本发明小波神经网络辨识过程:
1)网络训练初始,设定目标精度,设定最大迭代次数运用随机函数给权值、,节点伸缩参数、平移参数赋初值;
2)然后计算网络的所有输出,并与理想的参考信号进行做差比较求取误差函数;
3)判断误差函数的大小是否小于等于设定的目标精度,如果小于则执行步骤6),否则执行步骤4);
4)依据误差函数,计算新的参数的修正量,,,和对应的参数值、,、;
5)判断迭代次数是否达到设定的最大值,未到达则转移到步骤2),否则转移到步骤1)重新随机幅值进行训练;
6)停止训练,保存训练得到的参数值;
重复执行步骤1)~6),最终得到满足精度要求的网络权值、,节点伸缩参数、平移参数的取值,此时的小波神经网络结构的功能等同于DPI模型算子对应的权值的功能,即用已经训练好的小波神经网络代替权值矩阵,完成DPI模型权值的辨识。
本发明DPI模型与传统PI模型的实验结果对比,可以看出,DPI模型是与输入电压频率相关的动态模型,可以更好的描述压电陶瓷微定位平台迟滞非线性的动态特性,且适用于幅值有变化的任何驱动电压信号,为压电陶瓷微定位平台在精密定位控制中应用奠定了基础。
附图说明
图1压电陶瓷微定位平台实验原理图;
图2是压电陶瓷微定位平台实验装置图;
图3是单一幅值正弦输入电压信号的迟滞特性曲线;
图4是幅值递减正弦输入电压信号的迟滞特性曲线;
图5输入频率为0.1Hz时实际迟滞曲线与模型输出曲线对比图;
图6输入频率为1Hz时实际迟滞曲线与模型输出曲线对比图;
图7输入频率为5Hz时实际迟滞曲线与模型输出曲线对比图;
图8输入频率为10Hz时实际迟滞曲线与模型输出曲线对比图;
图9输入频率为0.1Hz时实际迟滞曲线与模型误差曲线对比图;
图10输入频率为1Hz时实际迟滞曲线与模型误差曲线对比图;
图11输入频率为5Hz时实际迟滞曲线与模型误差曲线对比图;
图12输入频率为10Hz时实际迟滞曲线与模型误差曲线对比图;
图13输入频率为0.1Hz时实际迟滞曲线与模型输出曲线对比图;
图14输入频率为1Hz时实际迟滞曲线与模型输出曲线对比图;
图15输入频率为5Hz时实际迟滞曲线与模型输出曲线对比图;
图16输入频率为10Hz时实际迟滞曲线与模型输出曲线对比图;
图17输入频率为0.1Hz时实际迟滞曲线与模型输出误差对比图;
图18输入频率为1Hz时实际迟滞曲线与模型输出误差对比图;
图19输入频率为5Hz时实际迟滞曲线与模型输出误差对比图;
图20输入频率为10Hz时实际迟滞曲线与模型输出误差对比图;
图21是传统PI迟滞算子结构图;
图22是PI模型结构图;
图23是DPI迟滞算子结构图;
图24是小波神经网络结构图。
具体实施方式
压电陶瓷微定位平台实验原理如图1,实验装置的真实硬件环境如图2。其中包括压电陶瓷微定位平台,型号为MPT-2MRL102A,集成精密定位控制台,型号为PPC-2CR0150,具有A/D,D/A转换的数据采集卡PCI1710,具有MATLAB实时工作空间RTW测试环境的PC机一台。在PC机上编写输出电压驱动信号,通过集成精密定位控制台施加到压电陶瓷微定位平台,使其产生位移,由集成精密定位控制台内部集成的位移传感器测量,通过采集卡实时采集位移在PC机上进行显示。驱动电压与位移的迟滞非线性表现分别用传统PI模型和DPI模型描述。
传统PI迟滞建模方法:PI迟滞模型选用Play迟滞算子作为其基本PI迟滞算子,该算子结构简单,只有一个未知参数。PI迟滞算子结构如下图21所示,数学描述形式为:
(1)
式中,是算子输入,代表着实际的输入电压信号;为PI模型迟滞算子对应的阈值;为PI迟滞算子当前t时刻输出值,为迟滞算子t-1刻输出值。
传统PI模型是多个PI算子的加权叠加,如图22所示,其表达式可以描述为:
(2)
式中,为模型的输出,为模型的输入,为play算子的个数,为第个算子的权值。
传统PI模型是一种静态模型,模型输出不受输入信号频率变化的影响。然而压电陶瓷微定位平台的迟滞非线性呈现的是一种动态特性,传统的PI模型无法准确描述平台的迟滞非线性特性。
本发明DPI迟滞动态建模方法:
根据一阶微分方程输入输出的类迟滞特性,对算子的表达式进行一阶微分处理构建DPI算子,DPI迟滞算子的最终表达式为:
(3)
式中,为DPI迟滞算子在t时刻的输出值,为DPI算子在t-1时刻的输出值;为DPI算子的阈值,决定了DPI算子的形状,在实际建模过程中需要求出各个阈值的大小,确定DPI算子;为惯性因子;由输入信号的频率决定,是模型具有动态特性的关键参数。改进后的DPI算子结构如图23所示。
可以看出改进后的DPI算子更接近于压电陶瓷微定位平台的迟滞非线性特性。传统PI算子的阈值一般采取等分取值方法确定,而DPI模型阈值采取非等分阈值确定,可以表述为:
(4)
式中,为算子输入,实际为输入的电压值,为常数,根据平台实际情况以及算子总数决定,近似等于DPI算子总数的三分之一。等式(4)表明:平台初载曲线中前半升程的算子个数为个,后半升程的算子个数为个,,分别代表前半升程和后半升程第个和第算子的序号。
DPI模型即为一系列DPI算子与对应的权值的乘积再求和,其表达式可以描述为:
(5)
式中,为DPI模型t时刻的输出值,为迟滞算子的个数,为DPI模型在t时刻第个迟滞算子输出,为DPI模型第个迟滞算子对应的权值。实际应用中需要运用优化算法辨识得到权值。
本发明DPI模型权值参数的求取,采用小波神经网络进行辨识:
小波神经网络辨识算法:
为了便于进行DPI动态迟滞模型与传统PI模型的精度比较,两种模型的权值参数都用网络结构完全相同的小波神经网络进行辨识。下面以辨识DPI动态迟滞模型的权值密度进行辨识过程的说明,辨识传统PI模型的权值参数过程相同不在赘述。
关于对DPI动态迟滞模型的权值参数的求取,采用小波神经网络进行辨识。本节选取的小波神经网络结构图如下图24所示。
将已求得的个DPI迟滞算子值作为小波神经网络的输入即,其中分别代表DPI动态迟滞模型中的第1个,第2个,第m个DPI算子的输出;定义网络输入层的第个输入到隐含层第个神经元间的权值为,隐含层第个神经元到输出层的第个输出之间的权值为,,,分别代表小波神经网络的输入层,隐含层,输出层的神经元序号;t时刻小波神经网络的第个输出为:
(6)
其中,代表隐含层神经元总数,代表隐含层的激励函数,通过公式(7)—(9)求得:
(7)
(8)
(9)
其中,式(9)为Mymorlet小波函数,代表输入层神经元的总数即DPI算子的个数;此处定义为小波函数的伸缩因子,为小波函数的平移因子,则公式(8)的含义为将输入的DPI算子乘以输入层与隐含层的权值后做平移和伸缩变换,在通过公式(7)计算求出激励函数;
由于在t时刻,网络的输出即DPI动态迟滞模型的输出与理想的参考位移不完全相等,存在着一定的位移误差,为了保证模型精度,小波神经网络需要根据此误差进行权值的调整,所以在辨识DPI权值参数时,小波神经网络的误差函数采取最小平方误差的形式,表示为:
(10)
式中,代表t时刻实验数据即为理想参考输出,代表着t时刻DPI模型的输出,小波神经网络权值、,节点伸缩参数、平移参数在每次训练过程中对应的修正量,,,分别按如下公式求得:
(11)
(12)
(13)
(14)
依照梯度下降法,小波神经网络t时刻的对应参数调整输出根据以下公式求得:
(15)
(16)
(17)
(18)
式中,为学习系数,为动量因子。
运用常规的辨识算法最终求得的DPI动态迟滞模型所有的权值参数构成一个维的数值矩阵,其中代表第1个算子对应的权值,代表第个算子对应的权值,代表第个算子对应的权值。但神经网络辨识是模拟人脑的神经元,神经元之间相互连,每个神经元接受数据,判断数据产生信号传递给下一个神经元,然后根据误差反向传递调整参数,最终达到辨识的目的。运用小波神经网络辨识出的DPI模型的权值参数并不是一个简单的数值,而是特定结构的小波神经网络,因其具有自学习和自调整能力,可以保证DPI动态迟滞模型具有更高的精度。
小波神经网络权值和阈值确定之后,对网络进行前向传播即可得到误差值,并通过式(11)-(18)对网络参数进行调整。反复训练,直到误差满足精度条件或者达到最大迭代次数时停止。
本发明小波神经网络辨识过程:
1)网络训练初始,设定目标精度,设定最大迭代次数运用随机函数给权值、,节点伸缩参数、平移参数赋初值;
2)然后计算网络的所有输出,并与理想的参考信号进行做差比较求取误差函数;
3)判断误差函数的大小是否小于等于设定的目标精度,如果小于则执行步骤6),否则执行步骤4);
4)依据误差函数,计算新的参数的修正量,,,和对应的参数值、,、;
5)判断迭代次数是否达到设定的最大值,未到达则转移到步骤2),否则转移到步骤1)重新随机幅值进行训练;
6)停止训练,保存训练得到的参数值;
重复执行步骤1)~6),最终得到满足精度要求的网络权值、,节点伸缩参数、平移参数的取值,此时的小波神经网络结构的功能等同于DPI模型算子对应的权值的功能,即用已经训练好的小波神经网络代替权值矩阵,完成DPI模型权值的辨识。
验证:
实际训练中,根据平台采样周期和算子规模大小,确定输入层神经元总数为m=60,隐含层神经元总数和输出层神经元总数。
为了验证两种建模方法对压电陶瓷微定位平台迟滞非线性动态特性描述的准确性,在Matlab环境下应用小波神经网络分别辨识出PI迟滞模型,DPI动态迟滞模型的权值参数。以下就前期的实验准备,建模的数据获得方法,DPI动态迟滞模型与传统PI模型精度比较和结论等方面进行详细描述。
压电陶瓷微定位平台实验的硬件设备主要包括:
①压电陶瓷微定位平台,型号为MPT-2MRL102A。压电陶瓷微定位平台是实验对象,通过外加载驱动电压信号,平台就会产生水平方向的输出位移。该实验对象可承受的最大外部电压为-20V-150V,最大的输出位移为0μm-到60μm。
②集成精密定位控制台,型号为PPC-2CR0150。内部集成了压电陶瓷驱动电源输出模块,可以与上位机通信输出稳定的电压以驱动微定位平台的产生位移。控制台内部也集成了位移测量的传感模块,可以精确的测量微定位平台的执行位移并输出。
③具有A/D,D/A转换功能的研华数据采集卡PCI1710U及外部端子。采集卡可以将数字的驱动电压信号通过D/A模块转换成模拟形式输出至集成精密定位控制台。也可以将通过集成精密定位控制台测的实际的位移新型通过A/D模块传送到上位机进行显示,便于检测与调试。
④具有MATLAB实时工作空间RTW测试环境的PC机一台。在PC机上编写输出电压驱动信号,通过集成精密定位控制台施加到压电陶瓷微定位平台,使其产生位移,由集成精密定位控制台内部集成的位移传感器测量,通过采集卡实时采集位移在PC机上进行显示。
⑤隔震台。整个实验过程都将在隔震台上完成,由于压电陶瓷微定位平台的输出位移是微纳米数量级的,非常同意受到外界的干扰如机器振动,人走动等影响,从而对微纳米定位产生严重坏影响,而隔震台可以给实验提供一个相对稳定的外界环境。
建模数据获得方法
实例中用的建模数据分为两类。第一种选幅值为80V的单一幅值的正弦电压信号作为压电陶瓷微定位平台的驱动信号,另一种选幅值按80V-70V-50V-40V-0V幅值递减的正弦电压信号作为压电陶瓷微定位平台的驱动信号。为了体现DPI迟滞模型的动态特性。两种信号分别在输入频率分别为0.1Hz,1Hz,5Hz,10Hz四种情况下,加载到微定位平台上并测得相应频率对应的位移输出。即在不同频率下对压电陶瓷微定位平台的迟滞非线性进行建模,从而对比DPI动态迟滞模型和传统的PI模型对迟滞非线性的动态特性的表述能力。
输入为单一幅值的正弦输入电压信号,频率分别为0.1Hz,1Hz,5Hz,10Hz时测的压电陶瓷微定位平台的实际迟滞特性曲线如图3所示,输入为幅值递减的正弦输入电压信号,频率分别为0.1Hz,1Hz,5Hz,10Hz时,测得的压电陶瓷微定位平台的实际迟滞特性曲线如图4所示。通过图3,图4可以看出随着频率的增大,迟滞环变宽,迟滞现象越严重,充分证明了压电陶瓷迟滞非线性表现是与输入电压信号频率相关的动态特性。为了便于实验进行和建模数据比较,将驱动电压和实际的驱动电压进行归一化到0-1之间。
为了验证DPI动态迟滞模型与传统PI模型对压电陶瓷微定位平台动态迟滞特性的表述能力,两种模型的算子总个数都取60个。DPI算子中非等分阈值运算过程中前半生成算子个数取40个,后半升程的算子个数取20个,DPI算子中的惯性因子在0.1Hz,1Hz,5Hz,10Hz四中不同输入信号频率下分别取0.50,0.75,0.83,0.92。对于DPI动态迟滞模型的权值和传统PI模型的权值的辨识采用结构相同的小波神经网络进行辨识。同时考虑精度要求和模型速度,经过多次实验最终确定小波神经网络为60×10×1的结构。即输入层有60个神经元分别对应着某个时刻60个独立的算子值,隐含层有10个神经元,输出层有一个神经元表示某一时刻网络只对应一个输出。辨识过程中,最大迭代次数设为,误差精度设为,学习因子选,动量因子。
动态迟滞模型与传统PI模型精度对比
在所有实验设备准备齐全,参数已确定情况下,分别用DPI动态迟滞模型与传统的PI模型对压电陶瓷微定位平台的迟滞特性的动态特性建模,并对其精度进行对比。
单一幅值正弦输入电压信号下建模结果对比:
图5-8为四种不同频率的单一幅值正弦输入信号下,压电陶瓷微定位平台实际迟滞曲线与传统PI模型,DPI模型输出迟滞曲线对比图。图中,黑实线为平台实际输入/输出迟滞曲线,·标记线为传统PI模型输入/输出曲线,*标记线为DPI模型输入/输出曲线。可以看出,在低频率使传统PI模型和DPI模型的输出基本与实验输出一致,当频率增大时,传统PI模型拟合度变差,表现出较大的建模误差,而DPI模型保持了较高的拟合程度。图9-12为平台实际位移输出与PI模型,DPI模型输出位移的误差曲线,图中,◇虚线为PI模型的位移误差,·实线为DPI模型误差。对比0.1HZ,1HZ,5HZ,10HZ不同输入频率下的误差曲线可以得出,DPI模型在多频率单一幅值驱动信号下,对比于传统PI模型具有很高的建模精度,且频率越大,DPI模型的优势越明显,这也充分证明了DPI模型是一种可以精确描述压电陶瓷微定位平台迟滞动态特性的动态模型。
让压电陶瓷微定位平台的实际输出位移与DPI动态迟滞模型和传统PI模型的输出分别做差,可以求得在不同频率下的最大的建模误差,根据误差数据可以求取不同频率下两种模型的均方根误差。表1为DPI模型和传统PI模型在输入信号频率分别为0.1HZ,1HZ,5HZ,10HZ时的建模误差和均方误差,对比最大建模误差和均方根误差同样可得出,对比与传统PI模型,DPI模型具有能够描述迟滞非线性动态特性的能力,是一种动态的迟滞模型。
表1 单一幅值信号输入下传统PI模型与DPI模型误差比较
。
幅值递减的正弦输入电压信号建模结果对比:
图13-16为0.1HZ,1HZ,5HZ,10HZ四种不同频率的变幅值正弦输入信号下,压电陶瓷微定位平台实际迟滞曲线与PI模型,DPI模型输出迟滞曲线对比图。黑实线为平台实际输入/输出迟滞曲线,·标记线为传统PI模型输入/输出曲线,*标记线为DPI模型输入/输出曲线。图17-20为四种不同频率实际输出与模型输出误差曲线,◇虚线为PI模型的位移误差,·实线为DPI模型误差。让压电陶瓷微定位平台的实际输出位移与DPI动态迟滞模型和传统PI模型的输出分别做差,同样可求得在在四种不同频率变幅值正弦信号下,DPI模型和传统PI模型的最大建模误差,根据误差数据求取不同频率下两种模型的均方根误差如表2所示。综合分析在幅值递减的正弦电压输入信号下,DPI模型和传统PI模型的输出对比图与误差对比图和表2中数据可知,与传统PI建模结果相比,在变幅值多频率输入信号下,采用DPI模型对具有内环的迟滞特性建模,建模精度较高,可以更好的描述压电陶瓷微定位平台迟滞非线性的动态特性。
表2幅值递减信号输入下传统PI模型与DPI模型误差比较
结论
总结实验可得出,无论是在描述幅值不变具有单一环的迟滞非线性特性,还是描述幅值变化的具有内迟滞环的迟滞非线性特性,DPI模型都可以很好的表述在低频和高频时的压电陶瓷微定位平台的迟滞动态特性,对比与传统的PI模型,具有更高的建模精度。这充分验证了DPI模型是动态迟滞模型,且有广泛的适用性。为压电陶瓷微定位平台在精密定位应用奠定了基础。
Claims (3)
1.一种基于DPI模型的压电陶瓷微定位平台建模方法,其特征在于:DPI模型的建立:
根据一阶微分方程输入输出的类迟滞特性,对算子的表达式进行一阶微分处理构建DPI算子,DPI迟滞算子的最终表达式为:
(3)
式中,为DPI迟滞算子在t时刻的输出值,为DPI算子在t-1时刻的输出值;为DPI算子的阈值,为惯性因子;
DPI模型阈值采取非等分阈值确定表述为:
(4)
式中,为DPI算子输入,为常数,,分别代表前半升程和后半升程第个和第算子的序号;
DPI模型为:
(5)
式中,为DPI模型t时刻的输出值,为迟滞算子的个数,为DPI模型在t时刻第个迟滞算子输出,为DPI模型第个迟滞算子对应的权值。
2.权利要求1所述的基于DPI模型的压电陶瓷微定位平台建模方法,其特征在于:DPI模型权值参数的求取:
将已求得的个DPI迟滞算子值作为小波神经网络的输入即,其中分别代表DPI动态迟滞模型中的第1个,第2个,第m个DPI算子的输出;定义网络输入层的第个输入到隐含层第个神经元间的权值为,隐含层第个神经元到输出层的第个输出之间的权值为,,,分别代表小波神经网络的输入层,隐含层,输出层的神经元序号;t时刻小波神经网络的第个输出为:
(6)
其中,代表隐含层神经元总数,代表隐含层的激励函数,通过公式(7)—(9)求得:
(7)
(8)
(9)
其中,式(9)为Mymorlet小波函数,代表输入层神经元的总数即DPI算子的个数;此处定义为小波函数的伸缩因子,为小波函数的平移因子,则公式(8)的含义为将输入的DPI算子乘以输入层与隐含层的权值后做平移和伸缩变换,在通过公式(7)计算求出激励函数;
在辨识DPI权值参数时,小波神经网络的误差函数采取最小平方误差的形式,表示为::
(10)
式中,代表t时刻实验数据即为理想参考输出,代表着t时刻DPI模型的输出,小波神经网络权值、,节点伸缩参数、平移参数在每次训练过程中对应的修正量,,,分别按如下公式求得:
(11)
(12)
(13)
(14)
依照梯度下降法,小波神经网络t时刻的对应参数调整输出根据以下公式求得:
(15)
(16)
(17)
(18)
式中,为学习系数,为动量因子;
最终求得的DPI动态迟滞模型所有的权值参数,其中代表第1个算子对应的权值,代表第个算子对应的权值,代表第个算子对应的权值。
3.根据权利要求2所述的基于DPI模型的压电陶瓷微定位平台建模方法,其特征在于:小波神经网络辨识过程:
1) 网络训练初始,设定目标精度,设定最大迭代次数运用随机函数给权值、,节点伸缩参数、平移参数赋初值;
2)然后计算网络的所有输出,并与理想的参考信号 进行做差比较求取误差函数;
3)判断误差函数的大小是否小于等于设定的目标精度,如果小于则执行步骤6),否则执行步骤4);
4)依据误差函数,计算新的参数的修正量,,,和对应的参数值、, 、;
5)判断迭代次数是否达到设定的最大值,未到达则转移到步骤2),否则转移到步骤1)重新随机幅值进行训练;
6)停止训练,保存训练得到的参数值;
重复执行步骤1)~6),最终得到满足精度要求的网络权值、,节点伸缩参数、平移参数的取值,此时的小波神经网络结构的功能等同于DPI模型算子对应的权值的功能,即用已经训练好的小波神经网络代替权值矩阵,完成DPI模型权值的辨识。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610583363.1A CN106125574B (zh) | 2016-07-22 | 2016-07-22 | 基于dpi模型的压电陶瓷微定位平台建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610583363.1A CN106125574B (zh) | 2016-07-22 | 2016-07-22 | 基于dpi模型的压电陶瓷微定位平台建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125574A true CN106125574A (zh) | 2016-11-16 |
CN106125574B CN106125574B (zh) | 2018-11-16 |
Family
ID=57290599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610583363.1A Expired - Fee Related CN106125574B (zh) | 2016-07-22 | 2016-07-22 | 基于dpi模型的压电陶瓷微定位平台建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125574B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873380A (zh) * | 2017-04-07 | 2017-06-20 | 哈尔滨理工大学 | 基于pi模型的压电陶瓷模糊pid控制方法 |
CN107367936A (zh) * | 2017-07-31 | 2017-11-21 | 广东工业大学 | 基于os‑elm的压电陶瓷驱动器建模、控制方法及系统 |
CN107578096A (zh) * | 2017-09-21 | 2018-01-12 | 胡明建 | 一种压频式择端人工神经元的设计方法 |
CN109670215A (zh) * | 2018-11-28 | 2019-04-23 | 上海工程技术大学 | 基于遗传算法的悬臂梁振动模型参数辨识方法及装置 |
CN110108443A (zh) * | 2019-05-05 | 2019-08-09 | 大连理工大学 | 一种基于神经网络的压电陶瓷驱动器输出控制方法 |
CN110245430A (zh) * | 2019-06-18 | 2019-09-17 | 吉林大学 | 改进Bouc-Wen模型迟滞建模方法 |
CN110470921A (zh) * | 2019-08-14 | 2019-11-19 | 上海卫星工程研究所 | 压电作动器输出力迟滞效应测试系统及测试方法 |
CN111142404A (zh) * | 2019-12-17 | 2020-05-12 | 吉林大学 | 一种基于压电陶瓷驱动的微定位平台及其建模与控制方法 |
CN111368400A (zh) * | 2020-02-17 | 2020-07-03 | 华南理工大学 | 基于pso算法的压电微驱动变频定位平台建模辨识方法 |
CN111459021A (zh) * | 2020-04-27 | 2020-07-28 | 沈阳建筑大学 | 一种基于分段pi模型的纳米定位平台补偿控制方法 |
CN111487862A (zh) * | 2020-04-27 | 2020-08-04 | 沈阳建筑大学 | 一种基于内环分离pi模型的压电驱动器迟滞补偿方法 |
CN111897211A (zh) * | 2020-05-31 | 2020-11-06 | 吉林大学 | 考虑约束条件的压电陶瓷微定位平台轨迹跟踪控制方法 |
CN111897210A (zh) * | 2020-05-24 | 2020-11-06 | 吉林大学 | 压电陶瓷微定位平台建模方法 |
CN111914981A (zh) * | 2020-05-31 | 2020-11-10 | 吉林大学 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
CN111931411A (zh) * | 2020-05-25 | 2020-11-13 | 吉林大学 | 压电驱动微定位平台Duhem动态迟滞建模方法 |
CN111930008A (zh) * | 2020-06-04 | 2020-11-13 | 吉林大学 | 基于数据驱动控制的压电微定位平台轨迹跟踪控制方法 |
CN113067497A (zh) * | 2021-03-26 | 2021-07-02 | 合肥工业大学 | 基于压电陶瓷驱动器的迟滞分割建模与补偿方法 |
CN113487016A (zh) * | 2021-07-07 | 2021-10-08 | 华中科技大学鄂州工业技术研究院 | 基于神经网络的扫描器件位移控制方法、设备及存储介质 |
CN115600480A (zh) * | 2022-06-13 | 2023-01-13 | 哈尔滨工业大学(Cn) | 一种针对压电换能器的全局线性化频散迟滞建模方法、装置、计算机设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1794116A (zh) * | 2005-12-22 | 2006-06-28 | 桂林电子工业学院 | 基于神经网络的迟滞特性建模方法 |
CN102280572A (zh) * | 2011-04-15 | 2011-12-14 | 重庆大学 | 压电陶瓷执行器迟滞特性的复合线性化控制方法及其实现电路 |
CN103853046A (zh) * | 2014-02-14 | 2014-06-11 | 广东工业大学 | 一种压电陶瓷驱动器的自适应学习控制方法 |
CN104678765A (zh) * | 2015-01-28 | 2015-06-03 | 浙江理工大学 | 压电陶瓷执行器迟滞模型及其控制方法 |
-
2016
- 2016-07-22 CN CN201610583363.1A patent/CN106125574B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1794116A (zh) * | 2005-12-22 | 2006-06-28 | 桂林电子工业学院 | 基于神经网络的迟滞特性建模方法 |
CN102280572A (zh) * | 2011-04-15 | 2011-12-14 | 重庆大学 | 压电陶瓷执行器迟滞特性的复合线性化控制方法及其实现电路 |
CN103853046A (zh) * | 2014-02-14 | 2014-06-11 | 广东工业大学 | 一种压电陶瓷驱动器的自适应学习控制方法 |
CN104678765A (zh) * | 2015-01-28 | 2015-06-03 | 浙江理工大学 | 压电陶瓷执行器迟滞模型及其控制方法 |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873380A (zh) * | 2017-04-07 | 2017-06-20 | 哈尔滨理工大学 | 基于pi模型的压电陶瓷模糊pid控制方法 |
CN107367936A (zh) * | 2017-07-31 | 2017-11-21 | 广东工业大学 | 基于os‑elm的压电陶瓷驱动器建模、控制方法及系统 |
CN107578096A (zh) * | 2017-09-21 | 2018-01-12 | 胡明建 | 一种压频式择端人工神经元的设计方法 |
CN109670215A (zh) * | 2018-11-28 | 2019-04-23 | 上海工程技术大学 | 基于遗传算法的悬臂梁振动模型参数辨识方法及装置 |
CN110108443A (zh) * | 2019-05-05 | 2019-08-09 | 大连理工大学 | 一种基于神经网络的压电陶瓷驱动器输出控制方法 |
CN110108443B (zh) * | 2019-05-05 | 2020-04-07 | 大连理工大学 | 一种基于神经网络的压电陶瓷驱动器输出控制方法 |
CN110245430A (zh) * | 2019-06-18 | 2019-09-17 | 吉林大学 | 改进Bouc-Wen模型迟滞建模方法 |
CN110470921A (zh) * | 2019-08-14 | 2019-11-19 | 上海卫星工程研究所 | 压电作动器输出力迟滞效应测试系统及测试方法 |
CN111142404A (zh) * | 2019-12-17 | 2020-05-12 | 吉林大学 | 一种基于压电陶瓷驱动的微定位平台及其建模与控制方法 |
CN111368400A (zh) * | 2020-02-17 | 2020-07-03 | 华南理工大学 | 基于pso算法的压电微驱动变频定位平台建模辨识方法 |
CN111459021A (zh) * | 2020-04-27 | 2020-07-28 | 沈阳建筑大学 | 一种基于分段pi模型的纳米定位平台补偿控制方法 |
CN111487862A (zh) * | 2020-04-27 | 2020-08-04 | 沈阳建筑大学 | 一种基于内环分离pi模型的压电驱动器迟滞补偿方法 |
CN111487862B (zh) * | 2020-04-27 | 2023-06-16 | 沈阳建筑大学 | 一种基于内环分离pi模型的压电驱动器迟滞补偿方法 |
CN111897210B (zh) * | 2020-05-24 | 2022-12-06 | 吉林大学 | 压电陶瓷微定位平台建模方法 |
CN111897210A (zh) * | 2020-05-24 | 2020-11-06 | 吉林大学 | 压电陶瓷微定位平台建模方法 |
CN111931411A (zh) * | 2020-05-25 | 2020-11-13 | 吉林大学 | 压电驱动微定位平台Duhem动态迟滞建模方法 |
CN111897211A (zh) * | 2020-05-31 | 2020-11-06 | 吉林大学 | 考虑约束条件的压电陶瓷微定位平台轨迹跟踪控制方法 |
CN111897211B (zh) * | 2020-05-31 | 2022-09-27 | 吉林大学 | 考虑约束条件的压电陶瓷微定位平台轨迹跟踪控制方法 |
CN111914981B (zh) * | 2020-05-31 | 2022-11-08 | 吉林大学 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
CN111914981A (zh) * | 2020-05-31 | 2020-11-10 | 吉林大学 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
CN111930008A (zh) * | 2020-06-04 | 2020-11-13 | 吉林大学 | 基于数据驱动控制的压电微定位平台轨迹跟踪控制方法 |
CN113067497A (zh) * | 2021-03-26 | 2021-07-02 | 合肥工业大学 | 基于压电陶瓷驱动器的迟滞分割建模与补偿方法 |
CN113067497B (zh) * | 2021-03-26 | 2022-12-06 | 合肥工业大学 | 基于压电陶瓷驱动器的迟滞分割建模与补偿方法 |
CN113487016A (zh) * | 2021-07-07 | 2021-10-08 | 华中科技大学鄂州工业技术研究院 | 基于神经网络的扫描器件位移控制方法、设备及存储介质 |
CN115600480A (zh) * | 2022-06-13 | 2023-01-13 | 哈尔滨工业大学(Cn) | 一种针对压电换能器的全局线性化频散迟滞建模方法、装置、计算机设备及存储介质 |
CN115600480B (zh) * | 2022-06-13 | 2023-06-27 | 哈尔滨工业大学 | 一种针对压电换能器的全局线性化频散迟滞建模方法、装置、计算机设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106125574B (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106125574A (zh) | 基于dpi模型的压电陶瓷微定位平台建模方法 | |
CN110245430B (zh) | 改进Bouc-Wen模型迟滞建模方法 | |
CN104991997B (zh) | 自适应差分进化算法优化的广义率相关p-i迟滞建模方法 | |
WO2003087746A2 (en) | Methods and systems for modeling material behavior | |
CN109948207A (zh) | 一种航空发动机高压转子装配误差预测方法 | |
CN109885916B (zh) | 一种基于lssvm的混合试验在线模型更新方法 | |
CN110007645A (zh) | 一种基于动力学及深度神经网络的进给系统混合建模方法 | |
CN110631792A (zh) | 基于卷积神经网络的抗震混合试验模型更新方法 | |
CN105259331A (zh) | 一种节理岩体单轴强度预测方法 | |
CN106568647B (zh) | 一种基于神经网络的混凝土强度预测方法 | |
CN115270239A (zh) | 基于动力特性和智能算法响应面法的桥梁可靠性预测方法 | |
Mohammadi | Groundwater table estimation using MODFLOW and artificial neural networks | |
Goktepe et al. | Shear strength estimation of plastic clays with statistical and neural approaches | |
Anton et al. | Identification of material parameters from full-field displacement data using physics-informed neural networks | |
CN109579967A (zh) | 智能动态称重方法与系统 | |
CN107314938B (zh) | 焊核区域材料塑性参数反演识别的实现方法 | |
Joghataie et al. | Designing high-precision fast nonlinear dam neuro-modelers and comparison with finite-element analysis | |
KR102480382B1 (ko) | 인공지능 기반 풍하중 산정 시스템 | |
CN114491730B (zh) | 一种高速铁路路基结构动力安定分析迭代方法及装置 | |
CN113973403B (zh) | 基于结构离散测点拓扑的温致应变场重分布智能感知方法 | |
CN114357878A (zh) | 一种水陆两栖飞机典型截面的着水面载荷预测方法 | |
Ahn et al. | Hysteresis modeling of magneto-rheological (MR) fluid damper by self tuning fuzzy control | |
Yang et al. | Statistical pattern recognition for structural health monitoring using esn feature extraction method | |
Rezania | Evolutionary polynomial regression based constitutive modelling and incorporation in finite element analysis | |
Zong-Yi et al. | Modelling and identification of electrohydraulic system and its application |
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 |
Granted publication date: 20181116 |