SylixOS平台下的火控实时解算与实现

第61卷第2期兵器装备工程学报2222年2月【装备理论与装备技术】dal： 2.1229/bqzhgcxX2222.2.725 SylixOS平台下的火控实时解算与实现
程文博、，屈艺2吴盘龙、，苏鼎4
(1.南京理工大学自动化学院，南京22094；2.中国人民解放军66483部队，北京20093)
摘要:为了研究SylmOS嵌入式实时操作系统平台下火控解算的实时性和可靠性，搭建了火控系统的总体框架;结合外弹道数学模型和命中问题模型设计了一种基于外弹道的火控算法，充分发挥了基于弹道微分方程解算通用性好、精度高的优点;基于Sy/OS操作系统，实现了外弹道解算和火控命中问题的仿真。

通过对匀速目标和匀加速目标模型的仿真数据分析，并与射表解算的结果比较，验证了火控算法的可行性，为国产实时操作系统平台下的火控系统研制提供了参考。

关键词：Sy/xOS实时操作系统;外弹道解算;命中问题;火控解算
本文引用格式:程文博，屈艺，吴盘龙，等.Syl/OS平台下的火控实时解算与实现[J].兵器装备工程学报,2020,41
(10)：29-34c
Citation format：CHENG Wendo,QU Yi,WU Panlong,et ai.Real-Time Solu/on and Imp/menta/on of Fire Controi Undci SylixOS P/Worm J J].Joarnai of Ordnance Equipment Engineering,2020,41(10)：29-30.
中图分类号:TJ33文献标识码：A文章编号：2092-2304(2420)10-0022-06
Real-Time Solution and Implementation
of Firc Contrvi Undec SylixOS Platform
CHENG Wendo1,QU Yi1,WU Panlong1,SU Ding4
(1.Schooi of Automation,Nanjing University of Science and Technoloy/,Nanjing210094,China;
2.The Noc66483w T uo/b PLA,Beijing140093,China)
Abstract:Ia ordvr to stuPy thv real-timv and reliaPi/tp of thv fire coatrot system undvr thv SylixOS embedded real-timv o/erating system platWrm,thv overali frameworI of thv fire coatrot system was estaP/s6er And ther,aa extemai ballistic fire was desiaded based oa thv extemai ballistic matUematicai model and thv hit problem model•Thv coatrot aleorithm maPvs fuli use of thv aPvontagvs of/oo/eeapn/e and high precisioa based oa thv ballistic differertiai epuatioa,fmn/ ,based oa thv SylixOS o/erating system,thv simulatioa of thv vxWrnn ballistic solutioa and ire coatroi hit was realized•Throaen thv analysis of thv simulatioa data of thv consWat velocity motioa and thv consWat acce/ratioa motioa model, and comparine with thv resa/s of thv s/o1taPiv solutioa,thv fpsidi/ty of thv fire caatroi aleorithm was 01X60,which providvs a refereace iu thv devv//mert of thv fire coatroi system undvr thv domestic1>1/ Umv o/erating system p/Worm.
Key words:SylixOS real-time o/erating system;exWrnn ballistic solutioa;hit problem;fire caatroi solnt hoa
收稿日期：2016-11-2;修回日期：2420-02-28
基金项目：国家自然科学基金项目（6273153）航空科学基金项目（2019ZC59002）
作者简介:程文博（294—）,男，硕士研究生，主要从事火控系统研究,E-m/C285400287@qq.am。

通讯作者:屈艺（282—）,男，副研究员，主要从事目标跟踪和信号处理研究,E-m/CInuhqy@。

