一种非球形气溶胶浓度谱探测分析方法、系统及设备与流程

1.本发明涉及气溶胶浓度分析技术领域，特别涉及一种非球形气溶胶浓度谱探测分析方法、系统及设备。


背景技术：

2.气溶胶是大气环境的重要组成部分，沙尘气溶胶等颗粒物的分布情况反映了大气的环境质量，大气环境质量直接影响到人类身体健康。人们目前对沙尘气溶胶的粒子分布进行监测和计算时存在很大误差。首先在计算的时候，通常将气溶胶粒子形状假设为球形。实际大气中的气溶胶粒子往往形态各异，例如沙尘气溶胶粒子是典型的非球形形态。因此计算非球形气溶胶粒子的粒子谱分布时，仍采用球形理论会造成计算不准确，影响气溶胶粒子谱分布计算准确性。如何解决非球形气溶胶粒子谱分布计算问题是当前的关键问题。


技术实现要素：

3.本发明的目的是提供一种非球形气溶胶浓度谱探测分析方法、系统及设备，以解决非球形气溶胶粒子谱分布计算不准确的问题。
4.为实现上述目的，本发明提供了如下方案：
5.一种非球形气溶胶浓度谱探测分析方法，包括：
6.利用dda算法计算沙尘气溶胶的光学特性参数；所述光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子；
7.基于所述光学特性参数计算核函数；
8.对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；所述光学数据包括消光系数和后散射系数；
9.基于所述核函数和所述光学数据构建沙尘气溶胶粒子谱分布反演函数；
10.根据所述沙尘气溶胶粒子谱分布反演函数确定确定大气中的沙尘气溶胶粒子谱分布；
11.利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。
12.可选地，利用dda算法计算沙尘气溶胶的光学特性参数，具体包括：
13.基于dda算法，将光学特性截面与等几何体截面之比定义为光学特性参数；
[0014][0015]
其中，q
ext
表示消光效率因子，q
abs
表示吸收效率因子，q
sca
表示散射效率因子，c
ext
消光截面，c
abs
表示吸收截面，c
sca
表示散射截面，r表示气溶胶粒子半径。
[0016]
可选地，所述核函数的计算公式如下：
[0017]
[0018]
其中，k表示核函数，q表示光学特性参数。
[0019]
可选地，基于所述核函数和所述光学数据构建气溶胶粒子谱分布反演函数，具体包括：
[0020]
利用混合优化算法求解最优拉格朗日乘子；所述混合优化算法包括粒子群算法和遗传算法；
[0021]
基于所述最优拉格朗日乘子和所述核函数确定权重因子；
[0022]
基于所述光学数据和所述权重因子构建沙尘气溶胶粒子谱分布反演函数。
[0023]
可选地，所述沙尘气溶胶粒子谱分布反演函数的表达式如下：
[0024][0025]
其中，g表示沙尘气溶胶光学数据，w表示权重因子，f(r)表示沙尘气溶胶粒子谱分布函数，r表示气溶胶粒子半径。
[0026]
可选地，还包括：
[0027]
计算偏振激光雷达探测大气中的沙尘气溶胶的回波信号的退偏振比；
[0028]
根据所述退偏振比确定沙尘气溶胶为非球形形态。
[0029]
本发明还提供了一种非球形气溶胶浓度谱探测分析系统，包括：
[0030]
光学特性参数计算模块，用于利用dda算法计算沙尘气溶胶的光学特性参数；所述光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子；
[0031]
核函数计算模块，用于基于所述光学特性参数计算核函数；
[0032]
沙尘气溶胶光学数据计算模块，用于对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；所述光学数据包括消光系数和后散射系数；
[0033]
沙尘气溶胶粒子谱分布反演函数构建模块，用于基于所述核函数和所述光学数据构建沙尘气溶胶粒子谱分布反演函数；
[0034]
沙尘气溶胶粒子谱分布确定模块，用于根据所述沙尘气溶胶粒子谱分布反演函数确定确定大气中的沙尘气溶胶粒子谱分布；
[0035]
近地面沙尘气溶胶粒子谱分布确定模块，用于利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。
[0036]
本发明还提供了一种电子设备，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备上述的非球形气溶胶浓度谱探测分析方法。
[0037]
本发明还提供了一种计算机可读存储介质，其存储有计算机程序，所述计算机程序被处理器执行时实现上述的非球形气溶胶浓度谱探测分析方法。
[0038]
根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0039]
(1)由于沙尘气溶胶是非球形粒子，利用dda算法解决沙尘气溶胶当做球形气溶胶计算的问题；
[0040]
(2)利用优化算法解决拉格朗日乘子计算效率低的问题；
[0041]
(3)利用粒径谱仪对近地面沙尘气溶胶进行探测，补充近地面数浓度谱数据缺失的问题，实现沙尘气溶胶在垂直方向上数量关系的监测，达到保护大气环境、监控环境污染的目的。
附图说明
[0042]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0043]
图1为本发明实施例一提供的非球形气溶胶浓度谱探测分析方法的流程图；
[0044]
图2为本发明实施例一提供的非球形气溶胶浓度谱探测分析方法的详细流程图；
[0045]
图3为采用混合优化算法求解拉格朗日乘子的流程图；
[0046]
图4为确定初始种群数量的示意图；
[0047]
图5为确定交叉概率的示意图；
[0048]
图6为确定变异概率的示意图；
[0049]
图7为退偏振比验证沙尘粒子非球形特性的示意图；
[0050]
图8为不同时间大气和地面沙尘气溶胶粒子谱分布示意图。
具体实施方式
[0051]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0052]
本发明的目的是提供一种非球形气溶胶浓度谱探测分析方法、系统及设备，以解决非球形气溶胶粒子谱分布计算不准确的问题。
[0053]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0054]
实施例一
[0055]
本发明实施例一提供了一种非球形气溶胶浓度谱探测分析方法，如图1-2所示，该方法包括以下步骤：
[0056]
s1：利用dda算法计算沙尘气溶胶的光学特性参数；光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子。
[0057]
离散偶极子近似法(dda)：dda方法是计算任意几何形状、方向和组成的粒子对电磁波的散射和吸收的一种通用方法，即一种体积积分法。dda方法是用一组特定间距的点偶极子按顺序排列取代了散射体。这些偶极子相互作用，并与入射场相互作用，形成一个线性方程组，求解此方程组后得到偶极子极化。dda逐渐成为目前求解粒子光散射特性的热门方法之一，该方法建模简单、适用于大粒径、多数量等粒子特性，在模拟非球形粒子光散射方面优于其他方法。
[0058]
用数量有限的偶极子近似模拟一个特定目标的几何形状，可以将目标体离散化为有限个相互关联且按照特定顺序排列的偶极子阵列，在此目标体的边界处用间隔紧密但散射较弱的偶极子来更好的模拟边界区域。对于一个给定的目标形状，可以使用以下算法来生成偶极子阵列。假设目标的方向对于坐标系是固定的：
[0059]
(1)生成一个由晶格间距d和坐标轴xyz定义的离原点最近的初始晶格体。
[0060]
(2)将这种晶格体内的分布定义为连续体目标的体积。
[0061]
(3)设置不同的晶格间距d和坐标xyz的值，避免偶极子点数太大导致计算机超负荷无法计算，最终得到最优的拟合参数。
[0062]
(4)在目前的假设中，并未区分出晶格体边界的偶极子和晶格体内部的偶极子。n个晶格点组成了晶格体，通过d＝(v/n)
1/3
来调整晶格点间距，所以可以将原始晶格体v表示为nd3。
[0063]
对每一个小的晶格点j分配一个偶极子αj。
[0064]
由于粒子散射时空间朝向是随机分布的，因此本文中光学参数用粒子在各个方向上的光学数值的平均值来表示，q为光学特性参数，其平均值记为《q》，表达式如下：
[0065][0066]
如果晶格体是各向同性的，则极化率张量是对角线的，各分量相等α
j,xx
＝α
j,yy
＝α
j,zz
，也就是说α可以看成一个标量。如果晶粒材料是各向异性的，极化率张量与单位矩阵不成正比。对于任何单个偶极子，可以选择一个笛卡尔坐标系使a对角化，对角元素为α
j,xx
，α
j,yy
，α
j,zz
。如果设定的晶格体是均匀的，那么介电函数e与晶粒内的位置无关，可以选择一个e对角化的坐标系，其中有对角元素e
||
，每个单独的偶极子极化率张量αj，在这个坐标系中也将是对角的，对角元素α
j,||
。在本发明中会受到这样一种限制：所选择的坐标系中，介电张量即使是非齐次的，也是对角的，则所有的偶极子极化张量都可以同时对角化。极化率的表达式如下：
[0067][0068][0069][0070]
b2＝0.1648469b3＝-1.7700004
[0071][0072]
其中ej为位置j处目标体的介电函数，和分别为入射方向的单位矢量和偏振方向的单位矢量。波在晶格上传播时，α
ldr
是最好的计算公式，所以可以假设此公式也可以用来计算有限偶极子的分布。利用dda方法尝试不同的极化率后，证明ldr方法最适合|m|kd＜1。
[0073]
通过要求每个点偶极子具有与直径为d的有限球体相同的极化率，单个偶极子将位于一个简单的立方晶格上。令坐标轴按照j＝1,...,n顺序排列已占据的格点，一个晶格点代表一队偶极子。每个偶极子通过和电场之间的相互作用产生偶极矩，所有偶极子电场的总和构成了目标体的散射特性。其中pj为偶极子j的瞬时电偶极矩，其中ej代表位于rj处
的电场。电场的值同时与入射波和n-1个偶极子的散射值相关。
[0074][0075]
在坐标轴j从1,...,n的分布中，αj代表具有位于rj的电偶极子的极化率，则电偶极矩pj的表达式为：
[0076][0077]
由于在偶极子所在j处的延迟效应和散射效应影响，j处电场为入射波处电场与其他影响之差。a
ij
中每个元素组成了3
×
3的系数矩阵：
[0078][0079]ajj
＝α
j-1
[0080]
其中k＝ω/c，r
jk
＝|r
j-rk|，i是3
×
3的单位矩阵。此时散射问题可以表示为n组非齐次线性方程组：
[0081][0082]
当电偶极矩pj已知时，则利用光学定理计算粒子的前向散射的消光截面：
[0083][0084][0085]
散射截面原则上可以由消光和吸收截面的差得到:csca＝cext-cabs。当吸收占主导地位时，这就要求cext和cabs的计算精度很高。由于计算方法是迭代的，在获得解向量p精确到必要的程度时，可能在数值上不可行的。当然，可以通过计算振荡偶极子阵列的辐射功率来直接计算散射截面：
[0086][0087]
消光效率因子、吸收效率因子和散射效率因子可定义为光学特性截面与等几何体截面之比。此公式为求解非球形粒子的光学特性提供了一种高效稳定的方法：
[0088][0089]
式中f为散射方向的单位矢量，dω为立体角元。
[0090]
s2：基于光学特性参数计算核函数。
[0091][0092]
其中，k表示核函数，q表示光学特性参数。
[0093]
s3：对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；光学数据包括消光系数和后散射系数。
[0094]
利用多波长偏振激光雷达系统对大气中的沙尘气溶胶进行探测，由于沙尘气溶胶具有非球形特性，因此利用传统的mie方法进行光学特行计算是不准确的，需要对雷达回波信号进行反演，根据下式计算出沙尘气溶胶的消光系数和后向散射系数，为反演沙尘气溶胶粒子谱分布做数据支撑。
[0095]
b＝cak[0096][0097]
式中rm为远端边界值，c为常数，a为消光系数，b为后向散射系数，k一般取0.76-1.0之间，r表示沙尘气溶胶粒子半径，
s(r)
表示取距离平方矫正信号的自然对数。消光系数与后向散射系数构成沙尘气溶胶光学数据g。
[0098]
s4：基于核函数和光学数据构建沙尘气溶胶粒子谱分布反演函数。
[0099]
通过第一类数学fredholm积分方程搭建气溶胶粒子谱分布反演函数：
[0100][0101]
求解大条件数的线性方程的不适定反问题通常用tikhonov正则化方法，大条件数代表着任何微小的误差都会给最终结果带来巨大影响。正则化方法通过数学约束条件逼近原问题的解，此方法为解决不适定问题提供一个有效方法。
[0102]
构造反演函数时，将消光系数和后向散射系数按照下式表示：
[0103][0104][0105]
其中ai和bi代表消光系数和后向散射系数，通过步骤s3可以计算得到，i代表在355nm，532nm，1064nm波长。ka和kb代表消光核函数和后向散射核函数，可以通过dda方法求解。m代表环境复折射率，r代表气溶胶粒子的粒径，λi代表入射波长，f(r)代表沙尘气溶胶粒子谱分布。
[0106][0107][0108]
其中g为探测得到的消光系数或后向散射系数，aj为权重矩阵，q
p
为效率因子，bj为b样条函数。因此对气溶胶粒子谱分布求解的重点在得出权重wj。
[0109]
g＝aw+error
p
[0110]
其中g为n行1列的矩阵，w为n行n列的矩阵。为解决解的过拟合情况，采用正则化方法时引入正则化参数解决矩阵方程的病态。表达式如下所示：
[0111]
wj＝(a
jtaj
+γh)-1ajtg[0112]
其中γ为拉格朗日乘子，h为平滑矩阵。上式表明拉格朗日乘子γ的求解对wn的数值结果至关重要，权重因子误差最小需要选取合适的γ。若拉格朗日乘子取值越大，则正则化操作产生作用越大算法抗噪能力越强，由于数据平滑所造成数据的损失也越多。若拉格朗日乘子取值偏小，则正则化操作产生的作用越小对最终反演结果的影响就越小。
[0113][0114]
其中k为数据个数，trace为矩阵的迹即对角元素求和，i为单位矩阵，g为激光雷达探测的回波信号反演后得到的消光系数和后向散射系数。
[0115]
利用混合优化算法求解拉格朗日乘子。如图3所示，首先利用粒子群算法全局寻优，再利用遗传算法局部寻优，求得优中最优的拉格朗日乘子，继而得到权重因子。
[0116]
1)首先确定初始种群数量，利用pso算法进行全局寻优。根据最佳迭代概率与迭代次数的关系，确定最适合计算机计算效率的迭代次数，图4为确定初始种群数量的示意图。
[0117]
2)利用ga算法进行局部寻优，并进行选择操作，确定选择概率。
[0118]
3)利用ga算法进行局部寻优，并进行交叉操作，确定交叉概率，图5为确定交叉概率的示意图。
[0119]
4)利用ga算法进行局部寻优，并进行变异操作，确定变异概率，图6为确定变异概率的示意图。
[0120]
5)将gcv函数作为适应度函数，最佳适应度值为此式的最小值。
[0121]
s5：根据沙尘气溶胶粒子谱分布反演函数确定确定大气中的沙尘气溶胶粒子谱分布。
[0122]
s6：利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。
[0123]
结合大气中的沙尘气溶胶粒子谱分布和近地面沙尘气溶胶粒子谱分布，完成沙尘气溶胶在垂直方向上数量关系的监测。图8为不同时间大气和地面沙尘气溶胶粒子谱分布示意图。
[0124]
此外，本发明在步骤s1之前，还包括：计算偏振激光雷达探测大气中的沙尘气溶胶的回波信号的退偏振比；根据退偏振比确定沙尘气溶胶为非球形形态。图7为退偏振比验证沙尘粒子非球形特性的示意图。
[0125]
退偏振比定义为垂直和平行于激光束的平面上的后向回波原始信号的比值，是非球形粒子的一个重要物理参量。退偏振比通常用于表示粒子的非球形程度，数值越大非球形程度越高，数值越趋于0越接近球形。由于本发明研究主要为沙尘气溶胶粒子，沙尘气溶胶粒子是典型的非球形粒子，退偏比能较好反映沙尘气溶胶粒子的非球形程度和对光的散射情况。
[0126]
对于某一特定高度处或特定天气条件下的沙尘气溶胶粒子，可将偏振通道测得的垂直回波信号强度与平行方向上的回波信号相比，则退偏比表示为：
[0127]
[0128]
p
rs
、p
rp
分别为范围z内垂直偏振模式和平行偏振模式下的接收强度。
[0129]
实施例二
[0130]
为了执行上述实施例一对应的方法，以实现相应的功能和技术效果，下面提供了一种非球形气溶胶浓度谱探测分析系统。
[0131]
该系统包括：
[0132]
光学特性参数计算模块，用于利用dda算法计算沙尘气溶胶的光学特性参数；光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子；
[0133]
核函数计算模块，用于基于光学特性参数计算核函数；
[0134]
沙尘气溶胶光学数据计算模块，用于对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；光学数据包括消光系数和后散射系数；
[0135]
沙尘气溶胶粒子谱分布反演函数构建模块，用于基于核函数和光学数据构建沙尘气溶胶粒子谱分布反演函数；
[0136]
沙尘气溶胶粒子谱分布确定模块，用于根据沙尘气溶胶粒子谱分布反演函数确定确定大气中的沙尘气溶胶粒子谱分布；
[0137]
近地面沙尘气溶胶粒子谱分布确定模块，用于利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。
[0138]
实施例三
[0139]
本发明实施例三提供了一种电子设备，包括存储器及处理器，该存储器用于存储计算机程序，该处理器运行计算机程序以使电子设备执行实施例一的非球形气溶胶浓度谱探测分析方法
[0140]
上述电子设备可以是服务器。
[0141]
实施例四
[0142]
本发明实施例四提供了一种计算机可读存储介质，其存储有计算机程序，该计算机程序被处理器执行时实现实施例一的非球形气溶胶浓度谱探测分析方法。
[0143]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0144]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

