一种pH计电压误差的双FIR测量方法与流程

本发明属于仪器测量技术领域，尤其涉及一种ph计电压误差的双fir测量方法，利用ph计检定仪进行ph计电压测量、以及利用fir滤波算法和fir平滑算法对测量数据进行融合处理，以实现ph计电压误差的精准测量。

背景技术：

本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。

ph计是利用原电池的原理来测定溶液酸碱度值的电子仪器，ph计的电压误差与测量结果直接相关，因此需要对ph计的电压误差进行测量和校准。ph计检定仪是一种标准直流电压信号输出装置，主要用于检定各种实验室ph计、便携式ph计和实验室通用离子计的电计部分。ph计检定仪通常由两部分组成，一是ph计检定仪的标准直流电压信号发生系统，二是ph计检定仪的高阻器，也称为ph计输入阻抗检测系统。当前ph计的电压误差测量仅采用ph计检定仪、通过多次测量求平均值的方法，不能有效的降低野值对ph计电压测量的影响，最终造成了ph计的电压误差不能准确校准。

技术实现要素：

为了解决上述问题，本发明提出了一种ph计电压误差的双fir测量方法，加入了双fir(有限脉冲响应)测量方法，以当前时刻的电压误差作为状态向量，将通过设备测量得到的电压误差值作为观测向量，通过fir滤波算法预估进行数据融合，在此基础上，利用fir平滑算法对fir滤波算法预估的电压误差值进行平滑，然后对平滑之后的电压误差值进行数值平均，最终得到ph计电压误差的最优测量值。本发明所采用的技术方案如下：

一种ph计电压误差的双fir测量方法，所述的双fir包括fir滤波算法和fir平滑算法两部分，具体包括以下步骤：

将k时刻的电压误差作为fir滤波算法的状态向量xk；

将ph计检定仪与被检的ph计连接，利用ph计检定仪向被检的ph计输入标准电压v1，读取被检的ph计的电压v2；将当前时刻从被检的ph计读出的电压v2与ph计检定仪输入的标准电压v1做差，得到电压误差v3作为fir滤波算法k时刻的观测向量yk，通过fir滤波算法对被检的ph计k时刻的电压误差进行预估；

对通过fir滤波算法预估得到的被检的ph计电压误差进行fir平滑，并对之前所有时刻经fir平滑过的被检的ph计电压误差取平均值，最终得到被检的ph计电压误差的最优测量值。

本发明的有益效果：

本发明以当前时刻的电压误差作为状态向量，构建相应的状态方程；将通过设备(ph计检定仪)测量得到的电压值作为观测向量，通过fir滤波算法预估进行数据融合，在此基础上，利用fir平滑算法对fir滤波算法预估的电压值进行平滑，最后进行数值平均，最终得到ph计电压误差的最优测量值。

本发明将fir滤波算法预估与fir平滑算法两部分有机结合，有效降低了野值对电压误差测量的影响，最终实现了ph计的电压误差的准确校准。

本发明可用于ph计的电压测量与校准。

附图说明

图1为本发明实施例的一种ph计电压误差的双fir测量方法的原理示意图。

具体实施方式

应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本发明使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

实施例一

在一个或多个实施方式中，公开了一种ph计电压误差的双fir测量方法，如图1所示。

一种ph计电压误差的双fir测量方法，所述的双fir包括fir滤波算法和fir平滑算法两部分，具体包括以下步骤：

将k时刻的电压误差作为fir滤波算法的状态向量xk；这个状态向量xk是一个理论值，是公式推导的技术问题，属于公知技术。

所述fir滤波算法的状态方程为：

xk|k-1＝xk-1+ωk；

其中，xk|k-1为由k-1时刻到k时刻的电压误差；xk-1为k-1时刻的电压误差，ωk为k时刻的系统噪声，其协方差阵为q。

将ph计检定仪与被检的ph计连接，利用ph计检定仪向被检的ph计输入标准电压v1，读取被检的ph计的电压v2；将当前时刻从被检的ph计读出的电压v2与ph计检定仪输入的标准电压v1做差，得到电压误差v3作为fir滤波算法k时刻的观测向量yk，通过fir滤波算法对被检的ph计k时刻的电压误差进行预估。

fir滤波算法的观测方程为：

yk＝xk|k-1+νk；

其中，yk为k时刻测量得到的电压误差；νk为观测噪声，其协方差阵为r。

作为优选的实施例，对于测量l次，利用fir滤波算法对ph计的电压误差值进行预估，具体为：

xk|k-1＝xk-1+ωk；

pk|k-1＝pk-1+q；

gk＝[(gk-1)^-1]^-1；

kk＝gk；

xk＝xk|k-1+kk(yk-xk|k-1)；

pk为k时刻的误差矩阵，pk|k-1为由k-1到k时刻的误差矩阵；gk为k时刻的广义噪声功率增益，kk为k时刻的误差增益，(·)^t为矩阵转置，其中，x0和p0需要用户提前设定。x0为xk的初始值，设置为0；p0为pk的初始值，设置为i；这里面i和t都是数学符号。

对通过fir滤波算法预估得到的被检的ph计电压误差进行fir平滑，并对之前所有时刻经fir平滑过的被检的ph计电压误差取平均值，最终得到被检的ph计电压误差的最优测量值。

首先通过获得电压误差的预估，从l次开始向前平滑，利用fir平滑算法对fir滤波算法预估的电压误差值进行平滑，具体为：

为由k到k-1时刻的平滑状态向量；为k时刻的平滑状态向量；为由k到k-1时刻的误差向量；为k时刻的平滑误差矩阵；为k-1时刻的广义平滑噪声功率增益；为k时刻的广义平滑噪声功率增益；为k-1时刻的平滑误差增益、为k-1时刻的平滑状态预估、为k-1时刻的平滑误差矩阵；

在此基础上，最优的平滑预估为：

kf为fir平滑算法的误差增益；xf为fir平滑算法的状态向量。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。