CN110346605B - 用于基于静压误差修正进行飞机空速校准的方法以及系统 - Google Patents
用于基于静压误差修正进行飞机空速校准的方法以及系统 Download PDFInfo
- Publication number
- CN110346605B CN110346605B CN201910708953.6A CN201910708953A CN110346605B CN 110346605 B CN110346605 B CN 110346605B CN 201910708953 A CN201910708953 A CN 201910708953A CN 110346605 B CN110346605 B CN 110346605B
- Authority
- CN
- China
- Prior art keywords
- static pressure
- mod
- error correction
- model
- pressure
- 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.)
- Active
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 207
- 238000012937 correction Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000012545 processing Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 14
- 239000003570 air Substances 0.000 description 12
- 238000004590 computer program Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 241000255925 Diptera Species 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000284 resting effect Effects 0.000 description 2
- 238000013515 script Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 241000084490 Esenbeckia delta Species 0.000 description 1
- 239000012080 ambient air Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
- G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本公开提供了一种用于基于静压误差修正模型进行飞机空速校准的方法,包括:获取所基本测量数据;建立静压误差修正模型;将静压误差修正模型的参数初始化;根据静压误差修正模型计算得到模型静温;计算得到模型静压;根据计算得到的模型静压与测量静压得到计算静压误差修正模型;确定所述静压误差修正模型;以及基于确定的静压误差修正模型来修正测量静压,从而根据修正后的测量静压获得远场静压来完成空速校准。本公开还提供了相应的系统以及计算机可读介质。
Description
技术领域
本发明涉及飞行数据的测量与校准,更具体地,涉及用于基于静压误差修正模型进行飞机空速校准的方法,以及相应的系统以及计算机存储介质。
背景技术
大气数据系统对于成功完成飞机飞行任务至关重要,并且可以根据飞机周围的飞行环境而获得。在大气数据系统飞行数据的测量和校准中,大气数据系统飞行数据一般包括所指示的和真实的空速、气压高度、环境空气温度、迎角和侧滑角、以及马赫数以及爬升率。通常,感测到的总压(皮托管压力)和静压会被转换(通过仪器本身中的机械装置)成为高度计、垂直速度指示器、空速指示器和马赫计的指示值。校准可以提供飞机的航空电子设备或飞行控制系统所需的其他参数,诸如真实空速。
飞机飞行中的空速计算关系,如下式:
其中Pt为总压,Ps为静压(测量静压或远场静压),Psea为海平面静压(即,101325Pa),asea为海平面音速(即,340.2940m/s)。通过上述公式,当Pt和Ps为真实飞行的总压和静压时,得到的速度v为真实空速,当Pt和Ps为飞机上传感器所测量的值时,则得到的速度v为指示空速。飞机在飞行中由于与气流的作用会导致测量仪器的输入误差,也即,飞机扰乱了它飞过的空气,从而也扰乱了飞行数据的测量。因此,指示空速存在误差。
在某种程度上,这些误差可以在风洞中进行研究,但风洞测量不能取代实际飞行中的测量。
为了安全飞行,有必要获得准确的飞行数据。其中,总压(皮托管压力)通常易于准确测量,并且总压的误差通常也可以忽略不计。然而静压的压力分布会随飞行条件变化而变化,局部测量的静压与环境静压之间的差异(取决于迎角、空速和飞机构型)被称为位置误差。另外,不受干扰环境下的外部大气温度(OAT)只能在飞机上以极低的速度直接测量。因此,空速校准的主要目的即为校准静压。
现有的技术中采用基于气压法来进行飞机空速校准。气压法原理是根据几何高度和大气压力之间的关系获得远场静压(大气压强)Psatmos。而通常飞机正常飞行时,总压和总温测量误差可忽略,即测量总温等于真实总温,测量总压等于真实总压,因而,空速的主要误差来源于感测得到的测量静压Psmeas与真实静压Pstrue之间存在的偏差-δPs(即,静压误差修正),该偏差会随着飞行高度的增加而不断增加,也即,根据测量静压计算得到的远场静压与真实的远场静压之间的误差是会累积增加的。因此,现有的气压法会降低空速校准中的精度。此外,现有的气压法只能用于低空飞行的飞机的空速校准。
发明内容
为了克服现有技术的缺陷,本发明提供一种用于基于静压误差修正模型进行飞机空速校准的方法,该方法提高了空速校准气压法的精度,并且可以对高空飞行的飞机进行空速校准。
根据各种实施例,提出了一种用于基于静压误差修正模型进行飞机空速校准的方法,该方法包括以下步骤:获取所在飞行高度处的测量总压、测量静压、测量总温、迎角以及在起飞处的场压;建立静压误差修正模型δPsmod,其中上述静压误差修正模型δPsmod是由上述测量总压、上述测量静压、上述迎角构造的三元回归方程;将上述静压误差修正模型δPsmod的参数初始化;根据上述测量总压、上述测量静压、上述测量总温和上述静压误差修正模型δPsmod计算得到模型静温Tsmod;根据计算得到的上述模型静温Tsmod,计算得到模型静压PsInteger mod;根据计算得到的上述模型静压PsInteger mod与上述测量静压得到计算静压误差修正模型δPsInteger mod;根据上述计算静压误差修正模型δPsInteger mod与上述静压误差修正模型函数δPsmod一致性关系回归计算,利用最小二乘法确定上述静压误差修正模型的参数,从而确定上述静压误差修正模型δPsmod;以及基于确定的上述静压误差修正模型δPsmod来修正上述测量静压,从而根据修正后的测量静压获得远场静压来完成空速校准。
根据各种实施例,上述方法还在于,上述静压误差修正模型δPsmod可以根据下式获得:
δPsmod=(a0+a1M+a2M2+a3M3+b1α+b2α2+b3α3+c1Mα+c2Mα2+c3M2α)q
其中,M为计算马赫数,α为迎角,q为动压,以及其中a0、a1、a2、a3、b1、b2、b3、c1、c2、c3为待求解参数。
根据各种实施例,上述方法还在于,上述计算静压误差修正模型δPsInteger mod可以根据下式获得:
δPsInteger mod=PsInteger mod-Psmeas
其中,PsInteger mod为模型静压,Psmeas为测量静压。
根据各种实施例,上述方法还在于,上述模型静温Tsmod可以根据下式获得:
其中,Ptmeas为测量总压,Psmeas为测量静压,Ttmeas为测量总温,δPsmod为静压误差修正模型。
根据各种实施例,上述方法还在于,上述动压可以根据下式获得:
q=0.7PsmeasM2
其中,M为计算马赫数,α为迎角,q为动压,Psmeas为测量静压。
根据各种实施例,上述方法还在于,上述计算马赫数M可以根据下式获得:
其中,Ptmeas为测量总压,Psmeas为测量静压。
根据各种实施例,上述方法还在于,上述模型静压PsInteger mod可以根据下式获得:
其中,Tsmod为模型静温,H为试飞高度、P0为飞机起飞时的场压、g0为重力加速度、R为大气的气体常数。
根据各种实施例,本发明还提供了一种用于基于静压误差修正模型进行飞机空速校准的系统,该系统包括:地面气压传感器,用于测量获得场压P0;差分GPS,用于测量获得几何高度H;总压传感器,用于获得测量总压Ptmeas;静压传感器,用于获得测量静压Psmeas;总温传感器,用于获得测量总温Ttmeas;以及静压误差修正装置,用于根据获得的上述场压P0、上述几何高度H、上述测量总压Ptmeas、上述测量静压Psmeas、上述测量总温Ttmeas来获得静压误差修正模型δPsmod,其中,上述地面气压传感器、上述差分GPS、上述总压传感器、上述静压传感器、上述总温传感器将各自所测量的值传输给上述静压误差修正装置,其中,上述静压误差修正装置被配置为执行本发明中任一项所述的方法。
根据各种实施例,本发明提供的系统还在于,上述总压传感器、上述静压传感器、上述总温传感器均布置于飞机上。
根据各种实施例,本发明提供的系统还在于,上述静压误差修正装置被配置为选自以下中的一个或多个:处理器、处理设备、微处理器。
根据各种实施例,本发明还提供了一种计算机可读介质,包括计算机可执行程序代码,上述计算机可执行程序代码在被执行时使数字处理设备执行本发明中任一项所述的方法。
本发明的积极进步效果在于:
根据本发明的用于基于静压误差修正模型进行飞机空速校准的方法,同传统的气压法相比,本发明专利采用静压修正模型,提高了传统气压法的精度,并且可以对高空飞行的飞机进行空速校准。
附图说明
图1为根据本发明优选实施例的用于基于静压误差修正模型进行飞机空速校准的系统的示意图。
图2为根据本发明优选实施例的用于基于静压误差修正模型进行飞机空速校准的方法200的流程图。
图3为根据本发明优选实施例的用于基于静压误差修正模型进行飞机空速校准的方法300的流程图。
图4为采用本发明优选实施例的用于基于静压误差修正模型进行飞机空速校准的方法与传统气压法积分值对比曲线图。
具体实施方式
下面结合说明书附图,进一步对本发明的优选实施例进行详细描述,以下的描述为示例性的,并非对本发明的限制,任何的其他类似情形也都将落入本发明的保护范围之中。
参考图1所示,本发明示例性地提供了一种用于基于静压误差修正模型进行飞机空速校准的系统100,该系统100包括:地面气压传感器101,用于测量获得场压P0;差分GPS102,用于测量获得几何高度H;总压传感器103,用于获得测量总压Ptmeas;静压传感器104,用于获得测量静压Psmeas;总温传感器105,用于获得测量总温Ttmeas;以及静压误差修正装置106,用于根据前述的场压P0、几何高度H、测量总压Ptmeas、测量静压Psmeas、测量总温Ttmeas来获得静压误差修正模型δPsmod。其中,可以使用当地气象站代替地面气压传感器101来测量场压P0。其中,总压传感器103、静压传感器104、总温传感器105均布置于飞机上,并可以根据下式来获得测量静温Tsmeas:
与之对应的,理想的真实环境静温Tsture为:
其中,Tttrue为真实总温;Pttrue为真实总压;Pstrue为真实静压。
当根据上述(1)或(2)式获得静温Ts后,可以根据下式获得在试飞高度处的远场静压:
其中,Psatmos(H)为试飞高度处远场静压;H为试飞高度;P0为飞机起飞时的场压(通过机场的气压站测量或通过飞机上的静压传感器测量);g0为重力加速度(一般地,g0=9.80665m/s2=32.17405ft/s2);R为大气的气体常数,(一般地,R=287.06(N·m)/(kg·K));Ts为环境静温。当Ts为真实的静温Tsture时,则得到的是真实的试飞高度处远场静压,也即修正期望逼近的值。
另外,上述静压误差修正装置106可以被配置为诸如各种处理器、处理设备、微处理器等。
优选地,上述静压误差修正装置106可以确定测量静压的修正关系,之后结合测量总温、测量总压、修正后的测量静压来获得修正后的测量静温,从而可以得到试飞高度处的远场静压。
优选地,上述静压误差修正装置106可以被配置为执行本发明所述的方法来得到试飞高度处的远场静压。
参考图2所示,根据本发明的优选实施例,提供了一种用于基于静压误差修正模型进行飞机空速校准的方法200,该方法200包括以下的步骤:
步骤202,获取基本测量数据,即,获取所在飞行高度处的测量总压、测量静压、测量总温、迎角以及在起飞处的场压。
步骤204,根据气动理论,建立静压误差修正模型δPsmod,如下式:
δPsmod=f(M,α,q) (4)
其中,M为计算马赫数,α为迎角,q为动压。其中,计算马赫数M可以根据下式获得:
其中,Ptmeas为测量总压,Psmeas为测量静压。
上述动压q可以根据下式获得:
q=0.7PsmeasM2 (6)
其中,Psmeas为测量静压,M为计算马赫数。
其中,迎角α由机载仪器测量获得。
其中,f(M,α,q)可以是标准回归方程;优先地,f(M,α,q)可以是三元线性回归方程。
步骤206,对上述静压误差修正模型初始化,也即令δPsmod=f(M,α,q)=0。
步骤208,获得模型静温Tsmod,也即利用上述初始化的静压误差修正模型δPsmod修正测量静压Psmeas,得到Psmeas+δPsmod,并根据测量总温Ttmeas和测量总压Ptmeas、并参照(1)式计算得到模型静温Tsmod,如下式:
步骤210,根据模型静温Tsmod,试飞高度H、飞机起飞时的场压P0、重力加速度g0、大气的气体常数R、并参照(3)式得到模型静压PsInteger mod,如下式:
步骤212,根据上述的模型静压PsInteger mod和测量静压Psmeas确定计算静压误差模型δPsInteger mod,即存在下式关系:
δPsInteger mod=PsInteger mod-Psmeas (9)
步骤214,根据回归计算获得静压误差修正模型δPsmod的模型参数,例如,其中δPsmod应与δPsInteger mod一致,即Ymod=δPsInteger mod-δPsmod→0。优选地,采用最小二乘法求解模型参数。
步骤216,根据求解后的静压误差修正模型δPsmod,以及上述(7)和(8)式,得到修正后的模型静压PsInteger mod,从而根据空速计算关系得到逼近真实空速的空速值。
同传统的方法相比,本发明的方法会提高传统空速校准气压法的精度。
根据本发明的优选实施例,提供了另一种用于基于静压误差修正模型进行飞机空速校准的方法300,该方法300包括以下的步骤:
步骤302,获取基本测量数据,即,获取所在飞行高度处的测量总压、测量静压、测量总温、迎角以及在起飞处的场压。
步骤304,类似方法200中的步骤204,建立静压误差修正模型δPsmod,如下式:
δPsmod=(a0+a1M+a2M2+a3M3+b1α+b2α2+b3α3+c1Mα+c2Mα2+c3M2α)q
其中,M为计算马赫数,α为迎角,q为动压。其中,计算马赫数M可以根据(5)式获得;动压q可以根据(6)式获得;迎角α由测量获得;并且其中a0、a1、a2、a3、b1、b2、b3、c1、c2、c3为待求解参数,也即回归方程的回归系数。
步骤306,类似方法200中的步骤206,对上述静压误差修正模型初始化,也即令上述待求解参数a0、a1、a2、a3、b1、b2、b3、c1、c2、c3均为0。
步骤308,类似方法200中的步骤208,获得模型静温Tsmod,也即利用上述初始化的静压误差修正模型δPsmod修正测量静压Psmeas,得到Psmeas+δPsmod,并根据测量总温Ttmeas和测量总压Ptmeas、并根据(7)式计算得到模型静温Tsmod。
步骤310,类似方法200中的步骤210,根据模型静温Tsmod,试飞高度H、飞机起飞时的场压P0、重力加速度g0、大气的气体常数R、并根据(8)式得到模型静压PsInteger mod。
步骤312,类似方法200中的步骤212,根据上述的模型静压PsInteger mod和测量静压Psmeas确定计算静压误差修正模型δPsInteger mod,即δPsInteger mod=PsInteger mod-Psmeas。
步骤314,类似方法200中的步骤214,根据回归计算获得模型参数,例如,其中δPsmod应与δPsInteger mod一致,即Ymod=δPsInteger mod-δPsmod→0。优选地,采用最小二乘法求解参数a0、a1、a2、a3、b1、b2、b3、c1、c2、c3。
步骤316,类似方法200中的步骤216,根据确定了参数的静压误差修正模型δPsmod,以及上述(7)和(8)式,得到修正后的模型静压PsInteger mod。
根据上述各个实施例,优选地可以在飞机上设置三套或三套以上的用于基于静压误差修正模型进行飞机空速校准的方法的系统,其中静压传感器为相互独立的,其他传感器或飞行数据收集装置可以共用或独立设置。优选地,上述系统中的静压传感器设置于机头与左右翼。优选地,上述系统中的静压传感器设置于尾翼与机体下部。采用三套上述系统不但可以验证本发明的正确性,也可以避免系统失效或个别部位对测量飞行数据带来的系统误差或偶然误差。
根据本发明优选实施例,还提供了一种静压误差修正的方法,该方法包括以下步骤:试飞数据准备,根据机载仪器或地面设备获取场压P0、几何高度H、测量总压Ptmeas、测量静压Psmeas、测量总温Ttmeas;选择静压误差修正模型δPsmod;将静压误差修正模型参数初始化;计算基于该模型的静温Tsmod;并根据计算的静温来计算模型静压PsInteger mod;计算静压误差修正模型δPsInteger mod=PsInteger mod-Psmod;根据静压误差回归函数Y=δPsInteger mod-δPsmod,利用最小二乘法确定静压误差修正模型参数;确定静压误差修正模型;对测量静压修正;从而根据修正后的测量静压获得远场静压来完成空速校准。
图4为采用本发明优选实施例的用于基于静压误差修正模型进行飞机空速校准的方法与传统气压法积分值对比曲线图。图中实线为采用本发明方法根据某型号飞机三套大气数据系统计算得到的静压真值,虚线为传统气压法得到的静压真值。理论上,分别采用三套的静压真值应该相等,即图中曲线理论值应为0。如图4所示,采用本发明的方法比传统气压法在精度上明显提高,并且传统气压法只能用于低空飞行空速校准,采用本发明的方法,可以将适用范围扩大到高空飞行中的空速校准。
此外,根据本发明的实施例,可以提供一种计算机可读介质,包括计算机可执行程序代码,所述计算机可执行程序代码在被执行时使数字处理设备执行上述实施例中所使用的方法以及步骤。
本公开内容的实施例和本文中提供的所有功能操作可以用数字电子电路、或者用计算机软件、固件或硬件,包括本说明书及其结构等同方案中所公开的结构、或者其中的一个或多个的组合来实现。本公开内容的实施例可以实施为一个或多个计算机程序产品,即在计算机可读介质上编码的计算机程序指令的一个或多个模块,这些指令由数据处理装置来执行或者用以控制数据处理装置的操作。该计算机可读介质可以是机器可读存储设备、机器可读存储基片、存储器设备、影响机器可读传播信号的组合物或者其中的一个或多个的组合。术语“数据处理装置”涵盖用于处理数据的所有装置、设备和机器,包括例如可编程处理器、计算机或者多个处理器或计算机。除了硬件之外,该装置可以包括为所讨论的计算机程序创建执行环境的代码,例如构成处理器固件、协议栈、数据库管理系统、操作系统或者其中的一个或多个的组合的代码。
计算机程序(也称为程序、软件、软件应用、脚本或代码)可以用任何形式的编程语言(包括编译语言或解释语言)来编写,并且计算机程序可以用任何形式来部署,包括作为独立程序或者作为模块、部件、子例程或者适合在计算环境中使用的其他单元。计算机程序并非必须对应于文件系统中的文件。程序可以存储在保持其他程序或数据(例如标记语言文档中所存储的一个或多个脚本)的文件的部分中,存储在专用于所讨论的程序的单个文件中,或者存储在多个协同文件(例如存储一个或多个模块、子程序或者代码的部分的文件)中。计算机程序可以被部署成在一个计算机上来执行,或者在位于一个站点处或分布在多个站点处且通过通信网络互连的多个计算机上来执行。
本公开内容中所描述的过程和逻辑流可以由执行一个或多个计算机程序的一个或多个可编程处理器来执行以通过操作输入数据并且生成输出来执行功能。该过程和逻辑流也可以由专用逻辑电路来执行,并且装置也可以实施为该专用逻辑电路,该专用逻辑电路例如为FPGA(现场可编程门阵列)或者ASIC(专用集成电路)。
适合执行计算机程序的处理器包括例如通用和专用微处理器二者、以及任何种类的数字计算机的任何一个或多个处理器。通常,处理器从只读存储器或者随机存取存储器或者二者接收指令和数据。计算机的元件可以包括用于执行指令的处理器以及用于存储指令和数据的一个或多个存储器设备。通常,计算机还将包括一个或多个海量存储设备以便存储数据,或者该计算机在操作上耦合以从海量存储设备接收或向海量存储设备传送数据或者二者,该海量存储设备例如是磁盘、磁光盘或者光盘。然而,计算机不需要具有这样的设备。此外,计算机可以嵌入在另一设备中,该另一设备例如为移动电话、个人数字助理(PDA)、移动音频播放器、全球定位系统(GPS)接收器等。适合存储计算机程序指令和数据的计算机可读介质包括所有形式的非易失性存储器、介质和存储器设备,包括例如:半导体存储器设备,如EPROM、EEPROM和闪存设备;磁盘,如内置硬盘或可移除盘;磁光盘;以及CD ROM和DVD-ROM盘。该处理器和存储器可以用专用逻辑电路来补充或者并入该专用逻辑电路中。
虽然本公开内容包括一些细节,然而不应当将这些细节理解为对本公开内容或者要求保护的内容的范围的限制,而是应当被理解为对本公开内容的示例实施例的特征的描述。本公开内容中在单独实施例的情境中描述的某些特征还可以与单个实施例组合来提供。相反地,在单个实施例的情境中描述的各个特征也可以分别在多个实施例中来提供或者在任何合适的子组合中来提供。此外,虽然以上可以将特征描述为以某种组合来执行并且甚至初始就要求这样保护,然而在一些情况下可以从组合中去掉来自要求保护的组合的一个或多个特征,并且要求保护的组合可以涉及子组合或子组合的变化。
类似地,虽然在附图中按照特定顺序来描绘操作,然而这不应当被理解为要求这样的操作按照所示的特定顺序或者按照相继顺序来执行,或者要求所有图示操作都被执行,以实现期望的结果。在一些境况下,多任务和并行处理可能是有利的。此外,以上描述的
实施例中的各种系统部件的分离不应当被理解为在所有实施例中都要求这样的分离,而且应当理解,所描述的程序部件和系统通常可以在单个软件产品中集成在一起或者被封装成多个软件产品。
因此,本公开内容的特定实施例已经被描述,并且其他实施例在以下权利要求的范围内。例如,权利要求中记载的动作可以按照不同的顺序来执行,并且这些动作仍然可以实现期望的结果。大量实施例已经被描述。然而,应当理解,可以在不偏离本公开内容的精神和范围的情况下做出各种修改。例如,可以使用以上示出的流程的各个形式,其中步骤可以被重新排序、添加或去除。因此,其他实施例在以下权利要求的范围内。
虽然以上描述了本发明的具体实施方式,但是本领域的技术人员应当理解,这些仅是举例说明,本发明的保护范围是由所附权利要求书限定的。本领域的技术人员在不背离本发明的原理和实质的前提下,可以对这些实施方式做出多种变更或修改,但这些变更和修改均落入本发明的保护范围。
Claims (7)
1.一种用于基于静压误差修正模型进行飞机空速校准的方法,所述方法包括以下步骤:
获取所在飞行高度处的测量总压、测量静压、测量总温、迎角以及在飞机起飞时的场压;
建立静压误差修正模型δPsmod,其中所述静压误差修正模型δPsmod是由所述测量总压、所述测量静压、所述迎角构造的三元回归方程;
将所述静压误差修正模型δPsmod的参数初始化;
根据所述测量总压、所述测量静压、所述测量总温和所述静压误差修正模型δPsmod计算得到模型静温Tsmod;
根据计算得到的所述模型静温Tsmod,计算得到模型静压PsInteger mod;
根据计算得到的所述模型静压PsInteger mod与所述测量静压得到计算静压误差修正模型δPsInteger mod;
根据所述计算静压误差修正模型δPsInteger mod与所述静压误差修正模型δPsmod一致性关系回归计算,利用最小二乘法确定所述静压误差修正模型的参数,从而确定所述静压误差修正模型δPsmod;以及
基于确定的所述静压误差修正模型δPsmod来修正所述测量静压,从而根据修正后的所述测量静压获得远场静压来完成空速校准;
其中,所述静压误差修正模型δPsmod可以根据下式获得:
δPsmod=(a0+a1M+a2M2+a3M3+b1α+b2α2+b3α3+c1Mα+c2Mα2+c3M2α)q
其中,M为计算马赫数,α为迎角,q为动压,以及其中a0、a1、a2、a3、b1、b2、b3、c1、c2、c3为待求解参数;
其中,所述计算静压误差修正模型δPsInteger mod可以根据下式获得:
δPsInteger mod=PsInteger mod-Psmeas
其中,PsInteger mod为模型静压,Psmeas为测量静压;
其中,所述模型静温Tsmod可以根据下式获得:
其中,Ptmeas为测量总压,Psmeas为测量静压,Ttmeas为测量总温,δPsmod为静压误差修正模型;
其中,所述模型静压PsInteger mod可以根据下式获得:
其中,Tsmod为模型静温,H为试飞高度、P0为飞机起飞时的场压、g0为重力加速度、R为大气的气体常数。
2.根据权利要求1所述的方法,其中,所述动压可以根据下式获得:
q=0.7PsmeasM2
其中,M为计算马赫数,α为迎角,q为动压,Psmeas为测量静压。
4.一种用于基于静压误差修正模型进行飞机空速校准的系统,所述系统包括:
地面气压传感器,用于测量获得场压P0;
差分GPS,用于测量获得几何高度H;
总压传感器,用于获得测量总压Ptmeas;
静压传感器,用于获得测量静压Psmeas;
总温传感器,用于获得测量总温Ttmeas;以及
静压误差修正装置,用于根据获得的所述场压P0、所述几何高度H、所述测量总压Ptmeas、所述测量静压Psmeas、所述测量总温Ttmeas来获得静压误差修正模型δPsmod,
其中,所述地面气压传感器、所述差分GPS、所述总压传感器、所述静压传感器、所述总温传感器将各自所测量的值传输给所述静压误差修正装置,
其中,所述静压误差修正装置被配置为能够执行根据权利要求1-3中任一项所述的方法。
5.根据权利要求4所述的系统,其中,所述总压传感器、所述静压传感器、所述总温传感器均布置于飞机上。
6.根据权利要求4或5所述的系统,其中,所述静压误差修正装置被配置为选自以下中的一个或多个:处理器、处理设备、微处理器。
7.一种计算机可读介质,包括计算机可执行程序代码,所述计算机可执行程序代码在被执行时使数字处理设备执行如权利要求1-3中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910708953.6A CN110346605B (zh) | 2019-08-01 | 2019-08-01 | 用于基于静压误差修正进行飞机空速校准的方法以及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910708953.6A CN110346605B (zh) | 2019-08-01 | 2019-08-01 | 用于基于静压误差修正进行飞机空速校准的方法以及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110346605A CN110346605A (zh) | 2019-10-18 |
CN110346605B true CN110346605B (zh) | 2021-05-07 |
Family
ID=68183714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910708953.6A Active CN110346605B (zh) | 2019-08-01 | 2019-08-01 | 用于基于静压误差修正进行飞机空速校准的方法以及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110346605B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111268165B (zh) * | 2019-11-15 | 2023-10-20 | 陕西飞机工业(集团)有限公司 | 一种中型飞机静压源高度误差的修正方法 |
EP3919869B1 (en) * | 2020-06-04 | 2023-08-02 | Stichting Nationaal Lucht- en Ruimtevaart Laboratorium | Method for determining a position error correction on a static pressure measurement at an aircraft |
CN112730887B (zh) * | 2020-12-29 | 2023-07-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种全压失效时马赫数的获取方法和系统 |
CN113465576B (zh) * | 2021-09-06 | 2021-11-19 | 中国商用飞机有限责任公司 | 基于飞行器的gnss高度来算出气压高度的方法和系统 |
CN113533784B (zh) * | 2021-09-07 | 2022-01-25 | 成都飞机工业(集团)有限责任公司 | 一种gps往返非等速平飞空速校准方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246078A (zh) * | 2008-03-18 | 2008-08-20 | 北京航空航天大学 | 应用于低速无人机上总静压受感器位置误差的修正方法 |
US7505145B2 (en) * | 2002-03-01 | 2009-03-17 | Michigan Aerospace Corporation | Optical air data system |
CN102050226A (zh) * | 2009-10-30 | 2011-05-11 | 航天科工惯性技术有限公司 | 一种航空应急仪表及其系统初始对准方法和组合导航算法 |
CN102360027A (zh) * | 2011-08-02 | 2012-02-22 | 中国科学院自动化研究所 | 差压式空速传感器测试系统 |
CN102419381A (zh) * | 2010-09-27 | 2012-04-18 | 波音公司 | 用于飞行器的空速感测系统 |
EP2810081A1 (en) * | 2012-02-03 | 2014-12-10 | Jason Thomas | Methods and systems for determining airspeed of an aircraft |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61155719A (ja) * | 1984-12-21 | 1986-07-15 | ザ ギヤレツト コーポレーシヨン | エラー補正装置 |
US9285387B2 (en) * | 2009-12-14 | 2016-03-15 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | In-flight pitot-static calibration |
CN101913427B (zh) * | 2010-08-04 | 2013-01-09 | 北京航空航天大学 | 一种适用于多用途无人飞行器的航空电子系统 |
CN103076013B (zh) * | 2012-12-27 | 2016-07-13 | 太原航空仪表有限公司 | 用于飞行导航的大气数据与姿态航向基准系统 |
CN103759884B (zh) * | 2014-01-29 | 2016-03-02 | 中国商用飞机有限责任公司 | 测量高压管路中高速流体的静压的装置和方法 |
CN103852398B (zh) * | 2014-03-12 | 2017-01-04 | 上海夏普电器有限公司 | 一种基于温度修正的集尘滤网静压差失效检测方法 |
CN105353761B (zh) * | 2015-12-02 | 2018-06-12 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种飞机的速度控制方法及装置 |
CN105973544B (zh) * | 2016-06-20 | 2018-04-10 | 中国民航大学 | 一种基于自适应pid方法的飞机全静压测试系统及方法 |
CN106324643A (zh) * | 2016-10-19 | 2017-01-11 | 山东科技大学 | 一种无人机空速估计和空速管故障检测方法 |
CN108090253B (zh) * | 2017-11-29 | 2019-02-26 | 中国直升机设计研究所 | 一种直升机数字式大气系统空速修正方法 |
-
2019
- 2019-08-01 CN CN201910708953.6A patent/CN110346605B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7505145B2 (en) * | 2002-03-01 | 2009-03-17 | Michigan Aerospace Corporation | Optical air data system |
CN101246078A (zh) * | 2008-03-18 | 2008-08-20 | 北京航空航天大学 | 应用于低速无人机上总静压受感器位置误差的修正方法 |
CN102050226A (zh) * | 2009-10-30 | 2011-05-11 | 航天科工惯性技术有限公司 | 一种航空应急仪表及其系统初始对准方法和组合导航算法 |
CN102419381A (zh) * | 2010-09-27 | 2012-04-18 | 波音公司 | 用于飞行器的空速感测系统 |
CN102360027A (zh) * | 2011-08-02 | 2012-02-22 | 中国科学院自动化研究所 | 差压式空速传感器测试系统 |
EP2810081A1 (en) * | 2012-02-03 | 2014-12-10 | Jason Thomas | Methods and systems for determining airspeed of an aircraft |
Also Published As
Publication number | Publication date |
---|---|
CN110346605A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110346605B (zh) | 用于基于静压误差修正进行飞机空速校准的方法以及系统 | |
US11085845B2 (en) | Systems and methods for correcting acoustic error in pressure sensors | |
US10352813B2 (en) | Calibration of aircraft instruments using a laser sensor | |
US9342988B2 (en) | Method and device for determining a linear terrain profile along a lateral approach trajectory of an airport | |
US9108745B2 (en) | Method for detecting a failure of at least one sensor onboard an aircraft implementing an anemo-inertial loop, and associated system | |
US9753052B2 (en) | Method and device for estimating the mach number of an aircraft | |
CN116380396B (zh) | 一种连续式风洞飞行器大气数据系统鉴定试验系统及方法 | |
Raab et al. | In-flight testing of MEMS pressure sensors for flight loads determination | |
US20050159857A1 (en) | System and method for determining aircraft tapeline altitude | |
US10921457B2 (en) | Method and device for detecting a fault of a barometric pressure measuring system arranged aboard a flying device | |
EP3919869B1 (en) | Method for determining a position error correction on a static pressure measurement at an aircraft | |
CN114778887A (zh) | 一种基于改进三角矢量模型的无人机测风方法和装置 | |
Cho et al. | Air data system calibration using GPS velocity information | |
EP1750134B1 (en) | Method using ratiometric characterization to improve accuracy of measurements in an air data system | |
Jain et al. | Calibration and accuracy determination of airdata system for a modern fighter aircraft | |
JPS61155719A (ja) | エラー補正装置 | |
Larin et al. | Investigation of the Accuracy of Barometric Pressure Sensors to Assessment of Their Possibility for UAVs Landing | |
Jain et al. | A novel methodology to update table in air data system of a high performance fighter aircraft | |
Efremova et al. | Study of Procedure Errors of the Vortex Air Data System of a Subsonic Aircraft | |
Zhubin et al. | Real-time estimation of altitude and airspeed under large airflow angles in flight test | |
CN117130072A (zh) | 一种用于探测临近空间大气环境的探空仪及方法 | |
Raab | Practical examples for the flight data compatibility check | |
Saraf et al. | Calibration and accuracy determination of airdata system for a modern fighter | |
Nguyen et al. | On the design of commercial off-the-shelf machmeter with CANaerospace interface | |
Raab | Check your Air Data-Setup and Practical Examples for the Flight Data Compatibility Check |
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 |