CN102216749B - 用于校正来自压力传感器的信号的漂移的方法 - Google Patents
用于校正来自压力传感器的信号的漂移的方法 Download PDFInfo
- Publication number
- CN102216749B CN102216749B CN200980145841.9A CN200980145841A CN102216749B CN 102216749 B CN102216749 B CN 102216749B CN 200980145841 A CN200980145841 A CN 200980145841A CN 102216749 B CN102216749 B CN 102216749B
- Authority
- CN
- China
- Prior art keywords
- signal
- spike
- pressure
- eps
- slope
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 15
- 238000002485 combustion reaction Methods 0.000 claims abstract description 12
- 238000005259 measurement Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 description 6
- 230000015654 memory Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005616 pyroelectricity Effects 0.000 description 2
- 206010038743 Restlessness Diseases 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L23/00—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid
- G01L23/08—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid operated electrically
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
一种用于校正来自测量内燃机的气缸中的压力的传感器的信号(Sb)的漂移的方法,其中该信号相当于等式的直线,压力尖峰覆盖在该信号上,所述校正方法包括如下步骤:I:使用快速卡尔曼滤波器用于检测属于压力尖峰的点,II:使用慢速卡尔曼滤波器用于确定斜率(A)和常量(B),III:对于每个点,根据这些点是否属于在步骤I期间确定的检测到的压力尖峰并且根据在步骤II期间确定的斜率(A)和常量(B)的值来校正信号漂移,其特征在于,在步骤I期间:使用所用的快速卡尔曼滤波器来估计信号的点上的预测误差(epsR),该预测误差的标准偏差(eps_sigma)被滤波并且被最大化,根据如下两个标准中的至少一个确定该点处的压力尖峰的开始和/或结束:该点上的预测误差(epsR)在尖峰开始阈值(delta1_up)以上,该点上的预测误差的经滤波的和被最大化的标准偏差(eps_sigma)在尖峰开始标准偏差阈值(eps_sigma_S1)以上。
Description
技术领域
本发明涉及一种用于校正来自压力传感器的信号的漂移的方法。本发明对于测量内燃机的气缸中的压力的压力传感器是特别有用的。
背景技术
在柴油机气缸的燃烧室中占优势的压力由位于例如预热塞(bougie de préchauffage)中的压力传感器测量。给出作为在发动机循环(进气、压缩、燃烧、排气)期间的时间的函数的压力的曲线呈现出理想地是直线的并且以零为中心的基本信号,压力尖峰周期性地覆盖在该基本信号上。这种类型的传感器通常配备有压电敏感元件。
顾及该传感器的操作环境,使该传感器遭受温度和压力变化。特别地,温度变化创建了传感器的压电敏感元件中的热电性(pyro-électricité),该热电性修改了该传感器递送的压力信号的值。因此,给出作为传感器输出时的时间的函数的压力的曲线与在气缸中占优势的真实压力的曲线表现不同。更具体地:
●基本信号不再以零为中心,也就是说,所测量的平均压力值偏移了常量B,
●基本信号不再平行于横坐标轴,也就是说,不再是水平的,而是呈现斜率A并且具有类型 ,
●斜率A和常量B在时间上不是固定的并且因此可以因发动机循环而异。
因此,可以将基本信号看作直线,其等式具有类型,压力尖峰周期性地覆盖在该基本信号上。
因此,必须对该基本信号进行处理,以便向发动机计算机(calculateur moteur)提供真实且可靠的压力测量结果,从而正确地回到为零(或者为预先定义的常量值)的中心位置并且没有时间漂移。
该信号的处理算法必须校正由传感器供给的信号,也就是说必须因此使得能:
●确定B,
●确定A,
●通过确定压力尖峰的位置和持续时间来相对于基本信号的漂移区别峰值。实际上,如果不独立于信号漂移对由于尖峰而引起的压力的突然增大进行处理,则使斜率A和常量B的确定失真,
●从传感器供给的信号减去所确定的基本信号(),以使其回到为零(或者为预先定义的常量值)的中心位置。
信号可以在捕获信号期间并且由传感器直接处理,或者在捕获信号之后由外部微处理器处理。后一种解决方案呈现出一旦捕获信号就执行处理的优点,其中在发动机控制计算机中有计算装置以及需要的时间可用。然而,这也呈现了持久地使计算机的存储容量过载的缺点。
由压力传感器直接处理呈现了诸多限制:该压力传感器必须是快速的、准确的并且受限于所使用的存储容量,因为该压力传感器被并入并不具有功能强大的配备有大存储器的内置计算机的传感器。从现有技术已知,信号的直接处理可以(根据最小二乘法)通过在包含N个点的点的滑动窗口上对线性模型进行估计来完成。
这样的处理的主要缺陷在于对于计算需要大存储容量。可以进行简化和近似以减小这种缺陷,这接着引起了计算精度和稳定性的问题。
可以设想用于处理信号的其它方法,诸如例如使用卡尔曼(Kalman)滤波器。该滤波器依赖于用于在信号及其按增益衰减的预测之间纠错的递归方法。信号的预测是基于在前一测量时刻被滤波和经校正的信号来计算的。然而,将这样的方法应用于包括具有规律间隔的尖峰的压力信号呈现了如下缺陷:
●如果校正过大,则斜率A和常量B的确定失真,因为其由于压力尖峰的存在而被过度估计,
●如果校正是低的,则经校正的信号并未考虑由于压力尖峰的存在造成的误差并且因而接近直线。因此,斜率A和常量B的值是正确的,但是压力尖峰被忽视了。
发明内容
本发明提出了在不要求相当大的计算存储容量的情况下以可靠的方式确定斜率A和常量B的值以及压力尖峰。
本发明的这些目的是借助于用于校正来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)的漂移的方法而实现的,其中该信号相当于形成了由斜率A和常量B的等式的直线表示的基本信号(Sa)的一连串点,压力尖峰覆盖在该基本信号上,所述校正方法包括如下步骤:
●I:使用快速卡尔曼滤波器,也就是说,使用包括斜率(KaR)和常量(KbR)的值接近于1的增益的快速卡尔曼滤波器,用于检测属于压力尖峰的点,
●II:使用慢速卡尔曼滤波器,也就是说,使用包括斜率(KaL)和常量(KbL)的值接近于0的增益的慢速卡尔曼滤波器,用于确定表示基本信号的直线的斜率(A)和常量(B),
●III:对于每个点,根据这些点是否属于在步骤I期间确定的检测到的压力尖峰并且根据在步骤II期间确定的斜率(A)和常量(B)的值来校正信号漂移,以便确定在气缸中占优势的压力的真实信号(Sr)。
所述方法在如下方面是值得注意的:在步骤I期间,
- 使用所用的快速卡尔曼滤波器来估计信号的点上的预测误差(epsR),
- 该预测误差的标准偏差(eps_sigma)被滤波并且被最大化,以估计该点相对于先前的点的稳定性,
- 根据如下两个标准中的至少一个来确定该点处的压力尖峰的开始和/或结束:
○该点上的预测误差(epsR)在尖峰开始阈值(delta1_up)以上,
○该点上的预测误差的经滤波的和被最大化的标准偏差(eps_sigma)在尖峰开始标准偏差阈值(eps_sigma_S1)以上。
优选地,尖峰开始标准偏差阈值等同于经滤波的和被最大化的标准偏差的最后的最小值乘以尖峰开始系数。
以互补的方式,根据如下两个标准中的至少一个在点处确定尖峰的结束:
●该点上的预测误差在尖峰结束阈值以下,
●该点上的误差的经滤波的和被最大化的标准偏差在尖峰结束标准偏差阈值以下。
有利地,尖峰结束标准偏差阈值等同于经滤波的和被最大化的标准偏差的最后的最大值乘以尖峰结束系数。
在另一实施例中,在步骤II期间:
●根据慢速卡尔曼滤波器估计直线的斜率和常量,
●在步骤I期间由快速卡尔曼滤波器确定的属于压力尖峰的点被替换为慢速卡尔曼滤波器通过使用先前估计的斜率和常量所预测的点。
在补充实施例中,在步骤III期间,在从传感器获得的信号中减去(在步骤II期间确定的)直线。
根据本发明,快速卡尔曼滤波器的斜率增益大于慢速卡尔曼滤波器的斜率增益,并且快速卡尔曼滤波器的恒定增益大于慢速卡尔曼滤波器的恒定增益。
明智地,快速卡尔曼滤波器的斜率增益小于快速卡尔曼滤波器的恒定增益,并且慢速卡尔曼滤波器的斜率增益小于慢速卡尔曼滤波器的恒定增益。
本发明还涉及实施呈现出前述特性中的任何一个的本方法的用于校正信号的任何设备,其中该信号可能是压力信号。
因此,本发明应用于包括根据本发明的用于校正压力信号的设备的任何压力信号传感器。
并且本发明还涉及包括根据本发明的用于校正压力信号的设备的任何电子计算机。
附图说明
通过阅读下列作为非限制实例给出的描述并且通过研究附图,本发明的其它特征和优点将变得显而易见,在附图中:
●图1a是内燃机的气缸中的真实压力随着时间推移、即在压缩期间的曲线的示意性表示,
●图1b是内燃机的气缸中的如由压力传感器递送的压力随着时间推移、即在压缩期间的曲线的示意性表示,
●图2a是在没有校正的情况下将卡尔曼滤波器应用于信号的示意性表示,
●图2b是在有校正的情况下将卡尔曼滤波器应用于信号的示意性表示,
●图3a是根据本发明将快速卡尔曼滤波器应用于压力信号的示意性表示,
●图3b是根据本发明将慢速卡尔曼滤波器应用于压力信号的示意性表示,
●图4是根据本发明将卡尔曼滤波器应用于压力尖峰的检测的示意性表示,
●图5是根据本发明的对压力信号进行处理的示意性图解。
具体实施方式
在图1a中表示了给出作为时间的函数的在气缸的燃烧室中占优势的真实压力Sr的变化的曲线。该曲线相当于以零为中心的直线,压力尖峰覆盖在该直线上。为了简化,在图1a中表示了单个压力尖峰。
图1b表示了如由压力传感器测量并供给的噪声信号Sb。
更具体地:
●基本信号Sa不再以零为中心,也就是说,测量的平均压力值偏移了常量B,
●基本信号Sa不再平行于横坐标轴,也就是说不再是水平的,而是呈现出如下类型的增加的斜率A:,
●斜率A和常量B在时间上不是固定的并且因此可以因循环而异。
因此,可以将基本信号Sa看作如下直线:对于该直线,等式具有类型:,压力尖峰周期性地覆盖在该基本信号上。
因此,为了获得表示在气缸中占优势的真实压力Sr的信号,对所测量的信号Sb的校正是必需的。为此,信号的处理从由传感器Sb测量并供给的信号的每个点中减去表示所述信号(就是说Sa)的漂移的直线,以便重新得到不受噪声影响且没有漂移的信号Sr。
图2a和2b图示了将卡尔曼滤波器应用于类型为的信号,其中x表示测量时刻t,以便确定斜率A和常量B。通过将该等式应用于每个测量时刻n并且通过假设斜率A保持恒定而且该斜率A在点n-1和n之间是相同的,获得了如下参考模型(图2a):
。
随后可以通过考虑在点n-1和n之间的时间间隔dt,根据在点n-1处的斜率和常量计算在点n处的常量B:
。
在点n+1处的信号的预测等同于:
(1)
因此,y_pred(n+1)表示根据在点n处确定的参数B和A对在点n+1处的信号的预测。
卡尔曼滤波器的目的在于在点n处将该预测与由传感器供给的受噪声影响的信号Sb的所测量的真实值y_meas(n)进行比较,并且随后校正在点n处的斜率A(n)和常量B(n),使得所预测的信号的值接近传感器所测量的信号Sb的值。
这样,在点n处的预测误差eps因此等同于:
(2)。
如果该误差非零,则在点n处的斜率A(n)并不等于在点n-1处的斜率A(n-1),(参见图2b),并且必须根据在点n处的斜率的预测误差eps(n)对在点n处的斜率A(n)进行校正。
该校正是使用增益Ka完成的,该增益Ka表示期望校正相对于所测量的误差的衰减。
(3)。
相似地,等同于所测量的预测误差的部分的校正被应用于具有增益Kb的常量B,其给出了:
(4)。
增益Ka和Kb的值介于0和1之间。在卡尔曼滤波器的应用中,增益Ka和Kb的调整使得能获得对相对于传感器所测量的值的动态性更大或更小的预测值的校正。因此,Ka和Kb越高(也就是说越接近于1),校正的动态性越大并且越接近所测量的值。另一方面,Ka和Kb越低(也就是说越接近于0),则校正越慢并且与所测量的值保持得越远。
在n+1处应用的用于下一点的预测公式(1)中使用如此校正过的参数A和B(图2b)。
本发明提出了针对压力信号使用该方法,以便可靠地由此确定压力尖峰、斜率A和常量B:
●在第一步骤I期间,第一对高增益KaR、KbR被用于所谓的“快速”滤波器的应用,以便获得接近于所测量的信号y_meas(n)的估计y_predR(n+1)。因此,通过斜率AR和常量BR的值的快速改变,快速地检测到由于尖峰引起的(相对于基本信号Sa)的斜率增加(图3a),
●在第二步骤II期间,(弱于用于“快速”滤波器的增益的)第二对增益KaL、KbL被用于所谓的“慢速”滤波器的应用,以便获得更接近于要被提取的基本信号的估计y_predL(n+1)。在这种情况下,对于信号的预测,未考虑由点y_meas(n)图示的斜率快速增加。斜率AL和常量BL的如此确定的值是基本信号Sa的信号的斜率和常量的正确的值,并且这些值并未因压力尖峰的存在而失真(图3b),
●最后,在最后的步骤III期间,一旦在步骤I中已确定压力尖峰,并且在步骤II期间已确定斜率A和常量B的值,就从压力传感器所测量的信号中减去直线,以便获得在燃烧室中占优势的压力的真实信号。
显然,步骤I和II同时运行,而随后立即在步骤III中进行校正。
增益KaR、KbR、KaL和KbL的值介于0和1之间,并且优选地斜率增益小于相应的恒定增益。
通过针对每个滤波器都应用等式(1)、(2)、(3)、(4),因此获得了如下等式,对于快速滤波器:
并且对于慢速滤波器:
应当注意,参数A(n)、B(n)、eps(n)和y_pred(n)对于每个滤波器都是特定的,因为后者并不提供相同水平的校正。如在图3a和3b中所图示的那样,经由这两个滤波器获得的经校正的点y_predR(n+1)和y_predL(n+1)是不同的。
根据本发明,并且如在图4的注释为4a和4b的曲线图中所图示的那样,快速卡尔曼滤波器被用于检测压力尖峰。先前由快速卡尔曼滤波器确定的预测误差epsR(n)提供了关于信号梯度的稳定性水平并且因此关于任何斜率的快速改变的指示。
受噪声影响的信号Sb的压力尖峰由压力上升、稳定、随后下降来表示。因此,对于压力尖峰,预测误差epsR是包括两个尖峰的信号,正尖峰表示压力尖峰的上升,并且负尖峰表示下降,(参见图4a)。
然而,背景噪声在压力尖峰中、在压力尖峰之前或在压力尖峰之后的出现也随着预测误差epsR信号的(正或负)尖峰而形成。因此,该信号是一连串正的和负的尖峰。因此,经由该预测误差信号确定压力尖峰整体的持续时间是不可能的。
在步骤I期间,本发明提出了如下附加步骤,以便仍然检测压力尖峰整体:
●使用预测误差epsR的平方,因此该信号epsR 2包括表示压力尖峰的两个正的尖峰,(参见图4a)。不幸地,利用该信号确定压力尖峰整体的持续时间仍是不可能的。实际上,该信号仍过零,从而使得不可能使用信号幅度标准,以便确定压力尖峰的持续时间。
●为了获得不过零的信号,预测误差的平方epsR 2(或者标准偏差)被滤波eps_sigma_filt:
为此,应用滤波系数Kys。例如Kys=0.5。
(5)
这具有如下效果:减慢epsR 2的第一尖峰的上升,在不过零的情况下与epsR 2的第二尖峰接合,随后减慢epsR 2的第二尖峰的下降,(图4b)。因此,所获得的信号eps_sigma_filt是并不过零的正的信号。
随后,为了确保快速地检测到尖峰的开始,该滤波器仅在尖峰一旦经过时才被应用,其总计产生如下信号:
●该信号包括针对尖峰的上升的预测误差的平方epsR 2,
●随后该信号包括针对尖峰的下降的预测误差eps_sigma_filt的平方的经滤波的信号。
为了产生该信号,取在该尖峰的持续时间上的这两个值中的最大值。
因此获得了被最大化的经滤波的信号eps_sigma,其等同于:
(6)。
在压力尖峰的持续时间上,被最大化的经滤波的标准偏差eps_sigma是并不通过零返回的正的信号(图4b),可以将与幅度相关的标准应用于该信号,以便确定压力尖峰的开始和结束。
因此,为了检测到压力尖峰的开始,应用如下两个条件中的至少一个:
●如果发现所测量的信号大于被添加常量的预测信号:,也就是说,如果由等式(2)定义的预测误差大于尖峰开始阈值,
●并且如果等式(6)定义的被最大化的经滤波的标准偏差eps_sigma大于尖峰开始标准偏差阈值eps_sigma_S1,,
则检测到压力尖峰的开始,否则压力尖峰尚未开始。
选择尖峰开始标准偏差阈值eps_sigma_S1,使得该尖峰开始标准偏差阈值eps_sigma_S1对于其值来说具有eps_sigma的最后的最小值(也就是说eps_sigma在压力尖峰开始时的值eps_sigma_min(参见图4b))乘以尖峰开始系数delta2_up。
也就是说:,
其中。
相似地,为了检测到尖峰的结束,应用如下两个条件中的至少一个:
●如果所测量的信号小于被添加常量的预测信号:,其等同于,也就是说,如果由等式(2)定义的预测误差在尖峰结束阈值delta1_down以下,
●并且如果由等式(6)定义的被最大化的经滤波的标准偏差eps_sigma小于尖峰结束标准偏差阈值eps_sigma_S2、即,则检测到压力尖峰的结束,否则压力尖峰未结束。
选择尖峰结束标准偏差阈值eps_sigma_S2,使得该尖峰结束标准偏差阈值eps_sigma_S2对于其值来说具有eps_sigma的最后的最大值(也就是说eps_sigma在压力尖峰顶部的值eps_sigma_max(参见图4b))乘以尖峰结束系数delta2_down。
也就是说:
其中。
因此,对于任何点n,可能确定该点n是否属于压力尖峰。
通常,尖峰开始系数delta2_up的值介于0和10之间,尖峰结束系数delta2_down的值介于0和1之间,并且尖峰开始阈值delta1_up和尖峰结束阈值delta1_down的值介于0和5伏特之间。
因此,对于每个点n,如果在该点处检测到尖峰,则该点y_predR(n)的值不能被选择用于估计斜率A和常量B,但是另一方面,如果在该点处未检测到尖峰,则该点的值可以被选择用于正确估计斜率A和常量B。
在步骤II期间,经由慢速卡尔曼滤波器执行根据本发明的斜率A和常量B的确定。实际上:
●通过具有低增益KaL和KbL的卡尔曼滤波器估计斜率和平均常量,以便获得慢校正,其中所述慢校正受信号改变的影响较小、也就是说相对独立于压力尖峰并且因此表示基本信号,
●随后通过使用先前在点n-1处计算的斜率AL和常量BL来预测点n,
●在其中先前经由快速卡尔曼滤波器在点n处已检测到尖峰的情况下,该点y_predR(n)的值被替换为由慢速卡尔曼滤波器预测的点y_predL(n)的值。因此,压力尖峰被替换为具有恒定斜率的信号,也就是说被替换为直线。该预测为了确定不因压力尖峰的存在而失真的斜率AL和常量BL是必需的。
为了提高斜率AL和常量BL的值的精度,该步骤II可以具有变型。实际上,当检测到尖峰开始时,y_predR(n)的值可以被替换为由慢速卡尔曼滤波器在点n-2处预测的值,也就是说被替换为y_predL(n-2)。这是为了对于小的尖峰开始增加不过度估计AL和BL的值。
相似地,为了避免在尖峰结束时对AL和BL的值估计不足,y_predR(n)的值被替换为由慢速卡尔曼滤波器在点n-1处预测的值,也就是说被替换为y_predL(n-1)。
在最后的步骤III期间,从传感器提供的信号中减去如此确定的、表示基本信号Sa的直线,以便重建在燃烧室中占优势的压力的真实曲线,也就是说包括以零为中心的直线的基本信号的真实曲线。
在图5中图示了根据本发明处理信号时的各种步骤,图5包括注释有5a、5b、5c、5d和5e的5个曲线图。
图5a表示由传感器测量的信号,其包括两个压力尖峰和类型为的基本信号的漂移。
图5b和5c图示了步骤I期间的信号处理:
●图5b表示所测量的信号的预测误差epsR,
●图5c表示被最大化的经滤波的预测误差的标准偏差eps_sigma,以及也表示值eps_sigma_min和eps_sigma_max。
在图5d中图示了步骤II。表示了由慢速卡尔曼滤波器获得的基本信号y_predL,其中由快速卡尔曼滤波器获得的压力尖峰被替换为直线。
图5e表示两个尖峰的检测区域D,并且也表示通过执行步骤III如此确定的基本信号y_predL。
因此,本发明使得能在不要求相当大的存储容量的情况下可靠地确定斜率A、常量B和压力尖峰,因为本方法从点n到点n+1一阶递归并且是推测性的,而且并不要求为了应用传统的最小二乘方公式而在多个点的长窗口上进行管理和存储。因此,该方法可以被并入气缸压力传感器中或者发动机计算机中。
显然,本发明并不限于所描述和所表示的实施例,该实施例仅被给出作为实例并且可以例如被应用于包括尖峰的任何测量信号。
Claims (13)
1.一种用于校正来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)的漂移的方法,其中所述来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)相当于形成基本信号(Sa)的一连串点,其中所述基本信号(Sa)由斜率A和常量B的等式 的直线来表示,其中x表示测量时刻,y表示信号,压力尖峰覆盖在所述基本信号上,所述校正方法包括如下步骤:
●I:使用快速卡尔曼滤波器,也就是说,使用其中斜率(KaR)和常量(KbR)的增益的值接近于1的快速卡尔曼滤波器,用于检测属于压力尖峰的点,
●II:使用慢速卡尔曼滤波器,也就是说,使用其中斜率(KaL)和常量(KbL)的增益的值接近于0的慢速卡尔曼滤波器,用于确定表示基本信号的直线的斜率A和常量B,
●III:对于每个点,根据所述点是否属于在步骤I期间确定的检测到的压力尖峰并且根据在步骤II期间确定的斜率A和常量B的值来校正信号漂移,以便确定在气缸中占优势的压力的真实信号(Sr),
其特征在于,在步骤I期间:
●使用所用的快速卡尔曼滤波器来估计信号的点上的预测误差(epsR),
●所述预测误差的标准偏差(eps_sigma)被滤波并且被最大化,以估计该点相对于先前的点的稳定性,
●根据如下两个标准中的至少一个确定该点处的压力尖峰的开始:
- 该点上的预测误差(epsR)在尖峰开始阈值(delta1_up)以上,
- 该点上的预测误差的经滤波的和被最大化的标准偏差(eps_sigma)在尖峰开始标准偏差阈值(eps_sigma_S1)以上。
2.根据权利要求1所述的方法,其特征在于,尖峰开始标准偏差阈值(eps_sigma_S1)等同于经滤波的和被最大化的标准偏差的最后的最小值(eps_sigma_min)乘以尖峰开始系数(delta2_up)。
3.根据权利要求2所述的方法,其特征在于,尖峰开始系数(delta2_up)的值介于0和10之间。
4.根据权利要求1所述的方法,其特征在于,在步骤I期间,根据如下两个标准中的至少一个在点处确定尖峰的结束:
●该点上的预测误差(epsR)在尖峰结束阈值(delta1_down)以下,
●该点上的预测误差的经滤波的和被最大化的标准偏差(eps_sigma)在尖峰结束标准偏差阈值(eps_sigma_S2)以下。
5.根据权利要求4所述的方法,其特征在于,尖峰结束标准偏差阈值(eps_sigma_S2)等同于经滤波的和被最大化的标准偏差的最后的最大值(eps_sigma_max)乘以尖峰结束系数(delta2_down)。
6.根据前述权利要求中的任一项所述的方法,其特征在于,在步骤II期间:
●根据慢速卡尔曼滤波器估计直线的斜率A和常量B,
●属于在步骤I期间由快速卡尔曼滤波器确定的压力尖峰的点被替换为慢速卡尔曼滤波器通过使用先前估计的斜率A和常量B预测的点。
7.根据权利要求6所述的方法,其特征在于,在步骤III期间,在从传感器获得的信号中减去在步骤II期间确定的所预测的直线,其中t表示测量时刻,y表示信号,其中AL是先前估计的斜率A,BL是先前估计的常量B。
8.根据权利要求1至5之一所述的方法,其特征在于,快速卡尔曼滤波器的斜率增益(KaR)大于慢速卡尔曼滤波器的斜率增益(KaL)。
9.根据权利要求1至5之一所述的方法,其特征在于,快速卡尔曼滤波器的恒定增益(KbR)大于慢速卡尔曼滤波器的恒定增益(KbL)。
10.根据权利要求1至5之一所述的方法,其特征在于,快速卡尔曼滤波器的斜率增益(KaR)小于快速卡尔曼滤波器的恒定增益(KbR)。
11.根据权利要求1至5之一所述的方法,其特征在于,慢速卡尔曼滤波器的斜率增益(KaL)小于慢速卡尔曼滤波器的恒定增益(KbL)。
12.一种用于校正来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)的漂移的设备,其中所述来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)相当于形成基本信号(Sa)的一连串点,其中所述基本信号(Sa)由斜率A和常量B的等式 的直线来表示,其中x表示测量时刻,y表示信号,压力尖峰覆盖在所述基本信号上,所述设备包括如下装置:
用于使用快速卡尔曼滤波器的装置,也就是说,使用其中斜率(KaR)和常量(KbR)的增益的值接近于1的快速卡尔曼滤波器,用于检测属于压力尖峰的点,
用于使用慢速卡尔曼滤波器的装置,也就是说,使用其中斜率(KaL)和常量(KbL)的增益的值接近于0的慢速卡尔曼滤波器,用于确定表示基本信号的直线的斜率A和常量B,
用于对于每个点根据所述点是否属于由用于使用快速卡尔曼滤波器的装置确定的检测到的压力尖峰并且根据由用于使用慢速卡尔曼滤波器的装置确定的斜率A和常量B的值来校正信号漂移的装置,以便确定在气缸中占优势的压力的真实信号(Sr),
其特征在于,所述设备还包括:
用于使用所用的快速卡尔曼滤波器来估计信号的点上的预测误差(epsR)的装置,
用于滤波并且最大化所述预测误差的标准偏差(eps_sigma)的装置,以估计该点相对于先前的点的稳定性,
用于根据如下两个标准中的至少一个确定该点处的压力尖峰的开始的装置:
- 该点上的预测误差(epsR)在尖峰开始阈值(delta1_up)以上,
- 该点上的预测误差的经滤波的和被最大化的标准偏差(eps_sigma)在尖峰开始标准偏差阈值(eps_sigma_S1)以上。
13.一种压力信号传感器,其包括根据权利要求12所述的用于校正来自测量内燃机的气缸中的压力的压力传感器的信号(Sb)的漂移的设备。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0806463A FR2938645B1 (fr) | 2008-11-19 | 2008-11-19 | Procede de correction de la derive du signal d'un capteur de pression |
FR0806463 | 2008-11-19 | ||
PCT/EP2009/007825 WO2010057571A1 (fr) | 2008-11-19 | 2009-11-02 | Procede de correction de la derive du signal d'un capteur de pression |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102216749A CN102216749A (zh) | 2011-10-12 |
CN102216749B true CN102216749B (zh) | 2015-02-25 |
Family
ID=40792626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980145841.9A Expired - Fee Related CN102216749B (zh) | 2008-11-19 | 2009-11-02 | 用于校正来自压力传感器的信号的漂移的方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20110264392A1 (zh) |
CN (1) | CN102216749B (zh) |
FR (1) | FR2938645B1 (zh) |
WO (1) | WO2010057571A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2964738B1 (fr) * | 2010-09-10 | 2013-06-21 | Yzatec | Procede et circuit de traitement d'un signal delivre par un capteur piezoelectrique et dispositif de mesure de pression pour machine a piston(s) |
CN102564486A (zh) * | 2011-12-21 | 2012-07-11 | 上海电机学院 | 一种传感器慢偏故障的校正方法 |
FR2995681B1 (fr) * | 2012-09-20 | 2014-09-05 | Continental Automotive France | Procede de traitement d'un signal d'un dispositif de mesure de pression au sein d'un moteur a combustion interne |
US9261026B2 (en) * | 2013-06-27 | 2016-02-16 | Pratt & Whitney Canada Corp. | System and method for conditioning noisy signals |
FR3011581B1 (fr) * | 2013-10-08 | 2018-08-24 | Continental Automotive France | Procede de compensation d'un signal d'un dispositif de mesure de pression au sein d'un moteur a combustion interne |
US10054511B2 (en) * | 2014-10-01 | 2018-08-21 | Sensata Technologies, Inc. | Pressure sensor with correction of offset drift in cyclic signal |
JP5923585B2 (ja) * | 2014-10-17 | 2016-05-24 | 日本写真印刷株式会社 | 圧力検出装置、圧力検出装置の制御方法、及びプログラム |
FR3028036B1 (fr) | 2014-11-03 | 2016-12-09 | Continental Automotive France | Procede de traitement d'un signal de tension relatif a la pression regnant dans une chambre de combustion d'un cylindre d'un moteur a combustion interne |
CN105043657B (zh) * | 2015-08-21 | 2017-09-05 | 麦克传感器股份有限公司 | 一种智能数显压力变送器手动清零的方法 |
FR3047072B1 (fr) * | 2016-01-21 | 2018-01-26 | Continental Automotive France | Procede et dispositif de traitement d'un signal fourni par un capteur de mesure de la pression regnant dans un cylindre |
CN107941417B (zh) * | 2017-11-10 | 2024-05-07 | 苏州华兴源创科技股份有限公司 | 一种压力传感器的输出校准装置及方法 |
CN109993088B (zh) * | 2019-03-22 | 2021-05-04 | 江南大学 | 一种无线传感网络数据漂移盲校准方法 |
CN112932475B (zh) * | 2021-02-01 | 2023-02-21 | 武汉泰利美信医疗科技有限公司 | 血氧饱和度的计算方法、装置、电子设备及存储介质 |
CN113640566B (zh) * | 2021-07-08 | 2024-04-26 | 国网江苏省电力有限公司电力科学研究院 | 一种foct漂移故障特征提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19900738C1 (de) * | 1999-01-12 | 2000-06-15 | Daimler Chrysler Ag | Verfahren und Vorrichtung zur Bestimmung eines Brennraumdruckverlaufs bei einer Brennkraftmaschine |
FR2872282A1 (fr) * | 2004-06-28 | 2005-12-30 | Renault Sas | Procede de traitement d'un signal de pression |
FR2878030A1 (fr) * | 2004-11-18 | 2006-05-19 | Renault Sas | Dispositif de filtrage d'un signal de mesure de pression |
CN1796756A (zh) * | 2004-12-27 | 2006-07-05 | 本田技研工业株式会社 | 缸内压的检测 |
DE102007045222A1 (de) * | 2007-09-21 | 2008-03-06 | Daimler Ag | Verfahren zur Korrektur eines Brennrauminnendrucksignals einer Verbrennungskraftmaschine |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7672797B2 (en) * | 2006-07-20 | 2010-03-02 | Ads Llc | Flow measurement in partially filled pipes using pulsed peak velocity doppler |
US7689343B2 (en) * | 2007-04-24 | 2010-03-30 | Gm Global Technology Operations, Inc. | Method and apparatus for enabling control of fuel injection for an engine operating in an auto-ignition mode |
US8072204B1 (en) * | 2007-05-09 | 2011-12-06 | Zilker Labs, Inc. | Control system optimization via digital diode emulation |
-
2008
- 2008-11-19 FR FR0806463A patent/FR2938645B1/fr not_active Expired - Fee Related
-
2009
- 2009-11-02 CN CN200980145841.9A patent/CN102216749B/zh not_active Expired - Fee Related
- 2009-11-02 US US13/129,736 patent/US20110264392A1/en not_active Abandoned
- 2009-11-02 WO PCT/EP2009/007825 patent/WO2010057571A1/fr active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19900738C1 (de) * | 1999-01-12 | 2000-06-15 | Daimler Chrysler Ag | Verfahren und Vorrichtung zur Bestimmung eines Brennraumdruckverlaufs bei einer Brennkraftmaschine |
FR2872282A1 (fr) * | 2004-06-28 | 2005-12-30 | Renault Sas | Procede de traitement d'un signal de pression |
FR2878030A1 (fr) * | 2004-11-18 | 2006-05-19 | Renault Sas | Dispositif de filtrage d'un signal de mesure de pression |
CN1796756A (zh) * | 2004-12-27 | 2006-07-05 | 本田技研工业株式会社 | 缸内压的检测 |
DE102007045222A1 (de) * | 2007-09-21 | 2008-03-06 | Daimler Ag | Verfahren zur Korrektur eines Brennrauminnendrucksignals einer Verbrennungskraftmaschine |
Also Published As
Publication number | Publication date |
---|---|
US20110264392A1 (en) | 2011-10-27 |
FR2938645A1 (fr) | 2010-05-21 |
FR2938645B1 (fr) | 2012-03-02 |
WO2010057571A1 (fr) | 2010-05-27 |
CN102216749A (zh) | 2011-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102216749B (zh) | 用于校正来自压力传感器的信号的漂移的方法 | |
US8342011B2 (en) | Method for determining a value representative of the pressure in a combustion chamber of an internal combustion engine | |
US7669583B2 (en) | Internal combustion engine control apparatus | |
CA2673216A1 (en) | Method and apparatus for reconstructing in-cylinder pressure and correcting for signal decay | |
US7124020B2 (en) | System for calibrating apparatus for acquiring the pressure in a motor vehicle diesel engine cylinder | |
US6711489B2 (en) | Method for estimating engine cylinder variables using second order sliding modes | |
US7194899B1 (en) | Method of estimating crack angles and rotation speeds of engine | |
US7295917B2 (en) | Method for determining a combustion chamber pressure | |
US10309854B2 (en) | Method for processing a signal of a pressure measuring device inside an internal combustion engine | |
US7454954B2 (en) | Knocking presence evaluation circuit for an internal combustion engine, knocking identification and control system and corresponding pressure signal processing method | |
JP2016151220A (ja) | 内燃機関の制御装置 | |
US8849591B2 (en) | Method for processing a signal from a flow meter for measuring a gas flow in an internal combustion engine | |
US20120160008A1 (en) | Apparatus for analyzing gas information | |
US7669456B2 (en) | Method and device for determining the rest time and/or for drift compensation in a combustion chamber pressure sensor | |
JP5257777B2 (ja) | 内燃機関の制御装置 | |
CN101688486B (zh) | 校准燃料喷射的方法、记录支持和装置 | |
JP2008520876A (ja) | 内燃機関の燃焼性制御を援助する加速度計信号を処理するためのシステム及び方法 | |
CN110242429A (zh) | 内燃机的控制装置 | |
JP6497283B2 (ja) | データ解析装置 | |
CN102192838A (zh) | 用于确定缸内压力传感器的压力偏移的方法 | |
EP3237871B1 (en) | A method of calibrating a pressure sensor and an internal combustion piston engine | |
JP4684944B2 (ja) | 内燃機関の制御装置 | |
JP6519520B2 (ja) | 内燃機関の制御装置 | |
WO2017091130A1 (en) | Method and device for determining in-cylinder pressure of a combustion engine | |
JP7055159B2 (ja) | 失火判定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221213 Address after: Toulouse, France Patentee after: WeiPai technology France Ltd. Address before: Toulouse, France Patentee before: CONTINENTAL AUTOMOTIVE FRANCE |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150225 |
|
CF01 | Termination of patent right due to non-payment of annual fee |