气动力多源数据融合和辨识方法、系统及介质技术方案

技术编号:25709238 阅读:46 留言:0更新日期:2020-09-23 02:55
本发明专利技术提供了一种气动力多源数据融合和辨识方法、系统及介质,包括:步骤1:建立全局相关的气动力数学模型;步骤2:根据现有计算数据、风洞试验数据和飞行试验数据,计算多种来源数据的标准差;步骤3:将现有计算数据、风洞试验数据和飞行试验数据作为气动力数学模型的数据样本,建立超定方程组;步骤4:采用加权最小二乘法求解超定方程组,权值由标准差计算获得,解出气动力数学模型中的系数,完成多源数据融合气动力数学模型的建立。本发明专利技术解决了气动力数学模型辨识和修正中不同来源数据可靠度不同及误差分解到样本点的问题,可对计算、风洞试验和飞行试验三个来源的数据进行综合求解超定方程组,一步完成辨识和模型修正。

【技术实现步骤摘要】
气动力多源数据融合和辨识方法、系统及介质
本专利技术涉及气动设计
,具体地,涉及一种气动力多源数据融合和辨识方法、系统及介质。
技术介绍
随着飞行器外形复杂度增加和飞行包线的扩大,沿袭传统的分段函数作为气动力数学模型的建模方法,难以科学地平衡不同来源数据的不同可靠度;同时,将飞行试验中非建模样本点数据误差分解到样本点上,存在分解难题,经常采用的贝叶斯方法被认为科学性不够。大量的工程师们对气动力数学模型的建模方法进行了持之以恒的努力和发展,以便形成连续函数形式的气动力数学模型的建模方法。专利文献CN108318203A(申请号:201810172976.5)公开了飞行器等价地面颤振综合试验方法,该方法通过分析不同温度、飞行速度、大气密度、气流环境对飞行器所受弹性气动力耦合的影响,建立非线性偏微分方程描述的飞行器弹性气动力模型;在此模型基础上,根据已有试验数据分析,优化信号激励/拾振点布局,采用偏微分方程反卷积方法,设计获取气动力输入信号理论值,采用多点分布激励力控制方法,有效跟踪输入信号理论值,精确控制多点激振力,实现多点激励力模拟气动力的效果;建立结构温度试验覆盖大范围温度的环境,采用专用多点非接触高速图像与接触式传感器一体化融合检测方法,解决了接触式传感器无法精确检测超低温、高温情况下结构振动幅度、频率的问题。
技术实现思路
针对现有技术中的缺陷,本专利技术的目的是提供一种气动力多源数据融合和辨识方法、系统及介质。根据本专利技术提供的气动力多源数据融合和辨识方法,包括:r>步骤1:建立全局相关的气动力数学模型;步骤2:根据现有计算数据、风洞试验数据和飞行试验数据,计算这三种来源数据的标准差;步骤3:将现有计算数据、风洞试验数据和飞行试验数据作为气动力数学模型的数据样本,建立超定方程组;步骤4:采用加权最小二乘法求解超定方程组,权值由标准差计算获得,解出气动力数学模型中的系数,完成多源数据融合气动力数学模型的建立。优选的,所述步骤1中的全局相关的气动力数学模型的表达式为:其中,f(Ma,αΦ,Φ,δP,δY,δR)表示气动力系数或气动力矩系数,包括法向力系数、俯仰力矩系数、侧向力系数、偏航力矩系数、滚转力矩系数和轴向力系数,轴向力系数不包括摩阻和底阻部分;Ma为来流马赫数;αΦ为合成攻角,速度矢量与体轴的夹角;Φ为气流滚转角,速度矢量在飞行器横截面的分量与选定子午面的夹角;δP为俯仰通道舵偏角;δY为偏航通道舵偏角;δR为滚转通道舵偏角;φ为初始相位角;表示对i、j、k、l、r和s所有取值范围进行求和;aijklrs为常系数,i、j、k、l均为0,1,2,…;Tr(x)为马赫数相关的切比雪夫多项式,Ts(y)为合成攻角相关的切比雪夫多项式。优选的,所述Tr(x)计算式为:Tr(x)=cos(rcos-1(x)),r=1,2,…;优选的,所述Ts(y)计算式为:Ts(y)=cos(scos-1(y)),s=1,2,…;其中下标max表示最大值,下标min表示最小值。优选的,对i、j、k、l、r和s进行组合,选取m项基函数E={e1,e2,…,em},全局相关的气动力数学模型写为:其中,ai表示数学模型各项系数;ei表示m项基函数E的第i项基。优选的,所述步骤3中的超定方程组表达式为:其中,n为现有计算数据、风洞试验数据和飞行试验数据三个来源的数据总量。优选的,写成矩阵形式为:fn×1=am×nen×1;即F=AE;F表示fn×1;A表示am×n;E表示en×1。优选的,所述步骤4中的加权最小二乘法求解超定方程组后,系数矩阵A为:A=(JTWJ)-1JTWF;其中,J为雅克比矩阵;W为权函数矩阵,公式为:wi为权函数,i=1,2,…,n,计算公式为:σi为现有计算数据、风洞试验数据和飞行试验数据三个来源数据的标准差,i=1,2,…,n;其中,上标T表示矩阵转置,上标-1表示矩阵求逆。根据本专利技术提供的气动力多源数据融合和辨识系统,包括:模块M1:建立全局相关的气动力数学模型;模块M2:根据现有计算数据、风洞试验数据和飞行试验数据,计算这三种来源数据的标准差;模块M3:将现有计算数据、风洞试验数据和飞行试验数据作为气动力数学模型的数据样本,建立超定方程组;模块M4:采用加权最小二乘法求解超定方程组,权值由标准差计算获得,解出气动力数学模型中的系数,完成多源数据融合气动力数学模型的建立。与现有技术相比,本专利技术具有如下的有益效果:1、本专利技术解决了气动力数学模型辨识和修正中不同来源数据可靠度不同及误差分解到样本点的问题,可对计算、风洞试验和飞行试验三个来源的数据进行综合求解超定方程组,一步完成辨识和模型修正;2、本专利技术随着研制的深入和飞行数据的积累,将新的飞行试验数据持续增加进来辨识数据、更新和修正气动力数学模型,从而不断提高数学模型精度。附图说明通过阅读参照以下附图对非限制性实施例所作的详细描述,本专利技术的其它特征、目的和优点将会变得更明显:图1为本专利技术方法的步骤流程图。具体实施方式下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本专利技术,但不以任何形式限制本专利技术。应当指出的是,对本领域的普通技术人员来说,在不脱离本专利技术构思的前提下,还可以做出若干变化和改进。这些都属于本专利技术的保护范围。本专利技术要解决的技术问题是:针对具有复杂气动特性飞行器的全局相关气动力数学模型,基于气动数据来源多样性的事实,将现有计算数据、风洞试验数据和飞行试验数据共同作为建模的样本点,建立超定方程组,采用加权最小二乘法求解,以权值平衡不同来源气动数据的可靠度,从而获得新的气动力数学模型。本专利技术可随着研制阶段的不断推进,将新的飞行试验数据持续增加进来辨识数据、更新和修正气动力数学模型,不断提高气动力数学模型的精度。针对具有复杂气动特性飞行器的全局相关气动力数学模型,基于气动数据来源多样性的事实,将现有计算数据、风洞试验数据和飞行试验数据共同作为建模的样本点,建立超定方程组,采用加权最小二乘法求解,以权值平衡不同来源气动数据的可靠度,从而获得多源数据融合气动力数学模型。在不增加飞行试验获取信息的基础上,本专利技术提供了一种全局相关、精确的气动力多源数据融合和辨识方法。根据本专利技术提供的气动力多源数据融合和辨识方法,如图1所示,包括如下步骤:S1、建立全局相关的气动力数学模型表达式:全局相关的气动力数学模型的一般表达式为:式中,f(Ma,αΦ,Φ,δP,δY,δR)为气动力系数或气动力矩系数,包括法向力系数、俯仰力矩系数、轴向力系数(不包括摩阻和底阻部分)、侧向力系数、偏航力矩系数和滚转力矩系数;Ma为来流马赫数本文档来自技高网...

