计算方法7 电阻网络

越发觉得自己物理学得烂了...

前置

前置的物理定律包括如下两条:

  1. 欧姆定律,即 ϕiϕj=Ui,j=Ii,jR。定义电导率为 w=R1,那么有 Ii,j=Ui,jw
  2. 基尔霍夫定律,即对于电阻网络的任意节点 v,流入的电流等于流出的电流 。

电路网络与 L

对于单位电阻组成的电路网络 G,其拉普拉斯矩阵 L 可以通过上面两条物理定律和电流联系起来。

考虑电路网络内部的节点 x,根据基尔霍夫定律有

bx=yN(x)Ix,y=yN(x)(ϕxϕy)w

为了方便讨论,一般会规定电源电势为 1V,或流入电路网络的电流总量为 1A,此处采用第二个约定。

对于单位电阻有 w=1,此时上述方程组可以写成

Lϕ=b

其中 bx 表示流入 x 的电流(流入为正,流出为负)。将电源接在图上任意两点间(不妨设为 s,tttt),其含义为向量 b 满足 bs=1,bt=1,其余均为 0

再考虑欧姆定律,有

Bϕ=I

其中 B 为图 Gn×m 关联矩阵,形如 [11],第 i 列表示边 ei 关联哪两条边,正负表示方向。m 维向量 I 表示每条边上电流的流量。

如果要考虑非单位电阻的矩阵,那么需要引入 m×m 的对角阵 W=diag\setwe1,we2wem 来分别建模每条边的电导率,在需要的时候乘上就行了。

考虑 L 的另一形式有

L=eE(G)Le=eE(G)webebe=BWB

因此还可以写成

Lϕ=BWBϕ=BWI=b

这也是很直观的,对边上的电流分布做一次图上的按邻居求和,就能得到一个节点上的全局电流分布 b

电路方程的解

定理:

Lϕ=b 有解当且仅当 b1

""

注意到 L 实对称,取一组由 1 扩充而来的正交基 \Setvi,则 ϕ=a11+i=2naivi

此时 Lϕ=L(a11+i=2naivi)=i=2naiλivi,根据正交基可知 Lϕb

直观含义为电阻网络流入的电流要等于流出的电流。

""

b1,那么 b=i=2nbivi

此时取 ϕ=i=2nbiλivi 即为一个解。

直观含义为对于一组外部电流的电势解,可以任意整体平移得到同方程的其余解(因为电流只和电势差有关)。在固定某个点电势为 0 的前提下,就能得到唯一解。

上面关于 "" 方向的证明用到了一个构造,实际上可以写成

ϕ=Lb

其中

L=i=2nλi1vivib=i=2nλivi

这意味着,当 b 是一个合法的电流(满足 b1)时,L 存在伪逆。并且 Lϕ=b 的解集为 \SetLb+k1kR

电势能和等效电阻

同样只考虑单位电阻。

考虑令 b 流入单位电流,电路网络为单位电阻,那么整个电路的等效电阻就是 s,t 间的电势差,即

Reff=ϕsϕt=bϕ=bLb

对于电势能同样可以通过等效电阻来算。注意到通的是单位电流,并且电阻为 Reff,因此电势能就是 Reff

另一种对每条边单独推导的方法如下:

E=eE(G)Ee=(x,y)E(G)(ϕxϕy)2

回忆关于 L 的二次型,有

E=ϕLϕ=Reff

并且有结论:对与任意的 s,t 流,其电势能不会比 Reff 更小。即这样的电势分布会最小化单位流的能量损耗,非常神奇的物理意义。

posted @   jjppp  阅读(746)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· Manus的开源复刻OpenManus初探
· AI 智能体引爆开源社区「GitHub 热点速览」
· 从HTTP原因短语缺失研究HTTP/2和HTTP/3的设计差异
· 三行代码完成国际化适配,妙~啊~
点击右上角即可分享
微信分享提示