一类任意m×n阶矩形网络的电特性*

谭志中 谭震

1) (南通大学物理系,南通 226019)

2) (南通大学信息科学技术学院,南通 226019)

任意矩形电路网络中的电位分布问题一直是科学研究的难题.本研究发展了研究电阻网络的递推-变换(RT)理论使之能够用于计算任意m×n阶电路网络模型.研究了一类含有任意边界的m×n阶矩形网络的电位分布及等效电阻,这是一个之前一直没有解决的深刻问题,因为之前的研究依赖于规则的边界条件或一个含有零电阻的边界条件.其他方法如格林函数技术和拉普拉斯矩阵方法计算电位函数比较困难,研究含有任意边界的电阻网络也是不可能的.电位函数问题是自然科学和工程技术领域研究的一个重要内容,如拉普拉斯方程的求解问题就是其中之一.本文给出了含有一条任意边界的m×n矩形电阻网络的节点电位函数解析式,并且得到了任意两节点间的等效电阻公式,同时导出了一些特殊情形下的结果.在对不同结果的比较研究时,得到了一个新的数学分式恒等式.

1 引 言

德国物理学家基尔霍夫(1824—1887)创立的节点电流定律和回路电压定律为人类研究电路网络模型提供了理论基础.现在电阻网络的研究已经不再局限于电路领域,已经成为一种基本模型被应用于许多学科领域,通过构建电阻网络模型进行计算机仿真研究已经成为解决一些科学问题的基本方法[1,2].复杂电阻网络的等效电阻解析式的获得通常比较困难,原因在于它是一个跨学科的科学难题,不仅需要电路理论知识,而且需要数学理论与方法的创新[3−6].在自然界中,石墨烯网络的研究、一些金属化合物晶体或非金属晶体结构的研究、纳米碳管结构的研究、星际空间的宇宙尺度研究、数学中的图论研究、聚合物中的高分子结构研究、统计物理学和复杂性科学研究等,都可能需要通过构建电阻网络模型进行模拟研究.

电位分布问题是自然科学和工程技术领域研究中的一个重要内容,如拉普拉斯方程是势函数的一种二阶偏微分方程,广泛应用于电学、磁学、力学及热学的研究与计算.除此之外,还有许多服从拉普拉斯方程的势场,尽管各种势场的物理本质不同,但如果具有相似的边界条件,则因拉普拉斯方程解的唯一性,任何一个势场的解或该势场模型中实验测绘的等热面或流线图,经过对应物理量的换算之后,都可以通用于其他的势场[6].求解复杂边界条件下的拉普拉斯方程往往比较困难.在实践研究中研究者发现可以采用模拟法通过建立电阻网络模型来研究温度场、扩散场等问题,这是一种很有意义的模拟法.寻找含有任意边界m×n阶电阻网络电位问题的精确解是重要的但却是困难的.以前的理论没有解决这样的问题,本研究所建立的研究理论与方法能够求出多种结构电阻网络模型中各个节点的电势,从而解出所要求解的函数方程.因此电阻网络模型的建立与研究具有重要的理论意义与实际应用价值.

纵观电阻网络模型的研究历史[1−43],电阻网络研究的重要进展可以概括为四个阶段.第一阶段是1845年德国物理学家基尔霍夫(1824—1887)创立了节点电流定律和回路电压定律,建立了电路分析的基本理论,推动了电路理论的研究与发展;第二阶段是2000年Cserti[3](匈牙利罗兰大学) 建立了计算无限矩形网络等效电阻的格林函数技术,该方法不适用于有限电阻网络,格林函数技术已经得到较好的应用与发展[4,5,7−12];第三阶段是国际著名理论物理学家Wu[13](伍法岳,美国东北大学教授)于2004年建立了计算等效电阻的Laplacian matrix方法,第一次系统地研究了多种拓扑结构的任意电阻网络的等效电阻,给出了一系列理想的结论.但是Laplacian matrix方法只适用于自由边界(规则边界)、周期边界,零电阻边界等特殊情形的网络,所给出的结论是由倍求和表达的结果[13−19],在后来的研究中Laplacian matrix方法在研究电阻方面已经得到较好的应用与发展[20−23];第四阶段是2011年谭志中(南通大学)建立了一种“递归-变换” (recursion-transform)方法[24],该方法适用于研究含有任意边界的有限和无限的电阻网络模型,弥补了 Laplacian matrix 方法的不足[25−30].“递归-变换”方法 (简称RT方法)的新发展是Tan[28−30]于 2015年发表的三篇论文,该方法不仅能够计算等效电阻而且能够研究复阻抗的电特性[31−34].特别需要说明的是,随着研究的不断深入与发展,谭志中的RT方法已经细化为RT-I方法和RT-V方法,其中RT-I方法是指基于支路电流参数建立矩阵方程的方法,而RT-V方法是指基于节点电压参数建立矩阵方程的方法[35,36],这两种方法构成了研究电阻网络模型的系统理论.

