一种基于数据滚动粒子滤波的燃油油量数据处理方法与流程

1.本发明属于航空燃油油量数据采集测量技术领域，具体地说，涉及一种基于数据滚动粒子滤波的燃油油量数据处理方法。


背景技术：

2.目前基于计算机技术的粒子式燃油测量系统，大大提高了油量测量的精度，却带来了测量稳定性不好的问题。在实际的应用中，输出的油量数据会频繁跳动，甚至输出和实际情况相去甚远的油量数据，这已严重影响到飞行员对飞机工作状态的正确判断。
3.现有技术通常采用硬件滤波器的方法，在测量部件前端对模拟式的燃油油量数据测量信号进行前端处理。硬件滤波器的不足之处是滤波参数不易更改，不能对频率很高或很低的信号进行滤波，且存在阻抗匹配问题。


技术实现要素：

4.本发明针对现有技术的上述缺陷，提出了一种基于数据滚动粒子滤波的燃油油量数据处理方法，通过在信号采集的过程中，设置一个对燃油油量信号进行后处理的滚动粒子软件滤波器，通过软件的滤波器在测量部件前端对模拟式的测量信号进行前端处理。本发明能够对频率很高或很低的信号进行滤波，同时，滤波参数非常容易更改。使用本发明的粒子滤波方法可以消除油面波动误差，对油量信号进行必要的平滑处理，使其更接近于真实的状态，克服飞机机动飞行带来油量测量的误差，提高了滤波的精度。
5.本发明具体实现内容如下：
6.本发明提出了一种基于数据滚动粒子滤波的燃油油量数据处理方法，设置一个对燃油油量信号进行后处理的滚动粒子软件滤波器，在滚动粒子软件滤波器中设置滤波队列地址、滤波队列长度设置单元；然后进行以下操作：
7.步骤1，通过滤波队列地址和滤波队列长度设置单元来设定滤波队列，并定义滤波极限δlimit、滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δlower；
8.步骤2，采用滤波极限δlimit对采样值si进行初次过滤；所述滤波极限δlimit为采样值si的有效范围，为正值；
9.步骤3，采用冒泡法对初次过滤后的采样值si进行排序；对排序后的采样值si，通过滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δlower对所有采样值si进行分类，根据采样值si的值的大小分为大权值粒子中值粒子和小权值粒子，得到大权值粒子集、中值粒子集和小权值粒子集及可信采样值序列vi；
10.步骤4，对于大权值粒子集、中值粒子集和小权值粒子集中的大权值粒子、中值粒子和小权值粒子进行重采样条件判断，对于符合重采样条件的粒子进行重采样，得到新的大权值粒子集、中值粒子集和小权值粒子集；
11.步骤5，将重采样后得到的大权值粒子集、中值粒子集和小权值粒子集中的粒子组成数据队列k；对数据队列k进行队列填充和k次平均；
12.步骤6，将经过k次取平均后的值作为油量输出值进行输出，并作为下一个滤波周期的基准油量值。
13.为了更好地实现本发明，进一步地，在所述步骤3中，对于经过初次过滤的采样值si，还要采用滤波队列上限δupper和滤波队列下限δlower构成的区间[δlower，δupper]进行再次过滤，将不再区间[δlower，δupper]内的采样值si进行剔除，然后再采用冒泡法进行排序。
[0014]
为了更好地实现本发明，进一步地，在所述步骤3中，对所有采样值si进行分类，并计算出可信采样值序列vi的中值δ；进行分类的的具体公式如下：
[0015][0016]
其中：u
i-1表示第i-1个周期的油量输出数据。
[0017]
为了更好地实现本发明，进一步地，得到的所述小权值粒子集为：
[0018][0019]
式中，为小权值粒子，0＜δlower≤si＜δ-为小权值粒子的取值范围。
[0020]
为了更好地实现本发明，进一步地，得到的所述中值粒子集为：
[0021][0022]
式中，为中值粒子，(δ-≤si≤δ

为中值粒子的取值范围。
[0023]
为了更好地实现本发明，进一步地，得到的大权值粒子集为：
[0024][0025]
式中，为大权值粒子，δ

＜si≤δupper为大权值粒子的取值范围。
[0026]
为了更好地实现本发明，进一步地，所述步骤4中，进行重采样条件判断的操作为：
[0027]
对于中值粒子集中的粒子，不进行重采样；
[0028]
对于大权值粒子集中的粒子，进行重采样的判断手段为：先求得滤波队列上限δupper，然后与预先设定的滤波阈值上限δ

进行比较；若在大权值粒子集中的某个粒子的
采样值si取值范围在δ

＜si≤δupper，则不进行重采样；若在大权值粒子集中的某个采样值si取值范围在δupper＜si，则进行重采样；
[0029]
对于小权值粒子集中的粒子，进行重采样的判断手段为：先求得滤波队列下限δlower，然后与预先设定的滤波阈值下限δ-进行比较；若小权值粒子集中的某个粒子的采样值si在δlower≤si＜δ-范围内，则不进行重采样；若若小权值粒子集中的某个粒子的采样值si在δloer＞si内，则进行重采样。
[0030]
为了更好地实现本发明，进一步地，在对大权值粒子集中的粒子进行了重采样后，还要使用线性方程式进行经验降解，将降解后得到的新的大权值粒子作为大权值粒子集中的粒子替换掉对应的旧的大权值粒子；具体经验降解的公式为：
[0031][0032]
为了更好地实现本发明，进一步地，在对小权值粒子集中的粒子进行重采样后，还要使用线性方程式进行经验补偿，将补偿得到新的小权值粒子作为小权值粒子集中的粒子替换掉对应的旧的小权值粒子；具体的经验补偿的公式为：
[0033][0034]
为了更好地实现本发明，进一步地，对大权值粒子集的粒子总数记为na，中值粒子集的粒子总数记为nb，小权值粒子集的粒子总数记为nc；得到数据队列k的值，具体计算为：
[0035]
na nb nc＝k。
[0036]
本发明与现有技术相比具有以下优点及有益效果：
[0037]
(1)本发明采用软件滤波器的方法对测量数据进行后处理，而不是采用硬件滤波器的方法，在测量部件前端对模拟式的测量信号进行前端处理。能够对频率很高或很低的信号进行滤波，同时，滤波参数非常容易更改。用粒子滤波方法消除油面波动误差，对油量信号进行必要的平滑处理，使其更接近于真实的状态，克服飞机机动飞行带来油量测量的误差，提高了滤波的精度。
[0038]
(2)本发明采用软件算法来实现滚动粒子滤波器，粒子滤波器无需硬件，只是一个计算过程，因此可靠性高，不存在阻抗匹配问题。且粒子滤波可以对频率很高或很低的信号进行滤波，这是模拟式硬件滤波器所不及的；用软件算法实现，可以使多个输入通道共用一个软件滤波器，降低了硬件成本；只要适当改变软件滤波器的滤波程序和运算参数，即可改变滤波特性，且便于更改。
附图说明
[0039]
图1为本发明实现的具体流程示意图。
具体实施方式
[0040]
为了更清楚地说明本发明实施例的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，应当理解，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例，因此不应被看作是对保护范围的限定。基于本发明中的实施例，本领域普通技术工作人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0041]
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“设置”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；也可以是直接相连，也可以是通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0042]
实施例1：
[0043]
本实施例提出了一种基于数据滚动粒子滤波的燃油油量数据处理方法，如图1所示，设置一个对燃油油量信号进行后处理的滚动粒子软件滤波器，在滚动粒子软件滤波器中设置滤波队列地址、滤波队列长度设置单元；然后进行以下操作：
[0044]
步骤1，通过滤波队列地址和滤波队列长度设置单元来设定滤波队列，并定义滤波极限δlimit、滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δlower；
[0045]
步骤2，采用滤波极限δlimit对采样值si进行初次过滤；所述滤波极限δlimit为采样值si的有效范围，为正值；
[0046]
步骤3，采用冒泡法对初次过滤后的采样值si进行排序；对排序后的采样值si，通过滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δower对所有采样值si进行分类，根据采样值si的值的大小分为大权值粒子、中值粒子和小权值粒子，得到大权值粒子集、中值粒子集和小权值粒子集及可信采样值序列vi；
[0047]
步骤4，对于大权值粒子集、中值粒子集和小权值粒子集中的大权值粒子、中值粒子和小权值粒子进行重采样条件判断，对于符合重采样条件的粒子进行重采样，得到新的大权值粒子集、中值粒子集和小权值粒子集；
[0048]
步骤5，将重采样后得到的大权值粒子集、中值粒子集和小权值粒子集中的粒子组成数据队列k；对数据队列k进行队列填充和k次平均；
[0049]
步骤6，将经过k次取平均后的值作为油量输出值进行输出，并作为下一个滤波周期的基准油量值。
[0050]
工作原理：在所述粒子软件滤波器中，设置有滤波队列地址1、滤波队列长度设定单元、滤波极限δlimit、滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δlower。通过滤波队列地址定义了采样数据的存放起始位置，方便我们找到滤波的初始值，滤波队列长度设定单元定义了滤波的响应速度和队列长度k，方便我们找到需要滤波的数据个数，滤波极限δlimit定义了采样值的取值范围，用以剔除明显虚假的数据，滤波队列上限δupper和波队列下限δlower定义了滤波序列可信的取值范围，滤波阈值上
限δ