30兵器装备工程学报http：//schn■qks.cquP /
软件国产化是我国软件产业领域的一个重要议题，其已经成为国家安全的重要组成部分,特别是军工领域对信息安全、信息自主化有着迫切的需求,而目前使用较多VxWo/s、Li/ux、F/eRTOS等操作系统都是非国产操作系统”
由翼辉信息公司研发的国产操作系统SylixOS是一款支持对称多处理器的大型开源实时操作系统，采用抢占式、多任务、硬实时的方式来设计,具有功能完善、性能卓越、可靠稳定的特点”SylixOS在设计思路上借鉴了众多实时操作系统的设计思想，使得SylixOS在功能和具体性能上达到或超过了众多实时操作系统的水平，已经在网络国防安全、工业自动化、航空航天等诸多领域有较多应用⑴，研究国产实时操作系统下的火控系统有着重要意义”
文献［2］利用射表解射击诸元，分析出弹丸飞行时间的计算精度对射击误差有较大影响”文献°］介绍了外弹道微分方程组的龙格库塔解法，分析了外弹道解算的实时性”文献°］对弹道方程的边值问题提出了迭代-修正法,能够满足精度和实时性要求”文献°］对高炮火控系统中的命中问题建模和仿真，通过MATLAB仿真具有良好的效果”文献°］设计了VxWuks平台下火控系统,设计并实现了各个火控子系统的功能模块及通信”
火控系统要求能够准确、快速的击中目标，传统的火控算法通过对射表数据多维插值或以连续解析函数逼近射表来解射击诸元⑺,实时性好,但精度低,而通过弹道微分方程解射击诸元需要多次迭代，精度高,但实时性差”国产Sylix­OS操作系统实时性强、稳定性好，本文基于SylixOS操作系统设计了基于外弹道微分方程实时解算射击诸元的近程防空火控系统”
1近程防空火控系统设计
分析火控系统的功能和需求，设计了近程防空火力控制系统，如图1所示，该系统包括信息模拟模块、火控解算模块、外弹道模拟模块”
图1近程防空火力控制系统框图
信息模拟模块产生火控解算所需要的目标位置、速度、加速度信息和风速、风向、气温等环境参数;火控解算模块根据目标当前状态，进行命中问题求解，计算命中点和射击诸元;外弹道模拟模块根据火控解算的射击诸元模拟完整弹丸轨迹,用来判断是否命中”各个模块通过UDP网络协议实现通信，通过GPIN接口读取同步信号”火控系统通过同步信号产生的中断执行数据发送任务，以保证系统的实时性
2基于外弹道的火控算法设计
基于外弹道的火控算法首先建立了合理的外弹道模型,在此基础上研究了命中点的求解过程和射击诸元求解的迭代算法”为了快速准确地解算射击诸元，选取了合适的迭代算法，并给出了迭代初值的选取方法”
22坐标系定义与转换
弹丸运动方程建立于弹道坐标系，用于弹丸在运动过程中的定位;而目标运动方程建立于东北天坐标系，火控解算时需要将弹道坐标系转换到东北天坐标系”
2.12弹道坐标系与东北天坐标系的定义
2）弹道坐标系©-XpYpZp”
q：火炮炮口为原点；
X:水平面内,以火炮发射方向为X p轴正向；
Y:垂直于地面、天顶方向为正；
乙:垂直于射面,右手准则确定正向”
4）东北天坐标系：02-X z Y2H z O
武器发射管回转中心；
X:水平面内正东方向；
Y:水平面内正北方向；
H：垂直于地面、天顶方向为正”
2.1.2坐标系的转换
弹道坐标系与东北天坐标系的转换见图4所示”
图4东北天坐标系与弹道坐标系转换
如图2所示，弹丸发射方位角（弹丸射向与Y轴的夹角）为0,弹道坐标系先绕Y轴顺时针转动n/2-0,再绕旋转后的X p轴顺时针转动n/2到东北天坐标系”因此坐标系转换矩阵为：
_sinO
cosf
-0
cosf0T「Xp
-sin/0Y
0l」-乙」
(2)
2.2外弹道模型
外弹道模型主要有质点弹道模型和刚体弹道模型，不同的弹道模型是在不同的假设条件下得到的⑻”弹道运动描述越复杂,就越接近真实弹道”现阶段刚体弹道模型难以满
程文博，等:Sy/xOS 平台下的火控实时解算与实现
31
足火控解算中实时性要求,而质点弹道模型在射程较近，飞
行过程中弹丸质量不变时,两者相差很小，能满足射表的精
度要求，因此采用质点弹道模型。

质点弹道模型将弹丸作质点处理,不考虑弹丸飞行姿态 对弹道的影响，弹丸只受重力、空气阻力的作用。

以弹丸飞
行时间/为自变量，考虑风速对弹道的影响,得到弹丸运动
微分方程组4］为：
1丫(0 )二 Y o
式(4)中：Y = 4%,5y,5 *z,5,s,5］T ； Y 2 为解算初始值。

