2024 年华数杯全国大学生数学建模竞赛题目B 题 VLSI 电路单元的自动布局完整思路 源代码 模型 结果分享

超大规模集成电路(VLSI，Very Large Scale Integration)将大量电路单元集成于单一芯片。随着设计复杂度增加，如今开展VLSI设计已离不开电子设计自动化(EDA，Electronic DesignAutomation)工具的支持。EDA 作为算法密集型产业，需要对数千种情境进行快速设计探索，是国家关键技术领域。其中，电路单元的自动布局是 EDA 研究的核心问题之一

电路单元的自动布局旨在矩形布局区域内确定所有电路单元位置,以最小化单元之间总连接线长并避免单元重叠。由于这是一个NP-难问题，通常分为全局布局和详细布局两个步骤。全局布局大致确定单元位置，允许单元重叠;详细布局则消除重叠并进一步优化。本问题聚焦于全局布局，将电路单元视为不同大小的矩形，矩形内分散有若干个连线接口，电路单元之间通过连线接口形成若干组连接关系。全局布局的目标是最小化总连接线长，同时满足单元密度约束。总连接线长等于每组有连接关系的电路单元的线长之和。由于布局阶段尚未实际布线每组线长通常可通过半周长线长(HPWL，Half-Perimeter Wirelength)或直线型斯坦纳最小树(RSMT，Rectilinear Steiner Minimal Tree)估计，要求连线水平或竖直。HPWL为连线接口外接矩形周长的一半，RSMT为通过插入斯坦纳点构建的线段长度之和。单元密度约束通过将矩形布局区域网格化后计算。每个网格的单元密度等于与网格重叠的电路单元面积和网格面积的比值,限制不超过特定阈值。附件1提供全局布局的中间状态，包括每组有连接关系的电路单元及其连线接口名称、连线接口坐标和对应的 HPWL和 RSMT 线长。附件2给出布局区域尺寸、网格划分粒度和密度阈值、电路单元的尺寸、坐标及其连线接口的基本信息。

请建立数学模型解决以下问题:

问题1 ：图2展示了3组具有不同连线接口数的 HPWL和RSMT 线长估计示意图。RSMT是布局阶段理想的线长表征，但是构建斯坦纳树是NP难问题。HPWL,简单有效，但对多连线接口情形估计偏小。根据附件1提供的信息，请设计一个与电路单元连线接口坐标相关的线长评估模型。该模型应满足:(1)每组估计线长与对应 RSMT 的差值尽可能小:(2)能应用于评估附件1中的总连接线长。

1. 问题背景
   超大规模集成电路(VLSI)设计涉及将大量电路单元集成到单一芯片中。随着设计复杂度的增加，电子设计自动化(EDA)工具在VLSI设计中的作用愈发重要。EDA工具需要快速探索数千种设计情景，而电路单元的自动布局是EDA研究的核心问题之一。电路单元的自动布局旨在矩形布局区域内确定所有电路单元的位置，以最小化单元之间的总连接线长度并避免单元重叠。由于这是一个NP难问题，通常分为全局布局和详细布局两个步骤。
   在全局布局阶段，电路单元视为不同大小的矩形，每个矩形内有若干个连线接口。电路单元之间通过连线接口形成若干组连接关系。全局布局的目标是最小化总连接线长度，同时满足单元密度约束。总连接线长度通常通过半周长线长(HPWL, Half-Perimeter Wirelength)或直线型斯坦纳最小树(RSMT, Rectilinear Steiner Minimal Tree)估计。HPWL为连线接口外接矩形周长的一半，RSMT为通过插入斯坦纳点构建的线段长度之和。
2. 数学建模过程

2.1 半周长线长模型（HPWL）

半周长线长模型是一种简单且广泛使用的线长估计方法。对于一组有连接关系的电路单元，其连线接口形成的外接矩形的宽度和高度分别为最大和最小横坐标之差、最大和最小纵坐标之差。具体计算公式如下：

设一组有连接关系的电路单元连线接口坐标为 (x1,y1),(x2,y2),…,(xn,yn)，则该组连接的半周长线长（HPWL）计算公式为：

其中 xi 和 yi分别表示第i个连线接口的横坐标和纵坐标。

2.2 直线型斯坦纳最小树（RSMT）

直线型斯坦纳最小树是另一种更精确的线长估计方法。通过在连接的电路单元间插入斯坦纳点来最小化总线长。尽管RSMT提供了更精确的线长估计，但计算复杂度较高，不适用于快速布局评估。因此，本模型主要采用HPWL估计线长，并将RSMT结果作为参考标准。

3.模型实现步骤

（1）数据解析与预处理：读取并解析电路单元连线接口的数据，提取每组连接的连线接口坐标、HPWL和RSMT值。

（2）HPWL计算：根据每组连接的连线接口坐标，计算其HPWL值。

