CN115512461B - 飞机巡航稳态的确定方法、装置以及计算机设备 - Google Patents

飞机巡航稳态的确定方法、装置以及计算机设备 Download PDF

Info

Publication number
CN115512461B
CN115512461B CN202211429957.9A CN202211429957A CN115512461B CN 115512461 B CN115512461 B CN 115512461B CN 202211429957 A CN202211429957 A CN 202211429957A CN 115512461 B CN115512461 B CN 115512461B
Authority
CN
China
Prior art keywords
sequence
data
sum
sub
prefix
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
Application number
CN202211429957.9A
Other languages
English (en)
Other versions
CN115512461A (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.)
South China Normal University
Original Assignee
South China Normal University
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 South China Normal University filed Critical South China Normal University
Priority to CN202211429957.9A priority Critical patent/CN115512461B/zh
Publication of CN115512461A publication Critical patent/CN115512461A/zh
Application granted granted Critical
Publication of CN115512461B publication Critical patent/CN115512461B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G07—CHECKING-DEVICES
    • G07C—TIME OR ATTENDANCE REGISTERS; REGISTERING OR INDICATING THE WORKING OF MACHINES; GENERATING RANDOM NUMBERS; VOTING OR LOTTERY APPARATUS; ARRANGEMENTS, SYSTEMS OR APPARATUS FOR CHECKING NOT PROVIDED FOR ELSEWHERE
    • G07C5/00—Registering or indicating the working of vehicles
    • G07C5/08—Registering or indicating performance data other than driving, working, idle, or waiting time, with or without registering driving, working, idle or waiting time
    • G07C5/0808—Diagnosing performance data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C23/00—Combined instruments indicating more than one navigational value, e.g. for aircraft; Combined measuring devices for measuring two or more variables of movement, e.g. distance, speed or acceleration

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请涉及一种飞机巡航稳态的确定方法、装置以及计算机设备，该方法包括：获取获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态，从而提高了计算效率，节约了计算算力。

Description