式(4)
的解法采用定步长四阶龙格库塔法46 ］，迭代公式如下：
I k 1 =f (\Y n )
k 2 = /(t n + h/2, Y n + h/2 * k 8 )
v k 3 =(t n + //2,Y n + //2 *k )
(5)
k 6 =(t n +/, +/ * k 3)
1 Y n + = Y n + h/6 * (k 8 + 0k
2 + 0k
3 + k 4)
式(5)中,4为积分步长。

通过式(5)可以解算任意射角下和
任意弹丸飞行时间下的弹道坐标。

2・3命中问题求解
命中点即目标和弹头的相遇点，通过联立目标运动方程
和弹丸运动方程求解。

如图3所示，目标从初始位置运动T m
时间到命中点,弹丸从发射到运动到该点的飞行时间为T , 则必有T d = T 。

建立目标命中方程为：I no %= -pMsc xon (勺4m
a,
=-p ()c A j 佥--n ,
二-p(y)^c xon ( J ,(， -
d
4m
o o d a U y a
d
l a
式(2)中为弹丸质心坐标弭%、八,弹丸质心速度沿三 个坐标轴的分量;；(y )为空气密度函数;m 为弹丸质量;S 为
弹丸横截面积是空气阻力系数，根据马赫数在阻力系数
表中线性插值得到;a 为音速;少%、5为风速;g 为重力加速 度;，为弹丸相对速度：
J 二槡(J — d )2 + J ( — d )2
(3)
方程的初始条件为：：=2,4%二J * COS 62，4y 二，2 * si 呛，
J =0,4二y 二z 二0。

％为弹丸初速，o 为火炮高低角。

将式(2)和初值条件记为：
(4)
广m 二^ ++；(©)
<
Y m 二 Y 。

+ 九(T )
(2)
二 ％ +U.T m )
建立弹丸命中方程为：
T d 二只乙血,％)
(7)
联立式(5)、式(6)，得命中方程为：
入二 E +(()
Y m = Y +((m )
v £ 二 Z 2 +(T m ) (8)
T d =(X m ，Y m ，H m)
1 T d 二 T m
为方便书写,可将式(8)过程简写为:T d-f(T m ) =0
图3命中问题示意图
T m 与T 具有非线性函数关系⑸，常用迭代算法有简单
迭代法、快速迭代法、牛顿法、弦截法,本文采用收敛速度较 快弦截法。

弦截法迭代方程为：
t 2 +1 t 2 _ f ( T m ) / rrk t 2 _ 1、 / 9、m _ m 二/(T) -(2_)) m _ m 丿
i 丿
通过不断迭代使t 2值越来越接近命中方程的解。

弦截
法迭代过程如下：
1)选取迭代的两个初值：
{T=0
(10)
L T 2 =-(匕X o + 匕丫2 + 孚/2)/(尤 + K + ^)
当目标正在靠近火炮时,解位于T ~ T 1之间。

4)计算(T m )(m )的值,通过式(9)计算T 2，T 4。

3)根据T d -T m 的符号更新解区间。

如果为正,则解在
区间T 4~T m ,否则在0〜T ；区间。

4)判断丨T m - T |是否满足| T m - T < e ,为迭代精度
值。

如果满足，则此时的T 就是式(8)的解。

如果不满足,
则在新的解区间重复上述过程。

2・4射击诸元求解
图3中M 点为命中点,也即弹丸终点，又已知弹丸起点
为原点，则求解射击诸元是两点边值问题，解决方法是将边
值问题转化为初值问题。

求解方法有迭代-修正法、二分 法、基于落点诸元、粒子群优化算法等41「2］。

迭代-修正法
具有较快的收敛速度，本文采用这种方法。

射击诸元解算的
32
兵器装备工程学报
http ：//schg. qks. j q ui , c d u , cn/
迭代-修正过程如图4所示。

图4 迭代-修正法解射击诸元过程框图3软件设计与实现
火控软件需要实现各个模块的通信来完成火控解算,为
了保证通信的实时性，还创建了中断服务子程序。

翼辉公司
开发了 Re/eveAAE 集成开发环境，方便用户可以在WBtows
系统下开发Sy/xOS 应用程序。

3.1 Reelevo-IDE 集成开发环境
Re/eveADE 运行于宿主机上，通过网络或串口与目标 机交互。

Re/eveADE 的交叉编译器在宿主机上编译生成可
以在目标机上运行的Sy/xOS 镜像文件，交叉调试器实现了 宿主机与目标机之间的前后台调试,火控应用程序的编译调
试过程如图5所示。

迭代初值为:
图5 火控程序通过RealeveAAE 在目标机上的
调试过程框图
3.2 SylixOS 驱动
Sy/xOS 设备驱动负责控制和监视硬件,并提供通用的
接口函数，包括字符设备驱动、块设备驱动、网络设备驱动、 总线子系统等。

应用程序调用驱动接口的过程如图5所示。

二七/-（弘/槡巧+尤）(11)
口9
式（K ）中：纬为外弹道解算的火炮高低角;A ）为弹道坐标 转换的方位角。

