二、频率调制光谱（FMS）技术——error信号的数学模型及其他推论

一些推论。

目录

1 理论基础：线形与频率调制机制

2 相位免调节

3 一些推论

3.1 AFM、DFM信号峰值与调制频率的公式

3.2 AFM、DFM信号峰值与浓度的公式

附录1 AFM、DFM信号峰值、峰位的公式推导

---

1 理论基础：线形与频率调制机制

分子间碰撞导致能级扰动造成的展宽称之为碰撞展宽，通常用于描述气体压力导致的谱线加宽，该函数分布满足洛伦兹分布，可以表述为：

其中，v0 为跃迁中心频率，ΔvL 为碰撞展宽（半高全宽）。引入频率调制谱线形后，可归纳出吸收和色散分量线形。

将R(ω) 看作从频域-归一化频域的转换函数。它的特点有：在ω=ωc 时，Rω 为0；在ω=ωc±ΔvL/2 处，Rω 为1。

a. Rω

b. δ(ωc) 线形示意图

Fig1.洛伦兹归一化线形示意图(ωc=1500nm , ΔvL=3GHz )

将载波、边带的频率代入以上公式，得到δ 以及ϕ 分量各3项。例如：

 ​​​根据前文的推导，当ωm 大于吸收线宽时，吸收峰可以被单个边带探测。为便于理解，假设ωm 小于吸收线宽时，吸收峰也被载波及边带单独探测，也就是说生成三个独立分量。

δ 和ϕ 的下标，-1等表示边带等级。根据导数的定义，可以近似地认为Δδ ，Δϕ 分别是δ 、ϕ 分量的一阶、二阶导数。(Δϕ=(ϕ1-ϕ0)-(ϕ0-ϕ-1))

Fig3.吸收及色散分量、吸收信号AFM及色散信号DFM生成示意图（ωm=0.1ΔvL ）

当ωm <1.6*ΔvL/2 时，吸收信号在吸收中心频率表现为一道斜坡，而色散信号在任何调制频率都如此，均适合用于激光器在吸收中心的频率锁定。

2 相位免调节

对得到的信号进行IQ解调，得到一对正交的信号，令I信号的相位比Q信号超前90°：

根据式（12），IQ信号均含有吸收信号和色散信号。令吸收信号AFM为纯实数，色散信号DFM为纯虚数，得到如下的时钟图：

Fig4. IQ信号示意图

根据式（10）、（11），AFM、DFM信号的表达式如下：

AFM对ω 求偏导，AFM的极点频率则在其导数的零点上。定义极点频率差fAP 是极点频率与ωc 之差，与线宽ΔvL 均可与调制频率ωm 共同作为已知条件求出另一值。极点频率差fAP 与线宽ΔvL 均有4个解方程，两两共轭，仅列举各一解。附录1中说明了推导过程。

 

 

为确定IQ的相位，首先选取一解调相位测得一对正交的IQ信号。通过测算AFM信号的极点频率，可通过式(16)计算出该条件下的线宽及线型参数，代入DFM计算公式(15)，得到D’FM信号。可与第一步中测得的Q信号对比，从而判断解调相位的正确与否，得到相位差，进而实现相位免调节。

以下列出了一些不同初始解调相位的D’FM与Q信号的对比图。

 

Fig5. 不同初始解调相位的D’FM与Q信号的对比图（实线：反演计算得到的D’FM信号 虚线：实际测得的Q信号）

通过计算相关系数和ρx,y 以及RMSE，可以得到与实际相位差的相关程度。计算公式如下：

   

以下是不同初始解调相位的D’FM与Q信号的相关系数图及绝对误差和示意图。（令正确解调相位为0）

 

Fig6. 不同初始解调相位的D’FM与Q信号的ρx,y 及RMSE（正确解调相位为0）

能够看出，解调相位正确时，ρx,y 趋近于1而RMSE趋近于0。单一ρx,y 或RMSE信息都无法判断相位，而综合二者即可。

3 一些推论

3.1 AFM、DFM信号峰值与调制频率的公式

Fig6. 仿真error峰值（1654nm,50000ppm甲烷）考虑解调相位变化

Fig7. 实际测得error峰值（1654nm,50000ppm甲烷）不改变解调相位

结论如下。

		公式
AFM	峰位
	峰值
DFM	峰位
	峰值	过长略

Fig.8 载波1654nm，线宽0.002nm，error峰值随调制频率变化图

结论：当调制频率增大时，AFM峰值不断趋近于载波信号峰值*M/2，也就是边带信号峰值。载波波长越长、线宽越小，拐点越小（仅与这2个变量相关）。

3.2 AFM、DFM信号峰值与浓度的公式

Fig.9 载波1654nm，AFM信号峰值随线宽变化图

Fig.9 载波1654nm，AFM信号峰值随调制频率变化图

结论：AFM信号峰值仅与浓度成正比，M仅改变其斜率，DFM信号峰值与浓度成正比且M影响线形。

附录1 AFM、DFM信号峰值、峰位的公式推导

已知线宽、调制频率、激光器中心频率、吸收峰位置、载波峰值，求AFM(吸收)与DFM(色散)信号的峰值。AFM (吸收)与DFM(色散)峰值比与哪些参数有关。

问题：载波峰值是否可认为是DAS

% 已知线宽、调制频率求峰位(A)
syms fAP wm dv
y1=1./(dv^2+(fAP-wm).^2);
y2=1./(dv^2+(fAP+wm).^2);
x=diff(y1-y2,fAP);
a=solve(x==0,fAP);
pretty(a)
输出结果：
 
% 已知线宽、调制频率求峰值 (A)
peakvalue=subs(y,a(1));
latex(peakvalue);
输出结果：
 

% 已知线宽、调制频率求峰位(D)
syms xRw wm dv
y1=1./(dv^2+(xRw-wm).^2);
y2=1./(dv^2+(xRw+wm).^2);
y3=xRw./(dv^2+(xRw).^2);
y=y1+y2-2*y3;

x=diff(y1-y2,xRw);
a=solve(x==0,xRw);
pretty(a)

输出结果：
 
%对吸收峰位公式进行验算
peak_loc1=real(eval(subs(peak_loc(1),[deltav ,wm ,wc, M, cpeak ],[0.0566, 8.8e8 ,1391 ,2 ,320 ])))
%线宽0.0566nm，调制频率880MHZ，吸收中心频率1391nm，调制比2，载波吸收峰320mv，数据来源于陈健论文
输出结果：
ans= 1390.9

PeakOfErrorA_1=real(eval(subs(PeakOfErrorA,[deltav ,wm ,wc, M, cpeak ],[0.0566, 8.8e8 ,1391 ,2 ,320 ])))  %求 AFM峰值
输出结果：
ans= 20.8029

验算结果表明推测的公式符合实验结果。

数据来源：陈健_频率调制光谱的燃烧产物CO检测技术研究 图3-6