本发明专利技术公开了一种基于可变数据长度最大信息量-可信度准则的飞行器建模方法,用于解决现有的最大信息量准则不考虑数据长度而导致飞行试验给出的气动模型和参数验证正确性差的技术问题。技术方案是通过在最大信息量准则考虑数据长度并且引入可信度参数,可以根据飞行器的不同飞行试验因素修正建模准则;对测量方差估计Rj和Rj+1的U-D分解,得到了标量模型选择和验证判别式。便于直接根据飞行试验数据建立飞行器气动力、力矩模型,避免了最大信息量准则未直接考虑数据长度导致用不同飞行试验数据建立和验证气动模型不正确的技术问题。
【技术实现步骤摘要】
本专利技术涉及一种飞行器建模方法,特别是涉及一种。
技术介绍
根据飞机气动模型和参数不仅可以确定飞机的操纵稳定性,还可为地面和空中仿真器提供正确的数学模型;验证飞机气动参数的风洞实验和理论计算结果;为飞机控制系统的设计和改进提供基本数据;鉴定定型飞机的飞行品质;研究高性能飞机的飞行品质; 进行飞机失事的事故分析等等;准确地建立飞机数学模型问题与通过基本定律、定理等机理建模的理论方法截然不同,主要根据实验所得的输入和输出数据建立模型,其基本理论依据为非线性系统辨识学和非线性飞行动力学;当飞机作小迎角小扰动飞行时,气动力和力矩模可以用台劳级数展开取一次项,即Bryan模型表示。当马赫数、高度一定时,这一模型是线性定常模型,此模型因为形式简单而一直沿用至今,成为气动数学模型的基石;采用这种模型,飞行器系统辨识就成了对已知数学模型的系统参数估计了 ;现代战斗机、战术导弹在作战时需要较大机动、过失速甚至尾旋,其迎角可以从十几度、几十度直至一百多度,已不能采用线性模型;飞机大迎角形成的脱体涡、分离涡所引起的非定常下洗流场、使得定常模型也不能再适用了。研究在大迎角下飞行器的非定常、非线性气动模型已成为当前飞机研制的迫切需要的问题。然而,非线性气动力的辨识异常复杂,它是一般的非线性系统辨识问题,输入量与状态之间的函数关系很难确定,需要对模型进行辨识;模型辨识的关键是建模判据和优选算法,对于给定的结构形式,应用建模判据来确定模型的最优阶数并从侯选模型中选出最优模型;由于实测数据含有噪声,建模判据不能仅仅考察对现有数据的拟合误差大小,而且综合考虑其它因素,否则将会使模型不正确;通常,建模判据应能使优选出的模型具有以下特点1.模型很好地拟合现有飞行数据;2.模型各项有明显的物理意义;3.模型能预测类似条件下的实测数据;4.在性能相当的条件下阶次最低;最常用的模型辨识方法是逐步回归法,其原理是逐项将影响显著性的预报因子选入,并将影响小的因子剔除,建立回归方程的方法;这一方法计算简单、比较实用;但这一方法有两个明显的缺点一是选择标准由人而定,而且没有给出结果的可信度;二是误差积累大,容易漏选和误选;为此,人们对要求较高的飞行器模型辨识问题常常采用最大信息量准则AIC方法,但是该方法处理速度慢,信噪比较小时模型辨识精度差;由于在非线性情况下,只能对飞机非线性方程进行数值积分,进行灵敏度矩阵计算和迭代计算,从而使计算的复杂程度和计算量比线性估计高得多,同时也使模型输出与实验数据之间的拟合变得更加困难,特别是当飞行数据长度不同时,现有AIC准则没有直接考虑不同的数据长度,常常会导致飞行试验给出的气动模型和参数验证不正确
技术实现思路
为了克服现有最大信息量准则不考虑数据长度而导致飞行试验给出的气动模型和参数验证正确性差的不足,本专利技术提供一种。该方法通过分析数据长度的影响,并且引入可信度参数对最大信息量准则进行了修正,得到了新的模型辨识判据,由新判据建立了指数建模,直接可以用于飞行器的飞行试验建模和模型验证,可以避免根据飞行试验建立和验证飞行器大迎角模型存在的技术问题。本专利技术解决其技术问题所采用的技术方案是一种,其特点是包括以下步骤步骤一、飞行试验待确定的飞行器候选模型的状态方程为x(t) = f{fo, fjx ⑴,,X ⑴], (I)观测方程为 「 n iy(0 = g = g{g0,gj,fi, t} > g{g0, g!,t}为表达式已知的待确定模型结构函数,fJxU),Q0]>g(l为根据物理概念必须选入的模型,fi >gi (i = 1,2, ···, q)为候选模型,ζ (tk)为在 tk 时亥Ij对y(tk)的测量值;Ω为未知维数的参数向量,Qci为已知维数的参数向量;V(k)为测量噪声,假定方差为Rk的零均值高斯白噪声;fi-1η/-,Ν,),Ω丨'/]>&⑶ In/.,式中,LJx⑴,Q1, t]、L2为取不同的参数向量和02的极大似然函数,S为给定正实数、表示可信度参数,mN ( N、f 2π λ I Ν<2 ItiJ 2 L KmNJ J 2tf⑷InZ—2-手 1η ^1 + 1,η2(1) = z(tk) ~ g,R] # U^ ’R2,Ω i和Ω2为不同的参数向量,N为数据长度,In为自然对数符号; N k=l步骤二、根据假定fQ、Ωα=Ω。已经通过优选算法选入模型,并由以下算法迭代计算得到令j=0,l,2, ···, q,假定 fj > gj 、Ω j 已经选入模型,按照以下方式选择其它候选模型求(4)式极大值,迭代计算ΔΩ; = A^bj-(5)以及RiIii W = z(4)—g(6) jY dy _(5),(6)式中 A bi = Σ (^7)Τ R/ > ZjJa^ ~ j — j, k=i o%£tj「ml * -V1/_ p Tp-In pT—|γ 办(ΟγΓ /^(^2)\Τ bA bA —P1 =diag, R}1, ··· R1J 「Ω, 当两个飞行试验数据长度分别为N、Μ、设Ω ,= J ,Θ j+1的选入或剔除及验证 LU条件为θ j+1的选入或剔除模型验证条件为当:nlR v(7<r)| Η,、。ο(人:〉谷及_£(υ)+1η^Μ^ <ε(7) |/ (!n2^+]) + lii|RA;(/:)| 2、 |R.w.⑷| ΓΩ, 时,Θ」+1、fj+1、gj+1选入模型,且Ω , = ;否则剔除 LtJ'+! _fj+1 、gJ+1 候选项,且 Ω J+1= Ω j ;(7)式中ε 为给定正数,Rjw,RaX間二古IX+iWC丨W ’Rmj =去Σ #用(k),R_+1)=去, 皿 k^lM k^iη, W = <tk) -g ,r\J+l{k) = z(tk) — g[x(4),Ω/+1 ,Ik);步骤三、由于飞行器测量向量y的维数m较大,采用Gram-Schmidt正交化方法对 Rn」·、Rmj·、Rn(j+i) 和Rm(w)进行U-D分解,Rnj·、Rmj·、^NCj+l)和RM(j+1)的U-D分解分另1J 为认 Nj - V RNjB RNj^ RNj ,— U mf} D R%{j U mff , R^v(j+1) — Ujivi 州)D^v(汁咖J+1),R = 11ι\ [1ivAfi j+1) ~ RM(j+l)^iai(j+l) RAf ( f+l),式中,Uenj-λ Uemj- λ UEN(j+1) λ Uem(j+1)为单位上二角阵;Denj--diag, Drn(j_+1)-diagEd·。·+。(I),dRN(j_+1) (2),…,c^RN (j+1) (m)],D—j-diag,DRM(j+i)_diag[dEM(j+1) (I) dRM(j+1) (2),…,c^RM (j+1) (m)];diag为对角符号;模型验证的最大信息量准则写成当权利要求1.一种,其特征在于包括以下步骤 步骤一、飞行试验待确定的飞行器候选模型的状态方程为全文摘要本专利技术公开了一种,用于解决现有的最大信息量准则不考虑数据长度而导致飞行试验给出的气动模型和参数验证正确本文档来自技高网...
【技术保护点】
一种基于可变数据长度最大信息量?可信度准则的飞行器建模方法,其特征在于包括以下步骤:步骤一、飞行试验待确定的飞行器候选模型的状态方程为x·(t)=f{f0[x(t),Ω0],f1[x(t),θ1],...,fq[θq,x(t)],t}---(1)观测方程为y(t)=g[x(t),Ω,t]=g{g0[x(t),Ω0],g1[x(t),θ1],...,gq[θq,x(t)],t}z(tk)=y(tk)+v(k)---(2)(1)、(2)式中,x(t)为n维状态向量;y(t)为m维观测向量;f{f0[x(t),Ω0],f1[x(t),θ1],...,fq[θq,x(t)],t}、g{g0[x(t),Ω0],g1[x(t),θ1],...,gq[θq,x(t)],t}为表达式已知的待确定模型结构函数,f0[x(t),Ω0]、g0[x(t),Ω0]为根据物理概念必须选入的模型,fi[x(t),θi]、gi[x(t),θi](i=1,2,…,q)为候选模型,z(tk)为在tk时刻对y(tk)的测量值;Ω为未知维数的参数向量,Ω0为已知维数的参数向量;v(k)为测量噪声,假定方差为Rk的零均值高斯白噪声;fi[x(t),θi]、gi[x(t),θi](i=1,2,…,q)是否在模型中出现及Ω0、θi(i=1,2,…,q)的取值需要辨识,q为已知的候选模型个数;由于对飞行器的模型结构准确度要求较高,给出以下最大信息量?可信度准则:|lnL2[x(t),Ω2,t]-lnL1[x(t),Ω1,t]lnL1[x(t),Ω1,t]|>δ,---(3)式中,L1[x(t),Ω1,t]、L2[x(t),Ω2,t]为取不同的参数向量Ω1和Ω2的极大似然函数,δ为给定正实数、表示可信度参数,lnL1[x(t),Ω1,t]=-mN2ln(Σk=1Nη‾1T(k)R1-1(k)η‾1(k))-mN2[ln(2πmN)+1]-12Σk=1Nln|R1(k)|lnL2[x(t),Ω2,t]=-mN2ln(Σk=1Nη‾2T(k)R2-1(k)η‾2(k))-mN2[ln(2πmN)+1]-12Σk=1Nln|R2(k)|---(4),η‾1(k)=z(tk)-g[x(tk),Ω1,tk],η‾2(k)=z(tk)-g[x(tk),Ω2,tk],R1=1NΣk=1Nη‾1(k)η‾1T(k),Ω1和Ω2为不同的参数向量,N为数据长度,1n为自然对数符号;步骤二、根据假定f0[x(t),Ω0]、g0[x(t),Ω0]、Ω0=Ω0已经通过优选算法选入模型, 并由以下算法迭代计算得到:令j=0,1,2,…,q,假定fj[x(t),θj]、gj[x(t),θj]、Ωj已经选入模型,按照以下方式选择其它候选模型:求(4)式极大值,迭代计算:ΔΩj=Aj-1bj---(5)以及Rj=1NΣk=1Nη‾j(k)η‾jT(k),η‾j(k)=z(tk)-g[x(tk),Ωj,tk]---(6)(5)、(6)式中:ΔΩj=Ωj-Ω^j,bj=Σk=1N(∂y∂ΩjT)TRj-1[z(tk)-y(tk)],Aj=Σk=1N(∂y∂ΩjT)TRj-1∂y∂ΩjT=BjTPj-1Bj,BjT=[(∂y(t1)∂ΩjT)T,(∂y(t2)∂ΩjT)T,...,(∂y(tN)∂ΩjT)T]Pj-1=diag[Rj-1,Rj-1,...Rj-1],当两个飞行试验数据长度...
【技术特征摘要】
【专利技术属性】
技术研发人员:史忠科,
申请(专利权)人:西北工业大学,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。