本发明属于纳米颗粒粒径检测,特别涉及一种面向动态光散射测量的双分散系粒径精确反演算法。
背景技术:
1、动态光散射技术是一种快速有效的测量纳米及亚微米级粒径颗粒的有效方法,被广泛应用在医学、化学、生物和高分子材料等领域。当激光器发出激光入射到样品池中正在做布朗运动的纳米颗粒时,颗粒的散射光信号围绕某一平均值随时间不断的波动。布朗运动的快慢与颗粒粒径有关,因此散射光信号的波动包含颗粒粒度信息。通过测量液体样品中颗粒散射光强随时间的变化曲线得到光强的自相关函数,并通过反演算法拟合得到颗粒的粒径信息。
2、从散射光相关函数获得颗粒粒度分布需要求解第一类fredholm积分方程,这是一个典型的不适定问题。由于该方程的病态解问题,颗粒粒度的反演一直是动态光散射法纳米颗粒测量技术中的难点。为了得到准确的粒径分布,先后提出了多种基于不同原理的粒径反演方法,如累积分析法、加权贝叶斯反演算法、非负约束最小二乘法等。2008年,刘桂强等在《光电工程》中公开了一种基于累积分析法的颗粒粒径反演方法,采取了自动斜率法确定基线值并对颗粒粒径进行反演。2020年,梁一卓等在《光子学报》第49卷第10期提出了加权贝叶斯反演算法,实现了单峰体系的宽分布及窄分布的粒度反演。2021年,申晋等在发明专利cn113552031a中利用非负最小二乘法求解得到颗粒单、双峰粒径分布。然而,加权贝叶斯反演算法存在迭代次数多、收敛速度较慢的问题;累积法只能准确检测单分散系纳米颗粒粒径,当溶液中存在两种及以上颗粒时,反演结果常常会产生较大误差;非负最小二乘法对于纳米颗粒双分散系样品反演结果不稳定,且迭代过程受初值影响较大;因此,在保证精度的前提下,如何通过算法得到双分散系粒径的快速、稳定、高精度反演结果,成为本领域亟待解决的问题。
技术实现思路
1、本发明要解决的技术问题是克服累积法只能准确检测单分散系纳米颗粒粒度以及非负最小二乘法对于双分散系纳米颗粒粒径精度较低的不足,提供一种面向动态光散射测量的双分散系粒径精确计算方法。该方法采用levenberg-marquardt算法,在迭代过程中引入阻尼因子参数,增强了全局收敛性;提出了基于一阶多项式拟合的迭代初值的选取方法,减少了迭代次数,降低了初值选取对拟合结果的影响,实现双分散系颗粒粒径稳定、高精度测量。
2、为了达到上述目的,本发明采用的技术方案为:
3、一种面向动态光散射测量的双分散系粒径精确计算方法。首先,通过对光电探测器采集到的散射光强数据进行自相关计算得到散射光强自相关函数;然后,对归一化光强自相关函数进行拟合处理,从而更加准确地获得衰减线宽系数;最后,通过stokes-einstein公式得到双分散系颗粒粒径测量结果。具体步骤如下:
4、第一步:散射光数据采集与处理
5、采集散射光强信号,设置光子采集频率f,相关时间τ、采样时间t,相关器通道数m,确定散射光强脉冲序列w=[w1,w2,w3,...wn]t,wn表示光子脉冲数量,n表示光子计数器在采样时间内的采样点数量。通过对脉冲序列w进行自相关运算得到数组τm表示通道m的相关时间,m=(1,2,…m),并对其进行归一化处理得到自相关数组g(τm)。
6、第二步:基于levenberg-marquardt算法的光强自相关函数反演初值
7、在动态光散射测量过程中,对满足高斯分布的光场,多分散颗粒系的归一化光强自相关函数为单指数加权之和。理论上双分散颗粒系归一化光强自相关函数g(τ)可表示为:
8、
9、其中,c1=g(γ1)、c2=g(γ2),表示两种不同粒径颗粒的衰减线宽函数;τ表示相关时间;γi=diq2表示颗粒的散射光衰减线宽,i=1,2表示不同粒径颗粒种类。其中:di表示颗粒在溶液中的扩散系数,表示散射矢量的大小,λ表示激光在真空中的波长,n0表示分散介质的折射率,θ表示散射光与入射光之间夹角。
10、令m表示通道数,两边取对数并进行移项得ln[g(τm)]=lnc0-γ0τm,获得迭代初值c1(0)=c2(0)=c0,γ1(0)=γ2(0)=γ0。
11、计算衰减线宽γ1、γ2的最优解:
12、首先,设未知参数矩阵b=[b1,b2,b3,b4]=[c1,c2,γ1,γ2],此时,归一化光强自相关函数可表示为g(τm,b),将获得的迭代初值b(0)=[c1(0),c2(0),γ1(0),γ2(0)]代入公式(1)可得到初始解g(τm,b(0))。将g(τm,b)在b(0)处按泰勒级数展开,略去二次及二次以上的项得到下式,
13、
14、其次,对式(2)采用最小二乘法处理,
15、
16、式中,m表示通道数;s表示归一化光强自相关函数g(τm,b)与数组g(τm)之间的2范数;μ(μ≥0)表示阻尼因子;j表示矩阵b中的未知参数数目。当函数s取得最小值时,可以得到函数g(τm,b)的最优解。因此,分别对未知参数b1、b2、b3、b4求解一阶偏导数,并使一阶偏导数为0,得到如下方程组的表达形式,
17、
18、其中,j、k均表示矩阵b中的未知参数数目;
19、最终,最优解与参数矩阵b及阻尼因子μ相关。设置通道数m,最大迭代次数xmax,参数初值b(0)、阻尼因子初值μ0、初始迭代步长ζ,迭代精度δ。在迭代过程中,若相邻参数的范数cx=||b(x)-b(x-1)||衰减,则减小阻尼因子μ值,反之增大μ值,直至满足cx≤δ,结束迭代。此时可以得到s取得最小值时的参数矩阵b,从而得到衰减线宽的最优解γ1=b3、γ2=b4。
20、第三步:计算双分散系颗粒粒径。
21、根据上述通过levenberg-marquardt算法对于光强自相关函数的反演结果得到衰减线宽的最优拟合解γ1、γ2,通过stokes-einstein公式即可得到双分散系颗粒粒径di,
22、
23、式中,kb表示玻尔兹曼常数,t表示开氏温度,η表示溶液粘度系数,di表示待测样品颗粒的直径,γi=diq2,其中:di表示颗粒在溶液中的扩散系数,表示散射矢量,λ表示激光在真空中的波长,n0表示分散介质的折射率,θ表示散射光与入射光之间夹角。
24、本发明所具有的有益效果是:
25、本发明的双分散系粒度分levenberg-marquardt动态光散射反演方法,克服累积法只能准确检测单分散系纳米颗粒粒度以及非负最小二乘法对于双分散系纳米颗粒粒径精度较低的问题。该方法引入了阻尼因子参数,提高了反演精度,同时在迭代初值的选取过程中更加接近衰减线宽系数的真实值,减少迭代次数,提高计算速度,从而得到稳定、高精度的反演结果。
1.一种面向动态光散射测量的双分散系粒径精确反演算法,其特征在于,首先,通过对光电探测器采集到的散射光强数据进行自相关计算得到散射光强自相关函数;然后,对归一化光强自相关函数进行拟合处理,从而更加准确地获得衰减线宽系数;最后,通过stokes-einstein公式得到双分散系颗粒粒径测量结果。
2.根据权利要求1所述的一种面向动态光散射测量的双分散系粒径精确反演算法,其特征在于,具体步骤如下:
3.根据权利要求2所述的一种面向动态光散射测量的双分散系粒径精确反演算法,其特征在于,所述的第二步中,计算衰减线宽γ1、γ2的最优解的具体步骤为: