Analysis of the Scalar and Vector Random Coupling Models For a Four Coupled-Core Fiber

The files with simulation results for ECOC 20223 submission "Analysis of the Scalar and Vector Random Coupling Models For a Four Coupled-Core Fiber". "4CCF_eigenvectorsPol" file is the Mathematica code which enables to calculate supermodes (eigenvectors of M(w)) and their propagation constants of 4-coupled-core fiber (4CCF). These results are uploaded to the python notebook "4CCF_modelingECOC" in order to plot them to get Fig. 2 in the paper. "TransferMatrix" is the python file with functions used for modeling, simulation and plotting. It is also uploaded in the python notebook "4CCF_modelingECOC", where all the calculations for figures in the paper are presented. ! UPD 25.09.2023: There is an error in the formula of birefringence calculation. It is in the function "CouplingCoefficients" in "TransferMatrix" file. There the variable "birefringence" has to be calculated according to the formula (19) [A. Ankiewicz, A. Snyder, and X.-H. Zheng, "Coupling between parallel optical fiber cores–critical examination", Journal of Lightwave Technology, vol. 4, no. 9,pp. 1317–1323, 1986]: (4*U**2*W*spec.k0(W)*spec.kn(2, W_)/(spec.k1(W)*V**4))*((spec.iv(1, W)/spec.k1(W))-(spec.iv(2, W)/spec.k0(W))) The correct formula gives almost the same result (the difference is 10^-5), but one has to use a correct formula anyway. ! UPD 9.12.2023: I have noticed that in the published version of the code I forgot to change the wavelength range for impulse response calculation. So instead of seeing the nice shape as in the paper you will see resolution limited shape. To solve that just change the range of wavelengths, you can add "wl = [1545e-9, 1548e-9]" in the first cell after "Total power impulse response". P.s. In case of any questions or suggestions you are welcome to write me an email [e-mail]本数据集包含提交至ECOC 2023（原文疑似笔误，原为20223）的论文《四耦合纤芯光纤标量与矢量随机耦合模型分析》的仿真结果文件。 "4CCF_eigenvectorsPol"文件为Mathematica代码，可用于计算四耦合纤芯光纤（Four Coupled-Core Fiber, 4CCF）的超模（即M(w)的特征向量）及其传输常数。上述结果已上传至Python Notebook脚本"4CCF_modelingECOC"，用于绘图以生成论文中的图2。"TransferMatrix"为包含建模、仿真与绘图所用函数的Python文件，同样上传至"4CCF_modelingECOC"脚本中，该脚本包含了论文所有配图的完整计算过程。 ! 2023年9月25日更新："TransferMatrix"文件中的"CouplingCoefficients"函数存在双折射计算公式错误。其中变量"birefringence"需按照公式(19)[A. Ankiewicz、A. Snyder与X.-H. Zheng，《并行光纤纤芯间的耦合：批判性审视》，《光波技术杂志》，第4卷第9期，第1317-1323页，1986年]进行计算： (4*U**2*W*spec.k0(W)*spec.kn(2, W_)/(spec.k1(W)*V**4))*((spec.iv(1, W)/spec.k1(W))-(spec.iv(2, W)/spec.k0(W))) 修正后的公式计算结果与原公式几乎一致（差值仅为10^-5量级），但仍需使用正确公式。 ! 2023年12月9日更新：我发现代码的发布版本中未修改脉冲响应计算的波长范围，因此你将看到受分辨率限制的波形，而非论文中的理想波形。解决方法为修改波长范围：在"Total power impulse response"后的第一个代码单元中添加"wl = [1545e-9, 1548e-9]"即可。 备注：如有任何疑问或建议，欢迎通过邮箱[e-mail]与我联系。