CN107607492B - 机动车尾气标准的检测方法及设备 - Google Patents

机动车尾气标准的检测方法及设备 Download PDF

Info

Publication number
CN107607492B
CN107607492B CN201710828659.XA CN201710828659A CN107607492B CN 107607492 B CN107607492 B CN 107607492B CN 201710828659 A CN201710828659 A CN 201710828659A CN 107607492 B CN107607492 B CN 107607492B
Authority
CN
China
Prior art keywords
curve
fitting
motor vehicle
standard
data
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.)
Active
Application number
CN201710828659.XA
Other languages
English (en)
Other versions
CN107607492A (zh
Inventor
陈文亮
王小蒙
刘洪涛
张金永
樊海春
张涛
汪锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TIANJIN TONGYANG TECHNOLOGY DEVELOPMENT CO LTD
Original Assignee
TIANJIN TONGYANG TECHNOLOGY DEVELOPMENT CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by TIANJIN TONGYANG TECHNOLOGY DEVELOPMENT CO LTD filed Critical TIANJIN TONGYANG TECHNOLOGY DEVELOPMENT CO LTD
Priority to CN201710828659.XA priority Critical patent/CN107607492B/zh
Publication of CN107607492A publication Critical patent/CN107607492A/zh
Application granted granted Critical
Publication of CN107607492B publication Critical patent/CN107607492B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明公开了一种机动车尾气标准的检测方法及设备。该方法包括:采集N组机动车尾气的烟度数据;采用不同类型的曲线函数,通过最小二乘法对所述N组数据进行曲线拟合,生成与实际数据拟合度最高的曲线方程;计算拟合曲线的曲边梯形与实际数据的直角梯形的面积比,得出尾气的不透光度;将所述尾气的不透光度与参比数据进行对比,判断尾气标准是否合格。所述检测设备包括数据采集设备、数据处理设备、数据分析设备和结果显示设备。本发明可以应用于机动车尾气检测系统中,用于快速判定机动车尾气是否合格。

Description

机动车尾气标准的检测方法及设备
技术领域
本发明属于机动车尾气检测技术领域,具体涉及一种通过曲线拟合计算不透光烟度对尾气标准进行检测的方法及设备。
背景技术
机动车行驶过程中将燃料转化为动力,同时向环境中排放污染气体。随着汽车工业的飞速发展,机动车尾气成为影响空气质量的重要因素之一。根据国家权威部门的统计数据表明,中国机动车造成的城市大气污染高达30%以上,机动车排放物中大量的PM、CO、NOx、HC对环境及人体系统造成极大的污染和伤害。在此背景下2016年1月1日国家颁布了中国历史上最严格的《大气污染防治法》,该法中针对机动车尾气治理的处罚规定得相当严格,这就意味着针对大气污染中主要的污染源机动车排放将会实施严格的管控措施。
目前,在对机动车尾气的监管上,国内外对机动车尾气标准的测量主要采用不透光烟度检测方法。对尾气不透光烟度检测室获取污染数据并实行相应措施的重要环节。目前为止,针对车辆尾气不透光烟度的检测一般采用静态检测方法,比如:滤纸烟度法、林格曼烟度法,以及烟度计法。然而,这些传统的检测方法存在以下缺点:
(1)滤纸烟度法测量效率低、数据不具备实时性;
(2)林格曼烟度法受操作人员的主观因素和环境因素的影响,数据可靠性低;
(3)烟度计法需要停车检测,因此检测效率低。
由于汽车运行速度很快,为达到实时、准确、高效率的检测结果,检测系统必须在100毫秒内完成尾气合格性的检测。然而现行的检测方法均无法实现这一目的。
发明内容
(一)要解决的技术问题
有鉴于此,本发明提供一种机动车尾气标准的检测方法及设备,以解决上述至少一项技术问题。
(二)技术方案
为实现上述目的,作为本发明的一个方面,本发明提出了一种机动车尾气标准的检测方法,包括:
步骤1:采集机动车尾气数据;
步骤2:对尾气数据进行曲线拟合,确定最优拟合曲线,所述曲线拟合采用的曲线拟合算法为最小二乘法;
步骤3:确定不透光烟度L,将不透光烟度与参比数据进行对比,得出尾气合格或不合格的结论。
在进一步实施方案中,所述尾气数据包括尾气烟度数据,该尾气烟度数据包含N组数据。
在进一步实施方案中,所述曲线拟合包括以下步骤:
步骤21:建立直角坐标系,将所述N组尾气烟度数据在直角坐标系中标注为对应的N个点(xi,yi),其中i取值为1~N;
步骤22:对所述N个点采用最小二乘法进行曲线拟合,得到拟合曲线。
步骤23:通过步骤2进行四次曲线拟合、五次曲线拟合、六次曲线拟合、对数函数拟合、指数函数拟合和幂函数拟合,分别得到第一拟合曲线、第二拟合曲线、第三拟合曲线、第四拟合曲线、第五拟合曲线和第六拟合曲线;
步骤24:对上述第一拟合曲线、第二拟合曲线、第三拟合曲线、第四拟合曲线、第五拟合曲线和第六拟合曲线进行拟合优度对比,从中选取拟合优度最接近1的作为最优拟合曲线。
在进一步实施方案中,所述最小二乘法进行曲线拟合步骤如下:
步骤221:设拟合多项式为:
Y=a0+a1x+...+akxk (a)
上述拟合多项式(a)对应拟合曲线Y。
步骤222:计算所述直角坐标系中的点(xi,yi)(i取1~N)到上述拟合曲线Y的距离之和,即偏差平方和如下:
Figure BDA0001408147090000031
步骤223:为了求得符合条件的系数值a0、a1…ak,将公式(b)等式左右两边分别对a0、a1、a2…ak求偏导数,依次得到:
对a0求偏导数:
Figure BDA0001408147090000032
对a1求偏导数:
Figure BDA0001408147090000033
对a2求偏导数:
Figure BDA0001408147090000034
.......
对ak求偏导数:
Figure BDA0001408147090000035
步骤224:将上述方程(c0)~(ck)进行变形,得到下面的等式:
Figure BDA0001408147090000036
Figure BDA0001408147090000037
Figure BDA0001408147090000038
.......
Figure BDA0001408147090000039
步骤225:将上述方程(d0)~(dk)等式右边的多项式进行拆分整理,得到下面的等式:
Figure BDA0001408147090000041
Figure BDA0001408147090000042
Figure BDA0001408147090000043
.......
Figure BDA0001408147090000044
步骤226:将上述方程(e1)~(ek)表示成如下所示的范德蒙德矩阵的形式:
Figure BDA0001408147090000045
步骤227:将上述范德蒙德矩阵(f)进行简化:
Figure BDA0001408147090000046
上述公式(g)实际上为矩阵方程:
X*A=Y (h)
步骤228:由方程(h)可得系数矩阵A:
Figure BDA0001408147090000051
由系数矩阵A就可以确定拟合曲线Y=a0+a1x+...+akxk
作为更进一步的实施方案,所述不透光烟度L的计算过程如下:
步骤31:预设一个标准值y0,将yi与y0进行对比:
yi≤y0时,保留对应的点(xi,yi);
yi>y0时,舍去对应的点(xi,yi);
步骤32:在保留的点中,选取其中横坐标xi最小的点作为上边界点A(xa,ya),选取其中横坐标xi最大的点作为下边界点B(xb,yb);所述上边界点A、下边界点B以及横坐标轴围成一个直角梯形;所述最优拟合曲线在上边界A点和下边界B点之间与横坐标轴围成一个曲边梯形;
步骤33:计算直角梯形的面积S1
S1=(ya+yb)(xb-xa)/2 (j)
步骤34:用程序模拟分割梯形求面积的算法计算曲边梯形的面积S2
步骤35:计算不透光烟度L:
L=S2/S1 (k)
作为更进一步的技术方案,判断机动车尾气标准是否合格的过程如下:
预设一个参比数据Ref;
将所述不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
在进一步实施方案中,所述结果显示单元对机动车尾气标准的判断结果进行如下显示:
接收到机动车尾气标准的判断结果为合格时,通过文字W1和/或颜色C1和/或语音V1等给出合格指示信息;
接收机动车尾气标准的判断结果为不合格时,通过文字W2和/或颜色C2和/或语音V2等给出不合格指示信息。
在更进一步实施方案中,所述文字W1和语音V1包括“合格”、“通过”、“Yes”和“Positive”等;所述文字W2和语音V2包括“不合格”、“不通过”、“No”和“Negative”等;所述颜色C1为区别于颜色C2的颜色,如绿色、蓝色等;所述颜色C2为区别于颜色C1的颜色,如红色、黄色等。
为实现上述技术目的,作为本发明的另一个方面,本发明提出了一种检测机动车尾气标准的设备,该设备包括:
数据采集单元,用于采集机动车尾气的烟度数据,数据样品数量为N组;
数据处理单元,对所述N组烟度数据建立直角坐标系,得到N个点(xi,yi)(i取1~N),并采用最小二乘法进行曲线拟合,从中选出最优拟合曲线;
数据分析单元对N个点(xi,yi)(i取1~N)进行筛选、计算不透光烟度L,并将不透光烟度L与预设的参比数据Ref进行对比,判断该机动车尾气标准是否合格;
结果显示单元,用于给出该机动车尾气标准是否合格的指示信息;
其中,数据采集单元和数据处理单元之间设置有数据上传串口;
数据处理单元和数据分析单元之间设置有数据传输串口;
数据分析单元和结果显示单元之间设置有数据输出串口。
在进一步实施方案中,所述最小二乘法进行曲线拟合步骤如下:
(1)设拟合多项式为:
Y=a0+a1x+...+akxk (a)
上述拟合多项式(a)对应拟合曲线Y。
(2)计算所述直角坐标系中的点(xi,yi)(i取1~N)到上述拟合曲线Y的距离之和,即偏差平方和如下:
Figure BDA0001408147090000061
(3)为了求得符合条件的系数值a0、a1…ak,将公式(b)等式左右两边分别对a0、a1、a2…ak求偏导数,依次得到:
对a0求偏导数:
Figure BDA0001408147090000071
对a1求偏导数:
Figure BDA0001408147090000072
对a2求偏导数:
Figure BDA0001408147090000073
.......
对ak求偏导数:
Figure BDA0001408147090000074
(4)将上述方程(c0)~(ck)进行变形,得到下面的等式:
Figure BDA0001408147090000075
Figure BDA0001408147090000076
Figure BDA0001408147090000077
.......
Figure BDA0001408147090000078
(5)将上述方程(d0)~(dk)等式右边的多项式进行拆分整理,得到下面的等式:
Figure BDA0001408147090000079
Figure BDA00014081470900000710
Figure BDA0001408147090000081
.......
Figure BDA0001408147090000082
(6)将上述方程(e1)~(ek)表示成如下所示的范德蒙德矩阵的形式:
Figure BDA0001408147090000083
(7)将上述范德蒙德矩阵(f)进行简化:
Figure BDA0001408147090000084
上述公式(g)实际上为矩阵方程:
X*A=Y (h)
(8)由方程(h)可得系数矩阵A:
Figure BDA0001408147090000085
由系数矩阵A就可以确定拟合曲线Y=a0+a1x+...+akxk
在进一步实施方案中,所述曲线拟合包括四次曲线拟合、五次曲线拟合、六次曲线拟合、对数函数拟合、指数函数拟合和幂函数拟合,分别得到第一拟合曲线、第二拟合曲线、第三拟合曲线、第四拟合曲线、第五拟合曲线和第六拟合曲线。
在进一步实施方案中,对上述第一拟合曲线、第二拟合曲线、第三拟合曲线、第四拟合曲线、第五拟合曲线和第六拟合曲线进行拟合优度对比,从中选取拟合优度最接近1的作为最优拟合曲线。
在进一步实施方案中,所述N个点的筛选过程如下:
(1)预设一个标准值y0
(2)将yi与y0进行对比:
yi≤y0时,保留对应的点(xi,yi);
yi>y0时,舍去对应的点(xi,yi)。
在进一步实施方案中,所述不透光烟度L的计算过程如下:
(1)如图4所示,在上述保留的点中,选取其中横坐标xi最小的点作为上边界点A(xa,ya),选取其中横坐标xi最大的点作为下边界点B(xb,yb),A、B两个点与横坐标轴围成一个直角梯形,所述最优拟合曲线在A、B两个点之间与横坐标轴围成一个曲边梯形;
(2)如图4所示,计算直角梯形的面积S1
Sl=(ya+yb)(xb-xa)/2 (j)
(3)如图4所示,计算曲边梯形的面积S2,以A、B两点为上下限,用程序模拟分隔梯形求面积的算法,计算出曲边梯形的面积S2
(4)计算不透光烟度L:
L=S2/S1 (k)
在进一步实施方案中,所述判断机动车尾气标准是否合格的过程如下:
(1)预设一个参比数据Ref;
(2)将所述不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
在进一步实施方案中,所述结果显示单元给出如下指示信息:
(1)接收到机动车尾气标准的判断结果为合格时,通过文字W1和/或颜色C1和/或语音V1等给出合格指示信息;
(2)接收机动车尾气标准的判断结果为不合格时,通过文字W2和/或颜色C2和/或语音V2等给出不合格指示信息。
在进一步实施方案中,所述文字W1和语音V1包括“合格”、“通过”、“Yes”和“Positive”等;所述文字W2和语音V2包括“不合格”、“不通过”、“No”和“Negative”等;所述颜色C1为区别于颜色C2的颜色,如绿色、蓝色等;所述颜色C2为区别于颜色C1的颜色,如红色、黄色等。
在进一步实施方案中,所述检测设备可以为便携式一体机或现场检测装置。
在进一步实施方案中,所述便携式一体机包括将数据采集单元、数据处理单元、数据分析单元以及结果显示单元通过数据上传串口、数据传输串口和数据输出串口联接在一起的终端,该终端具有适合随身携带的外形尺寸和重量。
所述现场检测装置包括将数据采集单元、数据处理单元、数据分析单元以及结果显示单元中的至少一个进行单独安装的终端,该终端的外形尺寸和重量适合现场监测。
(三)有益效果
(1)本发明通过最小二乘法拟合算法对多种曲线拟合结果进行对比,取最优解的方案,提高了测量结果的准确性。
(2)在实现计算曲边梯形面积时,通过编程实现了微积分中的无限分割梯形以近似计算曲边梯形面积的算法,在保证计算速度的前提下,尽可能使分割的份数更大,以提高计算精度。
(3)本发明提供的不透光烟度计算方法准确,运行速度快,可以实现快速、准确的检测,所有的运算可以在100毫秒内得出结果,能够满足机动车尾气检测的需求。
附图说明
图1为本发明机动车尾气标准检测方法的流程图。
图2为本发明实施例1中现场尾气标准检测设备的示意图。
图3为本发明实施例2中便携式尾气标准检测设备的示意图。
图4为本发明不透光烟度计算方法示意图。
具体实施方式
根据本发明总体上的发明构思,提供一种机动车尾气标准的检测方法,如图1所示,包括以下步骤:
步骤1:数据采集设备1采集机动车尾气数据D1并传给数据处理设备2;
步骤2:数据处理设备2对尾气数据D1进行曲线拟合CF,确定最优拟合曲线Y0,并输出给数据分析设备3;
步骤3:数据分析设备3确定不透光烟度L,将不透光烟度L与参比数据Ref进行对比,得出尾气标准的结论,并输出给结果显示设备4;
其中,所述曲线拟合CF采用的曲线拟合算法为最小二乘法。
进一步地,所述尾气数据D1包括N组尾气烟度数据;
进一步地,所述曲线拟合包括以下步骤:
步骤21:建立直角坐标系,将所述N组尾气烟度数据在直角坐标系中标注为对应的N个点(xi,yi),其中i取值为1~N;
步骤22:对所述N个点采用最小二乘法进行曲线拟合CF,得到拟合曲线Y;
步骤23:通过步骤2进行四次曲线拟合CF1、五次曲线拟合CF2、六次曲线拟合CF3、对数函数拟合CF4、指数函数拟合CF5和幂函数拟合CF6,分别得到第一拟合曲线Y1、第二拟合曲线Y2、第三拟合曲线Y3、第四拟合曲线Y4、第五拟合曲线Y5和第六拟合曲线Y6
步骤24:对上述第一拟合曲线Y1、第二拟合曲线Y2、第三拟合曲线Y3、第四拟合曲线Y4、第五拟合曲线Y5和第六拟合曲线Y6进行拟合优度对比,从中选取拟合优度最接近1的作为最优拟合曲线Y0
进一步地,所述最小二乘法进行曲线拟合CF步骤如下:
步骤221:设拟合多项式为:
Y=a0+a1x+...+akxk (a)
上述拟合多项式(a)对应拟合曲线Y。
步骤222:计算所述直角坐标系中的点(xi,yi)(i取1~N)到上述拟合曲线Y的距离之和,即偏差平方和如下:
Figure BDA0001408147090000121
步骤223:为了求得符合条件的系数值a0、a1…ak,将公式(b)等式左右两边分别对a0、a1、a2…ak求偏导数,依次得到:
对a0求偏导数:
Figure BDA0001408147090000122
对a1求偏导数:
Figure BDA0001408147090000123
对a2求偏导数:
Figure BDA0001408147090000124
.......
对ak求偏导数:
Figure BDA0001408147090000125
步骤224:将上述方程(c0)~(ck)进行变形,得到下面的等式:
Figure BDA0001408147090000126
Figure BDA0001408147090000127
Figure BDA0001408147090000128
.......
Figure BDA0001408147090000129
步骤225:将上述方程(d0)~(dk)等式右边的多项式进行拆分整理,得到下面的等式:
Figure BDA0001408147090000131
Figure BDA0001408147090000132
Figure BDA0001408147090000133
.......
Figure BDA0001408147090000134
步骤226:将上述方程(e1)~(ek)表示成如下所示的范德蒙德矩阵的形式:
Figure BDA0001408147090000135
步骤227:将上述范德蒙德矩阵(f)进行简化:
Figure BDA0001408147090000136
上述公式(g)实际上为矩阵方程:
X*A=Y (h)
步骤228:由方程(h)可得系数矩阵A:
Figure BDA0001408147090000137
由系数矩阵A就可以确定拟合曲线Y=a0+a1x+...+akxk
进一步地,所述不透光烟度L计算过程如下:
步骤31:预设一个标准值y0
步骤32:将yi与y0进行对比:
yi≤y0时,保留对应的点(xi,yi);
yi>y0时,舍去对应的点(xi,yi);
步骤33:在保留的点中,如图4所示,选取其中横坐标xi最小的点作为上边界A点(xa,ya),选取其中横坐标xi最大的点作为下边界B点(xb,yb),A、B两个点与横坐标轴围成一个直角梯形T1,所述最优拟合曲线Yo在A、B两个点之间与横坐标轴围成一个曲边梯形T2
步骤34:计算直角梯形T1的面积S1
Sl=(ya+yb)(xb-xa)/2 (j)
步骤35:计算曲边梯形T2的面积S2,以A、B两点为上下限,用程序模拟分隔梯形求面积的算法,计算出曲边梯形的面积S2
步骤36:计算不透光烟度L:
L=S2/S1 (k)
进一步地,所述数据分析设备3通过如下方式判断机动车尾气标准是否合格:
预设一个参比数据Ref;
将上述计算得到的不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
进一步地,所述结果显示单元U4进行如下显示:
接收到机动车尾气标准的判断结果为合格时,通过文字W1和/或颜色C1和/或语音V1等给出合格指示信息;
接收机动车尾气标准的判断结果为不合格时,通过文字W2和/或颜色C2和/或语音V2等给出不合格指示信息。
优选地,所述文字W1和语音V1包括“合格”、“通过”、“Yes”和“Positive”等;所述文字W2和语音V2包括“不合格”、“不通过”、“No”和“Negative”等;所述颜色C1为区别于颜色C2的颜色,如绿色、蓝色等;所述颜色C2为区别于颜色C1的颜色,如红色、黄色等。
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。应当指出的是,以下所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护范围。
实施例1:
图2为基于本发明检测方法的现场尾气标准检测设备的示意图。如图2所示,该检测设备包括数据采集单元U1、数据处理单元U2、数据分析单元U3和结果显示单元U4,并且上述4个结构单元中至少有一个单独安装;
所述数据采集单元U1和数据处理单元U2之间设置有数据上传串口SP1
所述数据处理单元U2和数据分析单元U3之间设置有数据传输串口SP2
所述数据分析单元U3和结果显示单元U4之间设置有数据输出串口SP3
该设备的外形尺寸和重量适合现场监测,其工作过程如下:
(1)数据采集单元U1进行数据采集:
所述数据采集单元U1包括一智能终端,该智能终端具有操作界面,并与硬件采集设备连接;
在所述操作界面上进行参数设定:如采集方式(手动/自动)、采集尾气数据的样本数N(此处N=24)、积分时间(采集N组尾气数据D1的间隔,可以根据现场车流量等实际情况进行设定)、切分份数(计算面积时切割出的梯形的个数)等;
参数设定完毕后开始进行数据采集,若要中止采集可点击停止按钮。
(2)数据处理单元U2进行数据处理:
数据单元U1采集的24组尾气数据通过数据上传串口SP1传输给数据处理单元U2
进一步地,处理单元U2将上述24组尾气数据构建为直角坐标系中的24个点(xi,yi)(i取1~24);
进一步地,处理单元U2对上述24个点进行曲线拟合CF;
进一步地,所述曲线拟合CF采用最小二乘法实现,步骤如下:
(1)设拟合多项式为:
Y=a0+a1x+...+akxk (a)
上述拟合多项式(a)对应拟合曲线Y。
(2)计算所述直角坐标系中的点(xi,yi)(i取1~24)到上述拟合曲线Y的距离之和,即偏差平方和如下:
Figure BDA0001408147090000161
(3)为了求得符合条件的系数值a0、a1…ak,将公式(b)等式左右两边分别对a0、a1、a2…ak求偏导数,依次得到:
对a0求偏导数:
Figure BDA0001408147090000162
对a1求偏导数:
Figure BDA0001408147090000163
对a2求偏导数:
Figure BDA0001408147090000164
.......
对ak求偏导数:
Figure BDA0001408147090000165
(4)将上述方程(c0)~(ck)进行变形,得到下面的等式:
Figure BDA0001408147090000166
Figure BDA0001408147090000167
Figure BDA0001408147090000171
.......
Figure BDA0001408147090000172
(5)将上述方程(d0)~(dk)等式右边的多项式进行拆分整理,得到下面的等式:
Figure BDA0001408147090000173
Figure BDA0001408147090000174
Figure BDA0001408147090000175
.......
Figure BDA0001408147090000176
(6)将上述方程(e1)~(ek)表示成如下所示的范德蒙德矩阵的形式:
Figure BDA0001408147090000177
(7)将上述范德蒙德矩阵(f)进行简化:
Figure BDA0001408147090000181
上述公式(g)实际上为矩阵方程:
X*A=Y (h)
(8)由方程(h)可得系数矩阵A:
Figure BDA0001408147090000182
由系数矩阵A就可以确定拟合曲线Y=a0+a1x+...+akxk
进一步地,所述曲线拟合CF包括包括四次曲线拟合CF1、五次曲线拟合CF2、六次曲线拟合CF3、对数函数拟合CF4、指数函数拟合CF5和幂函数拟合CF6,分别得到拟合曲线Y1、拟合曲线Y2、拟合曲线Y3、拟合曲线Y4、拟合曲线Y5和拟合曲线Y6,选取其中拟合优度最接近1的作为所述最优拟合曲线Yo
进一步地,所述数据传输串口SP2将上述24个点(xi,yi)(i取1~24)以及最优拟合曲线Yo传输给数据分析单元U3
进一步地,所述数据分析单元(U3)对24个点(xi,yi)(i取1~24)进行筛选;
进一步地,所述筛选过程如下:
预设一个标准值y0
将yi与y0进行对比,yi≤y0时,保留对应的点(xi,yi);yi>y0时,舍去对应的点(xi,yi);
进一步地,如图4所示,在上述保留的点中,选取其中横坐标xi最小的点作为上边界点A(xa,ya),选取其中横坐标xi最大的点作为下边界点B(xb,yb),A、B两个点与横坐标轴围成一个直角梯形T1,所述最优拟合曲线Yo在A、B两个点之间与横坐标轴围成一个曲边梯形T2
进一步地,计算直角梯形T1的面积S1
以ya和yb分别为直角梯形(T1)的上底和下底,以xb-xa为直角梯形(T1)的高,则有
Sl=(ya+yb)(xb-xa)/2 (j)
进一步地,计算曲边梯形T2的面积S2,以A、B两点为上下限,用程序模拟分割矩形求面积的算法,计算出曲边梯形的面积S2
更进一步地,通过计算曲边梯形T2的面积S2与直角梯形T1的面积(S1)的比值,得到不透光烟度L,即:
L=S2/S1 (k)
进一步地,所述数据分析单元U3判断机动车尾气标准是否合格的过程如下:
预设一个参比数据Ref;
将不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
进一步地,所述数据输出串口SP3将上述机动车尾气标准的判断结果传输给结果显示单元U4,该结果显示单元U4给出指示信息。
优选地,所述结果显示单元U4对机动车尾气标准的判断结果进行如下显示:
接收到机动车尾气标准的判断结果为合格时,通过文字W1和/或颜色C1等给出合格指示信息;
接收机动车尾气标准的判断结果为不合格时,通过文字W2和/或颜色C2等给出不合格指示信息。
优选地,所述文字W1和语音V1包括“合格”、“通过”、“Yes”和“Positive”等;所述文字W2和语音V2包括“不合格”、“不通过”、“No”和“Negative”等;所述颜色C1为区别于颜色C2的颜色,如绿色、蓝色等;所述颜色C2为区别于颜色C1的颜色,如红色、黄色等。
本发明通过多种曲线拟合结果的对比,取最优解的方案,提高了测量结果的准确性;采用数学中的最小二乘法拟合算法,通过编程实现拟合算法,具有科学的理论依据;在实现计算曲边梯形面积时,通过编程实现了微积分中的无限分割矩形以近似计算曲边梯形面积的算法,在保证计算速度的前提下,尽可能使分割的份数更大,以提高计算精度。
本发明提供的不透光烟度计算方法准确,运行速度快,可以实现快速、准确的检测,所有的运算可以在100毫秒内得出结果,能够满足机动车尾气检测的需求。
实施例2:
图3为基于本发明检测方法的便携式尾气标准检测设备的示意图。如图3所示,该便携式尾气标准检测设备包括数据采集单元U1、数据处理单元U2、数据分析单元U3和结果显示单元U4
其中,数据采集单元U1和数据处理单元U2之间设置有数据上传串口SP1
数据处理单元U2和数据分析单元U3之间设置有数据传输串口SP2
数据分析单元U3和结果显示单元U4之间设置有数据输出串口SP3
该便携式尾气标准检测设备的工作过程与实施例1相同,区别在于本实施例检测设备的各个结构单元(U1~U4)通过串口(SP1~SP3)联接并安装于同一台终端内,且该终端的体积和重量便于携带,操作起来更加灵活方便。
综上所述,本发明提供了一种机动车尾气标准的检测方法,通过对机动车尾气烟度数据进行采集和曲线模拟,计算尾气的不透光度,进而判断尾气排放是否合格。本发明还提供了基于上述检测方法的检测设备,该检测设备可以是便携式一体机,方便灵活;也可以是现场检测设备,通过在公路上建立机动车尾气遥测门站等方式,使机动车通过检测门站时,可以快速准确地实时采集数据并判定尾气是否合格,对于尾气不合格的机动车,进行抓拍并上传至机动车管理系统,以便追查。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种机动车尾气标准的检测方法,其特征在于,包括以下步骤:
步骤1:采集机动车的尾气数据D1,所述尾气数据D1含N组尾气烟度数据;
步骤2:建立直角坐标系,将尾气数据D1的N组尾气烟度数据在直角坐标系中标注为对应的N个点(xi,yi),其中i取值为1~N,对所述N个点进行曲线拟合CF,确定最优拟合曲线Y0,所述曲线拟合CF采用的曲线拟合算法为最小二乘法;
步骤3:确定不透光烟度L,将不透光烟度L与参比数据Ref进行对比,得出尾气合格或不合格的结论,其中所述不透光烟度L的计算过程如下:
预设一个标准值y0,将yi与y0进行对比:
yi≤y0时,保留对应的点(xi,yi);
yi>y0时,舍去对应的点(xi,yi);
在保留的点中,选取其中横坐标xi最小的点作为上边界A点(xa,ya),选取其中横坐标xi最大的点作为下边界B点(xb,yb);所述上边界A点、下边界B点以及横坐标轴围成一个直角梯形T1;所述最优拟合曲线Yo在上边界A点和下边界B点之间与横坐标轴围成一个曲边梯形T2
计算直角梯形T1的面积S1
用程序模拟分割梯形求面积的算法计算曲边梯形T2的面积S2
计算不透光烟度L:L=S2/S1
2.根据权利要求1所述的机动车尾气标准的检测方法,其特征在于,所述曲线拟合CF包括以下子步骤:
步骤21:建立直角坐标系,将所述N组尾气烟度数据在直角坐标系中标注为对应的N个点(xi,yi),其中i取值为1~N;
步骤22:对所述N个点采用最小二乘法进行曲线拟合CF,得到拟合曲线Y;
步骤23:通过步骤22进行四次曲线拟合CF1、五次曲线拟合CF2、六次曲线拟合CF3、对数函数拟合CF4、指数函数拟合CF5和幂函数拟合CF6,分别得到第一拟合曲线Y1、第二拟合曲线Y2、第三拟合曲线Y3、第四拟合曲线Y4、第五拟合曲线Y5和第六拟合曲线Y6
步骤24:对上述第一拟合曲线Y1、第二拟合曲线Y2、第三拟合曲线Y3、第四拟合曲线Y4、第五拟合曲线Y5和第六拟合曲线Y6进行拟合优度对比,从中选取拟合优度最接近1的作为最优拟合曲线Yo
3.根据权利要求1所述的机动车尾气标准的检测方法,其特征在于,判断机动车尾气标准是否合格的过程如下:
预设一个参比数据Ref;
将所述不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
4.一种机动车尾气标准的检测设备,包括:
数据采集单元U1,用于采集机动车尾气的烟度数据,样品数量为N组;
数据处理单元U2,用于对所述烟度数据建立直角坐标系并赋予点坐标(xi,yi),其中i取1~N;对上述N个点进行曲线拟合CF,并找到最优拟合曲线Y0;所述曲线拟合CF采用最小二乘法算法进行曲线拟合;
数据分析单元U3,用于对上述N个点进行筛选,计算不透光烟度L,并将不透光烟度L与预设的参比数据Ref进行对比,判断该机动车尾气标准是否合格,其中所述不透光烟度L的计算过程如下:
预设一个标准值y0,将yi与y0进行对比:
yi≤y0时,保留对应的点(xi,yi);
yi>y0时,舍去对应的点(xi,yi);
在保留的点中,选取其中横坐标xi最小的点作为上边界A点(xa,ya),选取其中横坐标xi最大的点作为下边界B点(xb,yb);所述上边界A点、下边界B点以及横坐标轴围成一个直角梯形T1;所述最优拟合曲线Yo在上边界A点和下边界B点之间与横坐标轴围成一个曲边梯形T2
计算直角梯形T1的面积S1
用程序模拟分割梯形求面积的算法计算曲边梯形T2的面积S2
计算不透光烟度L:L=S2/S1
结果显示单元U4,用于给出该机动车尾气标准是否合格的指示信息。
5.根据权利要求4所述的一种机动车尾气标准的检测设备,其特征在于:所述曲线拟合CF包括四次曲线拟合CF1、五次曲线拟合CF2、六次曲线拟合CF3、对数函数拟合CF4、指数函数拟合CF5和幂函数拟合CF6,分别得到拟合曲线Y1、拟合曲线Y2、拟合曲线Y3、拟合曲线Y4、拟合曲线Y5和拟合曲线Y6,从中选取拟合优度最接近1的作为最优拟合曲线Yo
6.根据权利要求4所述的一种机动车尾气标准的检测设备,其特征在于:所述数据分析单元U3用于将不透光烟度L与参比数据Ref进行对比:
L≤Ref时,该机动车尾气标准合格;
L>Ref时,该机动车尾气标准不合格。
7.根据权利要求4所述的一种机动车尾气标准的检测设备,其特征在于,所述结果显示单元U4给出如下指示信息:
机动车尾气标准的判断结果为合格时,通过文字W1,和/或颜色C1,和/或语音V1给出合格指示信息;
接收机动车尾气标准的判断结果为不合格时,通过文字W2,和/或颜色C2,和/或语音V2给出不合格指示信息。
8.根据权利要求4~7中任一项所述的机动车尾气标准的检测设备,其特征在于,该检测设备为便携式一体机或现场检测装置。
CN201710828659.XA 2017-09-14 2017-09-14 机动车尾气标准的检测方法及设备 Active CN107607492B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710828659.XA CN107607492B (zh) 2017-09-14 2017-09-14 机动车尾气标准的检测方法及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710828659.XA CN107607492B (zh) 2017-09-14 2017-09-14 机动车尾气标准的检测方法及设备

Publications (2)

Publication Number Publication Date
CN107607492A CN107607492A (zh) 2018-01-19
CN107607492B true CN107607492B (zh) 2020-10-27

Family

ID=61063545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710828659.XA Active CN107607492B (zh) 2017-09-14 2017-09-14 机动车尾气标准的检测方法及设备

Country Status (1)

Country Link
CN (1) CN107607492B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106769732A (zh) * 2016-12-31 2017-05-31 中国科学技术大学 一种垂直式柴油车烟度检测方法
CN109447010A (zh) * 2018-11-02 2019-03-08 天津同阳科技发展有限公司 基于过程数据的车辆分析装置和方法
CN110441481A (zh) * 2019-08-19 2019-11-12 郑州今迈衡器有限公司 一种用于监测机动车尾气的监测车

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB918525A (en) * 1958-04-29 1963-02-13 Hartridge Ltd Leslie Improvements in and relating to smokemeters
US4544273A (en) * 1983-07-29 1985-10-01 Particulate Instruments Smoke opacity meter
KR20140094727A (ko) * 2013-01-21 2014-07-31 박정환 적외선 센서를 이용한 굴뚝 배기 가스의 연기 농도 측정기
CN106680281B (zh) * 2016-12-31 2018-07-06 中国科学技术大学 一种基于深度残差学习网络的柴油车尾气烟度检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
两种烟度计在柴油机排气净化装置检测中的应用研究;郭鹏等;《内燃机》;20080229(第1期);第31-33页 *

Also Published As

Publication number Publication date
CN107607492A (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107607492B (zh) 机动车尾气标准的检测方法及设备
CN103632158B (zh) 森林防火监控方法及森林防火监控系统
JP6867076B1 (ja) 室内汚染と室外汚染を組み合わせた条件下での室内ホルムアルデヒド濃度の予測方法
CN112070393B (zh) 一种工农业生产排放的废水水质检测系统
CN109799195B (zh) 一种激光诱导击穿光谱的高精度定量分析方法
US6623975B1 (en) Method and system for vehicle emission testing
CN115908432B (zh) 一种物料输出品质检测系统及预测方法
CN115600919B (zh) 用于对实时无组织排放定位与园区排放总量计算的方法
CN115631138A (zh) 一种锆合金板材激光切割质量监检测方法与装置
CN112367518A (zh) 一种输电线路无人机巡检图像质量评价方法
CN112085482B (zh) 一种基于大数据的工程项目监理质量监测管理系统
CN115062675A (zh) 一种基于神经网络的全光谱污染溯源方法及云系统
CN116665057A (zh) 一种基于图像处理的河道水质监测方法及系统
CN113720914B (zh) 超声波探伤系统及超声波探伤方法
CN113884492A (zh) 一种用于机动车尾气的林格曼黑度校准和检测方法及装置
CN110134040A (zh) 工业设备的运行数据的处理方法及系统
CN116466368B (zh) 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法
CN113267601A (zh) 一种基于机器视觉和数据分析的工业生产环境远程实时监测云平台
CN116808922A (zh) 一种指定调胶设备稳定可靠型高效精准加料系统
CN115841478A (zh) 应用于车载空调软管生产管控的质量检测系统
CN114965860B (zh) 基于阻燃板阻燃测试的作业环境火灾消防安全监测系统
CN206146647U (zh) 一种汽车尾气流量测量控制系统
CN115060323A (zh) 一种智慧城市环境影响评估装置及方法
CN111307753B (zh) 一种基于太赫兹技术的纤维材料检测系统及方法
CN114720425A (zh) 一种基于图像识别的能见度监测系统及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant