CN104749396A - 一种基于波长调制-tdlas的高速气流流速检测系统 - Google Patents

Abstract

本发明公开了一种基于波长调制-TDLAS的高速气流流速检测系统，包括激光器模块、高速气流模块、信号采集和数据处理模块以及上位机，激光器模块发射的两路激光位于高速气流流速方向的同一侧，一路顺着高速气流流速方向，一路逆着高速气流流速方向，均和高速气流流速不成90°夹角，位于高速气流后方的信号采集和数据处理模块接收并处理高速气流的吸收光谱二次谐波信号，得到高速气流信息。本发明操作简单、测量精度高，能消除背景信号、激光强度波动和光电放大系数等因素的影响，解决了目前高速气流中测量需要直接接触气体源导致气流受到干扰的难题，能够用于高温、高流速、振动等恶劣环境的检测，从而拓宽了WMS技术的应用范围。

Description

一种基于波长调制-TDLAS的高速气流流速检测系统

技术领域

本发明涉及气流测量领域，特别涉及一种基于波长调制-TDLAS的高速气流流速检测系统。

背景技术

目前国内外吸气式发动机的流速测量方法主要有四种方法，1是利用皮托管基于流体力学的流量方程来测量流速，2是利用热线/热膜测速仪器利用放置在流场中具有加热电流的细金属丝(直径0.5～5μm)来测量流速的仪器，3是利用激光多普勒粒径分析仪利用流体中运动微粒散射光的多普勒频移来获得流体速度信息，利用4是利用使用粒子图像测速仪(PIV),皮托管对流场的干扰较大，热线风速仪有无法弥补的原理性缺陷,如热丝长细比、需要标定等,,对被测流场产生也会产生干扰，激光多普勒测速仪的最大优点在于测量的非接触性,但强调它是非接触式测量的同时必须说明其测量结果并不是流体流动的速度,而是示踪颗粒的速度采用机械式仪表进行直接测量，容易对流场造成干扰，测量精度也受到一定影响，PIV是目前发展最快的测量流速方法，但是其测量的结果仍然只是颗粒的运动速度,也需要示踪粒子的问题。

因此，迫切需求一种非接触式不干扰流场、无需示踪粒子的新型测量手段。

TDLAS(Tunable Diode Laser Absorption Spectroscopy)是近些年发展起来的、先进的、可直接测量气体参数的检测技术，该技术利用波长可调谐的窄带激光扫描气体分子特征吸收谱线，通过分析吸收光谱谐波信号中心频率的偏移，得到待测气体速度，因此具有极高的波长选择性、高灵敏度、系统通用性等优点。

TDLAS技术主要有直接吸收光谱技术和波长调制光谱技术为主的两种测量方法。直接吸收光谱技术通过入射光强和透射光强的比值直接拟合气体吸收率函数，进而通过吸收率函数测量气体温度，但直接吸收在测量中容易受到颗粒物浓度、激光强度波动等因素的影响而无法精确拟合气体吸收率函数，进而导致温度测量误差。而波长调制光谱技术(Wavelength Modulation Spectroscopy，WMS)在测量过程中通过对目标信号进行高频调制，而非目标信号由于没有经过调制在后续的谐波检测过程中被除去，因此可以有效地降低测量系统中背景信号的干扰，极大地提高温度测量精度和灵敏度，尤其适用于恶劣环境或弱吸收条件下气体温度的在线测量

基于波长调制TDLAS技术的流速测量方法，利用吸收气体的二次谐波信号来反演气体流速信息，能够提取微弱的吸收信号、具有很高的信噪比，适宜于高速气流气源环境的测量具有非接触测量、动态响应快、测量精度高、测量量程大、易小型化等特点。

发明内容

本发明的目的为了克服现有高速气体测量的缺陷，解决高速气流发生装置中气体流速测量需要直接接触气体或采用示踪粒子的问题，同时为提供一种可提高即时测定的信赖性的装置，提出一种非接触式，无需示踪粒子的一种基于波长调制-TDLAS的高速气流流速检测的简易系统。

本发明采用波长调制可调谐激光二极管吸收光谱技术，利用多普勒频移计算气体流速算法，可以对特定气体进行有选择性地测定的气体流速测定系统。

本发明的具体技术方案如下：

一种基于波长调制-TDLAS的高速气流流速检测系统，包括有激光器模块、高速气流模块、信号采集和数据处理模块以及上位机，其特征在于：所述激光器模块包括激光器、激光器驱动电路、控制激光器驱动电路的信号调制和扫描电路、加法电路；还包括设置在激光器输出激光光路上的分束器、两路准直器；所述信号采集和数据处理模块包括位于准直器焦点上的探测器、依次连接探测器的放大电路、信号采集电路和数据处理模块，还包括信号采集控制模块；所述信号调制和扫描电路、激光器驱动电路以及信号采集电路连接信号采集控制模块；所述数据处理模块位于上位机软件内；所述信号调制和扫描电路的输出信号经加法电路处理后作用到激光器驱动电路上，激光器驱动电路控制激光器输出激光，激光器输出激光经分束器和两路准直器器后，位于高速气流流速方向的同一侧，一路顺着高速气流流速方向，一路逆着高速气流流速方向穿过气流，均和高速气流流速不成90°夹角，之后被信号采集和数据处理模块接收，信号采集和数据处理模块接收并处理通过高速气流得到的高速气流吸收光谱二次谐波信号，得到高速气流信息。

还包括连接在采集电路的标准具，所述标准具对激光器扫描波长进行定标。

所述信号调制和扫面电路包括信号发生器和函数信号发生器，用于向加法电路输出选定频率、幅值、相位的正弦波调制信号和锯齿波扫描信号；所述加法电路将调制扫描电路的信号相加；所述激光器驱动电路包括温度控制电路和电流驱动电路，采用ITC102激光器驱动器，该激光驱动器利用加法电路输出信号控制激光器的输出激光波长与频率，所述激光器为Nanoplus公司的DFB激光器，其输出激光为探测氧气用的波长761nm激光，或探测水汽用的波长1392nm激光。

所述高速气流模块为超音速风洞，所述超音速风洞的试验段设置有一个观测窗口，所述观测窗口的两侧分别放置有准直器和探测器。

所述的分束器为两等分功率光分束器；

所述的采集电路采用德国Spectrum公司M2i.4021采集卡，该采集卡为12位位宽，最高采用速率20Msps。

信号采集和数据处理模块接收并处理高速气流吸收光谱的二次谐波信号，其工作流程如下

(1)信号采集模块采集探测器探测到的两路高速气流的吸收光谱信号，一路为待测信号，一路为参考信号，在上位机控制模块中设置信号采集模块的基本参数；

(2)数据处理模块对采集的高速气流的吸收光谱二次谐波信号进行预处理并结合数字锁相算法提取二次谐波信号，首先对采集的待测信号进行累加平均、小波滤波处理等去噪处理，再利用数字锁相算法提取二次谐波信号，数据采集控制模块中设置参考信号频率，幅值，初始相位，数字滤波器的参数；

(3)采用最小二乘拟合对提取得到二次谐波信号进行数据拟合；

(4)高速气流流速信息反演：计算拟合后的两路二次谐波的中心频率υ₁、υ₂及其差值|υ₁-υ₂|，根据多普勒频移测量流速原理，计算出高速气流信息。

本发明工作时，信号发生器产生调制和扫描两路信号，扫描信号为频率100Hz、峰峰值为1.2V的锯齿波；调制信号为频率50kHz、峰峰值为40mv的正弦波，调制信号和扫描信号叠加后，送入电流驱动电路，电流驱动电路输出信号接入DFB激光器，激光器输出波长根据调制和扫描信号进行变化，激光器输出信号通过两等分功率光分束器，接入到位于高速气流装置实验段的两路准直器中，准直器光路在观察窗口侧与气流成一定夹角，准直器输出的激光束透过气流进入探测器，探测器接收到吸收光谱信号后输出到放大电路，放大电路对探测器输出的电压信号放大到合适的观察范围，再利用数据采集卡对放大器输出信号进行采集，采集的数据信号送到上位机数据处理软件进行在线处理或者后续离线处理。

所述数据处理的软件分为前面板界面和后台程序。

本发明所述波长调制原理为：

当一束强度为I₀，频率为ν的激光，通过光程为L的气体池，输出光强I与I₀关系为：

I(w)＝I₀(w)exp(-σ(w)LN)                (1)

式中σ(w)为吸收系数，L为吸收光程长度，N为吸收气体的浓度。

当激光器的中心频率ω_L受到频率为ω_m的正弦波调制时，其瞬时频率可表示为：

ω＝ω_L+δωcos(ω_mt)               (2)

式中δω是调制幅度。激光束通过样吸收池后，其强度可以用I(ω_L)的余弦傅里叶级数来表示：

$I({\mathrm{\ω}}_L,t)=\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{A}_n\left({\mathrm{\ω}}_L\right)\cos \left(n{\mathrm{\ω}}_{m}t\right)---\left(3\right)$

I(ω_L)每个谐波分量用锁相放大器可以测得为：

$A_n\left({\mathrm{\ω}}_L\right)\frac{2}{\mathrm{\π}}{\∫}_0^{\mathrm{\π}}{I}_0({\mathrm{\ω}}_L+\mathrm{\δ\ω}\cos \mathrm{\θ})\exp \{-\mathrm{\σ}({\mathrm{\ω}}_L+\mathrm{\δ\ω}\cos \mathrm{\θ})\mathrm{LN}\}\cos \left(\mathrm{n\θ}\right)\mathrm{d\θ}---\left(4\right)$

式中θ＝ω_mt。

对于ppm量级的弱吸收气体而言，σLN<<1，不考虑光强调制效应,I₀(ω_L+δωcosθ)≈I(ω_L)≈I₀，可得到：

$A_n\left({\mathrm{\&omega;}}_L\right)\frac{2I_0\mathrm{NL}}{\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}-\mathrm{\&sigma;}({\mathrm{\&omega;}}_L+\mathrm{\&delta;\&omega;}\cos \mathrm{\&theta;})\mathrm{LN}\cos \left(\mathrm{n\&theta;}\right)\mathrm{d\&theta;}---\left(5\right)$

对上式中的σ(ω)在ω_L处进行泰勒级数展开,整理(5)式得到：

$A_n\left({\mathrm{\&omega;}}_L\right)=\frac{2^{1-n}{I}_0\mathrm{NL}}{n!}\mathrm{\&delta;}{\mathrm{\&omega;}}^n\frac{d^n\mathrm{\&sigma;}}{d{\mathrm{\&omega;}}^2}{|}_{\mathrm{\&omega;}={\mathrm{\&omega;}}_L}---\left(6\right)$

可以看出n次谐波的幅值正比于吸收系数σ(ω)的n次导数、原始光强I₀、吸收光程长度L以及吸收气体的浓度N。

任意次谐波分量与气体浓度成正比，因此都可以用来表示气体的浓度等特性，实际应用时，一般多使用二次谐波，二次谐波分量如式(7)所示：

$A_2\left({\mathrm{\&omega;}}_L\right)=\frac{I_0\mathrm{NL}}{4}\mathrm{\&delta;}{\mathrm{\&omega;}}^2\frac{d^2\mathrm{\&sigma;}}{d{\mathrm{\&omega;}}^2}{|}_{\mathrm{\&omega;}={\mathrm{\&omega;}}_L}---\left(7\right)$

本发明所述流速测量原理如下：

根据激光多普勒原理，当气流流速在激光传输方向上有速度分量时，激光接收频率与实际激光发生频率之间有一定的偏移，二者频率之差即为多普勒频移，其频移公式为：

$\mathrm{\&Delta;\&upsi;}=\frac{V_{\mathrm{gas}}\&times;\cos \mathrm{\&theta;}}{c}\&times;{\mathrm{\&upsi;}}_0---\left(8\right)$

其中，υ₀为无频移时气体吸收线中心频率，θ为激光束方向和气流速度方向的夹角，c为光速，V_gas为待测气体流速。

实际测量中，往往采用双光路方法，如附图3所示，利用吸收谱线的二次谐波的频移来反演气体流速。双光路可以抵消实验系统中由于环境因素变化导致的频移偏差。双光路测量气体流速计算公式为：

$V_{\mathrm{gas}}=\frac{\mathrm{c\&Delta;}{\mathrm{\&upsi;}}_{12}}{{\mathrm{\&upsi;}}_0(\cos {\mathrm{\&theta;}}_1-\cos {\mathrm{\&theta;}}_2)}---\left(9\right)$

其中θ₁，θ₂分别是两路激光束与气流的夹角。

本发明的有益效果体现在：

本系统能够实现PPM量级的高速气流的在线实时快速测量；同时，本系统在工作的时候采用的非接触测量，不会干扰高速气流气源；再者，本系统工作的时候无需示踪粒子，无需高速相机，大大节约了设备成本。

附图说明

附图1为本发明的结构示意图；

附图2为本发明的信号采集和数据处理模块工作流程图；

附图3为本发明的双路激光与高速气流作用示意图。

具体实施方式

如图1所示，一种基于波长调制-TDLAS的高速气流流速检测系统，包括有激光器模块、高速气流模块、信号采集和数据处理模块以及上位机，其特征在于：所述激光器模块包括激光器、激光器驱动电路、控制激光器驱动电路的信号调制和扫描电路、加法电路；还包括设置在激光器输出激光光路上的分束器、两路准直器；所述信号采集和数据处理模块包括位于准直器焦点上的探测器、依次连接探测器的放大电路、信号采集电路和数据处理模块，还包括信号采集控制模块；所述信号调制和扫描电路、激光器驱动电路以及信号采集电路连接信号采集控制模块；所述数据处理模块位于上位机软件内；所述信号调制和扫描电路的输出信号经加法电路处理后作用到激光器驱动电路上，激光器驱动电路控制激光器输出激光，激光器输出激光经分束器和两路准直器器后，如图3所示，位于高速气流流速方向的同一侧，一路顺着高速气流流速方向，一路逆着高速气流流速方向穿过气流，均和高速气流流速不成90°夹角，之后被信号采集和数据处理模块接收，信号采集和数据处理模块接收并处理通过高速气流得到的高速气流吸收光谱二次谐波信号，得到高速气流信息。

还包括连接在采集电路的标准具，所述标准具对激光器扫描波长进行定标。

所述信号调制和扫面电路包括信号发生器和函数信号发生器，用于向加法电路输出选定频率、幅值、相位的正弦波调制信号和锯齿波扫描信号；所述加法电路将调制扫描电路的信号相加；所述激光器驱动电路包括温度控制电路和电流驱动电路，采用ITC102激光器驱动器，该激光驱动器利用加法电路输出信号控制激光器的输出激光波长与频率，所述激光器为Nanoplus公司的DFB激光器，其输出激光为探测氧气用的波长761nm激光，或探测水汽用的波长1392nm激光。

所述高速气流模块为超音速风洞，所述超音速风洞的试验段设置有一个观测窗口，所述观测窗口的两侧分别放置有准直器和探测器。

所述的分束器为两等分功率光分束器；

所述的采集电路采用德国Spectrum公司M2i.4021采集卡，该采集卡为12位位宽，最高采用速率20Msps。

如图2所示，所述信号采集和数据处理模块接收并处理高速气流吸收光谱的二次谐波信号，其工作流程如下

(1)信号采集模块采集探测器探测到的两路高速气流的吸收光谱信号，一路为待测信号，一路为参考信号，在上位机控制模块中设置信号采集模块的基本参数；

(2)数据处理模块对采集的高速气流的吸收光谱二次谐波信号进行预处理和并结合数字锁相算法提取二次谐波信号，首先对采集的待测信号进行累加平均、小波滤波处理等去噪处理，再利用数字锁相算法提取二次谐波信号，数据采集控制模块中设置参考信号频率，幅值，初始相位，数字滤波器的参数。

(3)采用最小二乘拟合对提取得到二次谐波信号进行数据拟合；

(4)高速气流流速信息反演：计算拟合后的两路二次谐波的中心频率υ₁、υ₂及其差值|υ₁-υ₂|，根据多普勒频移测量流速原理，计算出高速气流信息。

本发明工作时，信号发生器产生调制和扫描两路信号，扫描信号为频率100Hz、峰峰值为1.2V的锯齿波；调制信号为频率50kHz、峰峰值为40mv的正弦波，调制信号和扫描信号叠加后，送入电流驱动电路，电流驱动电路输出信号接入DFB激光器，激光器输出波长根据调制和扫描信号进行变化，激光器输出信号通过两等分功率光分束器，接入到位于高速气流装置实验段的两路准直器中，准直器光路在观察窗口侧与气流成一定夹角，准直器输出的激光束透过气流进入探测器，探测器接收到吸收光谱信号后输出到放大电路，放大电路对探测器输出的电压信号放大到合适的观察范围，再利用数据采集卡对放大器输出信号进行采集，采集的数据信号送到上位机数据处理软件进行在线处理或者后续离线处理。

所述数据处理的软件分为前面板界面和后台程序。

本发明所述波长调制原理为：

当一束强度为I₀，频率为ν的激光，通过光程为L的气体池，输出光强I与I₀关系为：

I(w)＝I₀(w)exp(-σ(w)LN)              (1)

式中σ(w)为吸收系数，L为吸收光程长度，N为吸收气体的浓度。

当激光器的中心频率ω_L受到频率为ω_m的正弦波调制时，其瞬时频率可表示为：

ω＝ω_L+δωcos(ω_mt)                (2)

式中δω是调制幅度。激光束通过样吸收池后，其强度可以用I(ω_L)的余弦傅里叶级数来表示：

$I({\mathrm{\&omega;}}_L,t)=\underset{n=0}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{A}_n\left({\mathrm{\&omega;}}_L\right)\cos \left(n{\mathrm{\&omega;}}_{m}t\right)---\left(3\right)$

I(ω_L)每个谐波分量用锁相放大器可以测得为：

$A_n\left({\mathrm{\&omega;}}_L\right)\frac{2}{\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}{I}_0({\mathrm{\&omega;}}_L+\mathrm{\&delta;\&omega;}\cos \mathrm{\&theta;})\exp \{-\mathrm{\&sigma;}({\mathrm{\&omega;}}_L+\mathrm{\&delta;\&omega;}\cos \mathrm{\&theta;})\mathrm{LN}\}\cos \left(\mathrm{n\&theta;}\right)\mathrm{d\&theta;}---\left(4\right)$

式中θ＝ω_mt。

对于ppm量级的弱吸收气体而言，σLN<<1，不考虑光强调制效应,I₀(ω_L+δωcosθ)≈I(ω_L)≈I₀，可得到：

$A_n\left({\mathrm{\&omega;}}_L\right)\frac{2I_0\mathrm{NL}}{\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}-\mathrm{\&sigma;}({\mathrm{\&omega;}}_L+\mathrm{\&delta;\&omega;}\cos \mathrm{\&theta;})\mathrm{LN}\cos \left(\mathrm{n\&theta;}\right)\mathrm{d\&theta;}---\left(5\right)$

对上式中的σ(ω)在ω_L处进行泰勒级数展开,整理(5)式得到：

$A_n\left({\mathrm{\&omega;}}_L\right)=\frac{2^{1-n}{I}_0\mathrm{NL}}{n!}\mathrm{\&delta;}{\mathrm{\&omega;}}^n\frac{d^n\mathrm{\&sigma;}}{d{\mathrm{\&omega;}}^2}{|}_{\mathrm{\&omega;}={\mathrm{\&omega;}}_L}---\left(6\right)$

可以看出n次谐波的幅值正比于吸收系数σ(ω)的n次导数、原始光强I₀、吸收光程长度L以及吸收气体的浓度N。

任意次谐波分量与气体浓度成正比，因此都可以用来表示气体的浓度等特性，实际应用时，一般多使用二次谐波，二次谐波分量如式(7)所示：

$A_2\left({\mathrm{\&omega;}}_L\right)=\frac{I_0\mathrm{NL}}{4}\mathrm{\&delta;}{\mathrm{\&omega;}}^2\frac{d^2\mathrm{\&sigma;}}{d{\mathrm{\&omega;}}^2}{|}_{\mathrm{\&omega;}={\mathrm{\&omega;}}_L}---\left(7\right)$

本发明所述流速测量原理如下：

根据激光多普勒原理，当气流流速在激光传输方向上有速度分量时，激光接收频率与实际激光发生频率之间有一定的偏移，二者频率之差即为多普勒频移，其频移公式为：

$\mathrm{\&omega;\&upsi;}=\frac{V_{\mathrm{gas}}\&times;\cos \mathrm{\&theta;}}{c}\&times;{\mathrm{\&upsi;}}_0---\left(8\right)$

