发明名称 一种具有任意光功率分布的光纤光栅设计方法
摘要 本发明公开了一种具有任意光功率分布的光纤光栅设计方法,重构出光纤光栅的折射调制函数。从给定的沿光纤光栅长度的光功率分布出发,对光纤光栅进行分段处理,结合光纤光栅的边界条件(uN=1和v0=0)和能量守恒定律,推导得到谐振波长处光纤光栅内部的前向传输模、后向传输模的幅度。同时,根据传输矩阵法,针对均匀型、切趾型光纤光栅推导得到谐振波长处前向传输模、后向传输模的恒定相位。立足于已得到的前向传输模、后向传输模的幅度和相位,运用耦合模矩阵计算得到反射系数、耦合系数,并从中推算出光纤光栅的折射率调制函数,完成光纤光栅的重构设计。除适用于光纤光栅,本发明还适用于具有任意光功率分布的均匀型、切趾型波导光栅的设计。
申请公布号 CN103076650A 申请公布日期 2013.05.01
申请号 CN201310013694.8 申请日期 2013.01.15
申请人 西南交通大学 发明人 邹喜华;戈伟伟;周丰洲;周志;潘炜;罗斌;闫连山
分类号 G02B6/02(2006.01)I 主分类号 G02B6/02(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 一种具有任意光功率分布的光纤光栅设计方法,该方法包括如下步骤:a)首先将光纤光栅划分为N段,光栅段的编号为1,2,3…k‑1,k,…N,在每一段光栅中,折射率调制视为常量;归一化长度为0~1,在每一段的界面处,都存在前向传输模uk和后向传输模vk,则光纤光栅各分段点的光功率Pk、差分形式的光功率ΔPk为:Pk∝|uk+vk|2          (1)ΔPk=Pk‑Pk‑1         (2)b)定义一个归一化函数: <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>根据式(1)‑(3)可得出: <mrow> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>N</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>ΔPk=(PN‑P0)f(k)        (5)c)从任意给定的光功率分布Pk或者差分形式的光功率分布ΔPk出发,结合光纤光栅的边界条件uN=1和v0=0以及能量守恒定律,计算得到光纤光栅内部前向传输模uk和后向传输模vk:uN=1,v0=0            (6) <mrow> <msup> <mrow> <mo>|</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>v</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>&Proportional;</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <mo>=</mo> <mfrac> <mrow> <mi>&Delta;</mi> <msub> <mi>P</mi> <mi>k</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>arg(uk)=0,arg(vk)=π/2    (8)其中,arg(·)表示取相位;d)立足于已得到的前向、后向传输模,运用传输矩阵: <mrow> <msub> <mi>M</mi> <mi>k</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>&times;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mo>-</mo> <msubsup> <mi>&rho;</mi> <mi>k</mi> <mo>*</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&times;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j&delta;&Delta;</mi> <msub> <mi>z</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j&delta;&Delta;</mi> <msub> <mi>z</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msub> <mi>M</mi> <mi>k</mi> </msub> <mo>&times;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>其中δ为失谐量,Δzk为第k个光栅段的长度,在谐振波长处失谐量δ为0,得到每一光栅段的反射系数ρk: <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>&times;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mo>-</mo> <msubsup> <mi>&rho;</mi> <mi>k</mi> <mo>*</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&times;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Re</mi> <mrow> <mo>(</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>E</mi> <mo>-</mo> <mi>C</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>G</mi> <mo>+</mo> <mi>A</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>F</mi> <mo>-</mo> <mi>D</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>H</mi> <mo>+</mo> <mi>B</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mrow> <mo>(</mo> <msup> <mi>G</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>A</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Im</mi> <mrow> <mo>(</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>E</mi> <mo>-</mo> <mi>C</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>H</mi> <mo>-</mo> <mi>B</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>F</mi> <mo>-</mo> <mi>D</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>G</mi> <mo>-</mo> <mi>A</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msup> <mi>G</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>A</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> </mrow>其中Re(·)和Im(·)分别表示实部、虚部;A=Re(ukuk‑1),B=Im(ukuk‑1),C=Re(ukvk‑1),D=Im(ukvk‑1),E=Re(vkuk‑1),F=Im(vkuk‑1),G=Re(vkvk‑1),H=Im(vkvk‑1);e)由于反射系数ρk与耦合系数qk存在着如下对应关系: <mrow> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>=</mo> <mi>tanh</mi> <mrow> <mo>(</mo> <mo>|</mo> <msub> <mi>q</mi> <mi>k</mi> </msub> <mo>|</mo> <mi>&Delta;</mi> <msub> <mi>z</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> <mfrac> <msubsup> <mi>q</mi> <mi>k</mi> <mo>*</mo> </msubsup> <mrow> <mo>|</mo> <msub> <mi>q</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>进而计算出每一段光栅的耦合系数qk <mrow> <mo>|</mo> <msub> <mi>q</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>=</mo> <mo>[</mo> <mi>ln</mi> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mo>|</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mo>|</mo> <msub> <mi>&rho;</mi> <mi>k</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>]</mo> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <mi>&Delta;</mi> <msub> <mi>z</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>arg(qk)=π‑arg(ρk).          (16)f)最后,根据耦合系数qk、折射率调制函数Δnk的对应关系: <mrow> <mi>&delta;n</mi> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>&Delta;</mi> <msub> <mi>n</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </mfrac> <mi>exp</mi> <mo>{</mo> <mi>j</mi> <mo>[</mo> <msub> <mi>K</mi> <mn>0</mn> </msub> <mi>z</mi> <mo>+</mo> <msub> <mi>&theta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <mi>c</mi> <mo>.</mo> <mi>c</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>q</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mi>j</mi> <mfrac> <msub> <mi>K</mi> <mn>0</mn> </msub> <mrow> <mn>2</mn> <msub> <mi>n</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mi>&Delta;</mi> <msub> <mi>n</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>[</mo> <mo>-</mo> <mi>j</mi> <msub> <mi>&theta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>其中K0=2π/Λ,Λ为光栅周期,c.c表示复共轭,θk(z)=0,计算出折射率调制函数Δnk的幅度和相位: <mrow> <mo>|</mo> <msub> <mi>&Delta;n</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>=</mo> <mfrac> <mrow> <msub> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mn>0</mn> </msub> <mi>&Lambda;</mi> </mrow> <mi>&pi;</mi> </mfrac> <mo>|</mo> <msub> <mi>q</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>arg</mi> <mrow> <mo>(</mo> <mi>&Delta;</mi> <msub> <mi>n</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>arg</mi> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <mi>&pi;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>e)具有此折射率调制函数的光纤光栅即为具有任意光功率分布的光纤光栅。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处