一种步进式比相方法、装置、设备及存储介质与流程

1.本发明属于等离子体电子密度参数测量技术领域，具体涉及一种步进式比相方法、装置、设备及存储介质。


背景技术：

2.在磁约束核聚变等离子体研究中，电子密度作为重要的等离子体参数，不仅能反映等离子体的约束和输运状况，还关系到装置的密度反馈控制。目前，有多种诊断技术可以测量等离子体电子密度参数，包括激光和微波干涉仪，汤姆逊散射等。大部分微波/激光干涉仪和部分偏振仪通常使用迈克尔逊或马赫贞德式光路结构，通过调频式调制解调方法，将测量道和参考道信号调制为呈调制频率的正弦波，并将所测量的信息转化为这两个信号间的相位差。
3.为提取该相位差，计算出所求的测量数据，传统的相位比较方法是快速傅里叶比相。这种方法首先将测量道和参考道同时做傅里叶变换得到两个频域信号，然后将这两个频域信号后半区域全部置零，随后做傅里叶反变换得到两个复信号。这两个复信号的辐角之差即是所求相位差曲线。
4.为了适应未来高密度聚变装置的需求，激光干涉仪正朝波长更短的方向发展。色散干涉仪是一种测量范围大、可靠性高的新型干涉仪。这种干涉仪采用外差式调频式调制时，因为调制元器件的性能限制，调制频率通常会高达40mhz以上，相应的采样率也会更高，数据系统需要处理的数据量会非常大，而采用传统的相位比较方法进行数据处理耗时会达到数分钟，为此亟需开发一种速度更快的数据处理方法。


技术实现要素：