综上所述,目前能够有效地计算大规模n×m有限电阻网络的方法只有两种,一种是伍法岳的Laplacian matrix方法,另外一种是谭志中的RT方法,而能够用于研究含有任意边界电阻网络的方法只有唯一的RT方法.因此有必要进一步介绍RT方法曾经给出的研究.2013年Tan等[25]应用RT方法研究了一个3×n阶cobweb模型,并且提出了任意m×n阶cobweb等效电阻的一个猜想,该猜想引起了国际同行专家的高度兴趣[20],吸引了国际同行专家与我们开展合作研究[23,26,27].2014年以来我们采用RT方法发表了一系列学术论文[26−43],解决了一系列之前没有解决的电阻网络难题.RT方法的优越性在于应用此理论所给出的等效电阻结论是由单求和表达的简单结果 (注:Wu教授建立的Laplacian matrix方法所给出的结果是由倍求和表达),便于推广应用,并且所得到的等效电阻公式能够有效地应用于复阻抗网络研究[32,33,37−40].尽管前面阐述RT方法能够用于研究复杂的电阻网络模型,但是由于时间的限制(RT理论提出的时间不长)及电阻网络边界条件的复杂多样,还没有来得及解决所有的问题.例如之前的RT-I理论所研究的电阻网络都依赖于含有零电阻的边界[25−31,41−43]或只是研究了特殊节点间的等效电阻[32−34].当然,研究复杂边界的电阻网络问题还需要技术创新,本文正是获得了技术创新才解决了非零边界的电阻网络,解决了以前没有解决的复杂问题.

考虑图1所示的电阻网络,该模型不含有零电阻边界但是含有一个任意的右边界,研究任意节点的电位函数是一个以前从来没有解决的问题,这主要因为以前的相关理论不能解决这样的问题(如Laplacian matrix 方法[13],格林函数技术[7−12]等).任意右边界条件代表了网络的多功能特性,因为可以任意改变右边界而得到各种各样的几何结构,如当R1=0 时得到一个扇形网络模型,当R1=r0时得到规则的矩形网络模型,等等.文献[13]曾经研究了当R1=r0时的规则矩形电阻网络(r1=r0),为了便于本文的比较研究,不妨首先给出文献[13]给出的一个结论.

图1 含有一个任意边界的矩形电阻网络,右边界上的任意电阻为r1,其他在水平和竖直方向的电阻元素分别为r和r0Fig.1.An arbitrary rectangular m×n resistor network with an arbitrary boundary.Bonds in the horizontal and vertical directions are resistors r and r0 except for r1 on the right boundary.

考虑一类规则的m×n阶矩形电阻网络模型,如图1所示(r1=r0),其中水平和竖直网格上的电阻元素分别为r和 r0,水平和竖直方向上的网格数分别为n和m.假设矩形网络的底边为X轴,左边界为Y轴.利用Laplacian matrix方法得到任意m×n阶矩形网络二节点间的等效电阻公式:

(1)式的获得是一次重要创新,但美中不足的是该公式是由倍求和表达的结果,略显复杂.本研究拟研究图1网络模型的等效电阻与电位函数解析式,所给出的结论都是由单求和表达的结果,实现了理论和方法上的创新.

2 RT-I理论与应用

2.1 参数定义

为了方便研究与计算,以及下文各种表达式的简化,给出如下定义:

以上参数是本文所使用的重要参数,适用于整个文章.这些参数的定义参考了文献[29,36]中给出的定义,可以与国际学术研究所使用的符号定义相统一.

2.2 总体设计及计算思路

考虑图1所示的任意m×n阶矩形电阻网络模型,其中水平和竖直方向上的网格数分别为n和m.水平和竖直网格上的电阻元素分别为r和r0,右边界的任意电阻为 r1.假设矩形网络的底边为X轴,左边界为Y轴.本文的研究目标是给出任意节点的电位函数解析式及任意节点间的等效电阻公式.

