基于贝叶斯反演算建模的微纳米探头动态特性补偿方法

摘要

本发明公开了一种基于贝叶斯反演算建模的微纳米探头动态特性补偿方法，其特征是通过微纳米探头系统的阶跃响应实测数据，采用自适应模型辨识法得到探头系统的传递函数模型；根据系统传递函数模型参数推导确定贝叶斯反演算补偿模型的观测回归矩阵和状态转移矩阵，获得动态补偿模型的状态方程和观测方程；采用参考先验分析法确定状态向量的初始先验分布，利用贝叶斯递推算法估计状态向量实时分布参数，得到补偿后的实时输入待估计量。本发明用于实现微纳米探头系统的动态特性补偿，能够快速有效地减小探头阶跃输出的超调量，缩短探头响应时间，以低成本改善微纳米探头系统的动态特性。

说明书

技术领域

本发明涉及微纳米探头动态特性补偿技术领域，更具体地说是一种基于贝叶斯反演算建模的微纳米探头动态特性补偿方法。

背景技术

动态特性不理想是微纳米探头系统动态测量误差的重要来源，严重制约探头系统测量速度和精度的提升。微纳米探头系统结构复杂、样式多变，难以通过理论建模或结构改进实现动态特性提升。为了在提升系统测量速度同时保证测量精度，对其进行动态特性补偿是较为经济有效的办法，是在其测量输出端串联动态补偿器，要求动态补偿器模型与系统的动态特性模型互逆，从而实现动态特性补偿。

动态补偿器的设计依据主要有基于输出数据和基于系统模型两类。前者虽然无需辨识系统模型，但输出数据的误差会对动态补偿器的设计精度和稳健性产生较大的不良影响。因此，基于系统模型设计传感器的动态补偿器仍然是目前最常用的办法。为了获得动态补偿器，常用的补偿方法有：零极点配置法、最小二乘法和神经网络法等。其中，零极点配置法基于传感器的动态模型，能够通过输入输出数据准确地获知传感器的参数，但是对于动态特性复杂的传感器，很难获得高精度的参数模型。最小二乘法易于实现、应用广泛，但在辨识过程中存在不确定度干扰、数据饱和及参数偏差等问题。神经网络法收敛速度慢、容易陷入局部极值，而且需要训练网络，程序运行时间较长。现有技术中微纳米探头系统动态补偿方法普遍存在着实时性差和建模精度低的问题，不适用于对微纳米探头系统进行在线动态特性补偿。

发明内容

本发明是为避免上述现有技术所存在的不足，提供一种基于贝叶斯反演算建模的微纳米探头动态特性补偿方法，以实现快速、准确地建立探头系统的动态补偿模型，经济、有效地改善系统动态测试性能的发明目的。

为实现发明目的，本发明采用如下技术方案：

本发明基于贝叶斯反演算建模的微纳米探头动态特性补偿方法的特点是按如下步骤进行：

步骤1、由动态阶跃响应实验获得微纳米接触式探头系统的阶跃响应实测数据，利用所述阶跃响应实测数据，根据自适应模型参数辨识法建立探头系统的输入输出传递函数模型，获得由式(1)表征的n阶线性常系数微分方程：

式(1)中：

y

u

步骤2、令x

式(2)中：

y

y

以x

以

式(3)中：

u

x

以x

以x

由式(3)分别获得式(4)和式(5)：

y

式(4)表征测试序列的观测方程；式(5)表征状态序列方程组；

步骤3、采用线性逼近建模方法通过拟合获得未知的待估计量参数规律，是将u

观测方程：u

状态方程：

式(6)和式(7)中：

μ

μ

v

v

步骤4、令θ

观测方程：y

状态方程：

式(8)中：

以F

ν

式(9)中：以G表示状态转移矩阵，即：

ω

对式(8)和式(9)进行简化获得由式(10)表征的贝叶斯反演算补偿模型：

步骤5、采用参考先验分析法和贝叶斯递推算法确定θ

本发明基于贝叶斯反演算建模的微纳米探头动态特性补偿方法的特点也在于：

以N表征数据量，将k时刻的取值范围0到N-1划分为第一阶段和第二阶段，第一阶段k取0,1,…,n+2，第二阶段k取n+3,n+4,…,N-1；将所述步骤5中为获得状态参数μ

第一阶段：k＝0,1,…,n+2，根据探头系统最初测得的n+3个阶跃响应测试序列，按如下方式采用参考先验分析法确定状态向量θ

测量数据误差服从正态分布，状态噪声方差W

p(θ

式(11)中：

以p(θ

以V为观测方差，以V

以D

由式(12)表征观测数据的似然函数：

p(y

式(12)中：

以p(y

以V

以exp表示以无理数e为底的指数函数；

以D

结合式(12)并根据贝叶斯公式和探头系统模型，分别获得状态向量θ

式(13)和(14)中：

以p(θ

以p(θ

以H

H

γ

各初始值为零，即：H

由此获得由式(15)表征的θ

式(15)中：

(θ

式(16)中：

Γ[n

n

第二阶段：k＝n+3,n+4,…,N-1，采用折扣因子法确定由式(17)表征的递推算法中k+1时刻W

W

式(17)中：

C

ρ为折扣因子，折扣因子ρ是以观测数据预测方差最小为准则通过最优化搜索以确定；

根据探头系统输出y

p(θ

由式(18)获得由式(19)表征的θ

θ

与已有技术相比，本发明有益效果体现在：

1、本发明在传递函数或微分方程已知的基础上直接建立贝叶斯反演算动态补偿模型，实现了快速、准确地探头系统动态特性补偿。

2、本发明基于贝叶斯统计原理，融合历史先验信息和当前样本信息进行参数递推估计，不仅所需样本少，而且还保障了动态补偿输出结果的可靠性与稳定性。

3、本发明的动态补偿方法可由软件编程实现，无需改变探头系统的原有结构，只要将该方法串联在系统的测量数据处理环节中，就能经济、有效地提高探头系统的动态测试性能。

附图说明

图1是补偿前的微纳米接触式探头系统阶跃响应实测数据图。

图2是基于贝叶斯反演算补偿探头系统阶跃响应实测数据的效果图。

具体实施方式

本实施例关于基于贝叶斯反演算建模的微纳米探头动态特性补偿方法，其微纳米探头是指微纳米接触式探头系统，其动态特性补偿是指针对微纳米接触式探头系统的实测输出数据进行动态补偿，以满足快速测量的需求；其动态特性补偿方法是按如下步骤通过动态阶跃响应实验、模型参数辨识以及动态补偿器构造，实现微纳米探头动态特性补偿：

步骤1、由动态阶跃响应实验获得微纳米接触式探头系统的阶跃响应实测数据，利用阶跃响应实测数据，根据自适应模型参数辨识法建立探头系统的输入输出传递函数模型，获得由式(1)表征的n阶线性常系数微分方程：

式(1)中：

y

u

步骤2、为了得到贝叶斯反演算补偿模型，需要先把线性常系数微分方程转换成状态空间方程，然后再进行近似地离散化处理，获得测试序列的观测方程和状态序列方程组。因此，首先，假设x

式(2)中：

y

y

以x

考虑到探头系统在动态响应或扫描测量过程中，采样频率较高，采样时间间隔较短，以Δt表示探头系统在动态响应中的采样时间间隔，设定Δt→0，由前后相邻的k+1时刻和k时刻的状态值近似估计获得t时刻的一阶导数，将式(2)转化为由式(3)表征的方程组：

式(3)中：

u

x

以x

以x

由式(3)分别获得式(4)和式(5)：

y

式(4)表征测试序列的观测方程；式(5)表征状态序列方程组；

步骤3、贝叶斯反演算补偿的主要目的是获得输入量u

观测方程：u

状态方程：

式(6)和式(7)中：

μ

μ

v

v

步骤4、根据式(4)表征的测试序列的观测方程、式(5)表征的状态序列方程组以及输入量u

令θ

观测方程：y

状态方程：

式(8)中：

以F

ν

式(9)中：

以G表示状态转移矩阵，即：

ω

对式(8)和式(9)进行简化获得由式(10)表征的贝叶斯反演算补偿模型：

步骤5、采用参考先验分析法和贝叶斯递推算法确定θ

具体实施中，以N表征数据量，将k时刻的取值范围0到N-1划分为第一阶段和第二阶段，第一阶段k取0,1,…,n+2，第二阶段k取n+3,n+4,…,N-1；将步骤5中为获得状态参数μ

第一阶段：k＝0,1,…,n+2，根据探头系统最初测得的n+3个阶跃响应测试序列，按如下方式采用参考先验分析法确定状态向量θ

在参考先验分析中，探头系统共有n+3个待估参数，即θ

测量数据误差服从正态分布，状态噪声方差W

p(θ

式(11)中：

以p(θ

以V为观测方差，以V

以D

由式(12)表征观测数据的似然函数：

p(y

式(12)中：

以p(y

以V

以exp表示以无理数e为底的指数函数；

以D

结合式(12)并根据贝叶斯公式和探头系统模型，分别获得状态向量θ

式(13)和(14)中：

以p(θ

以p(θ

以H

H

γ

各初始值为零，即：H

由此获得由式(15)表征的θ

式(15)中：

(θ

T

式(16)中：

Γ[n

n

第二阶段：k＝n+3,n+4,…,N-1，采用折扣因子法确定由式(17)表征的递推算法中k+1时刻W

W

式(17)中：

C

ρ为折扣因子，折扣因子ρ是以观测数据预测方差最小为准则通过最优化搜索以确定；

按照参考先验分析法，根据探头系统输出y

p(θ

由式(18)获得由式(19)表征的θ

θ

利用本发明方法按如下过程针对公开号为CN104457613A的发明专利说明书中提出的一种三维微纳米接触触发式探头进行动态特性补偿，其结构简单、装调方便，分辨力达到1nm，重复性为8.2nm，允许触碰范围为±10μm。

首先对探头系统进行阶跃响应实验测试，采样间隔为2.5×10

根据阶跃响应实测数据，采用自适应模型辨识法确定探头系统为典型的二阶系统，即n取值为2，探头系统输入输出传递函数模型G(s)由式(21)所表征：

式(21)中，s为复频率。

由式(21)计算得到探头系统的阻尼比ζ为0.0099、固有频率ω

令x

将式(6)和式(7)代入到式(23)中，可得由式(24)表征的贝叶斯反演算模型：

第一阶段：设定k＝0,1,2,3,4，假定测量数据误差服从正态分布，状态噪声方差W

在参考先验分析中，二阶探头系统共有5个待估参数，即θ

按照参考先验分析法，根据贝叶斯公式、探头系统模型以及探头前5个动态测量数据，分别获得状态向量θ

第二阶段：设定k＝5,6,…,15999，采用折扣因子法确定由式(17)表征的递推算法中k+1时刻W

根据探头系统输出y

微纳米探头系统动态特性补偿效果展示：

微纳米接触式探头系统的阶跃响应实测数据补偿后效果如图2所示，补偿前后探头系统的动特性参数如表1所示。可见：补偿后探头系统阶跃响应的超调量从99.15％降低到40.05％，响应时间由2.4825s缩短到0.1720s。从应用的角度出发，补偿后探头阶跃输出的响应时间很短，在实际测量过程中，一般会在输出信号稳定之后开始采集有效数据，所以在响应时间之前的输出信号的超调量不会对最终结果产生影响。因此，该贝叶斯反演算补偿结果满足动态测量的需求，探头系统的动态特性得到了较大改善。

表1贝叶斯反演算补偿前后效果比较

尽管已经示出和描述了本发明的实施例，对于本领域的技术人员，可以在不脱离本发明的原理和精神的情况下对实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。