27.10.2014 Views

Вычислительная математика - ИСЭМ СО РАН

Вычислительная математика - ИСЭМ СО РАН

Вычислительная математика - ИСЭМ СО РАН

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!