CN105718727A - Stratospheric airship flight performance parameter estimation method and system - Google Patents

Stratospheric airship flight performance parameter estimation method and system Download PDF

Info

Publication number
CN105718727A
CN105718727A CN201610035229.8A CN201610035229A CN105718727A CN 105718727 A CN105718727 A CN 105718727A CN 201610035229 A CN201610035229 A CN 201610035229A CN 105718727 A CN105718727 A CN 105718727A
Authority
CN
China
Prior art keywords
airship
flight
circling
overbar
estimating
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
CN201610035229.8A
Other languages
Chinese (zh)
Other versions
CN105718727B (en
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.)
Academy of Opto Electronics of CAS
Original Assignee
Academy of Opto Electronics of CAS
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 Academy of Opto Electronics of CAS filed Critical Academy of Opto Electronics of CAS
Priority to CN201610035229.8A priority Critical patent/CN105718727B/en
Publication of CN105718727A publication Critical patent/CN105718727A/en
Application granted granted Critical
Publication of CN105718727B publication Critical patent/CN105718727B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及一种估计平流层飞艇飞行性能参数的方法和系统,包括:使飞艇进入稳定盘旋飞行状态;获取飞艇的地速、姿态信息和位置信息;根据飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。本发明提供的估计平流层飞艇飞行性能参数的方法和系统,可在高空稀薄大气环境下,利用飞艇内的导航系统给出的地速、位置、姿态角信息,计算出准确的水平风速、风向、空速、稳态盘旋半径、盘旋周期、盘旋角速率、侧滑角和攻角等飞行性能参数;且本发明过程简单、易于实现,只需固定主桨转速和偏航控制量,飞艇即可进入盘旋飞行状态,盘旋飞行使飞艇处于固定高度,风速和风向固定,使飞艇处于恒定风场下,减小了估计飞艇飞行性能参数的误差。

The invention relates to a method and system for estimating the flight performance parameters of a stratospheric airship, comprising: making the airship enter a stable hovering flight state; obtaining the ground speed, attitude information and position information of the airship; , to estimate the flight performance parameters of the airship. The method and system for estimating the flight performance parameters of a stratospheric airship provided by the present invention can calculate accurate horizontal wind speed and wind direction by using the ground speed, position, and attitude angle information given by the navigation system in the airship in a high-altitude thin atmosphere environment , airspeed, steady-state circling radius, circling period, circling angular rate, side-slip angle and angle of attack and other flight performance parameters; It can enter the circling flight state, and the circling flight keeps the airship at a fixed height, and the wind speed and wind direction are fixed, so that the airship is in a constant wind field, which reduces the error in estimating the flight performance parameters of the airship.

Description

一种估计平流层飞艇飞行性能参数的方法和系统A method and system for estimating flight performance parameters of a stratospheric airship

技术领域technical field

本发明涉及航空测控技术领域,特别涉及一种估计平流层飞艇飞行性能参数的方法和系统。The invention relates to the technical field of aviation measurement and control, in particular to a method and system for estimating flight performance parameters of a stratospheric airship.

背景技术Background technique

近年来在区域大气环境监测、防灾减灾、高分辨率实时监视、预警和导弹防御、反恐、区域通信等需求的驱动下,平流层飞艇引起了各国的普遍重视,美国、日本、俄罗斯、韩国、欧盟、中国等主要国家陆续启动了相关的研究计划,开始了较深入的研究开发工作。一些国家已经开始着手进行平流层飞艇样机的研制和飞行试验。In recent years, driven by the needs of regional atmospheric environment monitoring, disaster prevention and mitigation, high-resolution real-time monitoring, early warning and missile defense, anti-terrorism, and regional communications, stratospheric airships have attracted widespread attention from all countries. The United States, Japan, Russia, and South Korea Major countries such as China, the European Union, and China have successively launched relevant research programs and started more in-depth research and development work. Some countries have already started to develop and test the prototype of the stratospheric airship.

飞行试验是检验飞行性能指标是否满足要求、飞艇总体设计是否合理的重要依据。当前平流层飞艇尚处在探索阶段,开展一次试验飞行的周期较长,成本和代价也很大,因此往往希望能通过一次飞行试验评估出尽可能多的飞行特性参数。然而,现阶段从飞行试验中,获取平流层飞艇的飞行性能参数或飞行特性参数,存在几个根本性的困难:首先,相对风速或真风速无法测量。由于平流层环境的特殊性,目前还无法实时获取飞行时的相对风速信息(空速、攻角、侧滑角)或者真风速信息。这些信息对基础飞行特性参数的辨识至关重要,因为平流层飞艇的动力学特性与其相对于大气的移动速度相关,而非与相对地面的移动速度相关。其次,平流层飞艇处于风场环境。平流层风场并不恒定,风速/风向随高度的变化而变化,仅当高度变化很小时,才能近似看作恒定风场。由于平流层飞艇体积庞大,运动速度低,其动力学特性受风的影响很大。和低空飞艇可选择无风环境开展飞行试验不同,平流层飞艇的飞行试验无法隔离风的影响。在风速不可测的情况下,加剧了参数估计的困难。其三,平流层试验飞艇自身存在的一些不足,加剧了飞行工况保持的难度。例如,由于俯仰姿态调节困难,平流层飞艇保持定高飞行比较困难,一旦飞艇推力或空速发生变化,高度就会发生飘移,而高度飘移又往往导致艇体压力变化,产生充放气效应,使得飞艇的质量参数也发生改变。然而,保持特定的飞行工况,是通过试验方法辨识飞行特性参数的关键。The flight test is an important basis for checking whether the flight performance indicators meet the requirements and whether the overall design of the airship is reasonable. At present, the stratospheric airship is still in the exploratory stage, and the period of carrying out a test flight is long, and the cost and cost are also high. Therefore, it is often hoped that as many flight characteristic parameters as possible can be evaluated through a flight test. However, there are several fundamental difficulties in obtaining flight performance parameters or flight characteristic parameters of stratospheric airships from flight tests at this stage: First, relative wind speed or true wind speed cannot be measured. Due to the particularity of the stratospheric environment, it is currently impossible to obtain real-time relative wind speed information (airspeed, angle of attack, sideslip angle) or true wind speed information during flight. This information is crucial for the identification of basic flight characteristic parameters, since the dynamics of a stratospheric airship are related to its velocity relative to the atmosphere rather than relative to the ground. Second, the stratospheric airship is in a wind field environment. The stratospheric wind field is not constant, and the wind speed/wind direction changes with the change of height. Only when the height changes are small, it can be approximately regarded as a constant wind field. Due to the large size and low speed of the stratospheric airship, its dynamic characteristics are greatly affected by the wind. Unlike low-altitude airships, which can choose a windless environment to carry out flight tests, the flight tests of stratospheric airships cannot isolate the influence of wind. In the case of unpredictable wind speed, the difficulty of parameter estimation is exacerbated. Third, some deficiencies in the stratospheric test airship itself have exacerbated the difficulty of maintaining flight conditions. For example, due to the difficulty in adjusting the pitch attitude, it is difficult for a stratospheric airship to maintain a constant altitude flight. Once the thrust or airspeed of the airship changes, the altitude will drift, and the altitude drift will often lead to changes in the pressure of the hull, resulting in inflation and deflation effects. The quality parameters of the airship are also changed. However, maintaining specific flight conditions is the key to identifying flight characteristic parameters through experimental methods.

以风速估计为例,来说明通过试验方法评估飞行性能参数的困境。曾有人认为,可用下述方法提取风速:首先让飞艇处于无动力飘飞状态,此时导航系统,如GPS给出的地速即为风速。由于平流层风场具有较好的时间稳定性,在之后的飞行性能评估试验时,可认为风速保持不变。但实际上,平流层风场随高度分层变化,一旦动力开启,由于俯仰姿态调节的困难,几乎不可能使得飞艇保持在零攻角飞行状态,飞行高度会发生飘移,而高度变化一旦和充放气耦合,高度飘移很容易就会超过200m。此时,新高度上的风场已经不是之前飘飞时的风场。Wind speed estimation is used as an example to illustrate the dilemma of evaluating flight performance parameters by experimental methods. Someone once thought that the following method can be used to extract the wind speed: first let the airship be in a state of unpowered flight, and at this time, the ground speed given by the navigation system, such as GPS, is the wind speed. Since the stratospheric wind field has good time stability, the wind speed can be considered to remain unchanged in the subsequent flight performance evaluation tests. But in fact, the stratospheric wind field changes with the altitude layer. Once the power is turned on, it is almost impossible to keep the airship at zero angle of attack due to the difficulty of adjusting the pitch attitude, and the flight altitude will drift. With deflation coupling, the altitude drift can easily exceed 200m. At this time, the wind field at the new altitude is no longer the wind field when flying before.

数学上,平流层飞艇飞行特性参数的试验辨识问题,可归结为:在风场和空速未知的前提下,如何进行试验设计和辨识算法设计,使得仅利用相对于地面的移动轨迹、速度、姿态等其他信息,就可评估出所期望的飞行特性参数。目前在平流层及以上高度,还没有已见文献的方法能直接为平流层飞艇提供评价飞行性能所需要的实时水平风速风向、空速、稳态盘旋半径、盘旋周期、盘旋角速率、侧滑角和攻角信息。Mathematically, the experimental identification problem of the flight characteristic parameters of the stratospheric airship can be summarized as: how to design the experimental design and the identification algorithm under the premise that the wind field and airspeed are unknown, so that only the moving trajectory, speed, Attitude and other information can be used to evaluate the desired flight characteristic parameters. At present, in the stratosphere and above altitudes, there is no known method in the literature that can directly provide real-time horizontal wind speed and direction, airspeed, steady-state circling radius, circling period, circling angular rate, sideslip Angle and angle of attack information.

发明内容Contents of the invention

本发明所要解决的技术问题是如何准确估计平流层飞艇的飞行性能参数。The technical problem to be solved by the invention is how to accurately estimate the flight performance parameters of the stratospheric airship.

为此目的,本发明提出了一种估计平流层飞艇飞行性能参数的方法,包括:To this end, the present invention proposes a method for estimating the flight performance parameters of a stratospheric airship, comprising:

使所述飞艇进入稳定盘旋飞行状态;making the airship enter a stable circling flight state;

获取所述飞艇的地速、姿态信息和位置信息;Obtain the ground speed, attitude information and position information of the airship;

根据所述飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。Estimate the flight performance parameters of the airship according to the ground speed, attitude information and position information of the airship.

优选地,所述使所述飞艇进入稳定盘旋飞行状态包括:Preferably, said making the airship enter a stable circling flight state includes:

固定所述飞艇的主桨转速和偏航控制量,使所述飞艇进入盘旋飞行状态;当所述飞艇的飞行高度的波动值小于预设值时,即认为所述飞艇进入稳定盘旋飞行状态。Fixing the speed of the main propeller and the yaw control amount of the airship, so that the airship enters a circling flight state; when the fluctuation value of the flying height of the airship is less than a preset value, it is considered that the airship enters a stable circling flight state.

优选地,所述飞艇的飞行性能参数包括风速、空速、偏航操纵性能参数和纵、横向稳定性能参数。Preferably, the flight performance parameters of the airship include wind speed, airspeed, yaw control performance parameters and longitudinal and lateral stability performance parameters.

优选地,所述风速、空速的计算过程包括:Preferably, the calculating process of described wind speed, airspeed comprises:

根据所述风速、空速和地速的关系建立第一圆方程:Establish the first circle equation according to the relationship between the wind speed, air speed and ground speed:

(( VV xx -- WW xx )) 22 ++ (( VV ythe y -- WW ythe y )) 22 == VV aa 22

其中,(Vx,Vy)为艇载导航系统测出的地速,(Wx,Wy)为待估计的风速,Va为待估计的空速;Among them, (V x , V y ) is the ground speed measured by the boat navigation system, (W x , W y ) is the wind speed to be estimated, and Va is the air speed to be estimated ;

利用盘旋飞行时测出的多个地速(Vxi,Vyi),i=1,2,…,通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速VaUsing multiple ground speeds (V xi , V yi ) measured during circling flight, i=1, 2,..., through the circle fitting algorithm, calculate the wind speed to be estimated (W x , W y ) and the space to be estimated Velocity V a .

优选的,所述通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速Va的步骤包括:Preferably, the step of calculating the estimated wind speed (W x , W y ) and the estimated air speed V a through the circle fitting algorithm includes:

A、通过多个地速(Vxi,Vyi),i=1,2,…n,定义中间参数A. Define intermediate parameters through multiple ground speeds (V xi , V yi ), i=1, 2,...n

a=2Wx,b=2Wy,c=Va 2-(Wx 2+Wy 2)a=2W x , b=2W y , c=V a 2 −(W x 2 +W y 2 )

构造出最小二乘方程;Construct the least squares equation;

VV xx 11 VV ythe y 11 11 VV xx 22 VV ythe y 22 11 .. .. .. .. .. .. .. .. .. VV xx nno VV ythe y nno 11 aa bb cc == VV xx 11 22 ++ VV ythe y 11 22 VV xx 22 22 ++ VV ythe y 22 22 .. .. .. VV xx nno 22 ++ VV ythe y nno 22

B、求解最小二乘方程,计算出a,b,c;B. Solve the least squares equation and calculate a, b, c;

C、由a,b,c计算风速(Wx,Wy),空速VaC. Calculate wind speed (W x , W y ) and air speed V a from a, b, and c.

WW xx == aa 22 ,, WW ythe y == bb 22 ,, VV aa == cc ++ WW xx 22 ++ WW ythe y 22

优选地,所述偏航操纵性能参数包括:所述飞艇的盘旋半径、盘旋周期、盘旋角速度;所述偏航操纵性能参数的计算过程包括:Preferably, the yaw maneuvering performance parameters include: the circling radius, circling period, and circling angular velocity of the airship; the calculation process of the yaw maneuvering performance parameters includes:

根据所述飞艇的盘旋轨迹,建立第二圆方程According to the circling track of the airship, establish the second circle equation

[(x-Wxt)-x00]2+[(y-Wyt)-y00]2=R2 [(xW x t)-x 00 ] 2 +[(yW y t)-y 00 ] 2 =R 2

其中,(Wx,Wy)为风速,(x00,y00)为初始时刻圆心的位置,(x,y)为艇载导航系统测出的飞艇位置,R为待估计的飞艇的盘旋半径,t为时间;Among them, (W x , W y ) is the wind speed, (x 00 , y 00 ) is the position of the center of the circle at the initial moment, (x, y) is the position of the airship measured by the on-board navigation system, and R is the circling of the airship to be estimated Radius, t is time;

从所述飞艇的盘旋轨迹上选取一系列点(xi,yi,ti),i=1,2,…,通过圆拟合算法,计算待估计的飞艇的盘旋半径R。Select a series of points (x i , y i , t i ), i=1, 2, .

优选地,所述盘旋周期、盘旋角速度通过以下公式进行计算:Preferably, the spiral period and spiral angular velocity are calculated by the following formula:

ωω == VV aa RR

TT cc ii rr cc ll ee == 22 ππ ωω

其中,ω为盘旋角速率,Tcircle为盘旋周期。Among them, ω is the circle angular rate, and T circle is the circle period.

优选地,所述纵、横向稳定性能参数包括所述飞艇的平均侧滑角、瞬时侧滑角、平均攻角和瞬时攻角;Preferably, the longitudinal and lateral stability performance parameters include the average sideslip angle, instantaneous sideslip angle, average angle of attack and instantaneous angle of attack of the airship;

所述飞艇的平均侧滑角和瞬时侧滑角通过以下公式进行计算:The average sideslip angle and instantaneous sideslip angle of the airship are calculated by the following formula:

ββ ‾‾ == sinsin -- 11 (( VV ‾‾ aa ,, ythe y bb VV aa ))

ββ ii == sinsin -- 11 [[ (( VV aa ,, ythe y bb )) ii VV aa ]]

VV ‾‾ aa ,, xx bb VV ‾‾ aa ,, ythe y bb VV ‾‾ aa ,, zz bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, xx bb )) ii ΣΣ ii == 11 mm (( VV aa ,, ythe y bb )) ii ΣΣ ii == 11 mm (( VV aa ,, zz bb )) ii ii

VV aa ,, xx bb VV aa ,, ythe y bb VV aa ,, zz bb ii == AA (( φφ ii ,, θθ ii ,, ψψ ii )) VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii

VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii == VV xx VV ythe y VV zz ii -- WW xx WW ythe y 00