设 A0是直角坐标系的坐标原点.假设稳恒电流J从节点 d1(x1,y1) 输入再从节点 d2(x2,y2) 输出.用图2子网络表达所有网络中的分布电流及参数符号.设经过所有m+1行水平电阻r上的电流分别为经过所有n+1行竖直电阻r0上的电流分别为设任意节点d(x,y) 的电位函数为U(x,y),选择坐标原点 O (0,0) 为电位参考点(因为电位函数是一种相对标量值),计算图1中任意节点的电位分布函数的基本思路如下.

图2 含有电流参数和方向的部分电阻网络Fig.2.Segment of resistor network with current directions and parameters.

2.3 递推方程建立

矩形网络的子网络如图2所示.根据RT-I理论应用基尔霍夫定律分析该网络,利用图2中的2个网孔可以得到相关回路电压方程,根据图2中的6个节点可以得到6个节点电流方程.通过消元法消去水平方向的电流参数而得到如下方程:

其中 h=r/r0.如果考虑电流J从 d1(x1,y1) 输入和从 d2(x2,y2) 输出的情形,则上述方程可以写成矩阵形式:

其中 Ik和Hx分别为 m×1 阶列矩阵,表达为

其中[ ]T表示为矩阵的转置,(Hk)i是矩阵 Hx中的元素伴随着电流 J从 d1(x1,y1) 输入和从d2(x2,y2)输出,δ i,y 的含义是 δ i,y|y=i=1,δ i,y|y/=i=0,

接下来考虑建立左边界和右边界的电流条件方程,相似于方程(11)的建立方法,应用基尔霍夫定律分别对左边界和右边界进行网络分析,获得如下2组约束条件方程:

其中 h1=r1/r0,E 是一个 m×m 阶恒等矩阵,矩阵 Am由方程(14)给出.由于右边界含有任意电阻,因此关于右边界的约束方程(16)比左边界的约束方程(15)复杂得多.

以上方程(11)—(16)就是计算图1所示任意m×n网络的等效电阻所需要的全部方程(属于RT方法中的递推方法),所有电阻网络中的电位分布和等效电阻等等问题的解决都依赖于以上方程的解.然而,从方程 (11)—(16)中直接解出电流的精确表达式是比较困难的事情.为此这里对以上诸方程采用变换方法(属于RT方法中的变换技术)间接地解决该问题.

2.4 矩阵变换方法

如何获得矩阵方程(11)的通解是解决问题的关键.这里首先采用矩阵变换方法将复杂的矩阵方程转变成为一维的差分方程而间接地给出它的通解.具体方法是用一个待定的m×m阶矩阵Pm左乘矩阵方程(11),得到

通过下列2个恒等式建立对角化矩阵变换:

其中Tm=diag(t1,t1,···,tm)是一个对角化矩阵,并且 ti(i=1,2,···,m)是矩阵 Am的特征值.解(18)式得到

其中 θi=iπ/(m+1),将 (20)式代入 (19)式计算解得

一个简单的计算表明新矩阵 Pm是可逆的,其逆矩阵为(可以参考文献[33]的研究)

为了研究的需要,根据(17)式和(19)式给出如下定义:

其中Xm是一个m×1阶列矩阵,

因此,将方程(13)及(17)式应用于(11)式得到,

将矩阵(21)中的元素应用于方程(13)及(25)式计算得到(简记 ζ xk,i 为 ζ k,i)

显然方程(25)是一个简单易解的线性方程.

接下来实施与(17)式所做变换相同的变换方法,应用矩阵 Pm左乘矩阵(15)式和(16)式实施矩阵变换,并且应用(19)式和(23)式得到,

以上就是RT-I理论中的矩阵变换方法,有了以上矩阵变换就可以解出含有电流参数的变量,进而获得所有支路电流函数的解析式,当然实际的求解还需要相关的计算技巧.下文将采用RT理论中的解方程技术给出以上方程组的通解与特解.

2.5 矩阵方程的解

考虑电流J从 d1(x1,y1) 输入和从 d2(x2,y2) 输出,假设是方程 (25) 关于Xk的特征方程的2个根,那么根据(25)式的特征方程可以解得方程的特征根(4).解方程(25)可以得到如下分段函数解:

显而易见以上方程组的解是一个复杂的问题.为了计算方程组的解,首先采用消元法消去一些变量而降低方程组的数量.将(27)式代入(29)式消去得到

