luceda ipkiss教程 65：通过spiral来仿真波导的损耗

硅光设计

已于 2024-03-21 23:48:20 修改

阅读量298

点赞数 6

分类专栏：光子集成线路设计 ipkiss 文章标签： python 开发语言

于 2024-03-21 23:45:14 首次发布

本文链接：https://blog.csdn.net/qq_34316088/article/details/136924231

版权

光子集成线路设计同时被 2 个专栏收录

62 篇文章 12 订阅

订阅专栏

ipkiss

62 篇文章 8 订阅

订阅专栏

本文通过Python代码展示了如何使用SiFab库中的FixedLengthSpiralRounded类来模拟不同波长下波导的损耗，以SWG550和SWG450模板为例，计算并可视化了传输效率和相位变化。

摘要由CSDN通过智能技术生成

案例分享：通过设定固定长度的波导来仿真波导的损耗：
如：
在这里插入图片描述

所有代码如下：

from si_fab import all as pdk
import numpy as np
from matplotlib import pyplot as plt
from picazzo3.wg.spirals import FixedLengthSpiralRounded

trace_templates = [pdk.SWG550(), pdk.SWG450()]
wls_oband = np.linspace(1.26, 1.36, 10001)
wls_cband = np.linspace(1.5, 1.6, 10001)
wls_bands = [wls_cband, wls_oband]

for trace_template, wls_band in zip(trace_templates, wls_bands):
    spiral = FixedLengthSpiralRounded(total_length=10000, n_o_loops=6, trace_template=trace_template)
    spiral.Layout(
        incoupling_length=0.0,
        bend_radius=50.0,
        spacing=4,
        stub_direction="V",  # either H or V
        growth_direction="H",  # either H or V
    )
    spiral.Layout().visualize(annotate=True)
    cm = spiral.CircuitModel()
    S = cm.get_smatrix(wavelengths=wls_band)
    transmission = 10 * np.log10(np.abs(S["out", "in"]) ** 2)
    phase = np.unwrap(np.angle(S["out", "in"]))
    plt.plot(wls_band, transmission, label="in -> out,\nLoss: {:.3g} dB/cm".format(np.mean(transmission)))
    plt.xlabel("Wavelength [um]")
    plt.ylabel("Transmission [dB] ")
    plt.ylim([-0.5, 0])
    plt.title(
        "Trace_template: {}, Spiral length: {} cm".format(spiral.trace_template.name, spiral.total_length / 10000))
    plt.tight_layout()
    plt.legend(fontsize=14, loc=0)
    plt.show()