其中, V a , x V a , x V a , z i T 为空速在导航坐标系下的分量; V x V y V z i T 为地速在导航坐标系下的分量,[WxWy0]T为风速在导航坐标系下的分量, V a , x b V a , y b V a , z b i T 为空速在体轴系上的分量,A(φiii)为从导航坐标系转到体轴系的转换矩阵,m为盘旋轨迹上点的个数,i=1,2,…,m,为平均侧滑角,βi为ti时刻的侧滑角, V ‾ a , x b V ‾ a , y b V ‾ a , z b T 为空速在体轴系上的分量平均值,为空速在体轴系上的y轴方向分量平均值,(Va,yb)i为ti时刻空速在体轴系上的y轴方向分量。in, V a , x V a , x V a , z i T is the component of airspeed in the navigation coordinate system; V x V the y V z i T is the component of the ground speed in the navigation coordinate system, [W x W y 0] T is the component of the wind speed in the navigation coordinate system, V a , x b V a , the y b V a , z b i T is the component of airspeed on the body axis system, A(φ iii ) is the transformation matrix from the navigation coordinate system to the body axis system, m is the number of points on the circle trajectory, i=1, 2,...,m, is the average sideslip angle, β i is the sideslip angle at time t i , V ‾ a , x b V ‾ a , the y b V ‾ a , z b T is the average value of the airspeed components on the body axis, is the average value of the y-axis direction component of the space velocity on the body axis system, and (V a,yb ) i is the y-axis direction component of the space velocity on the body axis system at time t i .

优选地,所述飞艇的平均攻角和瞬时攻角通过以下公式进行计算:Preferably, the average angle of attack and the instantaneous angle of attack of the airship are calculated by the following formula:

α ‾ = tan - 1 ( V ‾ a , z b V ‾ a , x b ) α ‾ ≈ 1 m Σ i = 1 m θ i α ‾ = the tan - 1 ( V ‾ a , z b V ‾ a , x b ) or α ‾ ≈ 1 m Σ i = 1 m θ i

α i = tan - 1 ( V a , b z ) i ( V a , b x ) i α i ≈ θ i α i = the tan - 1 ( V a , b z ) i ( V a , b x ) i or α i ≈ θ i

其中,为平均攻角,为m个俯仰角的平均值,为空速在体轴系上的z轴方向分量平均值,为空速在体轴系上的x轴方向分量平均值,αi为ti时刻的瞬时攻角,θi为ti时刻的俯仰角,(Va,bz)i为ti时刻空速在体轴系上的z轴方向分量平均值,(Va,bx)i为ti时刻空速在体轴系上的x轴方向分量平均值。in, is the average angle of attack, is the average value of m pitch angles, is the average value of the z-axis direction component of airspeed on the body axis system, is the average value of airspeed components in the x-axis direction on the body axis, α i is the instantaneous angle of attack at time t i , θ i is the pitch angle at time t i , (V a,bz ) i is the airspeed at time t i The average value of the z-axis direction component on the body axis system, (V a,bx ) i is the average value of the x-axis direction component of airspeed on the body axis system at time t i .

另一方面,本发明还提供了一种估计平流层飞艇飞行性能参数的系统,包括:飞行单元、采集单元、估计单元;On the other hand, the present invention also provides a system for estimating flight performance parameters of a stratospheric airship, comprising: a flight unit, an acquisition unit, and an estimation unit;

所述飞行单元用于使所述飞艇进入稳定盘旋飞行状态;The flight unit is used to make the airship enter a stable circling flight state;

所述采集单元获取所述飞艇的地速、姿态信息和位置信息;The acquisition unit acquires ground speed, attitude information and position information of the airship;

所述估计单元用于根据所述飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。The estimation unit is used for estimating the flight performance parameters of the airship according to the ground speed, attitude information and position information of the airship.

本发明提供的估计平流层飞艇飞行性能参数的方法和系统,可在高空稀薄大气环境下,利用平流层飞艇内的导航系统给出的地速、位置、姿态角信息,计算出准确的水平风速、风向、空速、偏航操纵性能(稳态盘旋半径、盘旋周期、盘旋角速率)、纵横向稳定性(侧滑角和攻角)等飞行性能参数;且本发明过程简单、易于实现,只需固定主桨转速和偏航控制量,飞艇即可进入盘旋飞行状态,盘旋飞行使飞艇处于固定高度,风速和风向固定,使飞艇处于恒定风场下,减小了估计平流层飞艇飞行性能参数的误差。The method and system for estimating the flight performance parameters of a stratospheric airship provided by the present invention can calculate the accurate horizontal wind speed by using the ground speed, position, and attitude angle information given by the navigation system in the stratospheric airship in a high-altitude thin atmosphere environment , wind direction, airspeed, yaw control performance (steady-state circling radius, circling period, circling angular rate), longitudinal and lateral stability (sideslip angle and angle of attack) and other flight performance parameters; and the process of the present invention is simple and easy to implement, Only need to fix the rotation speed of the main propeller and the yaw control amount, the airship can enter the circling flight state, the circling flight makes the airship at a fixed altitude, the wind speed and wind direction are fixed, so that the airship is under a constant wind field, which reduces the estimated flight performance of the stratospheric airship parameter error.

附图说明Description of drawings

通过参考附图会更加清楚的理解本发明的特征和优点,附图是示意性的而不应理解为对本发明进行任何限制,在附图中:The features and advantages of the present invention will be more clearly understood by referring to the accompanying drawings, which are schematic and should not be construed as limiting the invention in any way. In the accompanying drawings:

图1示出了本发明一种估计平流层飞艇飞行性能参数的方法的流程示意图;Fig. 1 shows a schematic flow chart of a method for estimating the flight performance parameters of a stratospheric airship in the present invention;

图2示出了本发明一种估计平流层飞艇飞行性能参数的方法的一种实施方式的流程示意图;Fig. 2 shows a schematic flow chart of an embodiment of a method for estimating the flight performance parameters of a stratospheric airship in the present invention;

图3示出了在风场环境下地速、风速、空速矢量的关系示意图;Fig. 3 shows a schematic diagram of the relationship between ground speed, wind speed, and airspeed vectors in a wind field environment;

图4示出了风场环境下平流层飞艇的盘旋轨迹示意图;Fig. 4 shows the schematic diagram of the hovering trajectory of the stratospheric airship under the wind field environment;

图5为风场环境下盘旋轨迹摆线与圆的关系示意图;Fig. 5 is a schematic diagram of the relationship between the cycloid and the circle of the spiral trajectory in the wind field environment;

图6为侧滑角随时间的变化示意图;Figure 6 is a schematic diagram of the change of sideslip angle with time;

图7为攻角随时间的变化示意图。Figure 7 is a schematic diagram of the variation of the angle of attack with time.

具体实施方式detailed description

下面将结合附图对本发明的实施例进行详细描述。Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

如图1所示,本发明提供了一种估计平流层飞艇飞行性能参数的方法,As shown in Figure 1, the present invention provides a method for estimating the flight performance parameters of a stratospheric airship,

使飞艇进入稳定盘旋飞行状态;Make the airship enter a stable circling flight state;

获取飞艇的地速、姿态信息和位置信息;Obtain the ground speed, attitude information and position information of the airship;

根据飞艇的地速、姿态信息和位置信息,估计飞艇的飞行性能参数。According to the ground speed, attitude information and position information of the airship, the flight performance parameters of the airship are estimated.

如图2所示,其中较优的,使飞艇进入稳定盘旋飞行状态包括:固定飞艇的主桨转速和偏航控制量,使飞艇进入盘旋飞行状态;当飞艇的飞行高度的波动值小于预设值时,即认为进入稳定盘旋飞行状态。其中较优的,获取飞艇的地速、姿态信息和位置信息包括:通过飞艇的艇载导航系统获取所述飞艇的地速、姿态信息和位置信息。为了评价飞艇的飞行性能,需要获取的飞艇的飞行性能参数包括水平风速、风向、空速、偏航操纵性能参数和纵、横向稳定性能参数。其中,该方法适用于估计平流层飞艇的飞行性能参数。As shown in Figure 2, preferably, making the airship enter a stable circling flight state includes: fixing the main propeller speed and yaw control amount of the airship to make the airship enter a circling flight state; when the fluctuation value of the airship's flight altitude is less than the preset value, it is considered to enter a stable circling flight state. Preferably, acquiring the ground speed, attitude information and position information of the airship includes: acquiring the ground speed, attitude information and position information of the airship through an on-board navigation system of the airship. In order to evaluate the flight performance of the airship, the flight performance parameters of the airship that need to be obtained include horizontal wind speed, wind direction, airspeed, yaw control performance parameters and longitudinal and lateral stability performance parameters. Among others, the method is suitable for estimating flight performance parameters of stratospheric airships.

当飞艇进入稳态盘旋飞行时,空速大小保持不变,且会自动进入定高飞行状态,定高飞行时,可近似认为风场恒定,风速、风向保持不变。飞艇进入盘旋飞行后,其风速为(Wx,Wy),地速为(Vx,Vy),稳态盘旋飞行时可认为航速(空速)也是常数,其大小为Va,具体的,所述风速、空速的计算过程包括:如图3所示,根据所述风速、空速和地速的关系建立第一圆方程:When the airship enters steady-state circling flight, the airspeed remains unchanged, and it will automatically enter the state of constant-altitude flight. When the airship is in constant-altitude flight, it can be approximately considered that the wind field is constant, and the wind speed and wind direction remain unchanged. After the airship enters into a circling flight, its wind speed is (W x , W y ), and its ground speed is (V x , V y ). When the airship is in a steady circling flight, it can be considered that the airspeed (airspeed) is also a constant, and its magnitude is V a . The calculation process of described wind speed, air speed comprises: as shown in Figure 3, establishes the first circle equation according to the relation of described wind speed, air speed and ground speed:

(( VV xx -- WW xx )) 22 ++ (( VV ythe y -- WW ythe y )) 22 == VV aa 22 -- -- -- (( 11 ))

其中,(Vx,Vy)为艇载导航系统测出的地速,(Wx,Wy)为待估计的风速,Va为待估计的空速;Among them, (V x , V y ) is the ground speed measured by the boat navigation system, (W x , W y ) is the wind speed to be estimated, and Va is the air speed to be estimated ;

利用盘旋飞行时测出的多个地速(Vxi,Vyi),i=1,2,…,通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速VaUsing multiple ground speeds (V xi , V yi ) measured during circling flight, i=1, 2,..., through the circle fitting algorithm, calculate the wind speed to be estimated (W x , W y ) and the space to be estimated Velocity V a .

其中,圆拟合算法是通过圆上的轨迹点来估计圆心位置和半径的方法,一般圆的曲线方程为Among them, the circle fitting algorithm is a method of estimating the position and radius of the center of the circle through the trajectory points on the circle. The general curve equation of the circle is

(x-x0)2+(y-y0)2=R2(2)(xx 0 ) 2 +(yy 0 ) 2 =R 2 (2)

式中,x、y为平面上x方向和y方向的轨迹坐标,x0、y0为圆心的x坐标和y坐标,R为圆的半径。In the formula, x and y are the trajectory coordinates in the x and y directions on the plane, x 0 and y 0 are the x and y coordinates of the center of the circle, and R is the radius of the circle.

展开后为expands to

22 xx 00 xx ++ 22 ythe y 00 ythe y ++ RR 22 -- (( xx 00 22 ++ ythe y 00 22 )) == xx 22 ++ ythe y 22 -- -- -- (( 33 ))

make

aa == 22 xx 00 ,, bb == 22 ythe y 00 ,, cc == RR 22 -- (( xx 00 22 ++ ythe y 00 22 )) -- -- -- (( 44 ))

则有then there is

ax+by+c=x2+y2(5)ax+by+c=x 2 +y 2 (5)

已知轨迹点(xi,yi),i=1,2,…,n,由上式可构造出最小二乘方程Given the trajectory points ( xi ,y i ), i=1,2,…,n, the least squares equation can be constructed from the above formula

xx 11 ythe y 11 11 xx 22 ythe y 22 11 .. .. .. .. .. .. .. .. .. xx nno ythe y nno 11 aa bb cc == xx 11 22 ++ ythe y 11 22 xx 22 22 ++ ythe y 22 22 .. .. .. xx nno 22 ++ ythe y nno 22 -- -- -- (( 66 ))

估计出a,b,c后,即可得到圆曲线参数(圆心位置和半径)After estimating a, b, c, you can get the parameters of the circular curve (center position and radius)

xx 00 == aa 22 ,, ythe y 00 == bb 22 ,, RR == cc ++ xx 00 22 ++ ythe y 00 22 -- -- -- (( 77 ))

公式(1)所述的第一圆方程与公式(2)所述的一般圆的曲线方程,参数对应关系为,地速(Vx,Vy)对应轨迹坐标x、y,待估计的风速(Wx,Wy)对应圆心坐标x0、y0,待估计的空速Va对应圆半径R。则可按照上述步骤,通过地速(Vx,Vy)计算风速(Wx,Wy)和空速VaThe first circle equation described in the formula (1) and the curve equation of the general circle described in the formula (2), the parameter correspondence relationship is that the ground speed (V x , V y ) corresponds to the trajectory coordinates x, y, and the wind speed to be estimated (W x , W y ) corresponds to the circle center coordinates x 0 , y 0 , and the airspeed V a to be estimated corresponds to the circle radius R. Then the wind speed (W x , W y ) and the air speed V a can be calculated from the ground speed (V x , V y ) according to the above steps.

所以通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速Va的步骤包括:Therefore, through the circle fitting algorithm, the steps of calculating the estimated wind speed (W x , W y ) and the estimated air speed V a include:

A、通过多个地速(Vxi,Vyi),i=1,2,…n,定义中间参数A. Define intermediate parameters through multiple ground speeds (V xi , V yi ), i=1, 2,...n

aa == 22 WW xx ,, bb == 22 WW ythe y ,, cc == VV aa 22 -- (( WW xx 22 ++ WW ythe y 22 ))

构造出最小二乘方程;Construct the least squares equation;

VV xx 11 VV ythe y 11 11 VV xx 22 VV ythe y 22 11 .. .. .. .. .. .. .. .. .. VV xx nno VV ythe y nno 11 aa bb cc == VV xx 11 22 ++ VV ythe y 11 22 VV xx 22 22 ++ VV ythe y 22 22 .. .. .. VV xx nno 22 ++ VV ythe y nno 22

B、求解最小二乘方程,计算出a,b,c;B. Solve the least squares equation and calculate a, b, c;

C、由a,b,c计算风速(Wx,Wy),空速VaC. Calculate wind speed (W x , W y ) and air speed V a from a, b, and c.

WW xx == aa 22 ,, WW ythe y == bb 22 ,, VV aa == cc ++ WW xx 22 ++ WW ythe y 22

因此,通过测出的多个地速(Vxi,Vyi),i=1,2,…,以及第一圆方程可以估算出风速(Wx,Wy)和空速Va的值。此外,风速是一个矢量,根据风速的正负也可以确定风速的方向。Therefore, through the measured multiple ground speeds (V xi , V yi ), i=1, 2,..., and the first circle equation The values of wind speed (W x , W y ) and air speed V a can be estimated. In addition, the wind speed is a vector, and the direction of the wind speed can also be determined according to the positive or negative of the wind speed.

如图4和图5所示,在风场环境下,盘旋轨迹曲线为摆线,而非闭合圆曲线,摆线相当于轮子滚动时,轮缘所划出的轨线。若以轮心为参照,轮缘相对于轮心的轨迹仍然是个圆。其中较优的,偏航操纵性能参数包括:所述飞艇的盘旋半径、盘旋周期、盘旋角速度;所述偏航操纵性能参数的计算过程包括:As shown in Figure 4 and Figure 5, in the wind field environment, the spiral trajectory curve is a cycloid, not a closed circular curve, and the cycloid is equivalent to the trajectory drawn by the wheel rim when the wheel rolls. If the wheel center is used as a reference, the track of the rim relative to the wheel center is still a circle. Preferably, the yaw control performance parameters include: the circling radius, the circling period, and the circling angular velocity of the airship; the calculation process of the yaw control performance parameters includes:

根据所述飞艇的盘旋轨迹,建立第二圆方程According to the circling track of the airship, establish the second circle equation

[(x-Wxt)-x00]2+[(y-Wyt)-y00]2=R2(8)[(xW x t)-x 00 ] 2 +[(yW y t)-y 00 ] 2 =R 2 (8)

其中,(Wx,Wy)为风速,(x00,y00)为初始时刻圆心的位置,(x,y)为艇载导航系统测出的飞艇位置,R为待估计的飞艇的盘旋半径;Among them, (W x , W y ) is the wind speed, (x 00 , y 00 ) is the position of the center of the circle at the initial moment, (x, y) is the position of the airship measured by the on-board navigation system, and R is the circling of the airship to be estimated radius;

从所述飞艇的盘旋轨迹上选取一系列点(xi,yi,ti),i=1,2,…,通过圆拟合算法,计算待估计的飞艇的盘旋半径R。Select a series of points (x i , y i , t i ), i=1, 2, .

具体的,在估算出风速(Wx,Wy)后,假设初始时刻圆心的位置为(x00,y00),则时刻t,圆心的位置(x0,y0)为Specifically, after estimating the wind speed (W x , W y ), assuming that the position of the center of the circle at the initial moment is (x 00 , y 00 ), then at time t, the position of the center of the circle (x 0 , y 0 ) is

