线性方程约束下修正方法的可行性分析（演绎推理，一）

前言

    博客《一种线性方程约束下生成随机数修正的一般性方法（上）》、《一种线性方程约束下生成随机数修正的一般性方法（中）》、《一种线性方程约束下生成随机数修正的一般性方法（下）》给出了生成随机数不满足线性方程约束时的修正公式，该修正公式既适用于单线性方程约束，也适用于多线性方程约束，可以作为线性方程约束下生成随机数修正的一般性方法。但上述博客未能判断修正后的结果是否满足变量上下限约束和线性方程约束，即未能完成对修正方法所得结果的可行性分析。为进一步完善生成随机数修正方法的理论框架，本系列博客将首先给出单线性方程约束问题可行性定理和多线性方程约束问题可行性定理，作为分别判断线性方程约束问题的充分必要条件，然后利用上述定理研究修正结果与线性方程约束问题可行性之间的关系，通过修正结果的可行性判断线性方程约束问题的可行性，并基此设计线性约束下生成随机数修正方法的可执行程序。
    可行性分析系列博客：
    1. 《线性方程约束下修正方法的可行性分析（演绎推理，一）》
    2. 《线性方程约束下修正方法的可行性分析（演绎推理，二）》
    3. 《线性方程约束下修正方法的可行性分析（演绎推理，三）》
    4. 《线性方程约束下修正方法的可行性分析（演绎推理，四）》
    5. 《线性方程约束下修正方法的可行性分析（演绎推理，五）》
    6. 《线性方程约束下修正方法的可行性分析（演绎推理，六）》
    7. 《线性方程约束下修正方法的可行性分析（程序设计）》
    8. 《线性方程约束下修正方法的可行性分析（实例应用）》

单线性方程约束下修正方法的可行性分析

    单线性方程约束对应于博客《一种线性方程约束下生成随机数修正的一般性方法（上）》中的情况2.2.1，首先需要求证单线性方程约束问题可行性定理，并结合利用博客《一种线性方程约束下生成随机数修正的一般性方法（上）》所提修正方法得到的修正结果取值范围，探究其与单线性方程约束问题可行性定理之间的联系，进而达到通过修正解的结果即可判定方程可行性的目的。

单线性方程约束问题可行性定理

    定理1：对于博客《一种线性方程约束下生成随机数修正的一般性方法（上）》式(1)所示问题，若经系数矩阵行最简形变换后的常数向量${\mathbf{b}}_{\mathbf{RREF}}$自第$\left(R+1\right)$列到第$M$列不存在非零元素，且线性方程组对应系数矩阵的秩$R=1$，则问题有解的充分必要条件为$\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\le b_1\le \sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j$。

    证明1：首先求证必要性，由于博客《一种线性方程约束下生成随机数修正的一般性方法（上）》式(1)所示问题中，变量$x_j$存在上下界约束 x j ∈ [ x ‾ j , x ‾ j ] , j ∈ { 1 , ⋯   , N } x_j ∈ \left[ \underline{x}_j, \overline{x}_j \right], j ∈ \{ 1, \cdots, N \} xj​∈[x​j​,xj​],j∈{1,⋯,N}，因此欲使该问题有解，则所求的$x_j$必须在上下界范围内，将$x_j$乘以方程中对应的系数$a_{1,j}$，不难得到$a_{1,j}{x}_j$的取值范围为
