CN102183814A - 基于混合优化算法的光纤布拉格光栅的逆向设计方法 - Google Patents
基于混合优化算法的光纤布拉格光栅的逆向设计方法 Download PDFInfo
- Publication number
- CN102183814A CN102183814A CN201110140310XA CN201110140310A CN102183814A CN 102183814 A CN102183814 A CN 102183814A CN 201110140310X A CN201110140310X A CN 201110140310XA CN 201110140310 A CN201110140310 A CN 201110140310A CN 102183814 A CN102183814 A CN 102183814A
- Authority
- CN
- China
- Prior art keywords
- bragg grating
- fiber bragg
- index modulation
- grating fbg
- optimization algorithm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Optical Communication System (AREA)
- Diffracting Gratings Or Hologram Optical Elements (AREA)
Abstract
本发明涉及一种基于混合优化算法的光纤布拉格光栅的逆向设计方法,包括获得初始折射率调制的范围:设定光纤布拉格光栅的目标反射率,得到相应的耦合系数,通过耦合系数获得直流折射率调制的范围。重构折射率调制:以直流折射率调制的上边界和下边界作为初始条件,用量子粒子群优化算法去减小适应值函数,从而得到优化的直流折射率调制。设计光纤布拉格光栅:根据优化的直流折射率调制、光栅长度、周期参数设计制作光纤布拉格光栅。
Description
技术领域
本发明涉及一种光纤布拉格光栅的设计方法,尤其涉及一种基于混合优化算法的光纤布拉格光栅的逆向设计方法。
背景技术
光纤布拉格光栅(FBG,即,fiber Bragg gratings)具有体积小、波长选择性好、稳定性好、重复性好、极化不敏感、附加损耗小、耦合性好、外形可变, 适合埋入大型结构中等特点,这使得基于光纤布拉格光栅的光器件成为全光通信网中理想的关键器件,如光纤激光器、光纤光学滤波器、光纤传感器、色散补偿器等。
光纤布拉格光栅的设计主要包括正向设计和逆向设计。正向设计即在光栅参数已知的情况下,应用传输矩阵法、直接傅立叶变换法、或龙格-库塔法等来求得光栅反射率、色散等特性。正向设计的方法虽然简单直接,但是工艺实现起来很不方便,因为在设计过程中,必须首先要对光栅参数进行处理(如取样、切趾等),如果参数设置不合适的话,将会导致设计结果错误。
逆向设计即在光栅特性(目标反射谱)已知的情况下,可以直接重构出相应的光栅参数。因此,逆向设计方法引起了人们的关注。目前主要包括三类逆向设计方法:第一类方法是逆散射算法,主要是层剥算法(DLP,即,discrete layer peeling)。层剥算法是目前最简单高效的光栅逆向设计方法,但由它设计的光栅参数具有结构复杂、光栅长度较长,因此很难在工艺上实现;第二类方法是优化算法,主要包括基因算法、模拟退火算法、协方差矩阵自适应进化算法、粒子群算法等,这种方法往往能够实现满足应用要求的光栅参数,但是往往效率很低下;第三类方法是混合算法,即将逆散射算法与优化算法各取其优点。在申请号为[200910023407.5]的发明专利中,基于层剥算法和非线性最小二乘法的混合算法被应用来设计三角反射谱光纤布拉格光栅。在该发明中合成反射谱时所采用的光栅上的离散点较多,光栅较长,不利于工艺上实现,并且没有和目标反射谱进行对比,导致合成谱的精确性没有更明确地给出。
发明内容
本发明解决的技术问题是:涉及一种基于混合优化算法的光纤布拉格光栅的逆向设计方法,克服现有技术中光纤布拉格光栅设计方法运行效率低、设计出的光栅长度较长、结构复杂、反射谱主要为矩形等问题。
本发明的技术方案是:提供一种基于混合优化算法的光纤布拉格光栅的逆向设计方法,包含以下步骤:
获得折射率调制的初始范围:
重构折射率调制:
设计光纤布拉格光栅:
根据优化的直流折射率调制、光栅长度、周期参数设计制作光纤布拉格光栅。
本发明的进一步技术方案是:在设计光纤布拉格光栅步骤中,所述光纤布拉格光栅是具有三角形反射谱的布拉格光纤光栅。
本发明的进一步技术方案是:所述量子粒子群算法是带有变异因子的自适应量子粒子群算法。
本发明的进一步技术方案是:在获得折射率调制的初始范围步骤中,用离散层剥法计算得到相应的耦合系数。
本发明的进一步技术方案是:在获得折射率调制的初始范围步骤中,用离散层剥法计算相应的耦合系数时包括选择光栅长度,在重构折射率调制步骤中, 重构折射率调制时包括选择光栅长度,在重构折射率调制时所选的光栅长度较计算相应耦合系数时选择的光栅长度短。
本发明的进一步技术方案是:在所述基于混合优化算法的光纤布拉格光栅的逆向设计方法中,在设计光纤布拉格光栅后,还包括设计光纤布拉格光栅步骤中通过传输矩阵法验证的反射谱与目标反射谱进行比对。
本发明的技术效果是:本发明一种基于混合优化算法的光纤布拉格光栅的逆向设计方法,包括获得初始折射率调制范围:设定光纤布拉格光栅的目标反射率,得到相应的耦合系数,通过耦合系数获得直流折射率调制的初始范围上边界和下边界。重构折射率调制:以初始直流折射率调制的初始范围作为初始条件,用量子粒子群优化算法去减小适应值函数,从而得到直流折射率调制。设计光纤布拉格光栅:根据得到的直流折射率调制、光栅长度、周期参数设计制作光纤布拉格光栅。本发明一种基于混合优化算法的光纤布拉格光栅的逆向设计方法具体如下特点:
1、本发明所提供的布拉格光纤光栅的逆向设计方法,可以高效地设计出光栅制作所需要的光纤布拉格光栅参数。
2、本发明所提供的布拉格光纤光栅的逆向设计方法,可以根据层剥算法高效地重构出直流折射率调制的初始范围,所述以初始下边界和初始上边界,应用带有变异因子的自适应量子粒子群算法设计光纤布拉格光栅。该方法只有一个控制算法速度的参数需要根据算法的运行过程进行调整,克服了传统光纤布拉格光栅优化设计方法的不足。该方法可以设计具有三角形反射谱的布拉格光纤光栅,设计出的光栅具有长度短、没有啁啾、折射率调制较低易于制作、反射谱线性好等优点。
3、本发明所提供的布拉格光纤光栅的逆向设计方法,可以用来设计大带宽的三角谱光纤布拉格光栅,因此可以在传感系统中用于传感参数较大动态范围的光波长解读器件。
4、本发明所提供的布拉格光纤光栅的逆向设计方法,可以用来设计多信道的三角谱光纤布拉格光栅,因此可以在光传感系统中作为多参量光波长解读器件。
5、本发明所提供的布拉格光纤光栅的逆向设计方法设计出的三角谱光纤布拉格光栅具有结构简单、成本低、体积小、易于制作等特点。
附图说明
图1为本发明的流程图。
图2为本发明的折射率调制及相应的反射光谱图示。
图3为本发明大带宽单信道三角反射谱光纤布拉格光栅的折射率调制及相应的反射光谱图示。
图4为本发明三信道三角反射谱光纤布拉格光栅的折射率调制及相应的反射光谱图示。
具体实施方式
下面结合具体实施例,对本发明技术方案进一步说明。
如图1所示,本发明的具体实施方式如下:
具体过程如下:
复反射系数可以表示为
本发明算法提出一种符合三角反射谱(AS-FBGs,即,triangular-spectrum fiber Bragg gratings)设计要求的适应值函数:
(7)
其中: 是反射谱的采样点总数;是反射谱上第 j个采样点的波长; 和 分别指光栅的目标反射谱和应用发明算法所设计光栅的反射谱;权重参数 描述了反射谱中不同波长反射系数误差的平衡关系。是ASFBG的带宽; 和分别是内最大和最小的采样点的阶数。
带变异因子的自适应量子粒子群(AQPSO,即,adaptive quantum-behaved particle swarm optimization)算法是一种具有量子行为的优化算法,算法的设置简单、控制参数少、搜索能力强。算法中粒子的移动遵循以下表达式:.
(8)
(10)
其中:向量 是粒子群的吸引点; 是第个粒子在第d维空间中的最优解;是粒子群中所有的粒子在第d维空间中的最优解; 被定义为所有粒子最好位置的均值;和是均匀分布在 [0,1]上的随机数; M 是粒子群的粒子数;是控制算法收敛速度的收缩扩张系数,它是算法唯一可调整的参数; 是指算法的迭代代数。
(11)
(13)
带有变异因子的 AQPSO 算法可以归纳如下:
10.重复步骤3–9 直到获得设定的最好适应值。
步骤300:设计光纤布拉格光栅,即,根据得到的直流折射率调制、光栅长度、周期参数设计制作光纤布拉格光栅。
根据得到的直流折射率调制、光栅长度、周期等参数设计制作光纤布拉格光栅。具体实施例中,设计的窄带宽单信道三角谱光纤布拉格光栅的折射率调制如 图2(a)所示,图2(a)为设计的窄带宽单信道三角反射谱光纤布拉格光栅的折射率调制,其中,横轴为光栅长度,纵轴为折射率调制;图2(b)为图2(a)相应的光纤布拉格光栅的三角反射谱,其中带宽为0.2nm,横轴为光栅长度,纵轴为反射率;图2(c)为算法的适应值函数在优化过程中的收敛曲线,其中横轴为迭代次数,纵轴为适应值函数的变化曲线。设计的大带宽单信道三角谱光纤布拉格光栅的折射率调制如图3(a)所示,图3(b)为图3(a)相应的光纤布拉格光栅的三角反射谱,带宽为2.5 nm,横轴为光栅长度,纵轴为反射率。设计的三信道三角谱光纤布拉格光栅的折射率调制如图4(a)所示,图4(b)为图4(a)相应的光纤布拉格光栅的三角反射谱:横轴为光栅长度,纵轴为反射率。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (9)
3.根据权利要求2所述的基于混合优化算法的光纤布拉格光栅的逆向设计方法,其特征在于,在设计光纤布拉格光栅步骤中,所述光纤布拉格光栅是具有三角形反射谱的布拉格光纤光栅。
6.根据权利要求1所述基于混合优化算法的光纤布拉格光栅的逆向设计方法,其特征在于:所述量子粒子群算法是带有变异因子的自适应量子粒子群算法。
7.根据权利要求1所述的基于混合优化算法的光纤布拉格光栅的逆向设计方法,其特征在于:在获得折射率调制的初始范围步骤中,用离散层剥法计算得到相应的耦合系数。
8.根据权利要求1所述的基于混合优化算法的光纤布拉格光栅的逆向设计方法,其特征在于:在获得折射率调制的初始范围步骤中,用离散层剥法计算相应的耦合系数时包括选择光栅长度,在重构折射率调制步骤中, 重构折射率调制时包括选择光栅长度,在重构折射率调制时所选的光栅长度较计算相应耦合系数时选择的光栅长度短。
9.根据权利要求1所述的基于混合优化算法的光纤布拉格光栅的逆向设计方法,其特征在于:在所述基于混合优化算法的光纤布拉格光栅的逆向设计方法中,在设计光纤布拉格光栅后,还包括设计光纤布拉格光栅步骤中通过传输矩阵法验证的反射谱与目标反射谱进行比对。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110140310XA CN102183814B (zh) | 2011-05-27 | 2011-05-27 | 基于混合优化算法的光纤布拉格光栅的逆向设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110140310XA CN102183814B (zh) | 2011-05-27 | 2011-05-27 | 基于混合优化算法的光纤布拉格光栅的逆向设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102183814A true CN102183814A (zh) | 2011-09-14 |
CN102183814B CN102183814B (zh) | 2012-12-19 |
Family
ID=44570010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110140310XA Expired - Fee Related CN102183814B (zh) | 2011-05-27 | 2011-05-27 | 基于混合优化算法的光纤布拉格光栅的逆向设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102183814B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108227188A (zh) * | 2018-01-30 | 2018-06-29 | 山东大学 | 一种基于混沌粒子群算法的全随机硅基波导光栅设计方法 |
CN108254000A (zh) * | 2017-12-08 | 2018-07-06 | 上海电力学院 | 一种基于改进模拟退火法的fbg传感网光谱解调方法 |
CN110596894A (zh) * | 2019-10-25 | 2019-12-20 | 清华大学深圳国际研究生院 | 一种用于衍射光学元件设计的方法和系统 |
CN112033645A (zh) * | 2020-08-07 | 2020-12-04 | 深圳市速腾聚创科技有限公司 | 波导光栅结构参数确定方法、装置、存储介质及激光雷达 |
WO2023240772A1 (zh) * | 2022-06-16 | 2023-12-21 | 苏州大学 | 基于稀疏度计算的硅基光学微环滤波器逆向设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003510627A (ja) * | 1999-09-21 | 2003-03-18 | ザ・ユニバーシティ・オブ・シドニー | グレーティング設計 |
CN101963682A (zh) * | 2009-07-23 | 2011-02-02 | 中国科学院西安光学精密机械研究所 | 三角形反射谱光纤布拉格光栅的合成方法 |
-
2011
- 2011-05-27 CN CN201110140310XA patent/CN102183814B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003510627A (ja) * | 1999-09-21 | 2003-03-18 | ザ・ユニバーシティ・オブ・シドニー | グレーティング設計 |
CN101963682A (zh) * | 2009-07-23 | 2011-02-02 | 中国科学院西安光学精密机械研究所 | 三角形反射谱光纤布拉格光栅的合成方法 |
Non-Patent Citations (2)
Title |
---|
《光电工程》 20061231 余华清,姜向东 基于反向耦合模原理光栅的薄膜模型研究 全文 1-9 第33卷, 第12期 * |
《应用光学》 20090731 谢军华等 用粒子群优化算法设计光纤布拉格光栅 全文 1-9 第30卷, 第4期 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108254000A (zh) * | 2017-12-08 | 2018-07-06 | 上海电力学院 | 一种基于改进模拟退火法的fbg传感网光谱解调方法 |
CN108227188A (zh) * | 2018-01-30 | 2018-06-29 | 山东大学 | 一种基于混沌粒子群算法的全随机硅基波导光栅设计方法 |
CN110596894A (zh) * | 2019-10-25 | 2019-12-20 | 清华大学深圳国际研究生院 | 一种用于衍射光学元件设计的方法和系统 |
CN112033645A (zh) * | 2020-08-07 | 2020-12-04 | 深圳市速腾聚创科技有限公司 | 波导光栅结构参数确定方法、装置、存储介质及激光雷达 |
WO2023240772A1 (zh) * | 2022-06-16 | 2023-12-21 | 苏州大学 | 基于稀疏度计算的硅基光学微环滤波器逆向设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102183814B (zh) | 2012-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102183814B (zh) | 基于混合优化算法的光纤布拉格光栅的逆向设计方法 | |
GB2334594A (en) | Arrayed waveguide grating device | |
Turduev et al. | Ultracompact photonic structure design for strong light confinement and coupling into nanowaveguide | |
Khan et al. | Determination of the best apodization function and grating length of linearly chirped fiber Bragg grating for dispersion compensation. | |
Chen et al. | Group velocity in lossy periodic structured media | |
CN103424995B (zh) | 导模共振滤光片光刻胶层的优化方法 | |
KR20050013165A (ko) | 다수의 광파장 검출을 위한 방법 및 장치 | |
Guo et al. | Broad band beam splitter based on the double-groove fused silica grating | |
Faiyaz et al. | Optimization of apodization profile of chirped fiber Bragg grating for chromatic dispersion compensation: dispersion compensation using chirped apodized FBG | |
CN113673677B (zh) | 一种非线性激活函数relu的实现方法、设备及介质 | |
WO2023240772A1 (zh) | 基于稀疏度计算的硅基光学微环滤波器逆向设计方法 | |
CN113392965B (zh) | 一种哈达玛积的实现方法、设备及存储介质 | |
CN110187499B (zh) | 一种基于神经网络的片上集成光功率衰减器的设计方法 | |
Balaji et al. | Design and analysis of waveguide taper for photonic crystal demultiplexer | |
Raghuwanshi et al. | Dispersion and peak reflectivity analysis in a non-uniform FBG based sensors due to arbitrary refractive index profile | |
Sakurai et al. | Proposal of tunable hollow waveguide distributed Bragg reflectors | |
Wang et al. | Reverse Design of Multifunctional Demultiplexing Devices | |
Xiao et al. | Mode control in planar waveguide grating couplers with double surface corrugation | |
Fasihi et al. | Crosstalk investigation in channel-drop filters with coupled-cavity based wavelength-selective reflection feedbacks | |
CN103792755A (zh) | 一种啁啾结构低阈值双稳态的全光开关设计方法 | |
Gong et al. | The steady SRS analysis theory of DWDM transmission system in single-mode silica fiber | |
CN101963682A (zh) | 三角形反射谱光纤布拉格光栅的合成方法 | |
CN218122293U (zh) | 一种平顶波分复用器 | |
CN111999888B (zh) | 一种基于光致聚合物浓度扩散模型优化材料折射率调制度和衍射效率的方法 | |
Khan et al. | Performance evaluation of different apodization profiles of linearly chirped FBG for dispersion compensation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121219 Termination date: 20140527 |