简介面向光通信学习与科研人员的2×2光纤耦合器MATLAB仿真资源聚焦于利用数值计算理解光功率分配与合路原理适合正在学习光纤通信、波导光学或从事耦合器件设计的读者。压缩包共3个文件包含1个MATLAB源程序及2张仿真结果图整体仅97KB小巧易用脚本可用于定义光纤参数、计算耦合常数、建立耦合器模型并绘制输出功率分布图片则直观呈现不同Δβ条件下耦合长度与输出特性的变化。目前已有185人浏览学习对于希望快速上手光纤耦合器仿真、直观体会参数影响并在此基础上开展优化设计的读者是一份简明实用的参考资料。1. 2×2耦合器的分光比只是起点从ouheqi.rar说起拿到ouheqi.rar_2×2耦合器_coupler_fiber_光纤_光纤耦合这种命名风格的压缩包时我一般不会期待里面装着什么新原理大概率是某个工程目录的打包产物工艺记录、测试脚本、拉锥机的参数备份或者几张拍得不算清楚的光纤端面照片。真正有价值的是包里体现的工艺流程——2×2耦合器的核心问题从来不是“能不能分光”而是分光比怎么控、插损怎么压、偏振相关损耗怎么降。3dB耦合器要控在50±1%WDM型的要拉出特定波长响应这些靠同一台拉锥机、同一套监控逻辑差别只在停火时机和封装手法。这篇就按我自己做器件验证时会走的完整路线讲耦合模理论怎么指导参数设计、熔融拉锥怎么控制停拉点、成品怎么测均匀性和波长平坦度、分光比偏了从哪里排查最后落到一种不用高端设备也能判断耦合器封装质量的验证技巧。2. 从拍频周期到熔融拉锥2×2光纤耦合器的耦合模型论与制备参数2.1 两纤芯距离如何决定2×2耦合器的功率转移周期2×2耦合器的工作基础是两根平行放置的光纤纤芯之间的倏逝波耦合。光从一根纤芯注入后能量并不会停留在原纤芯内而是以拍频的形式在两根纤芯之间周期性地来回转移。耦合模方程可以写成dA1/dz -jκ·A2, dA2/dz -jκ·A1其中 A1、A2 是两根纤芯内的模场复振幅κ 是耦合系数z 是传播距离。求解后直通端和耦合端的功率分别是P1(z) cos²(κz), P2(z) sin²(κz)这里的 κz 就是拍频角。κ 的单位是 m⁻¹它由纤芯间距 d、光纤数值孔径和工作波长共同决定。两根纤芯拉得越近κ 越大能量交换得越快。实际拉锥过程中光纤被加热软化并拉伸纤芯不断变细、间距不断缩小κ 持续增大功率转移曲线也就从慢变快。熔融拉锥的“拉”不是一个匀速等效过程越往后 κ 增速越快这直接决定了拉锥机的停止位置必须精确控制——多拉 0.2 毫米分光比可能就从 50% 跳到 60% 了。耦合周期 Lc 定义为完成一个完整功率转移周期所需的长度Lc π/(2κ)。一个标准 3dB 耦合器实际上只用了这个周期的四分之一所以它的工作区对 κ 的变化特别敏感任何封装应力导致的微弯都会轻微改变 κ进而扰动分光比。理解了这一点后面的排错逻辑就清楚了分光比偏了先想是不是几何尺寸变了。2.2 熔融拉锥coupler的最小制备流程与停拉点控制常见的 2×2 熔融拉锥制备流程不算复杂但每个环节的容错空间都很小。我习惯的顺序是剥涂覆层、清洁纤芯、光纤对轴、固定、加热拉锥、监控停拉、封装固化。两根单模光纤各自剥掉约 30mm 涂覆层用无尘纸蘸无水乙醇擦拭裸纤区这一步直接决定拉锥后的强度。然后把两根裸纤平行放入拉锥机的夹具保证两纤芯在显微镜下对准间距误差过大时耦合效率会明显下降但对准本身需要经验我一般靠侧视图像判断两纤边缘是否贴合。拉锥的热源有两种常见选择丁烷焰和电阻丝加热。丁烷焰升温快、热区宽适合拉宽带耦合器电阻丝热区更窄、重复性更好适合做 WDM 型窄带器件。火焰宽度直接影响拉锥区的过渡曲线火焰太宽时锥区被拉长耦合器对波长更敏感拉出的宽谱平坦度反而变差。拉锥速度我一般设在 0.3~0.5 mm/s太慢了玻璃在热区停留过久纤芯过度扩散太快了过渡区过陡机械强度扛不住。拉锥机上监控停拉点的经典做法是从 Port1 注入工作波长的光Port3 和 Port4 各接一支光电探测器实时采集两路光功率 P_pass 和 P_coup。耦合比 C 定义为 P_coup / (P_pass P_coup)。当 C 到达目标值附近时就立即停火但考虑停火后残余热量的继续作用我一般会在比目标值低 1~2 个百分点时先关闭热源让余热把最后一段走完。2.3 拉锥停止判据为什么必须用双通道实时监控目标耦合比 (C)直通端功率占比耦合端功率占比经验停拉余量30%70%30%低于目标 1.0%50%50%50%低于目标 1.5%60%40%60%低于目标 2.0%只监控单通道存在一个隐患光源功率不稳定时单通道功率波动会被误判为分光比变化。双通道同时采集后用归一化比值替代绝对功率光源的慢漂移就被抵消掉了。停拉时机的判断依赖的是比值而非绝对功率这是器件制备和普通功率测量之间一个容易被忽略的差别。为了验证不同拉锥长度下的分光比走向可以用一个简单的 Python 脚本模拟拍频曲线import numpy as np def coupling_ratio(z, kappa): # z: 耦合区长度(mm), kappa: 耦合系数(1/mm) pass_power np.cos(kappa * z) ** 2 coup_power np.sin(kappa * z) ** 2 c coup_power / (pass_power coup_power) return pass_power, coup_power, c # 模拟耦合系数逐渐增大时同一长度下的分光比变化 kappa_values [0.5, 0.8, 1.0, 1.2] z_fixed 1.0 # 固定耦合区长度1mm for k in kappa_values: p, q, c coupling_ratio(z_fixed, k) print(fkappa{k}/mm: 直通{p:.2%}, 耦合{q:.2%}, 耦合比C{c:.2%})这段代码演示的核心是固定耦合区长度 zκ 从 0.5 增到 1.2 时耦合比 C 会快速上升。实际拉锥中 κ 的变化正是由纤芯间距缩小和纤芯直径变细共同驱动的所以代码里 κ 的增大等价于拉锥程度的加深。拉锥机上监控软件做的本质上是反过程——看着 C 变大把它到目标值时对应的电机位置换算成拉伸位移写进工艺卡作为下一次拉锥的预停位置参考。第一次做某个新规格时我会先在废纤上跑一遍裸拉记下 C10%, 30%, 50%, 70% 对应的电机坐标再在正式器件上分段逼近避免一次拉过头。3. 测试熔融拉锥coupler的插损、均匀性与波长平坦度3.1 用双波长测试台标定2×2耦合器的分光比与插损器件做出来后第一件事是测试。分光比在不同波长下并不一致对宽带耦合器来说1310nm 和 1550nm 的耦合比可能有几个百分点的差异对 WDM 型耦合器两个波长的分光比可能是 10/90 和 90/10 的差别。所以标定一台耦合器的第一步就是双波长测试。测试台结构不复杂可调激光源、光功率计、FC/PC 跳线若干、一支 1×2 光开关如果需要自动化。先把光源直接接功率计记录参考光功率 Pref再接上器件分别测 Port3 和 Port4 的输出功率。插损包括两部分附加损耗Excess LossEL和插入损耗Insertion LossIL。EL 定义为输入功率与两输出端总功率之差单位 dBIL 则是对某一输出端单独计算的。一只合格的单模 2×2 耦合器附加损耗一般在 0.05~0.2 dB 之间插入损耗则取决于分光比。50/50 耦合器的典型 IL 在 3.2~3.5 dB 附近均匀性两端口 IL 之差通常要求小于 0.5 dB。测试项测试条件合格判据插入损耗 IL1310nm/1550nmLD 稳定 30 分钟≤ 3.6 dB50/50 型附加损耗 EL双端口功率之和 vs 输入功率≤ 0.2 dB均匀性两端 IL 差值≤ 0.5 dB回波损耗输入端口反射≥ 50 dBPDL偏振态扫描插损极差≤ 0.1 dB3.2 参考光功率的取法与连接器带来的测量误差功率计的参考值最好在器件接入前现场采得而不是直接信任激光器的设置功率。原因是连接器的插入损耗会随着插拔次数和端面清洁程度变化参考值里包含了这些固定损耗后后续计算的插损才是可比的。采参考值时光源输出端接一支高质量的 FC/APC 跳线测得的功率作为 Pref。然后跳线接器件 Port1Port3 和 Port4 分别接功率计。如果没有双通道功率计就用 1×2 光开关轮流切换但开关的重复性误差通常 ±0.05 dB会混入测量所以能用双通道就尽量用双通道。另一个常见的坑是端口接反。2×2 耦合器从 Port1 注入和从 Port3 注入得到的分光比可能不同尤其拉锥区不对称时更明显。验收单上必须标注测试方向或者在测试脚本里对两个输入方向都测一遍取均匀性较差的作为最终结果。这个方向性差异本身就是拉锥区对称性的一个判据差异超过 0.3 dB 时基本可以断定拉锥时两纤加热不均匀。3.3 自动采集分光比数据的Python测试脚本批量化测试时手工记录功率计读数既慢又容易出错。常见做法是用 Python 通过 GPIB/USB 读取功率计数据。下面是一个简化的自动采集脚本骨架import pyvisa import time import csv rm pyvisa.ResourceManager() # 替换为实际功率计的VISA地址 power_meter rm.open_resource(GPIB0::14::INSTR) laser rm.open_resource(GPIB0::10::INSTR) def set_wavelength(wl_nm): laser.write(f:SOUR:WAV {wl_nm}NM) time.sleep(2) # 等待波长稳定 def read_power(channel1): # 读取指定通道功率单位dBm return float(power_meter.query(f:MEAS:POW? CH{channel})) results [] for wl in [1310, 1550]: set_wavelength(wl) set_power float(laser.query(:SOUR:POW?)) p_pass read_power(1) p_coup read_power(2) el set_power - 10 * np.log10(10**(p_pass/10) 10**(p_coup/10)) il_pass set_power - p_pass il_coup set_power - p_coup uniformity abs(il_pass - il_coup) results.append([wl, il_pass, il_coup, el, uniformity]) with open(coupler_test.csv, w, newline) as f: writer csv.writer(f) writer.writerow([wavelength_nm, IL_pass_dB, IL_coup_dB, EL_dB, uniformity_dB]) writer.writerows(results)脚本的关键在于每切换一次波长都留足稳定时间这直接决定分光比测量的重复性。DFB 激光器的波长切换后输出功率可能需要几百毫秒到几秒才稳定先查功率再切换下一个波长是更稳妥的顺序。EL 用 dBm 值先转线性再相减不能直接在 dBm 域相加这是新手最容易算错的地方。4. 分光比不准时的排错拉锥参数、封装应力与波长漂移4.1 拉锥速度与火焰宽度如何影响耦合比重复性分光比偏了习惯性动作是调拉锥长度但拉锥速度和火焰宽度可能才是真正的元凶。拉锥速度过快时玻璃在热区内来不及充分软化纤芯扩散不足κ 增速低于预期速度过慢时玻璃过度软化纤芯塌陷κ 增速又高于预期。同一台拉锥机速度从 0.3 mm/s 调到 0.6 mm/s同一目标位移对应的分光比可能差出 8~10 个百分点。判断该调哪个参数的方法是看整条拉锥曲线如果分光比在 10% 以内还挺准越往后偏得越多说明拉锥后期 κ 增速模型没对准优先调停拉余量如果从一开始就跟不上预期说明初始对准或火焰温度有问题。火焰宽度的影响主要体现在过渡区斜率和波长敏感性上火焰偏窄时拉出的耦合器在 1310nm 和 1550nm 两个波长的分光比差异会变大。4.2 封装应力导致的偏振相关损耗过高耦合器拉好之后要封装进石英管或钢套用胶固定。胶水固化时的收缩会对裸纤区施加不均匀应力导致双折射表现为 PDL 升高。很多耦合器测试台上测 PDL 没问题封装后再测就超了原因就出在胶水。规避方法有两个方向一是控制胶量让胶水只覆盖在涂覆层区域裸纤区的应力由石英管本身承载二是封装后做一次热循环-40℃ 到 85℃再测 PDL让应力充分释放和重分布把不合格品在这一步筛出去。另一个应力来源是尾纤的盘放方式。测试时尾纤大幅弯曲或扭转会改变局部双折射间接影响耦合区。测 PDL 时尾纤要保持自然盘绕状态不能用手拽着也不能压着重物。测完一次后换个盘绕方向再测一次两次 PDL 的差异如果超过 0.05 dB说明器件对尾纤应力敏感这种货不适合用在有振动或温度循环的环境中。4.3 窄带器件与宽带器件的区分波长平坦度从哪来异常现象可能原因优先检查项分光比整体偏高拉锥过度电机位移是否超限分光比整体偏低拉锥不足 / 火焰温度低停拉时机的耦合比值1310nm正常1550nm偏差大过渡区过陡 / 火焰过窄火焰宽度调整后重拉PDL 封装后升高胶水固化应力胶量控制或换胶均匀性差两纤对准偏差显微镜下重新对轴宽带耦合器需要的是过渡区内纤芯直径变化尽量缓慢这样不同波长看到的几何结构差异小。拉锥时如果火焰宽度不够过渡区短而陡等效波导变得更依赖于波长分光比随波长变化的斜率就会变大。做宽带器件时我一般会把火焰宽度调到 4~5mm拉锥速度降到 0.2 mm/s 左右代价是拉锥时间变长、总长变长额外的附加损耗会增加 0.1 dB 左右。用偏振控制器扫描输入偏振态记录插损的最大最小值之差就能快速估算 PDL。这个测试不需要专用仪器一支偏振控制器加一支功率计就能做。5. 无源光网络与双工收发2×2耦合器在光通信链路里的3个典型接法5.1 PON OLT侧分光与监控链路在 PON 系统里2×2 耦合器最常见的用途是在 OLT 侧做分光监控。OLT 的发射端和接收端共用一根光纤时需要分光器或 WDM 器件来分离上下行信号。1% 的分光比耦合器常用来抽取少量下行光做功率监测主链路保留 99% 光功率抽出的 1% 进光电探测器判断光模块输出是否正常。这种场景对分光比精度要求不高但对插入损耗的长期稳定性要求很高因为监控链路的功率变化会被当作告警信号触发。分光比 1% 的耦合器在 1550nm 下插损较大主链路损耗约 0.04~0.1 dB加上连接损耗后仍能满足 PON 的功率预算。关键是耦合器的回波损耗要足够高回光进入 OLT 激光器会造成光谱抖动所以这类应用普遍要求回损高于 50 dB。选型时优先挑 1×2 单模耦合器里标称回损 55 dB 的型号测试记录里回损低于这个值的直接退给厂家换货。5.2 双纤双向收发器里的合波与分波双纤双向BIDI光模块内部通常集成 WDM 型 2×2 耦合器把 1310nm 的上行光和 1490nm 或 1550nm 的下行光合到一根光纤上。这种 WDM 耦合器与普通 3dB 耦合器的差别在于拉锥长度和停拉点WDM 型要在拍频周期的一半附近停拉让两个波长对应的输出端口刚好处于互补状态。一个端口对 1310nm 呈现接近全通95% 以上通过、对 1550nm 呈现接近交叉95% 以上转移到另一端口的状态。这类器件的验收测试必须用双波长同时注入而不是单波长切换。单波长测到的是纯串扰双波长同时注入才会暴露非线性串扰和模式耦合的问题。测试时先在 1310nm 下测 Port3 和 Port4再在 1550nm 下测同一对端口两波长的隔离度隔离度 非目标端口的插损 - 目标端口的插损一般要求大于 15 dB。隔离度不够时上下行信号会互相干扰接收机灵敏度下降几个 dBm 是常有的事。5.3 实验室里搭一套分光比验证链路没有专用测试台时用实验室常见仪器也能验证一只 2×2 耦合器是否堪用。需要的东西一台可调激光源或者一台带 1310/1550 双波长的光源、两支光功率计、两支 FC/APC 跳线。把光源接到 Port1一支功率计接 Port3另一支接 Port4先测参考光功率再算分光比分别给两个波长做一遍。下面是快速计算分光比和链路预算的 Python 脚本def link_budget(p_tx_dbm, p_rx_dbm, split_ratio): # 分光比: 耦合端占比 0~1 ldB -10 * math.log10(split_ratio) # 耦合端理论损耗 il p_tx_dbm - p_rx_dbm - ldB # 实际插损 return il p_tx 2.0 # 光源输出 2dBm p_pass -21.3 # 实测直通端功率 dBm p_coup -18.5 # 实测耦合端功率 dBm total 10 ** ((p_pass - p_tx) / 10) 10 ** ((p_coup - p_tx) / 10) el 10 * math.log10(total) c 10 ** ((p_coup - p_tx) / 10) / total print(f附加损耗 EL{el:.2f}dB, 耦合比 C{c:.1%})注意这里的 EL 是功率总和与输入功率之比单位是 dB正常情况是负值习惯上取绝对值报告为“附加损耗 xx dB”。分光比则用耦合端与总功率的线性比值来表示两个公式必须分开不能用同一个 dB 值直接相减。6. 用回损值和偏振稳定性验证一只2×2耦合器是否合格常规测试做完了还有两个容易被忽略的检查项用来暴露封装质量和长期可靠性的隐患。第一个是回损Return LossRL。回损反映的是输入端口自身的反射光水平和插损之间没有直接关系。一只插损合格的耦合器回损可能只有 30dB这会严重干扰上游光模块。没有回损仪时可以用一支 3dB 耦合器和功率计搭简易回损测试光源经 3dB 耦合器接到待测器件 Port1功率计接在 3dB 耦合器的另一臂上测到的反射功率减去已知光路损耗就是 RL。这个方法精度不高但判断“有没有异常反射尖峰”足够了。第二个检查项是偏振稳定性而且要用“摸了才能知道”的方法测完插损后用手轻轻弯折靠近封装区的尾纤观察两路输出功率是否有明显波动。弯折会改变局部弯曲半径进而改变模场分布和双折射。如果弯折幅度不大但功率波动超过 0.2 dB说明耦合区对包层模泄漏和应力过于敏感封装质量堪忧。这个方法虽然粗暴但能在没有偏振控制器的情况下快速筛掉封装不牢固的器件。验证一只 2×2 耦合器是否合格我建议的顺序是单波长插损初测 → 双波长分光比 → 尾纤弯折观察 → 回损测量 → 最后再过一遍温度循环。温度循环可以省略弯折观察不能省。等温度循环做完再回来复测插损如果复测值和初测值差超过 0.3 dB那只耦合器大概率会在现场使用几个月后逐渐劣化——这种问题在出厂全检里基本测不出来只有双波长复测才能定位。没有温度循环设备时至少把器件放在热风枪出口 50cm 处吹一分钟等它冷却后再测一次插损能多筛掉一部分隐患。本文还有配套的精品资源点击获取