$$\begin{array}{c}\begin{array}{ll}{a}_{1,j}{\underline{x}}_j& \left(a_{1,j}>0\right)\\ 0& \left(a_{1,j}=0\right)\\ a_{1,j}{\overline{x}}_j& \left(a_{1,j}<0\right)\end{array}\}\le a_{1,j}{x}_j\le \{\begin{array}{ll}{a}_{1,j}{\overline{x}}_j& \left(a_{1,j}>0\right)\\ 0& \left(a_{1,j}=0\right)\\ a_{1,j}{\underline{x}}_j& \left(a_{1,j}0\right)\end{array}\end{array}$$
    将$a_{1,j}{x}_j$从$j=1$到$j=N$进行求和可得
$$\begin{array}{c}\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\le \sum _{j=1}^N{a}_{1,j}{x}_j=b_1\le \sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j\end{array}$$
此时必要性可证。
    然后求证充分性，由于该问题中的系数矩阵行最简形变换后的常数向量${\mathbf{b}}_{\mathbf{RREF}}$自第$\left(R+1\right)$列到第$M$列不存在非零元素，且线性方程组对应系数矩阵的秩$R=1$，因此根据线性方程组理论可知，该问题在$N$维实数空间内必然有解，此时需求证$\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\le b_1\le \sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j$时，必能得到解向量$\mathbf{x}={\left[\begin{array}{ccccc}{x}_1& \cdots & x_j& \cdots & x_N\end{array}\right]}^T$，使得对于 ∀ j ∈ { 1 , ⋯   , N } \forall j ∈ \{ 1, \cdots, N\} ∀j∈{1,⋯,N}，都有$x_j\in \left[{\underline{x}}_j,{\overline{x}}_j\right]$。
    以下采用数学归纳法进行进一步证明：
    当$N=1$时，由于$R=1$，因此$a_{1,1}\ne 0$。若$a_{1,1}>0$，则此时$a_{1,1}{\underline{x}}_1\le a_{1,1}{x}_1=b_1\le a_{1,1}{\overline{x}}_1$；若$a_{1,1}<0$，则此时$a_{1,1}{\overline{x}}_1\le a_{1,1}{x}_1=b_1\le a_{1,1}{\underline{x}}_1$。两种情况均可求得${\underline{x}}_1\le x_1=\frac{b_1}{a_{1,1}}\le {\overline{x}}_1$，故必能在$\left[{\underline{x}}_1,{\overline{x}}_1\right]$上求得$x_1=\frac{b_1}{a_{1,1}}$，使得待求量既满足该问题中的所有约束条件，问题的充分性成立。
    假设当$N=K\left(K>1\right)$时，命题的充分性成立。则当$N=K+1$时，有
$$\begin{array}{c}\sum _{j=1,a_{1,j}>0}^{K+1}{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^{K+1}{a}_{1,j}{\overline{x}}_j\le \sum _{j=1}^{K+1}{a}_{1,j}{x}_j=b_1\le \sum _{j=1,a_{1,j}>0}^{K+1}{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^{K+1}{a}_{1,j}{\underline{x}}_j\end{array}$$
    将式(3)所示和式方程部分分解为前$K$项和式与第$\left(K+1\right)$项相加的形式，则式(3)变为
$$\begin{array}{c}\sum _{j=1}^K{a}_{1,j}{x}_j+a_{1,K+1}{x}_{K+1}=b_1\end{array}$$
    结合式(4)和题设中的充分性条件可知，若$\sum _{j=1}^K{a}_{1,j}{x}_j$满足
$$\begin{array}{c}\{\begin{array}{ll}{b}_1-a_{1,K+1}{\overline{x}}_{K+1}\le \sum _{j=1}^K{a}_{1,j}{x}_j=b_1-a_{1,K+1}{x}_{K+1}\le b_1-a_{1,K+1}{\underline{x}}_{K+1}& \left(a_{1,K+1}>0\right)\\ b_1-a_{1,K+1}{\underline{x}}_{K+1}\le \sum _{j=1}^K{a}_{1,j}{x}_j=b_1-a_{1,K+1}{x}_{K+1}\le b_1-a_{1,K+1}{\overline{x}}_{K+1}& \left(a_{1,K+1}<0\right)\\ \sum _{j=1}^K{a}_{1,j}{x}_j=b_1& \left(a_{1,K+1}=0\right)\end{array}\end{array}$$
则必然有
$$\begin{array}{c}\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\overline{x}}_j\le \sum _{j=1}^K{a}_{1,j}{x}_j=b_1-a_{1,K+1}{x}_{K+1}\le \sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\end{array}$$
    根据前述在$N=K$时的假设，必能得到解向量${\mathbf{x}}_{1-K}={\left[\begin{array}{ccccc}{x}_1& \cdots & x_j& \cdots & x_K\end{array}\right]}^T$，使得对于 ∀ j ∈ { 1 , ⋯   , K } \forall j ∈ \{ 1, \cdots, K\} ∀j∈{1,⋯,K}，都有$x_j\in \left[{\underline{x}}_j,{\overline{x}}_j\right]$。将其代入式(4)，可得
$$\begin{array}{c}{b}_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\le a_{1,K+1}{x}_{K+1}=b_1-\sum _{j=1}^K{a}_{1,j}{x}_j\le b_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\end{array}$$

    1) 当$a_{1,K+1}>0$时，式(7)变为
$$\begin{array}{c}{a}_{1,K+1}{\underline{x}}_{K+1}\le b_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\le a_{1,K+1}{x}_{K+1}\le b_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\le a_{1,K+1}{\overline{x}}_{K+1}\end{array}$$
    结合式(6)可知，此时必能得到$x_{K+1}=\frac{b_1-\sum _{j=1}^K{a}_{1,j}{x}_j}{a_{1,K+1}}$，使得${\underline{x}}_{K+1}\le x_{K+1}=\frac{b_1-\sum _{j=1}^K{a}_{1,j}{x}_j}{a_{1,K+1}}\le {\overline{x}}_{K+1}$。

    2) 当$a_{1,K+1}<0$时，式(7)变为
$$\begin{array}{c}{a}_{1,K+1}{\overline{x}}_{K+1}\le b_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\le a_{1,K+1}{x}_{K+1}\le b_1-\left(\sum _{j=1,a_{1,j}>0}^K{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^K{a}_{1,j}{\underline{x}}_j\right)\le a_{1,K+1}{\underline{x}}_{K+1}\end{array}$$
    结合式(6)可知，此时必能得到$x_{K+1}=\frac{b_1-\sum _{j=1}^K{a}_{1,j}{x}_j}{a_{1,K+1}}$，使得${\underline{x}}_{K+1}\le x_{K+1}=\frac{b_1-\sum _{j=1}^K{a}_{1,j}{x}_j}{a_{1,K+1}}\le {\overline{x}}_{K+1}$。

    3) 当$a_{1,K+1}=0$时，结合式(6)可得，式(7)变为$0=0$的恒等式，此时可在区间$\left[{\underline{x}}_{K+1},{\overline{x}}_{K+1}\right]$上任意取值均可作为$x_{K+1}$的解。

    因此可证，当$N=K+1$时，必能得到解向量$\mathbf{x}={\left[\begin{array}{ccccc}{x}_1& \cdots & x_j& \cdots & x_{K+1}\end{array}\right]}^T$，使得对于 ∀ j ∈ { 1 , ⋯   , K + 1 } \forall j ∈ \{ 1, \cdots, K + 1\} ∀j∈{1,⋯,K+1}，都有$x_j\in \left[{\underline{x}}_j,{\overline{x}}_j\right]$，命题的充分性成立。综合对$R=1$时的讨论和对$R=K$时的假设可知，命题的充分性对于$\forall R$恒成立，命题的充分性得证。
    综上所述，命题的必要性和充分性均成立，命题得证，该命题可作为对修正方法进行可行性分析时的预备定理——单线性方程约束问题可行性定理。

修正结果与单线性方程约束问题可行性定理之间的联系

    当生成的随机向量${\mathbf{x}}_{\mathbf{gen}}$不满足单线性方程约束时，可通过对$\sum _{j=1}^N{a}_{1,j}{x}_{gen,j}>b_1$和$\sum _{j=1}^N{a}_{1,j}{x}_{gen,j}<b_1$的情况分类讨论，并根据博客《一种线性方程约束下生成随机数修正的一般性方法（上）》所提方法修正为${\mathbf{x}}_{\mathbf{cor}}$。通过下述讨论可知，单线性约束问题的可行性可通过修正后${\mathbf{x}}_{\mathbf{cor}}$是否在变量上下界约束范围内的特征直接体现：当单线性约束问题可行时，修正结果${\mathbf{x}}_{\mathbf{cor}}$必然在变量上下界约束范围内；当单线性约束问题不可行时，修正结果${\mathbf{x}}_{\mathbf{cor}}$必然超出变量上下界约束范围。

    情况1： 当$\sum _{j=1}^N{a}_{1,j}{x}_{gen,j}>b_1$时，
$\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)>0$，其修正解在实数空间上的分布可分如下情况进行讨论：

    情况1.1：当单线性方程约束问题有解时，根据单线性方程约束问题可行性定理可知$\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\le b_1\le \sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j$，将其代入修正结果，可得
