一种对半共格两相衍射曲线进行拆分的方法

1.本公开属于衍射技术领域，具体涉及一种对半共格两相衍射曲线进行拆分的方法。


背景技术：

2.半共格的双相结构材料通常表现出优异的物理性能、机械性能等，为了满足特定服役工况，工程领域往往设计具备半共格晶体学关系的两相结构材料。与此同时，半共格的特点也决定了材料两相的错配度极小，以第三代沉淀强化的镍基高温合金为例，室温下的错配度一般小于百分之一。性能由结构决定，为了不断优化半共格结构材料的各方面表现，其微观结构的精细表征必不可少，常用手段为衍射技术，例如x射线衍射技术，中子衍射技术等。困难地是，半共格两相的错配度极小，造成两相的衍射曲线高度重合，无法直接辨别，这对精细表征而言是极大的挑战。
3.现在已有的曲线拆分方法，是基于衍射曲线在强度峰值两侧的不对称性，获得半共格两相各自的衍射曲线。该方法提供了最基础的拆分思路，具备借鉴意义，但同时存在很大的局限性，其在于：该方法只考虑衍射曲线的不对称性，未考虑半共格两相在结构材料中的体积分数以及衍射结构因子(体积分数以及衍射结构因子会直接影响半共格两相的衍射强度积分之比，即衍射曲线的积分面积之比)。所以也必将导致通过该方法拆分得到的衍射曲线不能合理真实地反映出双相材料的本征结构特点，这与精细表征的初衷背道而驰，甚至会对结构的设计造成误导。因此，急需开发一种新颖的更加合理的曲线拆分方法，用来处理半共格两相结构材料的衍射曲线。


技术实现要素：

