首页 / 专利库 / 人工智能 / 半实物仿真 / 一种飞行器半实物仿真系统误差分析方法

一种飞行器半实物仿真系统误差分析方法

阅读:693发布:2020-05-17

专利汇可以提供一种飞行器半实物仿真系统误差分析方法专利检索,专利查询,专利分析的服务。并且本 发明 提供了一种 飞行器 半实物仿真 系统误差分析方法,以飞行器 飞行控制系统 数学方程、两轴转台传递函数、三轴转台传递函数为 基础 ,与飞行器与目标相对运动学方程和比例导引方程联立,基于协方差分析描述函数法建立均值和协方差方程,利用均值和协方差方程,求解飞行器半实物仿真系统的误差。解决了已有导弹等半实物仿真系统误差方法中没有包括飞行器的飞行控制系统数学方程的计算误差。本发明获得的误差包括飞行控制系统数学方程的计算误差,可以用于揭示仿真实验结果中的异常现象,或者修正仿真实验结果。,下面是一种飞行器半实物仿真系统误差分析方法专利的具体信息内容。

1.一种飞行器半实物仿真系统误差分析方法,其特征在于,步骤如下:
(a)忽略目标模拟器的误差;
 (b)建立仿真系统的状态方程
已有飞行器飞行控制数学方程含有状态变量,包括了目标状态变量为
飞行器状态变量为
飞行器与目标视线变量为 ,另外,
飞行器质量为 飞行器受到的推为 飞行器受到的阻力为 飞行器受到的升力为 飞行器受到的侧向力为 飞行器单位时间内质量消耗量为 ,其中:
分别表示目标在惯性坐标系 中的三维坐标, 表示目标在惯性坐
标系 中的运动速度, 表示目标在惯性坐标系 中的俯仰角, 表示目标在惯性坐标系 中的偏航角, 分别表示飞行器在惯性坐标系
中的三维坐标, 表示飞行器在惯性坐标系 中的运动速度, 表示飞行器在惯性坐标系 中的弹道倾角, 表示飞行器在惯性坐标系 中的弹道偏角, 表示飞行器在速度坐标系 中的攻角, 表示飞行器在速度坐标系 中的侧滑角, 表示
飞行器在速度坐标系 中的倾斜角, 分别飞行器沿着飞行器
动坐标系 三个坐标轴的旋转角速度, 表示飞行器在惯性坐标系 中的俯仰角, 表示飞行器在惯性坐标系 中的偏航角,表示飞行器在惯性坐标系 中的滚转角;
以高低视线角 、方位视线角 、飞行器旋转角速度 和飞行
姿态角 、 、 均为随机状态变量,其余状态变量为确定性变量,建立含有随机状态变量的状态方程为:
(3)
式中: 分别为飞行器在动坐标系 中的气动力矩;
分为飞行器在动坐标系 中的转动惯量; ,
, , , ;
在飞行器的飞行控制数学方程当中,其余含有确定性状态变量的状态方程是飞行控制力学的通用方程,见公式(3a);
(3a)
(c)利用描述函数将非线性的状态方程进行统计线性化,成为线性状态方程随机状态变量方程(3)中含有非线性函数,根据(3)式中的非线性函数的概率密度函数形式,求出一个拟线性表达式,用该拟线性表达式等效代替非线性函数,使得该式与非线性函数的均方差达到极小值;
得到单变量非线性函数 的统计线性化可用描述函数 表达为:
,式中: ; ; ;
; ;和 统称为描述函数;
得到双变量非线性函数 的统计线性化描述函数为:
,式中: ; ;
; ; ; ; ;

据此,随机状态变量方程(3)中的非线性函数可以拟合为线性表达式,例如:
用类似方法完成其余函数的统计线性化;
(d)利用协方差分析法,导出仿真系统的均值和协方差方程
将统计线性函数代入到(3)式中,整理得到均值方程(4)和协方差方程(5):
(4)
  (5)
(5)式中:


其中: ;




用类似方法可以求出其他各项;
在(4)式、(5)式中,协方差矩阵P的对角线元素分别是8个随机状态变量的均方差,非对角线元素是不同随机状态变量的协方差,该矩阵是关于主对角线对称的,N为拟线性系统动态矩阵;
(e)利用均值和协方差方程,求解得出飞行器半实物仿真系统的误差
将均值方程(4)和协方差方程(5)与飞行器动力学方程、运动学方程、目标运动学方程、飞行器与目标相对运动学方程和比例导引中确定性方程联立,并给定初始条件,就能求出各随机状态变量在每一时刻的均值和协方差,其中协方差矩阵P的对角线元素分别是8个随机状态变量的均方差,就是飞行器半实物仿真系统的误差。

说明书全文

一种飞行器半实物仿真系统误差分析方法

技术领域

[0001] 本发明涉及一种飞行器半实物仿真系统误差分析方法。

背景技术

[0002] 将飞行器联接到半实物仿真系统中,进行半实物仿真实验时,输出结果会存在误差,需要进行误差分析,用于揭示仿真实验结果中的异常现象,或者修正仿真实验结果。
[0003] 已有导弹等半实物仿真系统误差方法,文献[1]~[4]有记载。但是,已有方法中没有包括飞行器的飞行控制系统数学方程的计算误差,因此,需要发明一种飞行器半实物仿真系统误差分析方法,其中要包括飞行控制系统数学方程的计算误差,获得的误差用于揭示仿真实验结果中的异常现象,或者修正仿真实验结果。

发明内容

[0004] 为了解决已有技术存在的问题,本发明提供了一种飞行器半实物仿真系统误差分析方法。本发明以飞行器飞行控制系统数学方程、两轴转台传递函数、三轴转台传递函数为基础,与飞行器与目标相对运动学方程和比例导引方程联立,基于协方差分析描述函数法建立均值和协方差方程,利用均值和协方差方程,求解飞行器半实物仿真系统的误差。
[0005] 本发明提供了一种飞行器半实物仿真系统误差分析方法,首先,分析并给出飞行器半实物仿真系统产生的误差类型如下:(1)两轴转台的误差
两轴转台用于模拟目标相对于飞行器的视线速度,其误差包括动态误差和静态误差两类。动态误差是由于两个轴的传递函数不为1而产生的幅值误差 和相位误差,式中: 表示误差、 表示传递函数的幅值、 表示传递函数的相位、表示虚数的单位、 表示角频率;静态误差是指两个轴的安装角误差,即高低视线角误差和方位视线角误差。由于飞行器采用比例导引,导引信息是视线角速度,因此静态误差对仿真结果影响很小。
[0006] (2)目标模拟器的误差目标模拟器实时提供目标、干扰及背景的模拟信号,这些信号随着导引头瞄准方向、目标与飞行器之间的相对距离、大气传输条件而实时变化。目标模拟器误差对仿真结果的影响通常可以忽略。
[0007] (3)三轴转台的误差三轴转台根据仿真计算机发出的飞行器姿态控制指令,模拟飞行器的实际飞行姿态。
设 时刻转台的飞行器理论姿态为俯仰角 、偏航角 和倾斜角 ,则该时刻理论弹体坐标系 与惯性坐标系 的转换矩阵为 ,式中:表示转换矩阵。
[0008] 由于三轴转台存在垂直度误差、回转误差和位置误差、相邻两框回转轴线不共面而引起的相交度误差的影响,以及动态误差。在 时刻转台实际的飞行器姿态角为,此时实际弹体坐标系 与惯性系 的转换矩阵为 。飞行器实际姿态角与理论姿态角的差异就是弹体姿态角误差 :
式中: 为惯性系下的单位矢量。
[0009] 根据三轴转台最大运动范围,取转台垂直度误差、回转误差、位置误差以及相交度误差的最大值,并叠加各部分误差,弹体姿态角静态误差的最大值为 。而其数学模型为二阶系统,三轴转台的三个框的动态误差由二阶传递函数给出,动态幅值误差不超过10%,相位延迟不大于10°。
[0010] (4)负载模拟器的误差在负载模拟器的工作频带10Hz内,限定幅差和相差均为10%,计算出时间常数和阻尼系数,得到负载模拟器的矩传递函数为:
式中: 负载模拟器的力矩,为拉普拉斯自变量。
[0011] (5)仿真计算机的误差仿真计算机实时解算飞行器动力学方程、运动学方程、目标运动学方程和飞行器与目标相对运动学方程,解算运动学的误差为算法引入的截断误差和舍入误差,该误差很小,忽略不计;解算动力学方程的误差包括飞行器质量特性误差、作用力误差和力矩误差,因为很难精确建立的气动方程,所以仅能利用洞吹风试验数据,采用查表、插值法得到气动力和气动力矩;气动力和气动力矩的最大误差范围为±15%。
[0012] 本发明提供了一种飞行器半实物仿真系统误差分析方法的步骤如下:(a)忽略目标模拟器的误差
虽然有负载模拟器的误差,但在飞行器飞行控制系统数学模型中没有相应的加载力矩项,也无法将其与飞行控制系统数学方程联立,进行误差分析。本发明以飞行器飞行控制系统数学方程、两轴转台传递函数、三轴转台传递函数为基础,基于协方差分析描述函数法建立均值和协方差方程,计算出任意时刻的均值和协方差。
[0013] (b)建立仿真系统的状态方程已有飞行器飞行控制数学方程含有状态变量,文献[5]和[6]有记载。包括了目标状态变量为 、飞行器状态变量为

