CN105069311B - 一种远程火箭发射初态误差传播估计方法 - Google Patents
一种远程火箭发射初态误差传播估计方法 Download PDFInfo
- Publication number
- CN105069311B CN105069311B CN201510523229.8A CN201510523229A CN105069311B CN 105069311 B CN105069311 B CN 105069311B CN 201510523229 A CN201510523229 A CN 201510523229A CN 105069311 B CN105069311 B CN 105069311B
- Authority
- CN
- China
- Prior art keywords
- deviation
- vector
- acceleration
- coordinate system
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000001133 acceleration Effects 0.000 claims abstract description 109
- 230000005540 biological transmission Effects 0.000 claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims description 50
- 230000014509 gene expression Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 239000002184 metal Substances 0.000 claims description 3
- 238000003491 array Methods 0.000 claims description 2
- 239000007787 solid Substances 0.000 claims description 2
- 239000000758 substrate Substances 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 6
- 235000015842 Hesperis Nutrition 0.000 abstract 1
- 235000012633 Iberis amara Nutrition 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000012938 design process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
Landscapes
- Navigation (AREA)
- Complex Calculations (AREA)
Abstract
一种远程火箭发射初态误差传播估计方法,本发明涉及远程火箭发射初态误差传播估计方法。本发明的目的是为了解决现有发射初态误差引起的关机点位置偏差、速度偏差以及落点纵向偏差和横向偏差的计算效率低、无法充分分析出在弹道设计过程中发射初态误差的传播机理的问题。通过以下技术方案实现的:步骤一、建立动力学摄动方程;步骤二:求解远程火箭推力加速度偏差、气动加速度偏差、正常引力加速度偏差、科氏加速度偏差和离心加速度偏差;步骤三、根据步骤一和步骤二,得到远程火箭发射初态误差引起关机点位置偏差、速度偏差的近似解析解以及落点纵向偏差、横向偏差的近似解析解。本发明应用于远程火箭或运载火箭飞行动力学领域。
Description
技术领域
本发明涉及远程火箭发射初态误差传播估计方法。
背景技术
远程火箭在弹道设计和计算中一般建立在地面发射坐标系中,然而由于发射初态误差的存在,使得实际发射坐标系与标称发射坐标系存在偏差。发射初态误差包括初始定位误差(大地经纬偏差、高程偏差)和初始定向误差(垂线偏差、发射方位角偏差)。发射初态误差会带来发射坐标系中各加速度的偏差,即引力加速度偏差、推力加速度偏差、气动加速度偏差、科氏加速度偏差、离心加速度偏差,进而引起关机点状态偏差和落点偏差。根据资料表明:对于10000km射程的远程火箭,机动发射时初态误差造成的落点偏差可达到千米量级。因此,开展发射初态误差的研究对于分析远程火箭命中精度具有十分重要的意义。
虽然发射初态误差可以通过数值计算,但发射初态误差引起的关机点位置偏差、速度偏差以及落点纵向偏差和横向偏差的计算效率低,因此需要对发射初态误差的传播机理进行理论分析,而目前针对发射初态误差的研究大多集中在导航误差方面,无法充分分析出在弹道设计过程中发射初态误差的传播机理。
发明内容
本发明的目的是为了解决现有发射初态误差引起的关机点位置偏差、速度偏差以及落点纵向偏差和横向偏差的计算效率低、无法充分分析出在弹道设计过程中发射初态误差的传播机理的问题,而提出了远程火箭发射初态误差传播估计方法。
上述的发明目的是通过以下技术方案实现的:
步骤一、建立动力学摄动方程;
步骤二:根据步骤一得出的动力学摄动方程,求解远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差;
步骤三、根据步骤一得出的动力学摄动方程和步骤二中得出的远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差,得到远程火箭发射初态误差引起关机点位置偏差、速度偏差的近似解析解以及落点纵向偏差、横向偏差的近似解析解。
发明效果
本发明在弹道设计过程中在标称发射坐标系中采用小扰动假设建立动力学摄动方程,分析各发射初态误差的传播机理,获得发射初态误差的传播规律。
为便于计算效率的比较,仿真采用的计算机配置为:CPU为i3-2100,主频为3.1GHz,内存3GB。仿真结果表明:数值解耗时12.69s,利用动力学摄动方程求解出的近似解析解耗时0.12s。因此,近似解析解可以显著提高发射初态误差引起的关机点速度偏差和位置偏差以及落点纵向偏差和横向偏差的计算效率。
通过传播矩阵不仅可以得出发射初态误差引起的关机点位置偏差、速度偏差以及落点纵向偏差和横向偏差的影响量级,还可以更充分地分析出发射初态误差的传播机理。
从图3a、图3b、图4a、图4b、图5a和图5b中看出,利用本发明提出的关机点速度偏差和位置偏差近似解析解公式(56)和发射初态误差引起落点纵向偏差ΔL和横向偏差ΔZ的近似解析解公式(58)得出的结果能够与数值解(弹道求差法)吻合的很好,其产生的关机点速度偏差和位置偏差如表1所示,利用近似解析解得到的关机点速度偏差与数值解相差0.0036m/s,相对误差为3.1%,现有技术相对误差为10%,精度提高了6.9%;得到的关机点位置偏差与数值解相差0.786m,相对误差为0.4%,现有技术相对误差为1%,精度提高了0.6%。
这表明,采用本发明所提出的发射初态误差传播模型相比于数值解不仅可以提高计算效率,而且能够分析出发射初态误差的传播规律。
表2得出了远程火箭落点偏差、射程、射程角以及球面方位角。可以看出,通过数值解计算出发射初态误差产生的落点纵向偏差和横向偏差分别为1505.5m和-102.8m,利用传播矩阵CT求解出的纵向偏差和横向偏差分别为1624.3m和-85.6m,距离相差120.1m,相对误差为7.9%,现有技术相对误差为15%,精度提高了7.1%。
这表明对于射程接近10000km的远程火箭,本文所提出的发射初态误差传播模型得出的落点纵向偏差和横向偏差与数值解相比较,相对误差较小,精度较高。
附图说明
图1为本发明流程图;
图2为具体实施方式二中地面发射坐标系偏差示意图;
图3a为实施例中发射初态误差产生的主动段x方向位置偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为x轴方向位置偏差,单位为米;
图3b为实施例中发射初态误差产生的主动段x方向速度偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为x方向速度偏差,单位为米每秒;
图4a为实施例中发射初态误差产生的主动段y方向位置偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为y方向位置偏差,单位为米;
图4b为实施例中发射初态误差产生的主动段y方向速度偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为y方向速度偏差,单位为米每秒;
图5a为实施例中发射初态误差产生的主动段z方向位置偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为z方向位置偏差,单位为米;
图5b为实施例中发射初态误差产生的主动段z方向速度偏差图,横坐标为主动段飞行时间,单位为秒,纵坐标为z方向速度偏差,单位为米每秒。
具体实施方式
具体实施方式一:结合图1说明本实施方式,一种远程火箭发射初态误差传播估计方法,其特征在于,一种远程火箭发射初态误差传播估计方法具体是按以下步骤进行的:
步骤一、建立动力学摄动方程;
步骤二:根据步骤一得出的动力学摄动方程,求解远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差;
步骤三、根据步骤一得出的动力学摄动方程和步骤二中得出的远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差,得到远程火箭发射初态误差引起关机点位置偏差、速度偏差的近似解析解以及落点纵向偏差、横向偏差的近似解析解。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中建立动力学摄动方程;具体过程为:
设标称发射坐标系为o1-x1y1z1,o1为标称发射坐标系原点,x1为标称发射坐标系x轴,y1为标称发射坐标系y轴,z1为标称发射坐标系z轴,实际发射坐标系为o2-x2y2z2,o2为实际发射坐标系原点,x2为实际发射坐标系x轴,y2为实际发射坐标系y轴,z2为实际发射坐标系z轴,如图2所示;
标称发射坐标系和实际发射坐标系的差别反映了发射初态误差,发射初态误差包括初始定位误差和初始定向误差,其中,初始定位误差为大地经度偏差Δλ0、大地纬度偏差ΔB0和高程偏差ΔH0,初始定向误差为垂线偏差子午方向分量ξ、垂线偏差卯酉方向分量η和发射方位角偏差ΔA0;为了便于误差传播分析,将实际发射坐标系动力学方程与标称发射坐标系动力学方程作差后,得到动力学摄动方程在标称发射坐标系o1-x1y1z1中的表达式为
式中,为正常引力加速度偏差,为推力加速度偏差,为气动力加速度偏差,为科氏加速度偏差,为离心加速度偏差,为速度偏差对时间的导数,为位置偏差对时间的导数,为速度偏差向量;
和为发射初态误差Δλ0、ΔB0、ΔH0、ξ、η、ΔA0的函数,因此有
式中,为标称发射坐标系o1-x1y1z1中的位置向量;为标称发射坐标系o1-x1y1z1中的速度向量;为正常引力加速度向量;为位置偏差向量;λ0为发射点大地经度;Δλ0为发射点大地经度偏差;B0为发射点大地纬度;ΔB0为发射点大地纬度偏差;H0为发射点高程;ΔH0为发射点高程偏差;A0为发射方位角;ΔA0为发射方位角偏差;ξ为垂线偏差子午方向分量;η为垂线偏差卯酉方向分量;为科氏加速度向量;为标称发射坐标系o1-x1y1z1中速度向量;为速度偏差向量;为离心加速度向量。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中根据步骤一得出的动力学摄动方程,求解远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差;具体过程为:
(1)推力加速度偏差和气动加速度偏差
式中,ΔΓ21为实际发射坐标系o2-x2y2z2到标称发射坐标系o1-x1y1z1的转换矩阵偏差;为推力加速度向量;为气动加速度向量;为发射初态误差向量;GD为引起推力加速度偏差、气动加速度偏差的传递矩阵,分别表示为的函数;
的表达式为
式中,为的转置矩阵,aPx、aPy、aPz分别为的三个方向分量;
式中,为的转置矩阵,aRx、aRy、aRz分别为的三个方向分量。
(2)正常引力加速度偏差
正常引力加速度在标称地面发射坐标系中的矢量形式为
式中,为标称发射坐标系中的位置向量;为地心指向标称发射点的向量;为地球自转角速度向量;ωe为地球自转角速度大小;r为导弹所处的地心距大小,gr为引力加速度沿导弹地心矢径方向的分量;gω为引力加速度沿地球自转角速度方向的分量;正常引力加速度是一个3×1的向量;
其中,
式中:μ为地球引力常数;ae为地球长半轴,为ae=6378137m;J为地球扁率修正项,J2为带谐系数,为1.08263×10-3;φ为当前点的地心纬度;为导弹地心距矢量,为sinφ为的函数,用f表示:
和的表达式为
式中,为发射点地心矢径在地心固联坐标系中的向量;N0为酉半径;e2为第二偏心率,e2=6.7395018×10-3;GE为地心固联坐标系到发射坐标系的转换矩阵;
式(7)中酉半径N0的表达式为
由于e2sin2B0是一个小量,这里
因此,结合式(5)至式(8),可得正常引力加速度向量是λ0、B0、H0、A0、的函数,下面求解对各项的偏导数表达式;
1)
对中求偏导数为:
式中,为对的偏导数,I3为一个三维的单位阵,为r对的偏导数;
其中,
式中,μ为地球引力常数,μ=3.9860047×1014;
由于J项而导致的误差是微小的,因此可忽略J项的影响,此时有
定义弹道角速度nb为
则将公式(11)进一步简化整理为
2)
对中λ0求偏导数为
式中
其中
将公式(18)至(20)代入(17),结合(14),得到
3)
对中求偏导数为
式中,
其中
将(23)-(26)代入(22),并结合(14),得到
4)
对中H0求偏导数为
式中
将(29)代入(28),并结合(14),得到
5)
对中A0求偏导数为
式中
将(32)、(33)代入(31),并结合(14),得到
6)
对中ξ求偏导数为
7)
对中ηA0求偏导数为
因此,将1)~7)中各偏导数矩阵的表达式代入(2)中的即可得到正常引力加速度偏差的表达式。
(3)科氏加速度偏差
科氏加速度向量为
式中,为标称发射坐标系中的速度向量,ωe ×为的反对称矩阵,即
式中,ωex、ωey、ωez分别为的三个坐标分量;
因此
其中,
另外,对ξ的偏导数为
因此,将(39)至(41)各偏导数矩阵代入(2)中即可得到科氏加速度偏差的表达式;
(4)离心加速度向量为
式中,为相对速度向量;
由于因此,对偏导数为
对λ0偏导数为
对B0偏导数为
对H0偏导数为
对A0偏导数为
对ξ的偏导数为
对η的偏导数为
因此,将(43)至(49)中各偏导数矩阵代入(2)中即可得到离心加速度偏差的表达。
具体实施方式四:本实施方式与具体实施方式一、二或三不同的是:所述步骤三中根据步骤一得出的动力学摄动方程和步骤二中得出的远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差,得到远程火箭发射初态误差引起关机点位置偏差、速度偏差的近似解析解以及落点纵向偏差、横向偏差的近似解析解;具体过程为:
将公式(1)、(2)进行整理,写出状态方程的形式为
式中,Gv为速度偏导数矩阵,Gρ为位置偏导数矩阵,03×3为一个3×3零矩阵,GΔ为发射初态误差偏导数矩阵,03×6为一个3×6的零矩阵,为发射初态误差向量。
其中
从图1标称发射坐标系和实际发射坐标系得出,系统仅存在初始位置偏差,初始速度误差为零。因此,
系统状态的初始值为
式中,为初始速度向量,为初始位置向量,Lρ0为引起初始位置误差的转换矩阵,Lρ0表达式为
式中,Re为地球平均半径,Re=6371004m;
根据微分方程理论,公式(50)的解为
式中,为关机点速度偏差向量,为关机点位置偏差向量,Φ(tk)为发射时刻到关机点时刻tk的状态转移矩阵,Φ(tk,τ)为某时刻τ到关机点时刻tk的状态转移矩阵,tk为关机点时刻,τ为飞行中某时刻;
记
式中,CA为发射初态误差引起的关机点速度偏差和位置偏差传播矩阵,CΔ为GΔ引起的关机点速度偏差和位置偏差传播矩阵,CaP为推力加速度偏差引起的关机点速度偏差和位置偏差传播矩阵,CaR为气动加速度偏差引起的关机单速度偏差和位置偏差传播矩阵;
因此,远程火箭发射初态误差引起关机点速度偏差和位置偏差近似解析解为
公式(58)即在标称发射坐标系中建立了发射初态误差引起的关机点速度偏差和位置偏差
设主动段关机点在标称发射坐标系的状态量为则状态偏差量为因此,落点纵向偏差ΔL和横向偏差ΔZ与的关系为
式中,ΔL为落点纵向偏差,ΔZ为落点横向偏差,为关机点状态偏差量,T为转置,L为纵向射程,为关机点状态向量,Z为横向射程,CX为关机点状态偏差引起落点纵向偏差和横向偏差的传播矩阵,为纵向偏差对关机点状态偏差的偏导数,为横向偏差对关机点状态偏差的偏导数;
因此,远程火箭发射初态误差引起落点纵向偏差ΔL和横向偏差ΔZ的近似解析解为
式中,CT为发射初态误差引起落点纵向偏差ΔL和横向偏差ΔZ的传播矩阵。
实施例1:
对于射程为10000km的远程远程火箭,发射点大地经度为117.3°,大地纬度为39.9°,高程为10m,发射方位角为30°。选择大地经度偏差Δλ0为-5″、纬度偏差ΔB0为5″,高程偏差ΔH0为10m,垂线偏差子午方向分量ξ和卯酉方向分量均为-20″,发射方位角偏差ΔA0为10″。
具体实施方式一:对关机点状态的影响
利用公式(52)得到关机点状态偏差的传播矩阵CA、CΔ、CaP和CaR为
将发射初态误差通过传播矩阵CA、CΔ、CaP和CaR产生的关机点速度和位置偏差进行整理,如表1所示。然后将各传播矩阵得到的关机点状态偏差求和,得到了近似解析解产生的关机点状态偏差,并与数值解进行对比。
表1数值解、近似解析解以及各传播矩阵产生的关机点状态偏差
从传播矩阵和表1中看出,发射初态误差产生的关机点速度、位置偏差的数值解和近似解析解能够很好的吻合。从上面的分析得出如下结论:
(1)发射初态误差主要通过CA和CaP产生关机点状态偏差,而通过CΔ和CaR产生的关机点状态偏差较小,可以忽略。
(2)关机点速度偏差主要由CaP作用产生,通过CaP产生的关机点位置偏差是由于速度偏差导致产生的位置偏差。发射初态误差在CaP中会产生较大的影响因子,这表明由于发射初态误差的作用导致产生的推力加速度偏差是误差传播过程中一个很重要的影响因素。
(3)由于CA产生的关机点速度偏差较小,因此其产生的关机点位置偏差主要是初始定位误差部分。在CA产生的关机点偏差中,大地经度偏差Δλ0、大地纬度偏差ΔB0是造成关机点x、z方向位置偏差的主要因素,高程偏差ΔH0是造成y方向位置偏差的主要因素。而ξ、η和ΔA0主要通过CaP进行传播,在CA中的影响因子为零。
从图3a、图3b、图4a、图4b、图5a和图5b中看出,利用本文提出的发射初态误差传播模型在主动段飞行的过程中能够与数值解吻合的很好,其产生的关机点状态偏差如表1所示。利用近似解析解得到的关机点速度偏差与数值解相差0.0036m/s,相对误差为3.1%,现有技术相对误差为10%,精度提高了6.9%;得到的关机点位置偏差与数值解相差0.786m,相对误差为0.4%,现有技术相对误差为1%,精度提高了0.6%。这表明,采用本文所提出的发射初态误差传播模型的正确性。
具体实施方式二:对落点偏差的影响
由发射初态误差的传播模型得到了关机点的状态偏差,再通过公式(54)、(55)就可以得到传播权值CX、CT,即
将发射初态误差通过传播矩阵CT计算出的落点偏差进行整理,并与数值解进行对比,如表2所示。
表2导弹落点偏差参数
表2得出了远程火箭落点偏差、射程、射程角以及球面方位角。可以看出,通过数值解计算出发射初态误差产生的落点纵向偏差和横向偏差分别为1505.5m和-102.8m,利用传播矩阵CT求解出的纵向偏差和横向偏差分别为1624.3m和-85.6m,距离相差120.1m,相对误差为7.9%,现有技术相对误差为15%,精度提高了7.1%。
这表明对于射程接近10000km的远程火箭,本文所提出的发射初态误差传播模型得出的落点纵向偏差和横向偏差与数值解相对误差较小,精度较高。
Claims (1)
1.一种远程火箭发射初态误差传播估计方法,其特征在于,一种远程火箭发射初态误差传播估计方法具体是按照以下步骤进行的:
步骤一、建立动力学摄动方程;具体过程为:
设标称发射坐标系为o1-x1y1z1,o1为标称发射坐标系原点,x1为标称发射坐标系x轴,y1为标称发射坐标系y轴,z1为标称发射坐标系z轴,实际发射坐标系为o2-x2y2z2,o2为实际发射坐标系原点,x2为实际发射坐标系x轴,y2为实际发射坐标系y轴,z2为实际发射坐标系z轴;
标称发射坐标系和实际发射坐标系的差别反映了发射初态误差,发射初态误差包括初始定位误差和初始定向误差,其中,初始定位误差为大地经度偏差Δλ0、大地纬度偏差ΔB0和高程偏差ΔH0,初始定向误差为垂线偏差子午方向分量ξ、垂线偏差卯酉方向分量η和发射方位角偏差ΔA0;将实际发射坐标系动力学方程与标称发射坐标系动力学方程作差后,得到动力学摄动方程在标称发射坐标系o1-x1y1z1中的表达式为
式中,为正常引力加速度偏差,为推力加速度偏差,为气动力加速度偏差,为科氏加速度偏差,为离心加速度偏差,为速度偏差对时间的导数,为位置偏差对时间的导数,为速度偏差向量;
和为发射初态误差Δλ0、ΔB0、ΔH0、ξ、η、ΔA0的函数,因此有
式中,为标称发射坐标系o1-x1y1z1中的位置向量;为标称发射坐标系o1-x1y1z1中的速度向量;为正常引力加速度向量;为位置偏差向量;λ0为发射点大地经度;Δλ0为发射点大地经度偏差;B0为发射点大地纬度;ΔB0为发射点大地纬度偏差;H0为发射点高程;ΔH0为发射点高程偏差;A0为发射方位角;ΔA0为发射方位角偏差;ξ为垂线偏差子午方向分量;η为垂线偏差卯酉方向分量;为科氏加速度向量;为标称发射坐标系o1-x1y1z1中速度向量;为速度偏差向量;为离心加速度向量;
步骤二:根据步骤一得出的动力学摄动方程,求解远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差;具体过程为:
(1)远程火箭推力加速度偏差和远程火箭气动加速度偏差
为发射初态误差向量的函数,T为转置,可表示成
式中,ΔΓ21为实际发射坐标系o2-x2y2z2到标称发射坐标系o1-x1y1z1的转换矩阵偏差;为推力加速度向量;为气动加速度向量;为发射初态误差向量;GD为引起推力加速度偏差、气动加速度偏差的传递矩阵,分别表示为的函数;的表达式为
式中,为的转置矩阵,aPx、aPy、aPz分别为的三个方向分量;
式中,为的转置矩阵,aRx、aRy、aRz分别为的三个方向分量;
(2)远程火箭正常引力加速度偏差
正常引力加速度在标称地面发射坐标系中的矢量形式为
式中,为标称发射坐标系中的位置向量;为地心指向标称发射点的向量;为地球自转角速度向量;ωe为地球自转角速度大小;r为导弹所处的地心距大小,gr为引力加速度沿导弹地心矢径方向的分量;gω为引力加速度沿地球自转角速度方向的分量;
其中,
式中:μ为地球引力常数;ae为地球长半轴,为ae=6378137m;J为地球扁率修正项,J2为带谐系数,为1.08263×10-3;φ为当前点的地心纬度;为导弹地心距矢量,为sinφ为的函数,用f表示:
和的表达式为
式中,为发射点地心矢径在地心固联坐标系中的向量;N0为酉半径;e2为第二偏心率,e2=6.7395018×10-3;GE为地心固联坐标系到发射坐标系的转换矩阵;
式(7)中酉半径N0的表达式为
取N0的一阶近似,则N0的高阶项为
结合式(5)至式(8),可得正常引力加速度向量是λ0、B0、H0、A0、的函数,下面求解对各项的偏导数表达式;
1)
对中求偏导数为:
式中,为对的偏导数,I3为一个三维的单位阵,为r对的偏导数;
其中,
式中,μ为地球引力常数,μ=3.9860047×1014;
忽略J项的影响,此时有
定义弹道角速度nb为
则将公式(11)进一步简化整理为
2)
对中λ0求偏导数为
式中
其中
将公式(18)至(20)代入(17),结合(14),得到
3)
对中B0求偏导数为
式中,
其中
将(23)-(26)代入(22),并结合(14),得到
4)
对中H0求偏导数为
式中
将(29)代入(28),并结合(14),得到
5)
对中A0求偏导数为
式中
将(32)、(33)代入(31),并结合(14),得到
6)
对中ξ求偏导数为
7)
对中ηA0求偏导数为
因此,将1)~7)中各偏导数矩阵的表达式代入(2)中的即可得到正常引力加速度偏差的表达式;
(3)远程火箭科氏加速度偏差
科氏加速度向量为
式中,为标称发射坐标系中的速度向量,ωe ×为的反对称矩阵,即
式中,ωex、ωey、ωez分别为的三个坐标分量;
其中,
另外,对ξ的偏导数为
将(39)至(41)各偏导数矩阵代入(2)中即可得到科氏加速度偏差的表达式;
(4)远程火箭离心加速度向量为
式中,为相对速度向量;
由于因此,对偏导数为
对λ0偏导数为
对B0偏导数为
对H0偏导数为
对A0偏导数为
对ξ的偏导数为
对η的偏导数为
将(43)至(49)中各偏导数矩阵代入(2)中即可得到离心加速度偏差的表达;
步骤三、根据步骤一得出的动力学摄动方程和步骤二中得出的远程火箭推力加速度偏差、远程火箭气动加速度偏差、远程火箭正常引力加速度偏差、远程火箭科氏加速度偏差和远程火箭离心加速度偏差,得到远程火箭发射初态误差引起关机点位置偏差、速度偏差的近似解析解以及落点纵向偏差、横向偏差的近似解析解;具体过程为:
将公式(1)、(2)进行整理,写出状态方程的形式为
式中,Gv为速度偏导数矩阵,Gρ为位置偏导数矩阵,03×3为一个3×3零矩阵,GΔ为发射初态误差偏导数矩阵,03×6为一个3×6的零矩阵,为发射初态误差向量;
其中
系统状态的初始值为
式中,为初始速度向量,为初始位置向量,Lρ0为引起初始位置误差的转换矩阵,Lρ0表达式为
式中,Re为地球平均半径,Re=6371004m;
根据微分方程理论,公式(50)的解为
式中,为关机点速度偏差向量,为关机点位置偏差向量,Φ(tk)为发射时刻到关机点时刻tk的状态转移矩阵,Φ(tk,τ)为某时刻τ到关机点时刻tk的状态转移矩阵,tk为关机点时刻,τ为飞行中某时刻;
记
式中,CA为发射初态误差引起的关机点速度偏差和位置偏差传播矩阵,CΔ为GΔ引起的关机点速度偏差和位置偏差传播矩阵,CaP为推力加速度偏差引起的关机点速度偏差和位置偏差传播矩阵,CaR为气动加速度偏差引起的关机单速度偏差和位置偏差传播矩阵;
远程火箭发射初态误差引起关机点速度偏差和位置偏差近似解析解为
公式(58)即在标称发射坐标系中建立了发射初态误差引起的关机点速度偏差和位置偏差
设主动段关机点在标称发射坐标系的状态量为则状态偏差量为因此,落点纵向偏差ΔL和横向偏差ΔZ与的关系为
式中,ΔL为落点纵向偏差,ΔZ为落点横向偏差,为关机点状态偏差量,T为转置,L为纵向射程,为关机点状态向量,Z为横向射程,CX为关机点状态偏差引起落点纵向偏差和横向偏差的传播矩阵,为纵向偏差对关机点状态偏差的偏导数,为横向偏差对关机点状态偏差的偏导数;
因此,远程火箭发射初态误差引起落点纵向偏差ΔL和横向偏差ΔZ的近似解析解为
式中,CT为发射初态误差引起落点纵向偏差ΔL和横向偏差ΔZ的传播矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523229.8A CN105069311B (zh) | 2015-08-24 | 2015-08-24 | 一种远程火箭发射初态误差传播估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523229.8A CN105069311B (zh) | 2015-08-24 | 2015-08-24 | 一种远程火箭发射初态误差传播估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105069311A CN105069311A (zh) | 2015-11-18 |
CN105069311B true CN105069311B (zh) | 2018-06-12 |
Family
ID=54498676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510523229.8A Expired - Fee Related CN105069311B (zh) | 2015-08-24 | 2015-08-24 | 一种远程火箭发射初态误差传播估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105069311B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105501467B (zh) * | 2015-12-11 | 2017-03-08 | 周沁心 | 一种基于零射程线的运载火箭耗尽关机控制方法 |
CN105894126B (zh) * | 2016-04-25 | 2019-10-22 | 王景国 | 一种火箭残骸的落点预报模型生成及定位方法 |
CN107036626A (zh) * | 2016-11-08 | 2017-08-11 | 哈尔滨工业大学 | 一种远程火箭初始定位定向误差影响分析方法 |
CN106970643B (zh) * | 2017-04-27 | 2020-03-27 | 中国人民解放军国防科学技术大学 | 一种解析的卫星非线性相对运动偏差传播分析方法 |
CN107480347A (zh) * | 2017-07-24 | 2017-12-15 | 湖北航天技术研究院总体设计所 | 一种分离体散布特性预示方法 |
CN110989650B (zh) * | 2019-12-30 | 2023-10-27 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608920A (zh) * | 2008-06-18 | 2009-12-23 | 中国科学院国家天文台 | 一种组合式空间位姿精密动态测量装置及方法 |
-
2015
- 2015-08-24 CN CN201510523229.8A patent/CN105069311B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608920A (zh) * | 2008-06-18 | 2009-12-23 | 中国科学院国家天文台 | 一种组合式空间位姿精密动态测量装置及方法 |
Non-Patent Citations (3)
Title |
---|
定位定向误差对弹道导弹主动段运动的影响分析;郑伟等;《弹道学报》;20081231;第20卷(第4期);第5-7,31页 * |
弹道导弹定位定向误差传播机理分析及补偿方法研究;郑旭;《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅱ辑》;20150228;第2015年卷(第2期);第C032-49页 * |
机动发射的弹道导弹飞行诸元的快速计算;韦文书等;《哈尔滨工业大学学报》;20121130;第44卷(第11期);第7-12页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105069311A (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105069311B (zh) | 一种远程火箭发射初态误差传播估计方法 | |
US9702674B2 (en) | Method and apparatus for GPS-denied navigation of spin-stabilized projectiles | |
CN105184109B (zh) | 扰动引力作用下弹道助推段状态偏差解析方法 | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN107036626A (zh) | 一种远程火箭初始定位定向误差影响分析方法 | |
Zhu et al. | Impact time and angle control guidance independent of time-to-go prediction | |
CN111351481A (zh) | 一种基于发射惯性坐标系的传递对准方法 | |
CN107478110B (zh) | 一种基于状态观测器的旋转弹姿态角计算方法 | |
CN110044210B (zh) | 考虑任意阶地球非球型引力摄动的闭路制导在线补偿方法 | |
CN108534783A (zh) | 一种基于北斗导航技术的飞行器导航方法 | |
CN111552003A (zh) | 基于球卫星编队的小行星引力场全自主测量系统及方法 | |
CN105138808A (zh) | 基于摄动理论的滑翔弹道误差传播分析方法 | |
CN105180728A (zh) | 基于前数据的旋转制导炮弹快速空中对准方法 | |
CN105718660B (zh) | 临近空间大范围机动弹道三维包络计算方法 | |
de Celis et al. | Attitude determination algorithms through accelerometers, GNSS sensors, and gravity vector estimator | |
Si et al. | A new transfer alignment of airborne weapons based on relative navigation | |
CN114357807B (zh) | 一种用于大气层外拦截的最优中制导方法及装置 | |
CN108827345A (zh) | 一种基于杆臂挠曲变形补偿的机载武器传递对准方法 | |
RU2423658C2 (ru) | Способ управления и стабилизации подвижного носителя, интегрированная система, устройство приведения зеркала антенны в поворотное движение в двух взаимно перпендикулярных плоскостях и устройство приведения в действие дифференциальных аэродинамических рулей для его осуществления | |
CN113093776A (zh) | 一种航天器的离轨参数确定方法及装置 | |
CN111027137A (zh) | 基于遥测数据的航天器动力学模型高精度动态构建方法 | |
CN110006455A (zh) | 用于冗余惯导系统中加速度计误差参数的快速标定方法 | |
de Celis et al. | An estimator for UAV attitude determination based on accelerometers, GNSS sensors, and aerodynamic coefficients | |
Sun et al. | Piecewise attitude tracking control of a gravity gradient microsatellite for coplanar orbital transfer | |
Zhang et al. | Velocity-to-be-gained deorbit guidance law using state space perturbation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180612 |