摘要：在高速印刷电路板(pcb)设计中，逻辑门元器件速度的提高，使得pcb传输线效应成了电路正常工作的制约因素。对传输线做计算机仿真，可以找出影响信号传输性能的各种因素，优化信号的传输特性。采用全电荷格林函数法结合矩量法提取高速pcb传输线分布参数并建立等效时域网络模型，应用端接i/o缓冲器ibis瞬态行为模型，对实际pcb布线进行电气特性仿真，其结果与cadence公司的specctraquest软件仿真结果一致，且仿真效率得到提高。

　　关键词：传输线；全电荷格林函数法；高速印刷电路板；ibis模型
　　
　　随着技术的进步，目前高速集成电路的信号切换时间已经达到几百皮秒(ps)，时钟频率也已达到几百兆赫兹(mhz)，如此高的边沿速率导致印刷电路板上的大量互连线产生低速电路中所没有的传输线效应，使信号产生失真，严重影响信号的正确传输。因此有必要对高速印刷电路板(pcb)传输线建模并进行计算机仿真，这样一方面可以确定高速信号传输线的时域网络模型，另一方面也可以通过仿真找出影响信号传输性能的各种因素，以便采取措施，尽可能地优化信号的传输特性，保证系统的可靠及高性能工作。

　　当前电路工作频率不断提高，当其达到一定程度后，系统的波特性必然变得十分明显。在pcb设计中传输线的尺寸较大，其波特性应首先考虑。对传输线的分析必须采用l、c、r、g分布参数模型，这样系统的电特性分析和电磁场分析密切相关。基于这种模型，应对传输线的布局布线进行分析和仿真，由此来指导pcb设计。文中就是采用全电荷格林函数法结合矩量法提取传输线的电路元件参数模型，包括集总参数和分布参数(分布电容c、电感l、电阻r和电导g)，在建立传输线等效时域模型及提取参数的基础上进行电路分析，应用端接i/o缓冲器ibis瞬态行为模型，对实际pcb布线进行电气特性仿真。

　　采用全电荷格林函数法结合矩量法提取分布参数

　　对单根传输线，c、l、r、g是4个基本的分布参量，由此还可导出特性阻抗、相速或相位常数等参量。其中分布电阻r由所采用的导体材料和物理特性所决定的，而分布电容参数c是最重要的分布参数，因为一旦获知分布电容参数，除r以外的其他分布参数都可以通过公式转换得到。

　　为提取多根导体的分布电容矩阵，必须首先在给定导体电位的条件下求出各根导体的自由电荷电量其求解过程应求得格林函数，而多层介质下的格林函数之所以复杂，在于介质的不均匀。界面上的极化电荷会产生附加电位，其影响将叠加到格林函数的自由空间分量上。因此不妨将自由电荷和极化电荷都作为产生电位的场源，格林函数就可看成单位点电荷(三维)或单位线电荷(二维)在介质均匀的无限空间产生的电位。矩量法即是近似地将待解函数表示为n个相互正交的基函数求和展开式，每一基函数均乘以某一系数。

对于具有多根导体的系统内的分布电容，除了要考虑每一根导体自身的分布参量，还应考虑其与其他各导体之间耦合效应的互分布参量，如图1所示。其分布参量应表示为分布参量矩阵。对n根导体进行分析，其分布电量q与电位φ的关系如下：

图1　多导体系统的部分电容

从式(1)多导体线分布电容参数的定义可知，电容参数的提取必须求解给定导体电位的静电场，它是一个偏微分方程的边值问题。通过源区解法求解，其主要问题是积分方程中的核函数——格林函数求取问题。将全电荷格林函数积分方程结合矩量法[3]以数值的方法求解，选择脉冲基函数并采用点匹配，可得到方程组

　
其中n₁表示导体和介质的分界分块数，n～n₁表示介质和介质的分界分块数，总共有n个分块。前面n1个方程表示场点所在分块均在导体和介质的分界上，方程式左边的值p_m(m=1，2，…，n1)为分块中心点的电位；后面的n-n1个方程表示场点所在分块均在介质和介质的分界，方程左边的值p_m(m=n1+1，n1+2，…，n)应为零，方程式右边的分块脉冲基函数α_m(m=n1+1，n1+2，…，n)则代表各分块上的全电荷；系数矩阵1_mn(m，n=1，2，…，n)，由公式(3)表示。

其中m=1，2，…，n1；n=1，2，…，n，x′及y′为源点直角坐标。

假设导体数量为j1，根据分布电容矩阵的定义即式(1)，可依次对j1块导体中的每一块赋以单位正电荷，其余导体电位为零，解出式(2)，求得各分块的全电荷，然后将同一导体上的分块进行组合，可得到各导体上的