一种基于比例回归法的二阶振动测量系统的系统参数标定方法
【技术领域】
[0001] 本发明涉及二阶振动测量系统的系统参数标定方法。
【背景技术】
[0002] 小卫星微推进系统的微小推力测量通常采用的是质量-弹簧-阻尼二阶振动系 统，而且通常为欠阻尼系统。在测量系统设计和实现之后，需要通过实验对测量系统进行标 定，以获取测量系统的系统参数，包括：静态灵敏度、固有频率和阻尼比。系统参数标定是进 行推力测量的基础工作，是至关重要的前期步骤。
[0003]目前常用的系统参数标定方法包括：
[0004] 1)对于静态灵敏度，采用阶跃响应法，将标准量作为阶跃输入信号，获取稳态输出 信号，阶跃输入和稳态输出的线性拟合系数即为静态灵敏度；
[0005] 2)对于阻尼比，采用阶跃响应法，利用系统响应的最大超调量与阻尼比的关系获 取阻尼比，或者采用自由振动法，利用对数衰减率获取阻尼比；
[0006] 3)对于固有频率，采用阶跃响应法或自由振动法获取系统振动周期，利用固有频 率与振动周期的关系获取固有频率。
[0007]系统参数标定的精度对推力测量结果有直接的影响，目前多数微小推力二阶振动 测量系统基于上述方法进行系统参数标定，但没有充分利用系统的响应信息，没有给出系 统参数的置信区间，不能够有效地评定系统参数的标定精度。如何便捷、有效、高精度地进 行二阶振动系统的系统参数标定是亟待解决的技术问题。
【发明内容】
[0008] 本发明要解决的技术问题在于：针对二阶振动欠阻尼测量系统的系统参数高精度 标定问题，提出一种基于比例回归法的二阶振动系统的系统参数标定方法，充分利用测量 系统阶跃响应中极值点和稳态点信息，利用比例回归法对已知变量的成比例约束关系获取 静态灵敏度、固有频率和阻尼比三个系统参数，同时能够获得三个系统的置信区间。
[0009] 具体技术方案为：
[0010] 第一步，在二阶振动测量系统零初始条件下，构造测量系统在某一标准阶跃输入 下的阶跃响应，通过数据采集设备获取系统响应随时间变化的数据集；
[0011] 第二步，提取该数据集中极值点信息（tei，Y(tJ)和稳态点信息（tsj，Y(<-))，并 将极值点信息（tei，Y(tJ)转化为数据集合（i，Y(tei))，其中，i= 1，2,…，j= 1，2,…；
[0012] 第三步，根据比例回归法中比例回归方程yi =bxi，设定
Xi=-i1
通过比例回归确定系数b的估计值€和置信区间sb，得到二阶系 统阻尼比的估计值和置信区间：
[0013]
[0014] 弟四步，很据比例回！y'/云甲比例回!y万桎yi=bXi，设定yi=i，Xi=tei/3T，b= ?d，通过比例回归确定系数b的估计值石和置信区间Sb，得到二阶系统振动频率《d=b的 估计值和置信区间；
[0015] 第五步，在二阶振动测量系统零初始条件下，构造测量系统在一系列标准阶跃输 入xk下的阶跃响应，通过数据采集设备获取每一个系统响应随时间变化的数据集，并提取 每一个数据集中稳态值Yk，其中，k= 1，2,…；
[0016] 第六步，根据比例回归法中比例回归方程yk=bxk，设定yk=Yk，xk=Xk，b=K， 通过比例回归确定系数b的估计值石和置信区间Sb，得到二阶系统静态灵敏度K=b的估 计值和置信区间。
[0017] 比例回归法是针对存在正比或反比关系y=bx的已知变量（x，y)提出的，其中， b为比例回归确定系数，求解b的估计值^和置信区间Sb包括以下步骤：
[0018] 第一步，通过试验获得xJPy4勺一组样本数据（xi,yj，其中，i= 1，2,…，n,n>l;
[0019] 第二步，求解b的估计1
[0020] 第三步，求解b的置信区间
t^/dn-l)，其中
心-。/2(11-1)为给定概率1-a/2的自由度为n-1的t 分布的下侧分位数。
[0021] 与现有技术项目，采用本发明可以达到以下技术效果：
[0022] 1.本发明具有较好的通用性，适用于所有二阶振动测量系统的系统参数标定；
[0023] 2.本发明实施方便，利用一次阶跃响应，即能够获取阻尼比和固有频率的估计值， 又能够获取相应的置信区间；
[0024] 3.本发明利用比例回归法，增加了变量之间的比例约束，使得比例回归确定系数 估计值具有较高的精度。
【附图说明】
[0025]图1为本发明的二阶振动系统的系统参数标定流程图；
[0026] 图2为本发明的二阶振动系统阶跃响应极值点和稳态点示意图；
[0027] 图3为本发明的比例回归确定系数的求解步骤。 具体实施方案
[0028] 结合附图对基于比例回归法的二阶振动测量系统的系统参数标定方法做进一步 详细描述。图1为本发明的二阶振动系统的系统参数标定流程图。第一步，在二阶振动测 量系统零初始条件下，构造测量系统在某一标准阶跃输入下的阶跃响应，通过数据采集设 备获取系统响应随时间变化的数据集；
[0029] 第二步，提取该数据集中极值点信息（tei，Y(tJ)和稳态点信息（tsj，Y( <-))，并 将极值点信息（tei，Y(tJ)转化为数据集合（i，Y(tei))，其中，i= 1，2,…，j= 1，2,…，极 值点和稳态点如图2所示；
[0030] 第三步，根据比例回归法中比例回归方程yi=bXi，设定
Xi=-i，
'通过比例回归确定系数b的估计值石和置信区间Sb，得到二阶系 统阻尼比的估计值和置信区间：
[0031]
[0032] 第四步，根据比例回归法中比例回归方程yfbxi，设定yfi，xi=tei/Jr，b= ?d，通过比例回归确定系数b的估计值/〗和置信区间Sb，得到二阶系统振动频率《d=b的 估计值和置信区间；
[0033] 第五步，在二阶振动测量系统零初始条件下，构造测量系统在一系列标准阶跃输 入Xk下的阶跃响应，通过数据采集设备获取每一个系统响应随时间变化的数据集，并提取 每一个数据集中稳态值Yk，其中，k= 1，2,…；
[0034] 第六步，根据比例回归法中比例回归方程yk=bxk，设定yk=Yk，xk=Xk，b=K， 通过比例回归确定系数b的估计值^和置信区间Sb，得到二阶系统静态灵敏度K=b的估 计值和置信区间。
[0035] 图3为比例回归确定系数的求解步骤。比例回归法是针对存在正比或反比关系y =bx的已知变量（x，y)提出的，其中，b为比例回归确定系数，求解b的估计值g和置信区 间Sb包括以下步骤：
[0036] 第一步，通过试验获得xJPy4勺一组样本数据（xi,yj，其中，i= 1，2,…，n,n>l;
[0037] 第二步，求解b的估计值
[0038] 第三步，求解b的置信区丨1
_t^/dn-l)，其中
心-。/2(11-1)为给定概率1-a /2的自由度为n-1的t 分布的下侧分位数。
【主权项】
1. 一种基于比例回归法的二阶振动测量系统的系统参数标定方法，其特征在于包括以 下步骤： 第一步，在二阶振动测量系统零初始条件下，构造测量系统在某一标准阶跃输入下的 阶跃响应，通过数据采集设备获取系统响应随时间变化的数据集； 第二步，提取该数据集中极值点信息（tei，Y(tJ)和稳态点信息（％」，￥(<-))，并将极 值点信息（tei，Y(tei))转化为数据集合（i，Y(tei))，其中，i= 1，2,…，j= 1，2,…； 第三步，根据比例回归法中比例回归方程yi=bxi，设定Xi =-iJT，'通过比例回归确定系数b的估计值€和置信区间Sb，得到二阶系统 阻尼比的估计值和置信区间：第四步，根据比例回归法中比例回归方程yi=bxi，设定yi=i，xi=tei/ 31，b= ?d， 通过比例回归确定系数b的估计值/〗和置信区间Sb，得到二阶系统振动频率《d=b的估计 值和置信区间； 第五步，在二阶振动测量系统零初始条件下，构造测量系统在一系列标准阶跃输入Xk 下的阶跃响应，通过数据采集设备获取每一个系统响应随时间变化的数据集，并提取每一 个数据集中稳态值Yk，其中，k= 1，2,…； 第六步，根据比例回归法中比例回归方程yk=bxk，设定yk=Yk，xk=Xk，b=K，通过 比例回归确定系数b的估计值^和置信区间Sb，得到二阶系统静态灵敏度K=b的估计值 和置信区间。2. 如权利要求1所述的一种基于比例回归法的二阶振动测量系统的系统参数标定方 法，其特征在于：比例回归法是针对存在正比或反比关系y=bx的已知变量（x，y)提出的， 其中，b为比例回归确定系数，求解b的估计值$和置信区间Sb包括以下步骤： 第一步，通过试验获得xJPly^勺一组样本数据（xi,yj，其中，i= 1，2,…，n,n>l; 第二步，求解b的估计值第三步，求解b的置信区间&1/2(11-1)其中1为给定概率1-a/2的自由度为n-1的 t分布的下侧分位数。3.如权利要求1所述的一种基于比例回归法的二阶振动测量系统的系统参数标定方 法，其特征在于：该标定方法适合阻尼比为0. 1 <G< 0. 2的二阶振动系统。
【专利摘要】本发明公开了一种基于比例回归法的二阶振动测量系统的系统参数标定方法，目的是较便捷地获取二阶测量系统有效的系统参数标定结果。技术方案是构造二阶测量系统在标准阶跃输入下的阶跃响应，根据阶跃响应的极值点信息和稳态点信息，基于比例回归法，获取静态灵敏度、固有频率和阻尼比三个系统参数的估计值和标准差。采用本发明能够充分利用二阶系统的阶跃响应信息，利用比例回归法对变量之间成比例的约束关系，直接获取系统参数的估计值和置信区间，可以获得较好的标定结果。
【IPC分类】G01L25/00
【公开号】CN104990666
【申请号】CN201510279618
【发明人】金星, 周伟静, 叶继飞, 吴洁
【申请人】中国人民解放军装备学院
【公开日】2015年10月21日
【申请日】2015年5月27日