由方程(30)和(31)及在方程(34)中取k=x1-1,x1解得

再由方程(32),(33)及(35)解得(x2≤k＜n)

其中定义 qi=ti+(h1-2).

矩阵方程(38)并不复杂,经过适当的计算之后解得

方程(39)是一个关键的公式.将(39)式分别代入(34)式—(36)式及(28)式,经过一系列化简运算之后得到

其中 h1=r1/r0,并且 fx(i) 是一个新定义的函数,

2.6 计算分支电流

因为

将此代入(44)式得到

其中 fs(i) 定义在方程(42)中.

另外,应用基尔霍夫节点电流定律可以计算出沿着下边界轴线方向的电流,即

将(43)式代入(46)式计算得到

为了研究需要,需要给出(47)式的分段表达式.

当 0 ≤k≤x1时,将 (34)式代入 (47)式化简得到

进一步对(48)式中的k从 k=1 到x求和,得到

其中 fx(i) 由(42)式给出,以下与此同.

当X1≤k≤x2时,将 (35)式代入 (47)式化简得到

进一步利用 (49)式对 (50)式中的 k从 k=1 到x 求和,得到

当X2≤k≤n 时,将 (36)式代入 (47)式计算得到

进一步利用 (51)式对 (52)式中的 k从 k=1 到x 求和,得到

(45)式和(50)式—(53)式即为电路网络中的各支路电流分布理论公式.获得是计算电位函数的预备工作,下面将依据这些结果推导电位分布函数公式.

3 电位函数解析式

3.1 电位分布函数公式

设图1网络中的 d1(x1,y1)和d2(x2,y2) 为电流输入和输出的2个任意节点(不失一般性设x1≤x2).设一个任意节点 d (x,y) 的电位函数为U(x,y),选择坐标原点 d0(0,0) 的电位 U (0,0) 为电位参考点,则矩形网络中的任意节点的电位函数解析式为

特别地,当取 (x,y)=(0,0) 代入 (54) 式时,得到坐标原点 d0(0,0) 的电位

(57)式即为本文所选择的电位参考点,即该节点d0(0,0)位于坐标原点的参考电位为方程(57).显然(57)式的选择是科学合理的,并且与(54)式—(56)式的表达结构完全自洽.

3.2 电位函数证明

选择 d0(0,0) 为电位参考点,设应用 (45)式和 (49)式—(53)式,以及(8)式和(9)式推导计算电位公式.

a) 当 0 ≤x≤x1时,将 (45)式和 (49)式代入(8) 式,得到两节点 d0(0,0)和d (x,y) 之间的电位差:

b) 当X1≤x≤x2时,将 (45)式和 (51)式代入(8)式得到两节点 d0(0,0)和d (x,y) 之间的电位差:

c) 当X2≤x＜n 时,将 (45)式及 (53)式代入(8)式和(9)式得到两节点 d0(0,0)和d (x,y) 之间的电位差:

其中 fx(i) 由(42)式给出,并且以上三式利用了下列恒等式

因为电位分布函数是一个相对值,根据方程(57)的假设及(42)式的定义,则坐标原点 O (0,0) 的电位函数为

将(61)式分别代入 (58)式—(60)式并且利用(42)式化简即得到任意位置节点 d (x,y) 的电位分布函数解析式(54)式—(56)式.

(54)式—(56)式是本文第一次发现的关于矩形电阻网络的电位函数解析式,是之前研究者一直没有解决的难题,本文应用RT-I方法彻底解决了该问题,第一次获得了含有任意边界的m×n矩形电阻网络的电位函数解析式,为电阻网络研究提供了新的理论工具.(54)式—(56)式是三个一般性结论,为了帮助读者易于理解我们的结论,下面给出一些具体应用.

3.3 电位函数的应用

特别说明,下面各种应用中出现的网络均是基于图1结构的网络模型,并且选择位于坐标原点d0(0,0)的参考电位为方程(57).

应用1输入输出电流分别在左右边界上的情形.设电流从电路网络的左边界上的节点d1(0,y1)输入而电流从右边界的节点 d2(n,y2) 输出,当 0 ≤x≤n 时,由 (55) 式得到任意节点 d (x,y) 的电位分布函数

这是因为当时,有x1=0,x2=n

特别地,当 h1=1⇔r1=r0时,有则 (62)式简化成为一个简单结论

