[發明專利]一種機動車尾氣遙測設備校正方法有效
| 申請號: | 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,t
6)基于檢測值相關性進行噪聲方差估計;
所述步驟3)中的粗大誤差判斷如下:
即
H1:δt不服從正態分布由于gt遠遠大于QV,所以|δt|>
其中H0為未發生故障時,即不存在粗大誤差的情況,H1為發生故障時,即存在粗大誤差的情況,根據這兩個條件將量測值分為正常值和故障值兩大類,δt表示新息,是實際輸出與預測輸出的差值,即
所述步驟6)基于檢測值相關性進行噪聲方差估計如下:
構造觀測差值:Zt+1=yt+1-yt,
構造統計量a與b:
其中n為統計數據數目,
測量噪聲估計:
系統噪聲估計:
該專利技術資料僅供研究查看技術是否侵權等信息,商用須獲得專利權人授權。該專利全部權利屬于中國科學技術大學,未經中國科學技術大學許可,擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作,請聯系【客服】
本文鏈接:http://www.szxzyx.cn/pat/books/201611267912.0/1.html,轉載請聲明來源鉆瓜專利網。