技术特征：
1.一种非球形气溶胶浓度谱探测分析方法，其特征在于，包括：利用dda算法计算沙尘气溶胶的光学特性参数；所述光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子；基于所述光学特性参数计算核函数；对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；所述光学数据包括消光系数和后散射系数；基于所述核函数和所述光学数据构建沙尘气溶胶粒子谱分布反演函数；根据所述沙尘气溶胶粒子谱分布反演函数确定大气中的沙尘气溶胶粒子谱分布；利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。2.根据权利要求1所述的非球形气溶胶浓度谱探测分析方法，其特征在于，利用dda算法计算沙尘气溶胶的光学特性参数，具体包括：基于dda算法，将光学特性截面与等几何体截面之比定义为光学特性参数；其中，q
ext
表示消光效率因子，q
abs
表示吸收效率因子，q
sca
表示散射效率因子，c
ext
消光截面，c
abs
表示吸收截面，c
sca
表示散射截面，r表示气溶胶粒子半径。3.根据权利要求1所述的非球形气溶胶浓度谱探测分析方法，其特征在于，所述核函数的计算公式如下：其中，k表示核函数，q表示光学特性参数。4.根据权利要求1所述的非球形气溶胶浓度谱探测分析方法，其特征在于，基于所述核函数和所述光学数据构建气溶胶粒子谱分布反演函数，具体包括：利用混合优化算法求解最优拉格朗日乘子；所述混合优化算法包括粒子群算法和遗传算法；基于所述最优拉格朗日乘子和所述核函数确定权重因子；基于所述光学数据和所述权重因子构建沙尘气溶胶粒子谱分布反演函数。5.根据权利要求1所述的非球形气溶胶浓度谱探测分析方法，其特征在于，所述沙尘气溶胶粒子谱分布反演函数的表达式如下：其中，g表示沙尘气溶胶光学数据，w表示权重因子，f(r)表示沙尘气溶胶粒子谱分布函数，r表示气溶胶粒子半径。6.根据权利要求1所述的非球形气溶胶浓度谱探测分析方法，其特征在于，还包括：计算偏振激光雷达探测大气中的沙尘气溶胶的回波信号的退偏振比；根据所述退偏振比确定沙尘气溶胶为非球形形态。7.一种非球形气溶胶浓度谱探测分析系统，其特征在于，包括：光学特性参数计算模块，用于利用dda算法计算沙尘气溶胶的光学特性参数；所述光学特性参数包括：消光效率因子、吸收效率因子和散射效率因子；
核函数计算模块，用于基于所述光学特性参数计算核函数；沙尘气溶胶光学数据计算模块，用于对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；所述光学数据包括消光系数和后散射系数；沙尘气溶胶粒子谱分布反演函数构建模块，用于基于所述核函数和所述光学数据构建沙尘气溶胶粒子谱分布反演函数；沙尘气溶胶粒子谱分布确定模块，用于根据所述沙尘气溶胶粒子谱分布反演函数确定大气中的沙尘气溶胶粒子谱分布；近地面沙尘气溶胶粒子谱分布确定模块，用于利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。8.一种电子设备，其特征在于，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行权利要求1-6中任一项所述的非球形气溶胶浓度谱探测分析方法。9.一种计算机可读存储介质，其特征在于，其存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-6中任一项所述的非球形气溶胶浓度谱探测分析方法。

技术总结
本发明公开了一种非球形气溶胶浓度谱探测分析方法、系统及设备。该方法包括：利用DDA算法计算沙尘气溶胶的光学特性参数；基于所述光学特性参数计算核函数；对偏振激光雷达探测大气中的沙尘气溶胶的回波信号进行反演，得到沙尘气溶胶光学数据；基于所述核函数和所述光学数据构建沙尘气溶胶粒子谱分布反演函数；根据所述沙尘气溶胶粒子谱分布反演函数确定大气中的沙尘气溶胶粒子谱分布；利用粒径谱仪对近地面沙尘气溶胶进行探测，得到近地面沙尘气溶胶粒子谱分布。本发明能够解决非球形气溶胶粒子谱分布计算不准确的问题。粒子谱分布计算不准确的问题。粒子谱分布计算不准确的问题。


技术研发人员：乔泽 赵虎
受保护的技术使用者：乔泽
技术研发日：2023.06.28
技术公布日：2023/9/20