CN107145081A - 一种空馈式低频寻的半实物仿真试验方法及系统 - Google Patents

一种空馈式低频寻的半实物仿真试验方法及系统 Download PDF

Info

Publication number
CN107145081A
CN107145081A CN201710500115.0A CN201710500115A CN107145081A CN 107145081 A CN107145081 A CN 107145081A CN 201710500115 A CN201710500115 A CN 201710500115A CN 107145081 A CN107145081 A CN 107145081A
Authority
CN
China
Prior art keywords
mrow
msub
mtr
mtd
missile
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710500115.0A
Other languages
English (en)
Other versions
CN107145081B (zh
Inventor
王超磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Simulation Center
Original Assignee
Beijing Simulation Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Simulation Center filed Critical Beijing Simulation Center
Priority to CN201710500115.0A priority Critical patent/CN107145081B/zh
Publication of CN107145081A publication Critical patent/CN107145081A/zh
Application granted granted Critical
Publication of CN107145081B publication Critical patent/CN107145081B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种空馈式低频寻的半实物仿真试验方法,包括:搭建空馈式低频寻的半实物仿真系统;基于惯性坐标系、弹体坐标系和暗室内实验室坐标系的变换关系,确定天线阵列位置控制指令;确定转台控制角;控制天线阵列和三轴转台;采集姿态运动信息和目标信息。本发明还公开一种空馈式低频寻的半实物仿真试验系统。本发明中,通过采用转台外置、导引头固定的方式彻底消除转台对导引头探测精度的影响,解决了空馈低频的技术难题,准确、逼真的完成了空馈式低频寻的半实物仿真试验。

Description

一种空馈式低频寻的半实物仿真试验方法及系统
技术领域
本发明涉及半实物仿真试验技术领域。更具体地,涉及一种空馈式低频寻的半实物仿真试验方法及系统。
背景技术
半实物仿真是一种将控制器与在计算机上实现的控制对象的仿真模型联接在一起进行试验的技术。在这种试验中,控制器的动态特性、静态特性和非线性因素等都能真实地反映出来,因此它是一种更接近实际的仿真试验技术。这种仿真技术可用于修改控制器设计(即在控制器尚未安装到真实系统中之前,通过半实物仿真来验证控制器的设计性能,若系统性能指标不满足设计要求,则可调整控制器的参数,或修改控制器的设计),同时也广泛用于产品的修改定型、产品改型和出厂检验等方面。
半实物仿真具有以下特点:
(1)只能是实时仿真,即仿真模型的时间标尺和自然时间标尺相同。
(2)需要解决控制器与仿真计算机之间的接口问题。例如,在进行飞行器控制系统的半实物仿真时,在仿真计算机上解算得出的飞机姿态角、飞行高度、飞行速度等飞行动力学参数会被飞行控制器的传感器所感受,因而必须有信号接口或变换装置。这些装置例如是三自由度飞行仿真转台、动压-静压仿真器、负载力仿真器等。
(3)半实物仿真的实验结果比数学仿真更接近实际。
目前,三轴转台能为导引头和惯性测量组合提供角运动环境,是半实物仿真试验系统中必不可少的设备。由于转台由金属材料构成,在射频仿真暗室中必然产生电磁反射,对安装于内框上的无线电导引头产生干扰,因此暗室中的转台需要进行屏蔽处理。
低频寻的制导仿真试验中,导引头频率的降低给射频仿真试验系统带来了显著的问题。由于仿真暗室中转台在低频段对导引头探测精度影响显著,目前常采用路馈方式解决低频段仿真问题。路馈式由于信号直接注入,不能验证导引头中射频天线从接收无线电波到提取出有效信号的这部分功能,缺乏高可信度的导引头数学模型,因此,路馈式仿真无法对导引头进行全面检验。空馈式因为完全模拟了导引头真实的工作环境,因此能够全面验证导引头的性能。
因此,需要提供一种空馈式低频寻的半实物仿真试验方法及系统。
发明内容
本发明目的在于提供一种空馈式低频寻的半实物仿真试验方法,解决传统路馈注入式仿真试验方法无法对低频导引头进行全面验证的半实物仿真试验问题。
为达到上述目的,本发明采用下述技术方案:
一种空馈式低频寻的半实物仿真试验方法,包括:
S1:搭建空馈式低频寻的半实物仿真系统;空馈式低频寻的半实物仿真系统包括弹上计算机、惯性测量装置、低频探测导引头、三轴转台、低频暗室、低频信号生成分系统、天线阵列及馈电分系统和仿真计算机,其中,低频探测导引头和天线阵列及馈电分系统设置于低频暗室内部;
S2:基于惯性坐标系、弹体坐标系和暗室内实验室坐标系的变换关系,确定天线阵列位置控制指令;
S3:确定转台控制角;
S4:控制天线阵列和三轴转台;
S5:采集姿态运动信息和目标信息。
优选地,弹上计算机用于采集低频探测导引头和惯性测量装置的输出信息,进行导航和制导控制;
惯性测量装置用于测量三轴转台复现的姿态运动信息;
低频探测导引头用于测量天线阵列输出的低频辐射,向弹上计算机提供目标角度信息;
三轴转台用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境;
低频暗室用于提供无回波的自由空间环境;
低频信号生成分系统用于模拟雷达照射信号;
天线阵列及馈电系统用于模拟弹目的视线角运动;
仿真计算机用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
优选地,在低频暗室内采用无电磁反射或低电磁反射的支撑结构对导引头进行支撑。
优选地,惯性坐标系,x轴在水平面上朝前指向发射方向,y轴在垂直面内朝上,z轴符合右手定则;
弹体坐标系,x轴沿导弹中轴从弹尾指向弹头,y轴垂直朝上,z轴符合右手定则;
暗室内实验室坐标系,原点固定在导引头安装支架出,x轴朝前指向阵列,y轴垂直朝上,z轴符合右手定则。
进一步优选地,步骤S2具体包括以下步骤:
S201:在暗室中设置用于模拟射频目标的三元组天线阵列,被导引头探测后形成两轴目标视线角;
S202:计算惯性坐标系Axyz下目标的归一化位置:
其中,(x,y,z)为惯性坐标系下目标的归一化位置,qα为惯性坐标系下的俯仰视线角,qβ为惯性坐标系下的方位视线角;
S203:计算弹体坐标系Ox1y1z1下目标的归一化位置:
其中,(x1,y1,z1)为弹体坐标系下目标的归一化位置,qα1为弹体坐标系下的俯仰视线角,qβ1为弹体坐标系下的方位视线角,为惯性坐标系到弹体坐标系的转换矩阵,n表示发射系,1表示弹体系,ψ为惯性坐标系到弹体坐标系弹体的方位角的变换角度,θ为惯性坐标系到弹体坐标系弹体的俯仰角的变换角度,γ为惯性坐标系到弹体坐标系弹体的滚转角的变换角度;
S204:计算暗室内实验室坐标系Oxsyszs下目标的归一化位置:
其中,(xs,ys,zs)为暗室内实验室坐标系下目标的归一化位置,qαs为暗室内实验室坐标系下的俯仰视线角,qβs为暗室内实验室坐标系下的方位视线角,为弹体坐标系到暗室内实验室坐标系的转换矩阵,为弹体坐标系到暗室内实验室坐标系的方位角的变换角度,为弹体坐标系到暗室内实验室坐标系的俯仰角的变换角度,为弹体坐标系到暗室内实验室坐标系的滚转角的变换角度;
S205:计算天线阵列的控制角度:
进一步优选地,步骤S3具体包括以下步骤:
S301:建立转台坐标系,其中x轴为内框旋转轴,y轴朝上,z轴与x、y轴符合右手定则;
S302:转台坐标系与本体坐标系之间的变换关系为:
其中,(ψTTT)为转台控制角,(ψ,θ,γ)为弹体相对于惯性坐标系的姿态角,(ψ***)为转台坐标系和惯性坐标系之间的偏差角,为惯性坐标系到弹体坐标系的变换矩阵,为转台坐标系到惯性坐标系的变换矩阵,为转台坐标系到本体坐标系的变换矩阵;
S303:求解变换矩阵计算得到转台方位、俯仰和滚转三轴的控制角(ψTTT)。
进一步优选地,步骤S4具体包括以下步骤:
S401:仿真计算机根据S2得到的天线阵列位置控制指令,控制信号源和阵列天线;
S402:仿真计算机根据S3得到的转台方位、俯仰、滚转三轴的控制角(ψTTT),控制三轴转台运动到指定的位置。
优选地,弹上计算机采集惯性测量装置测量得到的姿态运动信息和低频探测导引头测量得到的目标角度信息,进行导航和制导控制。
本发明的另一个目的在于提供一种空馈式低频寻的半实物仿真试验系统,其特征在于,包括:弹上计算机、惯性测量装置、低频探测导引头、三轴转台、低频暗室、低频信号生成分系统、天线阵列及馈电分系统和仿真计算机,其中,低频探测导引头和天线阵列及馈电分系统设置于低频暗室内部,
弹上计算机用于采集低频探测导引头和惯性测量装置的输出信息,进行导航和制导控制;
惯性测量装置用于测量三轴转台复现的姿态运动信息;
低频探测导引头用于测量天线阵列输出的低频辐射,向弹上计算机提供目标角度信息;
三轴转台用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境;
低频暗室用于提供无回波的自由空间环境;
低频信号生成分系统用于模拟雷达照射信号;
天线阵列及馈电系统用于模拟弹目的视线角运动;
仿真计算机用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
本发明的有益效果如下:
本发明的优点在于采用转台外置、导引头固定的方式彻底消除转台对导引头探测精度的影响,解决了空馈低频的技术难题,准确、逼真的完成了空馈式低频寻的半实物仿真试验。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出空馈式低频寻的半实物仿真试验方法步骤图。
图2示出空馈式低频寻的半实物仿真系统组成示意图。
图3示出不同坐标系变换关系示意图。
图中:1、弹上计算机;2、惯性测量装置;3、低频探测导引头;4、三轴转台;5、低频暗室;6、低频信号生成分系统;7、天线阵列及馈电分系统;8、仿真计算机。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
空馈式仿真方法已经在2-18GHz的导引头仿真试验中得到了充分的验证,由于和高频段相比,暗室中模拟弹体姿态的三轴转台在低频段对导引头影响显著,所以该方法之前没有能够应用在低频(200M~2GHz)导引头仿真试验。
本发明针对转台对导引头的潜在影响,提出采用转台外置的空馈式低频制导半实物仿真试验方法。转台外置是指将安装惯性测量组合的三轴飞行转台从暗室内转到暗室外,而将导引头固定在暗室内。由于采用导引头固定安装的方式,在仿真试验过程中暗室内复现的弹目相对运动关系不再处于惯性空间,而是以导引头(弹体)为基准的非惯性空间。此时,外置于暗室的三轴飞行转台连同内框上的惯性测量组合仍处在惯性空间中来模拟弹体的姿态运动。
如图1所示,一种空馈式低频寻的半实物仿真试验方法的具体步骤为:
第一步搭建空馈式低频寻的半实物仿真系统
如图2所示,空馈式低频寻的半实物仿真系统,包括:弹上计算机1,惯性测量装置2,低频探测导引头3,三轴转台4,低频暗室5,低频信号生成分系统6,天线阵列及馈电分系统7,仿真计算机8。其中,弹上计算机1、惯性测量装置2和低频探测导引头3通称弹上实物。
其中,各部分功能介绍如下:
弹上计算机1用于采集低频探测导引头3和惯性测量装置2的输出信息,进行导航和制导控制。
惯性测量装置2用于测量三轴转台4复现的姿态运动信息。
低频探测导引头3用于测量天线阵列输出的低频辐射,向弹上计算机1提供目标角度信息。
三轴转台4用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境。
低频暗室5用于提供一个无回波的自由空间环境,防止外部的电磁信号对试验的影响。
低频信号生成分系统6用于模拟雷达照射信号。
天线阵列及馈电系统7用于模拟弹目的视线角运动。
仿真计算机8用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
和现有半实物仿真系统相比,本系统最大的区别是在暗室中采用无或低电磁反射的支撑结构代替了三轴飞行转台来对导引头进行支撑,而将三轴飞行转台置于射频暗室外,彻底消除低频信号下三轴飞行转台对导引头探测精度的影响。
第二步确定天线阵列位置控制指令
暗室内复现的弹目相对运动关系是以导引头(弹体)为基准的非惯性空间。
(1)在暗室中设置用于模拟射频目标的三元组天线阵列,被导引头探测后形成两轴目标视线角;
(2)计算惯性坐标系Axyz下目标的归一化位置:
其中,(x,y,z)为惯性坐标系下目标的归一化位置,qα为惯性坐标系下的俯仰视线角,qβ为惯性坐标系下的方位视线角;
(3)计算弹体坐标系Ox1y1z1下目标的归一化位置:
其中,(x1,y1,z1)为弹体坐标系下目标的归一化位置,qα1为弹体坐标系下的俯仰视线角,qβ1为弹体坐标系下的方位视线角,为惯性坐标系到弹体坐标系的转换矩阵,n表示发射系,1表示弹体系,ψ为惯性坐标系到弹体坐标系弹体的方位角的变换角度,θ为惯性坐标系到弹体坐标系弹体的俯仰角的变换角度,γ为惯性坐标系到弹体坐标系弹体的滚转角的变换角度;
应注意的是,通常两个坐标系的转换关系如下:
设坐标系a按照顺序分别经过的ψ,θ,γ角度旋转变换到坐标系b对应的旋转矩阵为则有:
(4)计算暗室内实验室坐标系Oxsyszs下目标的归一化位置:
其中,(xs,ys,zs)为暗室内实验室坐标系下目标的归一化位置,qαs为暗室内实验室坐标系下的俯仰视线角,qβs为暗室内实验室坐标系下的方位视线角,为弹体坐标系到暗室内实验室坐标系的转换矩阵,为弹体坐标系到暗室内实验室坐标系的方位角的变换角度,为弹体坐标系到暗室内实验室坐标系的俯仰角的变换角度,为弹体坐标系到暗室内实验室坐标系的滚转角的变换角度;
(5)计算天线阵列的控制角度:
三个坐标系的相互关系是用于在实验室坐标系下去控制目标阵列来模拟发射坐标系和弹体坐标系下的导弹-目标视线角。
第三步确定转台控制角
如图3所示,外置于暗室的三轴飞行转台连同内框上的惯性测量组合仍处在惯性空间中来模拟弹体的姿态运动。(ψTTT)为转台控制角,(ψ,θ,γ)为本体相对于赤道惯性系的姿态角,(ψ***)为转台坐标系和惯性坐标系之间的偏差角。为惯性坐标系到本体坐标系的变换矩阵,值得注意的是并不是赤道惯性系到本体坐标系的变换矩阵。为转台坐标系到惯性坐标系的变换矩阵,为转台坐标系到本体坐标系的变换矩阵。转台坐标系与本体坐标系之间的变换关系为:
通过求解变换矩阵可计算得到转台方位、俯仰、滚转三轴的控制角(ψTTT)。
第四步控制天线阵列和三轴转台4
仿真计算机8根据第二步得到的天线阵列位置控制指令,控制信号源和阵列天线。根据第三步得到的转台方位、俯仰、滚转三轴的控制角(ψTTT),控制三轴转台4运动到指定的位置。
第五步采集姿态运动信息和目标信息
弹上计算机1采集惯性测量装置2测量得到的姿态运动信息和低频探测导引头3测量得到的目标角度信息,进行导航和制导控制。
至此,实现了一种空馈式低频寻的半实物仿真试验。
本发明还公开了一种应用上述方法的空馈式低频寻的半实物仿真试验系统,包括:弹上计算机1、惯性测量装置2、低频探测导引头3、三轴转台4、低频暗室5、低频信号生成分系统6、天线阵列及馈电分系统7和仿真计算机8,其中,低频探测导引头3和天线阵列及馈电分系统7设置于低频暗室内部,
弹上计算机1用于采集低频探测导引头和惯性测量装置的输出信息,进行导航和制导控制;惯性测量装置2用于测量三轴转台复现的姿态运动信息;低频探测导引头3用于测量天线阵列输出的低频辐射,向弹上计算机提供目标角度信息;三轴转台4用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境;低频暗室5用于提供无回波的自由空间环境;低频信号生成分系统6用于模拟雷达照射信号;天线阵列及馈电系统7用于模拟弹目的视线角运动;仿真计算机8用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

Claims (9)

1.一种空馈式低频寻的半实物仿真试验方法,其特征在于,包括:
S1:搭建空馈式低频寻的半实物仿真系统;
所述空馈式低频寻的半实物仿真系统包括弹上计算机、惯性测量装置、低频探测导引头、三轴转台、低频暗室、低频信号生成分系统、天线阵列及馈电分系统和仿真计算机,其中,所述低频探测导引头和所述天线阵列及馈电分系统设置于所述低频暗室内部;
S2:基于惯性坐标系、弹体坐标系和暗室内实验室坐标系的变换关系,确定天线阵列位置控制指令;
S3:确定转台控制角;
S4:控制天线阵列和三轴转台;
S5:采集姿态运动信息和目标信息。
2.根据权利要求1所述的半实物仿真试验方法,其特征在于,
所述弹上计算机用于采集低频探测导引头和惯性测量装置的输出信息,进行导航和制导控制;
所述惯性测量装置用于测量三轴转台复现的姿态运动信息;
所述低频探测导引头用于测量天线阵列输出的低频辐射,向弹上计算机提供目标角度信息;
所述三轴转台用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境;
所述低频暗室用于提供无回波的自由空间环境;
所述低频信号生成分系统用于模拟雷达照射信号;
所述天线阵列及馈电系统用于模拟弹目的视线角运动;
所述仿真计算机用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
3.根据权利要求1所述的半实物仿真试验方法,其特征在于,在所述低频暗室内采用无电磁反射或低电磁反射的支撑结构对导引头进行支撑。
4.根据权利要求1所述的半实物仿真试验方法,其特征在于,
所述惯性坐标系,x轴在水平面上朝前指向发射方向,y轴在垂直面内朝上,z轴符合右手定则;
所述弹体坐标系,x轴沿导弹中轴从弹尾指向弹头,y轴垂直朝上,z轴符合右手定则;
所述暗室内实验室坐标系,原点固定在导引头安装支架出,x轴朝前指向阵列,y轴垂直朝上,z轴符合右手定则。
5.根据权利要求4所述的半实物仿真试验方法,其特征在于,所述步骤S2具体包括以下步骤:
S201:在暗室中设置用于模拟射频目标的三元组天线阵列,被导引头探测后形成两轴目标视线角;
S202:计算惯性坐标系Axyz下目标的归一化位置:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> <mtr> <mtd> <mi>z</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> 1
其中,(x,y,z)为惯性坐标系下目标的归一化位置,qα为惯性坐标系下的俯仰视线角,qβ为惯性坐标系下的方位视线角;
S203:计算弹体坐标系Ox1y1z1下目标的归一化位置:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mn>1</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msubsup> <mi>C</mi> <mi>n</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>&amp;psi;</mi> <mo>,</mo> <mi>&amp;theta;</mi> <mo>,</mo> <mi>&amp;gamma;</mi> <mo>)</mo> </mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> <mtr> <mtd> <mi>z</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msubsup> <mi>C</mi> <mi>n</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>&amp;psi;</mi> <mo>,</mo> <mi>&amp;theta;</mi> <mo>,</mo> <mi>&amp;gamma;</mi> <mo>)</mo> </mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;alpha;</mi> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mi>&amp;beta;</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,(x1,y1,z1)为弹体坐标系下目标的归一化位置,qα1为弹体坐标系下的俯仰视线角,qβ1为弹体坐标系下的方位视线角,为惯性坐标系到弹体坐标系的转换矩阵,n表示发射系,1表示弹体系,ψ为惯性坐标系到弹体坐标系弹体的方位角的变换角度,θ为惯性坐标系到弹体坐标系弹体的俯仰角的变换角度,γ为惯性坐标系到弹体坐标系弹体的滚转角的变换角度;
S204:计算暗室内实验室坐标系Oxsyszs下目标的归一化位置:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>s</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>s</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>s</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mi>s</mi> </mrow> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mi>s</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mi>s</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mi>s</mi> </mrow> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mi>s</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msubsup> <mi>C</mi> <mi>n</mi> <mi>s</mi> </msubsup> <mrow> <mo>(</mo> <mover> <mi>&amp;gamma;</mi> <mo>^</mo> </mover> <mo>,</mo> <mover> <mi>&amp;psi;</mi> <mo>^</mo> </mover> <mo>,</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mn>1</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msubsup> <mi>C</mi> <mi>n</mi> <mi>s</mi> </msubsup> <mrow> <mo>(</mo> <mover> <mi>&amp;gamma;</mi> <mo>^</mo> </mover> <mo>,</mo> <mover> <mi>&amp;psi;</mi> <mo>^</mo> </mover> <mo>,</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>cos</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mn>1</mn> </mrow> </msub> <mi>sin</mi> <mi> </mi> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,(xs,ys,zs)为暗室内实验室坐标系下目标的归一化位置,qαs为暗室内实验室坐标系下的俯仰视线角,qβs为暗室内实验室坐标系下的方位视线角,为弹体坐标系到暗室内实验室坐标系的转换矩阵,为弹体坐标系到暗室内实验室坐标系的方位角的变换角度,为弹体坐标系到暗室内实验室坐标系的俯仰角的变换角度,为弹体坐标系到暗室内实验室坐标系的滚转角的变换角度;
S205:计算天线阵列的控制角度:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mrow> <mi>&amp;alpha;</mi> <mi>s</mi> </mrow> </msub> <mo>=</mo> <mi>arcsin</mi> <mi> </mi> <msub> <mi>y</mi> <mi>s</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>q</mi> <mrow> <mi>&amp;beta;</mi> <mi>s</mi> </mrow> </msub> <mo>=</mo> <mi>arctan</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>z</mi> <mi>s</mi> </msub> </mrow> <msub> <mi>x</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
6.根据权利要求5所述的半实物仿真试验方法,其特征在于,所述步骤S3具体包括以下步骤:
S301:建立转台坐标系,其中x轴为内框旋转轴,y轴朝上,z轴与x、y轴符合右手定则;
S302:转台坐标系与本体坐标系之间的变换关系为:
<mrow> <msubsup> <mi>A</mi> <mi>T</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;psi;</mi> <mi>T</mi> </msub> <mo>,</mo> <msub> <mi>&amp;theta;</mi> <mi>T</mi> </msub> <mo>,</mo> <msub> <mi>&amp;gamma;</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>A</mi> <mi>i</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>&amp;psi;</mi> <mo>,</mo> <mi>&amp;theta;</mi> <mo>,</mo> <mi>&amp;gamma;</mi> <mo>)</mo> </mrow> <msubsup> <mi>A</mi> <mi>T</mi> <mi>i</mi> </msubsup> <mrow> <mo>(</mo> <msup> <mi>&amp;psi;</mi> <mo>*</mo> </msup> <mo>,</mo> <msup> <mi>s</mi> <mo>*</mo> </msup> <mo>,</mo> <msup> <mi>&amp;gamma;</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> </mrow>
其中,(ψTTT)为转台控制角,(ψ,θ,γ)为弹体相对于惯性坐标系的姿态角,(ψ***)为转台坐标系和惯性坐标系之间的偏差角,为惯性坐标系到弹体坐标系的变换矩阵,为转台坐标系到惯性坐标系的变换矩阵,为转台坐标系到本体坐标系的变换矩阵;
S303:求解变换矩阵计算得到转台方位、俯仰和滚转三轴的控制角(ψTTT)。
7.根据权利要求6所述的半实物仿真试验方法,其特征在于,所述步骤S4具体包括以下步骤:
S401:仿真计算机根据S2得到的天线阵列位置控制指令,控制信号源和阵列天线;
S402:仿真计算机根据S3得到的转台方位、俯仰、滚转三轴的控制角(ψTTT),控制三轴转台运动到指定的位置。
8.根据权利要求1所述的半实物仿真试验方法,其特征在于,弹上计算机采集惯性测量装置测量得到的姿态运动信息和低频探测导引头测量得到的目标角度信息,进行导航和制导控制。
9.应用权利要求1-8中任一项所述方法的空馈式低频寻的半实物仿真试验系统,其特征在于,包括:弹上计算机、惯性测量装置、低频探测导引头、三轴转台、低频暗室、低频信号生成分系统、天线阵列及馈电分系统和仿真计算机,其中,所述低频探测导引头和所述天线阵列及馈电分系统设置于所述低频暗室内部,
所述弹上计算机用于采集低频探测导引头和惯性测量装置的输出信息,进行导航和制导控制;
所述惯性测量装置用于测量三轴转台复现的姿态运动信息;
所述低频探测导引头用于测量天线阵列输出的低频辐射,向弹上计算机提供目标角度信息;
所述三轴转台用于模拟弹体的姿态运动,为惯性测量组合提供角运动环境;
所述低频暗室用于提供无回波的自由空间环境;
所述低频信号生成分系统用于模拟雷达照射信号;
所述天线阵列及馈电系统用于模拟弹目的视线角运动;
所述仿真计算机用于得到导弹的飞行弹道和目标的运动轨迹,形成各试验设备的控制参数。
CN201710500115.0A 2017-06-27 2017-06-27 一种空馈式低频寻的半实物仿真试验方法及系统 Active CN107145081B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710500115.0A CN107145081B (zh) 2017-06-27 2017-06-27 一种空馈式低频寻的半实物仿真试验方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710500115.0A CN107145081B (zh) 2017-06-27 2017-06-27 一种空馈式低频寻的半实物仿真试验方法及系统

Publications (2)

Publication Number Publication Date
CN107145081A true CN107145081A (zh) 2017-09-08
CN107145081B CN107145081B (zh) 2020-08-14

Family

ID=59784493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710500115.0A Active CN107145081B (zh) 2017-06-27 2017-06-27 一种空馈式低频寻的半实物仿真试验方法及系统

Country Status (1)

Country Link
CN (1) CN107145081B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966164A (zh) * 2017-11-28 2018-04-27 北京仿真中心 一种基于五轴转台的天球幕星光标定方法
CN108563144A (zh) * 2018-03-30 2018-09-21 西安电子科技大学 一种弹载雷达信号处理半实物仿真测试系统
CN110260714A (zh) * 2019-05-21 2019-09-20 中国人民解放军海军工程大学 制导弹药外弹道半实物仿真平台及方法
CN111258229A (zh) * 2019-12-18 2020-06-09 中国船舶重工集团有限公司第七一0研究所 一种航空自导深弹控制系统半实物仿真系统
CN111522252A (zh) * 2020-04-02 2020-08-11 北京仿真中心 一种半实物仿真方法及系统
CN114184210A (zh) * 2021-12-03 2022-03-15 江西洪都航空工业集团有限责任公司 一种基于卧式转台的仿真试验方法、装置和系统
CN114371633A (zh) * 2021-12-03 2022-04-19 北京仿真中心 一种半实物仿真装置、系统及其方法
CN114396964A (zh) * 2021-12-10 2022-04-26 北京仿真中心 一种定位器基础板的安装方法

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2055603U (zh) * 1989-10-10 1990-04-04 罗洪发 超窄带甚低频式探测器寻的头
CN101122780A (zh) * 2007-09-04 2008-02-13 北京控制工程研究所 月球软着陆制导、导航与控制半物理仿真试验系统
CN102116641A (zh) * 2009-12-31 2011-07-06 北京控制工程研究所 一种深空自主导航星敏感器的半物理仿真试验系统及方法
CN102147987A (zh) * 2011-04-18 2011-08-10 中国人民解放军国防科学技术大学 飞行器导航、制导与控制技术教学实验装置
CN102410785A (zh) * 2011-12-13 2012-04-11 林德福 激光末制导炮弹半实物仿真用的转台
CN102589350A (zh) * 2012-01-09 2012-07-18 林德福 激光末制导炮弹研制用的半实物仿真系统
CN102590794A (zh) * 2012-02-28 2012-07-18 北京航空航天大学 一种宽带相参雷达目标模拟器
CN103560336A (zh) * 2013-10-16 2014-02-05 北京航天福道高技术股份有限公司 Ku频段双频双极化微带平面反射阵天线
CN103759729A (zh) * 2014-01-10 2014-04-30 北京空间飞行器总体设计部 采用捷联惯导的月球软着陆地面试验用初始姿态获取方法
CN103810018A (zh) * 2014-01-28 2014-05-21 北京仿真中心 一种组件化、参数化仿真模型的设计方法
CN204116836U (zh) * 2014-10-08 2015-01-21 南京长峰航天电子科技有限公司 一种导引头姿态模拟装置
CN104516348A (zh) * 2014-12-29 2015-04-15 贵州航天计量测试技术研究所 一种用于射频寻的半实物仿真系统的校准装置及校准方法
CN104536291A (zh) * 2014-11-25 2015-04-22 北京理工大学 基于射频系统模拟弹性振动对导引头测量信号影响的方法
CN104596357A (zh) * 2013-11-01 2015-05-06 上海机电工程研究所 射频/光学半实物仿真复合目标模拟系统
US9059682B2 (en) * 2008-07-14 2015-06-16 Macdonald, Dettwilwe And Associates Corporation Orthomode junction assembly with associated filters for use in an antenna feed system
CN105423823A (zh) * 2015-10-29 2016-03-23 中国科学院长春光学精密机械与物理研究所 激光制导导弹的注入式半实物仿真系统
CN105487409A (zh) * 2016-01-29 2016-04-13 中国航空无线电电子研究所 无人机空域综合飞行安全管控演示验证平台
CN106441374A (zh) * 2016-12-15 2017-02-22 北京仿真中心 一种实时半实物仿真全天区星光导航模拟装置
CN106681170A (zh) * 2016-11-22 2017-05-17 北京润科通用技术有限公司 一种半实物制导仿真方法及仿真系统
CN106767161A (zh) * 2016-12-09 2017-05-31 北京仿真中心 一种红外与可见光复合制导半实物仿真方法

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2055603U (zh) * 1989-10-10 1990-04-04 罗洪发 超窄带甚低频式探测器寻的头
CN101122780A (zh) * 2007-09-04 2008-02-13 北京控制工程研究所 月球软着陆制导、导航与控制半物理仿真试验系统
US9059682B2 (en) * 2008-07-14 2015-06-16 Macdonald, Dettwilwe And Associates Corporation Orthomode junction assembly with associated filters for use in an antenna feed system
CN102116641A (zh) * 2009-12-31 2011-07-06 北京控制工程研究所 一种深空自主导航星敏感器的半物理仿真试验系统及方法
CN102147987A (zh) * 2011-04-18 2011-08-10 中国人民解放军国防科学技术大学 飞行器导航、制导与控制技术教学实验装置
CN102410785A (zh) * 2011-12-13 2012-04-11 林德福 激光末制导炮弹半实物仿真用的转台
CN102589350A (zh) * 2012-01-09 2012-07-18 林德福 激光末制导炮弹研制用的半实物仿真系统
CN102590794A (zh) * 2012-02-28 2012-07-18 北京航空航天大学 一种宽带相参雷达目标模拟器
CN103560336A (zh) * 2013-10-16 2014-02-05 北京航天福道高技术股份有限公司 Ku频段双频双极化微带平面反射阵天线
CN104596357A (zh) * 2013-11-01 2015-05-06 上海机电工程研究所 射频/光学半实物仿真复合目标模拟系统
CN103759729A (zh) * 2014-01-10 2014-04-30 北京空间飞行器总体设计部 采用捷联惯导的月球软着陆地面试验用初始姿态获取方法
CN103810018A (zh) * 2014-01-28 2014-05-21 北京仿真中心 一种组件化、参数化仿真模型的设计方法
CN204116836U (zh) * 2014-10-08 2015-01-21 南京长峰航天电子科技有限公司 一种导引头姿态模拟装置
CN104536291A (zh) * 2014-11-25 2015-04-22 北京理工大学 基于射频系统模拟弹性振动对导引头测量信号影响的方法
CN104516348A (zh) * 2014-12-29 2015-04-15 贵州航天计量测试技术研究所 一种用于射频寻的半实物仿真系统的校准装置及校准方法
CN105423823A (zh) * 2015-10-29 2016-03-23 中国科学院长春光学精密机械与物理研究所 激光制导导弹的注入式半实物仿真系统
CN105487409A (zh) * 2016-01-29 2016-04-13 中国航空无线电电子研究所 无人机空域综合飞行安全管控演示验证平台
CN106681170A (zh) * 2016-11-22 2017-05-17 北京润科通用技术有限公司 一种半实物制导仿真方法及仿真系统
CN106767161A (zh) * 2016-12-09 2017-05-31 北京仿真中心 一种红外与可见光复合制导半实物仿真方法
CN106441374A (zh) * 2016-12-15 2017-02-22 北京仿真中心 一种实时半实物仿真全天区星光导航模拟装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钟显江: "平面空馈阵列天线的设计与研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966164A (zh) * 2017-11-28 2018-04-27 北京仿真中心 一种基于五轴转台的天球幕星光标定方法
CN107966164B (zh) * 2017-11-28 2020-06-02 北京仿真中心 一种基于五轴转台的天球幕星光标定方法
CN108563144A (zh) * 2018-03-30 2018-09-21 西安电子科技大学 一种弹载雷达信号处理半实物仿真测试系统
CN110260714A (zh) * 2019-05-21 2019-09-20 中国人民解放军海军工程大学 制导弹药外弹道半实物仿真平台及方法
CN111258229A (zh) * 2019-12-18 2020-06-09 中国船舶重工集团有限公司第七一0研究所 一种航空自导深弹控制系统半实物仿真系统
CN111522252A (zh) * 2020-04-02 2020-08-11 北京仿真中心 一种半实物仿真方法及系统
CN114184210A (zh) * 2021-12-03 2022-03-15 江西洪都航空工业集团有限责任公司 一种基于卧式转台的仿真试验方法、装置和系统
CN114371633A (zh) * 2021-12-03 2022-04-19 北京仿真中心 一种半实物仿真装置、系统及其方法
CN114371633B (zh) * 2021-12-03 2023-09-26 北京仿真中心 一种半实物仿真装置、系统及其方法
CN114396964A (zh) * 2021-12-10 2022-04-26 北京仿真中心 一种定位器基础板的安装方法
CN114396964B (zh) * 2021-12-10 2023-10-20 北京仿真中心 一种定位器基础板的安装方法

Also Published As

Publication number Publication date
CN107145081B (zh) 2020-08-14

Similar Documents

Publication Publication Date Title
CN107145081A (zh) 一种空馈式低频寻的半实物仿真试验方法及系统
CN106681170B (zh) 一种半实物制导仿真方法及仿真系统
CN106767548B (zh) 空间三坐标法检测火炮身管射击状态下指向性装置及方法
CN102393213B (zh) 天基探测与跟踪成像系统测试装置及测试方法
CN103245364B (zh) 一种星敏感器动态性能测试方法
CN103954179A (zh) 捷联红外导引头隔离度寄生回路评估系统
CN110262282A (zh) 火箭子级落区控制三轴转台半实物仿真系统及方法
CN104155538A (zh) 一种小卫星天线综合试验测试系统
CN108061477B (zh) 一种导引头与惯组系统间相对安装误差校正方法
CN105466476A (zh) 一种卫星光学特性的实验室测量方法
CN104296908A (zh) 三自由度气浮台干扰力矩组成测量装置
CN109445310A (zh) 一种红外成像制导飞行器半实物仿真系统及方法
CN103868648A (zh) 三轴气浮仿真实验平台的质心测量方法
CN112698580B (zh) 适用于红外精确制导炸弹的半实物仿真系统及仿真方法
CN105716891A (zh) 一种卫星隐身性能地面测试及在轨一体化验证方法
CN112379603B (zh) 一种射频制导仿真中捷联导引头安装偏心的补偿系统及方法
CN112668107B (zh) 基于立式转台的航天运载器试验方法、装置及存储介质
CN100443859C (zh) 用于光电经纬仪的目标仿真方法
CN105627982A (zh) 一种远程运载器倾斜瞄准方法
CN206496739U (zh) 空间三坐标法检测火炮身管射击状态下指向性装置
KR101742126B1 (ko) 고도 측정 제어 장치 및 그 제어 방법
CN113916499A (zh) 动平台光学测量设备跟踪性能检测系统和方法
CN202329647U (zh) 天基探测与跟踪成像系统的测试装置
CN112013716A (zh) 一种电视导引头模拟系统捕捉快速移动目标的方法
CN109932925B (zh) 多侦察点仿真方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant