电磁力中的对偶

电、磁、力中的对偶

刘红

摘要：本文从对偶的角度解释了电、磁、力之间的关系，总结了高扬提出的用于全局优

化的典范对偶理论及利用它解决非线性非凸问题的主要思路和优点。

引言

电、磁、力三大物理分支存在对偶关系。透过它们之间的不同外部现象，抽象出数学

模型，看到他们的本质却是相同的。三大系统的物理量间又存在着对偶关系，这就是典范对偶理论。非线性的变量关系或非凸性的能量函数是造成系统复杂性的关键原因。典范对偶理论旨在利用非线性变换，凸化的手段，把原空间中不便于处理的问题转化到对偶空间中来处理。这就是把“不美”的东西转化为“美”的东西，然后处理“美”的东西，最后通过能量守恒的原理把处理的结果反馈回原空间中。而三个驻点对偶定理提供了能量在原空间和对偶空间中进行的最优化的理论基础。

本文先从最简单的线性电阻电路模型开始，表示出在线性情况下的典范对偶模型。描述这种电路的数学模型是线性方程组。解这类线性方程组等价于二次规划的最优解。线性模型对应线性算子，非线性模型对应非线性算子。通过非线性变换，以及利用任何函数都可以分解为凸函数之差的方法，可将非线性非凸问题转换为线性的凸的问题。这种转换，有别于泰勒展开后取线性部分近似。这里不是近似而是变换，所以能得到更准确的效果。

1. 线性电阻电路的数学描述

考虑如图1所示的电路。此电路中，节点为1，2，3，4。令[]1234U ,,,T

U U U U =为各节点的电位，假设节点4的电位为零，[]1234f=,,,T

f f f f 分别从节点1，2，3，4流进电路的电流，设网络除节点4外没有其它的接地点，所以40f =。[]12345I ,,,,T

I I I I I =为各支路的电流，[]12345V ,,,,T

V V V V V =为各支路电阻上的电压。 各支路上电阻的电压与电流取关联参考方向。

图 1 一个电路

该电路各变量之间的关系可由下列三式描述。 由基尔霍夫电压定律可得：

12

341100001100V U b 001100

101010016t U U U U -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥=Λ+=+-⎢⎥⎢⎥⎢⎥

-⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥-⎣⎦⎣⎦ （1）

由欧姆定律可得：

112

23

34

4551/000001/000I D V 0

01/000001/000

1/R V R V R V R V R V ⎡⎤⎡⎤⎢⎥⎢⎥⎢

⎥⎢⎥

⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎣

⎦⎣⎦

（２） 由基尔霍夫电流定律可得： 123451

000111010f I 011000

1

1

1T

t I I I I I ⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥

-⎢

⎥⎢⎥=Λ=⎢⎥-⎢⎥⎢⎥⎢⎥--⎣⎦

⎢⎥⎣⎦

（３） 其中，式（１）称为代数变换关系，将各节点的电位变换为各支路电阻元件上的电压降，即仿射变换U =U b t ΛΛ+。式（２）称为对偶关系，所对应的矩阵D 称为本构矩阵，反映系统的本质。很显然矩阵D 是正定矩阵，这个矩阵确定了电压和电流的一一对应关系。一对一的关系就是“美”的关系，它常常使问题变得简单。式（3）称为平衡方程，是能量守恒（功率平衡）的必然结果。

2. 电路变量间的对偶关系

将上述三式合成，可得：f D U D b T T t t t =ΛΛ+Λ。令K D T t t =ΛΛ，f =f D b T

t -Λ，则

可得f KU =。电路各变量间的对偶关系如图2所示。

U

V

I

f

t U=U+b

ΛΛD

T t

ΛD U f

T t ΛΛ=

图 2 电路变量的典范对偶图

图2中，上面一行是在原空间中的两个向量。原空间中的内积定义为：

n

1

U ,f :U ,f n

i i

i U f =<>=∀∈∑ 。 （4）

图2中，下面一行是在对偶空间中的两个向量。对偶空间中的内积定义为：

m

1

V ;I :V ,I m

i i

i V I =<>=∀∈∑

。 （5）

由Λ，D ，T

t Λ定义的三个变换代表的三组对偶关系称为电路的典范对偶关系。从下文可以看出，典范的含义就在于对于非凸的系统或者非线性的系统，通过选取合适的变换算子，

总可以化成凸的系统， 即典范化理解为标准化、凸化。

3. 功率平衡与能量最小化

若b=0，将会有U,f V;I <>=<>。从物理的角度，可理解为功率平衡（能量守恒）。不同的空间，只是选择了不同的坐标系，也就是说选择了不同的度量方式，但无论怎么度量，

能量是不变的。 从数学的角度，根据两个向量内积的定义及矩阵乘法的结合律，易知：

T

U ,D U U D U U ;D U V ;I T

T

t t t t t t <ΛΛ>=ΛΛ=<ΛΛ>=<>。

系统的内能定义为1V ;I 2

W =

<>，对应于动力系统的动能；系统的外能定义为

U,f F =<>，对应于外力对动力系统所作的功。系统在运动中，具有动能，外力要使系统

稳定，就要对系统作功，外力所做的反功就是在消耗系统的动能。 为此，定义系统的总能量（自由能）为：11U ,K U f ,U b,D b 2

2

P W F =-=<>-<>+

<>。

系统总能量为U 的二次型。令

K U f=0U

dP d =-，即得到了平衡方程。这就说明了解电路的平衡方程可等价为

求解一个二次规划。

4. 二次规划

二次规划可描述为：

()1m in 2

T

T

P u u A u f u =

- （６）

其中A 为对称阵。如果有约束，则可以通过lagrangian 乘子法松弛为无约束规划。这里总假设A 为对称阵，否则用

2

T

A A +代替它，因为2

T

T T T T

A A u Au u A u u

u +==。

下面先讨论A 是正定矩阵的情况。若A 是正定矩阵，()P u 是A 的凸函数，令偏导数为零可解出()P u 唯一的最小值点1u A f -=。

事实上，正定矩阵A 可分解为T A D =ΛΛ，其中D 为对角阵，对角线元素都是正数。令v u =Λ（代数变换方程），*v Dv =（对偶方程），*T f v =Λ（平衡方程）。这三个方程合在一起，就是所谓的三典范对偶。

通过典范对偶的转化，原二次规划问题可转为问题：

()()()*

***1*1min max ,2T T T

u v L u v u v v D v f u -⎧⎫=Λ--⎨⎬⎩⎭

（7） 其中，()()*

*

*1*1m ax 2T T v u v v D v -⎧⎫Λ-⎨⎬⎩⎭

在*v Dv D u ==Λ处取到。函数()*,L u v 关于*v 是凹函数，关于u 是线性函数。上述问题的最优解在()*,L u v 鞍点处取到。

u

v

*

v f

<,>

u f *<;v v >

Λ

D

T Λf

T D u ΛΛ=

图 3 二次规划的典范对偶图

若A 不是正定矩阵，原二次规划不是凸规划，如果直接在原空间中求解，问题会变得

麻烦。为此，可将A 分解为两个正定矩阵之差A B C =-（任何实数可以分解为两个正数之差，任何对称矩阵都可以分解为两个正定矩阵之差，任何函数都可以分解为两个凸函数之差）。这样，原问题变成为：()1

1m in 22T T T

P u u Bu u C u f u ⎛⎫=

-+ ⎪⎝⎭

。