、 、飞行器与目标视线角变量为
。另外,飞行器质量为 飞行器受到的推力为 飞行器受到的阻力为
飞行器受到的升力为 飞行器受到的侧向力为 飞行器单位时间内质量消耗
量为 。其中: 分别表示目标在惯性坐标系 中的三维坐标 ,表
示目标在惯性坐标系 中的运动速度,表示目标在惯性坐标系 中的俯仰角,表示目标在惯性坐标系 中的偏航角, 分别表示飞行器在惯性
坐标系 中的三维坐标, 表示飞行器在惯性坐标系 中的运动速度, 表示
飞行器在惯性坐标系 中的弹道倾角, 表示飞行器在惯性坐标系 中的弹道偏角, 表示飞行器在速度坐标系 中的攻角,表示飞行器在速度坐标系
中的侧滑角, 表示飞行器在速度坐标系 中的倾斜角,
分别飞行器沿着飞行器动坐标系 三个坐标轴的旋转角速度, 表示飞行器在惯性坐标系 中的俯仰角, 表示飞行器在惯性坐标系 中的偏航角, 表示飞行器在惯性坐标系 中的滚转角。
[0014] 以高低视线角 、方位视线角 、飞行器旋转角速度和飞行器姿态角 、 、 均为随机状态变量,其余状态变量为确定性变量。建立含有随机状态变量的状态方程为:
(3)
式中: 分别为飞行器在动坐标系 中的气动力矩;
分为飞行器在动坐标系 中的转动惯量; ,
, , , 。
[0015] 在飞行器的飞行控制数学方程当中,其余含有确定性状态变量的状态方程是飞行控制力学的通用方程,文献[5][6]有记载,见公式(3a)。
[0016] (3a)(c)利用描述函数将非线性的状态方程进行统计线性化,成为线性状态方程随机状态变量方程(3)中含有非线性函数,需要运用描述函数理论对其统计线性化。
[0017] 根据(3)式中的非线性函数的概率密度函数形式,求出一个拟线性表达式,用该拟线性表达式等效代替非线性函数,使得该式与非线性函数的均方差达到极小值。
[0018] 得到单变量非线性函数 的统计线性化可用描述函数 表达为:,式中: ; ; ;
; 和 统称为描述函数。
[0019] 得到双变量非线性函数 的统计线性化描述函数为:,式中: ; ;
; ; ; ; ;

[0020] 据此,随机状态变量方程(3)中的非线性函数可以拟合为线性表达式,例如:用类似方法完成其余函数的统计线性化。
[0021] (d)利用协方差分析法,导出仿真系统的均值和协方差方程将统计线性函数代入到(3)式中,整理得到均值方程(4)和协方差方程(5):
(4)
  (5)
(5)式中:


[0022] 其中: ;;



用类似方法可以求出其他各项。
[0023] 在(4)式、(5)式中,协方差矩阵 的对角线元素分别是8个随机状态变量的均方差,非对角线元素是不同随机状态变量的协方差,该矩阵是关于主对角线对称的, 为拟线性系统动态矩阵。
[0024] (e)利用均值和协方差方程,求解得出飞行器半实物仿真系统的误差将均值方程(4)和协方差方程(5)与飞行器动力学方程、运动学方程、目标运动学方程、飞行器与目标相对运动学方程和比例导引中确定性方程联立,并给定初始条件,就能求出各随机状态变量在每一时刻的均值和协方差,其中协方差矩阵 的对角线元素分别是8个随机状态变量的均方差,就是飞行器半实物仿真系统的误差。
[0025] 有益效果本发明提供了一种飞行器半实物仿真系统误差分析方法,以飞行器飞行控制系统数学方程、两轴转台传递函数、三轴转台传递函数为基础,与飞行器与目标相对运动学方程和比例导引方程联立,基于协方差分析描述函数法建立均值和协方差方程,利用均值和协方差方程,求解飞行器半实物仿真系统的误差。解决了已有导弹等半实物仿真系统误差方法中没有包括飞行器的飞行控制系统数学方程的计算误差。本发明获得的误差包括飞行控制系统数学方程的计算误差,可以用于揭示仿真实验结果中的异常现象,或者修正仿真实验结果。

具体实施方式

