27.10.2014 Views

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

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

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

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!