Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
f 13 = d 4<br />
(<br />
d 2<br />
√<br />
b14 2 2 3 3 4 1<br />
− d 5<br />
√<br />
b14 3 3 4 1 5 2<br />
) 2<br />
+ d2<br />
(d 4<br />
√<br />
b14 2 1 3 3 4 2<br />
− d 5<br />
√<br />
b14 2 1 3 3 5 2<br />
) 2+<br />
d 5<br />
(d 2<br />
√<br />
b14 2 2 3 3 5 1<br />
− d 4<br />
√<br />
b14 3 3 4 2 5 1<br />
) 2<br />
+ µ31 d 2 d 4 d 5 ;<br />
f 14 = d 3<br />
( √b14<br />
2 2 3 1 4 3<br />
d 2 − √ b 14 3 1 4 3 5 2<br />
d 5<br />
) 2<br />
+ d2<br />
( √b14<br />
2 1 3 2 4 3<br />
d 3 − √ b 14 2 1 4 3 5 2<br />
d 5<br />
) 2+<br />
d 5<br />
( √b14<br />
2 2 4 3 5 1<br />
d 2 − √ b 14 3 2 4 3 5 1<br />
d 3<br />
) 2<br />
+ µ41 d 2 d 3 d 5 ;<br />
f 15 = d 3<br />
( √b14<br />
2 2 3 1 5 3<br />
d 2 − √ b 14 3 1 4 2 5 3<br />
d 4<br />
) 2<br />
+ d2<br />
( √b14<br />
2 1 3 2 5 3<br />
d 3 − √ b 14 2 1 4 2 5 3<br />
d 4<br />
) 2+<br />
d 4<br />
( √b14<br />
2 2 4 1 5 3<br />
d 2 − √ b 14 3 2 4 1 5 3<br />
d 3<br />
) 2<br />
+ µ51 d 2 d 3 d 4 .<br />
(12)<br />
(13)<br />
(14)<br />
Другие значения f j и f jk могут быть получены из (10)-(14) соответствующей перестановкой<br />
индексов и субиндексов.Значения µ ij приведены в приложении (A.9). В каждом коэффициенте<br />
при d 4 i в общем случае выбором d j > 0 можно обратить в ноль одновременно все<br />
три скобки, поэтому требования µ ij 0 являются необходимыми для положительности Γ 4<br />
при любых d j > 0. Таких условий 20:<br />
µ ij = ( − A 2k,m,n a j,j<br />
)<br />
δij 0, ( k < m < n; i ≠ j ≠ k ≠ m ≠ n; i, j, k, m, n = 1, 5; δ ij = δ ji<br />
)<br />
,<br />
но они сводятся к десяти неравенствам вида:<br />
δ 12 = (√ −A 45 A 33,4 + √ −A 44 A 33,5 + √ −A 43 A 34,5<br />
) 2<br />
+ A23,4,5 ∆ 5 0,<br />
δ 13 = (√ −A 45 A 32,4 + √ −A 44 A 32,5 + √ −A 42 A 34,5<br />
) 2<br />
+ A22,4,5 ∆ 5 0,<br />
δ 14 = (√ −A 45 A 32,3 + √ −A 43 A 32,5 + √ −A 42 A 33,5<br />
) 2<br />
+ A22,3,5 ∆ 5 0,<br />
δ 15 = (√ −A 44 A 32,3 + √ −A 43 A 32,4 + √ −A 42 A 33,4<br />
) 2<br />
+ A22,3,4 ∆ 5 0,<br />
δ 23 = (√ −A 45 A 31,4 + √ −A 44 A 31,5 + √ ) 2<br />
−A 41 A 34,5 + A21,4,5 ∆ 5 0,<br />
δ 24 = (√ −A 45 A 31,3 + √ −A 43 A 31,5 + √ ) 2<br />
−A 41 A 33,5 + A21,3,5 ∆ 5 0,<br />
δ 25 = (√ −A 44 A 31,3 + √ −A 43 A 31,4 + √ ) 2<br />
−A 41 A 33,4 + A21,3,4 ∆ 5 0,<br />
δ 34 = (√ −A 45 A 31,2 + √ −A 42 A 31,5 + √ ) 2<br />
−A 41 A 32,5 + A21,2,5 ∆ 5 0,<br />
δ 35 = (√ −A 44 A 31,2 + √ −A 42 A 31,4 + √ ) 2<br />
−A 41 A 32,4 + A21,2,4 ∆ 5 0,<br />
δ 45 = (√ −A 43 A 31,2 + √ −A 42 A 31,3 + √ ) 2<br />
−A 41 A 32,3 + A21,2,3 ∆ 5 0.<br />
(15)<br />
Если посмотреть на Γ 4 (8) как на полином по любому d i , у которого коэффициенты представляют<br />
собой полиномы по оставшимся d j , (i ≠ j = 1, 5), то необходимо обеспечить<br />
также неотрицательность коэффициентов при нулевой степени d j в коэффициенте при<br />
d 4 i и в слагаемом, не содержащем переменной d i . Таким способом получим как необходимые<br />
условия положительности Γ 4 в положительном ортанте еще две группы неравенств,<br />
которые в терминах главных миноров имеют вид:<br />
β ij 0, (i, j = 1, 5, β ij = β ji ) (16)<br />
γ ij = (√ −A 3n,i A 2k,m,i + √ −A 3k,i A 2m,n,i + √ −A 3m,i A 2k,n,i<br />
) 2<br />
+ A4i a j,j 0,<br />
(<br />
i ≠ j ≠ k ≠ m ≠ n; i, j, k, m, n = 1, 5<br />
)<br />
.<br />
(17)<br />
36