（3）差值分析：比较每组连接的HPWL值和RSMT值，计算二者差值，并分析差值分布情况。

（4）总连接线长估计：计算所有组连接的HPWL总和，作为全局布局的总连接线长估计值，并与RSMT总和进行比较。

4.差值分析

比较每组连接的HPWL值和RSMT值，计算其差值：

5. 总连接线长估计

计算所有组连接的HPWL总和，作为全局布局的总连接线长估计值：

并与RSMT总和进行比较：

综合以上模型求解结果如下所示：

问题2：图3展示了单元密度计算示意图，请以此设计一个与电路单元坐标相关的网格密度评估模型。应用问题1构建的线长评估模型，整合密度计算，建立一个数学模型，目标为:(1)最小化总连接线长;(2)满足单元密度约束。根据附件1和附件2提供的信息，应用此模型完成全局布局，输出总连接线长(HPWL)，并可视化结果(电路单元的位置)。

问题二涉及到电路单元的全局布局，目标是最小化总连接线长（HPWL），同时满足单元密度约束。为了实现这个目标，我们需要构建一个数学模型，结合附件1提供的连线数据和附件2提供的电路单元信息，计算每个单元的连接线长，并评估其网格密度。

1.符号定义

2.目标函数

我们希望最小化电路单元的总连接线长。每个电路单元的连接线长计算公式为：

其中 Ci是第i个电路单元的连线接口坐标集合。总连接线长（HPWL）为所有电路单元连接线长的和：

3.密度约束

为了满足单元密度约束，我们需要计算每个网格的密度。第 (i,j)个网格的密度计算公式为：

其中 δi,j,k表示第k个电路单元是否覆盖网格(i,j)。具体计算方法如下：

为了确保网格密度不超过阈值 ρ，我们需要约束每个网格的密度：

我们需要通过优化算法来求解电路单元的位置，使得总连接线长最小，并且满足密度约束。优化变量为每个电路单元的左下角坐标 (xi,yi)，优化目标函数和约束如下：

4.约束条件：

电路单元不能超出布局区域的边界：

网格密度约束：

5.优化算法

使用scipy.optimize.minimize函数进行优化。具体步骤如下：

（1）初始值： 将电路单元的初始左下角坐标作为优化的初始值。

（2）目标函数： 定义为最小化总连接线长。

（3）约束条件： 包括边界约束和密度约束。

通过优化算法求解出以下结果：

Total HPWL: 322522.0

问题三：除了连接线长和单元密度，布线密度也是衡量布局质量的重要指标之一。分析图4所示的网格布线密度计算模型，找出其存在的问题。针对发现的问题，提出改进方案。应用改进后的布线密度模型，计算问题2中更新后的全局布局结果的布线密度，并对结果(网格布线密度)进行可视化。

问题四：除了最小化总连接线长和满足单元密度约束外，希望网格布线密度的最大值越小越好。请在问题3的基础上，修正问题2所建立的数学模型。根据附件1和附件2提供的信息，应用修正后的模型完成全局布局，输出总连接线长(HPWL)，并可视化结果(电路单元的位置和网格布线密度)。

可视化结果如下：

 成品文章内容：共计1w字，47页

超大规模集成电路自动布局优化模型研究

摘 要

随着集成电路设计的日益复杂化，传统的手工设计方法已无法满足现代超大规模集成电路(VLSI)设计的需求。电子设计自动化(EDA)工具在其中扮演了至关重要的角色。本文聚焦于VLSI设计中的自动布局问题，探讨了全局布局和详细布局的优化策略，以实现更高效的电路设计和更低的功耗。本研究建立了多个数学模型，对电路单元的布局进行了深入分析和优化。

针对问题一，我们建立了基于电路单元连线接口坐标的线长评估模型，解决了传统HPWL估计在多连线接口情况下偏小的问题，并通过与RSMT的比较，验证了模型的准确性和有效性。

针对问题二，我们开发了网格密度评估模型，并结合问题一中的线长评估模型，构建了一个综合数学模型。该模型最小化了总连接线长，同时满足了单元密度约束，并通过优化算法完成了全局布局，得到了电路单元位置的优化布局方案。

针对问题三，我们分析了现有网格布线密度计算模型的不足，并提出了改进方案。通过引入连线权重和改进网格划分，我们建立了一个新的布线密度模型，该模型能更准确地反映实际布线情况，并通过可视化手段展示了更新后的布线密度。

针对问题四，我们在问题三的基础上进一步优化了布局，修正了问题二中的数学模型，不仅最小化了总连接线长和满足了单元密度约束，还优化了网格布线密度的最大值，实现了更高效、更均衡的布局效果。

关键词：超大规模集成电路；自动布局；电子设计自动化；全局布局；布线密度；优化模型

代码实现：