、滤波阈值下限δ-定义了滤波的区域划分，用以进行粒子种类的区分。其工作原理是首先利用滤波极限δlimit将明显的虚假数据去除；用冒泡法进行排序，然后利用滤波队列上限δupper、波队列下限δlower以及滤波阈值上限δ

和滤波阈值下限δ-对采样数据进行分类，得到三类粒子集，对中值粒子集不进行重采样；最后对大权值粒子和小权值粒子组成的粒子集先进行判断，若符合重采样条件，对其使用线性组合方式进行重采样得到新粒子；最后进行队列填充和k次平均。
[0051]
实施例2：
[0052]
本实施例在上述实施例1的基础上，如图1所示，为了更好地实现本发明，进一步地，在所述步骤3中，对于经过初次过滤的采样值si，还要采用滤波队列上限δupper和滤波队列下限δlower构成的区间[δlower，δupper]进行再次过滤，将不再区间[δlower，δupper]内的采样值si进行剔除，然后再采用冒泡法进行排序。
[0053]
进一步地，在所述步骤3中，对所有采样值si进行分类，并计算出可信采样值序列vi的中值δ；进行分类的的具体公式如下：
[0054][0055]
其中：u
i-1表示第i-1个周期的油量输出数据。
[0056]
进一步地，得到的所述小权值粒子集为：
[0057][0058]
式中，为小权值粒子，0＜δlower≤si＜δ-为小权值粒子的取值范围。
[0059]
进一步地，得到的所述中值粒子集为：
[0060][0061]
式中，为中值粒子，(δ-≤si≤δ

为中值粒子的取值范围。
[0062]
进一步地，得到的大权值粒子集为：
[0063][0064]
式中，为大权值粒子，δ

＜si≤δupper为大权值粒子的取值范围。
[0065]
进一步地，所述步骤4中，进行重采样条件判断的操作为：
[0066]
对于中值粒子集中的粒子，不进行重采样；
[0067]
对于大权值粒子集中的粒子，进行重采样的判断手段为：先求得滤波队列上限δupper，然后与预先设定的滤波阈值上限δ

进行比较；若在大权值粒子集中的某个粒子的采样值si取值范围在δ

＜si≤δupper，则不进行重采样；若在大权值粒子集中的某个采样值si取值范围在δupper＜si，则进行重采样；
[0068]
对于小权值粒子集中的粒子，进行重采样的判断手段为：先求得滤波队列下限δlower，然后与预先设定的滤波阈值下限δ-进行比较；若小权值粒子集中的某个粒子的采样值si在δlower≤si＜δ-范围内，则不进行重采样；若若小权值粒子集中的某个粒子的采样值si在δlower＞si内，则进行重采样。
[0069]
进一步地，在对大权值粒子集中的粒子进行了重采样后，还要使用线性方程式进行经验降解，将降解后得到的新的大权值粒子作为大权值粒子集中的粒子替换掉对应的旧的大权值粒子；具体经验降解的公式为：
[0070][0071]
进一步地，在对小权值粒子集中的粒子进行重采样后，还要使用线性方程式进行经验补偿，将补偿得到新的小权值粒子作为小权值粒子集中的粒子替换掉对应的旧的小权值粒子；具体的经验补偿的公式为：
[0072][0073]
进一步地，对大权值粒子集的粒子总数记为na，中值粒子集的粒子总数记为nb，小权值粒子集的粒子总数记为nc；得到数据队列k的值，具体计算为：
[0074]
na nb nc＝k。
[0075]
工作原理：由于飞机的实际油量是按一个比较稳定的耗油速度减少的，所以这里引入油量可信采样值的概念：可信采样值序列vi在一定的范围内变化，当它是往下变化时，可信采样值vi相对基准油量最小不低于滤波队列下限δlower，也是可以往上变化的，但最大不能超过滤波队列上限δupper。
[0076]
粒子软件滤波器开始运行后，得到的采样值si，必须经过k个测量周期将数据队列填满，作为后续比较的基础。在队列填满的基础上，首先将采样值si与滤波极限δlimit比较可将明显的虚假数据拒之门外；用冒泡法进行排序，如果采样值si小于滤波极限δlimit且又在[δlower，δupper]区间，就认为它是可信的进入数据队列，然后利用滤波队列上限δupper和滤波队列下限δlower以及滤波阈值上限δ

、滤波阈值下限δ-对采样数据进行分类，得到三类粒子集，对中值粒子集不进行重采样；最后对大权值粒子和小权值粒子组成的粒子集先进行判断，若符合重采样条件，则对其使用线性组合方式进行重采样得到新粒子。根据新粒子组成的数据队列，取这k个数据的平均值，同时要将最早的数据从队列中推
出，为新的数据进入留出空间；这个平均值就是一个新的经过软件滤波后的油量值，同时作为下一个滤波周期的基准油量。这样就能在保证测量的跟随性的前提下提高测量的稳定性，比较好的反映出油量的真实变化情况。
[0077]
选择参数：各个参数对测量系统的跟随性和稳定性有直接的影响，恰当的选择参数是系统能否成功的关键。
[0078]
在上述所有的参数中，滤波队列长度k和滤波极限δlimit、滤波队列上限δupper、滤波阈值上限δ

、滤波阈值下限δ-、滤波队列下限δlower对系统的跟随性和稳定性影响最为明显，必须在满足规定的滤波响应时间的基础上，经过详细的分析试飞数据和理论推理才能得到合适的值。
[0079]
本实施例的其他部分与上述实施例1相同，故不再赘述。
[0080]
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。