本文内容接续博客 《一种多线性方程约束下的生成随机数修正方法(上)》
实例分析
再次转写成线性方程组,即为
{
x
1
−
209
367
x
5
+
44
367
x
6
=
177
367
x
2
+
1297
1101
x
5
−
1077
734
x
6
=
−
252
367
x
3
−
662
1101
x
5
+
975
734
x
6
=
1323
367
x
4
−
1219
1101
x
5
+
562
367
x
6
=
1043
367
\begin{equation}\tag{50} \begin{cases} x_1 - \dfrac{209}{367}x_5 + \dfrac{44}{367}x_6 = \dfrac{177}{367} \\ x_2 + \dfrac{1297}{1101}x_5 - \dfrac{1077}{734}x_6 = -\dfrac{252}{367} \\ x_3 - \dfrac{662}{1101}x_5 + \dfrac{975}{734}x_6 = \dfrac{1323}{367} \\ x_4 - \dfrac{1219}{1101}x_5 + \dfrac{562}{367}x_6 = \dfrac{1043}{367} \\ \end{cases} \end{equation}
⎩
⎨
⎧x1−367209x5+36744x6=367177x2+11011297x5−7341077x6=−367252x3−1101662x5+734975x6=3671323x4−11011219x5+367562x6=3671043(50)
根据线性方程组理论可知,若随机向量
x
\boldsymbol{x}
x满足式(46)中的等式约束的充要条件是
x
\boldsymbol{x}
x满足式(50)中的等式约束。将式(50)中的
x
1
x_1
x1、
x
2
x_2
x2、
x
3
x_3
x3、
x
4
x_4
x4用
x
5
x_5
x5和
x
6
x_6
x6线性表示,可得
{
x
1
=
177
367
+
209
367
x
5
−
44
367
x
6
x
2
=
−
252
367
−
1297
1101
x
5
+
1077
734
x
6
x
3
=
1323
367
+
662
1101
x
5
−
975
734
x
6
x
4
=
1043
367
+
1219
1101
x
5
−
562
367
x
6
\begin{equation}\tag{51} \begin{cases} x_1 = \dfrac{177}{367} + \dfrac{209}{367}x_5 - \dfrac{44}{367}x_6 \\ x_2 = -\dfrac{252}{367} - \dfrac{1297}{1101}x_5 + \dfrac{1077}{734}x_6 \\ x_3 = \dfrac{1323}{367} + \dfrac{662}{1101}x_5 - \dfrac{975}{734}x_6 \\ x_4 = \dfrac{1043}{367} + \dfrac{1219}{1101}x_5 - \dfrac{562}{367}x_6 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1=367177+367209x5−36744x6x2=−367252−11011297x5+7341077x6x3=3671323+1101662x5−734975x6x4=3671043+11011219x5−367562x6(51)
结合式(46)中的对
x
1
x_1
x1、
x
2
x_2
x2、
x
3
x_3
x3、
x
4
x_4
x4的上下限约束,不难看出要使变量
x
1
x_1
x1、
x
2
x_2
x2、
x
3
x_3
x3、
x
4
x_4
x4满足上下限约束条件,必然有
{
−
8
≤
x
1
=
177
367
+
209
367
x
5
−
44
367
x
6
≤
9
−
15
≤
x
2
=
−
252
367
−
1297
1101
x
5
+
1077
734
x
6
≤
7
−
2
≤
x
3
=
1323
367
+
662
1101
x
5
−
975
734
x
6
≤
11
1
≤
x
4
=
1043
367
+
1219
1101
x
5
−
562
367
x
6
≤
7
\begin{equation}\tag{52} \begin{cases} -8 ≤ x_1 = \dfrac{177}{367} + \dfrac{209}{367}x_5 - \dfrac{44}{367}x_6 ≤ 9 \\ -15 ≤ x_2 = -\dfrac{252}{367} - \dfrac{1297}{1101}x_5 + \dfrac{1077}{734}x_6 ≤ 7 \\ -2 ≤ x_3 = \dfrac{1323}{367} + \dfrac{662}{1101}x_5 - \dfrac{975}{734}x_6 ≤ 11 \\ 1 ≤ x_4 = \dfrac{1043}{367} + \dfrac{1219}{1101}x_5 - \dfrac{562}{367}x_6 ≤ 7 \\ \end{cases} \end{equation}
⎩
⎨
⎧−8≤x1=367177+367209x5−36744x6≤9−15≤x2=−367252−11011297x5+7341077x6≤7−2≤x3=3671323+1101662x5−734975x6≤111≤x4=3671043+11011219x5−367562x6≤7(52)
此时
x
5
x_5
x5和
x
6
x_6
x6满足
{
−
3113
≤
209
x
5
−
44
x
6
≤
3126
−
31518
≤
2594
x
5
+
3231
x
6
≤
16926
−
12342
≤
1324
x
5
−
2925
x
6
≤
16284
−
2028
≤
1219
x
5
−
1686
x
6
≤
4578
\begin{equation}\tag{53} \begin{cases} -3113 ≤ 209x_5 - 44x_6 ≤ 3126 \\ -31518 ≤ 2594x_5 + 3231x_6 ≤ 16926 \\ -12342 ≤ 1324x_5 - 2925x_6 ≤ 16284 \\ -2028 ≤ 1219x_5 - 1686x_6 ≤ 4578 \\ \end{cases} \end{equation}
⎩
⎨
⎧−3113≤209x5−44x6≤3126−31518≤2594x5+3231x6≤16926−12342≤1324x5−2925x6≤16284−2028≤1219x5−1686x6≤4578(53)
根据题设可知,
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen满足
{
−
3113
≤
209
x
5
g
e
n
−
44
x
6
g
e
n
=
−
49
≤
3126
−
31518
≤
2594
x
5
g
e
n
+
3231
x
6
g
e
n
=
−
18479
≤
16926
−
12342
≤
1324
x
5
g
e
n
−
2925
x
6
g
e
n
=
13301
≤
16284
1219
x
5
g
e
n
−
1686
x
6
g
e
n
=
7211
>
4578
\begin{equation}\tag{54} \begin{cases} -3113 ≤ 209x_5^{gen} - 44x_6^{gen} = -49 ≤ 3126 \\ -31518 ≤ 2594x_5^{gen} + 3231x_6^{gen} = -18479 ≤ 16926 \\ -12342 ≤ 1324x_5^{gen} - 2925x_6^{gen} = 13301 ≤ 16284 \\ 1219x_5^{gen} - 1686x_6^{gen} = 7211 > 4578 \\ \end{cases} \end{equation}
⎩
⎨
⎧−3113≤209x5gen−44x6gen=−49≤3126−31518≤2594x5gen+3231x6gen=−18479≤16926−12342≤1324x5gen−2925x6gen=13301≤162841219x5gen−1686x6gen=7211>4578(54)
因此
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen不满足式(53)所示的约束条件,须将
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen修正为
x
5
c
o
r
x_5^{cor}
x5cor和
x
6
c
o
r
x_6^{cor}
x6cor,使其能够满足式(53)所示的约束条件,进而使得修正后的随机向量
x
c
o
r
\boldsymbol{x^{cor}}
xcor满足式(46)中的所有约束条件。
令修正后的
x
5
c
o
r
x_5^{cor}
x5cor和
x
6
c
o
r
x_6^{cor}
x6cor满足
{
b
~
1
c
o
r
=
209
x
5
c
o
r
−
44
x
6
c
o
r
∈
[
−
3113
,
3126
]
b
~
2
c
o
r
=
2594
x
5
c
o
r
+
3231
x
6
c
o
r
∈
[
−
31518
,
16926
]
b
~
3
c
o
r
=
1324
x
5
c
o
r
−
2925
x
6
c
o
r
∈
[
−
12342
,
16284
]
b
~
4
c
o
r
=
1219
x
5
c
o
r
−
1686
x
6
c
o
r
∈
[
−
2028
,
4578
]
\begin{equation}\tag{55} \begin{cases} \tilde{b}_1^{cor} = 209x_5^{cor} - 44x_6^{cor} ∈ \left[ -3113, 3126 \right] \\ \tilde{b}_2^{cor} = 2594x_5^{cor} + 3231x_6^{cor} ∈ \left[ -31518, 16926 \right] \\ \tilde{b}_3^{cor} = 1324x_5^{cor} - 2925x_6^{cor} ∈ \left[ -12342, 16284 \right] \\ \tilde{b}_4^{cor} = 1219x_5^{cor} - 1686x_6^{cor} ∈ \left[ -2028, 4578 \right] \\ \end{cases} \end{equation}
⎩
⎨
⎧b~1cor=209x5cor−44x6cor∈[−3113,3126]b~2cor=2594x5cor+3231x6cor∈[−31518,16926]b~3cor=1324x5cor−2925x6cor∈[−12342,16284]b~4cor=1219x5cor−1686x6cor∈[−2028,4578](55)
可根据式(55)所示的关系将
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen修正为
x
5
c
o
r
x_5^{cor}
x5cor和
x
6
c
o
r
x_6^{cor}
x6cor,再根据以上修正结果将
x
1
g
e
n
x_1^{gen}
x1gen、
x
2
g
e
n
x_2^{gen}
x2gen、
x
3
g
e
n
x_3^{gen}
x3gen、
x
4
g
e
n
x_4^{gen}
x4gen修正为
x
1
c
o
r
x_1^{cor}
x1cor、
x
2
c
o
r
x_2^{cor}
x2cor、
x
3
c
o
r
x_3^{cor}
x3cor、
x
4
c
o
r
x_4^{cor}
x4cor。
根据式(54)可知,
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen不满足式(54)中的前三个约束条件但满足第四个约束条件,则此时可取第四个约束条件并采用博客《一种单线性方程约束下的生成随机数修正方法(结论与应用)》所提修正方法修正
x
5
g
e
n
x_5^{gen}
x5gen和
x
6
g
e
n
x_6^{gen}
x6gen,使得修正后的
x
5
c
o
r
x_5^{cor}
x5cor和
x
6
c
o
r
x_6^{cor}
x6cor满足式(54)中的第一个约束条件,并将其代入式(54)中的第二个约束条件,以判断其同时满足式(54)中所有约束的条件。
由于
1219
x
5
g
e
n
−
1686
x
6
g
e
n
=
7211
>
4578
1219x_5^{gen} - 1686x_6^{gen} = 7211 > 4578
1219x5gen−1686x6gen=7211>4578,则根据博客《一种单线性方程约束下的生成随机数修正方法(结论与应用)》中的式(5)可得
{
x
5
c
o
r
=
1
7120
b
~
4
c
o
r
−
14331
7120
x
6
c
o
r
=
−
7
14240
b
~
4
c
o
r
−
20723
14240
\begin{equation}\tag{56} \begin{cases} x_5^{cor} = \dfrac{1}{7120}\tilde{b}_4^{cor} - \dfrac{14331}{7120} \\ x_6^{cor} = -\dfrac{7}{14240}\tilde{b}_4^{cor} - \dfrac{20723}{14240} \\ \end{cases} \end{equation}
⎩
⎨
⎧x5cor=71201b~4cor−712014331x6cor=−142407b~4cor−1424020723(56)
将其代入式(54)中的前三个约束条件,则有
{
b
~
1
c
o
r
=
846
14240
b
~
4
c
o
r
−
6798266
14240
∈
[
−
3113
,
3126
]
b
~
2
c
o
r
=
−
17429
14240
b
~
4
c
o
r
−
160691241
14240
∈
[
−
31518
,
16926
]
b
~
3
c
o
r
=
23123
14240
b
~
4
c
o
r
+
40216287
14240
∈
[
−
12342
,
16284
]
\begin{equation}\tag{57} \begin{cases} \tilde{b}_1^{cor} = \dfrac{846}{14240}\tilde{b}_4^{cor} - \dfrac{6798266}{14240} ∈ \left[ -3113, 3126 \right] \\ \tilde{b}_2^{cor} = -\dfrac{17429}{14240}\tilde{b}_4^{cor} - \dfrac{160691241}{14240} ∈ \left[ -31518, 16926 \right] \\ \tilde{b}_3^{cor} = \dfrac{23123}{14240}\tilde{b}_4^{cor} + \dfrac{40216287}{14240} ∈ \left[ -12342, 16284 \right] \\ \end{cases} \end{equation}
⎩
⎨
⎧b~1cor=14240846b~4cor−142406798266∈[−3113,3126]b~2cor=−1424017429b~4cor−14240160691241∈[−31518,16926]b~3cor=1424023123b~4cor+1424040216287∈[−12342,16284](57)
联立并约简可得
b
~
4
c
o
r
∈
[
−
215966367
23123
,
191667873
23123
]
\begin{equation}\tag{58} \tilde{b}_4^{cor} ∈ \left[ -\dfrac{215966367}{23123}, \dfrac{191667873}{23123} \right] \end{equation}
b~4cor∈[−23123215966367,23123191667873](58)
结合式(55)中对
b
~
4
c
o
r
\tilde{b}_4^{cor}
b~4cor的约束可知,若使修正后的
x
5
c
o
r
x_5^{cor}
x5cor和
x
6
c
o
r
x_6^{cor}
x6cor满足式(54)中的所有约束,则
b
~
4
c
o
r
\tilde{b}_4^{cor}
b~4cor满足
b
~
4
c
o
r
∈
[
−
2028
,
4578
]
\begin{equation}\tag{59} \tilde{b}_4^{cor} ∈ \left[ -2028, 4578 \right] \end{equation}
b~4cor∈[−2028,4578](59)
修正方式3.1:令式(55)中的
b
~
4
c
o
r
=
−
2028
\tilde{b}_4^{cor} = -2028
b~4cor=−2028,则根据式(56)可得
{
x
5
c
o
r
=
−
16359
7120
x
6
c
o
r
=
−
6527
14240
\begin{equation}\tag{60} \begin{cases} x_5^{cor} = -\dfrac{16359}{7120} \\ x_6^{cor} = -\dfrac{6527}{14240} \\ \end{cases} \end{equation}
⎩
⎨
⎧x5cor=−712016359x6cor=−142406527(60)
将其代入式(51),则
{
x
1
c
o
r
=
−
5491
7120
x
2
c
o
r
=
7675
5696
x
3
c
o
r
=
80663
28480
x
4
c
o
r
=
1
\begin{equation}\tag{61} \begin{cases} x_1^{cor} = -\dfrac{5491}{7120} \\ x_2^{cor} = \dfrac{7675}{5696} \\ x_3^{cor} = \dfrac{80663}{28480} \\ x_4^{cor} = 1 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−71205491x2cor=56967675x3cor=2848080663x4cor=1(61)
经如式(62)所示的验算可知,所得
x
c
o
r
=
[
−
0.7712
,
1.3474
,
2.8323
,
1
,
−
2.2976
,
−
0.4584
]
T
\boldsymbol{x^{cor}} = \left[ -0.7712, 1.3474, 2.8323, 1, -2.2976, -0.4584 \right]^T
xcor=[−0.7712,1.3474,2.8323,1,−2.2976,−0.4584]T满足
{
x
1
c
o
r
=
−
5491
7120
∈
[
−
8
,
9
]
x
2
c
o
r
=
7675
5696
∈
[
−
15
,
7
]
x
3
c
o
r
=
80663
28480
∈
[
−
2
,
11
]
x
4
c
o
r
=
1
∈
[
1
,
7
]
x
5
c
o
r
=
−
16359
7120
∈
[
−
3
,
8
]
x
6
c
o
r
=
−
6527
14240
∈
[
−
10
,
2
]
4
x
1
c
o
r
+
3
x
2
c
o
r
+
5
x
3
c
o
r
−
7
x
4
c
o
r
+
6
x
5
c
o
r
−
8
x
6
c
o
r
=
−
2
−
7
x
1
c
o
r
−
4
x
2
c
o
r
+
8
x
3
c
o
r
−
5
x
4
c
o
r
+
8
x
6
c
o
r
=
14
10
x
1
c
o
r
+
3
x
2
c
o
r
−
3
x
3
c
o
r
+
6
x
4
c
o
r
−
7
x
5
c
o
r
+
2
x
6
c
o
r
=
9
2
x
1
c
o
r
−
3
x
2
c
o
r
−
7
x
3
c
o
r
+
5
x
4
c
o
r
−
6
x
5
c
o
r
+
3
x
6
c
o
r
=
−
8
\begin{equation}\tag{62} \begin{cases} x_1^{cor} = -\dfrac{5491}{7120} ∈ \left[ -8, 9 \right] \\ x_2^{cor} = \dfrac{7675}{5696} ∈ \left[ -15, 7 \right] \\ x_3^{cor} = \dfrac{80663}{28480}∈ \left[ -2, 11 \right] \\ x_4^{cor} = 1 ∈ \left[ 1, 7 \right] \\ x_5^{cor} = -\dfrac{16359}{7120} ∈ \left[ -3, 8 \right] \\ x_6^{cor} = -\dfrac{6527}{14240} ∈ \left[ -10, 2 \right] \\ 4x_1^{cor} + 3x_2^{cor} + 5x_3^{cor} - 7x_4^{cor} + 6x_5^{cor} - 8x_6^{cor} = -2 \\ -7x_1^{cor} - 4x_2^{cor} + 8x_3^{cor} - 5x_4^{cor} + 8x_6^{cor} = 14 \\ 10x_1^{cor} + 3x_2^{cor} - 3x_3^{cor} + 6x_4^{cor} - 7x_5^{cor} + 2x_6^{cor} = 9 \\ 2x_1^{cor} - 3x_2^{cor} - 7x_3^{cor} + 5x_4^{cor} - 6x_5^{cor} + 3x_6^{cor} = -8 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−71205491∈[−8,9]x2cor=56967675∈[−15,7]x3cor=2848080663∈[−2,11]x4cor=1∈[1,7]x5cor=−712016359∈[−3,8]x6cor=−142406527∈[−10,2]4x1cor+3x2cor+5x3cor−7x4cor+6x5cor−8x6cor=−2−7x1cor−4x2cor+8x3cor−5x4cor+8x6cor=1410x1cor+3x2cor−3x3cor+6x4cor−7x5cor+2x6cor=92x1cor−3x2cor−7x3cor+5x4cor−6x5cor+3x6cor=−8(62)
故而满足式(46)中的所有约束条件,
x
c
o
r
\boldsymbol{x^{cor}}
xcor可以作为式(46)所示问题的其中一个可行解。
修正方式3.2:令式(55)中的
b
~
4
c
o
r
=
2215
\tilde{b}_4^{cor} = 2215
b~4cor=2215,则根据式(56)可得
{
x
5
c
o
r
=
−
3029
1780
x
6
c
o
r
=
−
9057
3560
\begin{equation}\tag{63} \begin{cases} x_5^{cor} = -\dfrac{3029}{1780} \\ x_6^{cor} = -\dfrac{9057}{3560} \\ \end{cases} \end{equation}
⎩
⎨
⎧x5cor=−17803029x6cor=−35609057(63)
将其代入式(51),则
{
x
1
c
o
r
=
−
118747
653260
x
2
c
o
r
=
−
3786287
1567824
x
3
c
o
r
=
46730213
7839120
x
4
c
o
r
=
5344
1101
\begin{equation}\tag{64} \begin{cases} x_1^{cor} = -\dfrac{118747}{653260} \\ x_2^{cor} = -\dfrac{3786287}{1567824} \\ x_3^{cor} = \dfrac{46730213}{7839120} \\ x_4^{cor} = \dfrac{5344}{1101} \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−653260118747x2cor=−15678243786287x3cor=783912046730213x4cor=11015344(64)
经如式(65)所示的验算可知,所得
x
c
o
r
=
[
−
0.1818
,
−
2.4150
,
5.9612
,
4.8538
,
−
1.7017
,
−
2.5441
]
T
\boldsymbol{x^{cor}} = \left[ -0.1818, -2.4150, 5.9612, 4.8538, -1.7017, -2.5441 \right]^T
xcor=[−0.1818,−2.4150,5.9612,4.8538,−1.7017,−2.5441]T满足
{
x
1
c
o
r
=
−
118747
653260
∈
[
−
8
,
9
]
x
2
c
o
r
=
−
3786287
1567824
∈
[
−
15
,
7
]
x
3
c
o
r
=
46730213
7839120
∈
[
−
2
,
4
]
x
4
c
o
r
=
5344
1101
∈
[
0
,
5
]
x
5
c
o
r
=
−
3029
1780
∈
[
−
3
,
8
]
x
6
c
o
r
=
−
9057
3560
∈
[
−
10
,
2
]
4
x
1
c
o
r
+
3
x
2
c
o
r
+
5
x
3
c
o
r
−
7
x
4
c
o
r
+
6
x
5
c
o
r
−
8
x
6
c
o
r
=
−
2
−
7
x
1
c
o
r
−
4
x
2
c
o
r
+
8
x
3
c
o
r
−
5
x
4
c
o
r
+
8
x
6
c
o
r
=
14
10
x
1
c
o
r
+
3
x
2
c
o
r
−
3
x
3
c
o
r
+
6
x
4
c
o
r
−
7
x
5
c
o
r
+
2
x
6
c
o
r
=
9
2
x
1
c
o
r
−
3
x
2
c
o
r
−
7
x
3
c
o
r
+
5
x
4
c
o
r
−
6
x
5
c
o
r
+
3
x
6
c
o
r
=
−
8
\begin{equation}\tag{65} \begin{cases} x_1^{cor} = -\dfrac{118747}{653260} ∈ \left[ -8, 9 \right] \\ x_2^{cor} = -\dfrac{3786287}{1567824} ∈ \left[ -15, 7 \right] \\ x_3^{cor} = \dfrac{46730213}{7839120} ∈ \left[ -2, 4 \right] \\ x_4^{cor} = \dfrac{5344}{1101} ∈ \left[ 0, 5 \right] \\ x_5^{cor} = -\dfrac{3029}{1780} ∈ \left[ -3, 8 \right] \\ x_6^{cor} = -\dfrac{9057}{3560} ∈ \left[ -10, 2 \right] \\ 4x_1^{cor} + 3x_2^{cor} + 5x_3^{cor} - 7x_4^{cor} + 6x_5^{cor} - 8x_6^{cor} = -2 \\ -7x_1^{cor} - 4x_2^{cor} + 8x_3^{cor} - 5x_4^{cor} + 8x_6^{cor} = 14 \\ 10x_1^{cor} + 3x_2^{cor} - 3x_3^{cor} + 6x_4^{cor} - 7x_5^{cor} + 2x_6^{cor} = 9 \\ 2x_1^{cor} - 3x_2^{cor} - 7x_3^{cor} + 5x_4^{cor} - 6x_5^{cor} + 3x_6^{cor} = -8 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−653260118747∈[−8,9]x2cor=−15678243786287∈[−15,7]x3cor=783912046730213∈[−2,4]x4cor=11015344∈[0,5]x5cor=−17803029∈[−3,8]x6cor=−35609057∈[−10,2]4x1cor+3x2cor+5x3cor−7x4cor+6x5cor−8x6cor=−2−7x1cor−4x2cor+8x3cor−5x4cor+8x6cor=1410x1cor+3x2cor−3x3cor+6x4cor−7x5cor+2x6cor=92x1cor−3x2cor−7x3cor+5x4cor−6x5cor+3x6cor=−8(65)
故而满足式(46)中的所有约束条件,
x
c
o
r
\boldsymbol{x^{cor}}
xcor可以作为式(46)所示问题的其中一个可行解。
修正方式3.3:令式(55)中的
b
~
4
c
o
r
=
4396
\tilde{b}_4^{cor} = 4396
b~4cor=4396,则根据式(56)可得
{
x
5
c
o
r
=
−
1987
1424
x
6
c
o
r
=
−
10299
2848
\begin{equation}\tag{66} \begin{cases} x_5^{cor} = -\dfrac{1987}{1424} \\ x_6^{cor} = -\dfrac{10299}{2848} \\ \end{cases} \end{equation}
⎩
⎨
⎧x5cor=−14241987x6cor=−284810299(66)
将其代入式(51),则
{
x
1
c
o
r
=
−
63343
522608
x
2
c
o
r
=
−
27273689
6271296
x
3
c
o
r
=
47470423
6271296
x
4
c
o
r
=
7525
1101
\begin{equation}\tag{67} \begin{cases} x_1^{cor} = -\dfrac{63343}{522608} \\ x_2^{cor} = -\dfrac{27273689}{6271296} \\ x_3^{cor} = \dfrac{47470423}{6271296} \\ x_4^{cor} = \dfrac{7525}{1101} \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−52260863343x2cor=−627129627273689x3cor=627129647470423x4cor=11017525(67)
经如式(68)所示的验算可知,所得
x
c
o
r
=
[
0.1212
,
−
4.3490
,
7.5695
,
6.8347
,
−
1.3954
,
−
3.6162
]
T
\boldsymbol{x^{cor}} = \left[ 0.1212, -4.3490, 7.5695, 6.8347, -1.3954, -3.6162 \right]^T
xcor=[0.1212,−4.3490,7.5695,6.8347,−1.3954,−3.6162]T满足
{
x
1
c
o
r
=
−
63343
522608
∈
[
−
8
,
9
]
x
2
c
o
r
=
−
27273689
6271296
∈
[
−
15
,
7
]
x
3
c
o
r
=
47470423
6271296
∈
[
−
2
,
4
]
x
4
c
o
r
=
7525
1101
∈
[
0
,
5
]
x
5
c
o
r
=
−
1987
1424
∈
[
−
3
,
8
]
x
6
c
o
r
=
−
10299
2848
∈
[
−
10
,
2
]
4
x
1
c
o
r
+
3
x
2
c
o
r
+
5
x
3
c
o
r
−
7
x
4
c
o
r
+
6
x
5
c
o
r
−
8
x
6
c
o
r
=
−
2
−
7
x
1
c
o
r
−
4
x
2
c
o
r
+
8
x
3
c
o
r
−
5
x
4
c
o
r
+
8
x
6
c
o
r
=
14
10
x
1
c
o
r
+
3
x
2
c
o
r
−
3
x
3
c
o
r
+
6
x
4
c
o
r
−
7
x
5
c
o
r
+
2
x
6
c
o
r
=
9
2
x
1
c
o
r
−
3
x
2
c
o
r
−
7
x
3
c
o
r
+
5
x
4
c
o
r
−
6
x
5
c
o
r
+
3
x
6
c
o
r
=
−
8
\begin{equation}\tag{68} \begin{cases} x_1^{cor} = -\dfrac{63343}{522608}∈ \left[ -8, 9 \right] \\ x_2^{cor} = -\dfrac{27273689}{6271296} ∈ \left[ -15, 7 \right] \\ x_3^{cor} = \dfrac{47470423}{6271296} ∈ \left[ -2, 4 \right] \\ x_4^{cor} = \dfrac{7525}{1101} ∈ \left[ 0, 5 \right] \\ x_5^{cor} = -\dfrac{1987}{1424} ∈ \left[ -3, 8 \right] \\ x_6^{cor} = -\dfrac{10299}{2848} ∈ \left[ -10, 2 \right] \\ 4x_1^{cor} + 3x_2^{cor} + 5x_3^{cor} - 7x_4^{cor} + 6x_5^{cor} - 8x_6^{cor} = -2 \\ -7x_1^{cor} - 4x_2^{cor} + 8x_3^{cor} - 5x_4^{cor} + 8x_6^{cor} = 14 \\ 10x_1^{cor} + 3x_2^{cor} - 3x_3^{cor} + 6x_4^{cor} - 7x_5^{cor} + 2x_6^{cor} = 9 \\ 2x_1^{cor} - 3x_2^{cor} - 7x_3^{cor} + 5x_4^{cor} - 6x_5^{cor} + 3x_6^{cor} = -8 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−52260863343∈[−8,9]x2cor=−627129627273689∈[−15,7]x3cor=627129647470423∈[−2,4]x4cor=11017525∈[0,5]x5cor=−14241987∈[−3,8]x6cor=−284810299∈[−10,2]4x1cor+3x2cor+5x3cor−7x4cor+6x5cor−8x6cor=−2−7x1cor−4x2cor+8x3cor−5x4cor+8x6cor=1410x1cor+3x2cor−3x3cor+6x4cor−7x5cor+2x6cor=92x1cor−3x2cor−7x3cor+5x4cor−6x5cor+3x6cor=−8(68)
故而满足式(46)中的所有约束条件,
x
c
o
r
\boldsymbol{x^{cor}}
xcor可以作为式(46)所示问题的其中一个可行解。
修正方式3.4:令式(55)中的
b
~
4
c
o
r
=
4578
\tilde{b}_4^{cor} = 4578
b~4cor=4578,则根据式(56)可得
{
x
5
c
o
r
=
−
9753
7120
x
6
c
o
r
=
−
52769
14240
\begin{equation}\tag{69} \begin{cases} x_5^{cor} = -\dfrac{9753}{7120} \\ x_6^{cor} = -\dfrac{52769}{14240} \\ \end{cases} \end{equation}
⎩
⎨
⎧x5cor=−71209753x6cor=−1424052769(69)
将其代入式(51),则
{
x
1
c
o
r
=
−
1043
7120
x
2
c
o
r
=
−
25691
5696
x
3
c
o
r
=
219401
28480
x
4
c
o
r
=
7
\begin{equation}\tag{70} \begin{cases} x_1^{cor} = -\dfrac{1043}{7120} \\ x_2^{cor} = -\dfrac{25691}{5696} \\ x_3^{cor} = \dfrac{219401}{28480} \\ x_4^{cor} = 7 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−71201043x2cor=−569625691x3cor=28480219401x4cor=7(70)
经如式(71)所示的验算可知,所得
x
c
o
r
=
[
0.1465
,
−
4.5104
,
7.7037
,
7
,
−
1.3698
,
−
3.7057
]
T
\boldsymbol{x^{cor}} = \left[ 0.1465, -4.5104, 7.7037, 7, -1.3698, -3.7057 \right]^T
xcor=[0.1465,−4.5104,7.7037,7,−1.3698,−3.7057]T满足
{
x
1
c
o
r
=
−
1043
7120
∈
[
−
8
,
9
]
x
2
c
o
r
=
−
25691
5696
∈
[
−
15
,
7
]
x
3
c
o
r
=
219401
28480
∈
[
−
2
,
4
]
x
4
c
o
r
=
7
∈
[
0
,
5
]
x
5
c
o
r
=
−
9753
7120
∈
[
−
3
,
8
]
x
6
c
o
r
=
−
52769
14240
∈
[
−
10
,
2
]
4
x
1
c
o
r
+
3
x
2
c
o
r
+
5
x
3
c
o
r
−
7
x
4
c
o
r
+
6
x
5
c
o
r
−
8
x
6
c
o
r
=
−
2
−
7
x
1
c
o
r
−
4
x
2
c
o
r
+
8
x
3
c
o
r
−
5
x
4
c
o
r
+
8
x
6
c
o
r
=
14
10
x
1
c
o
r
+
3
x
2
c
o
r
−
3
x
3
c
o
r
+
6
x
4
c
o
r
−
7
x
5
c
o
r
+
2
x
6
c
o
r
=
9
2
x
1
c
o
r
−
3
x
2
c
o
r
−
7
x
3
c
o
r
+
5
x
4
c
o
r
−
6
x
5
c
o
r
+
3
x
6
c
o
r
=
−
8
\begin{equation}\tag{71} \begin{cases} x_1^{cor} = -\dfrac{1043}{7120} ∈ \left[ -8, 9 \right] \\ x_2^{cor} = -\dfrac{25691}{5696} ∈ \left[ -15, 7 \right] \\ x_3^{cor} = \dfrac{219401}{28480} ∈ \left[ -2, 4 \right] \\ x_4^{cor} = 7 ∈ \left[ 0, 5 \right] \\ x_5^{cor} = -\dfrac{9753}{7120} ∈ \left[ -3, 8 \right] \\ x_6^{cor} = -\dfrac{52769}{14240} ∈ \left[ -10, 2 \right] \\ 4x_1^{cor} + 3x_2^{cor} + 5x_3^{cor} - 7x_4^{cor} + 6x_5^{cor} - 8x_6^{cor} = -2 \\ -7x_1^{cor} - 4x_2^{cor} + 8x_3^{cor} - 5x_4^{cor} + 8x_6^{cor} = 14 \\ 10x_1^{cor} + 3x_2^{cor} - 3x_3^{cor} + 6x_4^{cor} - 7x_5^{cor} + 2x_6^{cor} = 9 \\ 2x_1^{cor} - 3x_2^{cor} - 7x_3^{cor} + 5x_4^{cor} - 6x_5^{cor} + 3x_6^{cor} = -8 \\ \end{cases} \end{equation}
⎩
⎨
⎧x1cor=−71201043∈[−8,9]x2cor=−569625691∈[−15,7]x3cor=28480219401∈[−2,4]x4cor=7∈[0,5]x5cor=−71209753∈[−3,8]x6cor=−1424052769∈[−10,2]4x1cor+3x2cor+5x3cor−7x4cor+6x5cor−8x6cor=−2−7x1cor−4x2cor+8x3cor−5x4cor+8x6cor=1410x1cor+3x2cor−3x3cor+6x4cor−7x5cor+2x6cor=92x1cor−3x2cor−7x3cor+5x4cor−6x5cor+3x6cor=−8(71)
故而满足式(46)中的所有约束条件,
x
c
o
r
\boldsymbol{x^{cor}}
xcor可以作为式(46)所示问题的其中一个可行解。
研究目标
(1) 根据以上具体实例的修正思路,探究满足线性约束的生成随机数的一般性修正方法;
(2) 探究满足某些非线性约束的生成随机数修正方法;
(3) 将生成随机数修正方法应用到基于启发式算法的优化问题求解中,使得启发式算法能够始终在优化问题的可行域中搜寻问题的解,从而加快启发式优化算法的收敛速度。