[0026] 一种飞行器半实物仿真系统误差分析方法,步骤如下:(a)忽略目标模拟器的误差
虽然有负载模拟器的误差,但在飞行器飞行控制系统数学模型中没有相应的加载力矩项,也无法将其与飞行控制系统数学方程联立,进行误差分析。本发明以飞行器飞行控制系统数学方程、两轴转台传递函数、三轴转台传递函数为基础,基于协方差分析描述函数法建立均值和协方差方程,计算出任意时刻的均值和协方差。
[0027] (b)建立仿真系统的状态方程已有飞行器飞行控制数学方程含有状态变量,文献[5]和[6]有记载。包括了目标状态变量为 、飞行器状态变量为
、 、
、飞行器与目标视线角变量为 。另外,飞行器质量为 飞
行器受到的推力为 飞行器受到的阻力为 飞行器受到的升力为 飞行器受到的侧向力为 飞行器单位时间内质量消耗量为 。其中: 分别表示
目标在惯性坐标系 中的三维坐标, 表示目标在惯性坐标系 中的运动速度,表示目标在惯性坐标系 中的俯仰角, 表示目标在惯性坐标系 中的偏航角,分别表示飞行器在惯性坐标系 中的三维坐标, 表示飞行器
在惯性坐标系 中的运动速度, 表示飞行器在惯性坐标系 中的弹道倾角, 表示飞行器在惯性坐标系 中的弹道偏角, 表示飞行器在速度坐标系 中的攻角, 表示飞行器在速度坐标系 中的侧滑角, 表示飞行器在速度坐标系
中的倾斜角, 分别飞行器沿着飞行器动坐标系 三个
坐标轴的旋转角速度, 表示飞行器在惯性坐标系 中的俯仰角, 表示飞行器在惯性坐标系 中的偏航角,表示飞行器在惯性坐标系 中的滚转角。
[0028] 以高低视线角 、方位视线角 、飞行器旋转角速度 和飞行器姿态角 、 、 均为随机状态变量,其余状态变量为确定性变量。建立含有随机状态变量的状态方程为:
  (3)
式中: 分别为飞行器在动坐标系 中的气动力矩;
分为飞行器在动坐标系 中的转动惯量; ,
, , , 。
[0029] 在飞行器的飞行控制数学方程当中,其余含有确定性状态变量的状态方程是飞行控制力学的通用方程,文献[5][6]有记载,见公式(3a)。
[0030]   (3a)(c)利用描述函数将非线性的状态方程进行统计线性化,成为线性状态方程随机状态变量方程(3)中含有非线性函数,需要运用描述函数理论对其统计线性化。
[0031] 根据(3)式中的非线性函数的概率密度函数形式,求出一个拟线性表达式,用该拟线性表达式等效代替非线性函数,使得该式与非线性函数的均方差达到极小值。
[0032] 得到单变量非线性函数 的统计线性化可用描述函数 、表达为: ,式中: ; ; ; ;和 统称为描述函数。
[0033] 得到双变量非线性函数 的统计线性化描述函数为:,式中: ; ;
; ; ; ;;

[0034] 据此,随机状态变量方程(3)中的非线性函数可以拟合为线性表达式,例如:用类似方法完成其余函数的统计线性化。
[0035] (d)利用协方差分析法,导出仿真系统的均值和协方差方程将统计线性函数代入到(3)式中,整理得到均值方程(4)和协方差方程(5):
 (4)
  (5)
(5)式中:


[0036] 其中: ;;



用类似方法可以求出其他各项。
[0037] 在(4)式、(5)式中,协方差矩阵 的对角线元素分别是8个随机状态变量的均方差,非对角线元素是不同随机状态变量的协方差,该矩阵是关于主对角线对称的, 为拟线性系统动态矩阵。
[0038] (e)利用均值和协方差方程,求解得出飞行器半实物仿真系统的误差将均值方程(4)和协方差方程(5)与飞行器动力学方程、运动学方程、目标运动学方程、飞行器与目标相对运动学方程和比例导引中确定性方程联立,并给定初始条件,就能求出各随机状态变量在每一时刻的均值和协方差,其中协方差矩阵 的对角线元素分别是8个随机状态变量的均方差,就是飞行器半实物仿真系统的误差。
[0039] 参考文献[1]万士正 等.某型导弹半实物仿真系统误差建模及影响分析.指挥控制与仿真,
2014,36(1):117-120.
[2]肖卫国 等.雷达寻的制导半实物仿真误差研究.计算机仿真,2007,24(5):259-
263.
[3]丁长明 等.半实物仿真系统误差研究.航空兵器,1996,2:32-36.
[4]康凤举 等.下航行器半实物仿真精度统计的CADFT法应用研究.系统仿真学报,
1999,11(6):423-425.
[5]袁子怀 等.有控飞行力学与计算机仿真[M].北京:国防工业出版社,2001,4:139-
171.
[6]徐明友.弹箭飞行动力学[M].北京:国防工业出版社,2003,1:7-20.
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