本文认为修改公式的标准共轭梯度(CG)技术,由李和福岛计划。一个新的标量参数
无约束最优化CG技术的计划。血统和全局收敛的条件属性建立如下使用强大的沃尔夫条件。我们的数值实验表明,新提出的算法更加稳定和经济相比,一些著名的标准CG方法。
1。介绍
共轭梯度(CG)策略包括一类非线性优化算法,需要低内存和强大的本地和全局收敛性属性(1 ,2 ]。通常,CG方法是为了解决大规模扩展非线性优化问题:
在理解中定义的函数形式
是光滑的非线性函数。重复的公式形式
这个公式是最重要的组成部分
步长,搜索方向
由
而
表示
和
表示积极的标量。的步长
有时选择满足绑定线搜索条件(3 ]。在这些搜索方向条件下,强劲的沃尔夫线搜索条件有时概述如下:
和
。 有许多不同的公式为共轭系数在接下来的来源,例如,Hestenes施蒂费尔,HS (4 ];弗莱彻和Reeves, FR (5 ];波兰人和Ribiere公关(6 ];共轭后裔,CD (7 ];李和福岛,低频(1 ];刘和故事,LS (8 为标量参数),对应于不同的选择
。
在本文的这部分,我们提出了一个新版本的参数
依靠修改后的高炉煤气的方法提出的李和福岛(1 ]。在蓄热方法,矩阵
更新以下公式(9 ]:
在哪里
。 此外,列出的正常sec关系符合后续公式:
研究者李和福岛提出适当的修改全局和超线性收敛的蓄热技术,即使在不需要凸目标函数。随后修改sec方程符合后续列出公式如下:
在哪里
和
> 0;
> 0,
是一个参数定义为
具体地说,把价值c 是常数,是大于零的。
有三种不同的情况下
: 案例1:如果
在这种情况下我们有负值的明确的矩阵的问题,所以李和福岛
公式如(9 )和开发相应的高炉煤气的公式如下:
此外,的形式
在(10 ),当使用max的值(0;在这种情况下零),不是选择。通过这个公式,研究人员证明了对称矩阵是正定(修改10 ]。 例2:如果
, 在这种情况下,我们可以肯定说,bfg更新矩阵是对称正定当应用这个公式(换句话说,当内应用不平等
在这个公式
, max = 0) (11 ]。
如果我们使用
刘和故事(LS),我们使用任何标量
; 然后,(3 )成为
在哪里
当任何积极的价值
大于一个,新的并行搜索方向吗
在方程(12 )是牛顿方向。因此,牛顿方向
因此,
使用方程(8 ),新的标量
就变成了
用方程(9 )和(10 ),以
(因为我们使用强大的沃尔夫行搜索在方程(10 ),收益率
。 因此,在新的搜索方向是新的标量
因此,我们得出结论从方程(17 ),新参数
是最好的,因为它是最新发现的价值y 时,我们也发现不同形式改变的价值c 我们会注意到部分的数值结果。
3所示。新定理(充分下降方向)
如果我们假定线搜索满足条件(4 )和(5 ),然后生成新的搜索方向从方程(12 )和(17 )可以充分下降方向。
证明。 从方程(12 )和(17 ),我们得到
通过使用鲍威尔重启方程(即
),
如果
> 0,第二不平等是正确的:
使用强大的沃尔夫行搜索条件(5)产量
这意味着后者方程
因此,我们要求是完整的。
3.1。概述了新CG-Algorithms
步骤1:选择初始点
, 并选择一些积极的值
和
。 然后,设置
并设置
。 第二步:测试停止标准。如果满意,然后停止;否则,继续。 步骤3:确定
沃尔夫条件,定义在方程(4 )和(5 )。 第四步:计算第二个迭代点
从方程(2 )。 第五步:计算标量参数
从方程(17 )。 第六步:计算新的搜索方向,即
第七步:测试鲍威尔重新启动判据,即,如果
, 然后重新启动新的搜索方向
。 第八步:设置下一次迭代k =k + 1,然后转到步骤2。
4所示。收敛性分析的新算法
在接下来的部分,我们倾向于讨论新算法的收敛性分析房地产。首先,我们提供一个假设收敛分析房地产的新算法。然后,我们的研究中提供了另一种著名的引理需要收敛分析性质。最后,我们有一种倾向,设立新的定理在证明,单位面积与新算法的收敛性分析。
假设。 (我)水平集
是有界的,存在一个常数z > 0,如(12 ]
(2) 在邻里N的年代,f 是连续可微的,其梯度李普希茨是连续的,也就是说,存在一个常数l > 0,等
假设(i)和(ii)f ,我们能够推断出存在
> 0,如
引理。 如果我们假设3 ,13 ] (1) 假设成立。 (2) 搜索方向
在标准的CG方法是一个下降方向。 (3) 优化步骤
通过方程计算(4 )和(5 )。 (4) 满足收敛条件,即。,如果
然后,
如果我们假设(1) 假设成立。 (2) 新的搜索方向
由方程(12 )和(17 )是一个下降方向。 (3) 最优的步骤
通过方程计算(4 )和(5 )。 (4) 目标函数f 一致凸的;然后,
证明。 考虑到新的方向在方程(12 )和方程的参数(17 )满足下一个绝对值条件:
好参数
此外,通过结合的结果,我们得到
我们得到了所需的证据。我们把相似的点前面的假设,但也有一些变化公式。
6。新定理(通用功能)
如果我们假设(1) 假设成立。 (2) 新的搜索方向
由方程(12 )和(17 )是一个下降方向。 (3) 最优的步骤
通过方程计算(4 )和(5 )。 (4) 目标函数f 通用功能;然后,
证明。 使用相同的风格的定理证明的事实的差异函数算法的一般函数,
然后,我们得到
因此,关于新定理的证明该算法的收敛性分析完成。
7所示。数值实验
在本节中,我们报道一些数值试验进行一组(60)无约束优化测试问题分析的效率
。 这些测试的细节问题,可以找到与给定的初始点,在14 ,15 ]。我们处理这些(60)通过添加1000为每个测试函数n 到达最大数量n 等于10000。在我们的实验中使用的终止准则
。
在比较下,我们采用以下算法:(我) LS:沃尔夫行搜索 (2) CD:沃尔夫行搜索 (3) 海关:沃尔夫行搜索 (iv) 公关:沃尔夫行搜索 (v) 新算法,利用方程(17 )和标量c = 0.1 (vi) 新算法,利用方程(17 )和标量c = 0.001
在表中1 和2 ,我们数值比较提出新的CG算法与其他著名的CG算法使用已知的比较工具来验证其性能等算法如下:陈列=的总数计算迭代的迭代 NOFG =函数和梯度计算的总数 时间=处理器执行所需的总CPU时间重心算法和达到所需的最小值函数最小化
概率。
新算法(C = 0.1)陈列/ NOFG /时间
新算法(C = 0.001)陈列/ NOFG /时间
LS 有陈列/ NOFG /时间
CD 有陈列/ NOFG /时间
1
99/252/0.14
101/263/0.14
323/7045/2.31
414/9608/2.64
2
408/884/2.02
405/880/1.97
412/888/1.98
395/878/2.08
3
853/2183/1.06
830/2171/0.98
823/2159/1.02
824/2154/0.97
4
123/308/0.14
123/308/0.15
114/300/0.17
121/314/0.15
5
100/388/0.15
100/388/0.20
274/5586/1.14
167/1918/0.77
6
585/995/1.56
585/995/1.54
6576/19689/4.54
6848/20496/5.70
7
30/80/0.15
30/80/0.13
40/100/0.21
40/100/0.16
8
1032/2705/1.58
1020/2668/1.42
1033/2768/1.49
1026/2972/1.52
9
2388/4881/2.40
2503/5048/3.06
3520/7182/7.49
3676/7485/8.83
10
477/871/2.66
478/858/2.62
18795/503441/7.60
19346/518146/9.14
11
182/423/1.03
174/418/1.09
8669/270951/5.83
8252/257171/6.83
12
113/302/0.16
113/302/0.20
318/3881/1.69
279/2792/5.80
13
80/226/0.11
78/222/0.08
73/241/0.10
73/241/0.09
14
61/131/0.50
61/131/0.56
124/2140/7.12
160/2833/3.54
15
460/991/0.62
452/969/0.59
456/953/0.61
467/997/0.64
16
66/132/0.03
66/132/0.04
67/134/0.14
68/136/0.06
17
70/160/0.11
70/160/0.06
69/158/0.10
71/162/0.08
18
753/1577/0.79
791/1668/0.85
691/1477/0.74
782/1687/0.82
19
74/158/0.33
74/158/0.42
123/1905/2.31
70/150/0.30
20.
110/349/0.41
109/349/0.37
114/334/0.43
111/329/0.41
21
806/3224/1.60
600/1665/1.26
429/1895/4.93
561/5284/5.69
22
72/275/0.51
72/285/0.49
84/366/0.70
2084/2322/9.54
23
4470/9572/1.33
5033/10748/2.42
20010/168275/4.95
19334/128670/5.02
24
62/201/0.42
62/201/0.40
546/15413/3.73
495/13595/3.48
25
459/1091/0.60
521/1192/0.78
494/1101/0.65
531/1154/0.76
26
56/153/0.08
65/373/0.13
56/155/0.06
66/328/0.12
27
85/203/0.11
85/203/0.12
80/190/0.07
80/190/0.09
28
534/1139/0.65
514/1073/0.68
493/1055/0.63
492/1075/0.61
29日
540/1274/0.62
537/1267/0.58
941/2237/1.03
943/2239/1.11
30.
576/1440/0.96
591/1504/1.05
19026/581537/2.57
18144/543276/2.03
31日
113/236/0.16
113/236/0.11
118/245/0.11
123/253/0.11
32
813/2181/3.28
735/1997/3.03
18140/560188/8.67
18320/564064/8.04
33
98/268/0.18
98/268/0.17
359/5355/2.21
387/6200/3.06
34
346/766/0.48
348/766/0.50
356/786/0.44
332/732/0.44
35
7635/12820/8.59
7554/12715/8.73
7463/12576/8.29
7183/12051/8.37
36
280/978/0.43
280/978/0.35
264/935/0.45
273/978/0.39
37
217/534/0.31
217/534/0.27
221/551/0.31
219/546/0.31
38
121/287/0.14
120/285/0.13
814/20969/8.83
786/20117/7.46
39
150/329/0.64
153/328/0.65
141/305/0.59
137/297/0.58
40
107/217/0.65
107/217/0.68
144/1165/3.67
124/476/1.35
41
120/330/0.23
120/330/0.25
100/290/0.17
100/290/0.18
42
3832/9998/2.51
3373/8928/7.18
3404/9046/8.47
3182/8618/7.82
43
40/80/0.11
40/80/0.08
40/80/0.12
40/80/0.12
44
50/110/0.10
50/110/0.05
50/110/0.09
50/110/0.05
45
43/184/0.08
43/184/0.11
3903/129109/7.10
7265/242193/4.78
46
427/1323/2.56
427/1320/2.60
413/1307/2.53
428/1252/2.44
47
64/249/0.09
64/249/0.06
118/447/0.22
114/425/0.19
48
308/798/1.39
293/772/1.28
773/7349/0.72
1731/23050/5.88
49
22/89/0.11
22/89/0.13
316/8561/9.47
298/8022/8.33
50
20/50/0.03
20/50/0.02
20/50/0.00
20/50/0.04
51
92/1755/2.25
42/136/0.17
361/9646/9.02
315/8202/8.61
52
107/418/0.16
107/418/0.15
107/418/0.16
107/418/0.16
53
6199/52426/6.75
6199/52426/6.78
8071/53291/6.20
8071/53291/9.30
54
51/151/0.20
51/151/0.19
51/151/0.21
51/151/0.18
55
60/140/0.19
60/140/0.19
60/140/0.18
60/140/0.20
56
70/140/0.20
70/140/0.17
70/140/0.19
70/140/0.16
57
79/158/0.22
79/158/0.22
79/158/0.24
74/148/0.22
58
143/570/0.25
143/570/0.25
143/570/0.23
143/570/0.21
59
188/498/0.29
172/453/0.27
176/449/0.27
177/461/0.22
60
83/236/0.13
83/236/0.11
985/27409/9.27
917/25010/8.38
总
37602/124886/55.54
37426/121643/59.2
131933/2455353/154.98
136949/1242169/166.56
概率。
新算法(C = 0.1)陈列/ NOFG /时间
新算法(C = 0.001)陈列/ NOFG /时间
海关 有陈列/ NOFG /时间
公关 有陈列/ NOFG /时间
1
99/252/0.14
101/263/0.14
5902/173484/0.72
11798/270815/6.07
2
408/884/2.02
405/880/1.97
362/637/1.83
416/720/2.07
3
853/2183/1.06
830/2171/0.98
789/1817/0.86
989/1979/1.11
4
123/308/0.14
123/308/0.15
141/281/0.19
254/430/0.26
5
100/388/0.15
100/388/0.20
653/17073/3.79
410/7948/2.73
6
585/995/1.56
585/995/1.54
9881/16243/3.07
20010/22091/4.04
7
30/80/0.15
30/80/0.13
40/90/0.18
40/90/0.18
8
1032/2705/1.58
1020/2668/1.42
997/2279/1.28
8780/10109/1.57
9
2388/4881/2.40
2503/5048/3.06
4658/7644/2.48
14945/16013/9.92
10
477/871/2.66
478/858/2.62
20010/98744/4.70
20010/292528/9.39
11
182/423/1.03
174/418/1.09
14053/39805/1.28
15593/489487/8.58
12
113/302/0.16
113/302/0.20
428/6427/2.19
609/11752/4.81
13
80/226/0.11
78/222/0.08
113/234/0.08
291/509/0.26
14
61/131/0.50
61/131/0.56
906/23792/1.66
318/6032/3.76
15
460/991/0.62
452/969/0.59
636/1006/0.73
964/1479/1.25
16
66/132/0.03
66/132/0.04
60/120/0.03
1043/1116/0.38
17
70/160/0.11
70/160/0.06
207/339/0.19
110/230/0.10
18
753/1577/0.79
791/1668/0.85
821/1545/0.78
3732/4630/3.99
19
74/158/0.33
74/158/0.42
108/1352/2.44
303/7194/2.66
20.
110/349/0.41
109/349/0.37
135/321/0.42
154/339/0.45
21
806/3224/1.60
600/1665/1.26
875/14122/5.89
929/10442/9.02
22
72/275/0.51
72/285/0.49
2104/2442/3.61
161/440/0.92
23
4470/9572/1.33
5033/10748/2.42
18912/38658/8.87
20010/25808/6.95
24
62/201/0.42
62/201/0.40
1853/6983/3.39
2527/76854/9.01
25
459/1091/0.60
521/1192/0.78
304/606/0.40
1199/1793/1.55
26
56/153/0.08
65/373/0.13
128/1103/0.37
2697/14700/7.09
27
85/203/0.11
85/203/0.12
91/193/0.13
132/264/0.16
28
534/1139/0.65
514/1073/0.68
288/558/0.35
556/925/0.60
29日
540/1274/0.62
537/1267/0.58
852/1783/0.98
1014/2180/1.14
30.
576/1440/0.96
591/1504/1.05
20010/98171/6.44
20010/317766/8.02
31日
113/236/0.16
113/236/0.11
79/168/0.10
147/287/0.13
32
813/2181/3.28
735/1997/3.03
20010/91480/9.47
20010/179051/5.16
33
98/268/0.18
98/268/0.17
631/11069/5.15
837/17999/8.99
34
346/766/0.48
348/766/0.50
716/1148/0.90
744/1213/0.93
35
7635/12820/8.59
7554/12715/8.73
8375/13146/9.98
8539/12513/15.46
36
280/978/0.43
280/978/0.35
330/695/0.37
401/868/0.47
37
217/534/0.31
217/534/0.27
610/6778/2.40
820/11112/5.50
38
121/287/0.14
120/285/0.13
1565/42467/4.89
1624/45918/6.43
39
150/329/0.64
153/328/0.65
174/290/0.54
193/323/0.69
40
107/217/0.65
107/217/0.68
253/426/1.03
4281/4461/11.17
41
120/330/0.23
120/330/0.25
118/286/0.19
124/298/0.20
42
3832/9998/2.51
3373/8928/7.18
3685/8619/8.30
17419/22467/5.83
43
40/80/0.11
40/80/0.08
99/119/0.19
99/119/0.19
44
50/110/0.10
50/110/0.05
50/110/0.02
70/282/0.13
45
43/184/0.08
43/184/0.11
12047/91618/9.93
15601/521624/4.14
46
427/1323/2.56
427/1320/2.60
409/1040/2.07
597/1232/2.51
47
64/249/0.09
64/249/0.06
133/380/0.15
143/393/0.40
48
308/798/1.39
293/772/1.28
2549/27198/2.61
15662/156182/6.10
49
22/89/0.11
22/89/0.13
1377/37674/8.58
1571/43940/8.14
50
20/50/0.03
20/50/0.02
20/50/0.02
20/50/0.03
51
92/1755/2.25
42/136/0.17
1620/5439/3.43
1680/46332/4.04
52
107/418/0.16
107/418/0.15
125/360/0.14
125/360/0.16
53
6199/52426/6.75
6199/52426/6.78
4160/9534/3.23
2233/8070/9.66
54
51/151/0.20
51/151/0.19
75/173/0.25
77/144/0.20
55
60/140/0.19
60/140/0.19
60/120/0.14
60/120/0.17
56
70/140/0.20
70/140/0.17
70/140/0.24
80/160/0.22
57
79/158/0.22
79/158/0.22
79/158/0.24
86/172/0.23
58
143/570/0.25
143/570/0.25
177/506/0.25
177/506/0.25
59
188/498/0.29
172/453/0.27
209/460/0.30
641/917/1.01
60
83/236/0.13
83/236/0.11
1615/43549/5.17
2005/57225/7.75
总
37602/124896/60.1
37426/121943/70.14
167737/10741186/139.61
271227/7136001/214.33
因此,在这些重心算法,新算法似乎产生最好的搜索方向。在表3 ,有一个明确的证据表明,新算法优于标准的LS和CD算法详细如下(当c = 0.1):(一) 100% LS算法:该算法提高了(71.5%),由(94.92%)NOFG改进,提高了(64.2%) (b) 100%的CD算法:该算法提高了(72.6%),由(89.95%)NOFG改进,提高了(66.7%)
算法
工具
LS (1991) (%)
CD (1987) (%)
当c = 0.1
随机过程
28.5
27.4
NOFG
5.08
10.05
时间
35.8
33.3
当c = 0.001
随机过程
28.3
27.3
NOFG
4.9
9.7
时间
38.1
35.5
(当c = 0.001):(c) 100% LS算法:该算法提高了(64.34%),由(16.99%)NOFG改进,提高了(16.08%) (d) 100%的CD算法:该算法提高了(52.60%),由(14.63%)NOFG改进,提高了(12.75%) 在表4 ,有一个清晰明显,新算法优于标准的商品和公关算法如下详细(当c = 0.1): (e) 100% HS算法:该算法提高了(77.6%),由(98.9%)NOFG改进,提高了(57%) (f) 100%公关算法:该算法提高了(86.2%),由(98.3%)NOFG改进,提高了(72%) (当c = 0.001) (g) 100% HS算法:该算法提高了(77.7%),由(98.9%)NOFG改进,提高了(49.8%) (h) 100%公关算法:该算法提高了(86.3%),由(98.3%)NOFG改进,提高了(67.3%)
算法
工具
海关(1952)(%)
公关(1969)(%)
当c = 0.1
随机过程
22.4
13.8
NOFG
1。1
1。7
时间
43
28
当c = 0.001
随机过程
22.3
13.7
NOFG
1。1
1。7
时间
50.2
32.7
从上面的表格我们可以推导出和实验总结如下:(我) 点(d)上面,我们新提出的算法CG-type方法领域的经济和健壮的标准相比LS和CD算法 (2) 上述点(e h)是我们的新提出的算法CG-type方法领域的经济和健壮的HS和PR算法比标准
所有这些比较都是使用性能的多兰和更多16 ),我们可以得出这样的结论:(1) 图1 说明了新算法和(LS、CD、HS和公关)的活动新算法在计算迭代的数量 (2) 图2 解释了新算法和(LS、CD、HS和公关)的活动新算法在计算功能和梯度的数量评估 (3) 图3 显示的算法需要多长时间到达解决方案(即。所需的CPU时间) (4) 图4 显示了执行的函数在两个不同的常数的新算法比基本算法(LS、CD、HS和公关)基于迭代的数量 (5) 图5 演示了许多的杰出的表现相比,新算法与两个不同的常数函数基本算法(LS、CD、HS和公关)基于函数梯度的数量评估
8。结论
在这项研究中,我们已经提交了两个新CG方法(通过改变的价值c )。提出的一个重要属性CG方法是保证充分下降方向。在温和的条件下,我们证明了新算法是全局收敛的一致凸和一般函数使用强大的沃尔夫行搜索条件。初步的数值结果表明,如果我们决定好的参数值c ,新算法能够很好地执行。然而,一个最优值的参数c 理论上可以处理(在未来研究)实现更多最好的数值结果。
数据可用性
使用的数据支持本研究的发现可以从相应的作者。
的利益冲突
作者宣称没有利益冲突有关的出版。
确认
支持的研究是计算机科学和数学学院大学的摩苏尔,伊拉克共和国,在项目没有。3615208。