其中，υ₀为无频移时气体吸收线中心频率，θ为激光束方向和气流速度方向的夹角，c为光速，V_gas为待测气体流速。

实际测量中，往往采用双光路方法，如附图3所示，利用吸收谱线的二次谐波的频移来反演气体流速。双光路可以抵消实验系统中由于环境因素变化导致的频移偏差。双光路测量气体流速计算公式为：

$V_{\mathrm{gas}}=\frac{\mathrm{c\&Delta;}{\mathrm{\&upsi;}}_{12}}{{\mathrm{\&upsi;}}_0(\cos {\mathrm{\&theta;}}_1-\cos {\mathrm{\&theta;}}_2)}$

其中θ₁，θ₂分别是两路激光束与气流的夹角。υ₁₂为两路二次谐波吸收信号的频移差值。

Claims (7)

1.一种基于波长调制-TDLAS的高速气流流速检测系统，包括有激光器模块、高速气流模块、信号采集和数据处理模块以及上位机，其特征在于：所述激光器模块包括激光器、激光器驱动电路、控制激光器驱动电路的信号调制和扫描电路、加法电路；还包括设置在激光器输出激光光路上的分束器、两路准直器；所述信号采集和数据处理模块包括位于准直器焦点上的探测器、依次连接探测器的放大电路、信号采集电路和数据处理模块，还包括信号采集控制模块；所述信号调制和扫描电路、激光器驱动电路以及信号采集电路连接信号采集控制模块；所述数据处理模块位于上位机软件内；所述信号调制和扫描电路的输出信号经加法电路处理后作用到激光器驱动电路上，激光器驱动电路控制激光器输出激光，激光器输出激光经分束器和两路准直器器后，位于高速气流流速方向的同一侧，一路顺着高速气流流速方向，一路逆着高速气流流速方向穿过气流，均和高速气流流速不成90°夹角，之后被信号采集和数据处理模块接收，信号采集和数据处理模块接收并处理通过高速气流得到的高速气流吸收光谱二次谐波信号，得到高速气流信息。

2.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速检测系统，其特征在于：还包括连接信号采集电路的标准具，对激光器扫描波长进行标定。

3.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速检测系统，其特征在于：所述信号调制和扫描电路包括函数信号发生器，用于向加法电路输出选定频率、幅值、相位的正弦波调制信号和锯齿波扫描信号；所述激光器驱动电路包括温度控制电路和电流驱动电路，采用ITC102激光器驱动器，该激光驱动器利用加法电路输出信号控制激光器的输出激光波长与频率，所述激光器为Nanoplus公司的DFB激光器，其输出激光为探测氧气用的波长761nm激光，或探测水汽用的波长1392nm激光。

4.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速检测系统，其特征在于：所述高速气流模块为超音速风洞，所述超音速风洞的试验段设置有一个观测窗口，所述观测窗口的两侧分别放置有准直器和探测器。

5.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速检测系统，其特征在于：所述的分束器为两等分功率光分束器。

6.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速检测系统，其特征在于：所述的信号采集电路采用德国Spectrum公司M2i.4021采集卡，该采集卡为12位位宽，最高采用速率20Msps。

7.根据权利要求1所述的基于波长调制-TDLAS的高速气流流速系统，其特征在于：所述信号采集和数据处理模块接收并处理高速气流吸收光谱的二次谐波信号的工作流程如下：

(1)信号采集模块采集探测器探测到的两路高速气流的吸收光谱信号，一路为待测信号，一路为参考信号，在上位机控制模块中设置信号采集模块的基本参数；

(2)数据处理模块对采集的高速气流的吸收光谱二次谐波信号进行预处理和并结合数字锁相算法提取二次谐波信号，首先对采集的待测信号进行累加平均、小波滤波处理等去噪处理，再利用数字锁相算法提取二次谐波信号，数据采集控制模块中设置参考信号频率，幅值，初始相位，数字滤波器的参数。

(3)采用最小二乘拟合对提取得到二次谐波信号进行数据拟合；

(4)高速气流流速信息反演：计算拟合后的两路二次谐波的中心频率υ₁、υ₂及其差值|υ₁-υ₂|，根据多普勒频移测量流速原理，计算出高速气流信息。