[发明专利]一种机动车尾气遥测设备校正方法有效
申请号: | 201611267912.0 | 申请日: | 2016-12-31 |
公开(公告)号: | CN106644951B | 公开(公告)日: | 2019-04-26 |
发明(设计)人: | 康宇;李泽瑞;陈国勇;谭小彬;曹洋 | 申请(专利权)人: | 中国科学技术大学 |
主分类号: | G01N21/17 | 分类号: | G01N21/17 |
代理公司: | 北京科迪生专利代理有限责任公司 11251 | 代理人: | 杨学明;顾炜 |
地址: | 230026 安*** | 国省代码: | 安徽;34 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种基于遥测数据的机动车尾气遥测设备校正方法,首先给出尾气遥测过程的状态空间模型,然后使用自调节阻差卡尔曼滤波方法进行估计。利用噪声统计估值器和粗大误差判断器,该滤波算法能够处理噪声统计特性时变以及存在测量粗大误差的情况,从而对尾气遥测设备进行校正,给出相对精确的尾气遥测结果。 | ||
搜索关键词: | 尾气遥测 校正 机动车尾气 粗大误差 遥测设备 噪声统计特性 状态空间模型 卡尔曼滤波 滤波算法 遥测数据 噪声统计 估值器 判断器 时变 测量 | ||
【主权项】:
1.一种机动车尾气遥测设备校正方法,其特征在于步骤如下:步骤1:将机动车尾气遥测设备布置在道路上,连续采集经过车辆的尾气中污染物浓度的数据;步骤2:建立机动车尾气遥测过程的状态空间模型;步骤3:使用自调节阻差卡尔曼滤波方法进行尾气遥测数据的校正;所述尾气遥测设备包括移动式、水平式和垂直式机动车尾气遥测设备;所述步骤2中,机动车尾气遥测过程的状态空间模型如下:系统状态方程:Xt+1=Xt+Wt未发生故障时的测量方程:yt+1=Xt+1+Vt+1发生故障时的测量方程:yt+1=Xt+1+Vt+1+gt+1其中:Xt+1为某一辆车尾气中污染物成分的真实浓度值,yt+1为尾气中污染物成分的检测值,用尾气遥测设备获得,Wt和Vt是均值为零的高斯白噪声,方差分别为QW和QV,Wt、Vt两两不相关,即认为t时刻的方差值与t+1时刻的值相等,gt为发生故障时的故障幅值,其幅值远远大于QV;所述步骤3中,自调节阻差卡尔曼滤波方法如下:初始条件:
其中:E表示数学期望,X0表示初始时尾气污染物成分的真实浓度值,
是X0的最优估计,P0是初始估计误差;具体步骤如下:1)状态预测:
其中,
是t+1时刻的预测估计,
是t时刻的最优估计;2)计算预测估计的方差:
其中,Pt,t是t时刻的最优估计方差,Pt+1,t是t+1时刻的预测估计方差,
是QW的估值;3)基于拉伊达准则进行粗大误差判断;4)如果不存在粗大误差,则:4.1)计算卡尔曼增益矩阵:
其中,Kt是t时刻的卡尔曼增益,
是QV的估值;4.2)最优估计:
其中,
是t+1时刻的最优估计;4.3)计算最优估计方差:Pt+1,t+1=[1‑Kt]Pt+1,t其中,Pt+1,t+1是最优估计方差;5)如果存在粗大方差5.1)最优估计:
5.2)计算最优估计方差:Pt+1,t+1=Pt+1,t6)基于检测值相关性进行噪声方差估计;所述步骤3)中的粗大误差判断如下:
即
H1:δt不服从正态分布
由于gt远远大于QV,所以|δt|>
其中H0为未发生故障时,即不存在粗大误差的情况,H1为发生故障时,即存在粗大误差的情况,根据这两个条件将量测值分为正常值和故障值两大类,δt表示新息,是实际输出与预测输出的差值,即
所述步骤6)基于检测值相关性进行噪声方差估计如下:构造观测差值:Zt+1=yt+1‑yt,构造统计量a与b:![]()
其中n为统计数据数目,测量噪声估计:
系统噪声估计:![]()
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学技术大学,未经中国科学技术大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201611267912.0/,转载请声明来源钻瓜专利网。