CN105043944B - 基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 - Google Patents
基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 Download PDFInfo
- Publication number
- CN105043944B CN105043944B CN201510374742.5A CN201510374742A CN105043944B CN 105043944 B CN105043944 B CN 105043944B CN 201510374742 A CN201510374742 A CN 201510374742A CN 105043944 B CN105043944 B CN 105043944B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- light
- photomultiplier
- mtd
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 100
- 238000001514 detection method Methods 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 73
- 239000013598 vector Substances 0.000 claims abstract description 59
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 230000010287 polarization Effects 0.000 claims description 56
- 239000011248 coating agent Substances 0.000 claims description 22
- 238000000576 coating method Methods 0.000 claims description 22
- 230000009471 action Effects 0.000 claims description 21
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical group Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 9
- 239000007888 film coating Substances 0.000 claims description 6
- 238000009501 film coating Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000012798 spherical particle Substances 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 abstract description 5
- 238000009776 industrial production Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 3
- 239000013618 particulate matter Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000000149 argon plasma sintering Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005250 beta ray Effects 0.000 description 1
- 206010006451 bronchitis Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 206010039083 rhinitis Diseases 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种基于日盲紫外光进行雾霾粒子检测的装置,包括用于发射紫外光的LED灯,沿紫外光传播方向依次设有空间滤波器、准直透镜、起偏器、λ/4波片、分子气体容器和光电探测器;起偏器、λ/4波片分别与步进电机控制器连接,计算机与步进电机控制器连接。本发明还公开了利用上述装置进行雾霾粒子检测的方法:首先,测量光电探测器的仪器矩阵A;随后,利用测得的仪器矩阵A,得偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout;最后,根据mie散射理论,得待测雾霾粒子的粒径。本发明方法克服暗盒测量空间的局限性,不受环境湿度和温度的影响,不受背景光干扰,可实现对雾霾的全天候开放式测量,同时降低设备成本,适用于工业生产。
Description
技术领域
本发明属于超快光学技术和非线性光学领域,具体涉及一种基于日盲紫外光进行雾霾粒子检测的装置;本发明还涉及利用该检测装置进行雾霾粒子检测的方法。
背景技术
近年来,在我国经济社会快速发展的同时,雾霾天气成为了突出的环境问题,特别是大中城市的上空难见“蓝天白云”,社会各界对加大雾霾治理力度的呼声日益强烈。2015年初的纪录片《穹顶之下》也深刻的告诉人们:雾霾治理已经刻不容缓。雾霾天气是一种大气污染状态,是对大气中各种悬浮颗粒物含量超标的笼统表述,其PM2.5(ParticulateMatter 2.5)被认为是造成雾霾天气的“元凶”。PM2.5是指空气中飘浮的直径小于等于2.5μm的颗粒物,是形成雾霾天气的主要原因之一。雾霾天气被公认为灾害性天气:一是影响通行安全,地面能见度降低交通事故多发;二是导致气象变异,区域气候反常;三是损害人体健康,PM2.5能被人体吸入并粘附在呼吸道和肺叶而引起鼻炎、支气管炎等病症。目前对TSP(Total Suspended Particulates),PM10,PM2.5的测量技术主要有:微量振荡天平法、滤膜称重法、压电晶体法、β射线吸收法以及光散射测量法等。这些测量技术都有各自的优缺点,没有一种方法能够同时满足检测设备在实用性、精确性、可靠性及便利性等方面的要求。
日盲紫外光散射传播具有全天候、无背景光干扰、宽视场接收等优点,同时也具有在线测量、采样周期短、受环境因素影响小、无需暗盒等优点。偏振是各矢量波的共有性质,紫外光也具有这种特性,紫外偏振光分为线偏振光、圆偏振光、椭圆偏振光和部分偏振光,它们的偏振状态不同,即Stokes矢量不同,更重要的是可以作为信息调制的载体。偏振状态因雾霾颗粒物对光的散射而改变,经过起偏器后的紫外光成为偏振紫外光在雾霾天气传播过程中,由于雾霾颗粒物形状、大小、性质各异,即粒子在不同散射角度的散射矩阵不同,经过不同颗粒物散射后的偏振状态不同,偏振紫外光发生散射的同时,就将待测粒子相关信息加载到了光束偏振态上,在雾霾检测中可以测得散射后的偏振状态,计算得到雾霾粒子的散射矩阵,用于分析雾霾粒子粒径。
因此,需要结合日盲紫外光偏振的特性,运用偏振态测量技术,来实现对大气雾霾粒子粒径的检测。
发明内容
本发明的目的是提供基于日盲紫外光进行雾霾粒子检测的装置,解决了现有技术条件下,雾霾检测技术的精度不高、受环境影响因素大、成本高,不能对大气雾霾颗粒物进行实时快速动态检测的问题。
本发明的另一目的是提供利用上述检测装置进行雾霾粒子检测的方法。
本发明所采用的技术方案是,基于日盲紫外光进行雾霾粒子检测的装置,包括用于发射紫外光的LED灯,沿紫外光传播方向依次设有空间滤波器、准直透镜、起偏器、λ/4波片、分子气体容器和光电探测器;起偏器、λ/4波片分别与步进电机控制器连接,计算机与步进电机控制器连接;
光电探测器包括沿光路依次设置的镀膜分光器、偏振分光棱镜a、偏振分光棱镜b、光电倍增管a、光电倍增管b、光电倍增管c和光电倍增管d;其中,镀膜分光器与入射光呈45°夹角放置,偏振分光棱镜a位于镀膜分光器后方与入射光共线放置,镀膜分光器、偏振分光棱镜a、光电倍增管a和光电倍增管b构成反射光路;偏振分光棱镜b位于镀膜分光器正上方与入射光垂直放置,镀膜分光器、偏振分光棱镜b、光电倍增管c和光电倍增管d构成透射光路。
本发明所采用的另一技术方案是,基于日盲紫外光进行雾霾粒子检测的方法,采用一种基于日盲紫外光进行雾霾粒子检测的装置,其结构为:包括用于发射紫外光的LED灯(1),沿紫外光传播方向依次设有空间滤波器(2)、准直透镜(3)、起偏器(4)、λ/4波片(5)、分子气体容器(6)和光电探测器(7);起偏器(4)、λ/4波片(5)分别与步进电机控制器(8)连接,计算机(9)与步进电机控制器(8)连接;
其中,光电探测器(7)包括沿光路依次设置的镀膜分光器(10)、偏振分光棱镜a(11)、偏振分光棱镜b(12)、光电倍增管a(13)、光电倍增管b(14)、光电倍增管c(15)和光电倍增管d(16);其中,镀膜分光器(7)与入射光呈45°夹角放置,偏振分光棱镜a(11)位于镀膜分光器(10)后方与入射光共线放置,镀膜分光器(10)、偏振分光棱镜a(11)、光电倍增管a(13)和光电倍增管b(14)构成反射光路;偏振分光棱镜b(12)位于镀膜分光器(10)正上方与入射光垂直放置,镀膜分光器(10)、偏振分光棱镜b(12)、光电倍增管c(15)和光电倍增管d(16)构成透射光路;
基于上述装置的检测方法按以下步骤进行:
步骤1,测量光电探测器的仪器矩阵A;
步骤2,利用步骤1测得的仪器矩阵A,得偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout;
步骤3,根据mie散射理论,得待测雾霾粒子的粒径。
本发明的特征还在于,
步骤1中,仪器矩阵A的测量过程为:
步骤1.1:打开LED灯,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14];
步骤1.2:步骤1经偏振作用后的紫外光穿过未放置雾霾粒子的分子气体容器,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ia1和Ib1;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ic1和Id1;令I1=[Ia1 Ib1 Ic1 Id1];
步骤1.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤1.2的方法测得Ia2、Ib2、Ic2和Id2,令I2=[Ia2 Ib2 Ic2 Id2];
步骤1.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤1.2的方法测得Ia3、Ib3、Ic3和Id3,令I3=[Ia3 Ib3 Ic3 Id3];
步骤1.5:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤1.2的方法测得Ia4、Ib4、Ic4和Id4,令I4=[Ia4 Ib4 Ic4 Id4];
步骤1.6:将步骤1.1~1.5所测得偏振紫外光的Sokes矢量及电流值带入Muller矩阵计算公式I=SA中,即得光电探测器的仪器矩阵A:
A=Sin -1Im,
其中,Sin=[Sin0;Sin1;Sin2;Sin3],Sin -1为Sin的逆矩阵;
Im=[I1;I2;I3;I4]。
步骤2中,偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout的获取过程为:
步骤2.1:开LED灯,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14],在偏振后的紫外光穿过放置雾霾粒子的分子气体容器,接收端的光电探测仪是放置于散射后某角度,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流I11和I12;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流I13和I14;令Iout1=[I11 I12 I13 I14];
步骤2.2:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤2.1的方法测得I21、I22、I23和I24,令Iout2=[I21 I22 I23 I24];
步骤2.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤2.1的方法测得I31、I32、I33和I34,令Iout3=[I31 I32 I33 I34];
步骤2.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤2.1的方法测得I41、I42、I43和I44,令Iout4=[I41 I42 I43 I44];
步骤2.5:将步骤1所得仪器矩阵A和步骤2.1~2.4所得电流值带入公式I=SA中,可得经雾霾粒子散射后偏振紫外光的Stokes矢量矩阵Sout:
Sout=IoutA-1,
Iout=[Iout1;Iout2;Iout3;Iout4];A-1光电探测器的仪器矩阵A的逆矩阵。
步骤3中,粒子的粒径的获取过程为:
步骤3.1:根据mie散射理论,雾霾球形粒子振幅函数表达式为:
其中,
式中,x1(θs),x2(θs)为振幅函数;
an,bn为Mie散射系数;
m为粒子折射率;
α为无因次粒径参量,α=πD/λ;
λ为入射光波长;
D为粒子的粒径;
θs为散射角;
ψn为n阶第一类Bessel函数;
ξn为n阶第二类Bessel函数;
Pn (1)为1阶n次的连带勒让德函数;
步骤3.2:将步骤2.5所得的Stokes矢量矩阵Sout带入公式Sout=SinRM(ψSA)M(θS)RM(ψPS)中,即得待测雾霾粒子的散射矩阵:
式中,M(θS)为待测雾霾粒子的散射矩阵;
RM(ψPS)为起偏参考面到散射参考面旋转矩阵;
RM(ψSA)为散射参考面到检偏参考面旋转矩阵;
Sin -1=[Sin0;Sin1;Sin2;Sin3]-1为输入的Stokes矢量矩阵的可逆矩阵;
Sout为雾霾粒子散射后偏振紫外光的Stokes矢量矩阵;
其中,M11=(|x1(θs)|2+|x2(θs)|2);
M12=(|x2(θs)|2-|x1(θs)|2);
M33=(x2(θs)*x1(θs)+x2(θs)x1(θs)*);
M34=(x1(θs)x2(θs)*-x2(θs)x1(θs)*);
x1(θs),x2(θs)为振幅函数;
步骤3.3:将步骤3.1和步骤3.2所得数值带入待下式中,即可求得粒子粒径,并:
式中,an,bn为Mie散射系数;
τn,πn为散射角θs的函数。
本发明的有益效果是,
1.克服暗盒测量空间的局限性,不受环境湿度和温度的影响,不受背景光干扰,可实现对雾霾的全天候开放式测量,同时降低设备成本,适用于工业生产。
2.提高了粒子形态等的精确度,可以分析雾霾的粒子的粒径,为各地区治理雾霾提供帮助。
附图说明
图1为本发明基于日盲紫外光进行雾霾粒子检测装置的结构示意图;
图2为本发明基于日盲紫外光进行雾霾粒子检测的装置中光电探测器的结构示意图;
图3为本发明基于日盲紫外光进行雾霾粒子检测方法中mie理论散射示意图;
图4为本发明基于日盲紫外光进行雾霾粒子检测方法中紫外光非视距散射原理示意图;
图5为本发明基于日盲紫外光进行雾霾粒子检测方法中所得M11值与雾霾粒子粒径关系变化示意图。
图1中,1.LED灯,2.空间滤波器,3.准直透镜,4.起偏器,5.λ/4波片,6.分子气体容器,7.光电探测器8.步进电机控制器,9.计算机;
图2中,10.镀膜分光器,11.偏振分光棱镜,12.偏振分光棱镜b,13.光电倍增管a,14.光电倍增管b,15.光电倍增管c,16.光电倍增管d。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
基于日盲紫外光进行雾霾粒子检测的装置,包括用于发射紫外光的LED灯1,沿紫外光传播方向依次设有空间滤波器2、准直透镜3、起偏器4、λ/4波片5、分子气体容器6和光电探测器7;起偏器4、λ/4波片5分别与步进电机控制器8连接,计算机9与步进电机控制器8连接;本发明基于日盲紫外光进行雾霾粒子检测的装置示意图如图1所示。
光电探测器7包括沿光路依次设置的镀膜分光器10、偏振分光棱镜a11、偏振分光棱镜b12、光电倍增管a13、光电倍增管b14、光电倍增管c15和光电倍增管d16;其中,镀膜分光器10与入射光呈45°夹角放置,偏振分光棱镜a11位于镀膜分光器10后方与入射光共线放置,镀膜分光器10、偏振分光棱镜a11、光电倍增管a13和光电倍增管b14构成反射光路;偏振分光棱镜b12位于镀膜分光器10正上方与入射光垂直放置,镀膜分光器10、偏振分光棱镜b12、光电倍增管c15和光电倍增管d16构成透射光路。本发明基于日盲紫外光进行雾霾粒子检测的装置中光电探测器的结构示意图如图2所示。
基于日盲紫外光进行雾霾粒子检测的方法,按以下步骤进行:
步骤1,测量光电探测器的仪器矩阵A:
步骤1.1:打开LED灯,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14];
步骤1.2:步骤1经偏振作用后的紫外光穿过未放置雾霾粒子的分子气体容器,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ia1和Ib1;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ic1和Id1;令I1=[Ia1 Ib1 Ic1 Id1];
步骤1.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤1.2的方法测得Ia2、Ib2、Ic2和Id2,令I2=[Ia2 Ib2 Ic2 Id2];
步骤1.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤1.2的方法测得Ia3、Ib3、Ic3和Id3,令I3=[Ia3 Ib3 Ic3 Id3];
步骤1.5:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤1.2的方法测得Ia4、Ib4、Ic4和Id4,令I4=[Ia4 Ib4 Ic4 Id4];
步骤1.6:将步骤1.1~1.5所测得偏振紫外光的Sokes矢量及电流值带入Muller矩阵计算公式I=SA中,即得光电探测器的仪器矩阵A:
A=Sin -1Im,
其中,Sin=[Sin0;Sin1;Sin2;Sin3],Sin -1为Sin的逆矩阵;
Im=[I1;I2;I3;I4]。
步骤2,利用步骤1测得的仪器矩阵A,得偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout:
步骤2.1:开LED灯,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14],在偏振后的紫外光穿过放置雾霾粒子的分子气体容器,接收端的光电探测仪是放置于散射后某角度,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流I11和I12;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流I13和I14;令Iout1=[I11 I12 I13 I14];
步骤2.2:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤2.1的方法测得I21、I22、I23和I24,令Iout2=[I21 I22 I23 I24];
步骤2.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤2.1的方法测得I31、I32、I33和I34,令Iout3=[I31 I32 I33 I34];
步骤2.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤2.1的方法测得I41、I42、I43和I44,令Iout4=[I41 I42 I43 I44];
步骤2.5:将步骤1所得仪器矩阵A和步骤2.1~2.4所得电流值带入公式I=SA中,可得经雾霾粒子散射后偏振紫外光的Stokes矢量矩阵Sout:
Sout=IoutA-1,
其中,Iout=[Iout1;Iout2;Iout3;Iout4];
A为光电探测器的仪器矩阵,A-1为其逆矩阵。
步骤3,根据mie散射理论,得待测雾霾粒子的粒径,本发明基于日盲紫外光进行雾霾粒子检测方法中mie理论散射示意图如图3所示,其中,θs入射光与散射光的夹角;本发明基于日盲紫外光进行雾霾粒子检测方法中紫外光非视距散射原理示意图如图4所示,图中,r为收发距离,βT收送端仰角,βR接收端仰角,θs为粒子的散射角度,θT发射端视场半角,θR接收端视场半角,αT为发射端旋转角,αR为接收端旋转角;
步骤3.1:根据mie散射理论,雾霾球形粒子振幅函数表达式为:
其中,
式中,x1(θs),x2(θs)为振幅函数;
an,bn为Mie散射系数;
m为粒子折射率;
α为无因次粒径参量,α=πD/λ;
λ为入射光波长;
D为粒子的粒径;
θs为散射角;ψn为n阶第一类Bessel函数;
ξn为n阶第二类Bessel函数;
Pn (1)为1阶n次的连带勒让德函数;
步骤3.2:将步骤2.5所得的Stokes矢量矩阵Sout带入公式Sout=SinRM(ψSA)M(θS)RM(ψPS)中,得待测雾霾粒子的散射矩阵:
式中,M(θS)为待测雾霾粒子的散射矩阵;
RM(ψPS)为起偏参考面到散射参考面旋转矩阵;
RM(ψSA)为散射参考面到检偏参考面旋转矩阵;
Sin -1=[Sin0;Sin1;Sin2;Sin3]-1为输入的Stokes矢量矩阵的可逆矩阵;
Sout为雾霾粒子散射后偏振紫外光的Stokes矢量矩阵;
其中,M11=(|x1(θs)|2+|x2(θs)|2);
M12=(|x2(θs)|2-|x1(θs)|2);
M33=(x2(θs)*x1(θs)+x2(θs)x1(θs)*);
M34=(x1(θs)x2(θs)*-x2(θs)x1(θs)*);
x1(θs),x2(θs)为振幅函数;
步骤3.3:将步骤3.1和步骤3.2所得数值带入待下式中,即可求得粒子粒径,并:
式中,an,bn为Mie散射系数;
τn,πn为散射角θs的函数。
利用本发明方法进行雾霾粒子检测的具体过程为:
步骤1.1:打开LED灯,入射波长为365nm,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12a13 a14]=[1 0 0.5 0];
步骤1.2:步骤1经偏振作用后的紫外光穿过未放置雾霾粒子的分子气体容器,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ia1和Ib1;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ic1和Id1;令I1=[Ia1 Ib1 Ic1 Id1]=[0.3352 0.0816 0.1466 0.0141];
步骤1.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24]=[1 1 0 0],重复步骤1.2的方法测得Ia2、Ib2、Ic2和Id2,令I2=[Ia2 Ib2 Ic2 Id2]=[0.3717 -0.2481 -0.0461 0.0051];
步骤1.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34]=[1 0 0 1],重复步骤1.2的方法测得Ia3、Ib3、Ic3和Id3,令I3=[Ia3 Ib3 Ic3 Id3]=[0.5783 -0.0865 0.2116 0.4194];
步骤1.5:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44]=[-1 0.5 0 0],重复步骤1.2的方法测得Ia4、Ib4、Ic4和Id4,令I4=[Ia4 Ib4 Ic4 Id4]=[-0.0426 0.0245 -0.1963 -0.0168];
步骤1.6:将步骤1.1~1.5所测得偏振紫外光的Sokes矢量及电流值带入Muller矩阵计算公式I=SA中,即得光电探测器的仪器矩阵A:
A=Sin -1Im
其中,Sin=[Sin0;Sin1;Sin2;Sin3]=[1 0 0.5 0
1 1 0 0
1 0 0 1
-1 0.5 0 0]
Im=[I1;I2;I3;I4]=[0.3352 0.0816 0.1466 0.0141
0.3717 -0.2481 -0.0461 0.0051
0.5783 -0.0865 0.2116 0.4194
-0.0426 0.0245 -0.1963 -0.0168];
即得光电探测器的仪器矩阵A=[0.1523 -0.099 0 0.1155 0.0129
0.2194 -0.1491 -0.1616 -0.0078
0.3658 0.3612 0.0622 0.0023
0.4260 0.0125 0.0961 0.4065]
步骤2.1:打开LED灯,紫外光依次通过空间滤波器、准直透镜、偏振器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14]=[1 0 0.50],穿过放置雾霾粒子的分子气体容器,接收端的光电探测仪是放置于散射后45°,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流I11和I12;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流I13和I14;令Iout1=[I11 I12 I13 I14]=1.0e-099*[0.0640 -0.2462 0.0020 0.0055];
步骤2.2:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24]=[1 1 0 0],重复步骤2.1的方法测得I21、I22、I23和I24,令Iout2=[I21 I22 I23 I24]=1.0e-099*[-0.1851 -0.29710.06890.0091];
步骤2.3:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34]=[1 0 0 1],重复步骤2.1的方法测得I31、I32、I33和I34,令Iout3=[I31 I32 I33 I34]=1.0e-099*[0.2999 -0.1338 0.09970.2438];
步骤2.4:计算机通过控制步进电机控制器调整偏振器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44]=[-1 0.5 0 0],重复步骤2.1的方法测得I41、I42、I43和I44,令Iout4=[I41 I42 I43 I44]=1.0e-099*[-0.1698 0.0631 -0.0311-0.0056];
步骤2.5:将步骤1所得仪器矩阵A和步骤2.1~2.4所得电流值带入公式I=SA中,可得经雾霾粒子散射后偏振紫外光的Stokes矢量矩阵Sout:
Sout=IoutA-1
Iout=[Iout1;Iout2;Iout3;Iout4]=1.0e-099*[0.0640 -0.2462 0.0020 0.0055
-0.1851 -0.2971 0.0689 0.0091
0.2999 -0.1338 0.0997 0.2438
-0.1698 0.0631 -0.0311 -0.0056]
带入光电探测器的仪器矩阵的逆矩阵A-1,得:
Sout=1.0e-099*[0.7077 0.3632 -0.3376 -0.0000
0.7215 -0.1796 -0.6988 0.0000
0.6185 0.1031 -0.1785 0.5831
-0.5669 -0.2444 -0.0816 0.0000]
步骤3.1:将步骤2中入射的Stokes矢量矩阵Sin以及步骤2.5所得的Stokes矢量矩阵Sout带入公式Sout=SinRM(ψSA)M(θS)RM(ψPS)中,即得待测雾霾粒子的散射矩阵:
即得M(θS)=1.0e-099*[0.6185 -0.2062 0.0000 0
-0.2062 0.6185 -0.0000 -0.0000
0.0000 0.0000 0.5831 0.0000
0.0000 -0.0000 -0.0000 0.5831]
步骤3.2:将步骤3.1所得数值带入待下式中:
其中,
λ=365nm;m=1.64;θs=45°;
通过仿真所得M11与粒子粒径关系如图5所示,经即可得粒子粒径为5μm。
Claims (5)
1.基于日盲紫外光进行雾霾粒子检测的装置,其特征在于,包括用于发射紫外光的LED灯(1),沿紫外光传播方向依次设有空间滤波器(2)、准直透镜(3)、起偏器(4)、λ/4波片(5)、分子气体容器(6)和光电探测器(7);起偏器(4)、λ/4波片(5)分别与步进电机控制器(8)连接,计算机(9)与步进电机控制器(8)连接;
所述光电探测器(7)包括沿光路依次设置的镀膜分光器(10)、偏振分光棱镜a(11)、偏振分光棱镜b(12)、光电倍增管a(13)、光电倍增管b(14)、光电倍增管c(15)和光电倍增管d(16);其中,镀膜分光器(10)与入射光呈45°夹角放置,偏振分光棱镜a(11)位于镀膜分光器(10)的正上方与入射光方向垂直放置,光电倍增管a(13)和光电倍增管b(14)位于偏振分光棱镜a(11)两侧,且关于入射光垂直方向对称,镀膜分光器(10)、偏振分光棱镜a(11)、光电倍增管a(13)和光电倍增管b(14)构成反射光路;偏振分光棱镜b(12)位于镀膜分光器(10)的后方与入射光方向共线放置,光电倍增管a(15)和光电倍增管b(16)位于偏振分光棱镜a(12)两侧,且关于入射光方向对称,镀膜分光器(10)、偏振分光棱镜b(12)、光电倍增管c(15)和光电倍增管d(16)构成透射光路。
2.基于日盲紫外光进行雾霾粒子检测的方法,其特征在于,采用一种基于日盲紫外光进行雾霾粒子检测的装置,其结构为:包括用于发射紫外光的LED灯(1),沿紫外光传播方向依次设有空间滤波器(2)、准直透镜(3)、起偏器(4)、λ/4波片(5)、分子气体容器(6)和光电探测器(7);起偏器(4)、λ/4波片(5)分别与步进电机控制器(8)连接,计算机(9)与步进电机控制器(8)连接;
所述光电探测器(7)包括沿光路依次设置的镀膜分光器(10)、偏振分光棱镜a(11)、偏振分光棱镜b(12)、光电倍增管a(13)、光电倍增管b(14)、光电倍增管c(15)和光电倍增管d(16);其中,镀膜分光器(10)与入射光呈45°夹角放置,偏振分光棱镜a(11)位于镀膜分光器(10)后方与入射光共线放置,镀膜分光器(10)、偏振分光棱镜a(11)、光电倍增管a(13)和光电倍增管b(14)构成反射光路;偏振分光棱镜b(12)位于镀膜分光器(10)正上方与入射光垂直放置,镀膜分光器(10)、偏振分光棱镜b(12)、光电倍增管c(15)和光电倍增管d(16)构成透射光路;
基于上述装置的检测方法按以下步骤进行:
步骤1,测量光电探测器的仪器矩阵A;
步骤2,利用步骤1测得的仪器矩阵A,得偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout;
步骤3,根据mie散射理论,得待测雾霾粒子的粒径。
3.根据权利要求2所述的基于日盲紫外光进行雾霾粒子检测的方法,其特征在于,步骤1中,仪器矩阵A的测量过程为:
步骤1.1:打开LED灯,紫外光依次通过空间滤波器、准直透镜、起偏器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14];
步骤1.2:步骤1.1经偏振作用后的紫外光穿过未放置雾霾粒子的分子气体容器,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ia1和Ib1;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流Ic1和Id1;令I1=[Ia1 Ib1 Ic1 Id1];
步骤1.3:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤1.2的方法测得Ia2、Ib2、Ic2和Id2,令I2=[Ia2 Ib2 Ic2 Id2];
步骤1.4:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤1.2的方法测得Ia3、Ib3、Ic3和Id3,令I3=[Ia3 Ib3 Ic3 Id3];
步骤1.5:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤1.2的方法测得Ia4、Ib4、Ic4和Id4,令I4=[Ia4 Ib4 Ic4 Id4];
步骤1.6:将步骤1.1~1.5所测得偏振紫外光的Sokes矢量及电流值带入Muller矩阵计算公式I=SA中,即得光电探测器的仪器矩阵A:
A=Sin -1Im,
其中,Sin=[Sin0;Sin1;Sin2;Sin3],Sin -1为Sin的逆矩阵;
Im=[I1;I2;I3;I4]。
4.根据权利要求2所述的基于日盲紫外光进行雾霾粒子检测的方法,其特征在于,步骤2中,偏振紫外光照射待测雾霾粒子后散射光的Stokes矢量矩阵Sout的获取过程为:
步骤2.1:开LED灯,紫外光依次通过空间滤波器、准直透镜、起偏器和λ/4波片,计算机获取并记录偏振作用后紫外光的Stokes矢量:Sin0=[a11 a12 a13 a14],在偏振后的紫外光穿过放置雾霾粒子的分子气体容器,接收端的光电探测器是放置于散射后某角度,再通过镀膜分光器作用分为反射光和透射光,反射光通过棱镜a的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管a和光电倍增管b,光电倍增管a和光电倍增管b将偏振光信号转换为电信号输出,用电流表分别测得输出电流I11和I12;透射光通过棱镜b的分割作用形成两束正交的偏振光,偏振光分别进入光电倍增管c和光电倍增管d,光电倍增管c和光电倍增管d将偏振光信号转换为电信号输出,用电流表分别测得输出电流I13和I14;令Iout1=[I11 I12 I13 I14];
步骤2.2:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin1=[a21 a22 a23 a24],重复步骤2.1的方法测得I21、I22、I23和I24,令Iout2=[I21 I22 I23 I24];
步骤2.3:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin2=[a31 a32 a33 a34],重复步骤2.1的方法测得I31、I32、I33和I34,令Iout3=[I31 I32 I33 I34];
步骤2.4:计算机通过控制步进电机控制器调整起偏器和λ/4波片的位置,使偏振紫外光的Sokes矢量为Sin3=[a41 a42 a43 a44],重复步骤2.1的方法测得I41、I42、I43和I44,令Iout4=[I41 I42 I43 I44];
步骤2.5:将步骤1所得仪器矩阵A和步骤2.1~2.4所得电流值带入公式I=SA中,可得经雾霾粒子散射后偏振紫外光的Stokes矢量矩阵Sout:
Sout=IoutA-1,
Iout=[Iout1;Iout2;Iout3;Iout4];A-1光电探测器的仪器矩阵A的逆矩阵。
5.根据权利要求4所述的基于日盲紫外光进行雾霾粒子检测的方法,其特征在于,步骤3中,粒子的粒径的获取过程为:
步骤3.1:根据mie散射理论,雾霾球形粒子振幅函数表达式为:
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mfrac>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mfrac>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中,
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>m&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&xi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>m&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&xi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>m&psi;</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>dP</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>cos&theta;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi> </mi>
<msub>
<mi>cos&theta;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>dP</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>cos&theta;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>d&theta;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
式中,x1(θs),x2(θs)为振幅函数;
an,bn为mie散射系数;
m为粒子折射率;
α为无因次粒径参量,α=πD/λ;
λ为入射光波长;
D为粒子的粒径;
θs为散射角;ψn为n阶第一类Bessel函数;
ξn为n阶第二类Bessel函数;
Pn (1)为1阶n次的连带勒让德函数;
步骤3.2:将步骤2.5所得的Stokes矢量矩阵Sout带入公式Sout=SinRM(ψSA)M(θS)RM(ψPS)中,即得待测雾霾粒子的散射矩阵:
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>R</mi>
<mi>M</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>S</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mi>M</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mn>11</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mn>12</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mn>12</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mn>11</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mn>33</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mn>34</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>M</mi>
<mn>34</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mn>33</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
式中,M(θS)为待测雾霾粒子的散射矩阵;
RM(ψPS)为起偏参考面到散射参考面旋转矩阵;
RM(ψSA)为散射参考面到检偏参考面旋转矩阵;
Sin -1=[Sin0;Sin1;Sin2;Sin3]-1为输入的Stokes矢量矩阵的可逆矩阵;
Sout为雾霾粒子散射后偏振紫外光的Stokes矢量矩阵;
其中,M11=(|x1(θs)|2+|x2(θs)|2);
M12=(|x2(θs)|2-|x1(θs)|2);
M33=(x2(θs)*x1(θs)+x2(θs)x1(θs)*);
M34=(x1(θs)x2(θs)*-x2(θs)x1(θs)*);
x1(θs),x2(θs)为振幅函数;
步骤3.3:将步骤3.1和步骤3.2所得数值带入下式中,即可求得粒子粒径:
<mrow>
<mi>M</mi>
<mn>11</mn>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mfrac>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>)</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<mfrac>
<mrow>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&pi;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>)</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
式中,an,bn为mie散射系数;
τn,πn为散射角θs的函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510374742.5A CN105043944B (zh) | 2015-06-30 | 2015-06-30 | 基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510374742.5A CN105043944B (zh) | 2015-06-30 | 2015-06-30 | 基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105043944A CN105043944A (zh) | 2015-11-11 |
CN105043944B true CN105043944B (zh) | 2017-09-22 |
Family
ID=54450662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510374742.5A Expired - Fee Related CN105043944B (zh) | 2015-06-30 | 2015-06-30 | 基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105043944B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106018197B (zh) * | 2016-07-01 | 2019-06-14 | 聚光科技(杭州)股份有限公司 | 基于单光源的颗粒物在线监测装置及方法 |
CN109883901A (zh) * | 2019-02-15 | 2019-06-14 | 西安理工大学 | 一种无人机机载的紫外光探测雾霾粒子系统及其探测方法 |
CN109883902A (zh) * | 2019-02-28 | 2019-06-14 | 西安理工大学 | 基于日盲紫外光圆偏振的雾霾粒子检测装置及其检测方法 |
CN110687022B (zh) * | 2019-11-05 | 2021-09-24 | 珠海真理光学仪器有限公司 | 一种利用散射光的偏振差异测量颗粒折射率的方法及系统 |
CN111929203A (zh) * | 2020-06-10 | 2020-11-13 | 西安理工大学 | 一种基于日盲紫外光线偏振度的非球形灰霾粒子区分方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101363918B (zh) * | 2008-08-12 | 2011-01-26 | 阎锋 | 应用日盲紫外光对物体进行定位识别的方法 |
US8835831B2 (en) * | 2010-06-22 | 2014-09-16 | Zena Technologies, Inc. | Polarized light detecting device and fabrication methods of the same |
CN103323758B (zh) * | 2013-06-21 | 2016-01-27 | 江苏紫峰光电科技有限公司 | 日盲紫外成像式测距装置 |
CN103868593B (zh) * | 2014-03-07 | 2015-08-19 | 太原理工大学 | 基于光子晶体滤波和量子点光谱转换的日盲紫外成像装置 |
CN104637234B (zh) * | 2015-02-11 | 2017-01-11 | 中煤科工集团重庆研究院有限公司 | 基于激光散射测量原理的烟雾探测器检定装置及标定方法 |
CN104697906A (zh) * | 2015-03-18 | 2015-06-10 | 东南大学 | 一种基于近场散射的颗粒粒度测量装置及方法 |
-
2015
- 2015-06-30 CN CN201510374742.5A patent/CN105043944B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105043944A (zh) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105043944B (zh) | 基于日盲紫外光进行雾霾粒子检测的装置及其检测方法 | |
CN102692393B (zh) | 一种基于石墨烯偏振效应的折射率实时测定方法和装置 | |
CN110687022B (zh) | 一种利用散射光的偏振差异测量颗粒折射率的方法及系统 | |
US4306809A (en) | Polarimeter | |
CN101303256B (zh) | 基于液晶的嵌入式偏振态测量仪 | |
CN106597414A (zh) | 一种定标偏振激光雷达增益比的方法 | |
CN109883902A (zh) | 基于日盲紫外光圆偏振的雾霾粒子检测装置及其检测方法 | |
CN101738369B (zh) | 相位差检测装置 | |
CN108956391B (zh) | 探测大气中雾滴和气溶胶粒径谱分布的探测仪及探测方法 | |
CN107894395A (zh) | 一种基于腔衰荡光谱技术的大气气溶胶消光系数测量方法 | |
CN106769731B (zh) | 颗粒物浓度的测量方法及装置 | |
US9097647B2 (en) | Method for using polarization gating to measure a scattering sample | |
CN206208753U (zh) | 一种基于石墨烯表面波的高灵敏度超快折射率探测装置 | |
CN205941337U (zh) | 基于硅光电池的无法拉第调制器的自动旋光仪光学系统 | |
CN210775177U (zh) | 一种运用于libs中调节激光光束能量的装置 | |
CN103278432B (zh) | 一种滤膜完整性的检测装置及其检测方法 | |
CN102692367A (zh) | 纳米粒子辨识系统装置及其识别方法 | |
CN116008190A (zh) | 一种基于红外偏振谐波的高精度溶液浓度检测装置 | |
CN113776997B (zh) | 环境干扰抑制的大气颗粒物穆勒矩阵测量方法 | |
CN107192679A (zh) | 一种基于导光毛细管的光度分析仪及其检测方法 | |
CN110426346B (zh) | 一种运用于libs中调节激光光束能量的装置 | |
Wang et al. | Research of polarization imaging detection method for water surface target in foggy weather | |
CN102243385A (zh) | 测量硅基液晶成像芯片光电特性曲线的装置 | |
CN204214777U (zh) | 一种多功能食品安全综合检测仪 | |
Raman et al. | Polarization signatures of airborne particulates |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170922 Termination date: 20210630 |