配平化学方程式的基本原理与方法.pdf

配平化学方程式

kkew3

2022年10月14日

1基本原理

将化学方程式转化为系数矩阵，配平即为求系数矩阵所蕴含的齐次方

程组的非平凡解。存在非平凡解的充要条件是系数矩阵零空间的维度大于

零。非平凡解若存在即为无穷多个，但线性无关的非平凡解的数目一定为零

空间的维数。要想求得线性无关的非平凡解，只需通过高斯消元法得到零空

间的一组基即可。

2简单解法

可以直接用高斯消元法（例如全选主元法）解出零空间的整数基。这适

用于几乎所有只有一种配平方式的方程；然而对于有多种配平方式的方程，

该方法似乎无能为力，因为基向量很可能含有非正数，且没有有效方法求出

满足锥条件（见下文）的基变换。

3进一步讨论

然而化学方程式配平有更多的要求：任意配平系数必须大于零；也就

是说，配平系数所构成的向量，即方程系数矩阵所蕴含的齐次方程组的解

（下称为配平解），存在于一个opencone中：x∈Rm，其中m为化学

++

方程式的项数。因此，零空间的维度大于零事实上只是存在至少一种配平

方式的必要不充分条件。如果存在多种配平解，则这几种配平解满足conic

independence。另一个比较理想（也许不是必须的）的条件是配平解各分量

均为整数：x∈Zm。

++

1

我们可以把配平问题建模成一系列优化问题。首先求解：

⊺

min1x0

x0

d(1)

s.t.{Mt|t∈R}=N(A),

Mz0=x0≻0

假设已解得j个x（即将解x），需求解：

ij

⊺

min1xj(2a)

xj

d

s.t.{Mt|t∈R}=N(A),(2b)

Mzj=xj≻0,(2c)

z⊺z∥z∥·∥z∥∀0≤ij(2d)

jiji

直到无可行解。其中，M是系数矩阵A的零空间的基矩阵，可用QR分解

求得，参见[1]第682页。如果加上配平解各分量均为整数的条件，那么可

参考这个presentation，在已知整数解上界2p的情况下，整数规划可以归

约为0-1规划问题，通过求解系列问题：

⊺

min1x0

x0

d

s.t.