当 h1=1⇔r1=r0时,输入和输出节点电流d1(0,y1)和 d2(n,y1) 位于同一水平线上的两端时,由方程(63)得到任意节点 d (x,y) 的电位函数,

应用2输入输出电流在对角线的情形.设电流从电路网络的左下角的节点 d1(0,0) 输入而从右上角的节点d2(n,m) 输 出,当 0 ≤x≤n 时,由(55)式得到任意节点 d (x,y) 的电位分布函数,

这是因为当 y1=0,y2=n 时,有则将此代入(55)式即得到(65)式.

应用3输入输出电流在同一竖直线上的情形.设电流从电路网络的同一竖直线上的节点d1(x1,y1)输入而从上方的节点 d2(x1,y2) 输出,当0≤x≤n 时,由(54)—(56)式得到任意节点d(x,y)的电位函数,

特别地,当输入和输出电流位于同一竖直线上的上下端点时,因为所以由 (66)式得到任意节点d(x,y)的电位分布函数,

应用4输入输出电流在左边界上的情形.设电流从电路网络的左边界上的节点d1(0,y1) 输入而电流从节点d2(0,y2)输出,当0≤x≤n时,在(56)式中取x1=x2=0 得到任意节点d (x,y) 的电位分布函数,

特别地,当输入和输出电流位于左边界的上下端点时,由 (68)式得到

应用5含有2个节点输出电流的情形.设有电流大小为J的电流从左边界上节点d1(0,y1) 输入,再分别有电流大小为J/2的电流分别从右边界上的节点d2(n,0)和d3(n,m)的输出.这种情形视为2个独立电流源的作用.那么,2个电流输入时,任意节点的电位是2个电位之和,应用(62)式得到节点d(x,y) 的电位分布函数,

其中利用了

应用6含有4个节点输入和输出电流的情形.设分别有电流大小为J的电流从左边界的上下端点输入,再分别有电流大小为J的电流从右边界上的上下端点输出.这种情形视为2组独立电流源的作用.那么,2 组电流输入时,任意节点的电位是2个电位之和,应用(62)式得到节点d(x,y) 的电位分布函数,

(71)式是由矩形网络四个角点输入与输出电流时的电位分布函数.

应用7右边界含有m+1个输出电流的情形.设有电流大小为J的电流从左边界节点d1(0,y1)输入,设U (0,0)=U0,定义在 (57)式中,并且电流分别从右边界上的m+1个节点dk(n,yk)输出,每个节点输出的电流均为J/(m+1).则任意节点d(x,y) 的电位分布函数为

证明当存在m+1个电流源时,由于电位是标量函数,则任意节点的电位是由m+1个电位之代数和组成,应用(62)式得到节点d(x,y) 的电位分布函数

应用8含有m+1个输入和输出电流的情形.设分别有电流大小为J/(m+1) 的电流从左边界上m+1个节点dk(0,yk)输入,重新定义U(0,0)=U0,U (n,0)=-U0,并且电流分别从右边界上的m+1个节点dk(n,yk) 输出,每个节点输出的电流均为J/(m+1).则任意节点d(x,y) 的电位分布函数为

证明当存在m+1个电流源时,由于电位是标量函数,则任意节点的电位是由m+1个电位之代数和组成,应用(62)式得到节点 d (x,y) 的电位分布函数

注意,在该命题条件下应用(57)式得到 U(0,0)=0,即(76)式是以 U(0,0)=0 为参考点为条件的结果.当考虑 U (0,0)=U0,U (n,0)=-U0时,(76)式应该改写成

考虑 U (n,0)=-U0,从而由 (77)式得到将此代入(77)即式得到(74)式.

注意:(74)式是一个与y无关,仅仅与x有关的函数.说明在命题条件下矩形网络在同一竖直轴线上的电位都相同,这是一个很有趣的问题.

4 矩形网络的等效电阻公式

4.1 一个总的等效电阻公式

在图1所示的含有一个任意右电阻边界条件的m×n电路网络中,其中水平和竖直方向上的网格数分别为n和m.水平和竖直网格上的电阻元素分别为r和 r0,右边界的任意电阻为 r1.设矩形网络的底边为X轴,左边界为Y轴.则任意两节点d1(x1,y1)和d2(x2,y2) 间的等效电阻公式为

(78)式的推导过程如下.

由电位函数公式计算图1电路网络中任意两节点d1(x1,y1)和d2(x2,y2) 之间的等效电阻,应用欧姆定律得到

