CN106707850A - 一种太阳能环境监测车 - Google Patents
一种太阳能环境监测车 Download PDFInfo
- Publication number
- CN106707850A CN106707850A CN201611111086.0A CN201611111086A CN106707850A CN 106707850 A CN106707850 A CN 106707850A CN 201611111086 A CN201611111086 A CN 201611111086A CN 106707850 A CN106707850 A CN 106707850A
- Authority
- CN
- China
- Prior art keywords
- sensor
- data
- module
- registration
- monitoring device
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 52
- 238000012806 monitoring device Methods 0.000 claims abstract description 45
- 238000012545 processing Methods 0.000 claims abstract description 31
- 239000002689 soil Substances 0.000 claims abstract description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 14
- 238000007405 data analysis Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims description 45
- 238000000034 method Methods 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000008569 process Effects 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 20
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 16
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 claims description 16
- 238000013213 extrapolation Methods 0.000 claims description 12
- 230000009897 systematic effect Effects 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 9
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 claims description 8
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 8
- 239000001569 carbon dioxide Substances 0.000 claims description 8
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 claims description 8
- 229910052799 carbon Inorganic materials 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 229910001385 heavy metal Inorganic materials 0.000 claims description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 5
- 229910052760 oxygen Inorganic materials 0.000 claims description 5
- 239000001301 oxygen Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 239000002957 persistent organic pollutant Substances 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000007613 environmental effect Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000005265 energy consumption Methods 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种太阳能环境监测车,包括监测车本体、环境监测装置和中央处理器,监测车本体的顶部固定安装有电动机;电动机的顶部通过转动轴与旋转平台连接;旋转平台上固定安装有气泵;气泵通过气动伸缩杆与太阳能电池板的另一端活动连接;环境监测装置由大气监测装置、水质监测装置和土壤监测装置构成;中央处理器的输入端分别与太阳能电池板、输入模块、光强传感器和GPS定位模块的输出端电性连接;中央处理器的输出端通过采集数据分析模块与数据处理模块的输入端电性连接,该发明的有益效果是监测参数多样化,能对大气、水质和土壤的基本参数进行监测、分析和处理,太阳能电池板角度可随意调节。
Description
技术领域
本发明属于环境监测技术领域,尤其涉及一种太阳能环境监测车。
背景技术
随着经济社会的发展,环境质量日益受到国家和人民的重视,针对环境事件进行快速监测成为了迫切需求,移动式环境监测成为了一个必要的手段。可以进行现场采样、分析、储存样品的移动式环境监测车成为了一个较好的解决方案。
目前存在的移动式环境监测车,在作业过程中,大部分使用车载柴油发电机作为电源供应,因此,不可避免的会对周围环境排放各类污染物,对环境监测结果数据准确性造成很大影响。尤其是对一些大气环境监测实验,车载发电机作为电源供应的环境监测车所造成的有机物、颗粒物、重金属污染对实验结果影响较大,从而导致了实验分析数据偏离结果,造成实验分析失误。少数的采用太阳能电池板进行发电,但是太阳能电池板的可活动角度为180°,不能进行全方位全角度的调节。另一方面现有的太阳能环境监测车器,监测手段和监测设备单一,不能对周围环境进行多方面的监测,对污染不能提供有效措施。
还有就是针对海量数据的采集与分析的研究有很多,例如,对所采集的能耗数据进行预处理或执行简单固定的采集策略等,目前环境监测系统中的采集终端大多采用固定简单的模式进行能耗数据采集,会产生大量重复或者无效的数据,这对数据的存储、挖掘和分析造成额外负担,降低了整个系统的运行效率,另一方面就是本地数据的存储具有一定的局限性。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种监测参数多样化,能对大气、水质和土壤的基本参数进行监测、分析和处理,太阳能电池板角度可随意调节的一种太阳能环境监测车。
本发明是这样实现的,一种太阳能环境监测车,包括监测车本体、环境监测装置和中央处理器,所述监测车本体的顶部固定安装有电动机;所述电动机的顶部通过转动轴与旋转平台连接;所述旋转平台上固定安装有气泵;所述气泵通过气动伸缩杆与太阳能电池板的另一端活动连接;所述环境监测装置由大气监测装置、水质监测装置和土壤监测装置构成;所述环境监测装置的输出端与采集数据控制模块的输入端电性连接;所述采集数据控制模块的输出端与输入模块的输入端电性连接;所述输入模块的输入端通过CCD图像传感器与CCD摄像机的输出端电性连接;所述中央处理器的输入端分别与太阳能电池板、输入模块、光强传感器和GPS定位模块的输出端电性连接;所述中央处理器的输出端分别与电动机、气泵和调光模块的输入端电性连接;所述中央处理器的输出端通过采集数据分析模块与数据处理模块的输入端电性连接;所述中央处理器分别与操作显示屏、RAM存储器、MRAM存储器、数据库和无线射频收发模块电性连接;所述无线射频收发模块通过GPRS网络分别与云端存储器和监控中心连接;第一摄像机、第二摄像机分别安装在监测车本体的内部和外部;与中央处理器有线连接;
所述旋转平台的另一侧通过合页与太阳能电池板活动连接;
所述大气监测装置由第一温度传感器、第一湿度传感器、PM2.5传感器、二氧化碳传感器、二氧化氮传感器、二氧化硫传感器、有机污染物监测仪构成;
所述水质监测装置由电导率传感器、pH传感器、浊度传感器和溶解氧传感器构成;
所述土壤监测装置由第二温度传感器、第二湿度传感器和重金属离子监测仪构成。
进一步,所述环境监测装置设置有传感器空间配准模块;所述传感器空间配准模块与第一温度传感器、第一湿度传感器、PM2.5传感器、二氧化碳传感器、二氧化氮传感器、二氧化硫传感器、第二温度传感器、第二湿度传感器有线连接;所述传感器空间配准模块的传感器空间配准的方法包括时间对准过程和传感器系统误差估计过程,且伪量测方程的建立过程仅与目标的位置相关,而与目标的速度等状态无关;
时间对准过程完成传感器数据之间在时间上的对准,传感器A、传感器B在本地直角坐标系下的量测数据分别为YA(ti)和YB(ti),且传感器A的采样频率大于传感器B的采样频率,则由传感器A向传感器B的采样时刻进行配准,具体为:
采用内插外推的时间配准算法将传感器A的采样数据向传感器B的数据进行配准,使得两个传感器在空间配准时刻对同一个目标有同步的量测数据,内插外推时间配准算法如下:
在同一时间片内将各传感器观测数据按测量精度进行增量排序,然后将传感器A的观测数据分别向传感器B的时间点内插、外推,以形成一系列等间隔的目标观测数据,采用常用的三点抛物线插值法的进行内插外推时间配准算法得传感器A在tBk时刻在本地直角坐标系下的量测值为:
其中,tBk为配准时刻,tk-1,tk,tk+1为传感器A距离配准时刻最近的三个采样时刻,YA(tk-1),YA(tk),YA(tk+1)分别为其对应的对目标的探测数据;
完成时间配准后,根据传感器A的配准数据与传感器B的采样数据,采用基于地心地固坐标系下的伪量测法实现传感器A和传感器B的系统误差的估计;基于ECEF的系统误差估计算法具体为:
k时刻目标在本地直角坐标系下真实位置为X'1(k)=[x'1(k),y'1(k),z'1(k)]T,极坐标系下对应的量测值为分别为距离、方位角、俯仰角;转换至本地直角坐标系下为X1(k)=[x1(k),y1(k),z1(k)]T;传感器系统偏差为分别为距离、方位角和俯仰角的系统误差;于是有:
其中表示观测噪声,均值为零、方差为
式(1)用一阶近似展开并写成矩阵形式为:
X'1(k)=X1(k)+C(k)[ξ(k)+n(k)]\*MERGEFORMAT(3)
其中,
设两部传感器A和B,则对于同一个公共目标,地心地固坐标系下为X'e=[x'e,y'e,z'e]T,可得:
X'e=XAs+BAX'A1(k)=XBs+BBX'B1(k)\*MERGEFORMAT(4)
BA,BB分别为目标在传感器A与传感器B本地坐标下的位置转换到ECEF坐标系下的位置时的转换矩阵;
定义伪量测为:
Z(k)=XAe(k)-XBe(k)\*MERGEFORMAT(5)
其中,XAe(k)=XAs+BAXA1(k);XBe(k)=XBs+BBXB1(k)
将式(2)、式(3)代入式(4)可以得到关于传感器偏差的伪测量方程:
Z(k)=H(k)β(k)+W(k)\*MERGEFORMAT(6)
其中,Z(k)为伪测量向量;H(k)为测量矩阵;β为传感器偏差向量;W(k)为测量噪声向量;由于nA(k),nB(k)为零均值、相互独立的高斯型随机变量,因此W(k)同样是零均值高斯型随机变量,其协方差矩阵为R(k)。
进一步,所述传感器的空间配准方法包括的步骤如下:
步骤一、确定低采样率传感器的数据时刻;
步骤二、采用内插外推时间配准算法,根据传感器A向传感器B进行配准的计算方法,计算高采样率传感器的时间配准结果;
步骤三、根据伪量测构建过程计算伪量测方程;
步骤四、根据所构建的状态方程与伪量测方程,采用卡尔曼滤波,估计传感器的系统误差。
进一步,传感器A的量测模型如下:
YA(tk-1)、YA(tk)、YA(tk+1)分别为传感器A对目标在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的量测值,分别为:
其中,Y'A(tk-1)、Y'A(tk)、Y'A(tk+1)分别为传感器A在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的真实位置;CA(t)为误差的变换矩阵;ξA(t)为传感器的系统误差;为系统噪声,假设为零均值、相互独立的高斯型随机变量,噪声协方差矩阵分别为RA(k-1)、RA(k)、RA(k+1)。
进一步,传感器A向传感器B进行配准的具体过程如下:
将式(7)、式(8)、式(9)带入式(1),可得:
其中:为传感器A的本地直角坐标系下目标的真实位置在tBk时刻的时间配准值;为系统误差造成的误差项;为随机噪声,假定tk-1、tk、tk+1时刻的噪声互不相关的零均值白噪声,则为均值为零,协方差矩阵为RA=a2RA(k-1)+b2RA(k)+c2RA(k+1)的白噪声,而a、b、c、分别为且a+b+c=1。
进一步,伪量测构建过程仅与目标的位置相关,而与目标的速度等状态无关,具体过程如下:
tBk时刻,传感器A的配准量测为传感器B的量测为:
Y'B(tBk)为传感器B的本地直角坐标系下目标的真实位置,ξB(tBk)为传感器B的系统误差,为传感器B的随机误差;
根据权利1中所述的基于ECEF的空间配准算法,定义伪量测为:
对于同一公共目标,在ECEF坐标系下的位置为Xe则有:
XAS(t)、XBS(t)分别为传感器A和传感器B在ECEF坐标系下的位置,
将式:
Z(k)=G(k)β(k)+W(k)\*MERGEFORMAT(15)
其中G(k)=[-aJA(tk-1) -bJA(tk) -cJA(tk+1) JB(t)]=λJ为量测矩阵,其中,λ=[-a -b -c 1],J=[JA(tk-1) JA(tk) JA(tk+1) JB(t)]T,JA(tk-1)=BA(tk-1)CA(tk-1),JA(tk)=BA(tk)CA(tk);JA(tk+1)=BA(tk+1)CA(tk+1),JB(tk+1)=BB(t)CB(t);β(k)=[ξA(tk-1)ξA(tk)ξA(tk+1)ξB(t)]T,为系统误差;为均值为零,协方差矩阵为:
RW(k)=a2RA(k-1)+b2RA(k)+c2RA(k+1)-RB(k)\*MERGEFORMAT(16)
符号T表示矩阵的转置运算。
进一步,所构建的状态方程与伪量测方程如下:
其中F(k+1|k)为状态方程的转移矩阵,取值与传感器的系统误差的变化规律相关,若传感器的系统误差是缓变的,则F(k+1|k)近似为单位矩阵,取为F(k+1|k)=0.99I,I为单位阵。
进一步,所述采集数据分析模块设置有收敛加密单元,所述收敛加密单元的收敛加密方法包括:
(1)数据拥有者首先计算数据m的哈希值h(m),其中h(·)是密码学中的一个强哈希函数,h(m)作为加密数据m的密钥;
(2)用h(m)加密数据,假设E是一个对称密钥加密函数,则产生的密文就是Eh(m)(m);
(3)用所有的授权用户的公钥加密h(m),密文是:(C,C′),其中C=Eh(m)(m),C′=FPK(h(m)),F是一个公钥加密函数,PK是公钥;
(4)数据只能被授权用户解密,合法用户首先用自己的私钥解密得到h(m),最后用h(m)来恢复m,表示如下:
本发明提供的太阳能环境监测车,采用太阳能电池板进行发电提供电源,节能环保;太阳能电池板通过电动机和气泵可进行全方位的角度调整,最大程度的接收太阳光,利用率高;监测参数多样化,可以周围大气、水质和土壤进行监测、分析和处理,采集的数据通过采集数据控制模块、中央处理器和采集数据分析模块后,将被数据处理模块进一步的处理,采集与处理迅速,且不会出现混乱和遗漏的问题,保证了环境监测参数的采集与处理能稳定同步进行,利用RAM存储器、MRAM存储器和数据库,可对监测数据的信息进行对比,采样和查询,更好地了解环境污染问题,监测的数据通过无线射频收发模块和GPRS网络存储到云端存储器中,保证数据的安全存储。
附图说明
图1是本发明实施例提供的太阳能环境监测车的原理框图。
图2是本发明实施例提供的太阳能电池板安装示意图。
图中:1、监测车本体;2、电动机;3、转动轴;4、旋转平台;5、气泵;6、合页;7、太阳能电池板;8、气动伸缩杆;9、环境监测装置;10、大气监测装置;11、水质监测装置;12、土壤监测装置;13、第一温度传感器;14、第一湿度传感器;15、PM2.5传感器;16、二氧化碳传感器;17、二氧化氮传感器;18、二氧化硫传感器;19、电导率传感器;20、PH传感器;21、浊度传感器;22、溶解氧传感器;23、第二温度传感器;24、第二湿度传感器;25、重金属离子监测仪;26、采集数据控制模块;27、输入模块;28、CCD图像传感器;29、CCD摄像机;30、中央处理器;31、光强传感器;32、GPS定位模块;33、调光模块;34、采集数据分析模块;35、数据处理模块;36、操作显示屏;37、RAM存储器;38、MRAM存储器;39、数据库;40、无线射频收发模块;41、GPRS网络;42、云端存储器;43、监控中心;44、第一摄像机;45、第二摄像机。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合图1和图2对本发明的结构作详细的描述:
本发明实施例提供的太阳能环境监测车,包括监测车本体1、环境监测装置9和中央处理器30,所述监测车本体1的顶部固定安装有电动机2;所述电动机2的顶部通过转动轴3与旋转平台4连接;所述旋转平台4上固定安装有气泵5;所述气泵5通过气动伸缩杆8与太阳能电池板7的另一端活动连接;所述环境监测装置9由大气监测装置10、水质监测装置11和土壤监测装置12构成;所述环境监测装置9的输出端与采集数据控制模块26的输入端电性连接;所述采集数据控制模块26的输出端与输入模块27的输入端电性连接;所述输入模块27的输入端通过CCD图像传感器28与CCD摄像机29的输出端电性连接;所述中央处理器30的输入端分别与太阳能电池板7、输入模块27、光强传感器31和GPS定位模块32的输出端电性连接;所述中央处理器30的输出端分别与电动机2、气泵5和调光模块33的输入端电性连接;所述中央处理器30的输出端通过采集数据分析模块34与数据处理模块35的输入端电性连接;所述中央处理器30分别与操作显示屏36、RAM存储器37、MRAM存储器38、数据库39和无线射频收发模块40电性连接;所述无线射频收发模块40通过GPRS网络41分别与云端存储器42和监控中心43连接;第一摄像机44、第二摄像机45分别安装在监测车本体1的内部和外部;与中央处理器30有线连接。
进一步,所述旋转平台4的另一侧通过合页6与太阳能电池板7活动连接。
进一步,所述大气监测装置10由第一温度传感器13、第一湿度传感器14、PM2.5传感器15、二氧化碳传感器16、二氧化氮传感器17、二氧化硫传感器18、有机污染物监测仪构成。
进一步,所述水质监测装置11由电导率传感器19、PH传感器20、浊度传感器21和溶解氧传感器22构成。
进一步,所述土壤监测装置12由第二温度传感器23、第二湿度传感器24和重金属离子监测仪25构成。
进一步,所述环境监测装置设置有传感器空间配准模块;所述传感器空间配准模块与第一温度传感器、第一湿度传感器、PM2.5传感器、二氧化碳传感器、二氧化氮传感器、二氧化硫传感器、第二温度传感器、第二湿度传感器有线连接;所述传感器空间配准模块的传感器空间配准的方法包括时间对准过程和传感器系统误差估计过程,且伪量测方程的建立过程仅与目标的位置相关,而与目标的速度等状态无关;
时间对准过程完成传感器数据之间在时间上的对准,传感器A、传感器B在本地直角坐标系下的量测数据分别为YA(ti)和YB(ti),且传感器A的采样频率大于传感器B的采样频率,则由传感器A向传感器B的采样时刻进行配准,具体为:
采用内插外推的时间配准算法将传感器A的采样数据向传感器B的数据进行配准,使得两个传感器在空间配准时刻对同一个目标有同步的量测数据,内插外推时间配准算法如下:
在同一时间片内将各传感器观测数据按测量精度进行增量排序,然后将传感器A的观测数据分别向传感器B的时间点内插、外推,以形成一系列等间隔的目标观测数据,采用常用的三点抛物线插值法的进行内插外推时间配准算法得传感器A在tBk时刻在本地直角坐标系下的量测值为:
其中,tBk为配准时刻,tk-1,tk,tk+1为传感器A距离配准时刻最近的三个采样时刻,YA(tk-1),YA(tk),YA(tk+1)分别为其对应的对目标的探测数据;
完成时间配准后,根据传感器A的配准数据与传感器B的采样数据,采用基于地心地固坐标系下的伪量测法实现传感器A和传感器B的系统误差的估计;基于ECEF的系统误差估计算法具体为:
k时刻目标在本地直角坐标系下真实位置为X'1(k)=[x'1(k),y'1(k),z'1(k)]T,极坐标系下对应的量测值为分别为距离、方位角、俯仰角;转换至本地直角坐标系下为X1(k)=[x1(k),y1(k),z1(k)]T;传感器系统偏差为分别为距离、方位角和俯仰角的系统误差;于是有:
其中表示观测噪声,均值为零、方差为
式(1)用一阶近似展开并写成矩阵形式为:
X'1(k)=X1(k)+C(k)[ξ(k)+n(k)]\*MERGEFORMAT(3)
其中,
设两部传感器A和B,则对于同一个公共目标,地心地固坐标系下为X'e=[x'e,y'e,z'e]T,可得:
X'e=XAs+BAX'A1(k)=XBs+BBX'B1(k)\*MERGEFORMAT(4)
BA,BB分别为目标在传感器A与传感器B本地坐标下的位置转换到ECEF坐标系下的位置时的转换矩阵;
定义伪量测为:
Z(k)=XAe(k)-XBe(k)\*MERGEFORMAT(5)
其中,XAe(k)=XAs+BAXA1(k);XBe(k)=XBs+BBXB1(k)
将式(2)、式(3)代入式(4)可以得到关于传感器偏差的伪测量方程:
Z(k)=H(k)β(k)+W(k)\*MERGEFORMAT(6)
其中,Z(k)为伪测量向量;H(k)为测量矩阵;β为传感器偏差向量;W(k)为测量噪声向量;由于nA(k),nB(k)为零均值、相互独立的高斯型随机变量,因此W(k)同样是零均值高斯型随机变量,其协方差矩阵为R(k)。
进一步,所述传感器的空间配准方法包括的步骤如下:
步骤一、确定低采样率传感器的数据时刻;
步骤二、采用内插外推时间配准算法,根据传感器A向传感器B进行配准的计算方法,计算高采样率传感器的时间配准结果;
步骤三、根据伪量测构建过程计算伪量测方程;
步骤四、根据所构建的状态方程与伪量测方程,采用卡尔曼滤波,估计传感器的系统误差。
进一步,传感器A的量测模型如下:
YA(tk-1)、YA(tk)、YA(tk+1)分别为传感器A对目标在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的量测值,分别为:
其中,Y'A(tk-1)、Y'A(tk)、Y'A(tk+1)分别为传感器A在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的真实位置;CA(t)为误差的变换矩阵;ξA(t)为传感器的系统误差;为系统噪声,假设为零均值、相互独立的高斯型随机变量,噪声协方差矩阵分别为RA(k-1)、RA(k)、RA(k+1)。
进一步,传感器A向传感器B进行配准的具体过程如下:
将式(7)、式(8)、式(9)带入式(1),可得:
其中:为传感器A的本地直角坐标系下目标的真实位置在tBk时刻的时间配准值;为系统误差造成的误差项;为随机噪声,假定tk-1、tk、tk+1时刻的噪声互不相关的零均值白噪声,则为均值为零,协方差矩阵为RA=a2RA(k-1)+b2RA(k)+c2RA(k+1)的白噪声,而a、b、c、分别为且a+b+c=1。
进一步,伪量测构建过程仅与目标的位置相关,而与目标的速度等状态无关,具体过程如下:
tBk时刻,传感器A的配准量测为传感器B的量测为:
Y'B(tBk)为传感器B的本地直角坐标系下目标的真实位置,ξB(tBk)为传感器B的系统误差,为传感器B的随机误差;
根据权利1中所述的基于ECEF的空间配准算法,定义伪量测为:
对于同一公共目标,在ECEF坐标系下的位置为Xe则有:
XAS(t)、XBS(t)分别为传感器A和传感器B在ECEF坐标系下的位置,
将式:
Z(k)=G(k)β(k)+W(k)\*MERGEFORMAT(15)
其中G(k)=[-aJA(tk-1) -bJA(tk) -cJA(tk+1) JB(t)]=λJ为量测矩阵,其中,λ=[-a -b -c 1],J=[JA(tk-1) JA(tk) JA(tk+1) JB(t)]T,JA(tk-1)=BA(tk-1)CA(tk-1),JA(tk)=BA(tk)CA(tk);JA(tk+1)=BA(tk+1)CA(tk+1),JB(tk+1)=BB(t)CB(t);β(k)=[ξA(tk-1) ξA(tk) ξA(tk+1) ξB(t)]T,为系统误差;为均值为零,协方差矩阵为:
RW(k)=a2RA(k-1)+b2RA(k)+c2RA(k+1)-RB(k)
\*MERGEFORMAT(16)
符号T表示矩阵的转置运算。
进一步,所构建的状态方程与伪量测方程如下:
其中F(k+1|k)为状态方程的转移矩阵,取值与传感器的系统误差的变化规律相关,若传感器的系统误差是缓变的,则F(k+1|k)近似为单位矩阵,取为F(k+1|k)=0.99I,I为单位阵。
进一步,所述采集数据分析模块设置有收敛加密单元,所述收敛加密单元的收敛加密方法包括:
(1)数据拥有者首先计算数据m的哈希值h(m),其中h(·)是密码学中的一个强哈希函数,h(m)作为加密数据m的密钥;
(2)用h(m)加密数据,假设E是一个对称密钥加密函数,则产生的密文就是Eh(m)(m);
(3)用所有的授权用户的公钥加密h(m),密文是:(C,C′),其中C=Eh(m)(m),C′=FPK(h(m)),F是一个公钥加密函数,PK是公钥;
(4)数据只能被授权用户解密,合法用户首先用自己的私钥解密得到h(m),最后用h(m)来恢复m,表示如下:
工作原理:该太阳能环境监测车,采用太阳能电池板7进行发电提供电源,节能环保,太阳能电池板7通过电动机2、转动轴3和旋转平台4可进行旋转,通过气泵5和气动伸缩杆8可进行倾角的调节,以此对太阳能电池板7进行全方位的角度调整,最大程度的接收太阳光,利用率高,监测参数多样化,设置有大气监测装置10、水质监测装置11和土壤监测装置12,大气监测装置10由第一温度传感器13、第一湿度传感器14、PM2.5传感器15、二氧化碳传感器16、二氧化氮传感器17和二氧化硫传感器18构成,水质监测装置11由电导率传感器19、PH传感器20、浊度传感器21和溶解氧传感器22,土壤监测装置12由第二温度传感器23、第二湿度传感器24和重金属离子监测仪25构成,可以对周围大气、水质和土壤进行全方面的监测,监测的数据通过采集数据控制模块26、输入模块27、中央处理器30和采集数据分析模块34处理后,将会被数据处理模块35进行进一步的处理,RAM存储器37、MRAM存储器38和数据库39的配合使用,可对监测数据的信息进行对比、采样和查询,在保证了监测数据的采集与处理的快速性的同时更好地了解环境污染问题,监测的数据通过无线射频收发模块40和GPRS网络41存储到云端存储器42中,保证数据的安全存储,监测的数据可实时通过无线射频收发模块40和GPRS网络41发送到监控中心43,无线射频收发模块40用于接收和发送无线网络信号,CCD摄像机29将周围的画面进行拍摄并通过CCD图像传感器28发送到输入模块27中,光强传感器31用于监测周围环境的光照强度,当光照强度满足不了需求时,中央处理器30控制调光模块33进行补光,操作显示屏36用于工作人员进行操控和实时显示监测参数,GPS定位模块32可对监测地点进行定位。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种太阳能环境监测车,包括监测车本体、环境监测装置和中央处理器,其特征在于,所述监测车本体的顶部固定安装有电动机;所述电动机的顶部通过转动轴与旋转平台连接;所述旋转平台上固定安装有气泵;所述气泵通过气动伸缩杆与太阳能电池板的另一端活动连接;所述环境监测装置由大气监测装置、水质监测装置和土壤监测装置构成;所述环境监测装置的输出端与采集数据控制模块的输入端电性连接;所述采集数据控制模块的输出端与输入模块的输入端电性连接;所述输入模块的输入端通过CCD图像传感器与CCD摄像机的输出端电性连接;所述中央处理器的输入端分别与太阳能电池板、输入模块、光强传感器和GPS定位模块的输出端电性连接;所述中央处理器的输出端分别与电动机、气泵和调光模块的输入端电性连接;所述中央处理器的输出端通过采集数据分析模块与数据处理模块的输入端电性连接;所述中央处理器分别与操作显示屏、RAM存储器、MRAM存储器、数据库和无线射频收发模块电性连接;所述无线射频收发模块通过GPRS网络分别与云端存储器和监控中心连接;第一摄像机、第二摄像机分别安装在监测车本体的内部和外部;与中央处理器有线连接;
所述旋转平台的另一侧通过合页与太阳能电池板活动连接;
所述大气监测装置由第一温度传感器、第一湿度传感器、PM2.5传感器、二氧化碳传感器、二氧化氮传感器、二氧化硫传感器、有机污染物监测仪构成;
所述水质监测装置由电导率传感器、pH传感器、浊度传感器和溶解氧传感器构成;
所述土壤监测装置由第二温度传感器、第二湿度传感器和重金属离子监测仪构成。
2.如权利要求1所述的太阳能环境监测车,其特征在于,所述环境监测装置设置有传感器空间配准模块;所述传感器空间配准模块与第一温度传感器、第一湿度传感器、PM2.5传感器、二氧化碳传感器、二氧化氮传感器、二氧化硫传感器、第二温度传感器、第二湿度传感器有线连接;所述传感器空间配准模块的传感器空间配准的方法包括时间对准过程和传感器系统误差估计过程,且伪量测方程的建立过程仅与目标的位置相关,而与目标的速度状态无关;
时间对准过程完成传感器数据之间在时间上的对准,传感器A、传感器B在本地直角坐标系下的量测数据分别为YA(ti)和YB(ti),且传感器A的采样频率大于传感器B的采样频率,则由传感器A向传感器B的采样时刻进行配准,具体为:
采用内插外推的时间配准算法将传感器A的采样数据向传感器B的数据进行配准,使得两个传感器在空间配准时刻对同一个目标有同步的量测数据,内插外推时间配准算法如下:
在同一时间片内将各传感器观测数据按测量精度进行增量排序,然后将传感器A的观测数据分别向传感器B的时间点内插、外推,以形成一系列等间隔的目标观测数据,采用常用的三点抛物线插值法的进行内插外推时间配准算法得传感器A在tBk时刻在本地直角坐标系下的量测值为:
其中,tBk为配准时刻,tk-1,tk,tk+1为传感器A距离配准时刻最近的三个采样时刻,YA(tk-1),YA(tk),YA(tk+1)分别为其对应的对目标的探测数据;
完成时间配准后,根据传感器A的配准数据与传感器B的采样数据,采用基于地心地固坐标系下的伪量测法实现传感器A和传感器B的系统误差的估计;基于ECEF的系统误差估计算法具体为:
k时刻目标在本地直角坐标系下真实位置为X'1(k)=[x'1(k),y'1(k),z'1(k)]T,极坐标系下对应的量测值为分别为距离、方位角、俯仰角;转换至本地直角坐标系下为X1(k)=[x1(k),y1(k),z1(k)]T;传感器系统偏差为分别为距离、方位角和俯仰角的系统误差;于是有:
其中表示观测噪声,均值为零、方差为
式(1)用一阶近似展开并写成矩阵形式为:
X'1(k)=X1(k)+C(k)[ξ(k)+n(k)]\*MERGEFORMAT (3)
其中,
设两部传感器A和B,则对于同一个公共目标,地心地固坐标系下为X'e=[x'e,y'e,z'e]T,可得:
X'e=XAs+BAX'A1(k)=XBs+BBX'B1(k)\*MERGEFORMAT (4)
BA,BB分别为目标在传感器A与传感器B本地坐标下的位置转换到ECEF坐标系下的位置时的转换矩阵;
定义伪量测为:
Z(k)=XAe(k)-XBe(k)\*MERGEFORMAT (5)
其中,XAe(k)=XAs+BAXA1(k);XBe(k)=XBs+BBXB1(k)
将式(2)、式(3)代入式(4)得到关于传感器偏差的伪测量方程:
Z(k)=H(k)β(k)+W(k)\*MERGEFORMAT (6)
其中,Z(k)为伪测量向量;H(k)为测量矩阵;β为传感器偏差向量;W(k)为测量噪声向量;由于nA(k),nB(k)为零均值、相互独立的高斯型随机变量,因此W(k)同样是零均值高斯型随机变量,其协方差矩阵为R(k)。
3.如权利要求2所述的太阳能环境监测车,其特征在于,所述传感器的空间配准方法包括的步骤如下:
步骤一、确定低采样率传感器的数据时刻;
步骤二、采用内插外推时间配准算法,根据传感器A向传感器B进行配准的计算方法,计算高采样率传感器的时间配准结果;
步骤三、根据伪量测构建过程计算伪量测方程;
步骤四、根据所构建的状态方程与伪量测方程,采用卡尔曼滤波,估计传感器的系统误差。
4.如权利要求2所述的太阳能环境监测车,其特征在于,传感器A的量测模型如下:
YA(tk-1)、YA(tk)、YA(tk+1)分别为传感器A对目标在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的量测值,分别为:
其中,Y'A(tk-1)、Y'A(tk)、Y'A(tk+1)分别为传感器A在tk-1,tk,tk+1时刻的本地笛卡尔坐标系下的真实位置;CA(t)为误差的变换矩阵;ξA(t)为传感器的系统误差;为系统噪声,假设为零均值、相互独立的高斯型随机变量,噪声协方差矩阵分别为RA(k-1)、RA(k)、RA(k+1)。
5.如权利要求3所述的太阳能环境监测车,其特征在于,传感器A向传感器B进行配准的具体过程如下:
将式(7)、式(8)、式(9)带入式(1),可得:
其中:为传感器A的本地直角坐标系下目标的真实位置在tBk时刻的时间配准值;为系统误差造成的误差项;为随机噪声,假定tk-1、tk、tk+1时刻的噪声互不相关的零均值白噪声,则为均值为零,协方差矩阵为RA=a2RA(k-1)+b2RA(k)+c2RA(k+1)的白噪声,而a、b、c、分别为且a+b+c=1。
6.如权利要求3所述的太阳能环境监测车,其特征在于,伪量测构建过程仅与目标的位置相关,而与目标的速度状态无关,具体过程如下:
tBk时刻,传感器A的配准量测为传感器B的量测为:
Y'B(tBk)为传感器B的本地直角坐标系下目标的真实位置,ξB(tBk)为传感器B的系统误差,为传感器B的随机误差;
根据权利1中所述的基于ECEF的空间配准算法,定义伪量测为:
对于同一公共目标,在ECEF坐标系下的位置为Xe则有:
XAS(t)、XBS(t)分别为传感器A和传感器B在ECEF坐标系下的位置,
将式:
带入式(12),则有:
Z(k)=G(k)β(k)+W(k)\*MERGEFORMAT (15)
其中G(k)=[-aJA(tk-1) -bJA(tk) -cJA(tk+1) JB(t)]=λJ为量测矩阵,其中,λ=[-a -b-c 1],J=[JA(tk-1) JA(tk) JA(tk+1) JB(t)]T,JA(tk-1)=BA(tk-1)CA(tk-1),JA(tk)=BA(tk)CA(tk);JA(tk+1)=BA(tk+1)CA(tk+1),JB(tk+1)=BB(t)CB(t);β(k)=[ξA(tk-1)ξA(tk)ξA(tk+1)ξB(t)]T,为系统误差;为均值为零,协方差矩阵为:
RW(k)=a2RA(k-1)+b2RA(k)+c2RA(k+1)-RB(k)
\*MERGEFORMAT (16)
符号T表示矩阵的转置运算。
7.如权利要求3所述的太阳能环境监测车,其特征在于,所构建的状态方程与伪量测方程如下:
其中F(k+1|k)为状态方程的转移矩阵,取值与传感器的系统误差的变化规律相关,若传感器的系统误差是缓变的,则F(k+1|k)近似为单位矩阵,取为F(k+1|k)=0.99I,I为单位阵。
8.如权利要求1所述的太阳能环境监测车,其特征在于,所述采集数据分析模块设置有收敛加密单元,所述收敛加密单元的收敛加密方法包括:
(1)数据拥有者首先计算数据m的哈希值h(m),其中h(·)是密码学中的一个强哈希函数,h(m)作为加密数据m的密钥;
(2)用h(m)加密数据,假设E是一个对称密钥加密函数,则产生的密文就是Eh(m)(m);
(3)用所有的授权用户的公钥加密h(m),密文是:(C,C′),其中C=Eh(m)(m),C′=FPK(h(m)),F是一个公钥加密函数,PK是公钥;
(4)数据只能被授权用户解密,合法用户首先用自己的私钥解密得到h(m),最后用h(m)来恢复m,表示如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611111086.0A CN106707850A (zh) | 2016-12-06 | 2016-12-06 | 一种太阳能环境监测车 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611111086.0A CN106707850A (zh) | 2016-12-06 | 2016-12-06 | 一种太阳能环境监测车 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106707850A true CN106707850A (zh) | 2017-05-24 |
Family
ID=58937630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611111086.0A Pending CN106707850A (zh) | 2016-12-06 | 2016-12-06 | 一种太阳能环境监测车 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106707850A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107146501A (zh) * | 2017-06-29 | 2017-09-08 | 广东雅嘉思科技有限公司 | 一种多功能教学实验车 |
CN107184145A (zh) * | 2017-06-21 | 2017-09-22 | 武汉洁美雅科技有限公司 | 一种智能节能环保吸尘器控制系统 |
CN107588847A (zh) * | 2017-08-19 | 2018-01-16 | 安徽泛亚环保科技有限公司 | 一种基于物联网的噪音检测系统 |
CN107608009A (zh) * | 2017-09-15 | 2018-01-19 | 深圳市卡普瑞环境科技有限公司 | 一种大气质量监测设备、处理终端以及服务器 |
CN108049494A (zh) * | 2017-12-15 | 2018-05-18 | 南昌工程学院 | 一种混凝土结构地下室抗裂防渗施工方法 |
CN108801343A (zh) * | 2018-05-25 | 2018-11-13 | 江苏师范大学 | 一种应用于大气环境监测的智能小车系统 |
CN109238362A (zh) * | 2018-09-26 | 2019-01-18 | 佛山市二加环保科技有限公司 | 区域环境监测控制方法 |
CN117274884A (zh) * | 2023-11-21 | 2023-12-22 | 赣江新区慧工科技有限公司 | 一种基于图像识别的施工扬尘污染事件检测方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102955120A (zh) * | 2012-10-29 | 2013-03-06 | 北京华盛集智新能源科技有限公司 | 可调式全方位太阳能电池室外测试系统 |
CN103448599A (zh) * | 2013-08-22 | 2013-12-18 | 哈尔滨金大环境工程有限公司 | 太阳能环境监测车 |
CN104344978A (zh) * | 2014-11-11 | 2015-02-11 | 浙江诚德检测研究有限公司 | 一种环境监测样本采集器及采集方法 |
CN104809326A (zh) * | 2014-06-23 | 2015-07-29 | 方洋旺 | 一种异步传感器空间配准算法 |
CN205131043U (zh) * | 2015-10-28 | 2016-04-06 | 太原海纳辰科仪器仪表有限公司 | 一种多功能太阳能环境监测车 |
CN205283614U (zh) * | 2015-12-31 | 2016-06-01 | 天津海沃特科技发展有限公司 | 一种水环境无线通讯讯号智能监测装置 |
-
2016
- 2016-12-06 CN CN201611111086.0A patent/CN106707850A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102955120A (zh) * | 2012-10-29 | 2013-03-06 | 北京华盛集智新能源科技有限公司 | 可调式全方位太阳能电池室外测试系统 |
CN103448599A (zh) * | 2013-08-22 | 2013-12-18 | 哈尔滨金大环境工程有限公司 | 太阳能环境监测车 |
CN104809326A (zh) * | 2014-06-23 | 2015-07-29 | 方洋旺 | 一种异步传感器空间配准算法 |
CN104344978A (zh) * | 2014-11-11 | 2015-02-11 | 浙江诚德检测研究有限公司 | 一种环境监测样本采集器及采集方法 |
CN205131043U (zh) * | 2015-10-28 | 2016-04-06 | 太原海纳辰科仪器仪表有限公司 | 一种多功能太阳能环境监测车 |
CN205283614U (zh) * | 2015-12-31 | 2016-06-01 | 天津海沃特科技发展有限公司 | 一种水环境无线通讯讯号智能监测装置 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107184145A (zh) * | 2017-06-21 | 2017-09-22 | 武汉洁美雅科技有限公司 | 一种智能节能环保吸尘器控制系统 |
CN107146501A (zh) * | 2017-06-29 | 2017-09-08 | 广东雅嘉思科技有限公司 | 一种多功能教学实验车 |
CN107146501B (zh) * | 2017-06-29 | 2019-10-15 | 广东雅嘉思科技有限公司 | 一种多功能教学实验车 |
CN107588847A (zh) * | 2017-08-19 | 2018-01-16 | 安徽泛亚环保科技有限公司 | 一种基于物联网的噪音检测系统 |
CN107608009A (zh) * | 2017-09-15 | 2018-01-19 | 深圳市卡普瑞环境科技有限公司 | 一种大气质量监测设备、处理终端以及服务器 |
CN108049494A (zh) * | 2017-12-15 | 2018-05-18 | 南昌工程学院 | 一种混凝土结构地下室抗裂防渗施工方法 |
CN108801343A (zh) * | 2018-05-25 | 2018-11-13 | 江苏师范大学 | 一种应用于大气环境监测的智能小车系统 |
CN109238362A (zh) * | 2018-09-26 | 2019-01-18 | 佛山市二加环保科技有限公司 | 区域环境监测控制方法 |
CN117274884A (zh) * | 2023-11-21 | 2023-12-22 | 赣江新区慧工科技有限公司 | 一种基于图像识别的施工扬尘污染事件检测方法及系统 |
CN117274884B (zh) * | 2023-11-21 | 2024-02-20 | 赣江新区慧工科技有限公司 | 一种基于图像识别的施工扬尘污染事件检测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106707850A (zh) | 一种太阳能环境监测车 | |
Zuo et al. | Multi-phase fingerprint map based on interpolation for indoor localization using iBeacons | |
Pankratius et al. | Mobile crowd sensing in space weather monitoring: the mahali project | |
CN108919367B (zh) | 基于电流磁场的交流输电线路反演方法 | |
CN104936287A (zh) | 基于矩阵补全的传感网室内指纹定位方法 | |
CN103542868B (zh) | 基于角度和强度的车载激光点云噪点自动去除方法 | |
CN104407338B (zh) | 基于海洋二号卫星微波散射计的极地海冰识别方法 | |
CN111398980B (zh) | 一种机载LiDAR数据处理的方法及装置 | |
CN107255566A (zh) | 一种车载式车辆行驶跑偏检测装置和检测方法 | |
CN109327797A (zh) | 基于WiFi网络信号的移动机器人室内定位系统 | |
WO2018032471A1 (zh) | 智能化跟踪系统及方法 | |
CN113375721A (zh) | 一种基于物联网的环境监测系统 | |
CN107437149B (zh) | 一种光伏电站出力的确定方法及系统 | |
CN106160659A (zh) | 一种光伏电站区域定向故障诊断方法 | |
CN103743469A (zh) | 压电水声换能器声辐射模态测量方法及系统 | |
CN111504380A (zh) | 一种基于风光互补发电的地理环境信息监测系统 | |
CN111770482A (zh) | 一种适用于智能设备的移动测距方法及应用 | |
CN114973443A (zh) | 一种基于巡检机器人的复杂气体环境巡检系统及其方法 | |
CN106767698A (zh) | 基于北斗ⅱ姿态测量的输电线路杆塔倾斜度监测系统及监测方法 | |
CN107656905B (zh) | 一种利用误差传递的空气质量数据实时校准方法 | |
CN102830391B (zh) | 一种红外搜索与跟踪系统准确性指标计算方法 | |
CN109752714A (zh) | 一种旋转雷达微变监测数据处理方法和雷达系统 | |
CN109100769A (zh) | 车联蓝牙钥匙的定位方法及系统 | |
CN107421904A (zh) | 基于差分吸收激光小车的六氟化硫气体泄漏在线监测系统 | |
CN113311826A (zh) | 一种基于环形红外阵列的自动跟随系统及其工作方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170524 |
|
WD01 | Invention patent application deemed withdrawn after publication |