xx 00 == xx 0000 ++ WW xx tt ythe y 00 == ythe y 0000 ++ WW ythe y tt -- -- -- (( 99 ))

飞艇相对于圆心的位置也满足一般的圆曲线方程The position of the airship relative to the center of the circle also satisfies the general circular curve equation

(x-x0)2+(y-y0)2=R2(10)(xx 0 ) 2 +(yy 0 ) 2 =R 2 (10)

式中,x、y为平面上x方向和y方向的轨迹坐标,x0、y0为圆心的x坐标和y坐标,R为圆的半径。In the formula, x and y are the trajectory coordinates in the x and y directions on the plane, x 0 and y 0 are the x and y coordinates of the center of the circle, and R is the radius of the circle.

将公式(9)代入上式,即可得到第二圆方程Substituting formula (9) into the above formula, we can get the second circle equation

[(x-Wxt)-x00]2+[(y-Wyt)-y00]2=R2(11)[(xW x t)-x 00 ] 2 +[(yW y t)-y 00 ] 2 =R 2 (11)

make

Xx == xx -- WW xx tt YY == ythe y -- WW ythe y tt -- -- -- (( 1212 ))

得到如下的辨识方程get the following identification equation

[X-x00]2+[Y-y00]2=R2(13)[Xx 00 ] 2 +[Yy 00 ] 2 =R 2 (13)

取飞艇摆线轨迹上的某点为初始点,从轨迹上选取一系列点(xi,yi,ti),i=1,2,…,计算出(Xi,Yi),i=1,2,…,利用上述的圆拟合算法即可得到盘旋半径R的估计。Take a point on the cycloid trajectory of the airship as the initial point, select a series of points ( xi , y i , t i ), i=1, 2,... from the trajectory, and calculate (X i ,Y i ), i = 1, 2, ..., using the above-mentioned circle fitting algorithm, the estimation of the circle radius R can be obtained.

所述盘旋周期、盘旋角速度可以通过以下公式进行计算:The circling period and circling angular velocity can be calculated by the following formula:

ωω == VV aa RR -- -- -- (( 1414 ))

TT cc ii rr cc ll ee == 22 ππ ωω -- -- -- (( 1515 ))

其中,ω为盘旋角速率,Tcircle为盘旋周期。Among them, ω is the circle angular rate, and T circle is the circle period.

其中较优的,所述纵、横向稳定性能参数包括所述飞艇的平均侧滑角、瞬时侧滑角、平均攻角和瞬时攻角;Preferably, the longitudinal and lateral stability performance parameters include the average sideslip angle, instantaneous sideslip angle, average angle of attack and instantaneous angle of attack of the airship;

所述飞艇的平均侧滑角和瞬时侧滑角通过以下公式进行计算:The average sideslip angle and instantaneous sideslip angle of the airship are calculated by the following formula:

ββ ‾‾ == sinsin -- 11 (( VV ‾‾ aa ,, ythe y bb VV aa )) -- -- -- (( 1616 ))

ββ ii == sinsin -- 11 [[ (( VV aa ,, ythe y bb )) ii VV aa ]] -- -- -- (( 1717 ))

VV ‾‾ aa ,, xx bb VV ‾‾ aa ,, ythe y bb VV ‾‾ aa ,, zz bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, xx bb )) ii ΣΣ ii == 11 mm (( VV aa ,, ythe y bb )) ii ΣΣ ii == 11 mm (( VV aa ,, zz bb )) ii ii -- -- -- (( 1818 ))

VV aa ,, xx bb VV aa ,, ythe y bb VV aa ,, zz bb ii == AA (( φφ ii ,, θθ ii ,, ψψ ii )) VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii -- -- -- (( 1919 ))

VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii == VV xx VV ythe y VV zz ii -- WW xx WW ythe y 00 -- -- -- (( 2020 ))

其中, V a , x V a , y V a , z i T 为空速在导航坐标系下的分量; V x V y V z i T 为地速在导航坐标系下的分量,[WxWy0]T为风速在导航坐标系下的分量, V a , x b V a , y b V a , z b i T 为空速在体轴系上的分量,A(φiii)为从导航坐标系转到体轴系的转换矩阵,m为盘旋轨迹上点的个数,i=1,2,…,m,为平均侧滑角,βi为ti时刻的侧滑角, V ‾ a , x b V ‾ a , y b V ‾ a , z b T 为空速在体轴系上的分量平均值,为空速在体轴系上的y轴方向分量平均值,(Va,yb)i为ti时刻空速在体轴系上的y轴方向分量。in, V a , x V a , the y V a , z i T is the component of airspeed in the navigation coordinate system; V x V the y V z i T is the component of the ground speed in the navigation coordinate system, [W x W y 0] T is the component of the wind speed in the navigation coordinate system, V a , x b V a , the y b V a , z b i T is the component of airspeed on the body axis system, A(φ iii ) is the transformation matrix from the navigation coordinate system to the body axis system, m is the number of points on the circle trajectory, i=1, 2,...,m, is the average sideslip angle, β i is the sideslip angle at time t i , V ‾ a , x b V ‾ a , the y b V ‾ a , z b T is the average value of the airspeed components on the body axis, is the average value of the y-axis direction component of the space velocity on the body axis system, and (V a,yb ) i is the y-axis direction component of the space velocity on the body axis system at time t i .

具体的,考虑到垂直方向的风速很小,所以假设为零,由地速和风速可得出空速在导航系上分量(Va,x,Va,y,Va,z)Specifically, considering that the wind speed in the vertical direction is very small, it is assumed to be zero, and the airspeed component in the navigation system can be obtained from the ground speed and wind speed (V a,x ,V a,y ,V a,z )

VV aa ,, xx VV aa ,, ythe y VV aa ,, zz == VV xx VV ythe y VV zz -- WW xx WW ythe y 00 -- -- -- (( 21twenty one ))

设导航系到体轴系的姿态角为(ψ,θ,φ),将空速由导航系变换到体轴系 V a , x b V a , y b V a , z b = A n b V a , x V a , y V a , z - - - ( 22 ) Let the attitude angle from the navigation system to the body axis be (ψ, θ, φ), and transform the airspeed from the navigation system to the body axis V a , x b V a , the y b V a , z b = A no b V a , x V a , the y V a , z - - - ( twenty two )

转换矩阵根据姿态角(ψ,θ,φ)通过下式计算得到transformation matrix According to the attitude angle (ψ, θ, φ) calculated by the following formula

对盘旋轨迹上的每一个点,按照式(21)计算出(Va,yb)i,i=1,2,…,m。For each point on the spiral trajectory, (V a,yb ) i is calculated according to formula (21), i=1,2,...,m.

当滚动角φ较小时,(Va,yb)i的计算可简化为When the rolling angle φ is small, the calculation of (V a,yb ) i can be simplified as

(Va,yb)i≈[-Vxsinψ+Vycosψ+Vasinθsinφ]i(24)(V a,yb ) i ≈[-V x sinψ+V y cosψ+V a sinθsinφ] i (24)

其中,Vx,Vy分别为地速x和y方向的分量,Va为平均空速;Among them, V x , V y are the components of the ground speed in the x and y directions respectively, and V a is the average air speed;

Pick

VV ‾‾ aa ,, ythe y bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, ythe y bb )) ii -- -- -- (( 2525 ))

则有平均侧滑角为Then the average sideslip angle is

ββ ‾‾ == sinsin -- 11 (( VV ‾‾ aa ,, ythe y bb VV aa )) -- -- -- (( 2626 ))

公式(26)为稳态盘旋时的平均侧滑角,相应地在每个ti时刻,可计算出瞬时侧滑角βi Equation (26) is the average sideslip angle in steady-state circling, correspondingly at each time t i , the instantaneous sideslip angle β i can be calculated

ββ ii == sinsin -- 11 [[ (( VV aa ,, ythe y bb )) ii VV aa ]] -- -- -- (( 2727 ))

如图6所示,为飞艇瞬时侧滑角和平均侧滑角的曲线图,飞艇瞬时侧滑角βi相对于平均侧滑角的波动量,反映了飞艇横向稳定性,波动量越小,飞艇横向稳定性越好。As shown in Figure 6, it is a graph of the airship's instantaneous sideslip angle and the average sideslip angle, and the airship's instantaneous sideslip angle β i is relative to the average sideslip angle The fluctuation amount of reflects the lateral stability of the airship, the smaller the fluctuation amount, the better the lateral stability of the airship.

其中较优的,所述飞艇的平均攻角和瞬时攻角通过以下公式进行计算:Wherein preferably, the average angle of attack and the instantaneous angle of attack of the airship are calculated by the following formula:

α ‾ = tan - 1 ( V ‾ a , z b V ‾ a , x b ) α ‾ ≈ 1 m Σ i = 1 m θ i - - - ( 28 ) α ‾ = the tan - 1 ( V ‾ a , z b V ‾ a , x b ) or α ‾ ≈ 1 m Σ i = 1 m θ i - - - ( 28 )

α i = tan - 1 ( V a , b z ) i ( V a , b x ) i 或αi≈θi(29) α i = the tan - 1 ( V a , b z ) i ( V a , b x ) i or α i ≈ θ i (29)

其中为平均攻角,为m个俯仰角的平均值,为空速在体轴系上的z轴方向分量平均值,为空速在体轴系上的x轴方向分量平均值,αi为ti时刻的瞬时攻角,θi为ti时刻的俯仰角,(Va,bz)i为ti时刻空速在体轴系上的z轴方向分量平均值,(Va,bx)i为ti时刻空速在体轴系上的x轴方向分量平均值。in is the average angle of attack, is the average value of m pitch angles, is the average value of the z-axis direction component of airspeed on the body axis system, is the average value of airspeed components in the x-axis direction on the body axis, α i is the instantaneous angle of attack at time t i , θ i is the pitch angle at time t i , (V a,bz ) i is the airspeed at time t i The average value of the z-axis direction component on the body axis system, (V a,bx ) i is the average value of the x-axis direction component of airspeed on the body axis system at time t i .

具体的,平均攻角的估计有两种方法。其一,对平飞俯仰角θ取平均,则有平均攻角为Specifically, the average angle of attack There are two methods of estimating . First, take the average of the level flight pitch angle θ, then the average attack angle is

αα ‾‾ ≈≈ 11 mm ΣΣ ii == 11 mm θθ ii -- -- -- (( 3030 ))

其二,利用式(19)给出的体轴系空速分量(Va,xb,Va,yb,Va,zb)i计算Second, use the body shaft space velocity component (V a,xb ,V a,yb ,V a,zb ) i given by formula (19) to calculate

VV ‾‾ aa ,, xx bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, xx bb )) ii ,, VV ‾‾ aa ,, zz bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, zz bb )) ii -- -- -- (( 3131 ))

则有平均攻角为then the average angle of attack is

αα ‾‾ == tanthe tan -- 11 (( VV ‾‾ aa ,, zz bb VV ‾‾ aa ,, xx bb )) -- -- -- (( 3232 ))

相应地,在每个ti时刻,可计算出瞬时攻角αi Correspondingly, at each time t i , the instantaneous angle of attack α i can be calculated

α i = tan - 1 ( V a , b z ) i ( V a , b x ) i , 或αi≈θi(33) α i = the tan - 1 ( V a , b z ) i ( V a , b x ) i , or α i ≈ θ i (33)

如图7所示,为飞艇瞬时攻角和平均攻角的曲线图。瞬时攻角αi相对于平均攻角的波动量,反映了飞艇纵向稳定性。波动量越小,飞艇纵向稳定性越好。As shown in Figure 7, it is a graph of the instantaneous attack angle and the average attack angle of the airship. The instantaneous angle of attack α i relative to the mean angle of attack The amount of fluctuation reflects the longitudinal stability of the airship. The smaller the amount of fluctuation, the better the longitudinal stability of the airship.

另一方面,采用上述估计平流层飞艇飞行性能参数的方法,本发明还提供了一种估计平流层飞艇飞行性能参数的系统,包括:飞行单元、采集单元、估计单元;On the other hand, by adopting the method for estimating the flight performance parameters of a stratospheric airship, the present invention also provides a system for estimating the flight performance parameters of a stratospheric airship, comprising: a flight unit, an acquisition unit, and an estimation unit;

所述飞行单元用于使所述飞艇进入稳定盘旋飞行状态;The flight unit is used to make the airship enter a stable circling flight state;

所述采集单元获取所述飞艇的地速、姿态信息和位置信息;The acquisition unit acquires ground speed, attitude information and position information of the airship;

所述估计单元用于根据所述飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。The estimation unit is used for estimating the flight performance parameters of the airship according to the ground speed, attitude information and position information of the airship.

本发明提供的估计平流层飞艇飞行性能参数的方法和系统,可在高空稀薄大气环境下,利用飞艇内的导航系统给出的地速、位置、姿态角信息,计算出准确的水平风速、风向、空速、偏航操纵性能(稳态盘旋半径、盘旋周期、盘旋角速率)、纵横向稳定性(侧滑角和攻角)等飞行性能参数;且本发明过程简单、易于实现,只需固定主桨转速和偏航控制量,飞艇即可进入盘旋飞行状态,盘旋飞行使飞艇处于固定高度,风速和风向固定,使飞艇处于恒定风场下,减小了估计平流层飞艇飞行性能参数的误差。The method and system for estimating the flight performance parameters of a stratospheric airship provided by the present invention can calculate accurate horizontal wind speed and wind direction by using the ground speed, position, and attitude angle information given by the navigation system in the airship in a high-altitude thin atmosphere environment , airspeed, yaw control performance (steady-state circling radius, circling period, circling angular rate), longitudinal and lateral stability (sideslip angle and angle of attack) and other flight performance parameters; and the process of the present invention is simple and easy to implement, only need By fixing the speed of the main propeller and the yaw control amount, the airship can enter the state of circling flight. The circling flight keeps the airship at a fixed altitude, and the wind speed and direction are fixed, so that the airship is under a constant wind field, which reduces the time required for estimating the flight performance parameters of the stratospheric airship. error.

虽然结合附图描述了本发明的实施方式,但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型,这样的修改和变型均落入由所附权利要求所限定的范围之内。Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present invention. within the bounds of the requirements.

Claims (10)

1.一种估计平流层飞艇飞行性能参数的方法,其特征在于,1. A method for estimating the flight performance parameters of a stratospheric airship, characterized in that, 使所述飞艇进入稳定盘旋飞行状态;making the airship enter a stable circling flight state; 获取所述飞艇的地速、姿态信息和位置信息;Obtain the ground speed, attitude information and position information of the airship; 根据所述飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。Estimate the flight performance parameters of the airship according to the ground speed, attitude information and position information of the airship. 2.根据权利要求1所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述使所述飞艇进入稳定盘旋飞行状态包括:2. a kind of method for estimating the flight performance parameter of stratospheric airship according to claim 1, is characterized in that, described making described airship enter stable hover flight state comprises: 固定所述飞艇的主桨转速和偏航控制量,使所述飞艇进入盘旋飞行状态;当所述飞艇的飞行高度的波动值小于预设值时,即认为所述飞艇进入稳定盘旋飞行状态。Fixing the speed of the main propeller and the yaw control amount of the airship, so that the airship enters a circling flight state; when the fluctuation value of the flying height of the airship is less than a preset value, it is considered that the airship enters a stable circling flight state. 3.根据权利要求1-2任意一项所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,3. according to a kind of method for estimating the flight performance parameter of stratospheric airship according to any one of claim 1-2, it is characterized in that, 所述飞艇的飞行性能参数包括风速、空速、偏航操纵性能参数和纵、横向稳定性能参数。The flight performance parameters of the airship include wind speed, air speed, yaw control performance parameters and longitudinal and lateral stability performance parameters. 4.根据权利要求3所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述风速、空速的计算过程包括:4. a kind of method of estimating stratospheric airship flight performance parameter according to claim 3, is characterized in that, the calculation process of described wind speed, airspeed comprises: 根据所述风速、空速和地速的关系建立第一圆方程:Establish the first circle equation according to the relationship between the wind speed, air speed and ground speed: (( VV xx -- WW xx )) 22 ++ (( VV ythe y -- WW ythe y )) 22 == VV aa 22 其中,(Vx,Vy)为艇载导航系统测出的地速,(Wx,Wy)为待估计的风速,Va为待估计的空速;Among them, (V x , V y ) is the ground speed measured by the boat navigation system, (W x , W y ) is the wind speed to be estimated, and Va is the air speed to be estimated ; 利用盘旋飞行时测出的多个地速(Vxi,Vyi),i=1,2,…n,通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速VaUsing multiple ground speeds (V xi , V yi ) measured during circling flight, i=1, 2,…n, through the circle fitting algorithm, calculate the estimated wind speed (W x , W y ) and estimated Air velocity V a . 5.根据权利要求4所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述通过圆拟合算法,计算待估计的风速(Wx,Wy)和待估计的空速Va的步骤包括:5. a kind of method for estimating the flight performance parameter of stratospheric airship according to claim 4, is characterized in that, described by circle fitting algorithm, calculate wind speed to be estimated (W x , W y ) and space to be estimated The steps of speeding up V a include: A、通过多个地速(Vxi,Vyi),i=1,2,…n,定义中间参数A. Define intermediate parameters through multiple ground speeds (V xi , V yi ), i=1, 2,...n a=2Wx,b=2Wy,c=Va 2-(Wx 2+Wy 2)a=2W x , b=2W y , c=V a 2 −(W x 2 +W y 2 ) 构造出最小二乘方程;Construct the least squares equation; VV xx 11 VV ythe y 11 11 VV xx 22 VV ythe y 22 11 .. .. .. .. .. .. .. .. .. VV xx nno VV ythe y nno 11 aa bb cc == VV xx 11 22 ++ VV ythe y 11 22 VV xx 22 22 ++ VV ythe y 22 22 .. .. .. VV xx nno 22 ++ VV ythe y nno 22 B、求解最小二乘方程,计算出a,b,c;B. Solve the least squares equation and calculate a, b, c; C、由a,b,c计算风速(Wx,Wy),空速VaC. Calculate wind speed (W x , W y ) and air speed V a from a, b, and c. WW xx == aa 22 ,, WW ythe y == bb 22 ,, VV aa == cc ++ WW xx 22 ++ WW ythe y 22 6.根据权利要求3所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述偏航操纵性能参数包括:所述飞艇的盘旋半径、盘旋周期、盘旋角速度;所述偏航操纵性能参数的计算过程包括:6. a kind of method for estimating the flight performance parameter of stratospheric airship according to claim 3, is characterized in that, described yaw maneuvering performance parameter comprises: the circling radius of described airship, circling period, circling angular velocity; The calculation process of flight control performance parameters includes: 根据所述飞艇的盘旋轨迹,建立第二圆方程According to the circling track of the airship, establish the second circle equation [(x-Wxt)-x00]2+[(y-Wyt)-y00]2=R2 [(xW x t)-x 00 ] 2 +[(yW y t)-y 00 ] 2 =R 2 其中,(Wx,Wy)为风速,(x00,y00)为初始时刻圆心的位置,(x,y)为艇载导航系统测出的飞艇位置,R为待估计的飞艇的盘旋半径,t为时间;Among them, (W x , W y ) is the wind speed, (x 00 , y 00 ) is the position of the center of the circle at the initial moment, (x, y) is the position of the airship measured by the on-board navigation system, and R is the circling of the airship to be estimated Radius, t is time; 从所述飞艇的盘旋轨迹上选取一系列点(xi,yi,ti),i=1,2,…,通过圆拟合算法,计算待估计的飞艇的盘旋半径R。Select a series of points (x i , y i , t i ), i=1, 2, . 7.根据权利要求6所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述盘旋周期、盘旋角速度通过以下公式进行计算:7. a kind of method for estimating the flight performance parameter of stratospheric airship according to claim 6, is characterized in that, described circle period, circle angular velocity are calculated by following formula: ωω == VV aa RR TT cc ii rr cc ll ee == 22 ππ ωω 其中,ω为盘旋角速率,Tcircle为盘旋周期。Among them, ω is the circle angular rate, and T circle is the circle period. 8.根据权利要求3所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述纵、横向稳定性能参数包括所述飞艇的平均侧滑角、瞬时侧滑角、平均攻角和瞬时攻角;8. a kind of method for estimating stratospheric airship flight performance parameters according to claim 3 is characterized in that, described longitudinal and lateral stability performance parameters comprise the average sideslip angle of described airship, instantaneous sideslip angle, average attack angle and instantaneous angle of attack; 所述飞艇的平均侧滑角和瞬时侧滑角通过以下公式进行计算:The average sideslip angle and instantaneous sideslip angle of the airship are calculated by the following formula: ββ ‾‾ == sinsin -- 11 (( VV aa ‾‾ ,, ythe y bb VV aa )) ββ ii == sinsin -- 11 [[ (( VV aa ,, ythe y bb )) ii VV aa ]] VV ‾‾ aa ,, xx bb VV ‾‾ aa ,, ythe y bb VV ‾‾ aa ,, zz bb == 11 mm ΣΣ ii == 11 mm (( VV aa ,, xx bb )) ii ΣΣ ii == 11 mm (( VV aa ,, ythe y bb )) ii ΣΣ ii == 11 mm (( VV aa ,, zz bb )) ii ii VV aa ,, xx bb VV aa ,, ythe y bb VV aa ,, zz bb ii == AA (( φφ ii ,, θθ ii ,, ψψ ii )) VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii VV aa ,, xx VV aa ,, ythe y VV aa ,, zz ii == VV xx VV ythe y VV zz -- WW xx WW ythe y 00 其中, V a , x V a , y V a , z i T 为空速在导航坐标系下的分量; V x V y V z i T 为地速在导航坐标系下的分量,[WxWy0]T为风速在导航坐标系下的分量, V a , x b V a , y b V a , z b i T 为空速在体轴系上的分量,A(φiii)为从导航坐标系转到体轴系的转换矩阵,m为盘旋轨迹上点的个数,i=1,2,…,m,为平均侧滑角,βi为ti时刻的侧滑角, V ‾ a , x b V ‾ a , y b V ‾ a , z b T 为空速在体轴系上的分量平均值,为空速在体轴系上的y轴方向分量平均值,(Va,yb)i为ti时刻空速在体轴系上的y轴方向分量。in, V a , x V a , the y V a , z i T is the component of airspeed in the navigation coordinate system; V x V the y V z i T is the component of the ground speed in the navigation coordinate system, [W x W y 0] T is the component of the wind speed in the navigation coordinate system, V a , x b V a , the y b V a , z b i T is the component of airspeed on the body axis system, A(φ iii ) is the transformation matrix from the navigation coordinate system to the body axis system, m is the number of points on the circle trajectory, i=1, 2,...,m, is the average sideslip angle, β i is the sideslip angle at time t i , V ‾ a , x b V ‾ a , the y b V ‾ a , z b T is the average value of the airspeed components on the body axis, is the average value of the y-axis direction component of the space velocity on the body axis system, and (V a,yb ) i is the y-axis direction component of the space velocity on the body axis system at time t i . 9.根据权利要求8所述的一种估计平流层飞艇飞行性能参数的方法,其特征在于,所述飞艇的平均攻角和瞬时攻角通过以下公式进行计算:9. a kind of method of estimating stratospheric airship flight performance parameter according to claim 8, is characterized in that, the mean angle of attack of described airship and instantaneous angle of attack are calculated by following formula: α ‾ = tan - 1 ( V ‾ a , z b V ‾ a , x b ) α ‾ ≈ 1 m Σ i = 1 m θ i α ‾ = the tan - 1 ( V ‾ a , z b V ‾ a , x b ) or α ‾ ≈ 1 m Σ i = 1 m θ i α i = tan - 1 ( V a , b z ) i ( V a , b x ) i 或αi≈θi α i = the tan - 1 ( V a , b z ) i ( V a , b x ) i or α i ≈ θ i 其中,为平均攻角,为m个俯仰角的平均值,为空速在体轴系上的z轴方向分量平均值,为空速在体轴系上的x轴方向分量平均值,αi为ti时刻的瞬时攻角,θi为ti时刻的俯仰角,(Va,bz)i为ti时刻空速在体轴系上的z轴方向分量平均值,(Va,bx)i为ti时刻空速在体轴系上的x轴方向分量平均值。in, is the average angle of attack, is the average value of m pitch angles, is the average value of the z-axis direction component of airspeed on the body axis system, is the average value of airspeed components in the x-axis direction on the body axis, α i is the instantaneous angle of attack at time t i , θ i is the pitch angle at time t i , (V a,bz ) i is the airspeed at time t i The average value of the z-axis direction component on the body axis system, (V a,bx ) i is the average value of the x-axis direction component of airspeed on the body axis system at time t i . 10.一种估计平流层飞艇飞行性能参数的系统,其特征在于,包括:飞行单元、采集单元、估计单元;10. A system for estimating the flight performance parameters of a stratospheric airship, comprising: a flight unit, an acquisition unit, and an estimation unit; 所述飞行单元用于使所述飞艇进入稳定盘旋飞行状态;The flight unit is used to make the airship enter a stable circling flight state; 所述采集单元获取所述飞艇的地速、姿态信息和位置信息;The acquisition unit acquires ground speed, attitude information and position information of the airship; 所述估计单元用于根据所述飞艇的地速、姿态信息和位置信息,估计所述飞艇的飞行性能参数。The estimation unit is used for estimating the flight performance parameters of the airship according to the ground speed, attitude information and position information of the airship.
CN201610035229.8A 2016-01-19 2016-01-19 A kind of method and system of estimation stratospheric airship flight performance parameter Active CN105718727B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610035229.8A CN105718727B (en) 2016-01-19 2016-01-19 A kind of method and system of estimation stratospheric airship flight performance parameter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610035229.8A CN105718727B (en) 2016-01-19 2016-01-19 A kind of method and system of estimation stratospheric airship flight performance parameter