在 (55)式中分别取 (x,y)=(x1,y1)和(x,y)=(x2,y2)得到

将方程(81)和(82)代入(80)式化简即得到(78)式.

(78)式是图1所示的含有一个任意右边界条件的m×n电路网络的一个总的等效电阻公式,由本文第一次获得.由于右边界是一个任意电阻,并且节点d1(x1,y1)和d2(x2,y2) 是两个任意节点,所以(78)式包含了一系列特殊情形的等效电阻,作为公式的应用,下面将给出(78)式的若干有趣的推论.

4.2 电阻公式(78)式的应用

Case 1一个半无穷网络的电阻公式.在图1所示的含有一个任意右电阻边界条件的m×n电路网络中,当n→∞,x1,x2→∞但x1-x2为有限值,并且竖直轴方向的网格数也为有限时(上下有限,左右无限),称该网络为一个半无穷矩形网络,其任意二节点d1(x1,y1)和d2(x2,y2) 间的等效电阻为

(83)式是根据(78)式取极限得到的.

Case 2当h1=1 (r1=r0)时,图1 的网络退化为一个规则的矩形网络,根据(78)式得到图1所示的标准边界的矩形网络任意二节点d1(x1,y1)和d2(x2,y2) 之间的等效电阻

Case 3当h1=0 时,图1的网络退化为扇形网络,如图3所示.此时右边的零电阻边界塌陷为一个节点,根据(78)式得到任意两节点d1(x1,y1)和d2(x2,y2) 之间的等效电阻公式

图3 一类任意m×n阶Fan电阻网络模型Fig.3.An arbitrary m×n Fan resistor network.

特别说明:文献[27−29]曾经研究了Fan网络的等效电阻,但是这里给出的电阻公式与他们的结果不同,因为文献[27−29]研究Fan网络时是沿着经度(半径)方向建立矩阵方程,而本文是沿着纬度(圆弧)方向建立矩阵方程.这说明沿着不同轴线建立方程可以得到不同的等效电阻表达式.

Case 4当两个节点都处在相同竖直轴线上时,根据(78)式得到图1网络中任意两节点d1(x,y1)和d2(x,y2) 之间的等效电阻公式

特别地,当节点 d1和d2分别在相同轴线上的上下边界上时,(86)式可以进一步简化为

说明:文献 [33]曾经专门研究了Case 4的问题,并且得到的结果与(86)式完全一致,但是这里给出的电阻公式是一般结果(78)式的特殊情形,这也就间接验证了本研究的正确性.

Case 5当两个节点在相同的水平轴线上时,根据(78)式得到图1网络水平轴线上任意两节点d1(x1,y)和d2(x2,y) 之间的等效电阻公式

Case 6当h1=1,节点 d1=(0,y1) 位于左边界上,并且节点 d2=(n,y2) 位于右边界上,根据(78)式得到图1网络中任意两节点d1(0,y1) 和d2(n,y2)之间的等效电阻公式

Case 7当d1=(x1,0)位于底边界上,而d2=(x2,m)位于上边界上,根据(78)式得到图1网络中任意两节点d1(x1,0)和d2(x2,m) 之间的等效电阻公式

Case 8当h1=1 时,并且节点d1=(0,0) 和d2=(n,m)分别位于一对对角线上时,根据(78)式得到图1网络中任意两节点d1(0,0)和d2(n,m) 之间的等效电阻公式

请注意(91)式是一个规则连接的m×n矩形网络最大分离节点间的等效电阻公式.这是一个很有趣的问题,该问题曾经被文献[15,16]深入研究,但是这里给出的(91)式比较简单,更加有利于进行渐进性问题的研究.

以上特例虽然比较好地解释了(78)式,但是这些公式仍然比较复杂而不容易理解.为了帮助读者进一步理解和验证所得结论的正确性,下面将给出(78)式在 m=1 时的一个最简单情形下的结果.

Case 9当m=1 时,图1 退化为一个含有任意右边界的 1×n 电阻网络,如图4 所示.如此情形下有θi=iπ/(m+1)=π/2,并且

设节点Ak和Bk是图4 中的任意节点,其中k表示从左向右计数的第k个网格,根据(78)式得到如下一系列结论.

当A0=(0,0),Bk=(k,1),根据 (78)式得到