飞机巡航稳态的确定方法、装置以及计算机设备
技术领域
本申请涉及航空飞行状态检测领域，特别是涉及一种飞机巡航稳态的确定方法、装置以及计算机设备。
背景技术
飞机完成一次飞行一般包括滑行、起飞、爬升、巡航、下降、进近等多个阶段。巡航稳态是指巡航时飞机飞行及外界因素处于平稳状态。巡航稳态判别对航空数据分析中起重要作用，例如，可通过巡航稳态下的飞机油耗分析给出节油策略，又例如可通过巡航稳态下飞行姿态分析为飞机舵面调整提供支持。
现有的飞机巡航状态确定方法，是对某一段时间窗口内飞机的机载记录数据进行直接计算，由于窗口随着飞机飞行而逐步前进，窗口内的机载记录数据也随之改变，需要对窗口内的机载记录数据进行多次遍历，从而存在数据重复计算，导致计算效率低，浪费计算算力。
发明内容
基于此，本申请的目的在于，提供一种飞机巡航稳态的确定方法、装置以及计算机设备，其可提高计算效率，节约计算算力。
根据本申请实施例的第一方面，提供一种飞机巡航稳态的确定方法，包括如下步骤：
获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；
将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；
根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；
将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
根据本申请实施例的第二方面，提供一种飞机巡航稳态的确定装置，包括：
序列数据获取模块，用于获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；
序列数据划分模块，用于将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；
样本标准差获得模块，用于根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；
巡航稳态确定模块，用于将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
根据本申请实施例的第三方面，提供一种计算机设备，包括：处理器和存储器；其中，所述存储器存储有计算机程序，所述计算机程序适于由所述处理器加载并执行如上述任意一项所述的种飞机巡航稳态的确定方法。
本申请实施例通过获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态，从而提高了计算效率，节约了计算算力。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。
为了更好地理解和实施，下面结合附图详细说明本发明。
附图说明
图1为本申请一个实施例提供的飞机巡航稳态的确定方法的流程示意图；
图2为本申请一个实施例提供的飞机巡航稳态的确定方法中步骤S10的流程示意图；
图3为本申请一个实施例提供的机巡航稳态的确定方法中步骤S30的流程示意图；
图4为本申请一个实施例提供的机巡航稳态的确定方法中步骤S40的流程示意图；
图5为本申请另一个实施例提供的机巡航稳态的确定方法的流程示意图；
图6为本申请另一个实施例提供的机巡航稳态的确定方法中步骤S50的流程示意图；
图7为本申请另一个实施例提供的机巡航稳态的确定方法中步骤S60的流程示意图；
图8为本申请另一个实施例提供的机巡航稳态的确定方法中步骤S70的流程示意图；
图9为本申请一个实施例提供的飞机巡航稳态的确定装置的结构框图；
图10为本申请一个实施例提供的电子设备的结构示意框图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施例方式作进一步地详细描述。
应当明确，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本申请保护的范围。
在本申请实施例使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请实施例。在本申请实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。在本申请的描述中，需要理解的是，术语“第一”、“第二”、“第三”等仅用于区别类似的对象，而不必用于描述特定的顺序或先后次序，也不能理解为指示或暗示相对重要性。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本申请中的具体含义。
此外，在本申请的描述中，除非另有说明，“多个”是指两个或两个以上。“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，A 和/或 B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
实施例1
请参阅图1，其为本申请一个实施例提供的飞机巡航稳态的确定方法的流程示意图。本申请实施例提供的飞机巡航稳态的确定方法，包括如下步骤：
S10：获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和。
其中，机载记录时间序列数据为一段时间窗口内的机载记录数据，例如，从10点10 分到10点16分的机载记录数据。序列前缀和为某个序列的前n项和。具体地，对于序列
Figure 65825DEST_PATH_IMAGE001
，序列前缀和表示为
Figure 249681DEST_PATH_IMAGE002
，
Figure 9696DEST_PATH_IMAGE003
。序列前缀平方和为某个序列的前n项的平方和，序列 前缀平方和表示为
Figure 900291DEST_PATH_IMAGE004
，
Figure 228505DEST_PATH_IMAGE005
。
在本申请实施例中，可以从快速存取记录仪(Quick Access Recorder，简称QAR)中获取机载记录时间序列数据。其中，QAR是一种带有保护装置的机载飞行记录仪，它用于监控、记录大量的时序飞行数据，涵盖了飞行运行品质分析的绝大部分巡航参数。机载记录数据包括飞机的飞行高度、经纬度、飞行加速度以及飞行速度等巡航参数。
S20：将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和。
在本申请实施例中，可以通过滑动窗口将机载记录时间序列数据划分为若干个子序列数据。其中，滑动窗口具有一定的时间长度。例如，通过5分钟的滑动窗口，将10点10分到10点16分的机载记录数据划分为10点10分到10点15分的子序列数据和10点11分到10点16分的子序列数据。
由于机载记录时间序列数据的序列前缀和以及序列前缀平方和是已知的，可以直接获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和。
S30：根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差。
其中，样本标准差用于说明样本数据的离散程度。
在本申请实施例中，在获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和之后，可以直接通过每个所述子序列数据的子序列前缀和以及子序列前缀平方和，计算获得每个所述子序列数据的样本标准差。
S40：将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
在本申请实施例中，将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，若每个所述子序列数据的样本标准差均小于预设的第一阈值，则确定飞机巡航进入稳态。若存在至少一个所述子序列数据的样本标准差大于或等于预设的第一阈值，则确定飞机巡航进入非稳态。
应用本申请实施例，通过获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态，从而提高了计算效率，节约了计算算力。
在一个可选的实施例中，请参阅图2，所述步骤S10计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和的步骤，包括步骤S11~S13，具体如下：
S11：遍历所述机载记录时间序列数据的每个机载记录数据，计算第一个机载记录数据到当前机载记录数据的序列前缀和以及序列前缀平方和，获得第一序列前缀和以及第一序列前缀平方和；
S12：遍历所述当前机载记录数据的下一个机载记录数据，将所述下一个机载记录数据与所述第一序列前缀和求和，获得第二序列前缀和；将所述下一个机载记录数据的平方与所述第一序列前缀平方和求和，获得第二序列前缀平方和；
S13：继续遍历下一个机载记录数据，直至计算所述第一个机载记录数据到最后一个机载记录数据的序列前缀和以及序列前缀平方和。
在本申请实施例中，以机载记录时间序列数据为10点10分到10点16分的机载记录 数据为例进行说明，10点10分到10点16分的机载记录数据包括50个记载记录数据，表示为：
Figure 509444DEST_PATH_IMAGE006
，10点10分到10点15分的子序列数据包括40个记载记录数 据，表示为：
Figure 214095DEST_PATH_IMAGE007
，10点11分到10点16分的子序列数据包括40个 记载记录数据，表示为：
Figure 959197DEST_PATH_IMAGE008
，从第1个机载记录数据开始，分 别计算第1个机载记录数据到第2个机载记录数据的序列前缀和
Figure 597594DEST_PATH_IMAGE009
以及序列前缀平方和
Figure 224884DEST_PATH_IMAGE010
，计算第1个机载记录数据到第3个机载记录数据的序列前缀和
Figure 343013DEST_PATH_IMAGE011
以及序列前缀平方 和
Figure 208201DEST_PATH_IMAGE012
，直至计算所述第1个机载记录数据到第50个机载记录数据的序列前缀和
Figure 878217DEST_PATH_IMAGE013
以及 序列前缀平方和
Figure 383016DEST_PATH_IMAGE014
，表示如下：
Figure 163890DEST_PATH_IMAGE015
Figure 555689DEST_PATH_IMAGE016
从而10点10分到10点15分的子序列数据的序列前缀和即为
Figure 599868DEST_PATH_IMAGE017
，10点10分到10点 15分的子序列数据的序列前缀平方和即为
Figure 201751DEST_PATH_IMAGE018
，10点11分到10点16分的子序列数据的序 列前缀和可以表示为：
Figure 176529DEST_PATH_IMAGE019
，10点11分到10点16分的子序列数据的序列前缀和可 以表示为：
Figure 547467DEST_PATH_IMAGE020
。
通过计算机载记录时间序列数据的序列前缀和以及序列前缀平方和，从而可以自动快捷地获取机载记录时间序列数据中任意数量的机载记录数据的序列前缀和以及序列前缀平方和，以便后续计算子序列数据的样本标准差。
在一个可选的实施例中，请参阅图3，所述步骤S30，包括步骤S31~S32，具体如下：
S31：获取每个所述子序列数据的序列长度；
S32：根据所述序列长度、所述子序列前缀和以及所述子序列前缀平方和，获得每个所述子序列数据的标准差；其中，计算每个所述子序列数据的样本标准差的公式为：
Figure 762548DEST_PATH_IMAGE021
Figure 727093DEST_PATH_IMAGE022
其中，
Figure 115349DEST_PATH_IMAGE023
表示第
Figure 465428DEST_PATH_IMAGE024
个元素到第
Figure 851410DEST_PATH_IMAGE025
个元素且序列长度为
Figure 693464DEST_PATH_IMAGE026
的子序列数据 的样本标准差，
Figure 495198DEST_PATH_IMAGE027
表示第1个元素到第
Figure 778411DEST_PATH_IMAGE028
个元素的子序列前缀和，
Figure 397611DEST_PATH_IMAGE029
表示第1 个元素到第
Figure 588946DEST_PATH_IMAGE028
个元素的子序列前缀平方和，
Figure 53425DEST_PATH_IMAGE030
表示第1个元素到第
Figure 456725DEST_PATH_IMAGE025
个元素的子序列 前缀和，
Figure 387772DEST_PATH_IMAGE031
表示第1个元素到第
Figure 407680DEST_PATH_IMAGE025
个元素的子序列前缀平方和。
在本申请实施例中，若要计算10点10分到10点15分的子序列数据的样本标准差， 获取10点10分到10点15分的子序列数据的序列长度40、子序列前缀和
Figure 144692DEST_PATH_IMAGE017
以及所述子序列 前缀平方和
Figure 323870DEST_PATH_IMAGE018
，从而获得样本标准差：
Figure 284872DEST_PATH_IMAGE032
其中，
Figure 933023DEST_PATH_IMAGE033
和
Figure 208146DEST_PATH_IMAGE034
均为0。
若要计算10点11分到10点16分的子序列数据的样本标准差，获取10点11分到10点 16分的子序列数据的序列长度40、子序列前缀和
Figure 117196DEST_PATH_IMAGE019
以及所述子序列前缀平方和
Figure 639313DEST_PATH_IMAGE020
，从而获得样本标准差：
Figure 633814DEST_PATH_IMAGE035
通过样本标准差计算公式，可以自动快捷地计算出每个子序列数据的样本标准差，提高了计算效率，节约了计算算力。
在一个可选的实施例中，请参阅图4，所述步骤S40，包括步骤S41~S42，具体如下：
S41：若每个所述子序列数据的样本标准差均小于所述预设的第一阈值，则确定飞机巡航进入稳态；
S42：若存在至少一个子序列数据的样本标准差大于或等于所述预设的第一阈值，则确定飞机巡航进入非稳态。
在本申请实施例中，将每个子序列数据的样本标准差与预设的第一阈值进行比较，可以确定每个子序列数据对应的时间窗口内飞机是否进入巡航稳态。具体地，若某个子序列数据的样本标准差小于预设的第一阈值，则该子序列数据对应的时间窗口内飞机巡航进入稳态。
进一步地，可以确定飞机的机载记录时间序列数据对应的时间窗口内飞机是否进入巡航稳态。例如，将飞机的机载记录时间序列数据划分为两个子序列数据，若两个子序列数据对应的时间窗口内飞机均进入巡航稳态，则飞机的机载记录时间序列数据对应的时间窗口内飞机进入巡航稳态。若两个子序列数据对应的时间窗口内存在至少一个时间窗口内飞机进入巡航非稳态，则飞机的机载记录时间序列数据对应的时间窗口内飞机进入巡航非稳态。
通过将每个子序列数据的样本标准差与预设的第一阈值进行比较，可以自动快捷地确定飞机的机载记录时间序列数据对应的时间窗口内飞机是否进入巡航稳态。
在一个可选的实施例中，请参阅图5，飞机巡航稳态的确定方法，还包括步骤S50~S70，具体如下：
S50：对所述机载记录时间序列数据进行预处理，获得最大值稀疏表和最小值稀疏表。
其中，稀疏表（Sparse Table，简称ST表）是一种数据结构，用来解决区间最值查询（Range Minimum Query，简称RMQ）问题。
在本申请实施例中，通过对机载记录时间序列数据进行预处理，生成最大值稀疏表和最小值稀疏表。
S60：根据区间最值查询方法，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值。
其中，区间最值查询方法是指对于长度为n的数列A，查询数列A中元素的最大值或最小值。
在本申请实施例中，通过区间最值查询方法，获得每个所述子序列数据的最大值和最小值。
S70：计算每个所述子序列数据的最大值和最小值的极差，将所述极差与预设的第二阈值比较，获得第二比较结果，根据所述第二比较结果确定飞机巡航稳态。
在本申请实施例中，将每个所述子序列数据的极差与预设的第二阈值进行比较，若每个所述子序列数据的极差均小于预设的第二阈值，则确定飞机巡航进入稳态。若存在至少一个所述子序列数据的极差大于或等于预设的第二阈值，则确定飞机巡航进入非稳态。可选的，将每个所述子序列数据的样本标准差和极差分别与预设的第一阈值和预设的第二阈值比较，若每个所述子序列数据的样本标准差均小于预设的第一阈值，且每个所述子序列数据的极差均小于预设的第二阈值，则确定飞机巡航进入稳态。若存在至少一个所述子序列数据的样本标准差大于或等于预设的第一阈值，或者存在至少一个所述子序列数据的极差大于或等于预设的第二阈值，则确定飞机巡航进入非稳态。
在一个可选的实施例中，请参阅图6，所述步骤S50，还包括步骤S51~S52，具体如下：
S51：根据状态转移方程，确定所述机载记录时间序列数据在预设区间内任意子区间的最大值和最小值；
S52：根据所有子区间的最大值和最小值，获得最大值稀疏表和最小值稀疏表；其中，所述状态转移方程的公式为：
Figure 712629DEST_PATH_IMAGE036
Figure 617131DEST_PATH_IMAGE037
其中，
Figure 919936DEST_PATH_IMAGE038
表示所述机载记录时间序列数据在区间
Figure 791946DEST_PATH_IMAGE039
的最 大值，
Figure 408872DEST_PATH_IMAGE040
表示所述机载记录时间序列数据在区间
Figure 292515DEST_PATH_IMAGE039
的最小值。
在本申请实施例中，以
Figure 376008DEST_PATH_IMAGE041
=1，
Figure 610680DEST_PATH_IMAGE042
=2为例进行说明，表示机载记录时间序列数据在区 间[1，4]的最大值，即序列数据
Figure 955599DEST_PATH_IMAGE043
的最大值，根据状态转移方程，可得
Figure 897010DEST_PATH_IMAGE044
，其中，
Figure 541618DEST_PATH_IMAGE045
表示机载记 录时间序列数据在区间[1，2]的最大值，即
Figure 138952DEST_PATH_IMAGE046
。
Figure 894419DEST_PATH_IMAGE047
表示机载记录时间序列数据在区间[3，4]的最大值，即
Figure 955916DEST_PATH_IMAGE048
。机 载记录时间序列数据在区间[1，4]的最小值的求解过程，与上述过程类似，在此不做赘述。
通过状态转移方程，可以自动快捷地任意子区间的最大值和最小值，从而获得最大值稀疏表和最小值稀疏表。
在一个可选的实施例中，请参阅图7，所述步骤S60，还包括步骤S61~S62，具体如下：
S61：获取每个所述子序列数据的序列长度、起始序列编号以及末尾序列编号；
S62：根据所述序列长度、所述起始序列编号以及所述末尾序列编号，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值；其中，查询获得每个所述子序列数据的最大值和最小值的公式为：
Figure 896059DEST_PATH_IMAGE049
Figure 574165DEST_PATH_IMAGE050
其中，
Figure 274267DEST_PATH_IMAGE051
表示起始序列编号为
Figure 987009DEST_PATH_IMAGE024
，末尾序列编号为
Figure 176681DEST_PATH_IMAGE025
且序列长度为
Figure 732297DEST_PATH_IMAGE026
的子序列数据的最大值，
Figure 95145DEST_PATH_IMAGE052
表示起始序列编号为
Figure 537759DEST_PATH_IMAGE024
，末尾序列编号为
Figure 695070DEST_PATH_IMAGE025
且序列 长度为
Figure 551031DEST_PATH_IMAGE026
的子序列数据的最小值，
Figure 779887DEST_PATH_IMAGE053
表示向下取整运算。
在本申请实施例中，以
Figure 670483DEST_PATH_IMAGE024
=1，
Figure 998696DEST_PATH_IMAGE025
=7为例进行说明，要查询子序列数据
Figure 531833DEST_PATH_IMAGE054
的最大值，即：
Figure 970904DEST_PATH_IMAGE055
其中，
Figure 388110DEST_PATH_IMAGE056
，
Figure 887225DEST_PATH_IMAGE057
。由于最 大值稀疏表包括这些最大值，从而可以直接从最大值稀疏表中查询获得子序列数据的最大 值。子序列数据
Figure 639149DEST_PATH_IMAGE054
的最小值的查询过程，与上述过程类似，在此 不再赘述。
通过最大值稀疏表和最小值稀疏表，可以自动快捷地查询子序列数据的最大值和最小值，以便后续计算子序列数据的极差。
在一个可选的实施例中，请参阅图8，所述步骤S70，还包括步骤S71~S72，具体如下：
S71：若每个所述子序列数据的极差均小于所述预设的第二阈值，则确定飞机巡航进入稳态；
S72：若存在至少一个所述子序列数据的极差大于或等于所述预设的第二阈值，则确定飞机巡航进入非稳态。
在本申请实施例中，将每个所述子序列数据的极差与预设的第二阈值进行比较，可以直接确定每个所述子序列数据对应的时间窗口内飞机巡航是否进入稳态。具体地，若某个子序列数据的极差小于预设的第二阈值，则该子序列数据对应的时间窗口内飞机巡航进入稳态。
进一步地，可以确定飞机的机载记录时间序列数据对应的时间窗口内飞机是否进入巡航稳态。例如，将飞机的机载记录时间序列数据划分为两个子序列数据，若两个子序列数据对应的时间窗口内飞机均进入巡航稳态，则飞机的机载记录时间序列数据对应的时间窗口内飞机进入巡航稳态。若两个子序列数据对应的时间窗口内存在至少一个时间窗口内飞机进入巡航非稳态，则飞机的机载记录时间序列数据对应的时间窗口内飞机进入巡航非稳态。
通过将每个子序列数据的极差与预设的第二阈值进行比较，可以自动快捷地确定飞机的机载记录时间序列数据对应的时间窗口内飞机是否进入巡航稳态。
实施例2
下述为本申请装置实施例，可以用于执行本申请实施例1中方法的内容。对于本申请装置实施例中未披露的细节，请参照本申请实施例1中方法的内容。
请参见图9，其示出了本申请实施例提供的飞机巡航稳态的确定装置的结构示意图。本申请实施例提供的飞机巡航稳态的确定8，包括：
序列数据获取模块81，用于获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；
序列数据划分模块82，用于将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；
样本标准差获得模块83，用于根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；
巡航稳态确定模块84，用于将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
可选的，所述序列数据获取模块，包括：
第一遍历单元，用于遍历所述机载记录时间序列数据的每个机载记录数据，计算第一个机载记录数据到当前机载记录数据的序列前缀和以及序列前缀平方和，获得第一序列前缀和以及第一序列前缀平方和；
第二遍历单元，用于遍历所述当前机载记录数据的下一个机载记录数据，将所述下一个机载记录数据与所述第一序列前缀和求和，获得第二序列前缀和；将所述下一个机载记录数据的平方与所述第一序列前缀平方和求和，获得第二序列前缀平方和；
第三遍历单元，用于继续遍历下一个机载记录数据，直至计算所述第一个机载记录数据到最后一个机载记录数据的序列前缀和以及序列前缀平方和。
可选的，所述样本标准差获得模块，包括：
第一序列长度获取单元，用于获取每个所述子序列数据的序列长度；
标准差获得单元，用于根据所述序列长度、所述子序列前缀和以及所述子序列前缀平方和，获得每个所述子序列数据的样本标准差。
可选的，所述巡航稳态确定模块，包括：
第一判断单元，用于若每个所述子序列数据的样本标准差均小于所述预设的第一阈值，则确定飞机巡航进入稳态；
第二判断单元，用于若存在至少一个子序列数据的样本标准差大于或等于所述预设的第一阈值，则确定飞机巡航进入非稳态。
可选的，所述飞机巡航稳态的确定装置，包括：
预处理模块，用于对所述机载记录时间序列数据进行预处理，获得最大值稀疏表和最小值稀疏表；
最值查询模块，用于根据区间最值查询方法，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值；
极差计算模块，用于计算每个所述子序列数据的最大值和最小值的极差，将所述极差与预设的第二阈值比较，获得第二比较结果，根据所述第二比较结果确定飞机巡航稳态。
可选的，预处理模块，包括：
最值确定单元，用于根据状态转移方程，确定所述机载记录时间序列数据在预设区间内任意子区间的最大值和最小值；
稀疏表获得单元，用于根据所有子区间的最大值和最小值，获得最大值稀疏表和最小值稀疏表。
可选的，最值查询模块，包括：
第二序列长度获取单元，用于获取每个所述子序列数据的序列长度、起始序列编号以及末尾序列编号；
最值查询单元，用于根据所述序列长度、所述起始序列编号以及所述末尾序列编号，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值。
可选的，极差计算模块，包括：
第三判断单元，用于若每个所述子序列数据的极差均小于所述预设的第二阈值，则确定飞机巡航进入稳态；
第四判断单元，用于若存在至少一个所述子序列数据的极差大于或等于所述预设的第二阈值，则确定飞机巡航进入非稳态。
应用本申请实施例，通过获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态，从而提高了计算效率，节约了计算算力。
实施例3
下述为本申请设备实施例，可以用于执行本申请实施例1中方法的内容。对于本申请设备实施例中未披露的细节，请参照本申请实施例1中方法的内容。
请参阅图10，本申请还提供一种电子设备300，电子设备可以具体为计算机、手机、平板电脑、交互平板等，在本申请的示例性实施例中，电子设备300为交互平板，交互平板可以包括：至少一个处理器301、至少一个存储器302，至少一个显示器，至少一个网络接口303，用户接口304以及至少一个通信总线305。
其中，用户接口304主要用于为用户提供输入的接口，获取用户输入的数据。可选的，用户接口还可以包括标准的有线接口、无线接口。
其中，网络接口303可选的可以包括标准的有线接口、无线接口（如WI-FI接口）。
其中，通信总线305用于实现这些组件之间的连接通信。
其中，处理器301可以包括一个或者多个处理核心。处理器利用各种接口和线路连接整个电子设备内的各个部分，通过运行或执行存储在存储器内的指令、程序、代码集或指令集，以及调用存储在存储器内的数据，执行电子设备的各种功能和处理数据。可选的，处理器可以采用数字信号处理（Digital Signal Processing，DSP）、现场可编程门阵列（Field-Programmable Gate Array，FPGA）、可编程逻辑阵列（Programmable Logic Array，PLA）中的至少一种硬件形式来实现。处理器可集成中央处理器（Central ProcessingUnit，CPU）、图像处理器（Graphics Processing Unit，GPU）和调制解调器等中的一种或几种的组合。其中，CPU主要处理操作系统、用户界面和应用程序等；GPU用于负责显示层所需要显示的内容的渲染和绘制；调制解调器用于处理无线通信。可以理解的是，上述调制解调器也可以不集成到处理器中，单独通过一块芯片进行实现。
其中，存储器302可以包括随机存储器（Random Access Memory，RAM），也可以包括只读存储器（Read-Only Memory）。可选的，该存储器包括非瞬时性计算机可读介质（non-transitory computer-readable storage medium）。存储器可用于存储指令、程序、代码、代码集或指令集。存储器可包括存储程序区和存储数据区，其中，存储程序区可存储用于实现操作系统的指令、用于至少一个功能的指令（比如触控功能、声音播放功能、图像播放功能等）、用于实现上述各个方法实施例的指令等；存储数据区可存储上面各个方法实施例中涉及到的数据等。存储器可选的还可以是至少一个位于远离前述处理器的存储装置。作为一种计算机存储介质的存储器中可以包括操作系统、网络通信模块、用户接口模块、操作应用程序。
处理器可以用于调用存储器中存储的视频分辨率调整方法的应用程序，并具体执行上述所示实施例1的方法步骤，具体执行过程可以参见实施例1所示的具体说明，在此不进行赘述。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (9)

