[发明专利]一种计算飞机升力曲线的方法及装置

权利要求书

1.一种计算飞机升力曲线的方法，其特征在于，所述方法包括：

以滑动时间窗的方式，实时采集飞机飞行时的预设时间段内的至少一组与升力曲线相关的第一参数；

分别将所述至少一组中每一组与升力曲线相关的第一参数代入第一公式中，计算与升力曲线相关的第二参数；

所述第二参数包括升力曲线协方差函数表达式中的参数σ、λ和n，以及均值函数表达式中的升力曲线水平偏移因子a、升力曲线竖直偏移因子b、升力曲线缩放比例因子c；

所述第一公式的表达式由下式所示：

其中，arg函数表示求最优解；m(A)表示迎角为A时的均值向量；k(A,A)表示迎角为A时的协方差矩阵；C为升力系数，其中，A和C均为已知数据；

根据所述第二参数和第一升力曲线函数表达式，计算飞机飞行时的预设时间段内的第二升力曲线函数表达式，其中，所述第一升力曲线函数表达式为飞机飞行状态为预设状态时的升力曲线函数表达式；

所述第二升力曲线函数表达式分别由以下公式表示：

均值函数表达式：

m_p(α)＝m(α)+k(α,A)k(A,A)^-1[C-m(A)]

协方差函数表达式：

k_p(α,α’)＝k(α,α’)-k(α,A)k(A,A)^-1k(A,α’)

其中，m_p(α)为所述第二升力曲线函数表达式中的均值函数表达式，m(α)为与第二升力曲线函数表达式对应的公式模型中的均值函数，k(α,A)k(A,A)^-1[C-m(A)]为校正参数；k_p(α,α’)为所述第二升力曲线函数表达式中的协方差函数表达式，k(α,α’)为与第二升力曲线函数表达式对应的公式模型中的协方差函数；

所述与第二升力曲线函数表达式对应的公式模型分别由以下两个公式表示：

均值函数表达式为：m(α)＝c*C_Lm(α+a)+b；

协方差函数表达式为：k(α,α’)＝σ*exp[(α-α’)²/(2λ²)]+δ_αα’n²；

其中，α和α’分别为所述预设时间段内不同时刻的迎角；C_Lm(α)为第一升力曲线函数表达式；δ_αα’为系数，且当α＝α’时，δ_αα’为1，否则，δ_αα’为0；

根据所述第二升力曲线函数表达式，生成第二升力曲线，其中所述第一公式为计算所述第二参数最优解的公式，所述至少一组中每一组与升力曲线相关的第一参数均包括与所述升力曲线对应的升力系数和迎角。

2.根据权利要求1所述的方法，其特征在于，所述根据所述第二参数和第一升力曲线函数表达式，计算飞机飞行时的预设时间段内的第二升力曲线函数表达式，具体包括：

根据所述第二参数和所述第一升力曲线函数表达式，获取与第二升力曲线函数表达式对应的公式模型；

根据所述第二参数和所述第二升力曲线函数表达式对应的公式模型，获取所述第二升力曲线函数表达式。

3.根据权利要求2所述的方法，其特征在于，所述第一升力曲线函数表达式、所述与第二升力曲线函数表达式以及与第二升力曲线函数表达式对应的公式模型均为高斯过程模型的表达式，所述高斯过程模型的表达式具体包括：均值函数表达式和协方差函数表达式。