图4 一类含有任意右边界的 1×n 阶电阻网络Fig.4.1×n resistor network with an arbitrary right boundary.

其中0≤k≤n,参数ΔFk及定义在方程 (6) 中,即αx=ΔFx+(h1-1)ΔFx-1.在h1=1 的条件下,根据(93)式和(94)式得到

特别地,当h1=1时,分别设k=0和n,根 据(95)式和(96)式得到

(97)式—(99)式的正确性可以从文献[24]中的结论得到验证.当然,当分别取 n=0,1,2 时可以直接验证以上诸公式的正确性.事实上,由于本文的推导过程和计算过程都是精确的和自洽的,因而所得结论必然是正确的.从以上 1×n 阶电阻网络的研究发现:一个简单的 1×n 电阻网络都能给出 7个不同情形的等效电阻公式,显然一个任意m×n电阻网络的等效电阻公式(78)式是一个多功能的普适公式,具有重要的应用价值.

Case 10当h1=1,n=0 时,图1 退化成为一个线性串联电阻结构.根据(78)式得到

其中θi=iπ/(m+1).众所周知线性串联电阻d1(0,y1)和d2(0,y2) 之间点间有|y2-y1| 个电阻元素,其等效电阻为

比较(100)式和(101)式得到

(102)式是一个很有趣的恒等式,该恒等式可以利用文献[13]的理论给予证明,但是这里是采用物理方法而不是数学方法得到的.

4.3 一个新的分式恒等式

(1)式是由文献[13]在R1=r0情形下给出的矩形网络的等效电阻,而本文给出了相同条件电阻网络相应的等效电阻公式(84)式,既然是同一个矩形电阻网络模型,在相同的坐标规定下两者的等效电阻必然相等.比较(1)式与(84)式得到

其中m,n,x1,x2及y都是自然数,并且以及

(103)式是一个很有趣的分式恒等式,本文是第一次发现,而且是采用物理方法而非数学方法得到的复杂分式恒等式,所以是一个很有意义的工作.这一恒等式的发现为数学工作者提供了新的数学工具.由于(103)式比较复杂,下面将给出它的一些特殊推论.

推论1当y1=y2=y 时,φj=jπ/(n+1),θi=iπ/(m+1).由 (103)式得到

推论2当x1=x2=x 时,φj=jπ/(n+1),θi=iπ/(m+1).由 (103) 式得到

推论3当m=1,y1=0,y2=1 时,得到θi=π/2,φj=jπ/(n+1),由 (103) 式得到

其中λi,由(104)式简化得到

推论4当m=1,y1=0,y2=1 并且 x1=x2=x 时,φj=jπ/(n+1),由 (107)式得到

另外,如果在(103)式中继续取一些特殊坐标(xi,yi)值,还可以推导出一系列有趣的简单恒等式,这里不再举例.

5 结 论

进一步发展了之前建立的研究电路网络的RT-I理论[28−30],借助于图1的任意矩形电路网络模型阐述了RT-I方法的基本过程.本文克服了之前的RT-I理论研究电阻网络都依赖于含有零电阻的边界[25−31]或者要求规则边界[32−34]的缺点,通过技术创新解决了非零边界的电阻网络问题,解决了以前从来没有解决的图1网络问题.关键创新点是建立了计算沿着水平方向电流的方法,给出了任意的水平电流公式(48)式—(53)式.事实上,任意边界的网络才是现实中的实际电路网络模型,因此本文建立的理论推广了之前的理论而成为更具一般性的普适方法.新的RT-I理论可以计算复杂电路网络中任意节点的电位解析式,能够自然地得到电路中的电流分布,自然地得到任意节点间的等效电阻公式.

首先利用RT-I方法导出了任意节点的电位公式(54)式—(56)式,然后得到了含有一个任意边界的m×n阶电阻网络的任意节点间的等效电阻公式(78)式.(78)式是一次理论上的创新,包含了一系列的特殊结论.对于规则的m×n阶矩形电阻网络,文献[13]给出的结论(1)是由倍求和表达的结果,而本文给出的规则矩形网络的等效电阻公式(84)式是由单求和表达的结果,根据两种不同方法对同一问题得到的结论的等价性,得到了一个有趣的分式恒等式(103)式,问题的有趣之处就在于这是通过物理方法而不是数学方法得到的,这为数学研究者提供了一个新的研究课题.作为电阻网络公式(78)式的应用及对公式的理解,本文给出了一系列有趣的特殊推论.