1.一种飞机巡航稳态的确定方法，其特征在于，所述方法包括如下步骤：
获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；
将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；
根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；其中，获取每个所述子序列数据的序列长度；根据所述序列长度、所述子序列前缀和以及所述子序列前缀平方和，获得每个所述子序列数据的样本标准差；其中，计算每个所述子序列数据的样本标准差的公式为：
Figure 337032DEST_PATH_IMAGE001
Figure 785331DEST_PATH_IMAGE002
其中，
Figure 237172DEST_PATH_IMAGE003
表示第
Figure 38906DEST_PATH_IMAGE004
个元素到第
Figure 384437DEST_PATH_IMAGE005
个元素且序列长度为
Figure 879003DEST_PATH_IMAGE006
的子序列数据的样 本标准差，
Figure 677195DEST_PATH_IMAGE007
表示第1个元素到第
Figure 781155DEST_PATH_IMAGE008
个元素的子序列前缀和，
Figure 981192DEST_PATH_IMAGE009
表示第1个元 素到第
Figure 381081DEST_PATH_IMAGE008
个元素的子序列前缀平方和，
Figure 807514DEST_PATH_IMAGE010
表示第1个元素到第
Figure 75684DEST_PATH_IMAGE005
个元素的子序列前缀 和，
Figure 5594DEST_PATH_IMAGE011
表示第1个元素到第
Figure 966597DEST_PATH_IMAGE005
个元素的子序列前缀平方和；
将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
2.根据权利要求1所述的飞机巡航稳态的确定方法，其特征在于：
所述计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和的步骤，包括：
遍历所述机载记录时间序列数据的每个机载记录数据，计算第一个机载记录数据到当前机载记录数据的序列前缀和以及序列前缀平方和，获得第一序列前缀和以及第一序列前缀平方和；
遍历所述当前机载记录数据的下一个机载记录数据，将所述下一个机载记录数据与所述第一序列前缀和求和，获得第二序列前缀和；将所述下一个机载记录数据的平方与所述第一序列前缀平方和求和，获得第二序列前缀平方和；
继续遍历下一个机载记录数据，直至计算所述第一个机载记录数据到最后一个机载记录数据的序列前缀和以及序列前缀平方和。
3.根据权利要求1所述的飞机巡航稳态的确定方法，其特征在于：
所述将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态的步骤，包括：
若每个所述子序列数据的样本标准差均小于所述预设的第一阈值，则确定飞机巡航进入稳态；
若存在至少一个子序列数据的样本标准差大于或等于所述预设的第一阈值，则确定飞机巡航进入非稳态。
4.根据权利要求1所述的飞机巡航稳态的确定方法，其特征在于，还包括：
对所述机载记录时间序列数据进行预处理，获得最大值稀疏表和最小值稀疏表；
根据区间最值查询方法，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值；
计算每个所述子序列数据的最大值和最小值的极差，将所述极差与预设的第二阈值比较，获得第二比较结果，根据所述第二比较结果确定飞机巡航稳态。
5.根据权利要求4所述的飞机巡航稳态的确定方法，其特征在于：
所述对所述机载记录时间序列数据进行预处理，获得最大值稀疏表和最小值稀疏表的步骤，包括：
根据状态转移方程，确定所述机载记录时间序列数据在预设区间内任意子区间的最大值和最小值；
根据所有子区间的最大值和最小值，获得最大值稀疏表和最小值稀疏表；其中，所述状态转移方程的公式为：
Figure 113282DEST_PATH_IMAGE012
Figure 185143DEST_PATH_IMAGE013
其中，
Figure 235139DEST_PATH_IMAGE014
表示所述机载记录时间序列数据在区间
Figure 101464DEST_PATH_IMAGE015
的最大值，
Figure 236910DEST_PATH_IMAGE016
表示所述机载记录时间序列数据在区间
Figure 846883DEST_PATH_IMAGE015
的最小值。
6.根据权利要求4所述的飞机巡航稳态的确定方法，其特征在于：
所述根据区间最值查询方法，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值的步骤，包括：
获取每个所述子序列数据的序列长度、起始序列编号以及末尾序列编号；
根据所述序列长度、所述起始序列编号以及所述末尾序列编号，从所述最大值稀疏表和所述最小值稀疏表中查询获得每个所述子序列数据的最大值和最小值；其中，查询获得每个所述子序列数据的最大值和最小值的公式为：
Figure 751385DEST_PATH_IMAGE017
Figure 54191DEST_PATH_IMAGE018
其中，
Figure 169609DEST_PATH_IMAGE019
表示起始序列编号为
Figure 458639DEST_PATH_IMAGE004
，末尾序列编号为
Figure 342281DEST_PATH_IMAGE005
且序列长度为
Figure 691354DEST_PATH_IMAGE006
的子 序列数据的最大值，
Figure 660447DEST_PATH_IMAGE020
表示起始序列编号为
Figure 753168DEST_PATH_IMAGE004
，末尾序列编号为
Figure 756896DEST_PATH_IMAGE005
且序列长度 为
Figure 775405DEST_PATH_IMAGE006
的子序列数据的最小值，
Figure 231794DEST_PATH_IMAGE021
表示向下取整运算。
7.根据权利要求4所述的飞机巡航稳态的确定方法，其特征在于：
所述计算每个所述子序列数据的最大值和最小值的极差，将所述极差与预设的第二阈值比较，获得第二比较结果，根据所述第二比较结果确定飞机巡航稳态的步骤，包括：
若每个所述子序列数据的极差均小于所述预设的第二阈值，则确定飞机巡航进入稳态；
若存在至少一个所述子序列数据的极差大于或等于所述预设的第二阈值，则确定飞机巡航进入非稳态。
8.一种飞机巡航稳态的确定装置，其特征在于，包括：
序列数据获取模块，用于获取飞机的机载记录时间序列数据，计算所述机载记录时间序列数据的序列前缀和以及序列前缀平方和；
序列数据划分模块，用于将所述机载记录时间序列数据划分为若干个子序列数据；根据所述机载记录时间序列数据的序列前缀和以及序列前缀平方和，获得每个所述子序列数据的子序列前缀和以及子序列前缀平方和；
样本标准差获得模块，用于根据每个所述子序列数据的子序列前缀和以及子序列前缀平方和，获得每个所述子序列数据的样本标准差；其中，获取每个所述子序列数据的序列长度；根据所述序列长度、所述子序列前缀和以及所述子序列前缀平方和，获得每个所述子序列数据的样本标准差；其中，计算每个所述子序列数据的样本标准差的公式为：
Figure 862627DEST_PATH_IMAGE001
Figure 455282DEST_PATH_IMAGE002
其中，
Figure 146158DEST_PATH_IMAGE003
表示第
Figure 89843DEST_PATH_IMAGE004
个元素到第
Figure 524367DEST_PATH_IMAGE005
个元素且序列长度为
Figure 237108DEST_PATH_IMAGE006
的子序列数据的样 本标准差，
Figure 331840DEST_PATH_IMAGE007
表示第1个元素到第
Figure 762822DEST_PATH_IMAGE008
个元素的子序列前缀和，
Figure 735457DEST_PATH_IMAGE009
表示第1个元 素到第
Figure 568284DEST_PATH_IMAGE008
个元素的子序列前缀平方和，
Figure 335383DEST_PATH_IMAGE010
表示第1个元素到第
Figure 394605DEST_PATH_IMAGE005
个元素的子序列前缀 和，
Figure 295565DEST_PATH_IMAGE011
表示第1个元素到第
Figure 91221DEST_PATH_IMAGE005
个元素的子序列前缀平方和；
巡航稳态确定模块，用于将每个所述子序列数据的样本标准差与预设的第一阈值进行比较，获得第一比较结果，根据所述第一比较结果确定飞机巡航稳态。
9.一种电子设备，其特征在于，包括：处理器和存储器；其中，所述存储器存储有计算机程序，所述计算机程序适于由所述处理器加载并执行如权利要求1至7中任意一项所述的飞机巡航稳态的确定方法。
CN202211429957.9A 2022-11-16 2022-11-16 飞机巡航稳态的确定方法、装置以及计算机设备 Active CN115512461B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211429957.9A CN115512461B (zh) 2022-11-16 2022-11-16 飞机巡航稳态的确定方法、装置以及计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211429957.9A CN115512461B (zh) 2022-11-16 2022-11-16 飞机巡航稳态的确定方法、装置以及计算机设备

