一种基于移动最小二乘法的弹道曲线拟合方法
【技术领域】
[0001] 本发明涉及一种基于移动最小二乘法的弹道曲线拟合方法，属于弹道实验数据处 理技术领域。
【背景技术】
[0002] 弹道处理是一种从实验数据获得发射时状态量与弹道运动量依赖关系的过程，在 工程实践中弹道处理的发射时状态量一般为载机的发射时刻，发射高度、发射角度等，而处 理的弹道运动量一般为射程、时间等，且依赖关系需要满足光滑性且曲率小范围变化等的 假设。目前，弹道的处理主要方法包括传统最小二乘法。传统最小二乘法在弹道处理过程 中，多项式的阶次需要试算，当多项式阶次过高时，容易产生高变化率一阶导数，阶次过低 时，原始实验数据逼近精度不足，而当数据量较大、原始数据非线性性质较强时定义域往往 需要进行分段拟合，从而限制了其工程应用的通用性及便捷性。
【发明内容】
[0003] 本发明的目的在于克服现有传统最小二乘法弹道拟合方法多项式阶次需要试算、 定义域需分段及原始数据逼近精度不高等的缺点，针对移动最小二乘法中的权重函数影响 域提出一种基于De Giorgi迭代的估计算法，估计权重函数在局部邻域中的影响半径，构造 权重函数的具体形式，获得一种依赖实验数据拟合弹道函数的新方法。
[0004] 本发明为解决上述技术问题提供了一种基于移动最小二乘法的弹道曲线拟合方 法，该方法的步骤如下：
[0005] 1)依据传统最小二乘法对弹道实验数据进行拟合得到的弹道函数；
[0006] 2)根据传统最小二乘法得到的弹道函数取最小误差、平均误差和实验点误差计算 权重函数影响半径；
[0007] 3)根据权重函数影响半径构建权重函数；
[0008] 4)根据所构造的权重函数采用移动最小二乘法对弹道数据进行拟合，以得到相应 的弹道曲线。
[0009] 所述权重函数影响半径为：
[0011] 其中Cl1表示实验点所对应的权重函数影响半径，S表示平均误差，A1表示实验点 误差， 八min表不最小误差，i = 1，2，为实验点数。
[0012] 所述的步骤3)中构建的权重函数：
[0014] 其中ri= |x-Xi I /di表示紧支集影响半径，Θ表示归一因子，CM=ki，ki表示实验点初 始权重。
[0015]所述权重函数满足归一性，即：
[0017] Θ值可根据上式进行计算。
[0018] 所述最小误差、平均误差和实验点误差满足下列关系：
[0020] 其中η表示实验点数，Al表示最小误差，Ari表示平均误差，A r2表示实验点误差， ^表示大于设定正常值的权重函数，△表示期望移动最小二乘法拟合函数。
[0021] 本发明的积极效果是:本发明基于De Giorgi迭代估计的弹道处理方法，通过选取 初次拟合函数获得权重函数误差夹逼估计，计算权重函数影响半径，构造权重函数，依据移 动最小二乘法计算拟合函数四个步骤完弹道数据的处理。相比传统最小二乘法将原有分段 拟合的方式改变为紧支集邻域半径计算，本发明不需要进行分段拟合。由权重函数紧支集 邻域半径的估计的有界性及权函数的归一性，可高精度逼近原始数据，能拟合数据剧烈变 化的情况。由目标泛函的连续可微性，得到的弹道函数有足够光滑性。由权重函数的权重分 配性，不需多项式阶次的试算。
【附图说明】
[0022] 图1是本发明实施例中四阶基底、Al= Amin, Ar= Ai移动最小二乘法拟合数据 图。
【具体实施方式】
[0023]下面结合附图对本发明的【具体实施方式】做进一步的说明。
[0024] 本发明根据最小二乘法在弹道处理过程中存在的上述问题，本发明提供了一种基 于移动最小二乘法的弹道曲线拟合方法，该方法首先依据传统最小二乘法对弹道实验数据 进行拟合得到的弹道函数，计算该拟合弹道函数与实验数据的最小误差、平均误差和实验 点误差;并根据得到最小误差、平均误差和实验点误差计算权重函数影响半径;然后根据权 重函数影响半径构建权重函数;最后根据所构造的权重函数采用移动最小二乘法对弹道数 据进行拟合，以得到相应的弹道曲线。
[0025] 其中移动最小二乘法(Moving least square method，MLSM)作为一种拟合算法是 由Shepard于1968年提出，并最初应用在固体力学领域，1992年Nayroles等人最早将移动最 小二乘法用于无网格法中求解偏微分方程，并逐渐在该领域内得到了广泛应用。MLSM的基 本原理是，首先对全特性曲线在整个求解区域内分区，然后在不同的区域上用最小二乘法 拟合，即采用分区局部拟合，这对拟合精度将会有较大改善，与理论值和最小二乘法拟合结 果对比，结果表明移动最小二乘法通用性强，原始数据逼近精度高，可适用于弹道拟合。
[0026] 下面以某模型炸弹的实际飞行采样的实际飞行高度、时间作为原始实验数据，对 本发明的具体实施过程进行详细描述。
[0027] I.依据传统最小二乘法对弹道实验数据进行拟合得到的弹道函数
[0028] 选取初次拟合函数获得权重函数误差夹逼估计指:选取一定阶次的基底，根据传 统最小二乘法拟合得出的弹道函数，通过计算拟合弹道函数与实验数据的最小误差，平均 误差，实验点误差，实验点误差由原始实验数据与经典二乘法依据L 2模求得，其具体公式为
[0030] 其中W表示实验数据，u表示由经典二乘法获得的数据，
[0031] 本实施例中采用传统最小二乘法进行拟合时选取的初次拟合函数的基底形式如 下：
[0032] R = ao+aiH+a2H2
[0033] 依据最小二乘法获得拟合时间，计算拟合时间与原始数据时间的最小误差、平均 误差和实验点误差。最小误差、平均误差和实验点误差满足下列关系：
[0035] 其中η表示实验点数，△遠示最小误差，Ar1表示平均误差，Ar2表示实验点误差， iPi表示大于设定正常值的权重函数，A表示期望移动最小二乘法拟合函数。
[0036] 2.根据最小误差、平均误差、实验点误差和实验点初始权重计算权重函数影响半 径di〇
[0037]权重函数影响半径的计算过程如下：
[0038] 第一步，引入De Giorgi迭代引理，设识(〇是定义在[k0,上的非负的单调不增 函数，且满足
[0040]其中 α>〇，β>1。则 [0041 ] (p{k{)+d) = Q
[0042] 其中
[0044]第二步，给出定义
[0045]将一个集合的测度记成|e|，同时定义
[0046] E(k) = {xED| co(x)>k}
[0047] Di={x^D|〇< I I χ-Xi I I <Ri