当未来点横坐标X .位于离散的弹道诸元（F-3, F ）之
间,可通过线性插值精确计算弹道终点（X ” ,,,乙）和弹丸飞
行时间。

修正后的+ + 3时刻的射击诸元为：
=a + T 9 -T e
-0+3 二 0+ + 人9 -人
式（12）中:T 、、。

分别为弹道终点对应的高低角、方位角，计
算方式如式（11）所示。

弹丸飞行时间为：
{ a+i (12)
d
r(3)-厂3(i)
(13)
式（13）中：为弹道微分方程解算的积分步长;.为解算的
积分步数。

通过迭代-修正法可以解出射击允许范围内目标轨迹
上任意一点的射击诸元。

图6 应用程序调用Sy/xOS 驱动过程框图
火控应用程序通过调用Sy/xOS 驱动提供的API 实现对
UDP 网络、GPIA 接口、RS232串口等设备的使用。

3.2 UDP 网络编程
火控系统既需要接收目标状态信息，也需要发送射击诸 元。

Sy/xOS 系统提供了 socket 通信的API ,帮助用户快速实
现服务器、客户端系统搭建。

创建UDP 通信需要的API 有：
1） bint 函数。

bBd 函数把本地协议地址赋予一个套接
字，用于绑定IL 地址或端口。

程文博，等/ylixOS平台下的火控实时解算与实现33
4)/cvf/m函数”/cvf/m函数用于接收客户端数据,存储在指定的缓冲区中,并返回读取到的字节个数”
3)sen/te函数”sen/te函数用于服务器端发送数据，只需要指定要发送的数据缓冲区的地址即可，返回读取到的字节个数”
火控系统的UDP回射模型如图7所示”
表2匀速运动目标射击诸元对比
弹飞时间/s高低角/mil方位角/mil 射表值3■292900359.4601622135.449060 MATLAB3■294376359.5239632135.219234 SylixOS 3.237032359.6424222135.274723
弹道轨迹和目标轨迹如图9所示”
图7火控系统的UDP通信流程框图
火控程序创建socbet套接字，通过bmn绑定本地地址和端口，通过/cvf/m^sen/te循环接收/发送数据，结束通信则关闭$0€0/”
3.2中断服务任务
火控系统需保证通信的实时性,采用中断的方式保证数据的及时发送”在接收到同步信号时，触发中断服务，使与中断向量绑定的数据发送任务优先级最高。

SylixOS系统中，一个中断服务函数对应一个中断描述符LW_CLASS_INTACT,该描述符会被加入到中断向量表的表项中”中断连接函数AP)_Inter_ector_on/eot将中断向量号与中断服务程序连接，中断释放函数APN_Inter/ectorkis-lie释放中断向量对应的中断服务函数”火控系统中以中断方式处理射击诸元发送任务的过程如图8所示”
图8数据发送任务在中断服务中的处理过程框图
4仿真
在RealevaDNE集成开发环境编写程序,通过网络或串口将编译后的镜像文件发送到SylixOS目标机”仿真的硬件平台采用0-4976CPU,主频3.60GHz”以某36mm榴弹为仿真弹丸，其初速为896m/s,质量为6■339ky,风速为0,积分步长为6.22、,射击诸元解算迭代误差值为0.C m,弦截法迭代精度为0.02、”
目标做匀速直线运动时，初始位置和速度为(2000叫2000m,1000m,-200m/s,-200m/s:-100m/、)”分别在MATLAB平台和SylixOS平台利用外弹道解算射击诸元,并与射表解算的结果比较，结果如表2所示”
1000y
800-
£600­
400...
200_
1000
图9弹道轨迹和匀速目标轨迹三维视图
弹道轨迹
目标轨迹
'跑
不同平台下解算的命中点坐标如表4所示”
表4MATLAB与SylixOS平台下的命中点坐标
MATLAB命中点SylixOS命中点距离/m2672.287OO2672.299310
偏流/m672.287OO671.299584
高度/m672.287OO672.299310
从表2中可以计算出:MATLAB平台和SylixOS平台解算的弹丸飞行时间与射表值的误差均小于0.02、,高低角和方位角的绝对误差小于0.2mii,表明了火控算法的可行性”两种平台下的解算结果误差为6■C mii,表明了SylixOS平台的准确性”通过200次蒙特卡洛仿真,SylixOS平台下火控解算平均耗时为14ms,表明了火控系统的快速性”从表2中可以计算出:命中点坐标误差分别为0・K8m、6.I385m、0.K8m,说明SylixOS解算的射击诸元能精确命中目标”
目标做匀加速运动时，初始位置和速度不变，加速度为(-50m/s2,-56m/s2,-50m/s2),其他设置不变,分别在MATLAB平台和SylixOS平台利用外弹道解算射击诸元，并与射表解算的结果比较，结果如表3所示”
表3匀加速运动目标射击诸元
弹飞时间/s射高/mil射向/mil 射表值 2.762925317.6543271180.532127 MATLAB 2.762423317.6582241130.O k496 SylixOS 2.757713317.C160041180.265
749
34
兵器装备工程学报
http ：//schg. qks. cqui. c d u . cn/
弹道轨迹和目标轨迹如图14所示。

