运载火箭弹性频率在线辨识的方法及存储介质与流程

1.本发明涉及航空航天技术领域，尤其涉及一种运载火箭弹性频率在线辨识的方法及存储介质。


背景技术：

2.近年来，全球商业航天产业规模持续增长，已成为带动全球航天经济发展的重要引擎，低成本、高可靠性的商业运载火箭研制和发射更是商业航天的核心基础。为提高运载能力，运载火箭要尽量提高发动机质量比、减小结构质量，箭体的结构刚度减小，会表现出较为明显的弹性特征。
3.由于弹性作用，飞行中的箭体在外力作用下可能会产生显著的弯曲形变，弹性的弯曲变形将造成附加角速率和角偏差。而惯测组合敏感装置(速率陀螺和加速度计)会敏感检测到这些弯曲形变，并会将弹性信息反馈至姿态控制回路，从而影响姿态控制系统工作，严重时甚至直接导致姿态控制发散，引起箭体的损毁。为避免弹性振动信息的影响，侧喷流姿态控制回路设计时，会根据箭体的弹性固有频率(弹性频率)设计对应的滤波器，实现对特定弹性频率噪声衰减和抑制，降低箭体弹性对姿态控制的影响。
4.传统技术中，箭体弹性频率都是依赖环境专业数学建模计算和模态试验结合来确定，但这种方式由于数学模型偏差以及地面模态试验环境与箭体在实际飞行时的环境的不一致性，环境专业计算得到的弹性频率会存在较大的偏差，即是现有技术获得的箭体弹性频率的精度较低，不利于箭体的稳定飞行控制。


技术实现要素：

