CN109470252A - 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 - Google Patents
一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 Download PDFInfo
- Publication number
- CN109470252A CN109470252A CN201811235049.XA CN201811235049A CN109470252A CN 109470252 A CN109470252 A CN 109470252A CN 201811235049 A CN201811235049 A CN 201811235049A CN 109470252 A CN109470252 A CN 109470252A
- Authority
- CN
- China
- Prior art keywords
- follows
- landing
- constraint
- optimization
- section
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000005457 optimization Methods 0.000 title claims abstract description 62
- 238000001228 spectrum Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000005484 gravity Effects 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 125000004432 carbon atom Chemical group C* 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 239000012855 volatile organic compound Substances 0.000 claims description 3
- 206010048669 Terminal state Diseases 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 3
- 230000003416 augmentation Effects 0.000 description 2
- 101150036841 minJ gene Proteins 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,属于飞行器轨迹优化与制导技术领域。所述方法包括:步骤一:建立垂直起降重复使用运载器返回着陆轨迹优化模型;步骤二、对所述返回着陆轨迹优化模型进行离散化处理;步骤三:对离散的轨迹优化问题进行凸化处理;步骤四:建立保证收敛性能的序列凸化算法。本发明提出的方法具有高精度,高准确度和求解快速性的特点。
Description
技术领域
本发明涉及一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,属于飞行器轨迹优化与制导技术领域。
背景技术
精确软着陆是一个经典问题,在国内外有着较长的研究和实践历史。近年来随着美国太空探索技术公司的猎鹰9(Falcon-9)火箭多次发射与子级回收成功,垂直起降可重复使用运载火箭(Vertical Takeoff Vertical Landing Reusable Launch Vehicle,VTVL-RLV)子级返回技术引起了国内外学者的广泛关注。子级返回制导是VTVL-RLV的一项关键技术,子级返回过程中,需要利用有限的控制能力实现大范围减速、满足过程约束和苛刻的定点垂直着陆终端约束,同时要尽量减少返回过程中的燃料消耗以提高火箭运载能力并应对突发情况。此问题可描述为强约束下的轨迹优化问题。在复杂的飞行环境下,传统的离线轨迹优化无法满足返回制导的需求。
发明内容
本发明为了解决传统的离线轨迹优化无法满足返回制导的需求的技术问题,在基于在线优化的闭环制导方法框架下,解决垂直起降重复使用运载器高精度返回着陆制导问题,提出了一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法。所采取的技术方案如下:
一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,所述优化方法包括:
步骤一:建立垂直起降重复使用运载器返回着陆轨迹优化模型;
步骤二、对所述返回着陆轨迹优化模型进行离散化处理;
步骤三:对离散的轨迹优化问题进行凸化处理;
步骤四:建立保证收敛性能的序列凸化算法。
进一步地,步骤一所述建立垂直起降重复使用运载器返回着陆轨迹优化模型的具体过程如下:针对垂直起降重复使用运载器返回着陆问题,建立以燃料最省为性能指标的轨迹优化模型,所述轨迹优化模型为:
minJ=-mf (1)
其中,V为速度,γ为航迹倾角,ψ为航向角,φ为纬度,θ为经度,r为地心距,m为质量,Pi为火箭发动机推力在箭体坐标系下的分量,D为子级受到的空气阻力,fM为地球引力常数,Ω为地球自转角速度,PTotal为发动机总推力,Isp为发动机比冲,g0为海平面重力加速度,εH、εV、εγ、εφ、εθ分别为终端高度、速度、倾角和经纬度容许误差,mdry为子级干重;mdry为子级干重;J为性能指标函数;mf为终端质量;ρ为大气密度;qmax和分别表示容许最大动压和热流密度;kQ为与飞行器外形相关的常数;Hf、Vf、γf、φf、θf、φd和θd分别表示终端状态值。
进一步地,步骤二所述离散化处理过程中,火箭子级的滑行段飞行和动力段飞行在基于伪谱法的离散优化模型建立过程中分为两个阶段,在两个阶段内分别对问题进行伪谱离散并设置阶段间连接条件,所述离散化处理的具体过程包括:
第一步:根据VTVL-RLV子级动力学模型,滑行段的动力学flip-Radau伪谱离散方程为:
其中,D为flip-Radau微分矩阵,f(x,u)为微分方程的右端函数,τi,(i=1,...,N1)为(-1,1]区间内的配点,N1为配点数量,τ0为在-1处的离散点,x1为滑行段的状态变量,u1为滑行段的控制变量,由于滑行段发动机不开机,u1为0,和分别为滑行段的起止时间,
第二步:根据VTVL-RLV子级动力学模型,动力段的动力学flip-Radau伪谱离散方程为:
并且,两阶段间的连接条件为x1_f=x2_0;
第三步:在滑行段发动机控制量为零,仅考虑动力段的动压和热流约束,则滑行段的离散形式为:
动力段发动机推力约束的离散形式为:
其中,Cq、分别表示离散化的约束函数;r2、V2分别表示动力段的地心距与速度。
进一步地,步骤三所述对离散的轨迹优化问题进行凸化处理的具体过程包括:
步骤1:引入新的控制变量Γ对原控制约束进行松弛,将问题中控制变量约束做如下变换:
Pmin≤Γ(τi)≤Pmax,(i=1,...,N2) (10)
其中,Pmin和Pmax分别表示最小和最大推力值;
步骤2:将动力段动力学约束重新列写,所述动力段动力学约束如下:
其中,farg表示增广的有端动力学函数,yi表示增广的优化变量;
步骤3:对公式(11)中的farg进行一阶泰勒展开,所述farg的一阶泰勒展开式如下:
其中,下标i表示在第i个配点处的约束,xk和为第k次迭代时泰勒展开的参考点,h(x)是动力学中与地球自转相关的项;H.O.T.表示高阶小项;和分别为线性化系统的状态矩阵和控制矩阵,且有:
上式中,
步骤4:所述farg的一阶泰勒展开式得出线性化的动力段动力学约束方程,所述线性化的动力段动力学约束方程如下:
步骤5:对过程状态约束进行线性化处理,获得线性化的过程状态约束方程如下:
其中,
c11=ρ0e-κ·HV,
进一步地,步骤四所述序列凸化算法的具体过程为:
Step1:输入给定初始优化变量x1,给定迭代停止条件以及最大迭代次数kmax,其中,最大迭代次数kmax的取值范围为10至20次。
Step2:求解上述凸化子问题,得到最优解Δxk;如果求解结束,返回成功标志与解xk;如果k>kmax,结束并返回失败标志;否则,进入step3; 和分别表示约束残差和相应的收敛条件。
Step3:置xk+1:=xk+Δx k:=k+1,返回Step 1。
进一步地,步骤四所述序列凸化算法中,第一次迭代的参考点由粗选的初值提供,随后的每次迭代过程中参考点取为上一次迭代得到的最优解,并且,在所述凸化处理加入信赖域约束,所述信赖域约束为:
其中,δx和分别为针对各个状态和控制变量设计的信赖域值。
本发明有益效果:
本发明提出的一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,首先针对垂直起降重复使用运载器返回着陆问题,建立了相应的轨迹优化模型。随后,采用flip-Radau伪谱法对问题进行了离散化,并提出了一种序列凸化算法对问题进行凸化合求解。基于伪谱法的离散化以及凸优化方法收敛,本发明提出的方法具有高精度,高准确度和求解快速性的特点,适合发展为在线轨迹优化方法。且为后续开展基于在线轨迹优化方法的闭环制导方法打下了基础。对于未来的垂直起降重复使用运载器返回着陆制导方法具有借鉴和参考价值。
具体实施方式
下面结合具体实施例对本发明做进一步说明,但本发明不受实施例的限制。
实施例1:
一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,所述优化方法包括:
步骤一:建立垂直起降重复使用运载器返回着陆轨迹优化模型;
步骤二、对所述返回着陆轨迹优化模型进行离散化处理;
步骤三:对离散的轨迹优化问题进行凸化处理;
步骤四:建立保证收敛性能的序列凸化算法。
其中,步骤一所述建立垂直起降重复使用运载器返回着陆轨迹优化模型的具体过程如下:针对垂直起降重复使用运载器返回着陆问题,建立以燃料最省为性能指标的轨迹优化模型,所述轨迹优化模型为:
minJ=-mf (1)
其中,V为速度,γ为航迹倾角,ψ为航向角,φ为纬度,θ为经度,r为地心距,m为质量,Pi为火箭发动机推力在箭体坐标系下的分量,D为子级受到的空气阻力,fM为地球引力常数,Ω为地球自转角速度,PTotal为发动机总推力,Isp为发动机比冲,g0为海平面重力加速度,εH、εV、εγ、εφ、εθ分别为终端高度、速度、倾角和经纬度容许误差,mdry为子级干重。
步骤二所述离散化处理过程中,火箭子级的滑行段飞行和动力段飞行在基于伪谱法的离散优化模型建立过程中分为两个阶段,在两个阶段内分别对问题进行伪谱离散并设置阶段间连接条件,所述离散化处理的具体过程包括:
第一步:根据VTVL-RLV子级动力学模型,滑行段的动力学flip-Radau伪谱离散方程为:
其中,D为flip-Radau微分矩阵,f(x,u)为微分方程的右端函数,τi,(i=1,...,N1)为(-1,1]区间内的配点,N1为配点数量,τ0为在-1处的离散点,x1为滑行段的状态变量,u1为滑行段的控制变量,由于滑行段发动机不开机,u1为0,和分别为滑行段的起止时间,
第二步:根据VTVL-RLV子级动力学模型,动力段的动力学flip-Radau伪谱离散方程为:
并且,两阶段间的连接条件为x1_f=x2_0;
第三步:在滑行段发动机控制量为零,仅考虑动力段的动压和热流约束,则滑行段的离散形式为:
动力段发动机推力约束的离散形式为:
问题的另外一个非凸性来源于问题非线性动力学模型和过程状态约束。本实施例利用序列凸化来解决非线性约束的凸化问题。序列凸化是将非线性约束进行线性化近似,进而得到近似凸问题的最优解,并通过序列迭代求解的方式使近似解不断逼近原非线性非凸问题的最优解。经过线性化处理,算法每一步迭代求解一个凸的子问题,并将每一次迭代得到的最优解作为下一次迭代的线性化参考点,直至最优解收敛。
由于火箭子级滑行段与动力段动力学模型相似,且动力段动力学较为复杂,首先给出动力段的动力学约束线性化形式,然后滑行段线性化动力学约束可类比得到。将动力段动力学约束重新列写如下,为表述简洁在此忽略表示不同飞行阶段的下标。步骤三所述对离散的轨迹优化问题进行凸化处理的具体过程包括:
步骤1:引入新的控制变量Γ对原控制约束进行松弛,将问题中控制变量约束做如下变换:
Pmin≤Γ(τi)≤Pmax,(i=1,...,N2) (10)
步骤2:将动力段动力学约束重新列写,所述动力段动力学约束如下:
步骤3:公式(11)中的非线性相关项集中在增广的动力学方程右端函数farg中,对公式(11)中的farg进行一阶泰勒展开,所述farg的一阶泰勒展开式如下:
其中,下标i表示在第i个配点处的约束,xk和为第k次迭代时泰勒展开的参考点,h(x)是动力学中与地球自转相关的项;和分别为线性化系统的状态矩阵和控制矩阵,且有:
上式中,
矩阵中其他元素与上述元素类似,为增广动力学方程右端函数对各个状态变量和增广控制变量求一阶偏导的结果;
步骤4:所述farg的一阶泰勒展开式得出线性化的动力段动力学约束方程,所述线性化的动力段动力学约束方程如下:
步骤5:对过程状态约束进行线性化处理,获得线性化的过程状态约束方程如下:
其中,
c11=ρ0e-κ·HV,
上述线性化过程中一个重要的环节是泰勒展开参考点的选择,在序列凸化算法中,第一次迭代的参考点由粗选的初值提供,随后的迭代过程中参考点取为上一次迭代得到的最优解。
根据泰勒展开的特点,在序列凸化迭代过程中,只有优化变量在参考点附近取值时,线性化的动力学和状态约束才是对原非线性形式的良好近似。因此,在序列凸化算法中加入如下信赖域约束:
式中δx和分别为针对各个状态和控制变量设计的信赖域值。
步骤四所述序列凸化算法的具体过程为:
Step1:输入给定初始优化变量x1,给定迭代停止条件以及最大迭代次数kmax,其中,最大迭代次数kmax的取值范围为10至20次。
Step2:求解上述凸化子问题,得到最优解Δxk;如果求解结束,返回成功标志与解xk;如果k>kmax,结束并返回失败标志;否则,进入step3;
Step3:置xk+1:=xk+Δx k:=k+1,返回Step 1。
虽然本发明已以较佳的实施例公开如上,但其并非用以限定本发明,任何熟悉此技术的人,在不脱离本发明的精神和范围内,都可以做各种改动和修饰,因此本发明的保护范围应该以权利要求书所界定的为准。
Claims (6)
1.一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,所述优化方法包括:
步骤一:建立垂直起降重复使用运载器返回着陆轨迹优化模型;
步骤二、对所述返回着陆轨迹优化模型进行离散化处理;
步骤三:对离散的轨迹优化问题进行凸化处理;
步骤四:建立保证收敛性能的序列凸化算法。
2.根据权利要求1所述基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,步骤一所述建立垂直起降重复使用运载器返回着陆轨迹优化模型的具体过程如下:针对垂直起降重复使用运载器返回着陆问题,建立以燃料最省为性能指标的轨迹优化模型,所述轨迹优化模型为:
min J=-mf (1)
其中,V为速度,γ为航迹倾角,ψ为航向角,φ为纬度,θ为经度,r为地心距,m为质量,Px、Py和Pz为火箭发动机推力在箭体坐标系下的分量,D为子级受到的空气阻力,fM为地球引力常数,Ω为地球自转角速度,PTotal为发动机总推力,Isp为发动机比冲,g0为海平面重力加速度,εH、εV、εγ、εφ、εθ分别为终端高度、速度、倾角和经纬度容许误差,mdry为子级干重;J为性能指标函数;mf为终端质量;ρ为大气密度;qmax和分别表示容许最大动压和热流密度;kQ为与飞行器外形相关的常数;Hf、Vf、γf、φf、θf、φd和θd分别表示终端状态值。
3.根据权利要求1所述基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,步骤二所述离散化处理过程中,火箭子级的滑行段飞行和动力段飞行在基于伪谱法的离散优化模型建立过程中分为两个阶段,在两个阶段内分别对问题进行伪谱离散并设置阶段间连接条件,所述离散化处理的具体过程包括:
第一步:根据VTVL-RLV子级动力学模型,滑行段的动力学flip-Radau伪谱离散方程为:
其中,D为flip-Radau微分矩阵,f(x,u)为微分方程的右端函数,τi,(i=1,...,N1)为(-1,1]区间内的配点,N1为配点数量,τ0为在-1处的离散点,x1为滑行段的状态变量,u1为滑行段的控制变量,由于滑行段发动机不开机,u1为0,和分别为滑行段的起止时间,
第二步:根据VTVL-RLV子级动力学模型,动力段的动力学flip-Radau伪谱离散方程为:
并且,两阶段间的连接条件为x1_f=x2_0;
第三步:在滑行段发动机控制量为零,仅考虑动力段的动压和热流约束,则滑行段的离散形式为:
动力段发动机推力约束的离散形式为:
其中,Cq、分别表示离散化的约束函数;r2、V2分别表示动力段的地心距与速度。
4.根据权利要求1所述基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,步骤三所述对离散的轨迹优化问题进行凸化处理的具体过程包括:
步骤1:引入新的控制变量Γ对原控制约束进行松弛,将问题中控制变量约束做如下变换:
Pmin≤Γ(τi)≤Pmax,(i=1,...,N2) (10)
其中,Pmin和Pmax分别表示最小和最大推力值;
步骤2:将动力段动力学约束重新列写,所述动力段动力学约束如下:
其中,farg表示增广的有端动力学函数,yi表示增广的优化变量;
步骤3:对公式(11)中的farg进行一阶泰勒展开,所述farg的一阶泰勒展开式如下:
其中,下标i表示在第i个配点处的约束,xk和为第k次迭代时泰勒展开的参考点,h(x)是动力学中与地球自转相关的项;H.O.T.表示高阶小项;和分别为线性化系统的状态矩阵和控制矩阵,且有:
上式中,
步骤4:所述farg的一阶泰勒展开式得出线性化的动力段动力学约束方程,所述线性化的动力段动力学约束方程如下:
步骤5:对过程状态约束进行线性化处理,获得线性化的过程状态约束方程如下:
其中,
5.根据权利要求1所述基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,步骤四所述序列凸化算法的具体过程为:
Step1:输入给定初始优化变量x1,给定迭代停止条件以及最大迭代次数kmax,其中,最大迭代次数kmax的取值范围为10至20次。
Step2:求解上述凸化子问题,得到最优解Δxk;如果求解结束,返回成功标志与解xk;如果k>kmax,结束并返回失败标志;否则,进入step3; 和分别表示约束残差和相应的收敛条件;
Step3:置xk+1:=xk+Δx k:=k+1,返回Step 1。
6.根据权利要求1或5所述基于凸优化的垂直起降重复使用运载器快速轨迹优化方法,其特征在于,步骤四所述序列凸化算法中,第一次迭代的参考点由粗选的初值提供,随后的每次迭代过程中参考点取为上一次迭代得到的最优解,并且,在所述凸化处理加入信赖域约束,所述信赖域约束为:
其中,δx和分别为针对各个状态和控制变量设计的信赖域值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811235049.XA CN109470252A (zh) | 2018-10-23 | 2018-10-23 | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811235049.XA CN109470252A (zh) | 2018-10-23 | 2018-10-23 | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109470252A true CN109470252A (zh) | 2019-03-15 |
Family
ID=65665976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811235049.XA Pending CN109470252A (zh) | 2018-10-23 | 2018-10-23 | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109470252A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109974538A (zh) * | 2019-03-25 | 2019-07-05 | 哈尔滨工业大学 | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 |
CN110989644A (zh) * | 2019-11-29 | 2020-04-10 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN111210668A (zh) * | 2019-12-30 | 2020-05-29 | 四川函钛科技有限公司 | 一种基于时序qar参数的着陆阶段飞行轨迹偏移修正方法 |
CN112093078A (zh) * | 2020-08-21 | 2020-12-18 | 哈尔滨工业大学 | 一种高精度高效率飞行器轨迹优化与制导切换方法 |
CN112149225A (zh) * | 2020-09-14 | 2020-12-29 | 南京航空航天大学 | 基于高精度离散格式的火星探测器定点着陆轨迹凸优化方法 |
CN112395689A (zh) * | 2020-11-19 | 2021-02-23 | 清华大学 | 基于凸优化的火箭故障后在线重构方法 |
CN112580188A (zh) * | 2020-11-12 | 2021-03-30 | 北京航天自动控制研究所 | 一种运载火箭圆轨道在线规划方法 |
CN112629339A (zh) * | 2020-12-15 | 2021-04-09 | 北京航天自动控制研究所 | 一种基于直接法的火箭软着陆轨迹规划方法 |
CN112693631A (zh) * | 2020-11-27 | 2021-04-23 | 中国人民解放军国防科技大学 | 飞行器在线序列凸优化中的初始轨迹生成方法及系统 |
CN113467498A (zh) * | 2021-07-14 | 2021-10-01 | 西北工业大学 | 一种基于Bezier-凸优化的运载火箭上升段轨迹规划方法 |
CN117826616A (zh) * | 2024-03-04 | 2024-04-05 | 西北工业大学 | 一种基于序列凸优化的飞行器快速轨迹规划方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120090458A1 (en) * | 2009-05-01 | 2012-04-19 | Floyd Brian A | Gas Strut Separation For Staged Rocket |
CN107092731A (zh) * | 2017-04-05 | 2017-08-25 | 北京临近空间飞行器系统工程研究所 | 亚轨道运载火箭级间比推力弹道一体化优化方法 |
CN108287476A (zh) * | 2018-02-06 | 2018-07-17 | 哈尔滨工业大学 | 基于高阶滑模控制和扰动观测器的空间翻滚非合作目标自主交会制导方法 |
-
2018
- 2018-10-23 CN CN201811235049.XA patent/CN109470252A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120090458A1 (en) * | 2009-05-01 | 2012-04-19 | Floyd Brian A | Gas Strut Separation For Staged Rocket |
CN107092731A (zh) * | 2017-04-05 | 2017-08-25 | 北京临近空间飞行器系统工程研究所 | 亚轨道运载火箭级间比推力弹道一体化优化方法 |
CN108287476A (zh) * | 2018-02-06 | 2018-07-17 | 哈尔滨工业大学 | 基于高阶滑模控制和扰动观测器的空间翻滚非合作目标自主交会制导方法 |
Non-Patent Citations (3)
Title |
---|
崔乃刚等: "垂直起降可重复使用运载器发展现状与关键技术分析", 《宇航总体技术》 * |
张志国等: "火箭垂直回收着陆段在线制导凸优化方法", 《弹道学报》 * |
王劲博等: "火箭返回着陆问题高精度快速轨迹优化算法", 《控制理论与应用》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109974538A (zh) * | 2019-03-25 | 2019-07-05 | 哈尔滨工业大学 | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 |
CN110989644A (zh) * | 2019-11-29 | 2020-04-10 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN110989644B (zh) * | 2019-11-29 | 2021-04-23 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN111210668A (zh) * | 2019-12-30 | 2020-05-29 | 四川函钛科技有限公司 | 一种基于时序qar参数的着陆阶段飞行轨迹偏移修正方法 |
CN112093078B (zh) * | 2020-08-21 | 2022-03-22 | 哈尔滨工业大学 | 一种高精度高效率飞行器轨迹优化与制导切换方法 |
CN112093078A (zh) * | 2020-08-21 | 2020-12-18 | 哈尔滨工业大学 | 一种高精度高效率飞行器轨迹优化与制导切换方法 |
CN112149225A (zh) * | 2020-09-14 | 2020-12-29 | 南京航空航天大学 | 基于高精度离散格式的火星探测器定点着陆轨迹凸优化方法 |
CN112149225B (zh) * | 2020-09-14 | 2024-04-30 | 南京航空航天大学 | 基于高精度离散格式的火星探测器定点着陆轨迹凸优化方法 |
CN112580188A (zh) * | 2020-11-12 | 2021-03-30 | 北京航天自动控制研究所 | 一种运载火箭圆轨道在线规划方法 |
CN112580188B (zh) * | 2020-11-12 | 2023-12-12 | 北京航天自动控制研究所 | 一种运载火箭圆轨道在线规划方法 |
CN112395689B (zh) * | 2020-11-19 | 2022-05-06 | 清华大学 | 基于凸优化的火箭故障后在线重构方法 |
CN112395689A (zh) * | 2020-11-19 | 2021-02-23 | 清华大学 | 基于凸优化的火箭故障后在线重构方法 |
CN112693631A (zh) * | 2020-11-27 | 2021-04-23 | 中国人民解放军国防科技大学 | 飞行器在线序列凸优化中的初始轨迹生成方法及系统 |
CN112693631B (zh) * | 2020-11-27 | 2022-07-29 | 中国人民解放军国防科技大学 | 飞行器在线序列凸优化中的初始轨迹生成方法及系统 |
CN112629339A (zh) * | 2020-12-15 | 2021-04-09 | 北京航天自动控制研究所 | 一种基于直接法的火箭软着陆轨迹规划方法 |
CN113467498A (zh) * | 2021-07-14 | 2021-10-01 | 西北工业大学 | 一种基于Bezier-凸优化的运载火箭上升段轨迹规划方法 |
CN113467498B (zh) * | 2021-07-14 | 2022-07-01 | 西北工业大学 | 一种基于Bezier-凸优化的运载火箭上升段轨迹规划方法 |
CN117826616A (zh) * | 2024-03-04 | 2024-04-05 | 西北工业大学 | 一种基于序列凸优化的飞行器快速轨迹规划方法及装置 |
CN117826616B (zh) * | 2024-03-04 | 2024-05-10 | 西北工业大学 | 一种基于序列凸优化的飞行器快速轨迹规划方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109470252A (zh) | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 | |
CN111428372B (zh) | 基于凸规划和自适应迭代的火箭动力故障降级入轨制导方法 | |
CN107544262B (zh) | 一种运载火箭自适应精确回收控制方法 | |
CN110989644B (zh) | 一种考虑目标点多终端约束的飞行器轨迹规划方法 | |
CN108180910B (zh) | 一种基于气动参数不确定的飞行器快速高精度制导方法 | |
CN101937233B (zh) | 近空间高超声速飞行器非线性自适应控制方法 | |
CN108646554B (zh) | 一种基于指定性能的飞行器快速抗干扰纵向制导方法 | |
CN108536020A (zh) | 一种针对垂直起降重复使用运载器的模型参考自适应滑模控制方法 | |
CN111306989A (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
CN112395689B (zh) | 基于凸优化的火箭故障后在线重构方法 | |
CN110736470A (zh) | 一种不规则形状小天体附近连续推力轨道混合搜索方法 | |
CN111191368B (zh) | 一种连续小推力行星际转移轨道优化方法和装置 | |
CN113900448B (zh) | 一种基于滑模干扰观测器的飞行器预测校正复合制导方法 | |
Bialy et al. | An adaptive backstepping controller for a hypersonic air-breathing missile | |
CN104567545A (zh) | Rlv大气层内主动段的制导方法 | |
Xu et al. | Multi-constrained ascent guidance for solid propellant launch vehicles | |
Peng et al. | Free return orbit design and characteristics analysis for manned lunar mission | |
Tan et al. | Optimal maneuver trajectory for hypersonic missiles in dive phase using inverted flight | |
CN113504723B (zh) | 一种基于逆强化学习的运载火箭减载控制方法 | |
Martin et al. | Saturn V guidance, navigation, and targeting. | |
CN105354380A (zh) | 面向摄动因素影响补偿的滑翔弹道快速修正方法 | |
CN113093776A (zh) | 一种航天器的离轨参数确定方法及装置 | |
CN109062044B (zh) | 一种终端迭代学习对接控制方法 | |
Zhou et al. | Ascent trajectory optimization for air‐breathing vehicles in consideration of launch window |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190315 |