800、•600、-400、-200…
01 000
1 000 T
图14 弹道轨迹和匀加速目标轨迹三维视图
----弹道轨迹----自标轨迹♦命中点
不同平台下解算的命中点坐标如表4所示。

表4 MATLAB 与Sy/xOS 平台下的命中点坐标
MATLAB 命中点
Sy/xOS 命中点
距离/m
1 532.959 8371 533.454 778
偏流/m
532.959 837533.554 778高度/m
532.959 837
533.554 778
从表3和表4计算出Sy/xOS 平台下的弹丸飞行时间与 射表值的绝对误差为0, 414 s :高低角、方位角绝对误差小于
0.3 mil,与MATLAB 相比，命中点的位置误差分别0, 00%、
0.18%、0. 18%，说明Sy/xOS 系统下对于匀加速目标的火控 解算也有较高的准确性。

通过100次蒙特卡洛仿真,得到该
平台下火控解算平均耗时为10 ms,表明实时操作系统下的 火控解算实时性较强。

5结论
在国产实时操作系统Sy/xOS 上实现了一种基于质点外 弹道解算射击诸元的火控算法，并对匀速目标和匀加速目标
进行仿真。

通过对仿真数据的分析，表明了弹丸能够准确、
快速的击中目标。

在20 ms 同步信号周期下持续解算，都满 足精度和实时性要求。

实时操作系统Sy/xOS 为火控系统的设计和开发提供了
一个较为理想的平台，在该系统下编程、调试简单，开发、维
护方便,性能稳定,实时性好,具有良好的发展前景。

本文中
目标运动在理想环境下,下一步模拟杂波环境下的目标跟踪
和火控解算，使系统更接近实际的火控系统。

参考文献：
[]焦进星.Sy/xOS 的来龙去脉[].软件与集成电路,2013
(7) ：63 -69.
[] 孟灿，毛征，孟博，等.弹丸飞行时间与射击诸元关系分
析[J ],兵工自动化,2016,35(11):24-27.
[3] 李丹•四阶龙格库塔法在火控解算中的应用[J],微计算
机信息,2011,27(3):192-193.
[4] 周启煌•战车火控外弹道实时解算的研究[J].火力与指
挥控制,2001,26(4)：15-13.
[]李强•基于外弹道的高炮火控算法与仿真火力与指
挥控制,2014,39(5):131 -134.
[]高金艳.VxWom (平台下的火控系统设计与仿真[].哈
尔滨工业大学学报,2015,47(9) ：42 -45.
[]郝伟，侯朝桢，王小艺.一种火炮快速解命中问题的新方
法[J ].火力与指挥控制,2007,32(10) ：65 -67.
J] LI D,GUO L H, The application of digitot signet processing
(DSP ) fos the reel time solving of artillero firo control exte ­riOS Wajectoro [ C ]//3m ILEE IntemationoI Conference on SoUwaN Engineering & Service Science. ILEE ,2012.
[9] EBAID A, Analysis of projectile motion in view of fractioxct
calcnlus [ J ]. Applieh MathematiccI Motelling , 2011,35
(3) ： 1231 -1239.
[10 ] CORTELL R, Application of the fourth-3mes Ruuge-Kutto
methob fos the solution of high-3mer generot inidot value
problems - J ], Computers & Structures : 1293,49 ( 5 ) ： 897
-900.
[11]秦鹏飞•基于大口径火炮的实时弹道解算方法与研究
[J].火炮发射与控制学报,2015,36(1):23 -70.
[10] 胡星，马国梁，葛敬飞•多旋翼无人机火控指令解算方法
研究[],兵器装备工程学报,2018,39(09) ：33 -33.
[13]崔静，邓方，方浩.基于改进粒子群优化算法的弹道求解
方法[]•东南大学学报(自然科学版),2013,43(S1)：
015 -218.
科学编辑 迟圣威博士(西安应用光学研究所高级工程师)
责任编辑周江川。