
- •Экономико-математические методы и модели
- •П.И. Кузьмин
- •Рекомендуется к печати научно-методическим советом (кафедрой) Международного института экономики, менеджмента и информационных систем Алтайского государственного университета
- •Содержание
- •Введение
- •Глава 1. Основные понятия и этапы математического моделирования
- •1.1. Основные задачи экономико-математического моделирования.
- •1.2. Использование экономико-математических методов при принятии управленческих решений.
- •Глава 2. Функции многих переменных. Квадратичные формы.
- •2.1. Квадратичные формы. Критерий Сильвестра. Достаточное условие оптимальности.
- •2.2. Представление прибыли как функции многих переменных.
- •Задание к лабораторной работе №1 по теме «Квадратичные формы. Функции многих переменных».
Глава 2. Функции многих переменных. Квадратичные формы.
2.1. Квадратичные формы. Критерий Сильвестра. Достаточное условие оптимальности.
С линейными функциями многих переменных студентам пришлось иметь дело в задачах линейного программирования.
При решении различных прикладных задач часто приходится исследовать квадратичные формы [Мальцев, Кремер].
Определение 1. Квадратичной формой L(xl,x2,...,xn) от п переменных называется сумма, каждый член которой является либо квадратом одной из переменных, либо произведением двух разных переменных, взятых с некоторым коэффициентом:
, (1)
Предполагаем, что
коэффициенты квадратичной формы аij
действительные
числа, причем аij
= аji.
Матрицу А
всегда можно предполагать симметрической
(показать это на практике). Действительно,
значение L(X)
не измениться, если каждый из пары
коэффициентов
и
заменить на
.
Матрица А
= (аij)
(i, j = =
1, 2, ..., n),
составленная
из этих коэффициентов,
называется
матрицей
квадратичной формы.
В матричной записи квадратичная форма имеет вид:
L(X) = X'AX, (2)
где X = (xl, х2,..., хn )' — матрица-столбец переменных. В самом деле
.
и эквивалентность формул (1) и (2) установлена.
Пример 2.1. Дана квадратичная форма L(x1,x2,x3) = 2x12 –12x1x2 - 10x1x3 +4x22-5x32 . Записать ее в матричном виде.
Решение. Найдем матрицу квадратичной формы. Ее диагональные элементы равны коэффициентам при квадратах переменных, т.е. 2, 4, -5, а другие элементы половинам соответствующих коэффициентов квадратичной формы. Поэтому
.
Выясним, как изменяется квадратичная форма при невырожденном линейном преобразовании переменных.
Пусть вектор-столбцы переменных X' = (x1, x2,...,xn)' и Y'= (y1, y2, …., yn)' связаны линейным преобразованием X = CY, где C=(cij) (i,j = 1, 2, ...,n) есть некоторая невырожденная матрица n-го порядка. Тогда квадратичная форма после применения преобразования имеет вид
L = X'AX = (CY)'A(CY) = (Y'С')A(СY) = Y'(C'AC)Y.
Итак, при невырожденном линейном преобразовании X = CY матрица квадратичной формы принимает вид:
В = С'АС. (3)
Пример 2.2. Дана квадратичная форма L(x1,x2) = 2x12 –4x1x2 – 3x22. Найти квадратичную форму L(у1,y2), полученную из данной линейным преобразованием xl = 2у1 - Зу2, х2 = у1 + у2
Решение.
Матрица данной квадратичной формы A
=
,
а матрица линейного
преобразования C
=
Следовательно, по (3) матрица искомой квадратичной формы
=
а квадратичная форма имеет вид L(y1, у2) = 13у12 – 34y1y2 + Зу22. Следует отметить, что при некоторых удачно выбранных линейных преобразованиях вид квадратичной формы можно существенно упростить.
Определение 2. Квадратичная форма
называется канонической (или имеет к анонический или диагональный вид), если все ее коэффициенты аij = 0 при ij:
.
а ее матрица является диагональной. Справедлива следующая теорема.
Теорема 1. Любая квадратичная форма с помощью невырожденного линейного преобразования переменных может быть приведена к каноническому виду.
Пример 2.3. Привести к каноническому виду квадратичную форму методом выделения полного квадрата
L(x1,x2,x3) = x12 – 3x1x2 + 4x1x3 + 2x2x3 +x32
Решение. Вначале выделим полный квадрат при переменной х1, коэффициент при квадрате которой отличен от нуля:
L(x1, x2, x3) = [x12 – 2x1((3·(x2 – 4x3)/2)+((3x2 – 4x3)/2)2] –
– ((3x2 – 4x3)/2)2+2x2x3 + x32 =
= (x12 – 3·x2/2+ 2·x3)2 – 9·x22/4 + 6·x2·x3 – 4·x32 + 2·x2·x3 + x32 =
= (x12 – 3·x2/2+ 2·x3)2 – 9·x22/4 + 8·x2·x3 – 3·x32.
Теперь выделяем полный квадрат при переменной х2, коэффициент при которой отличен от нуля:
L(x1,
x2,
x3)
=
.
'Итак, невырожденное линейное преобразование
приводит исходную квадратичную форму к каноническому виду
L(y, y2, y3) = y12 – 9y22/4 + 37·y32/9.
Канонический вид квадратичной формы не является однозначно определенным, так как одна и та же квадратичная форма может быть приведена к каноническому виду многими способами, однако полученные различными способами канонические формы обладают рядом общих свойств. Одно из этих свойств формулируем в виде теоремы.
Теорема 2 (закон инерции квадратичных форм). Число слагаемых с положительными (отрицательными) коэффициентами квадратичной формы в каноническом виде не зависит от способа приведения формы к этому виду.
Например, квадратичную форму L в примере 2.3 можно было привести к виду
L(y, y2, y3) = 37y12/4 + y22 y32,
применив невырожденное линейное преобразование
.
Как видим, число положительных и отрицательных коэффициентов (соответственно два и один) сохранилось.
Следует отметить, что ранг матрицы квадратичной формы, называемый рангом квадратичной формы, равен числу отличных от нуля коэффициентов канонической формы и не меняется при линейных преобразованиях.
Определение 3. Квадратичная форма L(х1,х2,...,хn) называется положительно (отрицательно) определенной, если при всех значениях переменных, из которых хотя бы одно отлично от нуля,
L(х1,х2,...,хn) > 0 (L(х1,х2,...,хn) < 0).
Так, например, квадратичная форма
L1 = 3x12 + 4x12 + 9x32
является положительно определенной, а форма
L2 = 10x12 + 2x1x2 5x22
— отрицательно определенной.
Нарисуем трехмерный график (пространственную форму) отрицательно-определенной квадратичной формы.
Рис.1. График отрицательно-определенной формы.
Теорема 3. Для того чтобы квадратичная форма L = Х'АХ была положительно (отрицательно) определенной, необходимо и достаточно, чтобы все собственные значения i, матрицы А были положительны (отрицательны).
В ряде случаев для установления знакоопределенности квадратичной формы удобнее бывает применить критерий Сильвестра.
Теорема 4 (критерий Сильвестра). Для того чтобы квадратичная форма была положительно определенной, необходимо и достаточно, чтобы все главные миноры матрицы этой формы были положительны, т.е 1> 0, 2> 0,…, n>0, где
.
Следует отметить, что для отрицательно определенных квадратичных форм знаки главных миноров чередуются, начиная со знака "минус" для минора первого порядка.
Доказательство. При n=1 теорема верна, так как форма имеет в этом случае вид ах2 и поэтому положительно определена тогда и только тогда, если а>0. Будем поэтому доказывать теорему для случая n неизвестных, предполагая, что для квадратичных форм от n 1 неизвестных она уже доказана.
Сделаем сначала следующее замечание:
Если квадратичная форма L с действительными коэффициентами, составляющими матрицу A подвергается невырожденному линейному преобразованию с действительной матрицей Q, то знак определителя квадратичной формы (т.е. определителя ее матрицы) не меняется.
Действительно, после преобразования мы получаем квадратичную форму с матрицей QТAQ, однако, ввиду |QТ| = |Q|,
QТAQ = QТAQ=AQ2,
т. е. определитель |А| умножается на положительное число.
Пусть теперь дана квадратичная форма
Ее можно записать в виде
,
(4)
где будет квадратичной формой от n1 неизвестных, составленной из тех членов формы L, в которые не входит неизвестное хn.
Главные миноры формы совпадают, очевидно, со всеми, кроме последнего, главными минорами формы L.
Необходимость. Пусть форма L положительно определена. Форма также будет в этом случае положительно определенной: если бы существовали такие значения неизвестных х1, х2, …, хn-1, не все равные нулю, при которых форма получает не строго положительное значение, то, полагая дополнительно хn = 0, мы получили бы, ввиду (4), также не строго положительное значение формы L, хотя не все значения неизвестных x1, х2, ... , xn-1, хn равны нулю. Поэтому, по индуктивному предположению, все главные миноры формы , т. е. все главные миноры формы L, кроме последнего, строго положительны. Что же касается последнего главного минора формы L, т. е. определителя самой матрицы А, то его положительность вытекает из следующих соображений: форма L, ввиду ее положительной определенности, невырожденным линейным преобразованием приводится к нормальному виду, состоящему из n положительных квадратов. Определитель этого нормального вида строго положителен, а поэтому ввиду сделанного выше замечания положителен и определитель самой формы L.
Достаточность. Пусть теперь строго положительны все главные миноры формы L. Отсюда вытекает положительность всех главных миноров формы , т. е., по индуктивному предположению, положительная определенность этой формы. Существует, следовательно, такое невырожденное линейное преобразование неизвестных х1, х2, ... , хn-1, которое приводит форму к виду суммы n-1 положительных квадратов от новых неизвестных у1, у2, ... , уn-1. Это линейное преобразование можно дополнить до (невырожденного) линейного преобразования всех неизвестных xl, x2, ..., хn, полагая хn= уn. Ввиду (4) форма L приводится указанным преобразованием к виду
;
(5)
точные выражения коэффициентов bin для нас несущественны. Так как
,
то невырожденное линейное преобразование
,
,
приводит, ввиду (5), форму L к каноническому виду
.
(6)
Для доказательства положительной определенности формы L остается доказать положительность числа с. Определитель формы, стоящей в правой части равенства (6), равен с. Этот определитель должен, однако, быть положительным, так как правая часть равенства (6) получена из формы L двумя невырожденными линейными преобразованиями, а определитель формы L был, как последний из главных миноров этой формы, положительным.
Доказательство теоремы закончено.
Пример 2.4. Доказать, что квадратичная форма L = 13x12 6x1x2 + 5x22 является положительно определенной.
Решение. Первый способ. Матрица А квадратичной формы имеет вид
A
=
Для матрицы А характеристическое уравнение
.
Решая уравнение, найдем 1 = 14, 2 = 4. Так как корни характеристического уравнения матрицы А положительны, то на основании приведенной теоремы квадратичная форма L — положительно определенная.
Второй способ. Так как главные миноры матрицы А
положительны, то по критерию Сильвестра данная квадратичная форма L положительно определенная.
Нарисуем трехмерный график (пространственную форму) нашей квадратичной формы L = 13x12 6x1x2 + 5x22
Рис.2. График положительно-определенной формы.
Таким образом, график положительно-определенной квадратичной формы имеет такой вид. Он ограничен областью определения, имеющей вид прямоугольника. На самом деле трехмерный график неограничен.
Определение 4.
Квадратичная форма называется:
положительно определенной, если L(X)>0 для всех X0;
положительно полуопределенной, если L(X)0, для X и X0 такое, что L(X)=0.
отрицательно определенной, если (–L(X)) есть положительно определенная квадратичная форма;
отрицательно полуопределенной, если (–L(X)) – положительно полуопределенная квадратичная форма;
неопределенной – в остальных случаях.
Теорема 5 (Необходимые и достаточные условия положительной (отрицательной) определенности квадратичной формы):
Квадратичная форма L(X) является положительно определенной , когда значения всех угловых главных миноров матрицы А положительны (матрица А – называется положительно определенной матрицей).
Квадратичная форма L(X) является отрицательно определенной , когда значение к-го углового главного минора матрицы А имеет знак
(А–отрицательно определенная матрица).
Квадратичная форма L(X) является положительно полуопределенной , когда А – вырожденная матрица и все ее главные миноры неотрицательны (А – положительно полуопределенная матрица).