5.为了适应高测量范围、高可靠性要求的新型干涉仪对数据处理数据的需求，本发明提供了一种步进式比相方法，本发明能够快速准确解算出干涉仪所测量的相位曲线，抗干扰能力强。
6.本发明通过下述技术方案实现：
7.一种步进式比相方法，包括：
8.获取干涉仪系统的测量道信号和参考道信号；
9.对参考道信号进行预设幅度的延迟，得到第二参考信号；
10.分别将测量道信号与参考道信号、测量道信号与第二参考信号进行相关运算，并进行滤波处理得到两相关信号；
11.将两相关信号逐点相除并做反正切计算得到相位曲线，去除基线后即可解算出密度相位。
12.优选的，本发明的在获取干涉仪系统的测量道信号和参考道信号步骤之后还包括：
13.对获取的测量道信号和参考道信号进行滤波处理，具体采用带通滤波器以滤除噪
声信号，所述带通滤波器的通带中心频率为调制频率。
14.优选的，本发明的分别将测量道信号与参考道信号、测量道信号与第二参考信号进行相关运算，具体为：
15.将测量道信号与参考道信号进行逐点相乘得到第一相关信号；
16.将测量道信号与第二参考信号进行逐点相乘得到第二相关信号。
17.优选的，本发明的进行滤波处理得到两相关信号，具体为：
18.通过数字滤波或平滑处理滤除逐点相乘后的第一相关信号和第二相关信号中角频率高于调制频率的成分，即可得到两相关信号。
19.优选的，本发明的平滑处理具体采用n点平滑处理，n取一个调制周期的采样点数量。
20.优选的，本发明的将两相关信号逐点相除并做反正切计算得到相位曲线，具体为：
21.将两相关信号逐点相除并进行反正切运算，得到原始相位曲线；
22.根据测量道信号与参考道信号的相关信号，所述原始相位曲线进行象限修正；
23.对经象限修正后的曲线进行条纹修正，即可得到完整的相位曲线。
24.优选的，本发明的象限修正具体为将提取的相位差取值范围从(-π/2，+π/2)扩充到(-3π/2，+π/2)，即一整个条纹的区间内。
25.第二方面，本发明提出了一种步进式比相装置，包括：
26.信号获取模块，用于获取干涉仪系统的测量道信号和参考道信号；
27.延迟模块，用于对所述参考道信号进行预设幅度的延迟，得到第二参考信号；
28.相关运算模块，分别将测量道信号与参考道信号、测量道信号与第二参考信号进行相关运算，并进行滤波处理得到两相关信号；
29.解算模块，用于将两相关信号逐点相除并做反正切计算得到相位曲线，去除基线后即可解算出密度相位。
30.第三方面，本发明提出了一种电子设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现本发明所述方法的步骤。
31.第四方面，本发明提出了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现本发明所述方法的步骤。
32.本发明具有如下的优点和有益效果：
33.1、本发明可快速比较干涉仪测量信号和参考信号的相位差，从而解算出干涉仪所测量的的相位曲线，解算速度快，相位计算精确，抗干扰能力强。
34.2、本发明适用于各种调频调制的外差式激光干涉仪，尤其适用于高密度等离子体装置的外差式色散干涉仪。本发明已在软件模拟实验和中国环流器2号a(hl-2a)托卡马克装置甲酸激光(hcooh，波长为432.5μm)偏振/干涉仪上实现。
附图说明
35.此处所说明的附图用来提供对本发明实施例的进一步理解，构成本技术的一部分，并不构成对本发明实施例的限定。在附图中：
36.图1为本发明实施例的比相方法流程示意图。
37.图2为本发明实施例的电子设备结构示意图。
38.图3为本发明实施例的比相装置原理框图。
39.图4为采用本发明测量线性变化的密度时各信号变化状况。
40.图5为采用本发明在极低性噪比下的相位测量模拟实验结果。
41.图6为采用本发明和传统的快速傅里叶比相的相位测量模拟实验对比结果。
42.图7为采用本发明和传统的快速傅里叶比相在hl-2a托卡马克装置第35188号放电实验中，对甲酸激光干涉仪测量数据的相位提取对比结果。
具体实施方式
43.为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
44.实施例1
45.传统的快速傅里叶比相技术，处理速度无法满足更高性能的激光干涉仪的需求，基于此，本实施例提出了一种步进式比相方法，该方法能够快速进行正弦信号的相位比较，适用于各种跳频调制的外差式激光干涉仪，尤其适用于高密度等离子体装置的外差式色散干涉仪。
46.具体如图1所示，本实施例的步进式比相方法具体包括：
47.步骤s1，获取干涉仪系统的测量道信号s和参考道信号r。
48.本实施例中，在干涉仪系统中，探测器测量得到的干涉信号(即测量道信号s)，其相位可表示为：
49.i
∝
cos(ωmt+φ
p
+ψ)
50.其中，ωm为调制频率，φ
p
为等离子体产生的密度相移，ψ为干涉信号的初始相位。
51.参考道信号r，其相位可表示为：
52.i
r1
∝
cos(ωmt+ψ2)
53.其中，ψ2为参考信号的初始相位。
54.本实施例还可对获取的信号进行滤波处理，以进一步提高处理精度，具体为采用带通滤波器以滤除常规的噪声信号，通带中心频率为调制频率ωm。
55.步骤s2，对参考道信号r进行预设幅度的延迟，得到第二参考信号r2。
56.其中，本实施例中的延迟幅度为(k
±
1/4)个调制周期，其中，k为整数；本实施例优选为1/4个调制周期。
57.第二参考信号r2，其相位可表示为：
58.i
r2
∝
sin(ωmt+ψ2)
59.步骤s3，分别将测量道信号s与参考道信号r、测量道信号s与第二参考信号r2进行相关运算，并进行数字滤波或平滑处理得到两相关信号i1和i2。
60.本实施例将两个参考信号和测量信号进行相关运算(逐点相乘)，得到两相关信号，其相位表示为：
61.ii
r1
∝
[cos(2ωmt+φ
p
+ψ+ψ2)+cos(φ
p
+ψ-ψ2)]
[0062]
ii
r2
∝
[sin(2ωmt+φ
p
+ψ+ψ2)-sin(φ
p
+ψ-ψ2)]
[0063]
通过数字滤波或平滑处理滤掉角频率高于ωm的成分，就可以得到两相关信号i1和
i2，其相位表示为：
[0064][0065][0066]
其中，本实施例的平滑处理具体为：进行n点平滑，n取一个调制周期的采样点数量。
[0067]
步骤s4，将两相关信号i1和i2逐点相除并做反正切计算得到相位曲线，去除基线后就可以解算出密度相位φ
p
。
[0068]
本实施例中，将两相关信号i1和i2逐点相除并做反正切计算得到相位曲线具体为：
[0069]
将两相关信号逐点相除并进行反正切运算，得到提取的原始相位曲线；
[0070]
根据i1的正负符号，对原始相位曲线做象限修正，将提取的相位差取值范围从(-π/2，+π/2)扩充到(-3π/2，+π/2)，即一整个条纹的区间内。
[0071]
对象限修正后的曲线进行条纹修正即可得到完整的提取相位曲线。
[0072]
本实施例在获得密度相位后，即可计算出线密度积分
[0073]
采用上述方法测量线性变化的密度时各信号变化情况如图2所示。其中，图2(a)为参与干涉的两探测束电矢量曲线；图2(b)为两探测束干涉后由探测器得到的信号曲线(即测量道信号和参考道信号)，该曲线主要频率成分和参考信号cos2ωmt一致。图2(c)为计算得到的两相关曲线经滤波后的结果，滤波去除了频率超过2ωm的信号成分。图2(d)为最终得到的相位曲线，该曲线与预设密度曲线完全平行，且信噪比良好，可以证明本实施例的方法能够正确解算出所测量得到的相位曲线。
[0074]
采用上述方法在极低信噪比下的相位测量模拟实验结果如图3所示，测量密度为恒定的常数，测量信号包含极高的干扰成分，信噪比仅有1.8db，原始信号已难以辨认，但采用本实施例提出的方法仍能计算出相位曲线，且不会出现跳变或无法回归基线等测量错误，此条件下相位分辨率降低值6
°
。可以证明，本实施例的方法的适用范围更广，抗干扰能力强。
[0075]
本实施例还提出了一种电子设备(计算机设备)，用于执行本实施例的上述方法。
[0076]
具体如图4所示，本实施例的电子设备包括处理器、内存储器和系统总线；内存储器和处理器在内的各种设备组件连接到系统总线上。处理器是一个用来通过计算机系统中基本的算术和逻辑运算来执行计算机程序指令的硬件。内存储器是一个用于临时或永久性存储计算程序或数据(例如，程序状态信息)的物理设备。系统总线可以为以下几种类型的总线结构中的任意一种，包括存储器总线或存储控制器、外设总线和局部总线。处理器和内存储器可以通过系统总线进行数据通信。其中内存储器包括只读存储器(rom)或闪存(图中未示出)，以及随机存取存储器(ram)，ram通常是指加载了操作系统和计算机程序的主存储器。
[0077]
电子设备一般包括一个外存储设备。外存储设备可以从多种计算机可读介质中选择，计算机可读介质是指可以通过计算机设备访问的任何可利用的介质，包括移动的和固定的两种介质。例如，计算机可读介质包括但不限于，闪速存储器(微型sd卡)，cd-rom，数字通用光盘(dvd)或其它光盘存储、磁带盒、磁带、磁盘存储或其它磁存储设备，或者可用于存
储所需信息并可由计算机设备访问的任何其它介质。
[0078]
电子设备可在网络环境中与一个或者多个网络终端进行逻辑连接。网络终端可以是个人电脑、服务器、路由器、智能电话、平板电脑或者其它公共网络节点。计算机设备通过网络接口(局域网lan接口)与网络终端相连接。局域网(lan)是指在有限区域内，例如家庭、学校、计算机实验室、或者使用网络媒体的办公楼，互联组成的计算机网络。wifi和双绞线布线以太网是最常用的构建局域网的两种技术。
[0079]
应当指出的是，其它包括比电子设备更多或更少的子系统的计算机系统也能适用于发明。
[0080]
如上面详细描述的，适用于本实施例的电子设备能执行步进式比相方法的指定操作。电子设备通过处理器运行在计算机可读介质中的软件指令的形式来执行这些操作。这些软件指令可以从存储设备或者通过局域网接口从另一设备读入到存储器中。存储在存储器中的软件指令使得处理器执行上述的群成员信息的处理方法。此外，通过硬件电路或者硬件电路结合软件指令也能同样实现本发明。因此，实现本实施例并不限于任何特定硬件电路和软件的组合。
[0081]
实施例2
[0082]
本实施例提出了一种步进式比相装置，具体如图5所示，本实施例的装置包括：
[0083]
信号获取模块，用于获取或读入测量道信号s和参考道信号r。
[0084]
本实施例的信号获取模块还能够对获取的信号进行滤波处理，具体采用带通滤波器以除去常规的噪声信号，通带中心频率为调制频率ωm。
[0085]
延迟模块，用于对获取的参考道信号r进行预设幅度的延迟，得到第二参考信号r2。
[0086]
相关运算模块，用于分别将测量道信号s与参考道信号r、测量道信号s与第二参考信号r2进行相关运算，并进行数字滤波或平滑处理得到两相关信号i1和i2。
[0087]
解算模块，将两相关信号i1和i2逐点相除并做反正切计算得到提取相位曲线，去除基线后就可以解算出密度相位φ
p
。
[0088]
本实施例在获得密度相位φ
p
后，即可计算出线密度积分
[0089]
实施例3
[0090]
本实施例为验证上述实施例提出的步进式比相技术相对于传统傅里叶快速比相的优势，分别进行了软件模拟对比实验和干涉仪实际测量对比实验。
[0091]
具体为：
[0092]
模拟实验中，使用两种不同的比相算法对干涉信号i
∝
cos(2ωmt+3k/2ω+ψ)进行解算。干涉信号总长度为0.1s，采样率f＝200mhz，调制频率ωm＝40mhz，信号中加入信噪比为50db的白噪声。模拟实验的测量曲线见图6。步进式比相耗时约为3.38s，快速傅里叶比相耗时约为19.30s，在计算速度方面步进式比相有明显优势。步进式比相的相位分辨率为0.21
°
，快速傅里叶比相的相位分辨率为0.88
°
，步进式比相的相位分辨率比起传统比相方式，在高噪声环境下也具有优势，抗噪能力更强。
[0093]
干涉仪测量对比实验中，使用甲酸激光干涉仪在hl-2a托卡马克装置(该装置采用双激光偏振/干涉仪光路结构设计，包含一套完整的激光干涉测量光路，可通过两台甲酸激光器的腔长调节，由锁频系统获得可调的稳定差频)#35188号放电的测量数据进行相位提
取运算，原始信号的采样率为6.25mhz，调制频率ωm约为1.46mhz。测量实验的结果见图7。步进式比相耗时约为1.49s，快速傅里叶比相耗时约为7.89s，步进式比相同样体现出了明显的优势。通过测量曲线的局部放大，可以确认步进式比相能够正确测量密度的快速扰动，测量的时间分辨率不低于传统的比相方式。
[0094]
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。