Publications (2)

Publication Number Publication Date
CN115512461A CN115512461A (zh) 2022-12-23
CN115512461B true CN115512461B (zh) 2023-02-07

Family

ID=84514310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211429957.9A Active CN115512461B (zh) 2022-11-16 2022-11-16 飞机巡航稳态的确定方法、装置以及计算机设备

Country Status (1)

Country Link
CN (1) CN115512461B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1479744A (en) * 1973-11-16 1977-07-13 Tvi Television Ind Sa Device for the automatic reading and evaluation of recorded curves
WO2007115140A2 (en) * 2006-03-31 2007-10-11 Alaka'i Technologies Aircraft-engine trend monitoring methods and systems
EP2709328A2 (en) * 2012-09-14 2014-03-19 Fujitsu Limited Coefficient determining apparatus, equalizer, receiver and transmitter
CN103778104A (zh) * 2012-10-22 2014-05-07 富士通株式会社 信息处理装置、信息处理方法以及电子设备
CA2953251A1 (fr) * 2014-06-25 2015-12-30 Safran Aircraft Engines Procede de surveillance d'une degradation d'un dispositif embarque d'un aeronef incluant la determination d'un seuil de comptage
CN107767612A (zh) * 2016-08-22 2018-03-06 波音公司 用于飞行器发动机的热事件指示器
CN114641061A (zh) * 2022-03-08 2022-06-17 重庆邮电大学 一种空对地随机接入级联长前导序列检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1479744A (en) * 1973-11-16 1977-07-13 Tvi Television Ind Sa Device for the automatic reading and evaluation of recorded curves
WO2007115140A2 (en) * 2006-03-31 2007-10-11 Alaka'i Technologies Aircraft-engine trend monitoring methods and systems
EP2709328A2 (en) * 2012-09-14 2014-03-19 Fujitsu Limited Coefficient determining apparatus, equalizer, receiver and transmitter
CN103778104A (zh) * 2012-10-22 2014-05-07 富士通株式会社 信息处理装置、信息处理方法以及电子设备
CA2953251A1 (fr) * 2014-06-25 2015-12-30 Safran Aircraft Engines Procede de surveillance d'une degradation d'un dispositif embarque d'un aeronef incluant la determination d'un seuil de comptage
CN107767612A (zh) * 2016-08-22 2018-03-06 波音公司 用于飞行器发动机的热事件指示器
CN114641061A (zh) * 2022-03-08 2022-06-17 重庆邮电大学 一种空对地随机接入级联长前导序列检测方法

Also Published As

Publication number Publication date
CN115512461A (zh) 2022-12-23

Similar Documents

Publication Publication Date Title
CN109669663B (zh) 音区幅值获取方法、装置、电子设备及存储介质
US20160125053A1 (en) System and method for storage and analysis of time-based data
CN106326928B (zh) 一种目标识别方法及设备
CN105608113A (zh) 判断文本中poi数据的方法及装置
JP2022003388A (ja) 車載機器の応答速度を試験する方法、装置、機器及び記憶媒体
CN111079428B (zh) 一种分词和行业词典构建方法、装置以及可读存储介质
US6480770B1 (en) Par system for analyzing aircraft flight data
WO2024192998A1 (zh) 确定目标位置的方法、装置、系统及存储介质
CN113962420A (zh) 一种飞行器燃油消耗的预测方法、计算设备及存储介质
CN111540344B (zh) 声学网络模型训练方法、装置及电子设备
CN115512461B (zh) 飞机巡航稳态的确定方法、装置以及计算机设备
US20210390802A1 (en) Method, Computer Program And Device For Processing Signals
CN114117255A (zh) 车辆道路滑行阻力的测试方法、设备及存储介质
CN106970998B (zh) 新闻数据的更新方法及装置
CN115372988B (zh) 一种飞机尾涡的识别定位方法、装置及介质
CN111221874A (zh) 盘旋航段的识别方法、装置、电子设备和存储介质
CN111240361A (zh) 一种基于四维数据的飞行控制方法、装置、设备及存储介质
CN104615008B (zh) 一种自动紧急刹车aeb系统的测试结果分析方法和系统
CN113554072B (zh) 一种飞行动作划分方法、计算机可读存储介质及终端设备
CN114861321A (zh) 交通流仿真的问题场景提取方法、装置、设备及介质
CN111627422B (zh) 语音加速检测方法、装置、设备及可读存储介质
CN111159172B (zh) 一种数据处理方法、装置及电子设备
CN112308285A (zh) 信息处理方法、装置、电子设备和计算机可读存储介质
CN117831352B (zh) 飞机起飞机场预测方法、系统、电子设备及存储介质
CN111104374A (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