$$\begin{array}{c}{x}_{cor,j}=\{\begin{array}{ll}\frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\underline{x}}_j\right)+{\underline{x}}_j\in \left[{\underline{x}}_j,{\overline{x}}_j\right]& \left(a_{1,j}>0\right)\\ x_{gen,j}\in \left[{\underline{x}}_j,{\overline{x}}_j\right]& \left(a_{1,j}=0\right)\\ \frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\overline{x}}_j\right)+{\overline{x}}_j\in \left[{\underline{x}}_j,{\overline{x}}_j\right]& \left(a_{1,j}<0\right)\end{array}\end{array}$$

    情况1.2：当单线性方程约束问题无解时，根据单线性方程约束问题可行性定理可知$b_1<\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j$或$b_1>\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j$，将其代入修正结果，可得

    情况1.2.1：当$b_1<\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j$时
$$\begin{array}{c}{x}_{cor,j}=\{\begin{array}{ll}\frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\underline{x}}_j\right)+{\underline{x}}_j<{\underline{x}}_j& \left(a_{1,j}>0\right)\\ x_{gen,j}\in \left[{\underline{x}}_j,{\overline{x}}_j\right]& \left(a_{1,j}=0\right)\\ \frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\overline{x}}_j\right)+{\overline{x}}_j>{\overline{x}}_j& \left(a_{1,j}<0\right)\end{array}\end{array}$$

    情况1.2.2：当$b_1>\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\overline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\underline{x}}_j$时
$$\begin{array}{c}{x}_{cor,j}=\{\begin{array}{ll}\frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\underline{x}}_j\right)+{\underline{x}}_j>{\overline{x}}_j& \left(a_{1,j}>0\right)\\ x_{gen,j}\in \left[{\underline{x}}_j,{\overline{x}}_j\right]& \left(a_{1,j}=0\right)\\ \frac{b_1-\left(\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}{\underline{x}}_j+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}{\overline{x}}_j\right)}{\sum _{j=1,a_{1,j}>0}^N{a}_{1,j}\left(x_{gen,j}-{\underline{x}}_j\right)+\sum _{j=1,a_{1,j}<0}^N{a}_{1,j}\left(x_{gen,j}-{\overline{x}}_j\right)}\left(x_{gen,j}-{\overline{x}}_j\right)+{\overline{x}}_j<{\underline{x}}_j& \left(a_{1,j}<0\right)\end{array}\end{array}$$

    剩余内容参见博客 《线性方程约束下修正方法的可行性分析（演绎推理，二）》