基于MATLAB的光纤布拉格光栅反射谱特性仿真系统

资 源 简 介

本项目基于MATLAB平台开发，旨在通过数值计算模拟光纤布拉格光栅（Fiber Bragg Grating, FBG）的反射光谱特性。项目包含两段独立但互补的核心代码，采用耦合模理论（Coupled-mode Theory）与传输矩阵法（Transfer Matrix Method, TMM）作为数学模型。第一段代码专注于均匀光纤光栅的理论仿真，用于计算理想状态下的反射率、透射率与波长的关系，分析光栅长度、有效折射率调制深度及光栅周期对主峰带宽和峰值反射率的影响；第二段代码进阶模拟非均匀光栅或环境变化下的光谱响应，例如切趾（Apodization）光纤光栅以抑制旁瓣，或模拟在不同温度和应变加载条件下布拉格中心波长的漂移情况。系统最终通过图形化界面或二维绘图直观展示波长与反射率的非线性曲线，能够精确计算3dB带宽、零点带宽及旁瓣抑制比，适用于光通信滤波器设计、光纤传感机理研究及光电子学教学演示。

详 情 说 明

基于MATLAB的光纤布拉格光栅(FBG)反射谱特性仿真系统

项目介绍

本项目是一个基于MATLAB平台开发的光学仿真系统，专门用于数值模拟光纤布拉格光栅（FBG）的反射光谱特性。核心算法结合了耦合模理论（Coupled-mode Theory）与传输矩阵法（Transfer Matrix Method, TMM），能够精确计算光在光栅结构中的传播行为。

系统通过图形化界面直观展示了均匀光栅与切趾光栅的光谱差异，并模拟了环境温度变化引起的中心波长漂移现象，最终自动计算传感灵敏度，适用于光通信器件设计与光纤传感原理的教学演示。

功能特性

• 均匀光栅光谱仿真：计算并绘制理想均匀FBG的反射率曲线，自动分析峰值反射率、中心波长及FWHM（半高全宽/3dB带宽）。
• 变迹（切趾）效应分析：支持高斯（Gaussian）变迹模拟，用于对比分析切趾技术对旁瓣（Side-lobes）的抑制效果。
• 环境传感特性模拟：基于热光效应和热膨胀效应，模拟不同温度条件下（25°C, 50°C, 75°C）布拉格中心波长的漂移。
• 线性度与灵敏度分析：自动追踪不同环境下的波长峰值，进行线性拟合计算，并输出温度传感灵敏度（pm/°C）。
• 可视化数据展示：在一个综合图表中通过四个子图展示反射谱、对比谱、漂移谱及拟合曲线。

系统要求

• MATLAB R2016a 或更高版本
• 无需额外工具箱（仅使用基础数学与绘图函数）

核心算法与实现细节

本项目采用严格的物理模型进行仿真，而非简单的公式拟合。

1. 传输矩阵法 (TMM)

为了解决非均匀或分段光栅的计算问题，代码将整个光栅长度（L = 10mm）在轴向z上分割为100个微小段（M=100）。对于每一个微小段，系统构建一个2x2的传输矩阵，该矩阵描述了前向传播模和后向传播模之间的耦合关系。

• 耦合系数计算：根据局部折射率调制深度计算直流耦合系数（$sigma$）和交流耦合系数（$kappa$）。
• 色散方程：对每一段计算传播常数，代码中包含了对实数根（阻带）和虚数根（通带）的完整处理，确保数学运算的稳定性。
• 矩阵级联：通过循环将所有微小段的矩阵相乘，得到表达整体光栅特性的总传输矩阵。
• 反射率提取：假设输出端边界条件，利用总矩阵元素计算输入端的反射系数，进而得到反射率。

2. 变迹函数实现

仿真器内置了多种折射率分布模型的支持逻辑（在子函数中实现）：

• Uniform：沿光栅轴向折射率调制深度保持恒定。
• Gaussian：折射率调制深度沿轴向呈高斯分布，代码中通过参数控制高斯窗的宽度，有效平滑了光栅两端的折射率突变，从而模拟旁瓣抑制效果。
• Raised Cosine：代码逻辑中包含了升余弦变迹的计算公式（虽然主程序主要演示了高斯变迹）。

3. 温度传感机理

在模拟温度影响时，代码并非简单平移光谱，而是更新了底层的物理参数传入TMM算法重新计算：

• 热光效应：根据温度变化量和热光系数更新光纤的有效折射率（$n_{eff}$）。
• 热膨胀效应：根据热膨胀系数更新光栅的物理周期（$Lambda$）。
• 这意味着仿真结果真实的反映了物理参数改变对光谱形状和位置的综合影响。

运行与输出说明

运行主程序后，系统将生成包含四个子图的窗口：

子图 1：均匀FBG反射谱

• 展示基于默认参数（L=10mm, $delta n=2e-4$）计算出的反射光谱。
• 图中标注了峰值反射率的具体数值和计算出的半高全宽（FWHM）。

子图 2：均匀光栅 vs 高斯切趾光栅

• 将均匀光栅（蓝色虚线）与高斯变迹光栅（红色实线）叠加显示。
• 用户可以直观观察到变迹技术如何显著降低中心主峰两侧的旁瓣幅度。

子图 3：温度传感特性模拟

• 绘制在25°C、50°C和75°C三个不同温度点下的反射光谱。
• 随着温度升高，由于折射率增大和光栅周期伸长，布拉格波长向长波方向（红移）移动。

子图 4：布拉格中心波长随温度漂移曲线

• 提取子图3中各个光谱的峰值波长。
• 绘制“温度-波长”关系图。
• 进行线性回归分析，并在图中显示计算出的灵敏度（例如：约 10.9 pm/°C）。

使用方法

1. 下载源代码文件。
2. 在MATLAB环境中打开文件。
3. 直接运行函数，无需输入参数。
4. 观察弹出的图形窗口及命令行窗口的进度提示。