【技术保护点】
1.一种气动力多源数据融合和辨识方法,其特征在于,包括:/n步骤1:建立全局相关的气动力数学模型;/n步骤2:根据现有计算数据、风洞试验数据和飞行试验数据,计算这三种来源数据的标准差;/n步骤3:将现有计算数据、风洞试验数据和飞行试验数据作为气动力数学模型的数据样本,建立超定方程组;/n步骤4:采用加权最小二乘法求解超定方程组,权值由标准差计算获得,解出气动力数学模型中的系数,完成多源数据融合气动力数学模型的建立。/n

【技术特征摘要】
1.一种气动力多源数据融合和辨识方法,其特征在于,包括:
步骤1:建立全局相关的气动力数学模型;
步骤2:根据现有计算数据、风洞试验数据和飞行试验数据,计算这三种来源数据的标准差;
步骤3:将现有计算数据、风洞试验数据和飞行试验数据作为气动力数学模型的数据样本,建立超定方程组;
步骤4:采用加权最小二乘法求解超定方程组,权值由标准差计算获得,解出气动力数学模型中的系数,完成多源数据融合气动力数学模型的建立。


2.根据权利要求1所述的气动力多源数据融合和辨识方法,其特征在于,所述步骤1中的全局相关的气动力数学模型的表达式为:



其中,f(Ma,αΦ,Φ,δP,δY,δR)表示气动力系数或气动力矩系数,包括法向力系数、俯仰力矩系数、侧向力系数、偏航力矩系数、滚转力矩系数和轴向力系数,轴向力系数不包括摩阻和底阻部分;Ma为来流马赫数;αΦ为合成攻角,速度矢量与体轴的夹角;Φ为气流滚转角,速度矢量在飞行器横截面的分量与选定子午面的夹角;δP为俯仰通道舵偏角;δY为偏航通道舵偏角;δR为滚转通道舵偏角;φ为初始相位角;表示对i、j、k、l、r和s所有取值范围进行求和;aijklrs为常系数,i、j、k、l均为0,1,2,…;Tr(x)为马赫数相关的切比雪夫多项式,Ts(y)为合成攻角相关的切比雪夫多项式。


3.根据权利要求2所述的气动力多源数据融合和辨识方法,其特征在于,所述Tr(x)计算式为:
Tr(x)=cos(rcos-1(x)),r=1,2,…;





4.根据权利要求2所述的气动力多源数据融合和辨识方法,其特征在于,所述Ts(y)计算式为:
Ts(y)=cos(scos-1(y)),s=1,2,…;



其中下标max表示最大值,下标min表示最小值。


5.根据权利要求2所述的气动力多源数据融合和辨识方法,...

【专利技术属性】
技术研发人员:傅建明李欣益唐海敏伍彬梁伟李小林段旭吴潜施振兴施臣钢
申请(专利权)人:上海机电工程研究所
类型:发明
国别省市:上海;31

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1