Publications (2)

Publication Number Publication Date
CN105718727A true CN105718727A (en) 2016-06-29
CN105718727B CN105718727B (en) 2018-09-28

Family

ID=56147246

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610035229.8A Active CN105718727B (en) 2016-01-19 2016-01-19 A kind of method and system of estimation stratospheric airship flight performance parameter

Country Status (1)

Country Link
CN (1) CN105718727B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106227060A (en) * 2016-07-19 2016-12-14 中国航天员科研训练中心 A kind of stratospheric airship load cabin electronic equipment Thermal characteristic simulation and simulation and prediction System and method for
CN106250581A (en) * 2016-07-13 2016-12-21 南京航空航天大学 A kind of stratosphere aerostatics multi-objective method based on NSGA II
CN106681337A (en) * 2017-01-05 2017-05-17 烟台南山学院 Height keeping flight control method for stratospheric airship based on odd-order sliding mode
CN106781706A (en) * 2016-12-16 2017-05-31 北京航空航天大学 Air traffic Track Design method based on wind field distribution
CN108195558A (en) * 2017-12-06 2018-06-22 太原航空仪表有限公司 A kind of ascending air measuring method for weather modification aircraft
CN108388263A (en) * 2018-03-01 2018-08-10 中国航空无线电电子研究所 Circular course autoflight system
CN109917493A (en) * 2019-04-03 2019-06-21 西安邮电大学 Device and method for measuring high-altitude wind energy resources
CN110548288A (en) * 2019-09-05 2019-12-10 腾讯科技(深圳)有限公司 Virtual object hit prompting method and device, terminal and storage medium
CN113515901A (en) * 2021-04-08 2021-10-19 中国科学院空天信息创新研究院 Stability evaluation method and device for stratospheric airship, electronic equipment and medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565451A (en) * 2011-12-28 2012-07-11 中国航空工业集团公司第六三一研究所 Method for measuring and calculating navigation wind direction and wind speed of general-purpose airplane
CN103852081A (en) * 2014-03-20 2014-06-11 南京航空航天大学 Vacuum speed resolving method for air data/serial inertial navigation combined navigation system
CN103994748A (en) * 2014-05-27 2014-08-20 中国航天空气动力技术研究院 Method for estimating trim incidence angle of unmanned aerial vehicle by using flight and wind tunnel test data
CN104459193A (en) * 2014-12-05 2015-03-25 中国航天空气动力技术研究院 Crosswind information estimation method based on unmanned aerial vehicle crabbing method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565451A (en) * 2011-12-28 2012-07-11 中国航空工业集团公司第六三一研究所 Method for measuring and calculating navigation wind direction and wind speed of general-purpose airplane
CN103852081A (en) * 2014-03-20 2014-06-11 南京航空航天大学 Vacuum speed resolving method for air data/serial inertial navigation combined navigation system
CN103994748A (en) * 2014-05-27 2014-08-20 中国航天空气动力技术研究院 Method for estimating trim incidence angle of unmanned aerial vehicle by using flight and wind tunnel test data
CN104459193A (en) * 2014-12-05 2015-03-25 中国航天空气动力技术研究院 Crosswind information estimation method based on unmanned aerial vehicle crabbing method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
刘让贤,晏初宏主编: "《航空概论》", 31 July 2013 *
李智斌等: "平流层飞艇动力学与控制研究进展", 《力学进展》 *
苗景刚等: "基于Matlab/Simulink的飞艇全数字仿真系统", 《系统仿真学报》 *
苗景刚等: "风场环境下平流层飞艇运动建模", 《系统科学与数学》 *
蒋坤培等: "复杂大气环境下大型飞艇建模研究", 《计算机仿真》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106250581A (en) * 2016-07-13 2016-12-21 南京航空航天大学 A kind of stratosphere aerostatics multi-objective method based on NSGA II
CN106250581B (en) * 2016-07-13 2019-06-04 南京航空航天大学 A multi-objective optimization decision-making method for stratospheric aerostat based on NSGAⅡ
CN106227060A (en) * 2016-07-19 2016-12-14 中国航天员科研训练中心 A kind of stratospheric airship load cabin electronic equipment Thermal characteristic simulation and simulation and prediction System and method for
CN106781706A (en) * 2016-12-16 2017-05-31 北京航空航天大学 Air traffic Track Design method based on wind field distribution
CN106781706B (en) * 2016-12-16 2019-10-18 北京航空航天大学 Air Traffic Route Design Method Based on Wind Field Distribution
CN106681337A (en) * 2017-01-05 2017-05-17 烟台南山学院 Height keeping flight control method for stratospheric airship based on odd-order sliding mode
CN108195558A (en) * 2017-12-06 2018-06-22 太原航空仪表有限公司 A kind of ascending air measuring method for weather modification aircraft
CN108195558B (en) * 2017-12-06 2019-09-13 太原航空仪表有限公司 A kind of ascending air measurement method for weather modification aircraft
CN108388263A (en) * 2018-03-01 2018-08-10 中国航空无线电电子研究所 Circular course autoflight system
CN109917493A (en) * 2019-04-03 2019-06-21 西安邮电大学 Device and method for measuring high-altitude wind energy resources
CN110548288A (en) * 2019-09-05 2019-12-10 腾讯科技(深圳)有限公司 Virtual object hit prompting method and device, terminal and storage medium
CN113515901A (en) * 2021-04-08 2021-10-19 中国科学院空天信息创新研究院 Stability evaluation method and device for stratospheric airship, electronic equipment and medium
CN113515901B (en) * 2021-04-08 2025-01-07 中国科学院空天信息创新研究院 Stability assessment method, device, electronic equipment and medium for stratospheric airship

Also Published As

Publication number Publication date
CN105718727B (en) 2018-09-28

Similar Documents

Publication Publication Date Title
CN105718727B (en) A kind of method and system of estimation stratospheric airship flight performance parameter
CN103852081B (en) Vacuum speed resolving method for air data/serial inertial navigation combined navigation system
CN108152529A (en) A kind of method based on flight parameter calculation of wind speed and wind direction
CN109612676B (en) Pneumatic parameter back calculation method based on flight test data
CN107544532B (en) Long-range high-altitude flight mission planning method for low-altitude airship
Karlgaard et al. Mars entry atmospheric data system trajectory reconstruction algorithms and flight results
CN109240323B (en) A real-time analytical structure for reentry guidance method of aerospace vehicle
CN103303495A (en) Method for estimating disturbance moment in power decreasing process
CN102902204B (en) A kind of dirigible remote control is taken off/landing method for early warning
Park Wind and airspeed error estimation with GPS and pitot-static system for small UAV
CN106705996A (en) Aircraft navigation information correcting method based on atmospheric feature parameters
CN108592911A (en) A kind of quadrotor kinetic model/airborne sensor Combinated navigation method
CN105005099A (en) Atmospheric parameter calculation method based on strapdown inertial navigation and flight control system
CN108225323B (en) Method, medium, and apparatus for determining landing zone boundaries based on deviation influence direction combinations
US10876920B1 (en) Auxiliary aerial vehicles for flow characterization
Marzouk Coupled differential-algebraic equations framework for modeling six-degree-of-freedom flight dynamics of asymmetric fixed-wing aircraft
CN103034125A (en) Pneumatic and jet stream control moment parameter identification method for airship re-entry capsule
CN115683544A (en) Unmanned aerial vehicle rotor disturbance correction method and device
CN109492237A (en) A kind of Aerodynamic Coefficient preparation method
CN104677634B (en) Crosswind ground experiment method on aero-engine machine
CN116227013A (en) Aircraft energy consumption calculation method considering wind field influence
Ostler et al. Flight testing of small, electric powered unmanned aerial vehicles
RU2744772C1 (en) Method for determining wind speed and direction using unmanned aircraft
Qu et al. Real-time wind vector estimation for a micro UAV
Cho et al. Air data system calibration using GPS velocity information

Legal Events

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