luceda ipkiss教程 51：画宽度随长度变化的波导

硅光设计

已于 2023-12-22 20:57:18 修改

阅读量457

点赞数 6

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

于 2023-12-22 20:56:13 首次发布

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

版权

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

62 篇文章 12 订阅

订阅专栏

ipkiss

62 篇文章 8 订阅

订阅专栏

本文介绍了一种使用Python和SiFab库在i3.PCell中创建宽度随长度变化的波导结构的方法，通过`width_change_with_length`函数计算边界位置，利用i3.Boundary功能绘制出具有自适应宽度的光学波导。

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

通过i3.Boundary也可以来定义宽度随长度变化的波导，
如：
在这里插入图片描述
全部代码如下：

from si_fab import all as pdk
from ipkiss3 import all as i3
import numpy as np


def width_change_with_length(L, a, m, w0):
    x = np.linspace(0, L, 101)
    y1 = (a * (L - x) ** m + w0) / 2
    y2 = -y1[::-1]
    x_list = np.concatenate((x, x[::-1]))
    y_list = np.concatenate((y1, y2))
    return list(zip(x_list, y_list))


class waveguide(i3.PCell):
    L = i3.PositiveNumberProperty(default=1)
    a = i3.PositiveNumberProperty(default=0.1)
    m = i3.IntProperty(default=4)
    w0 = i3.PositiveNumberProperty(default=0.2)

    class Layout(i3.LayoutView):
        def _generate_elements(self, elems):
            position = width_change_with_length(L=self.L, a=self.a, m=self.m, w0=self.w0)
            elems += i3.Boundary(
                layer=i3.TECH.PPLAYER.SI,
                shape=position)
            return elems

        def _generate_ports(self, ports):
            ports += i3.OpticalPort(name="in",
                                    position=(0, 0),
                                    angle=180.0,
                                    layer=i3.TECH.PPLAYER.SI)
            ports += i3.OpticalPort(name="out",
                                    position=(self.L, 0),
                                    angle=0.0,
                                    layer=i3.TECH.PPLAYER.SI)

            return ports

    class Netlist(i3.NetlistFromLayout):
        pass


if __name__ == '__main__':
    waveguide().Layout().visualize(annotate=True)

通过：

def width_change_with_length(L, a, m, w0):
    x = np.linspace(0, L, 101)
    y1 = (a * (L - x) ** m + w0) / 2
    y2 = -y1[::-1]
    x_list = np.concatenate((x, x[::-1]))
    y_list = np.concatenate((y1, y2))
    return list(zip(x_list, y_list))

可以得到波导边界的坐标点，再通过i3.Boundary可以画出波导