​​T=Dimensions[qw];​​delays={6,12,18,24};​​maxdelay=Max[delays];ranges=Range[maxdelay+1,T[[2]]];X={};For[d=1,d<maxdelay+1,d+=6,AppendTo[X,qw[[2,ranges-d]]]];t=qw〚2,ranges〛;​​(*Y//MatrixForm*)​​X=X;​​Dimensions[t];​​Dimensions[X][[1]];
Out[]=
{176}
​
a=0.5;n=4;K[u_,v_]:=
n
∏
i=1
Cos[1.75(u[[i]]-v[[i]])/a]Exp-
2
(u[[i]]-v[[i]])
2
2
a
;
ϵ=0.0025;c=200.;m=Length[X];
In[]:=
M=(Table[N[K[X[[i]],X[[j]]]],{i,1,m},{j,1,m}]+(1/c)IdentityMatrix[m]);
In[]:=
W=
m
∑
i=1
α
i
t[[i]]-ϵ
m
∑
i=1
Abs[
α
i
]-(1/2)
m
∑
i=1
m
∑
j=1

α
i
α
j
M[[i,j]];
In[]:=
q=Apply[And,Join[Table[-c<
α
i
≤c,{i,1,m}],{(
m
∑
i=1
α
i
==0)}]];
In[]:=
vars=Table[
α
i
,{i,1,m}];
sol=NMaximize[{W,q},vars,AccuracyGoal9,PrecisionGoal8,WorkingPrecision8,MaxIterations10,Method{"DifferentialEvolution"}]
NMaximize
:Failed to converge to the requested accuracy or precision within 10 iterations.
Out[]=
{1144.2448,{
α
1
-6.0700,
α
2
-7.4024114,
α
3
-5.7326631,
α
4
-4.6291851,
α
5
-11.786311,
α
6
-2.9104057,
α
7
-1.9592662,
α
8
-3.7530973,
α
9
-6.0099946,
α
10
-5.7345274,
α
11
-3.4237032,
α
12
-3.4486595,
α
13
-5.8556800,
α
14
-6.2094720,
α
15
-10.800526,
α
16
-5.6820036,
α
17
-3.8399903,
α
18
-8.0607139,
α
19
-6.3768424,
α
20
-10.691821,
α
21
-5.4710934,
α
22
-5.0760439,
α
23
-11.585661,
α
24
-11.738912,
α
25
-7.9214204,
α
26
-9.8717513,
α
27
-6.8719561,
α
28
-10.075113,
α
29
-3.1992382,
α
30
-8.6485663,
α
31
-5.6561623,
α
32
-7.6426037,
α
33
-4.8648042,
α
34
-4.4304463,
α
35
-5.2717089,
α
36
-2.5420508,
α
37
-2.5779218,
α
38
-1.2886771,
α
39
0.32649711,
α
40
-0.31100478,
α
41
-1.0551780,
α
42
-4.0034025,
α
43
-4.0996268,
α
44
-5.4967829,
α
45
-0.73879925,
α
46
-0.81033068,
α
47
-1.5929827,
α
48
-2.0208546,
α
49
-1.7840498,
α
50
-5.6269995,
α
51
-5.1464830,
α
52
-3.3854392,
α
53
-3.8577302,
α
54
-7.0771051,
α
55
-4.5628406,
α
56
-13.099354,
α
57
-8.6143327,
α
58
-3.8968887,
α
59
-2.9263743,
α
60
-0.99577529,
α
61
0.75226440,
α
62
-6.4953370,
α
63
3.0061613,
α
64
-11.202525,
α
65
-1.7077698,
α
66
-15.305809,
α
67
-9.0412460,
α
68
0.91062393,
α
69
4.5561567,
α
70
-7.2194150,
α
71
-0.29889095,
α
72
7.5314084,
α
73
-0.32559145,
α
74
-6.3202603,
α
75
3.9364789,
α
76
0.37912606,
α
77
3.7555046,
α
78
3.8115240,
α
79
3.7908179,
α
80
-1.5902621,
α
81
0.81599287,
α
82
4.8925599,
α
83
7.6202917,
α
84
1.8924110,
α
85
4.1145162,
α
86
3.5132104,
α
87
3.9260567,
α
88
16.140014,
α
89
-2.5039513,
α
90
3.6725529,
α
91
21.285196,
α
92
-6.5647872,
α
93
2.8866314,
α
94
-8.9676051,
α
95
9.4777536,
α
96
7.6862937,
α
97
38.887996,
α
98
53.409206,
α
99
4.1827239,
α
100
-12.980137,
α
101
3.1224981,
α
102
0.26287974,
α
103
72.828057,
α
104
-58.761255,
α
105
-2.3989716,
α
106
-26.066234,
α
107
2.7358108,
α
108
0.59531752,
α
109
-3.9839312,
α
110
29.404325,
α
111
-4.7624919,
α
112
25.646572,
α
113
16.662786,
α
114
-6.5311482,
α
115
9.4909358,
α
116
7.1531043,
α
117
8.0936983,
α
118
-8.5287916,
α
119
2.0317210,
α
120
4.3570097,
α
121
27.138310,
α
122
32.466272,
α
123
12.405534,
α
124
29.160774,
α
125
14.889394,
α
126
5.6039125,
α
127
11.195291,
α
128
-5.6021723,
α
129
3.8147558,
α
130
2.9431486,
α
131
1.6867925,
α
132
-1.6833426,
α
133
-0.37688247,
α
134
3.7610989,
α
135
4.6697241,
α
136
3.8732275,
α
137
6.2357974,
α
138
-2.6598779,
α
139
2.5822290,
α
140
-3.3320785,
α
141
-6.1449289,
α
142
4.6978221,
α
143
0.21099083,
α
144
-2.2270633,
α
145
-4.9318945,
α
146
-0.41422318,
α
147
-2.8155336,
α
148
4.3707475,
α
149
-3.3190804,
α
150
0.088618572,
α
151
3.2267303,
α
152
1.6302329,
α
153
-0.14157510,
α
154
0.49520176,
α
155
-0.59107908,
α
156
-0.077752602,
α
157
-2.9628083,
α
158
-2.4542100,
α
159
3.9202894,
α
160
-4.6026854,
α
161
-2.0516871,
α
162
-7.4330803,
α
163
-5.4603897,
α
164
-2.8948489,
α
165
-3.2569607,
α
166
-3.4275100,
α
167
-2.0227517,
α
168
-2.8821555,
α
169
2.8395081,
α
170
1.9390556,
α
171
1.1343286,
α
172
0.75382804,
α
173
2.7805756,
α
174
4.6071561,
α
175
-1.6815581,
α
176
-1.4477215}}