航天器抓捕目标后组合体无模型快速消旋稳定控制方法技术

本发明专利技术公开了一种航天器抓捕目标后组合体无模型快速消旋稳定控制方法，在不需要组合体动力学模型信息的基础上，实现对组合体的低复杂度鲁棒控制器设计。同时，组合体的瞬态与稳态性能能够得到先验设计和保障。本发明专利技术在不需要对系统未知动力学模型在线辨识以及未知外界干扰观测条件下，就可以实现对服务星抓捕后组合体的快速消旋与稳定控制，控制系统设计的复杂度大大降低；本发明专利技术实现了对组合体控制系统的瞬态态与稳态性能的先验设计，且组合体控制系统的性能不依赖于大量繁复的控制器调参过程，控制系统设计简单，易于实际工程应用。

【技术实现步骤摘要】

本专利技术属于航天器控制
，涉及一种航天器抓捕目标后组合体无模型快速消旋稳定控制方法。
技术介绍
近几十年，在轨服务技术取得了突飞猛进的发展，但是随着空间卫星发射数量的不断增加，由于各种故障和燃料的消耗殆尽，失效报废的卫星数量也越来越多。考虑到失效报废卫星上仍存在大量可用的部件，为了延长其寿命，通常借助于服务卫星，对失效报废卫星进行接管控制。在对失效报废卫星的接管控制中，存在两个亟待解决的问题：一是服务卫星在对失效卫星抓捕后形成组合体的未知动力学模型的辨识；二是抓捕后组合体的快速消旋与稳定控制。其中第一个问题往往是第二个问题中控制器设计的关键。传统针对组合体动力学模型辨识的方法有基于组合体动力学的最小二乘拟合方法、基于神经网络和模糊系统的智能辨识方法。但是辨识算法的可靠性以及复杂度影响到后续控制系统设计的时效性，因此传统针对未知组合体动力学模型辨识的方法对于后续控制系统设计具有较大的影响。对于组合体快速消旋和稳定控制方法上，传统方法有基于模型的状态依赖Riccati方法、基于神经网络和模糊辨识的自适应控制方法、基于模型的预测控制方法等,～。传统方法多依赖于组合体的动力学模型，因此在模型不精确或者模型未知情况下，控制系统的性能受限较大。与此同时，组合体系统的瞬态与稳态性能难以得到保障。针对传统组合体控制系统设计存在的缺点，需要新的控制策略和控制方法，在不需要对组合体未知动力学模型辨识基础上，对组合体进行快速消旋与稳定控制，从而为实现后续组合体的其他任务提供充足的时间和良好的操作条件。
技术实现思路
本专利技术的目的在于解决上述现有技术中的问题，提供一种航天器抓捕目标后组合体无模型快速消旋稳定控制方法，该方法针对服务星抓捕失效空间目标后形成的组合体，在不需要对包含未知动力学模型的组合体辨识基础上，设计低复杂度的鲁棒控制器，实现对组合体的快速消旋与稳定控制。与此同时，组合体控制系统的瞬态与稳态性能能够得到先验设计。为达到上述目的，本专利技术采用以下技术方案予以实现：航天器抓捕目标后组合体无模型快速消旋稳定控制方法，包括以下步骤：步骤一：组合体姿态动力学模型针对抓捕后形成的组合体姿态动力学模型为：其中：姿态采用修正罗格里格斯参数—MRPs，σ＝[σ1,σ2,σ3]T∈R3,ω∈R3分别为在MRPs表示下的姿态角和角速度；J∈R3×3为组合体惯性坐标系的转动惯量；Hw＝diag(Jw1,Jw2,Jw3,Jw4),Ω＝[Ω1,Ω2,Ω3,Ω4]T分别为四个反作用轮的转动惯量和转速；C∈R3×4为四个反作用轮的安装矩阵；τext∈R3为未知干扰力矩；σ×是反对称矩阵，具体形式为：步骤二：姿态动力学模型转换定义χ＝Γ(σ)ω，对式(1)进行简化得：其中：v＝[σT,χT]T；Π是对角正定矩阵；uc:＝-Cuw为虚拟控制力；f(ω)＝-(J-CHwCT)-1ω×(Jω+CHwΩ),G(v)＝Γ(σ)g,g＝(J-CHwCT)-1,其中参数F(·),G(·),d(·)为未知参数；假定输出跟踪轨迹为yr，则输出误差为ε＝[ε1,ε2,ε3]T＝y-yr∈R3，定义如下性能指标：-δi,1μi(t)＜εi(t)＜δi,2μi(t)(4)其中：δi,1,δi,2为待设计常值参数；μi(t)＞0为严格递减函数，并且取为μi0＞μi∞＞0,κi＞0；在预设性能(4)下，为了降低控制系统设计复杂度，定义εi(t):＝μi(t)P(zi)，且P(zi)取为：则有对新定义的转化误差zi求导得其中：z＝[z1,z2,z3]T.μ＝[μ1,μ2,μ3]T,ξ＝diag(ξ1,ξ2,ξ3),Λ＝diag(Λ1,Λ2,Λ3).步骤三：无模型鲁棒控制器设计设计的无模型鲁棒控制器为其中：k＝diag{k1,k2,k3本文档来自技高网...

【技术保护点】
航天器抓捕目标后组合体无模型快速消旋稳定控制方法，其特征在于，包括以下步骤：步骤一：组合体姿态动力学模型针对抓捕后形成的组合体姿态动力学模型为：σ·=Γ(σ)ωω·=-(J-CHwCT)-1[ω×(Jω+CHwΩ)+Cuw-τext]Ω·=Hw-1uw-CTω·---(1)]]>其中：姿态采用修正罗格里格斯参数—MRPs，σ＝[σ1,σ2,σ3]T∈R3,ω∈R3分别为在MRPs表示下的姿态角和角速度；J∈R3×3为组合体惯性坐标系的转动惯量；Hw＝diag(Jw1,Jw2,Jw3,Jw4),Ω＝[Ω1,Ω2,Ω3,Ω4]T分别为四个反作用轮的转动惯量和转速；C∈R3×4为四个反作用轮的安装矩阵；τext∈R3为未知干扰力矩；Γ(σ)＝1/4[(1‑σTσ)I3+2σ×+2σσT].σ×是反对称矩阵，具体形式为：σ×=0-σ3σ2σ30-σ1-σ2σ10---(2)]]>步骤二：姿态动力学模型转换定义χ＝Γ(σ)ω，对式(1)进行简化得：σ·=χχ·=F(v)+G(v)uc+dy=Πσ---(3)]]>其中：v＝[σT,χT]T；Π是对角正定矩阵；uc:＝‑Cuw为虚拟控制力；G(v)＝Γ(σ)g,g＝(J‑CHwCT)‑1,其中参数F(·),G(·),d(·)为未知参数；假定输出跟踪轨迹为yr，则输出误差为ε＝[ε1,ε2,ε3]T＝y‑yr∈R3，定义如下性能指标：‑δi,1μi(t)＜εi(t)＜δi,2μi(t)   (4)其中：δi,1,δi,2为待设计常值参数；μi(t)＞0为严格递减函数，并且取为μi0＞μi∞＞0,κi＞0；在预设性能(4)下，为了降低控制系统设计复杂度，定义εi(t):＝μi(t)P(zi)，且P(zi)取为：P(zi)=δi,2ezi-δi,1e-ziezi+e-zi---(5)]]>则有对新定义的转化误差zi求导得z·=ξ(Πσ·-y·r-Λμ·)---(6)]]>其中：z＝[z1,z2,z3]T.μ＝[μ1,μ2,μ3]T,ξ＝diag(ξ1,ξ2,ξ3),Λ＝diag(Λ1,Λ2,Λ3).步骤三：无模型鲁棒控制器设计设计的无模型鲁棒控制器为其中：k＝diag{k1,k2,k3},η＝diag{η1,η2,η3}是正定对角矩阵；其中i＝1,2,3；步骤四：控制力矩鲁棒分配虚拟控制力与四个反作用轮之间的关系为uc＝‑Cuw＝‑(C0+ΔC0)uw   (8)其中：C0,ΔC0∈R3×4分别为反作用轮标称安装矩阵和偏差矩阵；则控制力矩鲁棒分配的优化问题为：其中：θmin,θmax,c#分别为反作用轮的控制上下界，以及偏差矩阵的范数上界；为Lagrangian乘子；Q为待设计正定矩阵；对(9)式进行内外最值进行优化，则实现控制力矩的鲁棒分配。...

