延长捷联惯组稳定期的方法研究
第28卷第6期2007年11月
宇航学报
V01.28November
No.62007
JournalofAstronautics
延长捷联惯组稳定期的方法研究
徐军辉,汪立新,钱培贤
(第二炮兵工程学院304教研室.西安710025)
擒要:针对捷联惯性测量组合(简称捷联惯组)稳定期较短的问题,在工具误差对导弹落点精度影响分析、稳定性分析以及统计分析的基础上,提出通过验前信息建模预测与射前标定相结合的方法来延长捷联惯组的稳定期。弹道仿真结果表明这一方法能够在导弹命中精度不受影响的情况下,延长捷联惯组稳定期,减少惯组的测试次数,提高其使用性能,具有良好的工程应用前景。
关键词:捷联惯组;稳定期;最小二乘支持向量机;建模预测;射前标定
中图分类号:研65.2
文献标识码:A
文章编号:1000.1328(2007)06-1559.06
0引言
差对落点的影响占整个落点偏差的70%一∞%,而在
捷联惯组已经广泛应用于航天和国防等诸多领域。但是受多种因数影响,目前捷联惯组的稳定期较短,为了保证其使用精度,必须采用定期循环测试方案。由于捷联惯组的使用寿命(由惯性器件的累积工作时间表示)有限,反复测试会给惯组的寿命带来损耗,同时也会给使用单位带来很大的工作负担。这些都会直接影响和制约捷联惯组的使用性能。捷联惯组稳定期主要由硬件性能所决定,但是要在硬件和工艺方面的改进面I临着投入大、周期长、见效慢等问题。因而这一问题多年来一直难以找到合适的解决途径。随着捷联惯组使用数量的不断增加,积累了大量的测试数据,因而借鉴前苏联和美国的经验…,即通过建模补偿技术来研究和解决这一问题成为一种新的、可行的途径。参考文献[2]的统计分析表明有一部分捷联惯组的误差系数长期稳定性和短期稳定性都比较好;而参考文献[3]的研究表明,在捷联惯组常用的误差系数中,只有一部分系数对武器的命中精度影响比较大。因此本文提出通过验前信息建模和射前标定相结合的方法在不影响捷联惯组使用精度的情况下延长捷联惯组的稳定期。1误差系数对导弹命中精度影响分析
以上研究结果为延长捷联惯组稳定期提供了一
制导误差是影响导弹命中精度的一个重要因素,
个可行的方案:即对那些比较稳定而又对导弹命中导弹落点偏差中有绝大部分是由于惯性测量组合的测量误差引起,因此要提高导弹的射击精度,就必须分离工具误差,求出各项误差对导弹落点的影响。参考文献E3]建立了弹道仿真系统,在建立的仿真系统基础上进行干扰弹道解算,计算出各项误差系数对导弹落点的影响。弹道仿真结果表明,某型捷联惯组常用的误差系数中,影响导弹落点精度50米以上的误差系数有11个,这其中有6个误差系数对落点精度的影响达到100米以上。而其余的误差系数对导弹命中精度的影响相对要小的多。
参考文献[2]表明,这些误差系数中,大部分误差系数的长期稳定性和短期稳性都比较好,这其中包含有9个对导弹命中精度影响较大的系数,影响100米以上的有4个,影响50米以上的有5个。统计分析表明,在没有经过长途运输、没有更换仪器和测试环境没有出现异常变化的情况下,这些误差系数具有较好的正态性,其变化比较稳定。另外有6个误差系数的短期稳定性较好,但是长期稳定性相对较差,这其中包含有两个对命中精度影响较大的系数,影响都在100米以上。另外4个对命中精度影响很小。
对于采用捷联式惯性制导的导弹来说,其制导工具误
收稿日期:2006-07.25;修回日期:2007-01.25
万方数据
1560宇航学报
第28卷
精度影响较小的误差系数利用历次测试结果建立其预测模型,来预测惯组误差系数在未来一段时间内的变化情况,并用其部分的代替测试结果,而对于另外一些对导弹命中精度影响大,或者不稳定的误差系数则采用射前标定方案,在临射前利用单次通电特性对其进行标定。由于捷联惯组目前利用的是逐次通电特性(即装订的参数是上次通电标定结果),能的因素很多,比如,测试环境的差异,温度的变化,运输的影响等等。由于我们并不清楚惯组误差系数的变化趋势和影响因子之间的关系,为了更有效的进行建模预测,需要利用相空间重构的思想,将其转化成矩阵形式以获得数据间的关联关系,把其中蕴藏的信息充分的显露出来,挖掘到尽可能大的信息量。以便提高建模预测的精度。对于给定的一组捷联惯组
而单次通电特性的精度比逐次通电特性的精度要提高3至5倍H],所以完全能够满足这些误差系数标定的精度要求。从以上分析可以提出如下方案:对比较稳定的误差系数进行验前信息建模预测,而对另外6各误差系数进行射前标定,这6个误差系数对导弹命中精度的影响都很大,包括4个陀螺的误
差系数(D。,,D。:,D…D。:)和2个加速度计的误差
系数(Kh,Ko,),其中K.,,‰。的长期稳定性较差。2捷联惯组历次测试数据的预处理2.1非等间隔测试数据的等间隔化
捷联惯组的测试应该严格按照稳定期要求进行测试。但从得到的数据看,一般不可能完全满足这个要求,因而不能保证测试结果的等间隔性。所以首先要解决的是测试数据等间隔化。另外由于惯组测试次数较少,样本容量也无法满足模型训练所需样本量的要求。参考文献[2]提出了二次修正插值法来实现测试数据的等间隔化和样本容量扩大。所谓二次修正插值法,即通过两次使用样条函数插值法来产生能够满足模型训练要求的训练样本。第一次插值的目的是实现测试数据的等间隔化,即以原始测试数据为基础,以稳定期为间隔,采用样条函数
插值法,对惯组历次测试数据进行第一次插值。第
二次插值的目的是为了扩大样本容量,解决小样本建模容量不足的问题,即以第一次插值所得等间隔序列的元素作为基本点,在相邻的两个基本点之间进行等点数插值,插值点的个数可根据建模所需样本容量来选取,一般最小间隔取为月的整数倍。2.2样本维数的扩展
捷联惯组是一个非常复杂的系统,在实际测试中得到的只是其历次测试数据,反映的主要是时间的影响。而实际上由于惯组是一种精度要求高,结构复杂,对测试环境要求比较苛刻的系统,因而影响其性
万
方数据Y
5
对于嵌入维数的选取尚无严格意义上的理论依pre—
error,FPE)评价模型的预测误差,并根据误差
SeE(k)=筹寺:
(2)
口:=弛n,)2寺。萎。
^々
1
~一“
【Xt一(∑(口。一口:)_|}(”互。))+6】2
数据建模预测
捷联惯组使用寿命有限,测试次数较少,因而捷历次测试数据|sl,…,s。},经过相空间变换后,得到
用于预测器学习的样本如(1)式所示。
据。文献E5]提出采用最终误差预报准则(finM
diction
大小选取嵌入维数m。
式中:
式中:厅。——用于训练的数据;m——需要确定的嵌入维数。从(1)式可看出,当m值增大时,残差仃:
将减少,因此总可以找到一个最优值m使得FPE达到最小。
3基于最小二乘支持向量机的捷联惯组历次测试
联惯组历次测试数据建模预测属于小样本建模问题。由于支持向量机具有任意逼近非线性映射的能力,采用的结构风险最小化准则,在最小化样本点误差的同时,缩小模型泛化误差的上界,即最小化模型的结构风险,从而提高了模型的泛化能力,这一优点在小样本学习中更为突出。因此本文采用最小二乘支持向量机建立捷联惯组历次测试数据的预测模型。
Il
第6期徐军辉等:延长捷联惯组稳定期的方法研究
156l
3.1最小二乘支持向量机
最小二乘支持向量机估计算法在优化目标中的损
失函数为误差毫的二范数n。17。,故最优化问题为:
.1
曾圳∞II2+iC台"I:2;
s.t.
Y‘=OJ・妒(并‘)一6+£,i=1,…,Z
咖愚知)=丢御‘∞+虿C善1£一
瓦aL=o,硒aL=o,琵aL=o,瓦aL=o
∽童一之量/c]…{:]
㈤
万
方数据则有一步预测模型为
n一靠
Z+。=∑(白;)K(Z・Z+,)+矿,
l=1
.『=1,2,…,P
(9)
式中:Z+。=(‘,‘一。,…,s:一川)
一般的可得到z步预测模型为:
Z“=∑(口;)K(∥・矗+。)+矿
(10)
式中:Z+。=(《小。,《小:,…,一儿训),根据式(10)
可以实现对捷联惯组的P个误差系数的多步预测。4射前标定
射前标定就是在不动弹上控制系统仪器与不改变测试发射工作流程的情况下,在发射阵地,导弹临
发射之前,利用导弹起竖前后和起竖过程中的信息,
标定捷联惯组误差系数(D。,,D∽D…D。,和K。,,
Ko,)。捷联惯组误差系数标定要求精度高,因此必须遵循两个原则:一是必须有效激励要标定的参数;二是尽量消除或减小干扰因素对要标定参数的影响。根据上述原则,利用导弹起竖过程的不同状态,综合考虑标定精度与导弹发射快速性与可操作性,射前标定过程如图1所示。
I
区啐正卜{驷
臣夏回
I臣扫l鉴查!.!垄墅!
l
臣圆
I‘
图1射前标定过程
Fig.1
Theprocedureofpre-launch
状态1时,导弹基本处于水平状态;状态3时,为导弹起竖到45。位置;状态5为导弹起竖到900位置。坐标转换关系如图2所示。此时惯组坐标系(OX。y6z。)的O邑基本处于水平状态。取发射坐标
1562宇航学报第28卷
系为基准坐标系(OX,YIZ,),且ox,与北向的夹角为A,则基准坐标系到惯组坐标系的变换矩阵为:
“P:釜鬻蓦
sinAsin妒‘coscp‘)+ccJzCOS7‘sin9isinAcos妒‘sin7tcosq,i)+∞zcos7icosq)‘
∞i=ct,_lvcosAsinl9bjcos7i一∞zsin7‘+∞ⅣsinAcos妒‘cos7'‘
喜耋嘲
删州删
埘瞄叩9
cc,五=∞Ⅳ(cosAcos妒‘costpi+cosAsin‘b‘sin7‘8in妒‘一sinAcos驴‘sin7‘sin5p‘+
c|’一2(cJ_lv(cosAsin妒‘sin7‘cos妒‘一cosAcos驴‘sin妒‘一sinAsin妒‘sinl9v‘一
i=1,2,3
(12)
将(12)、(13)式代人到(14)、(15)式,解算方程即可获得相应系数的射前标定值。详细情况见参考文献[14]。5标准弹道仿真
延长稳定期的目的为了减少惯组的测试次数,提高其使用性能。但是有一个前提,就是在延长稳定期的同时不能影响导弹的命中精度。为了检验本文所提出的方法对导弹命中精度的影响,设计了一种标准弹道仿真方案。弹道仿真原理如图3所示。其原理是:将由标准弹道方程组仿真出的角速率和
图2坐标转换
Fig.2
Conversionofcoordinates
、
视加速度(舀,,舀,,a,和痧,,谚,,影;)代人到惯组误差模型中,利用标定的误差系数可以仿真出陀螺
加速度计通道的输入分别为:
仪和加速度计输出脉冲(』、r&。,JI、r研。,』\r&。,JIv舡。,』、r肿。,
(13)
丙捌)。然后将预测的误差系数+射前标定的误差系数代人到误差模型中,利用仿真的脉冲再计算出
f
A一2
gcos7tsin纯
gcos7‘cos纯i=1,2,3
{A,‘2
LAi=一gsin7j
其中A由瞄准信息给出,y,9由加速度计输出计算,而9则可在瞄准结束后,通过反相导航法获得,详细情况见参考文献[14]。
陀螺仪通道的动、静态简化误差模型为:』V盘/El。=叫,+E。‘I,r+J艮∞:+Do,+
DI,A,+D2,A,+D3,A:
Ⅳ毋/E,,=E,yto;+叫y+E可to;+Do,+
角速率和视加速度(旁,,旁,,口;和叽,咖,,吮)。在
相同的诸元条件下,以标准弹道为基础,进行制导和导航计算,计算出落点偏差△L和AH。6实例分析
某型捷联惯组共有14次历次测试结果,用前13次测试结果作为建模预测用,以第14次测试作为当
D14。+D24,+D3一,
Ⅳ盘/EI:=E盥叫。+E,:tor+∞:+Do,+
(14)前测试,用于和预测结果进行比对。对D。,,D。:,
D…D。:和墨,,K,进行射前标定。在支持向量预
测模型中,训练样本的维数可取为4。核函数采用高斯径向基函数,口2可根据参考文献[2]中的误差系数合格性分析标准来分别确定;e可根据参考文
(15)
献[2]中的误差系数准确性分析标准来分别确定;常数C取为500。标准弹道仿真结果如表1所示。
.Dl;A,+D24,+D3:A:
加速度计通道的简化静态误差模型为:
fⅣ血/KI,=Ko,+A。+K二^y+疋A:
\N姆|Ktr=Kor+K吖AI+Ar+K玎Az
L^k/KI;=Ko:+K譬A,+K,zA,+A,
万方数据
第6期
徐军辉等:延长捷联惯组稳定期的方法研究
1563
图3弹道仿真
Fig.3
Standard
trajectoryemulation
表1标准弹道仿真结果
Table1
Theresultsofthestandard
traieetorr
CEP(横向)(m)CEP(纵向)(m)
6.45892
相对于导弹命中精度的要求,这个偏差可以忽略。
7结论
本文提出了通过验前信息建模和射前标定相结合的方法在不影响导弹命中精度的情况下延长捷联惯组的稳定期。标准弹道仿真结果表明这一方法在理论上是正确的,在工程应用上是可行性。但是要用于实际工程应用,还需要作进一步的工作:(1)规范对惯组的测试,以得到完整、准确的测试数据,这是保证建模预测准确性的前提;(2)建立各种补偿模型(运输的影响、累计通电时间的影响、主要环境因数的影响),使本文所提出的方案能够适应不同的情况以及不同的环境要求;(3)由于条件有限,目前只在实验室状态下进行了射前标定方案的验证实验,因此需要将其与具体应用背景相结合,研究在发射阵地的射前标定技术,主要是动基座条件下的数据采集和处理技术。分析表明,这一方法可以将捷联惯组的稳定期延长两倍以上。
万
方数据参考文献:
李庆中.国外提高导弹、运载火箭精度的一些措施[C].航天科技情报研究报告系列文集(一),1991[LIQing-曲m唱.so∞
Bletlsures
abouti=provi,lgthelw,eurneyof
miBsile
andcalllJeltrocketabroad[R].Re.archReport(on6)of
nstxonautiesScientific
and
TechnicalInformation,1991
[2]徐军辉.联惯性测量组合数据分析方法研究:学位论文[D].第二炮兵工程学院,2003[XlJJun.hui.Dnalysi8
MethodStudy
0f
Sll~llJ
Data.Master’aegree
Paper[D].The
SecondArtillery
Ellgi—
neerinlgCollege,.2003J
[3]鲜勇,张毅,制导工具误差对导弹射击精度的影响分析[J].飞行力学,1999:17(2)[XIAN
Yong,ZHAI、IGYi.Analysis
On
effect
of
a
mi68ile’8
guicla№insta'ument
erl'ol"Oil
fire
accuracy[J].niO,t
Dynanlic8.1999:17(2)]
[4]钟万登.液浮惯性器件[M].宇航出版社,1992[ZItOI、I(;w醣
de唱.LiqⅡ试InertiⅡm旆№[M].At打tmsutien,bli出%}hm.1992][5】
尉询楷,李应红,张朴,路建明.基于支持向量机的时间序列预测模型分析与应用[J].系统工程与电子技术,2005,03[WEI
Xim-kai,U
Ying-hong,ZttAN(;Pu,LU
Jian・面嵋.AMlysis
and
applicationsof
time・耐es
fot'ee.uting
model
viasupport
ved膻nsa-
ehi,,es[J].Sy,,ms
Engineering
and口∞咖i∞,2005,03]
[6]
l~lollerKR,SmolaAJ,RIttsehG.el11.Plt,edietinlgdme
6d∞with
illlploOi't
Vl鹭tOl"machines【C],,Pr∞∞击n庐ofIntemntionalCDⅡfe-即∞
m
Artifieinl.NaIr丑lNetworks。1997."999—1004
[7】s,,yke∞J^K,Vnndewalle
J.k武“iuares目Ip州YeCtOl"mehil"
elmifiell[CJ.Neural
Ptoee商nlgLetters,1999,9(3)..293—300
[8]
sI|蛔J
AK,Vnndewalle
J.Mnltiel,柚s
leastsquaressupport
VeCtOr
mehineB【JJ.International
JointCollfel'ence
oB
Neural
NeI啊讪,
Wnshington。USA,1999,2:∞O一903
[9]
&呵k傩J
A
K,Vnndewnile
J.R嘲m哪J耐squnrms,蛔_,姗,47(7)
support删加
珊Id血I∞[J].IEEETrnlmetiom
GEt
CircuitsEl
[10]
su,jke,,,J
AK。YandewalleJ,De
Moor
B.Optimal劬0lhl嘲曩
IIcluresHIppcnw咖machine,Neural
Network,2001,14
s,,ykens
JA
K.Nonlinenr
modeling
mt叨n】嵋I州孵tor
machines
【CJ,,Phceedin粤B
ofthe18thIEEE
In嘲n皿明舶枷蚰“l~lemmmem
Technology
Cor面el'ence,2001,1(1)
[12]Vap.1kV.Thel、IntureofStati硝iealⅧ玛嘲[M].N_York:
Springer-Verlag,1999
[13]
邓乃扬,田荚杰.数据挖掘中的新方法支持向量机[M].科学出版社,2004[DENCl、lai—yanlg。TIAI、IYing-jie.TheNewMethod
inDalaNining
SVI~I[M].SeientitlePubliBlaingtiott∞,2004]
[14]
肖正林,徐军辉,钱培贤.惯性组合误差射前标定方法研究[J].弹箭与制导学报,2004:24(4):362—366[XIAOZheng-lin,
XIJ
Jim-hui,QI^NPei-xian.R∞eareh
011thecalibrationdmethod
lh
IlVlu
er呻before
tnnehi,lcEJ】.Joumalof
Projectiles,Rod【咖,
Missilesand
Guidanee,2004。24(4):362—366]
[15]张毅.杨辉耀.李俊莉.弹道导弹弹道学[M].国防大学出版社,1999[ZItAI、IG
Yi,YAI、IGI-lni—yao,ILl
Jlm-11.B-Ⅱ谴icMi∞丑e
Ballistics[M].1YatioMlDdbM
IJnivenIi竹№,1999
作者简介:徐军辉(1974一),男,博士,主要从事惯性技术及小样本建模的研究。通信地址:第二炮兵工程学院304教研室(7100"Z5)
电话:(029)S474463S
E-mail:∞em—xjh@163.om
1564
宇航学报第28卷
The
Method
StudyofIncreasetheStabilization
XUJun—hui,WANG
Period
of
SIMU
Li-xin,QIAN
Pei-xian
(Lab304,theSecond
ArtilleryEngineeringCollege,Xi’¨710025.China)
Abstract:Aiming
at
theshortstabilizationpm40dofShpdown
InertialMeasurement
to
Unit(SIMU),a
newmethod
combining
on
thepriorinformationmodelingwithpnv-launchingcalibrationispresentedbilityanalysis,analysisofimpactthismethodfeet
On
Call
liberalizethestabilizationperiodofSIMUbased
sis-
accuracy
affectedbytool-errorandstatisticanalysis.Thestandard
trajectoryemulationshowsthat
no
liberalizedthestabilization
accuracy
period
r
reducedthetesttimeandimpmvetheperformanceofSIMUwith
a
obviousef-
theimp∞t
ofmissile.Thismet.hodhas
great
valueinengineeringapplication.
vector
Keywords:SIMU;Stabilizationperiod;Leastsquaressupport
bration
machines;Modelingandforecasting;Pre—hunchcall—
螬_蜘寺喇”寺“却峥o.却e“却oo.细eo.知寺o.杈峥o.和e“却々“知e“知峥“却寺畦mo“和e“却e“却々o.却e“酗专“*睁d知々喇"々州细々“}“良
(上接第1519页)
FindingandAnalysisofNewChaoticAttractorsinEuler’S
DynamicalEquation
ZHOUFeng—qi,KONGLing-yun
(Cdlnge
of
Astronautics.NorthwestemPolytechn/calUniversity.Xi’缸710072,China)
Abstract:Nonlinear
bodyand
chacacttristieofEnler’Bdynamicalequationisanalyzedinthispaper.Periodicmotionofthefreer!igid
ate
multifariouschaoticmotionsofthissystem
attractor
are
investigated.Some
conclusionswhich
are
differentfrom
system
trajectorymanifol—
I由rem
systems
dstheoryaboutfamily.Periodic
educedbyanalyzingthe
attractors
or
especialexampleofthissystemNewton-Leipnik
a
and
attractom.quasi-periodic
andseriesofnew
chaoticattractorsa托foundfromcontinuousdynamicalsys・
andthe
tem.Someconclusionsabouttheeduced.
convergent
emanativemechanismofsystemmotion
configurationof
attractorsare
Keywords:Euler’8equation;Nonlinear;Chaoticattractor;Newum-Leipniksystem;Lorenz
systems
万方数据