光纤法布里珀罗传感器游标效应增敏（Matlab代码实现）

游标效应增敏原理

由于FP腔可以产生周期性的余弦光谱，当再经过若干级次之后又会重新对齐。因此，级联FP腔的光谱同样具有周期性，而且级联FP腔光谱的自由光谱范围和两个FP腔的自由光谱范围相关。当其中一个FP腔的光谱发生微弱漂移而另一个FP腔的光谱没有发生变化时，通过检测原始对齐的波峰（或波谷）和当前对齐的波峰（或波谷）的位置，就可以将微小的变化进行放大，则相比于单个FP腔，级联FP腔的输出光谱图会发生更大的变化。

游标效应实现方法

将两个腔长近似的FP型干涉仪级联或并联，其中一个腔为传感腔，用来测量待测压力。另一个腔为参考腔，腔长与传感腔长相似，用作参考光谱。

公式表达

令入射电场强度为1，总的反射强度可以表示为从两个反射面反射回来的电场强度的和，表达式为：

其中，a1是第一个反射面的透射损耗；

          𝑅1 是第一个反射面的反射率；

          𝑅2是 第二个面的反射率；

          θ=2𝜋𝑛𝐿/𝜆是腔内光束传播的相位移动。

因为反射光强是电场强度振幅的平方，所以对应的反射光强根据公式表示为:

为了方便matlab代码编写，参考腔的反射光强表示为：

传感腔的反射光强表示为： 

总反射光强为：

干涉光谱

设置传感腔长为300um，辅助腔长为280um。通过matlab仿真等到的干涉叠加光谱如下图1，对光谱中极大值组成的包络进行拟合，得到的包络曲线如下图2，通过记录包络曲线的漂移量可分析出传感量的变化。

游标效应增敏分析

通过对比单个传感器腔长变化时干涉光谱前后发生的变化验证游标效应增敏效果，当传感腔腔长从300um 改变为301um，单腔的腔长变化前后的干涉光谱（上）和双腔仿真的干涉光谱（下）。

其中，传感腔长增加1um后，单腔的仿真干涉光谱红移了2nm，双腔仿真干涉光谱红移了约20nm。

结论

本F-P压力传感器应用游标效应增敏可放大传感器灵敏度理论值约10倍左右。（根据传感腔长和参考腔长的调整可改变放大的具体数值）

游标效应Matlab代码

单腔传感matlab代码

liamud = linspace(1450,1460,10000); %波长
R1=0.04; %第一个反射面的反射率
R2=0.08; %第二个反射面反射率
a1 = 0.01; %第一个反射面的投射损耗
n1 = 1;  %折射率
L1 = 300000; %腔长 默认单位nm
h1 = (2.*pi.*n1.*L1)./liamud;
I1 = R1 + (1-R1)^2.*(1-a1)^2.*R2 + 2.*(1-a1).*(1-R1).*(R1.*R2)^0.5.*cos(2.*h1);
%传感腔发生形变后
L3 = 301000;
h3 = (2.*pi.*n1.*L3)./liamud;
I3 = R1 + (1-R1)^2.*(1-a1)^2.*R2 + 2.*(1-a1).*(1-R1).*(R1.*R2)^0.5.*cos(2.*h3);
figure;  
plot(liamud, I1);          %游标效应干涉光谱
hold on;
plot(liamud, I3);       %传感腔形变后，干涉光谱的变化情况
% title('干涉光谱');     
xlabel('波长'); % 可以根据实际情况修改x轴标签  
ylabel('反射光强');

双腔游标matlab代码

liamud = linspace(1350,1550,10000); %波长
R1=0.04; %参考腔：第一个反射面的反射率
R2=0.08; %        第二个反射面反射率
R3=0.04; %传感腔：第一个反射面的反射率
R4=0.08; %        第二个反射面的反射率
a1 = 0.01; %参考腔：第一个反射面的投射损耗
a2 = 0.01; %传感腔：第一个反射面的投射损耗
n1 = 1;  %参考腔折射率
n2 = 1;  %传感腔折射率
L1 = 280000; %参考腔长 默认单位nm
L2 = 300000; %传感腔长 默认单位nm
%辅助腔的反射光强
h1 = (2.*pi.*n1.*L1)./liamud;
I1 = R1 + (1-R1)^2.*(1-a1)^2.*R2 + 2.*(1-a1).*(1-R1).*(R1.*R2)^0.5.*cos(2.*h1);
%传感腔的反射光强
h2 = (2.*pi.*n2.*L2)./liamud;
I2 = R3 + (1-R3)^2.*(1-a2)^2.*R4 + 2.*(1-a2).*(1-R3).*(R3.*R4)^0.5.*cos(2.*h2);
%总反射光强
IR = I1+I2;
%传感腔发生形变后
L3 = 301000;
h3 = (2.*pi.*n2.*L3)./liamud;
I3 = R3 + (1-R3)^2.*(1-a2)^2.*R4 + 2.*(1-a2).*(1-R3).*(R3.*R4)^0.5.*cos(2.*h3);
IR2 = I1+I3;
% 找到IR的极大值点  
[peaks,locs] = findpeaks(IR);  
% 生成函数图像  
figure;  
plot(liamud, IR, liamud(locs) ,peaks ,'linewidth',1);  %干涉光谱包络函数
% plot(liamud, IR);          %游标效应干涉光谱
hold on;
% plot(liamud, IR2);       %传感腔形变后，干涉光谱的变化情况
% title('干涉光谱');     
xlabel('波长'); % 可以根据实际情况修改x轴标签  
ylabel('反射光强');