4.针对现有技术中的不足，本公开的目的在于提供一种对半共格两相衍射曲线进行拆分的方法，本方法综合考虑了半共格两相的相对体积分数以及衍射结构因子，相比仅依靠曲线特征对半共格两相衍射曲线进行拆分的方法，本方法得到的半共格两相的衍射曲线能够合理真实地反应出双相材料的本征特点。
5.为实现上述目的，本公开提供以下技术方案：
6.一种对半共格两相衍射曲线进行拆分的方法，包括如下步骤：
7.s100：通过对金属材料的微观组织进行形貌表征，获得金属材料中半共格两相a和b的相对体积分数va和vb，且va+vb＝1；
8.s200：分别计算a相和b相的衍射结构因子fa和fb；
9.s300：通过衍射获得a相和b相的衍射矢量模长-衍射强度曲线xq-yi；
10.s400：对所述衍射矢量模长-衍射强度曲线xq-yi进行拟合，以获得衍射矢量模长-衍射强度拟合曲线xq-yi；
11.s500：根据镜面法则对所述衍射矢量模长-衍射强度拟合曲线xq-yi进行拆分，以获得若干组衍射曲线；
12.s600：根据a相和b相的相对体积分数va和vb以及衍射结构因子fa和fb，从若干组衍射曲线中选出a相和b相的衍射曲线。
13.优选的，步骤s200中，所述衍射结构因子fa通过下式计算：
[0014][0015]
所述衍射结构因子fb通过下式计算：
[0016][0017]
其中，n是基本单元所占据的不同位置，n＝1，2，
…
，n，n是组成晶胞的基本单元所占据的位置总个数；π是圆周率；i是数学中的虚数符号；hkl是米勒指数，表示发生衍射的晶面为(hkl)，且h、k、l的取值范围均为区间[0，∞)内的整数；un、vn、wn分别是基本单元所占据的位置n的横坐标、纵坐标和竖坐标，且un、vn、wn的取值范围均为区间[0，1]内的实数；fn是基本单元的散射因子。
[0018]
优选的，所述散射因子fn通过下式计算：
[0019][0020]
其中，m是构成不同位置n的最小单元，且m＝1，2，
…
，m，m是最小单元的个数；f
n，m
是最小单元m的散射因子；是最小单元m在不同位置n出现的概率，且满足关系式
[0021]
优选的，步骤s500包括如下步骤：
[0022]
s501：将所述衍射矢量模长-衍射强度拟合曲线xq-yi中所有数值的横坐标组成一维等差向量xq，xq相邻元素的差值记为δxq；
[0023]
s502：以所述衍射矢量模长-衍射强度拟合曲线xq-yi的峰值的横坐标xq
max
为中心，以l为半径，生成具有2l+1个元素的一维等差向量xql，xql相邻元素的差值等于δxq；
[0024]
s503：在所述衍射矢量模长-衍射强度拟合曲线xq-yi上，以向量xql的第l个元素xq
l
为横坐标对应的数据点(xq
l
，yi
l
)及(xq
l
，0)所在的线段为镜面，对所述衍射矢量模长-衍射强度拟合曲线xq-yi位于(xq
l
，yi
l
)右侧的部分进行镜面对称操作，即得到初次对称曲线xq
l_rr-yi
l_rr
；
[0025]
s504：将所述衍射矢量模长-衍射强度拟合曲线xq-yi减去初次对称曲线xq
l_rr-yi
l_rr
，即得到初次差值曲线xq
l_r-‑
yi
l_r-；
[0026]
s505：以初次差值曲线xq
l_r-‑
yi
l_r-的峰值以及峰值在xq轴的投影点所在线段为镜面，对初次差值曲线xq
l_r-‑
yi
l_r-位于峰值左侧的部分进行镜面对称操作，即得到二次对称曲线xq
l_ll-yi
l_ll
；
[0027]
s506：将所述衍射矢量模长-衍射强度拟合曲线xq-yi减去二次对称曲线xq
l_ll-yi
l_ll
，即得到二次差值曲线xq
l_l-‑
yi
l_l-，二次对称曲线xq
l_ll-yi
l_ll
和二次差值曲线xq
l_l-‑
yi
l_l-即为所述衍射矢量模长-衍射强度拟合曲线xq-yi的第l组衍射曲线；
[0028]
s507：对一维等差向量xql中的2l+1个元素重复执行步骤s503到步骤s506，即对所述衍射矢量模长-衍射强度拟合曲线xq-yi进行2l+1种拆分，即得到2l+1组二次对称曲线xq
l_ll-yi
l_ll
和二次差值曲线xq
l_l-‑
yi
l_l-。
[0029]
优选的，步骤s600包括如下步骤：
[0030]
s601：对2l+1组二次对称曲线xq
l_ll-yi
l_ll
进行积分运算，得到2l+1组积分结果，记为itg
l_ll
；对2l+1组二次差值曲线xq
l_l-‑
yi
l_l-进行积分运算，得到2l+1组积分结果，记为itg
l_l-；
[0031]
s602：根据itg
l_ll
和itg
l_l-计算得到2l+1组偏差值δ
l
；
[0032]
s603：选择2l+1组偏差值δ
l
中的最小值，记为δ
l
′
，则δ
l
′
所对应的二次对称曲线xq
l
′
_ll-yi
l
′
_ll
和二次不对称曲线xq
l
′
_l-‑
yi
l
′
_l-即为a相和b相的衍射曲线，分别记为xq
a-yia和xq
b-yib。
[0033]
优选的，步骤s602中，通过下式计算偏差值δ
l
：
[0034][0035]
其中，va和vb分别是a相和b相的相对体积分数，fa和fb分别是a相和b相的衍射结构因子。
[0036]
优选的，步骤s400中，采用包括多项式函数p(x)，gauss函数g(x)，lorentz函数l(x)和voigt函数v(x)在内的任一函数对所述衍射矢量模长-衍射强度曲线xq-yi进行拟合，其中，
[0037]
所述多项式函数p(x)表示为：
[0038]
p(x)＝anxn+a
n-1
x
n-1
+
…
+a2x2+a1x+a0[0039]
所述gauss函数g(x)表示为：
[0040][0041]
所述lorentz函数l(x)表示为：
[0042][0043]
所述voigt函数v(x)表示为：
[0044][0045]
其中，a0，a1，a2，a
n-1
，an，a，b，c，a
′
，b
′
，c
′
，σ，γ是函数中的系数，e是自然数，π是圆周率。
[0046]
优选的，本公开还提供一种计算机，包括：
[0047]
存储器，用于存储计算机指令；
[0048]
处理器，用于执行计算机指令以实现前述任一方法。
[0049]
与现有技术相比，本公开带来的有益效果为：
[0050]
1、在实施曲线拆分之前，通过对半共格两相的衍射矢量模长-衍射强度曲线进行拟合得到原则上分辨率无限高的衍射矢量模长-衍射强度拟合曲线，能够提高曲线拆分的
精确度；
[0051]
2、通过综合考虑半共格两相的相对体积分数以及衍射结构因子，能够合理真实地反应出双相材料的本征特点，从而能够克服仅仅依靠曲线特征进行拆分的不足。
附图说明
[0052]
图1是本公开一个实施例提供的一种对半共格两相衍射曲线进行拆分的方法的流程图；
[0053]
图2是本公开另一个实施例提供的通过扫描电子显微镜观察获得的镍基高温合金半共格两相a和b的形貌表征示意图；
[0054]
图3是本公开另一个实施例提供的通过x射线衍射获得的镍基高温合金半共格两相a和b的衍射矢量模长-衍射强度曲线xq-yi；
[0055]
图4是本公开另一个实施例提供的利用gauss函数对曲线xq-yi进行拟合获得的衍射矢量模长-衍射强度拟合曲线xq-yi；
[0056]
图5(a)、图5(b)是本公开另一个实施例提供的根据镜面法则对衍射矢量模长-衍射强度拟合曲线xq-yi进行拆分获得的衍射曲线，其中，图5(a)对应于l＝1201，图5(b)对应于l＝2001；
[0057]
图6是本公开另一个实施例提供的二次对称曲线xq
l_ll-yi
l_ll
对应的积分值itg
l_ll
以及二次差值曲线xq
l_l-‑
yi
l_l-对应的积分值itg
l_l-的结果示意图；
[0058]
图7是本公开另一个实施例提供的偏差值δl的结果示意图；
[0059]
图8是本公开另一个实施例提供的初次对称曲线xq
1417_rr-yi
1417_rr
和初次差值曲线xq
1417_r-‑
yi
1417_r-的结果示意图；
[0060]
图9是本公开另一个实施例提供的镍基高温合金半共格两相a和b的衍射曲线xq
a-yia和xq
b-yib；
[0061]
图10是本公开另一个实施例提供的利用衍射峰的不对称性对曲线进行拆分，所获得的镍基高温合金半共格两相a和b的衍射曲线。
具体实施方式
[0062]
下面将参照附图1至图10详细地描述本公开的具体实施例。虽然附图中显示了本公开的具体实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。
[0063]
需要说明的是，在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可以理解，技术人员可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名词的差异作为区分组件的方式，而是以组件在功能上的差异作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语，故应解释成“包含但不限定于”。说明书后续描述为实施本公开的较佳实施方式，然所述描述乃以说明书的一般原则为目的，并非用以限定本公开的范围。本公开的保护范围当视所附权利要求所界定者为准。
[0064]
为便于对本公开实施例的理解，下面将结合附图以具体实施例为例做进一步的解
释说明，且各个附图并不构成对本公开实施例的限定。
[0065]
一个实施例中，如图1所示，本公开提供一种对半共格两相衍射曲线进行拆分的方法，包括如下步骤：
[0066]
s100：通过对金属材料的微观组织进行形貌表征，获得金属材料中半共格两相a和b的相对体积分数va和vb，且va+vb＝1；
[0067]
s200：分别计算a相和b相的衍射结构因子fa和fb；
[0068]
s300：通过衍射获得a相和b相的衍射矢量模长-衍射强度曲线xq-yi；
[0069]
s400：对所述衍射矢量模长-衍射强度曲线xq-yi进行拟合，以获得衍射矢量模长-衍射强度拟合曲线xq-yi；
[0070]
s500：根据镜面法则对所述衍射矢量模长-衍射强度拟合曲线xq-yi进行拆分，以获得若干组衍射曲线；
[0071]
s600：根据a相和b相的相对体积分数va和vb以及衍射结构因子fa和fb，从若干组衍射曲线中选出a相和b相的衍射曲线。
[0072]
本实施例中，在曲线拆分之前，先对半共格两相的衍射矢量模长-衍射强度曲线进行拟合得到衍射矢量模长-衍射强度拟合曲线，以此提高曲线拆分的精确度。因为在衍射实验的过程中，为了达到节省时间的主观目的，或者由于实验条件限制的客观因素，采集到的衍射数据往往分辨率不够，影响最终的曲线拆分准确性。因此，先通过合适的数学函数对实验数据点即衍射矢量模长-衍射强度曲线进行拟合，得到原则上分辨率无限高的衍射矢量模长-衍射强度拟合曲线，然后再实施曲线拆分，最终大大地提高了曲线拆分的精确度。此外，本实施例在确定半共格两相的衍射曲线时，综合考虑了半共格两相的相对体积分数以及衍射结构因子，充分利用了两相的微观组织以及晶体结构和组成成分，因此得到的半共格两相的衍射曲线能够合理真实地反应出双相材料的本征特点，进而达到精细表征以及指导材料设计的目的。
[0073]
另一个实施例中，步骤s200中，所述衍射结构因子fa通过下式计算：
[0074][0075]
所述衍射结构因子fb通过下式计算：
[0076][0077]
其中，n是基本单元所占据的不同位置，n＝1，2，
…
，n，n是组成晶胞的基本单元所占据的位置总个数；π是圆周率；i是数学中的虚数符号；hkl是米勒指数，代表发生衍射的晶面为(hkl)，且h、k、l的取值范围均为区间[0，∞)内的整数；un、vn、wn分别是基本单元所占据的位置n的横坐标、纵坐标和竖坐标，且un、vn、wn的取值范围均为区间[0，1]内的实数；fn是基本单元的散射因子。
[0078]
另一个实施例中，所述散射因子fn通过下式计算：
[0079]
[0080]
其中，m是构成不同位置n的最小单元，且m＝1，2，
…
，m，m是最小单元的个数；f
n，m
是最小单元m的散射因子；是最小单元m在不同位置n出现的概率，且满足关系式
[0081]
另一个实施例中，步骤s500包括如下步骤：
[0082]
s501：将所述衍射矢量模长-衍射强度拟合曲线xq-yi中所有数值的横坐标组成一维等差向量xq，xq相邻元素的差值记为δxq；
[0083]
s502：以所述衍射矢量模长-衍射强度拟合曲线xq-yi的峰值的横坐标xq
max
为中心，以l为半径，生成具有2l+1个元素的一维等差向量xql，xql相邻元素的差值等于δxq；
[0084]
s503：在所述衍射矢量模长-衍射强度拟合曲线xq-yi上，以向量xql的第l个元素xq
l
为横坐标对应的数据点(xq
l
，yi
l
)及(xq
l
，0)所在的线段为镜面，对所述衍射矢量模长-衍射强度拟合曲线xq-yi位于(xq
l
，yi
l
)右侧的部分进行镜面对称操作，即得到初次对称曲线xq
l_rr-yi
l_rr
；
[0085]
s504：将所述衍射矢量模长-衍射强度拟合曲线xq-yi减去初次对称曲线xq
l_rr-yi
l_rr
，以得到初次差值曲线xq
l_r-‑
yi
l_r-；
[0086]
s505：以初次差值曲线xq
l_r-‑
yi
l_r-的峰值以及峰值在xq轴的投影点所在线段为镜面，对初次差值曲线xq
l_r-‑
yi
l_r-位于峰值左侧的部分进行镜面对称操作，以得到二次对称曲线xq
l_ll-yi
l_ll
；
[0087]
s506：将所述衍射矢量模长-衍射强度拟合曲线xq-yi减去二次对称曲线xq
l_ll-yi
l_ll
，以得到二次差值曲线xq
l_l-‑
yi
l_l-，二次对称曲线xq
l_ll-yi
l_ll
和二次差值曲线xq
l_l-‑
yi
l_l-即为所述衍射矢量模长-衍射强度拟合曲线xq-yi的第l组衍射曲线；
[0088]
s507：对一维等差向量xql中的2l+1个元素重复执行步骤s503到步骤s506，即对所述衍射矢量模长-衍射强度拟合曲线xq-yi进行2l+1种拆分，即得到2l+1组二次对称曲线xq
l_ll-yi
l_ll
和二次差值曲线xq
l_l-‑
yi
l_l-。
[0089]
另一个实施例中，步骤s600包括如下步骤：
[0090]
s601：对2l+1组二次对称曲线xq
l_ll-yi
l_ll
进行积分运算，得到2l+1组积分结果，记为itg
l_ll
；对2l+1组二次差值曲线xq
l_l-‑
yi
l_l-进行积分运算，得到2l+1组积分结果，记为itg
l_l-；
[0091]
s602：根据itg
l_ll
和itg
l_l-计算得到2l+1组偏差值δ
l
，所述偏差值δ
l
通过下式计算：
[0092][0093]
其中，va和vb分别是半共格两相a和b的相对体积分数，fa和fb分别是半共格两相a和b的衍射结构因子。
[0094]
s603：选择2l+1组偏差值δ
l
中的最小值，记为δ
l
′
，下角标l
′
是位于区间[1，2l+1]的整数，则δ
l
′
所对应的二次对称曲线xq
l
′
_ll-yi
l
′
_ll
和二次不对称曲线xq
l
′
_l-‑
yi
l
′
_l-即为a相和b相的衍射曲线，分别记为xq
a-yia和xq
b-yib。
[0095]
下面，以镍基高温合金的x射线衍射为例，对本公开的技术方案作具体描述。
[0096]
1、针对镍基高温合金的微观组织，电化学腐蚀后利用扫描电子显微镜获得如图2所示的半共格两相中a相和b相的形貌表征，通过观察可获得半共格两相中a相和b相的相对体积分数分别为：va＝30％，vb＝70％，且va+vb＝30％+70％＝1。
[0097]
2、分别计算a相和b相的衍射结构因子fa和fb，其中，
[0098]
a相的衍射结构因子fa的计算公式为：
[0099][0100]
在镍基高温合金中，a相具有面心立方晶体结构：晶胞的基本单元所占据的位置的总个数为4，即n＝4；基本单元所占据的位置n的坐标为(un，vn，wn)，当n＝1时，(u1，v1，w1)＝(0，0，0)，当n＝2时，(u2，v2，w2)＝(0，1/2，1/2)，当n＝3时，(u3，v3，w3)＝(1/2，0，1/2)，当n＝4时，(u4，v4，w4)＝(1/2，1/2，0)；发生衍射的晶面指数(hkl)是(800)，即h＝8，k＝0，l＝0，其中，hkl是米勒指数，表示发生衍射的晶面为(hkl)，h，k，l的取值范围均是区间[0，∞)内的整数；f1，f2，f3，f4分别是占据n＝1，n＝2，n＝3，n＝4四个位置的基本单元的散射因子fn，fn的计算公式表示为：
[0101][0102]
在镍基高温合金中，a相的fn满足f1＝f2＝f3＝f4。占据位置n的基本单元由总计三种最小单元构成，即m＝3；m＝1的最小单元为ni原子，f
n，1
＝f
ni
＝28.3372+0.7725i，m＝2的最小单元为co原子，f
n，2
＝f
co
＝27.3187+0.6608i，m＝3的最小单元为cr原子，f
n，3
＝f
cr
＝24.2626+0.4175i；是最小单元m在位置n出现的概率，f
ni
，f
co
，f
cr
分别是ni原子的散射因子，co原子的散射因子，cr原子的散射因子。
[0103]
将以上变量带入镍基高温合金中a相的衍射结构因子计算公式：
[0104][0105]
最终得到fa＝107.69+2.57i。
[0106]
同理，b相的衍射结构因子fb的计算公式为：
[0107][0108]
在镍基高温合金中，b相具有l12晶体结构：晶胞的基本单元所占据的位置的总个数为4，即n＝4；基本单元所占据的位置n的坐标(un，vn，wn)，当n＝1时，(u1，v1，w1)＝(0，0，
0)，当n＝2时，(u2，v2，w2)＝(0，1/2，1/2)，当n＝3时，(u3，v3，w3)＝(1/2，0，1/2)，当n＝4时，(u4，v4，w4)＝(1/2，1/2，0)；发生衍射的晶面指数(hkl)是(800)；f1，f2，f3，f4分别是占据n＝1，n＝2，n＝3，n＝4四个位置的基本单元的散射因子fn，fn计算公式为：
[0109][0110]
在镍基高温合金中，b相的fn满足f1≠f2＝f3＝f4。占据位置n＝1的基本单元总计两种最小单元构成，即m＝2，占据位置n＝2，3，4的基本单元仅由一种最小单元构成；当n＝1，m＝1时，最小单元为al原子，f
1，1
＝f
al
＝13.0435+0.0337i，当n＝1，m＝2时，最小单元为ti原子，f
1，2
＝f
ti
＝22.2247+0.3060i，当n＝2，3，4时，占据位置n的最小单元为ni原子，f2＝f3＝f4＝f
ni
＝28.3372+0.7725i；是最小单元m在位置n出现的概率，f
al
是al原子的散射因子，f
ti
是ti原子的散射因子。
[0111]
将以上变量带入镍基高温合金中b相的衍射结构因子计算公式：
[0112][0113]
最终得到fb＝100.99+2.44i。
[0114]
3、通过包括但不限于x射线衍射、中子衍射和电子衍射在内的任一x射线衍射手段，获得如图3所示的a相和b相在晶面(800)的衍射矢量模长-衍射强度曲线xq-yi，其中，曲线xq-yi的横坐标范围是[139.51，141.31]，精度为0.03，单位为nm-1
。
[0115]
4、利用包括多项式函数p(x)、gauss函数g(x)、lorentz函数l(x)和voigt函数v(x)在内的任一函数对衍射矢量模长-衍射强度曲线xq-yi进行拟合，获得如图4所示的衍射矢量模长-衍射强度拟合曲线xq-yi，其中，
[0116]
多项式函数p红)表示为：
[0117]
p(x)＝anxn+a
n-1
x
n-1
+...+a2x2+a1x+a0[0118]
所述gauss函数g(x)表示为：
[0119][0120]
所述lorentz函数l(x)表示为：
[0121][0122]
所述voigt函数v(x)表示为：
[0123]
[0124]
其中，a0，a1，a2，a
n-1
，an，a，b，c，a
′
，b
′
，c
′
，σ，γ是函数中的系数，e是自然数，π是圆周率。
[0125]
下面，以gauss函数g(x)为例对衍射矢量模长-衍射强度曲线xq-yi进行拟合：
[0126][0127]
上式中，a，b，c为待定系数，等式两边取自然对数，得到：
[0128][0129]
将衍射矢量模长-衍射强度曲线xq-yi上的所有数据点(xj，yj)带入上述等式，得到：
[0130][0131]
根据广义最小二乘解，得到系数a，b，c，函数g(x)的表达式由此确定。其中，s是衍射矢量模长-衍射强度曲线xq-yi上所有数据点的总个数，j是区间为[1，s]内的整数。
[0132]
5、根据镜面法则对衍射矢量模长-衍射强度拟合曲线xq-yi进行拆分，具体步骤如下：
[0133]
5.1、将衍射矢量模长-衍射强度拟合曲线xq-yi中所有数值的横坐标组成一维等差向量xq，xq＝[139.5100：0.0001：141.3100]，向量xq相邻元素的差值δxq＝0.0001(δxq＝0.0001是人为规定的步长)，单位为nm-1
；
[0134]
5.2、以衍射矢量模长-衍射强度拟合曲线xq-yi的峰值的横坐标xq
max
＝140.4033nm-1
为中心，以l＝1000为半径(l是人为规定的一个范围，代表以曲线峰值为基准，向左取l个数据点，向右取l个数据点，因此称l为半径，且l的取值范围小于所述xq长度二分之一的正整数)，生成具有2l+1＝2001个元素的一维等差向量xql，且向量xql相邻元素之间的差值等于δxq＝0.0001nm-1
，即向量xql的第l个元素xql(l)＝xq
max-(l-1-l)
×
δxq＝140.4033-(l-1001)
×
0.0001，l是位于区间[1，2001]的整数，xql＝[140.3033：0.0001：140.5033]，单位为nm-1
；
[0135]
5.3、在所述衍射矢量模长-衍射强度拟合曲线xq-yi上，以向量xql的第l个元素xq
l
为横坐标对应的数据点(xq
l
，yi
l
)及(xq
l
，0)所在的线段为镜面，对所述衍射矢量模长-衍射强度拟合曲线xq-yi位于(xq
l
，yi
l
)右侧的部分进行镜面对称操作，即得到初次对称曲线xq
l_rr-yi
l_
rr；
[0136]
5.4、将衍射矢量模长-衍射强度拟合曲线xq-yi减去初次对称曲线xq
l_rr-yi
l_rr
，即得到初次差值曲线xq
l_r-‑
yi
l_r-；
[0137]
5.5、以初次差值曲线xq
l_r-‑
yi
l_r-的峰值以及峰值在xq轴的投影点所在线段为镜面，对初次差值曲线xq
l_r-‑
yi
l_r-位于峰值左侧的部分进行镜面对称操作，即得到二次对称曲线xq
l_ll-yi
l_ll
；
[0138]
5.6、将所述衍射矢量模长-衍射强度拟合曲线xq-yi减去二次对称曲线xq
l_ll-yi
l_ll
，即得到二次差值曲线xq
l_l-‑
yi
l_l-，二次对称曲线xq
l_ll-yi
l_ll
和二次差值曲线xq
l_l-‑
yib的积分面积之比则为48.71％。理论的半共格两相a和b的衍射曲线积分面积之比为显而易见，根据本公开所提供的方法进行曲线拆分时，所获得的半共格两相a和b的衍射曲线积分面积之比更接近理论值，因此也更能反应双相材料的本征结构特点。
[0149]
另一个实施例中，本公开还提供一种计算机，包括：
[0150]
存储器，用于存储计算机指令；
[0151]
处理器，用于执行计算机指令以实现前述实施例任一方法。
[0152]
以上结合具体实施例对本公开所提供的技术方案进行了详细介绍，同时，上述实施例的说明只是用于帮助理解本公开的方法及其核心思想，对于本领域的一般技术人员，依据本公开的思想，在具体实施方式及应用范围上均会有改变之处。因此，本说明书内容不应理解为对本公开的限制。