一种多干扰因素耦合的光栅尺误差补偿方法
【专利摘要】本发明公开了一种多干扰因素耦合的光栅尺误差补偿方法，包括以下步骤：在光栅尺测量过程中采用多个传感器测量获得多种干扰因素的作用强度值，进而与误差补偿数据库进行匹配，获得该组作用强度值对应的最优误差补偿量；采用该最优误差补偿量对光栅尺系统进行补偿。本发明通过测量多种干扰因素的作用强度值后，与误差补偿数据库进行匹配获得最优误差补偿量后对光栅尺系统进行补偿，操作简单、成本较低，而且补偿效果好，可实现对光栅尺系统的有效补偿，可广泛应用于光栅尺测量行业中。
【专利说明】
一种多干扰因素耦合的光栅尺误差补偿方法
技术领域
[0001] 本发明涉及光栅的误差测量领域，特别是涉及一种多干扰因素耦合的光栅尺误差 补偿方法。
【背景技术】
[0002] 名词解释：
[0003] 经验模式分解:Empirical Mode Decomposition，简称EMD，一种信号分析处理的 算法，算法思想为:将复杂信号分解为有限个本征模函数（Intrinsic Mode Function，简称 頂F)分量，所分解出来的各頂F分量包含了原信号的不同时间尺度的局部特征信号；
[0004] IMF: Intrinsic Mode Function，本征模态函数，信号经过EMD分解后所产生的满 足一定条件的信号，通常EMD分解对信号进行分解后，会产生若干个頂F分量和一个残余分 量；
[0005] HMS:Hilbert marginal spectrum，希尔伯特边缘谱，一种频谱图。
[0006] 光栅尺作为一种高精度的位置传感器己经广泛应用于各种位移测量系统中，是全 闭环数控机床最为常用的位置检测设备。光栅尺的测量精度直接影响数控机床的加工精 度。光栅尺的光电系统误差、振动误差、安装误差和温度误差是影响其测量精度的关键因 素，误差补偿法是一种既有效又经济的提高光栅尺测量精度的手段。
[0007] 从误差源的角度分析，光栅尺的误差因素主要来源于四个方面：1)光电系统误差， 光栅尺光电系统误差主要来源于莫尔条纹信号质量及细分，同时在绝对式光栅尺中，还需 另外考虑绝对码道编码图像的图像信号处理方法和相关的像素细分技术;2)制造误差，光 栅尺制造误差包括标尺光栅刻划误差及机械安装误差等；3)运行过程中产生的误差，主要 包括运行过程中的温度误差、振动误差、尺身污染等;4)电子器件带来的误差，如电子器件 在光栅尺采集、位置解码中的时序延时、电子元器件的非线性等。
[0008] 自光栅尺得到应用以来，国内外相继展开了光栅尺误差补偿的研究，从误差补偿 的过程来说，光栅尺误差补偿的研究大致分为单误差因素抑制及补偿以及多干扰因素耦合 下的多误差元素综合补偿两个方面。单误差因素抑制及补偿主要是针对光栅尺误差的多个 主要来源分别进行抑制和补偿，采用不同的补偿方法对多个来源的光栅尺误差进行抑制， 从而整体上达到误差补偿的目的。而在进行单误差因素抑制补偿的基础上，研究发现多个 误差元素之间存在耦合关系，有些误差元素在实际加工中会相互抵消，有些误差元素会相 互叠加，这样单一误差元素的补偿会产生过补偿或欠补偿的问题。因此，出现了多误差元素 综合补偿的补偿方式，通过多个误差元素之间的耦合关系进行综合补偿。但是，目前技术 中，多误差元素综合补偿的补偿方法主要是通过对单项误差近似线性叠加的方式进行补 偿，由于多干扰因素在测量误差中的增减性方向不一致，通过简单的线性叠加进行补偿的 方式并不准确，难以进行有效的补偿。而如果通过分析每种误差元素对总的测量误差的影 响后，将所有的误差元素进行矢量累加来获得实际误差再进行补偿，这种方式虽然理论上 可以有效地进行误差补偿，但是操作繁琐、困难、操作效率低，而且需要耗费较多的人力物 力。因此，总的来说，目前技术中，难以对光栅尺测量过程中的多干扰因素造成的误差进行 耦合补偿。
【发明内容】
[0009] 为了解决上述的技术问题，本发明的目的是提供一种多干扰因素耦合的光栅尺误 差补偿方法。
[0010] 本发明解决其技术问题所采用的技术方案是：
[0011] -种多干扰因素耦合的光栅尺误差补偿方法，包括以下步骤：
[0012] 在光栅尺测量过程中采用多个传感器测量获得多种干扰因素的作用强度值，进而 与误差补偿数据库进行匹配，获得该组作用强度值对应的最优误差补偿量；
[0013] 采用该最优误差补偿量对光栅尺系统进行补偿。
[0014] 进一步，还包括以下步骤：
[0015] 依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每组作用强度值 下的多组误差数据；
[0016] 基于经验模式分解算法计算在每组作用强度值下的最优误差补偿量，进而对应建 立误差补偿数据库。
[0017] 进一步，所述依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每 组作用强度值下的多组误差数据的步骤，其具体为：
[0018] 通过多个误差调节器依次调节多种干扰因素的作用强度，并采用传感器测量获得 调节过程中的每组作用强度值，同时在每组作用强度值下，采用光栅尺系统和激光干涉仪 对待测对象进行多次测量后计算获得多组误差数据。
[0019] 进一步，所述基于经验模式分解算法计算在每组作用强度值下的最优误差补偿 量，进而对应建立误差补偿数据库的步骤，包括步骤：
[0020] A1、将每组作用强度下的多组误差数据平均分成两个集合，分别为分析数据集和 测试数据集；
[0021] A2、基于经验模式分解算法分别对分析数据集的每组误差数据进行分解计算后， 获得每组误差数据所对应的趋势误差分量；
[0022] A3、分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进 行补偿效果统计，计算获得每个趋势误差分量对应的多个补偿余量的平均值；
[0023] A4、将平均值最小的多个补偿余量对应的趋势误差分量作为该组作用强度下的最 优误差补偿量。
[0024] 进一步，所述步骤A2,包括：
[0025] A21、针对分析数据集的每组误差数据，采用经验模式分解算法将其分解成多个 IMF分量和一个残余分量；
[0026] A22、分别对每个IMF分量进行希尔伯特变换，并计算获得每个IMF分量的希尔伯特 谱和希尔伯特边缘谱；
[0027] A23、根据预设筛选条件对所有IMF分量进行分析后，筛选出符合条件的IMF分量；
[0028] A24、将残余分量和所有符合条件的IMF分量累加后作为误差信号的趋势误差分 量。
[0029] 进一步，所述步骤A23,包括：
[0030] A231、将每个頂F分量的希尔伯特边缘谱按照频率顺序均分为低频区域、中频区域 和尚频区域；
[0031] A232、分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而 分别计算每个区域的归一化权重因子；
[0032] A233、筛选出低频区域的归一化权重因子大于预设阈值的MF分量作为符合条件 的IMF分量。
[0033] 进一步，所述步骤A232,其具体为：
[0034]分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而根据下 式分别计算每个区域的归一化权重因子：
[0036]上式中，Mj, i、Mj, 2和Mj,3依次表示第j个MF分量的低频区域、中频区域和高频区域 的归一化权重因子，和叫,3分别表示第j个頂F分量的低频区区域、中频区域和高频 区域的权重因子，j为自然数。
[0037] 进一步，所述步骤A233中所述预设阈值为0.5。
[0038]进一步，所述步骤A3,其具体为：
[0039]分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进行补 偿效果统计，根据下式计算获得每个趋势误差分量对应的多个补偿余量的平均值ε1:
[0041] ￡l表示第i组误差数据的趋势误差分量对应的多个补偿余量的平均值，e2j表示测 试数据集中的第j组误差数据，su表示分析数据集中的第i组误差数据对应的趋势误差分 量，i、j均为自然数，η表示测试数据集的误差数据的总组数。
[0042]进一步，所述步骤Α21，其具体为：
[0043]针对分析数据集的每组误差数据，根据下式，采用经验模式分解算法将其分解成 多个IMF分量和一个残余分量和一个残余分量：
[0045]上式中，X(t)表示误差数据，i为自然数，m表示分解的MF分量的总个数，Cl(t)表 示经验模式分解后的第i个IMF分量，rm( t)表示残余分量。
[0046] 进一步，所述步骤A22,其具体为：
[0047] 分别对每个IMF分量进行希尔伯特变换后，计算获得每个IMF分量的瞬时幅值、瞬 时相位和瞬时频率，进而计算获得每个IMF分量的希尔伯特边缘谱。
[0048] 进一步，所述步骤A22,包括：
[0049] A221、根据下式，分别对每个IMF分量进行希尔伯特变换：
[0051 ] 上式中，P表不奇异积分的柯西主值，c(t)表不IMF分量，CH(t)表不IMF分量c(t)的 希尔伯特变换；
[0052] A222、根据下式，计算获得每个IMF分量的瞬时幅值、瞬时相位和瞬时频率：
[0054] 上式中，a(t)表示MF分量c(t)的瞬时幅值，Φ (t)表示MF分量c(t)的瞬时相位， ω (t)表示IMF分量c(t)的瞬时频率；
[0055] A223、根据下式计算获得每个IMF分量的希尔伯特边缘谱：
[0057] 上式中，Η( ω，t)表示IMF分量c(t)的希尔伯特谱，h( ω )表示IMF分量c(t)的希尔 伯特边缘谱，T表示积分周期。
[0058]本发明的有益效果是:本发明的一种多干扰因素耦合的光栅尺误差补偿方法，包 括以下步骤:在光栅尺测量过程中采用多个传感器测量获得多种干扰因素的作用强度值， 进而与误差补偿数据库进行匹配，获得该组作用强度值对应的最优误差补偿量;采用该最 优误差补偿量对光栅尺系统进行补偿。本方法通过测量多种干扰因素的作用强度值后，与 误差补偿数据库进行匹配获得最优误差补偿量后对光栅尺系统进行补偿，操作简单、成本 较低，而且补偿效果好，可实现对光栅尺系统的有效补偿。
【附图说明】
[0059] 下面结合附图和实施例对本发明作进一步说明。
[0060] 图1是本发明的一种多干扰因素耦合的光栅尺误差补偿方法中如何建立误差补偿 数据库的流程图。
【具体实施方式】
[0061] 本发明提供了一种多干扰因素耦合的光栅尺误差补偿方法，包括以下步骤：
[0062] 在光栅尺测量过程中采用多个传感器测量获得多种干扰因素的作用强度值，进而 与误差补偿数据库进行匹配，获得该组作用强度值对应的最优误差补偿量；
[0063] 采用该最优误差补偿量对光栅尺系统进行补偿。
[0064] 进一步作为优选的实施方式，还包括以下步骤：
[0065] 依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每组作用强度值 下的多组误差数据；
[0066] 基于经验模式分解算法计算在每组作用强度值下的最优误差补偿量，进而对应建 立误差补偿数据库。
[0067] 进一步作为优选的实施方式，所述依次调节多种干扰因素的作用强度，测量获得 在多种干扰因素的每组作用强度值下的多组误差数据的步骤，其具体为：
[0068] 通过多个误差调节器依次调节多种干扰因素的作用强度，并采用传感器测量获得 调节过程中的每组作用强度值，同时在每组作用强度值下，采用光栅尺系统和激光干涉仪 对待测对象进行多次测量后计算获得多组误差数据。误差调节器用于调节各干扰因素的作 用强度，即调节各干扰因素所带来的误差大小；传感器用于测量各干扰因素所带来的作用 强度值，即误差值。
[0069] 进一步作为优选的实施方式，参照图1，所述基于经验模式分解算法计算在每组作 用强度值下的最优误差补偿量，进而对应建立误差补偿数据库的步骤，包括步骤：
[0070] A1、将每组作用强度下的多组误差数据平均分成两个集合，分别为分析数据集和 测试数据集；
[0071] A2、基于经验模式分解算法分别对分析数据集的每组误差数据进行分解计算后， 获得每组误差数据所对应的趋势误差分量；
[0072] A3、分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进 行补偿效果统计，计算获得每个趋势误差分量对应的多个补偿余量的平均值；
[0073] A4、将平均值最小的多个补偿余量对应的趋势误差分量作为该组作用强度下的最 优误差补偿量。
[0074] 进一步作为优选的实施方式，所述步骤A2,包括：
[0075] A21、针对分析数据集的每组误差数据，采用经验模式分解算法将其分解成多个 IMF分量和一个残余分量；
[0076] A22、分别对每个IMF分量进行希尔伯特变换，并计算获得每个IMF分量的希尔伯特 谱和希尔伯特边缘谱；
[0077] A23、根据预设筛选条件对所有IMF分量进行分析后，筛选出符合条件的IMF分量；
[0078] A24、将残余分量和所有符合条件的IMF分量累加后作为误差信号的趋势误差分 量。
[0079] 进一步作为优选的实施方式，所述步骤A23,包括：
[0080] A231、将每个頂F分量的希尔伯特边缘谱按照频率顺序均分为低频区域、中频区域 和尚频区域；
[0081] A232、分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而 分别计算每个区域的归一化权重因子；
[0082] A233、筛选出低频区域的归一化权重因子大于预设阈值的MF分量作为符合条件 的IMF分量。
[0083]进一步作为优选的实施方式，所述步骤A232，其具体为：
[0084]分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而根据下 式分别计算每个区域的归一化权重因子：
[0086]上式中，Mj, i、Mj, 2和Mj,3依次表示第j个MF分量的低频区域、中频区域和高频区域 的归一化权重因子，和叫,3分别表示第j个頂F分量的低频区区域、中频区域和高频 区域的权重因子，j为自然数。
[0087] 进一步作为优选的实施方式，所述步骤A233中所述预设阈值为0.5。
[0088] 进一步作为优选的实施方式，所述步骤A3，其具体为：
[0089] 分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进行补 偿效果统计，根据下式计算获得每个趋势误差分量对应的多个补偿余量的平均值ε 1:
[0091] ει表示第i组误差数据的趋势误差分量对应的多个补偿余量的平均值，e2j表示测 试数据集中的第j组误差数据，SU表示分析数据集中的第i组误差数据对应的趋势误差分 量，i、j均为自然数，η表示测试数据集的误差数据的总组数。
[0092]进一步作为优选的实施方式，所述步骤Α21，其具体为：
[0093]针对分析数据集的每组误差数据，根据下式，采用经验模式分解算法将其分解成 多个IMF分量和一个残余分量和一个残余分量：
[0095]上式中，X(t)表示误差数据，i为自然数，m表示分解的頂F分量的总个数，Cl(t)表 示经验模式分解后的第i个IMF分量，rm( t)表示残余分量。
[0096]进一步作为优选的实施方式，所述步骤A22，其具体为：
[0097] 分别对每个IMF分量进行希尔伯特变换后，计算获得每个IMF分量的瞬时幅值、瞬 时相位和瞬时频率，进而计算获得每个IMF分量的希尔伯特边缘谱。
[0098] 进一步作为优选的实施方式，所述步骤A22,包括：
[0099] A221、根据下式，分别对每个IMF分量进行希尔伯特变换：
[0?0? ] 上式中，P表不奇异积分的柯西主值，c(t)表不IMF分量，CH(t)表不IMF分量c(t)的 希尔伯特变换；
[0102] A222、根据下式，计算获得每个IMF分量的瞬时幅值、瞬时相位和瞬时频率：
[0104] 上式中，a(t)表示MF分量c(t)的瞬时幅值，Φ (t)表示MF分量c(t)的瞬时相位， ω (t)表示IMF分量c(t)的瞬时频率；
[0105] A223、根据下式计算获得每个IMF分量的希尔伯特边缘谱：
[0107] 上式中，Η( ω，t)表示IMF分量c(t)的希尔伯特谱，h( ω )表示IMF分量c(t)的希尔 伯特边缘谱，T表示积分周期。
[0108] 以下结合一具体实施例对本发明做详细说明。
[0109] 实施例一
[0110] -种多干扰因素耦合的光栅尺误差补偿方法，包括以下步骤：
[0111] 误差数据库建立步骤，包括步骤1和步骤2:
[0112] 步骤1、依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每组作用 强度值下的多组误差数据；
[0113] 步骤2、基于经验模式分解算法计算在每组作用强度值下的最优误差补偿量，进而 对应建立误差补偿数据库。
[0114] 测量补偿步骤，包括步骤3和步骤4:
[0115] 步骤3、在光栅尺测量过程中采用多个传感器测量获得多种干扰因素的作用强度 值，进而与误差补偿数据库进行匹配，获得该组作用强度值对应的最优误差补偿量;本步骤 相当于通过多种干扰因素的作用强度值的组合，进行查表，获得对应的最优误差补偿量。
[0116] 步骤4、采用该最优误差补偿量对光栅尺系统进行补偿。补偿时，将光栅尺系统的 测量值减去最优误差补偿量即为补偿后的测量值。
[0117]步骤1具体为：
[0118] 通过多个误差调节器依次调节多种干扰因素的作用强度，并采用传感器测量获得 调节过程中的每组作用强度值，同时在每组作用强度值下，采用光栅尺系统和激光干涉仪 对待测对象进行多次测量后计算获得多组误差数据。误差调节器用于调节各干扰因素的作 用强度，即调节各干扰因素所带来的误差大小；传感器用于测量各干扰因素所带来的作用 强度值，即误差值。具体的，测量过程中，激光干涉仪的测量数据作为标准测量数据，将光栅 尺系统的测量数据与激光干涉仪的测量数据相减获得的差值即为误差数据。
[0119] 参照图1，步骤2具体包括步骤Α1~Α4:
[0120] Α1、将每组作用强度下的多组误差数据平均分成两个集合，分别为分析数据集和 测试数据集；因此，本实施例中，多组误差数据的数量为2η组，分析数据集和测试数据集中 每个集合的组数为η组。
[0121] Α2、基于经验模式分解算法分别对分析数据集的每组误差数据进行分解计算后， 获得每组误差数据所对应的趋势误差分量；
[0122] A3、分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进 行补偿效果统计，根据下式计算获得每个趋势误差分量对应的多个补偿余量的平均值ε 1:
[0124] ^表示第i组误差数据的趋势误差分量对应的多个补偿余量的平均值，e2j表示测 试数据集中的第j组误差数据，su表示分析数据集中的第i组误差数据对应的趋势误差分 量，i、j均为自然数，η表示测试数据集的误差数据的总组数。
[0125] e2j-Sll表示每次误差补偿操作得到的补偿余量，显然的，补偿余量越小，意味着补 偿效果越好。当对η组测试数据进行补偿后，获得η个补偿余量，对η个补偿余量求平均值获 得平均值^。当分析数据集的η组误差数据的趋势误差分量都完成补偿操作后，即获得η组 平均值^，其中最小值对应的趋势误差分量即为最优的趋势，可以作为最优的补偿量，通过 该最优的补偿量进行光栅尺误差补偿，可以获得最好的补偿效果。
[0126] Α4、将平均值最小的多个补偿余量对应的趋势误差分量作为该组作用强度下的最 优误差补偿量。
[0127] 其中，步骤Α2包括Α21~Α24:
[0128] Α21、针对分析数据集的每组误差数据，采用经验模式分解算法将其分解成多个 IMF分量和一个残余分量；
[0129] A22、分别对每个IMF分量进行希尔伯特变换，并计算获得每个IMF分量的希尔伯特 谱和希尔伯特边缘谱；
[0130] A23、根据预设筛选条件对所有IMF分量进行分析后，筛选出符合条件的IMF分量；
[0131] A24、将残余分量和所有符合条件的IMF分量累加后作为误差信号的趋势误差分 量。
[0132] 具体的，步骤A23，包括A231~A233:
[0133] A231、将每个頂F分量的希尔伯特边缘谱按照频率顺序均分为低频区域、中频区域 和尚频区域；
[0134] A232、分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而 根据下式分别计算每个区域的归一化权重因子：
[0136] 上式中，Mj, i、Mj, 2和Mj,3依次表示第j个MF分量的低频区域、中频区域和高频区域 的归一化权重因子，和叫,3分别表示第j个頂F分量的低频区区域、中频区域和高频 区域的权重因子，j为自然数。本实施例中，预设阈值优选为0.5。
[0137] A233、筛选出低频区域的归一化权重因子大于预设阈值的MF分量作为符合条件 的IMF分量。
[0138] 步骤A231~A233中，通过计算各MF分量的希尔伯特边缘谱的低频区域、中频区域 和高频区域的归一化权重因子后，根据低频区域的归一化权重因子是否大于预设阈值来筛 选该IMF分量是否符合条件。本方法通过计算归一化权重因子，根据归一化情况来筛选判 断，而不是人为的设定某参数的阈值来进行判断，从而可以使得筛选过程更为科学、严谨， 最后提取获得的趋势误差分量更为接近实际情况，即使得补偿效果更好。
[0139] 步骤A21具体为：
[0140] 针对分析数据集的每组误差数据，根据下式，采用经验模式分解算法将其分解成 多个IMF分量和一个残余分量和一个残余分量：
[0142]上式中，X(t)表示误差数据，i为自然数，m表示分解的頂F分量的总个数，Cl(t)表 示经验模式分解后的第i个IMF分量，rm( t)表示残余分量。
[0143] 步骤A22的具体计算过程为:分别对每个IMF分量进行希尔伯特变换后，计算获得 每个IMF分量的瞬时幅值、瞬时相位和瞬时频率，进而计算获得每个IMF分量的希尔伯特边 缘谱。
[0144] 具体的，步骤A22包括A221~A223:
[0145] A221、根据下式，分别对每个IMF分量进行希尔伯特变换：
[0147] 上式中，P表不奇异积分的柯西主值，c(t)表不IMF分量，CH(t)表不IMF分量c(t)的 希尔伯特变换；
[0148] A222、根据下式，计算获得每个IMF分量的瞬时幅值、瞬时相位和瞬时频率：
[0150] 上式中，a(t)表示MF分量c(t)的瞬时幅值，Φ (t)表示MF分量c(t)的瞬时相位， ω (t)表示IMF分量c(t)的瞬时频率；
[0151] A223、根据下式计算获得每个IMF分量的希尔伯特边缘谱：
[0153] 上式中，Η( ω，t)表示IMF分量c(t)的希尔伯特谱，h( ω )表示IMF分量c(t)的希尔 伯特边缘谱，T表示积分周期。
[0154]以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施 例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替 换，这些等同的变型或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，包括W下步骤:在光栅尺 测量过程中采用多个传感器测量获得多种干扰因素的作用强度值，进而与误差补偿数据库 进行匹配，获得该组作用强度值对应的最优误差补偿量； 采用该最优误差补偿量对光栅尺系统进行补偿。2. 根据权利要求1所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，还 包括W下步骤： 依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每组作用强度值下的 多组误差数据； 基于经验模式分解算法计算在每组作用强度值下的最优误差补偿量，进而对应建立误 差补偿数据库。3. 根据权利要求2所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述依次调节多种干扰因素的作用强度，测量获得在多种干扰因素的每组作用强度值下的多 组误差数据的步骤，其具体为： 通过多个误差调节器依次调节多种干扰因素的作用强度，并采用传感器测量获得调节 过程中的每组作用强度值，同时在每组作用强度值下，采用光栅尺系统和激光干设仪对待 测对象进行多次测量后计算获得多组误差数据。4. 根据权利要求2所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述基于经验模式分解算法计算在每组作用强度值下的最优误差补偿量，进而对应建立误差 补偿数据库的步骤，包括步骤： A1、将每组作用强度下的多组误差数据平均分成两个集合，分别为分析数据集和测试 数据集； A2、基于经验模式分解算法分别对分析数据集的每组误差数据进行分解计算后，获得 每组误差数据所对应的趋势误差分量； A3、分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进行补 偿效果统计，计算获得每个趋势误差分量对应的多个补偿余量的平均值； A4、将平均值最小的多个补偿余量对应的趋势误差分量作为该组作用强度下的最优误 差补偿量。5. 根据权利要求4所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述步骤A2,包括： A21、针对分析数据集的每组误差数据，采用经验模式分解算法将其分解成多个IMF分 量和一个残余分量； A22、分别对每个IMF分量进行希尔伯特变换，并计算获得每个IMF分量的希尔伯特谱和 希尔伯特边缘谱； A23、根据预设筛选条件对所有IMF分量进行分析后，筛选出符合条件的IMF分量； A24、将残余分量和所有符合条件的IMF分量累加后作为误差信号的趋势误差分量。6. 根据权利要求5所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述步骤A23,包括： A231、将每个IMF分量的希尔伯特边缘谱按照频率顺序均分为低频区域、中频区域和高 频区域； A232、分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而分别 计算每个区域的归一化权重因子； A233、筛选出低频区域的归一化权重因子大于预设阔值的IMF分量作为符合条件的IMF 分量。7. 根据权利要求6所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述步骤A232,其具体为： 分别将每个区域中的全部极大值进行累加后作为该区域的权重因子，进而根据下式分 别计算每个区域的归一化权重因子：上式中，My、Mw和Mw依次表示第j个IMF分量的低频区域、中频区域和高频区域的归 一化权重因子，mw、mw和叫,3分别表示第j个IMF分量的低频区区域、中频区域和高频区域 的权重因子，j为自然数。8. 根据权利要求6所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述步骤A233中所述预设阔值为0.5。9. 根据权利要求4所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于，所 述步骤A3,其具体为： 分别根据每个趋势误差分量对测试数据集中的每组误差数据进行补偿，并进行补偿效 果统计，根据下式计算获得每个趋势误差分量对应的多个补偿余量的平均值ει:ει表示第i组误差数据的趋势误差分量对应的多个补偿余量的平均值，e2j表示测试数 据集中的第j组误差数据，S11表示分析数据集中的第i组误差数据对应的趋势误差分量，i、j 均为自然数，η表示测试数据集的误差数据的总组数。10. 根据权利要求5所述的一种多干扰因素禪合的光栅尺误差补偿方法，其特征在于， 所述步骤Α21，其具体为： 针对分析数据集的每组误差数据，根据下式，采用经验模式分解算法将其分解成多个 IMF分量和一个残余分量和一个残余分量：上式中，X(t)表示误差数据，i为自然数，m表示分解的IMF分量的总个数，ci(t)表示经验 模式分解后的第i个IMF分量，；Tm( t)表示残余分量。
【文档编号】G06F17/14GK106091925SQ201610436459
【公开日】2016年11月9日
【申请日】2016年6月16日
【发明人】蔡念, 谢伟, 彭红霞, 王晶, 王晗, 陈新度, 陈新
【申请人】广东工业大学