5.本技术实施例通过提供一种运载火箭弹性频率在线辨识的方法，解决了现有技术中弹性频率计算的精度较低等技术问题，提升了弹性频率辨识的精度，进而有利于控制箭体的稳定飞行。
6.一方面，本技术通过本技术的一实施例提供一种运载火箭弹性频率在线辨识的方法，所述方法包括如下步骤：
7.s1、根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；
8.s2、根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；
9.s3、根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率；
10.s4、对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率；
11.s5、从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角
速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
12.可选地，所述方法还包括：
13.当所述运载火箭进入下一个飞行时段时，获取并更新所述运载火箭的当前飞行时段，重复执行所述s2～s5的步骤。
14.可选地，所述s1具体包括：
15.根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1；
16.基于所述带阻滤波器的倒数同频谐振放大原理，计算每个所述带阻滤波器的倒数环节所对应的滤波频率。
17.可选地，所述带阻滤波器的倒数环节个数n1为所述每个所述带阻滤波器的倒数环节所对应的滤波频率为f(i)；
18.其中，所述弹性频率范围为(f0～f1)，所述辨识频率间隔为δf，f(i)∈[f0,f1]，i为不超过n1的正整数。
[0019]
可选地，所述飞行运行拍数n2为：
[0020][0021]
其中，所述当前飞行时段为(t0～t1)，所述飞行控制周期为t。
[0022]
可选地，所述s3具体包括：
[0023]
采集所述运载火箭在所述当前飞行时段中n2个所述飞行运行拍数所对应时刻的俯仰角速率和偏航角速率；
[0024]
分别将所述n2个不同时刻的俯仰角速率和偏航角速率输入n1个所述带阻滤波器的倒数环节中，计算得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵的大小均为n1×
n2。
[0025]
可选地，所述s4具体包括：
[0026]
对所述俯仰方向和偏航方向角速率矩阵中对应的元素进行矢量合成，然后进行行向量绝对值求和处理，得到求和角速率列向量，所述求和角速率列向量中包括n1个求和角速率；
[0027]
从所述求和角速率列向量中的n1个求和角速率中，选取数值最大的目标角速率，并确定所述目标角速率对应的倒数环节序号。
[0028]
可选地，所述求和角速率列向量为：
[0029][0030]
其中，δθ
f
(i)为所述求和角速率列向量，为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的仰俯角速率，δψ
f
(i,j)为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的偏航角速率，i为不超过n1的正整数，j为
不超过n2的正整数。
[0031]
可选地，若所述弹性频率范围为(f0～f1)、所述目标角速率对应的倒数环节序号为i
flag
，则所述目标辨识频率为f(i
flag
)；其中，f(i
flag
)∈(f0～f1)。
[0032]
另一方面，本技术通过本技术的一实施例提供一种运载火箭弹性频率的辨识装置，所述装置包括计算模块、处理模块及频率辨识模块，其中：
[0033]
所述计算模块，用于根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；
[0034]
所述计算模块，还用于根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；
[0035]
所述处理模块，用于根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率；
[0036]
所述处理模块，还用于对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率；
[0037]
所述频率辨识模块，用于从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
[0038]
关于本技术未阐述的内容可对应参见前述方法实施例中的相关介绍，这里不再赘述。
[0039]
另一方面，本技术通过本技术的一实施例提供一种终端设备，所述终端设备包括处理器、存储器、通信接口和总线；所述处理器、所述存储器和所述通信接口通过所述总线连接并完成相互间的通信；所述存储器存储可执行程序代码；所述处理器通过读取所述存储器中存储的可执行程序代码来运行与所述可执行程序代码对应的程序，以用于如上所提供的一种运载火箭弹性频率在线辨识的方法。
[0040]
另一方面，本技术通过本技术的一实施例提供一种计算机可读存储介质，所述计算机可读存储介质包括计算机指令，当所述计算机指令在终端设备上运行时，使得所述终端设备执行如上所述的运载火箭弹性频率在线辨识的方法。
[0041]
本技术实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：本技术通过根据预设的弹性频率范围及辨识频率间隔计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；根据所述运载火箭的当前飞行时段和预设的飞行控制周期计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵；进一步对所述角速率矩阵进行筛选处理得到目标角速率对应的倒数环节序号；最后从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。这样既能够解决现有技术中弹性频率计算的精度较低等技术问题，又能提升弹性频率辨识的精度，进而有利于控制箭体的稳定飞行。
附图说明
[0042]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0043]
图1是本技术实施例提供的一种带阻滤波器倒数环节构成的单输入
‑
单输出系统的结构示意图。
[0044]
图2是本技术实施例提供的一种带阻滤波器倒数环节g(s)波特图的示意图。
[0045]
图3是本技术实施例提供的一种运载火箭弹性频率在线辨识的方法的流程示意图。
[0046]
图4是本技术实施例提供的一种运载火箭弹性频率的辨识装置的结构示意图。
[0047]
图5是本技术实施例提供的一种终端设备的结构示意图。
具体实施方式
[0048]
申请人在提出本技术的过程中还发现：传统技术中箭体弹性(固有)频率都是依赖环境专业数学建模计算和模态试验结合来确定，这种方式存在以下两大不足之处：
[0049]
一是，由于数学模型偏差以及地面模态试验环境与箭体在实际飞行时的环境的不一致性，环境专业计算得到的弹性频率会有较大偏差，降低了姿态控制的回路裕度，不利于运载火箭的稳定飞行。
[0050]
二是，商业运载火箭要求低成本、研制周期短，箭体振动模态试验耗时长、费用昂贵，卫星载荷不一样时，同一型号每一发箭都开展模态试验，无法满足商业运载火箭快速研制的需求。因此，低成本、高可靠性的弹性频率辨识方法具有重要的工程应用价值。
[0051]
目前，一般地箭体弹性频率辨识方法(离线辨识或在线辨识)均采用傅里叶变换的方式，即对箭体角速度信号进行傅里叶变换，识别出对应的频率点。但这种方式算法比较复杂、计算量较大、对计算机性能要求高、可靠性低，一旦辨识错误会导致飞行姿态失稳。
[0052]
为解决上述问题，本技术实施例通过提供一种运载火箭弹性频率在线辨识的方法。本技术实施例的技术方案为解决上述技术问题，总体思路如下：s1、根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；s2、根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；s3、根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率；
[0053]
s4、对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率；s5、从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
[0054]
为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
[0055]
首先说明，本文中出现的术语“和/或”，仅仅是一种描述关联对象的关联关系，表
示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
[0056]
本技术旨在提供一种运载火箭弹性频率在线辨识的方法、装置、终端设备及计算机可读存储介质。本技术利用带阻滤波器的倒数同频谐振放大的原理，解决运载火箭飞行过程中侧喷流控制段弹性频率在线辨识的应用需求。
[0057]
首先，介绍本技术涉及的带阻滤波器倒数同频谐振放大原理。请参见图1是本技术实施例示出的一种带阻滤波器单倒数环节构成的输入
‑
单输出系统的结构示意图。如图1中，记带阻滤波器倒数环节g(s)的函数表达式为如下公式(1)所示：
[0058][0059]
其中，ωn为谐振频率对应的角速率，ζ1、ζ2为带阻滤波器自身的阻尼系数。ω
n
>0,0<ζ2<ζ1<1。
[0060]
需要说明的是，带阻滤波器倒数环节是一个单输入
‑
单输出系统，不同频率信号经过该系统后会出现幅值和相位的改变。当输入信号频率等于系统谐振频率ωn时，在对应频率点处会出现最大谐振峰值即在该频率点处对输入信号具有放大作用，放大的谐振峰值由阻尼系数ζ1、ζ2确定。因此，本技术可根据幅值最大的输出信号就可以反向确定输入信号频率，实现频率辨识。
[0061]
请参见图2示出一种带阻滤波器倒数环节g(s)的波特图。如图2中，以ω
n
＝120rad/s,k＝ζ1/ζ2＝60为例，当与ω
n
＝120rad/s同频的角速率信号进入该系统后，幅值被放大到35.6db。系统根据被放大的信号35.6db即可确定输入信号的频率ω
n
＝120rad/s，从而实现频率辨识。
[0062]
接着，介绍本技术涉及的方法实施例。请参见图3，是本技术实施例提供的一种运载火箭弹性频率在线辨识的方法的流程示意图。如图3所示的方法包括如下实施步骤：
[0063]
s1、根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率。
[0064]
本技术可根据环境专业数学模型计算结果和箭体地面模态首次试验结果，初始化确定带阻滤波器的倒数环节个数及其谐振频率点(也可称为滤波频率)。环境专业会根据总体输入(包含弹道、气动数据等)，通过数学建模仿真和模态试验给出带较大偏差的弹性频率。本技术考虑到箭体弹性频率及偏差设置辨识的弹性频率范围为(f0～f1)hz，辨识频率间隔为δf hz。
[0065]
相应地，本技术根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1，其具体如下公式(2)所示：
[0066][0067]
进一步本技术还可根据带阻滤波器倒数同频谐振放大原理，计算每个带阻滤波器倒数环节所对应的频率点ωn或滤波频率f(i)。其中，ω
n
(i)＝2πf(i),f(i)∈[f0,f1],i＝1～n1。即i为不超过n1的正整数。可选地，本技术可记k＝ζ1/ζ2，根据辨识区别分段k一般取值为30～100左右。
[0068]
举例来说，本技术设计弹性频率范围为(f0～f1)，其首尾偏差为4hz，以0.01hz为辨识频率间隔。则设计的带阻滤波器倒数环节个数为对应振荡环节参数设置为ω
n
(i)＝2πf(i),k＝60,i＝1～41，其中f(i)∈[f0,f1]。
[0069]
s2、根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2。
[0070]
本技术可初始化时间段，即初始化当前飞行时段。本技术记运载火箭的飞行控制周期为t，运载火箭飞行过程中，短时间内弹性频率的变化较小，可近似认为不变。初始化当前飞行时段为(t0～t1)，则这段飞行时段内运载火箭的飞行运行拍数为n2，其具体如下公式(3)所示：
[0071][0072]
举例来说，本技术取当前飞行时间段(t0～t1)为1s，飞行控制周期t＝0.01s，则其计算的飞行运行拍数
[0073]
s3、根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率。
[0074]
本技术可通过速率陀螺仪采集按照飞行控制周期t在当前飞行时段(t0～t1)内不同时刻的俯仰角速率和偏航角速率δψ，即采集n2个所述飞行运行拍数所对应的俯仰角速率和偏航角速率δψ
f
(i,j)。i＝1～n1,j＝1～n2。进一步可分别将这n2个不同时刻的俯仰角速率和偏航角速率输入n1个所述带阻滤波器的倒数环节中，计算得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵的大小均为n1×
n2。
[0075]
具体地在t0时刻，分别将速率陀螺仪采集输出的俯仰角速率和偏航角速率(信号)输入到已经预置的n1个带阻滤波器倒数环节中，得到对应计算的角速率δψ
f
(i,1)，i＝1～n1。下一个周期t中，继续将速率陀螺仪输出的俯仰角速率和偏航角速率(信
号)输入至n1个带阻滤波器倒数环节中，得到对应的角速率δψ
f
(i,2)，i＝1～n1；依次类推，
……
，直至t1时刻得到δψ
f
(i,n2),i＝1～n1。
[0076]
举例来说，引用步骤s3例子在1s内，输入的俯仰和偏航角速率通过带阻滤波器的振荡环节放大得到的对应角速率矩阵为：δψ
f
(i,j),i＝1～41,j＝1～101。
[0077]
s4、对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率。
[0078]
本技术可对俯仰方向和偏航方向角速率矩阵中对应的元素进行矢量合成，然后进行行向量绝对值求和处理，得到求和角速率列向量，该列向量中包括n1个求和角速率；从列向量的n1个求和角速率中选取数值最大的目标角速率，并确定出该目标角速率对应的倒数环节序号。
[0079]
在一具体实施例中，所述求和角速率列向量为如下公式(4)所示：
[0080][0081]
其中，δθ
f
(i)为所述求和角速率列向量，为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的仰俯角速率，δψ
f
(i,j)为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的偏航角速率，i为不超过n1的正整数，j为不超过n2的正整数。
[0082]
s5、从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
[0083]
具体地本技术根据(t0～t1)当前飞行时段内得到的角速率矩阵δψ
f
(i,j),i＝1～n1,j＝1～n2进行对应的元素矢量合成，然后进行行向量绝对值求和处理，得列向量记(t0～t1)当前飞行时段内列向量δθ
f
(i),i＝1～n1中的最大值为δθ
max
，与之对应的倒数环节序号为i
flag
，i∈[1,n1]。则相应地在辨识频率范围(f0～f1)内对应的目标辨识频率为f(i
flag
)。
[0084]
举例来说，引用前述步骤s1和s2的例子，本技术根据(t0～t1)时间段内得到的角速率矩阵δψ
f
(i,j),i＝1～41,j＝1～101，进行对应的元素矢量合成，然后进行行向量绝对值求和处理，得到列向量记(t0～
t1)时间段内列向量δθ
f
(i),i＝1～41中最大值为δθ
max
，与之对应的序号为i
flag
∈[1,41]，在(f0～f1)内对应的目标辨识频率为f(i
flag
)。
[0085]
在可选实施例中，当运载火箭进入下一个飞行时段时，可获取并更新所述运载火箭的当前飞行时段(t0'～t1')，其时长仍可为1s。重复上述步骤s2～s5，得到当前飞行时段的目标辨识频率f
′
(i
′
flag
)。循环往复，从而达到侧喷流控制段箭体弹性频率在线辨识的目标，同时保证辨识频率在(f0～f1)范围内，不会出现辨识频率偏差过大的情况。
[0086]
本技术提出一种运载火箭弹性频率在线辨识的方法，其具有设计原理简单、算法简洁、计算量小、可靠性高、辨识精度高等优点，同时融合了环境专业提供的首次模态试验数据信息，不会出现很大的辨识偏差，辨识结果可靠性高，适合箭体弹性频率快速在线辨识。在商业运载火箭使用上，可以极大地减少模态试验次数，缩短研制试验周期，减低研制成本，提升火箭竞争能力，具有很高的工程应用价值。
[0087]
通过实施本技术，通过根据预设的弹性频率范围及辨识频率间隔计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；根据所述运载火箭的当前飞行时段和预设的飞行控制周期计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵；进一步对所述角速率矩阵进行筛选处理得到目标角速率对应的倒数环节序号；最后从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。这样既能够解决现有技术中弹性频率计算的精度较低等技术问题，又能提升弹性频率辨识的精度，进而有利于控制箭体的稳定飞行。
[0088]
基于同一发明构思，本技术另一实施例提供一种实施本技术实施例中所述方法对应的装置及终端设备。请参见图4，是本技术实施例提供的一种运载火箭弹性频率的辨识装置的结构示意图。如图4所示的装置400包括：计算模块401、处理模块402及频率辨识模块403，其中：
[0089]
所述计算模块401，用于根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；
[0090]
所述计算模块401，还用于根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；
[0091]
所述处理模块402，用于根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率；
[0092]
所述处理模块402，还用于对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率；
[0093]
所述频率辨识模块403，用于从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
[0094]
可选地，所述装置还包括更新模块404，
[0095]
所述更新模块404，用于当所述运载火箭进入下一个飞行时段时，获取并更新所述
运载火箭的当前飞行时段，重复执行所述计算模块401、处理模块402及频率辨识模块403所执行的步骤。
[0096]
可选地，所述计算模块401具体用于：
[0097]
根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1；
[0098]
基于所述带阻滤波器的倒数同频谐振放大原理，计算每个所述带阻滤波器的倒数环节所对应的滤波频率。
[0099]
可选地，所述带阻滤波器的倒数环节个数n1为所述每个所述带阻滤波器的倒数环节所对应的滤波频率为f(i)；
[0100]
其中，所述弹性频率范围为(f0～f1)，所述辨识频率间隔为δf，f(i)∈[f0,f1]，i为不超过n1的正整数。
[0101]
可选地，所述飞行运行拍数n2为：
[0102][0103]
其中，所述当前飞行时段为(t0～t1)，所述飞行控制周期为t。
[0104]
可选地，所述处理模块402具体用于：
[0105]
采集所述运载火箭在所述当前飞行时段中n2个所述飞行运行拍数所对应时刻的俯仰角速率和偏航角速率；
[0106]
分别将所述n2个不同时刻的俯仰角速率和偏航角速率输入n1个所述带阻滤波器的倒数环节中，计算得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵的大小均为n1×
n2。
[0107]
可选地，所述处理模块402具体用于：
[0108]
对所述俯仰方向和偏航方向角速率矩阵中对应的元素进行矢量合成，然后进行行向量绝对值求和处理，得到求和角速率列向量，所述求和角速率列向量中包括n1个求和角速率；
[0109]
从所述求和角速率列向量中的n1个求和角速率中，选取数值最大的目标角速率，并确定所述目标角速率对应的倒数环节序号。
[0110]
可选地，所述求和角速率列向量为：
[0111][0112]
其中，δθ
f
(i)为所述求和角速率列向量，为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的仰俯角速率，δψ
f
(i,j)为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的偏航角速率，i为不超过n1的正整数，j为不超过n2的正整数。
[0113]
可选地，若所述弹性频率范围为(f0～f1)、所述目标角速率对应的倒数环节序号为
i
flag
，则所述目标辨识频率为f(i
flag
)；其中，f(i
flag
)∈(f0～f1)。
[0114]
请一并参见图5，是本技术实施例提供的一种终端设备的结构示意图。如图5所示的终端设备500包括：至少一个处理器501、通信接口502、用户接口503和存储器504，处理器501、通信接口502、用户接口503和存储器504可通过总线或者其它方式连接，本发明实施例以通过总线505连接为例。其中，
[0115]
处理器501可以是通用处理器，例如中央处理器(central processing unit，cpu)。
[0116]
通信接口502可以为有线接口(例如以太网接口)或无线接口(例如蜂窝网络接口或使用无线局域网接口)，用于与其他终端或网站进行通信。用户接口503具体可为触控面板，包括触摸屏和触控屏，用于检测触控面板上的操作指令，用户接口503也可以是物理按键或者鼠标。用户接口503还可以为显示屏，用于输出、显示图像或数据。
[0117]
存储器504可以包括易失性存储器(volatile memory)，例如随机存取存储器(random access memory，ram)；存储器也可以包括非易失性存储器(non
‑
volatile memory)，例如只读存储器(read
‑
only memory，rom)、快闪存储器(flash memory)、硬盘(hard disk drive，hdd)或固态硬盘(solid
‑
state drive，ssd)；存储器504还可以包括上述种类的存储器的组合。存储器504用于存储一组程序代码，处理器501用于调用存储器504中存储的程序代码，执行如下操作：
[0118]
s1、根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；
[0119]
s2、根据所述运载火箭的当前飞行时段和预设的飞行控制周期，计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；
[0120]
s3、根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率和偏航角速率进行谐振放大处理，得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵均包括n1×
n2个角速率；
[0121]
s4、对所述角速率矩阵进行筛选处理，得到目标角速率对应的倒数环节序号，所述目标角速率为对所述角速率矩阵中的角速率进行处理后的角速率；
[0122]
s5、从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。
[0123]
可选地，所述处理器501还用于当所述运载火箭进入下一个飞行时段时，获取并更新所述运载火箭的当前飞行时段，重复执行所述s2～s5的步骤。
[0124]
可选地，所述s1具体包括：
[0125]
根据预设的弹性频率范围及辨识频率间隔，计算带阻滤波器的倒数环节个数n1；
[0126]
基于所述带阻滤波器的倒数同频谐振放大原理，计算每个所述带阻滤波器的倒数环节所对应的滤波频率。
[0127]
可选地，所述带阻滤波器的倒数环节个数n1为所述每个所述带阻滤波器的倒数环节所对应的滤波频率为f(i)；
[0128]
其中，所述弹性频率范围为(f0～f1)，所述辨识频率间隔为δf，f(i)∈[f0,f1]，i
为不超过n1的正整数。
[0129]
可选地，所述飞行运行拍数n2为：
[0130][0131]
其中，所述当前飞行时段为(t0～t1)，所述飞行控制周期为t。
[0132]
可选地，所述s3具体包括：
[0133]
采集所述运载火箭在所述当前飞行时段中n2个所述飞行运行拍数所对应时刻的俯仰角速率和偏航角速率；
[0134]
分别将所述n2个不同时刻的俯仰角速率和偏航角速率输入n1个所述带阻滤波器的倒数环节中，计算得到对应的角速率矩阵，所述俯仰方向和偏航方向角速率矩阵的大小均为n1×
n2。
[0135]
可选地，所述s4具体包括：
[0136]
对所述俯仰方向和偏航方向角速率矩阵中对应的元素进行矢量合成，然后进行行向量绝对值求和处理，得到求和角速率列向量，所述求和角速率列向量中包括n1个求和角速率；
[0137]
从所述求和角速率列向量中的n1个求和角速率中，选取数值最大的目标角速率，并确定所述目标角速率对应的倒数环节序号。
[0138]
可选地，所述求和角速率列向量为：
[0139][0140]
其中，δθ
f
(i)为所述求和角速率列向量，为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的仰俯角速率，δψ
f
(i,j)为所述角速率矩阵中处于第i行第j列所述飞行运行拍数所对应时刻的偏航角速率，i为不超过n1的正整数，j为不超过n2的正整数。
[0141]
可选地，若所述弹性频率范围为(f0～f1)、所述目标角速率对应的倒数环节序号为i
flag
，则所述目标辨识频率为f(i
flag
)；其中，f(i
flag
)∈(f0～f1)。
[0142]
由于本实施例所介绍的终端子设备为实施本技术实施例中所述方法所采用的终端设备，故而基于本技术实施例中所介绍的方法实施例，本领域所属技术人员能够了解本实施例的终端设备的具体实施方式以及其各种变化形式，所以在此对于该终端设备如何实现本技术实施例中的方法不再详细介绍。只要本领域所属技术人员实施本技术实施例中方法所采用的终端设备，都属于本技术所欲保护的范围。
[0143]
上述本技术实施例中的技术方案，至少具有如下的技术效果或优点：本技术通过根据预设的弹性频率范围及辨识频率间隔计算带阻滤波器的倒数环节个数n1及每个所述带阻滤波器的倒数环节所对应的滤波频率；根据所述运载火箭的当前飞行时段和预设的飞行控制周期计算所述当前飞行时段内所述运载火箭的飞行运行拍数n2；根据所述飞行运行拍数n2和所述倒数环节个数n1，分别对采集的发生弹性弯曲时所述运载火箭的俯仰角速率
和偏航角速率进行谐振放大处理，得到对应的角速率矩阵；进一步对所述角速率矩阵进行筛选处理得到目标角速率对应的倒数环节序号；最后从每个所述带阻滤波器的倒数环节所对应的滤波频率中，查找与所述目标角速率对应的倒数环节序号所对应的目标频率，以作为所述运载火箭的目标辨识频率。这样既能够解决现有技术中弹性频率计算的精度较低等技术问题，又能提升弹性频率辨识的精度，进而有利于控制箭体的稳定飞行。
[0144]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd
‑
rom、光学存储器等)上实施的计算机程序产品的形式。
[0145]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0146]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0147]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0148]
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
[0149]
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。