Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
умноженную на любое вещественное число α. В дальнейшем будем исследовать на знакоопределенность<br />
квадратичную форму<br />
V (z, α) = V 1 (z) + α V 2 (z) . (1.1)<br />
Приводя каким-либо способом [6, 7] квадратичную форму к сумме квадратов, по коэффициентам<br />
при этих квадратах можно сделать определенные заключения о знакоопределенности<br />
V (z, α). Выполним это самым простым методом последовательного исключения [6]<br />
и получим<br />
V (z, α) = (z 1 + α z 3 ) 2 + J 2 (α) [z 2 + gz 3 /2J 2 (α)] 2 + ∆(α) z 2 3 ,<br />
где обозначено: ∆(α) = J 3 (α)/J 2 (α); J 2 (α) = p − 2 α; J 3 (α) = 2 α 3 − p α 2 − r α +<br />
p r − g 2 /4.<br />
При этом переменные z 1 , z 2 , z 3 удовлетворяют условию связи вида V 2 (z) = 0, которую<br />
можно рассматривать ограничением на область определения переменных: z 1 ≥ 0; z 3 ≥ 0.<br />
Проведем вначале анализ квадратичной формы (1.1), затем установим ее свзь с формой<br />
(1). По критерию Сильвестра [6, 7] для положительной определенности V (z, α) достаточно,<br />
чтобы существовала область D ⊂ R (D ≠ 0) значений параметра α, в которой<br />
выполняется система неравенств:<br />
{<br />
J2 (α) > 0<br />
J 3 (α) > 0<br />
(1.2)<br />
Область D будем выбирать с целью приблизить достаточные условия знакоопределенности<br />
квадратичной формы к необходимым. Для этого исходя из условия J 2 (α) > 0,<br />
потребуем положительным наибольшее значение локального экстремума J 3 (α). Таким<br />
образом, поставлена задача условного экстремума:<br />
max J 3(α) > 0 .<br />
J 2 (α)>0<br />
Исследование условного экстремума проще провести простым способом: определить точку<br />
α 1 локального максимума функции f(α) = J 3 (α)/2, и затем проверить условие J 2 (α 1 ) ><br />
0. Искомое значение α 1 находится среди корней уравнения f ′ (α) = 3 α 2 − p α −<br />
r = 0. Локальному максимуму здесь соответствует меньший корень α 1 = (p −<br />
√<br />
p2 + 12 r)/6 < 0. Вычисления дают значения: f(α 1 ) = [36 p r − p 3 + √ (p 2 + 12 r) 3 −<br />
27 g 2 /2]/108; J 2 (α 1 ) = (2p + √ p 2 + 12 r)/3. Рассматривая полученные значения при<br />
разных знаках p, сведем требования (1.2) к системе<br />
{<br />
1. p > −2 √ r ,<br />
2. q 2 < 8 pr/3 − 2 p 3 /27 + 2 √ (p 2 + 12r) 3 /27<br />
(1.3)<br />
Любое отклонение от α 1 нарушает второе условие (1.3), тем самым сужает условия<br />
знакоопределенности (1.1). Следовательно, (1.3) являются необходимыми и достаточными<br />
условиями положительной определенности (1.1).<br />
135