【技术特征摘要】
1.航天器抓捕目标后组合体无模型快速消旋稳定控制方法，其特征在于，包括以下步骤：步骤一：组合体姿态动力学模型针对抓捕后形成的组合体姿态动力学模型为：σ·=Γ(σ)ωω·=-(J-CHwCT)-1[ω×(Jω+CHwΩ)+Cuw-τext]Ω·=Hw-1uw-CTω·---(1)]]>其中：姿态采用修正罗格里格斯参数—MRPs，σ＝[σ1,σ2,σ3]T∈R3,ω∈R3分别为在MRPs表示下的姿态角和角速度；J∈R3×3为组合体惯性坐标系的转动惯量；Hw＝diag(Jw1,Jw2,Jw3,Jw4),Ω＝[Ω1,Ω2,Ω3,Ω4]T分别为四个反作用轮的转动惯量和转速；C∈R3×4为四个反作用轮的安装矩阵；τext∈R3为未知干扰力矩；Γ(σ)＝1/4[(1-σTσ)I3+2σ×+2σσT].σ×是反对称矩阵，具体形式为：σ×=0-σ3σ2σ30-σ1-σ2σ10---(2)]]>步骤二：姿态动力学模型转换定义χ＝Γ(σ)ω，对式(1)进行简化得：σ&...

【专利技术属性】
技术研发人员：罗建军，魏才盛，袁建平，王明明，朱战霞，殷泽阳，夏鹏程，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61