Mathematical Programming Textbook 9965296308 Compress
Mathematical Programming Textbook 9965296308 Compress
MATHEMATICAL
PROGRAMMING
Textbook
Almaty
"Kаzаkh university"
2011
1
ББК 22. 1
А 36
Рекомендовано к изданию
Ученым советом механико-математического факультета
и РИСО КазНУ им. аль-Фараби
(Протокол №1 от 19 октября 2010 г.)
Reviewers:
doctor of physical and mathematical sciences, professor M.Akhmet;
candidate of philological sciences, professor A.A.Moldagalieva
ББК 22. 1
Lecture 1
THE MAIN DETERMINATIONS.
STATEMENT OF THE PROBLEM
J (u ) J (u1 , u 2 , u 3 , u 4 , u 5 ) c1u1 c 2 u 2 c3 u 3 c 4 u 4 c5 u 5
3u1 5u 2 10u3 u 4 8u5 max (1)
u1 0, u 2 0, u3 0, u 4 0, u5 0, (3)
5
c1 u1
c2 u2 b1 g1
c c3 ; u u3 ; b b2 ; g g 2 ,
c4 u4 b g
3 3
u
c5 5
g (u ) Au b; ( 2' )
u 0, (3' )
u U u E 5 u U 0 , g (u ) 0 E 5 . (5)
The problem (1) – (3) (or (1' ) – (3' ) , or (4), (5)) belongs to type
so called problem of the linear programming, since J (u ) is linear
with respect to u function; the vector-function g (u ) also is linear
with respect to u .
In the case above vector u has five components, vector-function
g (u ) – three and matrix A has an order 3 5 was considered. In the
case vector u has dimensionality n , g (u ) – m -a measured
function, but matrix A has an order m n and there are restrictions
of the equality type (for instance g1 (u ) A1u b1 , where A1 – a
6
matrix of the order m1 n , b1 E 1 ) problem (4), (5) is written:
m
u U ,
n
U u E u U 0 , g ( u ) Au b 0 , g1( u ) A1u b1 0 , (7)
where set U 0 u E n u 0 , c E n . Problem (6), (7) belongs
to type of the general problem of the linear programming.
We suppose that J (u ) is convex function determined on the
convex set U 0 (unnecessary linear); g (u ) is a convex function
determined on the convex set U 0 . Now problem (6), (7) is written as
J ( u ) max; (8)
u U u E n u U 0 , g ( u ) 0 , g1( u ) A1u b1 0 . (9)
J (u ) max; (10)
u U u E n u U 0 , g (u ) 0, g 1 (u ) 0. (11)
J (u ) max; u U , U En.
7
Finally, we note that problem of the function maximization J (u )
on a set U tantamount to the problem of the function minimization –
J (u ) on set U. So further it is possible to limit by consideration of the
problem
J (u ) min; u U , U E n , (12)
value of the function J (u ) on set U.
The set U * u* U J (u* ) min J (u ) contains all minimum
uU
points of the function J (u ) on set U.
It follows from the definition, that the global (or absolute)
minimum of the function J (u ) on the set U is reached on the set
U * U . We remind that a point u** U is called the local
minimum point of the function J (u ) , if inequality J (u** ) J (u )
is valid under all u o(u** , ) U , where set
o(u** , ) u E n | u u** | E n
a
i 1
2
i
.
8
, and U u E 1 1 2 u 1.
Example 1. Let J (u ) cos 2
u
Then set U * 2 3 ; b) U u E 1 1 4 u 2 , then set
U * 2 / 7, 2 / 5, 2 / 3, 2 ; c) U u E 1 / 2 u ,
then set
U * , where - empty set.
Example 2. The function J (u ) ln u , set U u E 1 0 u 1.
The set U * .
Example 3. The function J (u ) J 1 (u ) c, c const , where
function
| u a |, если u a;
J 1 (u ) u b, если u b;
c, если a u b,
min J u does not always exist. Since lower bound J * for function
uU
J u inf, u U , U E n . (13)
11
Lecture 2
WEIERSTRASS’S THEOREM
J u inf, u U , U E n . (1)
12
Definition 1. The set U E n is identified compact, if any sequence
u k U has at least one limiting point v, moreover v U .
It is easy make sure the definition is equally to known from the
course of the mathematical analysis statement about any bounded
and closed set is compactly. In fact, according to Bolzano and
Weierstrass theorem any bounded sequence has at least one limiting
point (the set U is bounded ), but from inclusion v U follows the
property of being closed of the set U .
Half-continuously from below. Let u k E n be a sequence.
Then J k J uk is the number sequence. We notice that: a)
numerical set J k is bounded from below, if the number exists,
such that J k , k 1,2,3,... ; b) numerical set J k is not bounded
from below, if J km subsequence exists such the limit lim J km .
m
Definition 2. By the lower limit of the bounded from below numeric
sequence J k is identified the value а denoted a lim J k if: 1)
m
subsequence J km exists for which lim J km a ; 2) all other limiting
m
U * u* E n u* U , J (u* ) min J (u )
uU
is inempty, compactly and any minimizing sequence converges to set
U* .
Proof. Let u k U be any minimizing sequence, i.e. the limit
of the numeric set lim J u k J * . We notice, that such minimizing
k
m
J * J u* lim J u km lim J u k J * .
k
(2)
M (v ) u E n u U , J (u ) J (v )
16
is bounded . Then J * , set U * is inempty, compact and any
minimizing sequence u k M (v) converges to set U * .
Proof. Since all condition of the lemma are executed, set M (v )
is closed. From restrictedness and closeness of the set M (v ) follows
its compactness. The set U M (v) M 1 (v) , where set M 1 (v)
u E n u U , J (u ) J (v) , moreover on set M 1 (v) function
J (u ) does not reach its lower bound J * . Hereinafter proof of the
theorem 1 is repeated with change set U on compact set M (v ) .
Theorem is proved.
Theorem 3. Let function J (u ) be determined, finite and half-
continuous from below on inempty closed set U E n . Let for any
sequence vk U , vk , under k the equality
lim J (v k ) is executed. Then J * , set U * is inempty,
k
J (v) J * exists. We enter Lebesgue’s set M (v) u E n u U ,
J (u ) J (v ) . With consideration of lemma set M (v ) is closed. It is
easy to show that set M (v ) is bounded. In fact, if set M (v ) is not
bounded, then the sequence wk M (v) , such that wk under
k exists. By condition of the theorem for such sequences the
value J ( wk ) under k . Since J ( wk ) J (v) it is not
possible. Finally, set M (v ) is bounded and closed, consequently, it is
compact.
Hereinafter proof of the theorem 1 is repeated for set M (v ) .
Theorem is proved.
17
Chapter I. CONVEX ROGRAMMING.
ELEMENTS OF THE CONVEX
ANALYSIS
Lecture 3
CONVEX SETS
J (u ) inf,
u U u E n u U 0 , g i (u ) 0, i 1, m, (1)
g i (u ) a i , u bi 0, i m 1, s ,
18
Definition 1. Set U E n called by convex, if for any u U ,
v U and under all , 0 1 , the point u u (1 )v
v (u v) U .
Example1. We show that closed sphere
S (u 0 , R ) u E n u u 0 R
Is a convex set. Let the points u S (u 0 , R ) , S (u 0 , R ) , i. e.
norms u u 0 R , v u 0 R . We take a number , 0, 1 ,
and define the point u u (1 )v . Norm
u u 0 u (1 )v u 0 (u u 0 ) (1 )(v u 0 )
u u 0 (1 ) v u 0 R (1 ) R R.
u E n c, u
is a convex set, where с E n is the vector, is the number. Let the
points u , v . Consequently, scalar products c, u ,
c, v . Let u u (1 )v , 0, 1 . Then scalar product
c, u c, u (1 )v c, u
(1 ) c, v (1 ) .
19
vector which is convex. For the point u u (1 )v , u М ,
v М , 0, 1 we have
nr
M u E n u u 0 v, v L , L u E n v i u i
i 1
20
points a b c, b B, c C . The set A D , where is a
real number, if it contains the points a d , d D.
Theorem 1. If the sets A1 , A2 ,..., Am , B, C , D are convex, so
the sets A A1 A2 ... Am , A B C , A D are convex.
m m
Proof. Let the points be a ai A, e ei A,
i 1 i 1
U u E 1 0 u 1; 2 u 3 E 1 .
For the example the sets
int U u E 1 0 u 1; 2 u 3 .
affU u E 1 0 u 3, riU u E 1 0 u 1; 2 u 3.
the following statements are faithful:
1) If А is a convex set, so closing A is convex too. In fact, if
a, b A , so subsequences ak A, bk A, such that ak a,
bk b under k exists. With consideration of convexity of the
set A the point a k 1 bk A, k 1,2,3,... . Then limiting
point a lim a k 1 bk a 1 b A under all
k
v v u 0 v int U , u 0 inf U , v U , 0 1 .
22
To prove.
3) If U is a convex set, so int U is convex too. To prove.
4) If U is a convex inempty set, so ri U and ri U is
convex. To prove.
5) If U is a convex set and ri U , so the point
v v v v u 0 v riU , u 0 riU , v U , 0 1
To prove.
Definition 5. The point u E n called the convex combination of
the points u 1 , u 2 ,..., u m from E n , if it is represented in the manner
m
of u u ,
i 1
i
i
where the numbers i 0, i 1, m but their sum
1 2 ... m 1 .
Theorem 2. The set U is convex if and only if, when it contains
all convex combinations of any finite number their own points.
Proof. Necessity. Let U be a convex set. We show, that it
contains all convex combinations of the finite number of their own
points. We use the method of induction to proof of the theorem.
From definition of the convex set follows that statement faithfully for
any two points from U. We suppose, that set U contains the convex
combinations m 1 of their own points, i.e. the point
m 1
v i u i U , i 0, i 1, m 1, 1 ... m1 1 . We prove,
i 1
that it contains the convex combinations т their own points. In fact,
m 1
i
expression u 1u 1 ... m u m 1 m 1 u i mu m .
i 1 m
v i
i
, where i i 0, i 1, m, i 1 i 0, i 1, p
i 1
m p
moreover the sum
i 1
i
i 1
i 1 . Thence follows that u W ,
24
inclusion Co U W exists. From inclusion W Co U , Co U W
follows that Co U W . Theorem is proved.
Theorem 4. Any point u Co U can be presented as convex
combination no more than n 1 points from U .
Proof. As follows from theorem 3, the point u Co U is
m m
represented in the manner of u ,
i 1
i
i
i 0, i 1
i 1 . We
i 1
i 0 , that such i* 0 always exists. Now the point
m m m
u i u i i u i , i 0, i 1. (2)
i 1 i 1 i 1
i i* i i*
25
u 1
u 0 , i 1, m are linear independent and it is denoted by S m .
The points u 0 , u1 ,..., u m are called the tops of the simplex.
Set S m is a convex polyhedron with dimension m and by
theorem 3 it is represented in the manner of
m m
S m u E n u i u i , i 0, i 0, m, i 1.
i 1 i 0
26
Lecture 4
CONVEX FUNCTIONS
J u 1 v J u 1 J v . (1)
J u 1 v J u 1 J v 1 | u v |2 . (2)
27
the point u u 1 v U on the strength of convexity of the
set U . Then
J u c, u c, u 1 c, v J u 1 J v .
J u u 1 v u 1 v, u 1 v
2
2 | u |2 2 1 u, v 1 | v |2 . (3)
2
u v, u v | u v | 2 | u | 2 2 u, v | v |2 .
J u | u |2 1 | v |2 1 | u v |2
J u 1 J v 1 | u v |2 .
28
J u u 1 v | u | 2 1 | v | 2 J u 1 J v ,
2
29
Sufficiency. Let for function J u C 1 (U ) , U be a convex set
inequality (4) is executed. We show, that J (u ) is convex on U.
Since the set U is convex, the point u u 1 v U ,
u , v U , 0,1 . Then from inequality (4) follows, that
J u J u J ' u , u u , J v J u J ' u , v u ,
u , u , v U . We multiply the first inequality on , but the
second inequality - on the number 1 and add them. As a result
we give J u 1 J v J u . Thence follows the convexity
to functions J (u ) on set U. Theorem is proved.
Theorem 2. In order the function J (u ) C 1 (U ) to be convex
on convex set U necessary and sufficiently executing of the inequality
u , u h U .
30
The first equality follows from formula of the finite increments,
but the second - from theorem about average value, 0 1 1 Then
difference
J u 1 J u J u J u J u
1
1 J v J u J ' u t u u , u u dt
0
1
1 J ' u t v u , v u dt.
0
Let z1 u t u u u t 1 u v , z 2 u
t v u u t v u . Then z1 z 2 t u v , but
differences u u 1 z1 z 2 / t , v u z1 z 2 / t .
Now previous inequality is written:
1
1
J u 1 J v J u 1 J ' z1 J ' z 2 z1 z 2 dt.
0
t
J " (u ) , 0, E n , u U . (6)
J ' (u ) J ' (u ), 0, E n , | | 0 .
32
2 J u / u12 2 J u / u1u 2 ... 2 J u / u1u n
2
J u / u 2 u1 2 J u / u 22 ... 2 J u / u 2 u n –
J " u
. . . . . . . . . . . . . . . . . . . . . . . .
2 J u / u u 2 J u / u u ... 2 J u / u n2
n 1 n 2
a scalar product.
J " u , , ' J " u .
The theorems 1 - 3 for strongly convex function J (u ) on U are
formulated in the manner specified below and its are proved by the
similar way.
Theorem 4. In order the function J (u ) C 1 (U ) to be strongly
convex on the convex set U necessary and sufficiently executing of
the inequality
J " u , | | 2 , E n , u U ,
(9)
k const 0.
The formulas (4) - (9) can be applied for defining of the
convexity and strong convexity of the even functions J (u )
determined on the convex set U E n .
For studying of convergence of the successive approximation
methods are useful the following lemma.
33
Definition 3. It is spoken that gradient J ' u to function
J (u ) C 1 (U ) satisfies to Lipschitz’s condition on the set U, if
| J u J v J ' v , u v | L | u v |2 / 2, (11)
follows that
1
J u J v J ' v , u v J ' v t u v J ' v u v dt.
0
34
2) If J i u , i I is a certain family of the convex functions on
the convex set U, then function J u sup J i u , i I is convex on
set U .
In fact, by the determination of the upper boundary index i* I
and the number 0, such that J (u ) J i* (u ),
i* i* , I are found. Thence we have
J (u ) J i* (u ) 1 J i* (v ), u u 1 v.
m m m
J i u i i J ui , i 0, i 1.
i 1 i 1 i 1
J u f F u f F u 1 F v
f F u 1 f F v J u 1 J v .
f F u 1 F v f t1 1 t 2
f t1 1 f t 2 , t1 , t 2 [a, b] ,
M (C ) u E n u U , J (u ) C
u U u E n / u U 0 , g i (u ) 0, i 1, m;
g i (u ) ai , u bi 0, i m 1, s , (3)
U i u i E n /u U 0 , g i (u ) 0 , i 1, m,
U m1 u E n / ai , u bi 0, i m 1, s .
U m 1 u E n Au b u E n ai , u bi 0, i m 1, s
where А is a matrix of the order s m n ; the vectors
ai , i m 1, s are the lines of the matrix А; b bm 1 ,, bs E n s
is affine set.
Now the problem (2), (3) possible write as
m 1
J (u ) inf, u U U i U 0 . (4)
i 0
uU
J (u* ), u u* 0, u U . (5)
40
(refer to section "Convexity criteria of the even functions"), we have
J (u ) J (u* ) J (u* ), u u* 0 under all u U . This implies
that J (u* ) J (u ), u U . Consequently, the point u* U * .
Theorem is proved.
Consequence. If J (u ) C 1 (U ) , U is a convex set, U *
and the point u* U * , u* int U , then the equality J (u* ) 0
necessary is executed.
In fact, if u* int U , then the number 0 0 such that for any
e E n point u u* eU under all , 0 . is found. Then
from (5) follows J (u* ), e 0, e E n , under all , 0 .
Thence we have J (u* ) 0 .
Formula (5) as precondition of the optimality for nonlinear
programming problem and as necessary and sufficiently condition of
the optimality for convex programming problem will find applying
in the following lectures.
The point projection on set. As application to theorems 1, 2 we
consider the point projection on convex closed set.
Definition 2. Let U be certain set from E n , but the point E n .
The point w U is called the projection of the point E n on set U
, if norm w inf u , and it is denoted by w PU ( )
uU
Example 1. Let the set be
U S (u0 , R) u E n / u u0 R ,
but the point E n . The point w PU ( ) u 0 R( u 0 ) u 0 ,
that follows from geometric interpretation.
Example 2. If the set
U u E n / c, u ,
41
w PU ( ) c, c c . Norm of the second composed
2
From inequality (6) and (8) follows the relationship (7). Theorem
is proved.
42
Example 3. Let the set U u E n / u 0 , the point E n .
Then projection w PU (v) (v1 ,, vn ), where vi max(0, vi ),
i 1, n . To prove.
Example 4. Let the set
U u (u1 ,, u n ) E n / i u i i , i 1, n –
n-dimensional parallelepiped, but the point E n . Then components
of the vector w PU ( ) are defined by formula
i , если vi i
wi i , если vi i
vi , если i vi i , i 1, n.
n
In fact, from formula (6) follows that sum (w v )(u
i 1
i i i wi ) 0
43
Lecture 6, 7
SEPARABILITY OF THE
CONVEX SETS
v U .
Theorem 2. If convex sets A and B from E n have a no points
in common, that hyperplane c, u separating set А and В, as
well as their closures A and B , in the event of point y A B ,
45
number c, y exists.
Proof. We denote U A B . As it is proved earlier (refer to
theorem 1 from lecture 3), set U is convex. Since the intersection
A B , that 0 U . Then with consideration of theorem 1
hyperplane c, u , c, 0 0, separating set U and the
point 0, i.e. c, u 0, u U exists. Thence with consideration of
u a b U , a A, b B , we get c, a c, b, a A,
b B . Finally, hyperplane c, u , where the number
satisfies to inequality inf c, a sup c, b separates the sets А
a A bB
c0 c1 cm 0; (3)
0 1 m 0; (4)
exist.
Proof. We enter the set A A0 A1 Am - a direct product of
the sets Ai , i 0, m with elements a (a0 , a1 ,, a m ) A , where
ai Ai , i 0, m; . We notice that set A E , where E E n ( m1) . It
is easy to show, that set А is convex. In fact, if a1
( a01 , a11 , , a1m ) A, a 2 ( a 02 , a12 , , a m2 ) A , then under all ,
0,1 , the point a a 1 (1 ) a 2 ( a 0 , a1 , , a m ) A ,
since ai ai (1 )ai Ai , i 0, m . We enter the diagonal set
1 2
B b (b0 , b1 ,, bm ) E b0 b1 bm b, bi E n , i 0, m .
Set B is convex. It is not difficult to make sure, the intersection
A0 A1 Am , if and only if, when intersection
A B .
Finally, under performing the condition of theorem the sets А and
B are convex and A B , i.e all conditions of theorem 2 are
47
executed. Consequently, the vector c (c0 , c1 , , c m ) E , c 0
such that c, a c, b, a A, b B exists. The inequality can be
written in the manner of
m m m
i 0
ci , ai ci , bi
i 0
c , b
i 0
i , ai A, bi E n . (5)
m
As follows from inequality (5) linear function J (b) c ,b
i 0
i ,
c ,a
i 0
i i 0, ai Ai , i 0, m . We fix the vectors ai ai Ai ,
m
i k , then ck , ak ci , ai const, a k Ak . Finally, the value
i 0
ck , a k inf ck , a k k , a k Ak , k 1, m . We denote
ak Ak
m
0 ( 1 2 m ) . Then c0 , a0 inf
ak Ak
c
k 1
k , ak
m
k 0 .
k 1
Thereby, values ci , ai i a i Ai , i 0, m ,
0 1 m 0 . Faithfulness of the relationships (2), (4) is
proved. Theorem is proved.
Convex cones. Theorems 1 - 5 are formulated and proved for
convex sets from E n . Separability theorems are often applied in the
extreme problem theories for events when convex sets are convex
cones.
48
Definition 2. Set K from E n is called by cone with top in zero,
if it together with any its point u K contains the points u K
under all 0 . If the set K is convex, that it is called the convex
cone, if K is closed, that it is called by closed cone, if K is open,
that it is called by open cone.
Example 1. The set K u E n / u 0 -a convex closed cone.
In fact if u K , that point u K under all 0, u 0 .
Example 2. The set K u E n / a, u 0 - a convex closed
cone. Since for any 0 , scalar product a, u a, u 0,
u K , consequently, vector u K .
Example 3. The set K u E n / a, u 0 - an open convex
cone, K u E / a, u 0 - closed convex cone, K u E n E
n
E - a convex cone.
n
K * c E n / c, u 0, u K . (6)
49
dual cones K * c E n / c , 0, K * c E n / c
, 0, K * 0 accordingly.
Theorem 5. If the intersection of the inempty convex with vertex
in zero cones K 0 , K1 ,, K m - empty set, that necessary the vectors
ci K i* , i 0, m not all are equal to zero and such that
c0 c1 cm 0 exist.
Proof. Since all condition of theorem 4 are executed, that the
relationship (2) - (4) are faithful. We notice, that K i , i 0, m -
cones, then from ai K i , i 0, m follows that vectors
ai K i , i 0, m under all 0 . Then inequality (2) is written
so: ci , ai i , a i K i , i 0, m , for any 0 . Thence we
have ci , ai i / , i 0, m . In particular, when , we
get ci , ai 0 , i 0, m, ai Ai . On the other hand,
i inf ci , ai 0 , i 0, m . From condition ci , ai 0,
ai A
50
are executed. Since sum c0 c1 cm 0 and ci , w 0,
i 0, m , then from equality c0 c1 cm , w 0 follows that
ci , w 0, i 0, m . By condition of the theorem not all
ci , i 0, m are equal to zero, consequently, at least, two vectors
ci , c j , i j differenced from zero exist. Since all cones are open,
except can be one, that in particular, it is possible to suppose that
cone K i - an open set. Then K i contains the set o( w, 2 )
u K i / u w 2 , 0. In particular, the point u w
w ci / ci K i . Since ci , u 0, u K i , ci K i* , that we have
c , w c / c c , w c c 0,
i i i i i i
where 0,
ci 0 . It is impossible. Consequently, the intersection
K 0 K1 K m . Necessity follows from theorem 5 directly.
Theorem is proved.
51
Lecture 8
LAGRANGE’S FUNCTION.
SADDLE POINT
u U u E n u U 0 , g i (u ) 0, i 1, m;
g i (u ) 0, i m 1, s , (2)
L u* , L u* , * L u , * , u U 0 , 0 . (4)
L u* , * L u , * , u U 0 ; (5)
Proof. Necessity. Let pair u* , * U 0 0 - a saddle point.
We show, that the condition (5). (6) are executed. Since pair
u* , * U 0 0 - saddle point to Lagrange’s function (3), that
inequalities (4) are executed. Then from the right inequality follows
the condition (5). It is enough to prove fairness of the equality (6).
The left inequality from (4) is written so:
s s
J (u* ) i g i (u* ) J (u* ) *i g i (u* ) ,
i 1 i 1
(1 ,, s ) 0 E s .
Consequently, the inequality exists
s
(
i 1
*
i i ) g i (u* ) 0, 0 . (7)
53
At first, we show that u* U , i.e. g i (u* ) 0, i 1, m,
g i (u* ) 0, i m 1, s . It is easy to make sure in that vector
i *i , i 1, s, i j ,
(1 , , s ) (8)
j j 1, for j of 1 j m,
*
belongs to the set E s / 1 0,, m 0 .
Substituting value 0 from (8) to the inequality (7) we get
(1) g i (u* ) 0 . Thence follows that g i (u* ) 0, i 1, m (since it
is faithful for any j of 1 j m ). Similarly the vector
i *i , i 1, s, i j ,
(1 ,, s )
j j g j (u* ), для j из m 1 j s,
*
*i , i 1, s, i j ,
(1 ,, s ) i
j 0, для j из 1 j m.
s
J (u* , * ) J (u* ) g (u ) J (u ) , since the point
i 1
*
i i * * u* U .
Now inequality (5) is written so:
s
J (u* ) J (u ) *i g i (u ), u U 0 . (9)
i 1
56
s
L (u,λ ) 0 J (u ) i g i (u* ), u U 0 ,
i 1 (11)
0 (1 ,, s ) E s 1 0 0,, m 0.
57
It can turn out to be that source problem J (u ) inf, u U ,
U * Ø , has a solution, i.e. the point u* U * exists, however
Lagrange’s function for the given problem has not saddle point.
Example. Let function J (u ) u 1 , but set U u E 1 / 0
0 u 1; (u 1) 2 0 . The function g (u ) (u 1) 2 , set U 0
uE 1
/ 0 u 1, moreover functions J (u ) and g (u ) are
convex on set U 0 . Since set U consists of single element, i.e.
U {1} , that set U * {1} . Consequently, the point u* 1 .
Lagrange’s function L u , (u 1) (u 1) 2 , 0, u U 0
for the problem has not saddle point. In fact, from formula (4)
follows
L u* , 0 L u* , * 0 L u, * (u 1) * (u 1) 2 ,
0, 0 u 1.
u U u E n u U 0 , g i (u ) 0, i 1, m , (13)
i 0, i 1, k , 1 k 1 all restrictions g j (u ) 0, j 1, m,
are regular. In fact,
k k
g j (u ) g j i u i i g j (u i ) 0, j 1, m.
i 1 i 1
59
Lectures 9, 10
KUHN-TUCKER’S THEOREM
A a ( a 0 , a1 , , a m ) E m 1 a 0 J (u ),
a1 g1 (u ) , am g m (u ), u U 0 ; (15)
B b (b0 , b1 , , bm ) E m 1 b0 J * , b1 0 , bm 0 , (16)
a (1 )d a g (u ) (1 ) g (v ) g i (u ),
i i i i i
i 1, m ,
with consideration of function convexity g i (u ), i 1, m on set U 0 .
Finally, the point a (a0 , a1 , , a m ), a0 J (u ), ai
g i (u ), i 1, m, moreover u U 0 . This implies that point
a A under all , 0 1 . Consequently set А is convex. By
the similar way it is possible to prove the convexity of the set B.
c) Since sets А and В are convex and A B Ø , that, by
theorem 2 (the lecture 6), hyperplane c, w , where normal
m 1
vector c (*0 , 1* , *m ) E m 1 , c 0, w E , separating sets
А and В, as well as its closures A A , B b (b0 ,
b1 , , bm ) E m 1
b0 J * , b1 0 , bm 0 exists. Consequently,
inequality are executed
m m
c, b *i bi c, a *i ai , a A, b B . (17)
i 0 i 0
61
We notice that if the point u* U * , that J (u* ) J * ,
g i (u* ) 0, i 1, m , consequently, the vector y ( J * ,0,,0)
A B and the value c, y *0 J * .
Now inequalities (17) are written as
m m
*0 b0 *i bi *0 J * *0 a 0 *i ai ,
i 0 i 0
(18)
a A, b B .
i 1
m
*i g i (u* ) 0, i 1, m , that J * *i g i (u ) L u* , *
i 1
L u , * , L u , * , u U 0 . From the inequality and from
equality g i (u* ) 0, i 1, m, u* U , where *i 0, i 1, m
*
i
follows that pair u* , * U 0 0 - saddle point to Lagrange’s
function. Theorem is proved.
2-nd case. Now we consider the next problem of the convex
programming:
J (u ) inf , (19)
u U u E n u j 0, j I , g i (u ) ai , u bi 0, i 1, m;
g i (u ) ai , u bi 0, i m 1, s , (20)
(1 , , s ) 0 E s 1 0, , m 0.
We notice that, if function J (u ) cu is linear, the task (19),
(20) is called the general problem of the linear programming.
We suppose, that set U * . It is appeared that for task of the
convex programming (19), (20) Lagrange’s function (21) always has
saddle point without some additional requirements on convex set U .
For proving of this it is necessary the following lemma.
n
Lemma. If a1 , , a p - a finite number of the vectors from E
and the numbers i 0, i 1, p , that set
p
Q a E n a i ai , i 0, i 1, p (22)
i 1
Is convex closed cone.
Proof. We show, that set Q is a cone. In fact, if a Q and
p p
0 - an arbitrary number, that a i ai i ai , where
i 1 i 1
a (1 )b i ai (1 ) i ai
p
i 1
p p
( i (1 ) i )ai i ai , i i (1 ) i 0,
i 1 i 1
p
i 1, m , where 0,1, b i ai , i 0, i 1, p .
i 1
64
Thence we have a (1 )b Q under all , 0 1 .
Consequently, set Q - a convex cone.
We show, that Q is convex closed cone. We prove this by the
method of mathematical induction. For values p 1 we have
Q a E n a a1 , 0 - a half-line. Therefore, Q - closed
p
set. We suppose, that set Q p 1 a E n a i ai , i 0 is
i 1
closed. We prove, that set Q Q p 1 a p , 0 is closed. Let
c E n is limiting point of the set Q. Consequently, the sequence
cm Q exists, moreover cm c, m . The sequence cm
is represented in the manner of cm bm m a p , where
bm Q p 1 , m - a numeric sequences. As far as set Q p 1 is
closed, that bm b, m , under b Q p 1 . It is remains to
prove that numeric sequence m is bounded. We suppose opposite
i.e. mk under k . Since bmk / mk c mk / mk a p and
bmk
/ mk Q p 1 the sequence {c mk } is bounded, so far as
c mk c, k , then under k we have a p Q p 1 . Since
set Q p 1 Q (in the event of 0, Q p 1 Q ), that vector
a p Q . It is impossible. Therefore, the sequence m is
bounded. Consequently, vector c b a p Q , where
mk , 0 under k . Theorem is proved.
We formulate and prove Farkas’s theorem having important
significant in theories of the convex analysis previously than prove
the theorem about existence of saddle point to Lagrange’s function
for task (19), (20).
Farkas’ theorem. If cone K with vertex in zero is determined by
inequalities
65
K e E n ci , e 0, i 1, m; ci , e 0, i m 1, p;
ci , e 0, i p 1, s, ci E n , i 1, s , (23)
the dual to it cone K * c E n c, e 0, e K has the type
s
K * c E n c i c i , 1 0, , p 0 . (24)
i 1
s
Q c E n c i ci , 1 0, , p 0 (25)
i 1
complies with K * c E n c, e 0, e K . We represent
i i i , i 0, i 0, i p 1, s . Then set Q from (25) is
written as
p s s
Q с E n c i (ci ) i (ci ) i (ci ),
i 1 i p 1 i p 1
i 0, i 1, p, i 0, i p 1, s, i 0, i p 1, s . (26)
As follows from expressions (26) and (22) the set Q is convex
closed cone generated by the vectors c1 , , c p , c p 1 ,....,
cs , c p 1 , , cs (refer to the lemma).
s
We show, that Q K * . In fact, if c Q , i.e. c c ,
i 1
i i
66
s
1 0,, p 0 , that product c, e i ci , e 0 with
i 1
K * c E n c, e 0, e K .
This implies that Q K * .
We show, that K * Q . We suppose opposite, i.e. vector a K ,
*
1 0,, p 0 . (27)
0, i j , i 1, s
( 1 , , s ) i
j t , t 0, j из 1 j p.
i 0, i j , i 1, s,
( 1 ,, s )
j t c j , d , t 0, j из p 1 j s.
67
2
From inequality (27) we have t c j , d d , a . We divide into
t 0 and, directing t , we get c j , d 0, j p 1, s . Finally,
for the vector d E inequalities c j , d 0,
n
j 1, p, c j , d 0,
j p 1, s are executed. Then from (23) follows that d K ,
where K - a closure of the set K .
Since the vector a K * , that the inequality a , e 0, e K
exists. Thence, in particular, for e d K follows a, d 0 .
However from (27) under i 0, i 1, s , we have a, d 0 . We
have got the contradiction. Consequently, vector a K * belongs to
the set Q. From inclusions Q K * , K * Q follows that K * Q .
Theorem is proved.
Now we formulate theorem about existence of the saddle point to
Lagrange’s function (21) for problem of the convex programming
(19), (20).
Theorem 2. If J (u ) is a convex function on convex set
U 0 , J (u ) C 1 (U 0 ) and set U * for problems (19), (20), then
for each point u * U * necessary the Lagrangian multipliers
* (1* , , *s ) 0 E s / 1 0, , m 0 exist, such
that pair u* , U 0 0 forms a saddle point to Lagrange’s
*
u *j e j 0, j I ; ai , u* e bi 0, i 1, m ;
(28)
ai , u* e bi 0, i m 1, p, , 0 0 .
If the sets of the indexes I 1 j u *j 0, j I , I 2 j ai ,
u * bi , 1 i m are entered, that conditions (28) are written so:
e j 0, j I1 ; ai , e 0, i I 2 ;
ai , e 0, i m 1, s. (29)
K e E n e j e j , e 0, j I 1 ; ai , e 0, i I 2 ,
ai , e 0, i m 1, s , (30)
69
s
J (u* ) *j e j *i ai a *
i i (31)
jI1 iI 2 i m 1
exist.
Let *i 0 for values i 1,2, , m \ I 2 . Then expression (31)
is written as
s
J (u* ) *i ai *j e j . (32)
i 1 jI1
s
L u , * L u* , * J (u ) J (u* ) *i ai , u u* . (33)
i 1
s
L u , * L u* , * J (u* ) *i a i , u u*
i 1
jI1
*
j e , u u* (u j u *j )
j
jI1
*
j u
jI1
*
j j 0,
u U u E n u U 0 , g i (u ) 0, i 1, m; g i (u ) ai , u
bi 0 , i m 1, p,
g i (u ) ai , u bi 0 , i p 1, s , (35)
that the pair u* , * U 0 0 forms a saddle point to Lagrange’s
function (36) on set U 0 0 .
Proof of theorem requires of using more fine separability
theorems of the convex set, than theorems 1 - 6 (the lectures 6, 7).
For studying more general separability theorems of the convex set
and other Kuhn-Tucker’s theorems is recommended the book:
Rocafellar R. Convex analysis. Moskow, 1973.
We note the following:
10. The theorems 1-3 give the sufficient conditions of existence
of the saddle point for the convex programming problem.
Example 1. Let function be J (u) 1 u , but set U u
71
E 1 0 u 1; (1 u ) 2 0. For this example the functions
J (u ), g (u ) (1 u ) 2 are convex on convex set U 0 u E 1 0
u 1 . The set U 1 , consequently U * 1, set U * 1. The
point u* 1 and the value J (u * ) 0 are solution of the problem
J (u ) inf, u U . Lagrange’s function L u, (1 u ) (1
1 u ) 2 , u U 0 , 0 . The pair u* 1, * 0 - saddle point to
Lagrange’s function, since L u , 0 L u, (1 u)
*
* *
* (1 u ) 2 , 0 u 1, * g (u * ) 0 .
We notice, that for the example neither Sleyter’s theorem from
theorem 1, nor condition of theorem 3 is not executed.
2°. In general event Lagrangian multipliers for the point u* U *
are defined ambiguous. In specified example the pair u* 1, * 0
is the saddle point to Lagrange’s function under any 0 .
*
72
s
30. To form Lagrange’s function L u , J (u ) g (u ) i i
i 1
0 E s 1 0, , m 0 .
40. To find a saddle point u* , * U 0 0 to Lagrange’s
function from the conditions (the main lemma from lecture 8):
L u* , * L u , * , u U 0 ,
(37)
*i g i (u* ) 0, i 1, s, 1* 0, , *m 0.
s
where L u u * , * J (u * ) g (u ) .
i 1
*
i i * Now condition (37) are
written as
L u u* , * , u u* 0 , u U 0 ,
(38)
*i g i (u* ) 0, i 1, s, 1* 0, , *m 0.
b) If we except J (u ) C 1 (U 0 ), g i (u ) C 1 (U 0 ), i 1, m ,
set U 0 E n , then according to optimality criteria (lecture 5) the
conditions (38) possible present in the manner of
73
L u u* , * 0, *i g i (u* ) 0, i 1, s,
1* 0, , *m 0 . (39)
74
Chapter II. NONLINEAR
PROGRAMMING
75
Lectures 11, 12.
STATEMENT OF THE PROBLEM.
NECESSARY CONDITIONS OF THE
OPTIMALITY
u U u E n / u U 0 , g i (u ) 0, i 1, m; g i (u ) 0
0, i m 1, s , (2)
U i u E n g i (u ) 0 , i 1, m,
(3)
U m 1 u E n g i (u ) 0 , i m 1, s ,
ensemble U possible to present in the manner of
U i u E n g i (u ) 0 , i 1, m,
(3)
U m 1 u E n g i (u ) 0 , i m 1, s ,
Now problem (1), (2) can be written in such form:
J (u ) inf, u U
L u (u * , * ), u u*
s
*0 J (u* ) *i g i(u* ), u u* 0, u U 0 , (7)
i 1
U u E 1 0 u 1; g (u ) (u 1) 2 0
Generalized Lagrange’s function for the problem has the form
L (u , ) 0 J (u ) (u 1) 2 , u U 0 u E 1 0 u 1, 0 0,
0 . Ensemble U 1 , consequently, ensemble U * 1
moreover J (u* ) 0, g (u* ) 0 . The condition (7) is written as
(*0 J (u* ) * g (u* )) (u u* ) 0, u U 0 . Thence we have
(*0 )(u 1) 0, u U 0 , or ( *0 )(1 u ) 0, 0 u 1 . Since
expression (1 u ) 0 under 0 u 1 , then for executing the
inequality necessary that value *0 0 , so as *0 0 by the
condition (6). Finally, the source problem is degenerate. The
78
condition (8) is executed for any 0 . Thereby, all conditions of
*
U (u1 , u2 ) E 2 g (u1 ) u1 0 .
Ensemble U 0 E 2 . Generalized Lagrange’s function
is executed.
We denote through K y K y (u* ) the ensemble of all directions
of the function decrease J (u ) in the point u* . Finally, the ensemble
K y e E n e e , J (u* e ) J (u* ),
0 0 . (11)
K y e E n J (u* ), e 0 (12)
- open convex cone. By Farkas’s theorem the dual cone to cone (12)
is defined on formula
80
K y* c y E n c y *0 J (u* ), *0 0 . (13)
K bi e E n e e , u* e U i ,
, 0 0 . (14)
U i u E n g i (u ) 0 , i 1, m
K bi e E n g i ( u* e ) 0, , 0 0 . (15)
K bi e E n g i (u* ), e 0 , i 1, m. (16)
K b*i ci E n ci *i g i (u* ), *i 0 , i 1, m. (17)
K m 1 e E n u* e r ( ) U m 1 , , 0 0 ;
r (0) 0, r ( ) / 0 при 0.
According to formula (3) the ensemble
U m1 u E n g i (u ) 0, i m 1, s ,
consequently,
82
K m1 e E n g i (u* e r ( )) 0, i m 1, s,
, 0 0 , (18)
K m 1 e E n g i (u* ), e 0, i m 1, s . (19)
s
K m* 1 c E n c *i g i (u* ) . (20)
i m 1
Finally, we define K 0 - an ensemble internal directions of the
convex ensemble U 0 in the point u* U U 0 . It is easy to make
sure in that, if u* intU 0 , that K 0 E n , consequently, K * 0.
If u* pU 0 , that
K 0 e E n e (u u* ), u int U 0 , 0 (21)
K 0 с 0 E n с 0 , u с 0 , u * при всех u U 0 . (22)
K y K 0 K 1 K m K m 1 . (23)
e K y K 0 K 1 K m K m 1 .
u* e u* e r ( ) U 0 U 1 U m U m 1 U .
84
impossible, since the point u* U U 0 is solution of the problem
(1), (2). The lemma is proved.
Previously than transfer to proof of the theorem, we note the
following:
1) If J (u* ) 0, that under *0 1, 1* *2 ... *s 0 all
condition of the theorem 1 are fulfilled. In fact, norm * 1 0 is
a scalar product
L u (u * , * ), u u * *0 J (u * ), u u * 0,
u U 0 , *i g i (u* ) 0, i 1, s.
s
L u ( u * , * ), u u *
i m 1
*
i g i ( u * ) , u u * 0 ,
u U 0 , *i g i ( u * ) 0 , i 1, s ,
since g i (u* ), i m 1, s .
From items 1 - 3 results that theorem 1 should prove for event,
85
when J (u* ) 0, g i (u* ) 0 under all i 1, m , and vectors
g (u ), i m 1, s are linearly independent.
i *
Proof of the theorem 1. By the data of the theorem the ensemble
u* U . Then, as follows from proved lemma, the intersection of the
convex cones
K y K 0 K 1 K m K m 1 (24)
s
c0 , u u 0 *0 J (u* ) *i g i (u* ), u u 0
i 1
L u (u* , * ), u u 0 0, u U 0 .
86
Thence follows fairness of the correlations (7). The condition (6)
follows from that not all c y , c 0 , c1 , , c m , are equal to zero. We
notice, that if g i (u* ) 0 for a certain i from 1 i m , that cone
Ki E n , consequently, K i* 0. It means that
ci *i g i (u* ) 0, g i (u* ) 0 . Thence follows that *i 0 .
Thereby, the products i g i (u* ) 0, i 1, s , i.e. the condition (8)
*
87
Lecture 13
SOLUTION ALGORITHM OF THE
NONLINEAR PROGRAMMING PROBLEM
u U u E n u U 0 , g i (u ) 0, i 1, m;
g i (u ) 0, i m 1, s , (2)
where functions J (u ) C (U 01 ), g i (u ) C (U 01 ), i 1, s, U 01
1 1
88
s
L (u, ) 0 J (u ) i g i (u ), u U 0 ;
i 1
(0 , 1 , , s ) 0 E s 1 / 0 0, 1 0, , m 0.
L u (u * , * ), u u *
s
*0 J (u* ) *i g i(u* ), u u* 0, u U 0 , (4)
i 1
0 E s / 1 0, , m 0,
u U u E n g i (u ) 0, g i (u ) C 1 ( E n ), i 1, m . (8)
The problem (7), (8) is a particular event of the problem (1), (2).
The ensemble U 0 E n . The point in the problem (7), (8) is named
by normal point of the minimum, if the vectors g i (u* ), i 1, m are
linearly independent. We notice, that if u* U - normal point, that
problem (7), (8) are nondegenerate. In fact, for problem (7), (8) the
equality
m
*0 J (u * ) *i g i (u * ) 0 (9)
i 1
exists.
If *0 0 , that on the strength of linearly independence of the
vectors g i (u* ), i 1, m we would get i 0, i 1, m . Then the
*
90
Theorem. Let functions J (u ), g i (u ), i 1, m are determined and
twice continuously differentiable in neighborhoods of the point u* U .
In order the point u* U to be a point of the local minimum J (u ) on
ensemble U , i.e. J (u* ) J (u ), u, u o(u* , ) U sufficiently
2 L (u* , * )
that quadratic form y y to be positively determined on
u 2
hyperplane
y 0. (10)
u
g u / u g m u* / u n y n
m * 1
2 L (u* , * )
y y
y 2
2 L (u* , * ) / u12 2 L (u* , * ) / u1un
2
L (u* , * ) / u2 u1 2 L (u* , * ) / u2 un
( y1 , , yn )
2 L (u , * ) / u u 2 L (u* , * ) / un2
* n 1
y1
y
2 0, y 0,
y
n
91
function minimum J (u ) on ensemble U. We notice, that u* - a
normal point and it is determined from the conditions
L u (u* , * ) 0, *i g i (u* ) 0, i 1, m , i.e. g i (u* ) 0, i 1, m .
Finally, pair (u* , * ) is known.
Since functions
2 L (u* , * )
min y y, y V
u 2
We enter the ensemble
A u E n u u* y, y 1, (g u* / u )* y 0 , (11)
1
L (u, * ) L (u* , * ) L u (u* , * ), u u* (u u* )
2
L uu (u* , )(u u* ) o u u*
* 2
2
2
o( 2 )
92
o( 2 )
2 2 2 , 0 1 , (13)
2 4
B u E n u u* , g i (u ) 0, i 1, m U . (14)
u A : g i (u ) g i (u * ) g i (u ) g i (u * ), u u *
1 2
(u u* ) g iuu (u* ) (u u* ) oi u u*
2
1
2 y g iuu (u * ) y oi ( ), y 1; i 1, m; (15)
2
1
u~ B : g i (u~ ) g i (u* ) 0 g i(u* ), u~ u* (u~ u* )
2
g iuu (u * )(u~ u * ) oi u~ u * g i (u* ), u~ u u u *
2
1 ~
(u u u u * ) g iuu (u* )(u~ u u u * ) oi ( 2 )
2
1 1
g i (u * ), u~ u (u~ u ) g iuu (u * )(u~ u * ) (u u * )
2 2
~
g (u )(u u ) (u u ) g (u )(u u )
iuu * * iuu * *
oi ( 2 ), i 1, m; (16)
93
From (11), (14) and (15), (16) follows that norm u~ u K ,
2
u~ B , u A .
Since the function L u (u , * ) continuously differentiable by u
in neighborhoods of the point u* , and derivative L u (u* , * ) 0 ,
that difference L u (u, ) L u (u* , ) L u (u , * ) L uu (u* , )(×
* * *
×(u u* ) o u u* in neighborhoods of the point u* . This
implies that, in particular, u A , norm L (u , * ) K 1 , if
0 1 , 1 0 - sufficiently small number.
For the points u A , u~ B the difference L (u ~, * ) L
1
L (u, * ) L u (u , * ), u~ u (u~ u )L uu (u* , * )(u~ u )
2
o u~ u ,
*
2
consequently, norm L (u~, * ) L (u, * )
u
g (u* ) 0 : u1*
2 * 2
2 1,
1) * 3 / 2, u1* 5 / 6 , u 2* 1 / 6 ;
2) * 3 / 2, u1* 5 / 6 , u 2* 1 / 6 ;
3) * 3 / 2, u1* 5 / 6 , u 2* 1 / 6 ;
4) * 3 / 2, u1* 5 / 6 , u 2* 1 / 6 .
2 2* 5 2u *
L uu (u* , * ) , g u (u* ) g (u* ) 1* .
2u
5 2 2* 2
I i i m 1, s, g i (u* ) 0, 1 i m .
If the problem (1), (2) is nondegenerate, that vectors g i (u* ), i I
are linearly independent. According to specified above theorem from item
40 the point u* U is the point of the local minimum to functions J (u )
on U , if quadratic form y L uu (u* , ) y 0, y 0 on hyperplane
*
g i u* / u * y 0, i I .
96
Lecture 14
DUALITY THEORY
u U u E n u U 0 , g i (u ) 0, i 1, m;
g i (u ) 0, i m 1, s , (2)
The main task. We enter the function [refer to. the formula (3)]:
97
X (u ) sup L (u, ), u U 0 . (4)
0
J (u ) при всех u U ,
X (u ) (5)
при всех u U 0 \ U .
In fact, if u U , that g i (u ) 0, i 1, m, g i (u ) 0,
i m 1, s , consequently,
s
X (u ) sup J (u ) i g i (u )
0 i 1
m
sup J (u ) i g i (u ) J (u ) ,
0 i 1
m
since 1 0,, m 0, g (u ) 0
i 1
i i under all u U , moreover
X (u ) inf, u U 0 , (6)
98
U * u* U min J (u ) J (u* ) J *
uU
( ) sup, 0 , (8)
0
* J * X (u ), u U 0 , 0 . (9)
99
on the strength of correlations (4). From correlations (10)
transferring to lower bound on u we get * inf X (u ) J * .
uU 0
It means that pair u* , * U * * - saddle point. Moreover
ensemble U * belongs to the ensemble of the saddle points to
*
100
Lagrange’s function (3). We show that correlations (11) are fulfilled.
As follows from determination of the saddle point which has the
form (13) inequality L (u* , ) L (u* , * ), 0 faithfully.
Consequently,
is equitable.
2°. If u* , * , a* , b * U 0 0 are saddle points of the
Lagrange’s function (3) on U 0 0 , that u* , b* , a* , * - also
saddle points to function (3) on U 0 0 , moreover L (u* , b* ) L
101
L (a* , * ) L (u* , * ) L (a* , b* ) (* ) X (u* ) J * * .
However inverse statement, i.e. that from L (u* , * ) L (a, b)
follows a, b U 0 0 - saddle point in general event untrue.
3°. Dual problem (8) it is possible to write as
( ) inf, 0 . (16)
J (u ) c, u inf,
(17)
u U u E n u 0, Au b 0 ,
where c E n , b E m are the vectors; А is the matrix of the order
m n ; the ensemble
U 0 u E n u u1 0, , u n 0 0
L (u, ) c, u , Au b c A* , u b, ,
(18)
u U 0 , (1 , , m ) 0 E m 1 0, , m 0 .
As follows from formulas (17), (18), function
102
b, , если c A* 0
inf L (u, )
uU 0
, если ( c A* ) i 0, 0 .
We notice, that under c A* 0 lower bound which is equal
to b, is got under u 0 U 0 . If ( c A* ) i 0 , that it is
possible to choose ui , but all u j 0, j 1, n, i j ,
herewith , * 0 . Finally, the dual task to task (18)
has the form
b, inf,
(19)
E m 0, c A* 0.
The dual task to task (19) complies with source task (17).
By introducing of the additional variables u n i 0, i 1, m ,
optimization problem (17) can be written in the following type:
J (u ) inf
. (20)
Au i bi u ni 0, u 0, u n i 0, i 1, m
J (u ) c, u inf,
(21)
u U u E n u j 0, j I , Au b 0, A u b 0 ,
where c E n , b E m , b E s are the vectors; A, A are the
matrixes accordingly to the orders m n, s n ; index ensemble
I 1,2,, n. The ensemble
U 0 u u1 ,, u n E n u j 0, j I
Lagrange’s function for task (21) is written as
103
L (u , ) c, u , Au b , A u b
c A* A * , u b, b , ,
u U 0 , ( , ) 0 ( , ) E m E s 0 .
Function
inf L(u, )
uU 0
b, b , , if
c A* A * i 0, i I ;
c A* A * j 0, j I ;
under rest.
b, b , inf; c A* A *
i 0, i I ;
(22)
c A A
* *
j 0, j I ; ( , ) E n E s , 0.
J (u ) c, u inf,
(24)
u U u E n u 0, Au b 0 ,
104
where c E n , b E s are the vectors; А is the matrix of the order
s n ; the ensemble
U 0 u E n u u1 , , u n 0
L (u , ) c, u , Au b c A* , u b, ,
u U 0 , 0 E s .
Function
b, , если c A* 0,
( ) inf L (u, )
uU 0
, если c A* 0.
Then the dual problem to the problem (24) has the form
( ) b, inf; c A* 0, E s . (25)
It is easy to make sure in the dual task to task (25) complies with
task (24). Finally, we note that main and the general task of the linear
programming by the way of introducing some additional variables
are reduced to the canonical tasks of the linear programming [refer to
formulas (20), (23)].
105
Chapter III. LINEAR
PROGRAMMING
Lectures 15, 16
STATEMENT OF THE PROBLEM.
SIMPLEX-METHOD
J (u ) c, u inf, u U u E n u 0, Au b 0 , (1)
where c E n , b E m are the vectors; А is the matrix of the order
m n . Matrix A aij , i 1, m, j 1, n can be represented in the
manner of
106
a1 j
a2 j
2 n
A a , a ,..., a , a , j 1, n .
1 j
...
a
mj
The vectors a j , j 1, n are identified by the condition vectors, but
vector b E m - a vector of the restrictions. Now the equation
Au b can be written in the manner of a1u1 a 2 u 2 ... a n u n b
. Since ensemble U 0 u E n u 0 and U u E n Au b
- affine ensemble which are convex, that (1) is a problem of the
convex programming. We notice that if ensemble
U * u* E n u* U , J (u* ) min c, u ,
uU
that Lagrange’s function for problem (1) always has saddle point,
any point of the local minimum simultaneously is the point of the
global minimum and necessary and sufficiently condition of the
optimality is written as J ' (u* ), u u* c, u u* 0, u U 0 .
We suppose, that ensemble U * . It is necessary to find the point
u* U * and value J * inf J (u ) J (u* ) min J (u ) .
uU uU
Simplex-method. For the first time solution of the problem (1)
was considered on simplex
n
U u E n u 0, ui 1
i 1
so solution method of such linear programming problem was called
by simplex-method. Then method was generalized for event of
ensemble U specified in the problem (1), although initial name of
the method is kept.
Definition 1. The point u U is identified by extreme (or
angular), if it is not represented in the manner of
107
u u 1 1 u 2 , 0 1, u 1 , u 2 U . From given definition
follows that extreme point is not an internal point of any segment
belonging to ensemble U.
Lemma 1. Extreme point u U has not more m positive
coordinates.
Proof. Not derogating generalities, hereinafter we consider that
first components k , k m of the extreme point are positive, since
by the way of recalling the variables always possible to provide the
given condition.
We suppose opposite, i.e. that extreme point u U has m 1
positive coordinates
u u1 0, u 2 0, ..., u m1 0, 0,...,0 .
We compose the matrix A1 (a 1 , a 2 ,..., a m 1 ) of the order
m m 1 from condition vectors corresponding to positive
coordinates of the extreme point. We consider the homogeneous
linear equation A1 z 0 for vector z E m 1 . The equation has a
nonzero solution ~ z 0 . We define n-vector u~ ~
z, ~ z ,0 and
consider two vectors: u 1 u u~, u 2 u u~ , where 0 -
sufficiently small number. We notice that vectors u 1 , u 2 U under
0 1 , where 1 0 is sufficiently small number. In fact,
Au 1 Au A u~ Au A1 ~ z Au b, u 1 u u~ 0 under
sufficiently small 1 0 similarly Au 2 b, u 2 0 . Then the
extreme point u 1 / 2 u 1 1 / 2 u 2 , 1 / 2 . It opposites to the
definition of the extreme point. The lemma is proved.
Lemma 2. Condition vectors corresponding to positive
coordinates of the extreme point are linear independent
Proof. Let u u1 0, u 2 0, ..., u k 0, 0,...,0 U , k m
be an extreme point. We show that vectors a1 , a 2 ,..., a k , k m are
linear independent. We suppose opposite, i.e. that there are the
numbers 1 ,..., k not all equal to zero such that
1a1 2 a 2 ...
108
k ak 0 (the vectors a1 ,..., a k are linear dependent). From
inclusion u U follows that a1u1 a 2 u 2 ... a k u k b,
ui 0, i 1, k . We multiply the first equality on 0 and add
(subtract) from the second equality as a result we get
a 1 u1 1 a 2 u2 2 ... a k u k k b . We denote
by u 1 u1 1 ,..., u k k ,0,...,0 E n , u 2 u1 1 ,...,
k ,0,...,0 E n . There is a number 1 0 such that u 1 U ,
u 2 U under all , 0 1 . Then vector u 1 / 2u 1
1 / 2u 2 U , u 1 U , u 2 U . We obtained the contradiction in
that u U - extreme point. Lemma is proved.
From lemmas 1, 2 follows that:
a) The number of the extreme points of the ensemble U is finite
m
and it does not exceed the sum C
k 1
k
n , where C nk - a combinations
109
problem the feasible vector has exactly m positive coordinates, that
u - an extreme point of the ensemble U.
Proof. Let the feasible vector u u1 0, u m 0, 0,...,0 U
has exactly m positive coordinates. We show, that u is an extreme
point of the ensemble U .
We suppose opposite i.e. there are the points
u , u 2 U , u 1 u 2 , and number , 0 1 such that
1
U u (u1 , u 2 , u 3 , u 4 ) E 4 u j 0, j 1,4;
3u1 u 2 u 3 u 4 3, u1 u 2 2u 3 u 4 1.
3 1 1 1
A a1 , a 2 , a 3 , a 4 , rangA 2.
1 1 2 1
110
By the extreme points of the ensemble U are u 1 1,0,0,0 ,
u 2 0,5 / 3, 4 / 3, 0, u 3 0,1,0,2 . Here u 2 , u 3 are nondegenerate
extreme points, u 1 is degenerate extreme point. Since there is
feasible vector u 1 , number of the positive coordinates less than
rangA , the problem of the linear programming (1) is not
nondegenerate. The number of the extreme points are equal to 3, that
does not exceed the amount c14 c42 10 ; the vectors corresponding
to the positive coordinates of the extreme point u 2 a 2 ,a 3 ,
u a ,a
3
2 4
are linear independent.
Example 2. Let ensemble be
U u (u1 , u 2 , u 3 , u 4 ) E 4 u j 0, j 1,4;
3u1 u 2 u 3 u 4 3, u1 u 2 2u 3 u 4 1.
k 1
k 1, u1 , u 2 ,..., u s - an extreme points of the ensemble U .
k 1
k 1 . It is remains to show, the points u1 , u 2 ,..., u s1 are extreme
112
s1
on the strength of proved a k v k , k 0, k 1, s1 ,
k 1
s1 s2 s2
k 1
k 1, b k w , k 0,
k 1
k
k 1, s 2 ,
k 1
k 1, where
s1 s1
the point u a 1 b k v k 1 k wk , moreover
k 1 k 1
k k 0, k 1, s1 , k 1 k 0, k 1, s 2 ,
s1 s2
k k 1 1 . Lemma is proved.
k 1 k 1
s
i 0, i 1 , where u1 ,..., u s - extreme points of the ensemble
i 1
s s
U . The value J (u* ) c, u* i c, u i i J i , where
i 1 i 1
J i J (u ) c, u , i 1, s.
i i
113
s
Let J 0 min J i J (u i0 ) . Then J (u* ) J 0 i J 0
1i s
i 1
u * u1* , u 2* ,..., u m* ,0,...,0 , u1* 0, u 2* 0,..., u m* 0 .
115
under all 0 0 . We present the vector l E n in the manner of
l lБ , lH , where lБ l1 ,..., lm E m , lH lm1 ,..., ln E n m .
From inclusion u* l U follows that u Б* l Б 0, u H* l H
l H 0, Au* l AБ u Б* l Б AH l H b . Since u* U
, that Au* AuБ* b consequently from the last equality we have
AБ lБ AH lH 0 . Thence follows the vector lБ АБ1 АН lH . Since
under sufficiency small 0 inequality u Б* lБ 0 is executed,
the feasible directions in the point u* are defined by the correlations
l H 0, l Б AБ1 AH l H . (3)
u U .
Ba
sis
с с1 … сj … с j0 … сn
AБ Сб b
а1 aj …
a j0 …
an
а1 с1 u11 … u1 j … u1 j0 … u1n
u1*
… … … … … … … … … … …
a i0
ci0 u *
i0
ui0 1 … ui0 j … ui0 j0 … ui0 n 0
… … … … … … … … … … …
ai ci ui* ui1 … uij … uij0 … uin i
… … … … … … … … … … …
m
cm *
um1 … u mj … umj0 … umn
a u m
z z1 … zj … z j0 … zn
z c 0 … zj cj … zj0 cj0 … zn cn
117
The condition vectors a1 , a 2 ,..., a m corresponding to positive
coordinates of the extreme point u* u1* ,..., um* , 0,..., 0 are leaded
in the first column of the table. Matrix АБ a ,..., a .
1 m
zc z c Б
, z c H 0, z c H ,
118
m
Since the values z j c u
i 1
i ij cБ u j cБ АБ1а j , j 1, n , that
z j c j c Б А a c j , j 1, n . Consequently, vector ( z c )
1
Б
j
(0, c Б АБ1 AH c H ) .
Comparing the correlation with optimality criterion (2), we make
sure in the extreme point u* U to be a solution of the problem
necessary and sufficiency that values z j c j 0, j 1, n . Finally,
by sign of the values in the last line of the simplex-table it is possible
to define whether the point u* U a solution of the nondegenerate
problem (1).
119
Lecture 17
DIRECTION CHOICE.
NEW SIMPLEX-TABLE CONSRUCTION.
THE INITIAL EXTREME POINT
CONSRUCTION
J (u ) c, u inf, u U u E n u 0, Au b 0 . (1)
120
Direction choice. Since u U is an extreme point, not
derogating generalities, it is possible to consider that the first m
components of the vector u are positive, i.e. u (u1 , u 2 ,..., u m ,
0,...,0), u i 0, i 1, m . According to formula (2) (the lecture 16),
feasible directions in the point u are defined from correlations
lH 0, lБ АБ1 АН lH . Derivative of function J (u ) on feasible
direction l in the point u is equal to
J u / l c, l c Б , l Б c H , l H c H' c Б' AБ1 AH l H
z c j l j , l j 0, j m 1, n .
n
j
j m 1
l 0 L, l 0 l Б0 , l H0 , l H0 0,...,0, 1, 0,...0,
121
l Б0 AБ1 AH l H0 AБ1a j0 u j0
(u1 j0 , u 2 j0 ,...,u m j0 ) (2)
0 min i min ui uij0 i0 0, i0 I 2 .
iI 2 iI 2
122
J (u ) 0 ( z j0 c j0 ). Thence we have J (u ( 0 )) J (u )
0 ( z j0 c j0 ) 0. Then J (u ( 0 )) J (u ) J (u ) , i.e. in the
extreme point u 0 u value J u less, than in extreme point u .
We note the following:
1) If ensemble of the indexes I 2 , i.e. uij0 0, i 1, m , then
for any 0 the point u U , moreover J (u ( )) J (u )
( z j0 c j0 ) under . In this case source problem
has a no solution. However, if ensemble U is compact, that it isn’t
possible.
2) It can turn out so that 0 i0 i0 , i0 I 2 , i0 I 2 . In this
case vector u 0 has m 1 positive coordinates. It means that
source problem is degenerate. In such events possible appearance
"thread", i.e. through determined number steps once again render in
the extreme point u 0 . There are the different methods from
"recirculations". We recommend the following books on this
questions: Gabasov R., Kirillova F. M. Methods of the optimization.
Minsk: BGU, 1975, Karmanov V.G. Mathematical programming.
M.: Science, 1975; Moiseev N.N., Ivanilov U.P., Stolyarova E.M.
Methods of the optimization. M.: Science, 1978.
To make sure in the extreme point u 0 u U is a solution
of the problem (1), need to build the simplex-table for the extreme
point u with aim to check the optimality criterion (2) (lecture 16) in
the point u .
Construction of the new simplex-table. The simplex-table built
for the extreme point u* U in the previous lecture, in particular for
the extreme point u U is true. In the simplex-table column j0
and line i0 and values i , i I 2 are indicated.
We build the simplex-table for the extreme point u 0 u U
on base of simplex-table for the point u* u . We notice, that
123
u U by the base vectors were condition vectors a1 ,..., a m for the
point corresponding to positive coordinates of the extreme point u .
By the base vectors will be
a1 , a 2 ,..., a i0 1 , a j0 , a i0 1 , a m , (2*)
have
m m uij0 1
a j0 a i uij0 a i0 ui0 j0 , a i0 a i a j0 . (3)
i 1 i 1 ui0 j0 ui0 j0
i i0 i i0
m
a j a i ui j a i0 ui0 j
i 1
i i0
124
m ui j ui j ui0 j j
a i ui j 0 0 a 0, j 1, n. (3*)
ui 0 j 0 ui j
i 1
i i0
0 0
J (u ) c, u inf, [ Au ]i u n i bi , i 1, m, (5)
J (u ) c , u inf,
u U u E n m u 0, A u b . (6)
For problem (6) initial extreme point u 0,..., 0, b1 ,..., bm
Є
ЄEn+m , since condition vectors a n 1 ,..., a n m (the columns of the
unit matrix I m ) corresponding to positive coordinates of the
extreme point u~ , b 0 are linear independent.
2. Danzig’s method (two-phase method). We consider canonical
problem of the linear programming
J (u ) c, u inf, u U u E n u 0, Au b 0 , (7)
where c E n , b E n ; А - a matrix of the order m n . We suppose
that vector b 0 (if bi 0, i 1, m , that being multiplied by the
first line, where bi 0 , always possible to provide to bi 0, i 1, m ).
126
The following canonical problem on the first stage (phase) is
solved:
u
i 1
ni inf,
(8)
Au i u ni bi , i 1, m, u 0, u n i 0, i 1, m.
Au i uni bi , u 0, uni 0, i 1, m ,
where M 0 - feasible large number. For M-problem (9) the initial
extreme point u 0,..., 0, b1 , b2 ,..., bm E n m , b 0 , and it is
solved by the simplex-method. If the source problem (7) has a
solution, that М- problem has such solution: u~* u~1* ,..., u~n* ,0,...,0 ,
127
where components un*i 0, i 1, m . Vector u* u~1* , u~2* ,..., u~n* E n
is the solution of the problem (7).
We notice, that for М-problem values z j c j j M j ,
j 1, n . So in the simplex-table for lines z j , z j c j are conducted
two lines: one for coefficients j , other for j . Index j0 , where
z j0 c j0 max ( z j c j ), z j c j 0 is defined by the value of
j
the coefficients j , j 0 .
128
REFERENCES
129
Appendix 1
130
Theorem 2. If function J (u ) C
2
E n , then in the point u 0 En
the following conditions: J u 0 0, J u 0 0 (necessary condition
of the second order) are executed.
Theorem 3. In order to the point u0 E n be a point of the local
minimum to function J (u ) C
2
E n , necessary and sufficiency
J u0 0, J u0 0 .
1 u1 u 2
6. J (u1 ,u 2 ) , u (u1 , u 2 ) E 2 .
1 u1 u1
2 2
u 22 u 32 2
7. J (u1 ,u 2 ,u 3 ) u1 , u1 0, u 2 0, u 3 0.
4u1 u 2 u 3
8. J (u1 ,u 2 ) u12 u1u 2 u 22 2u1 u 2 , u (u1 ,u 2 ) E 2 .
131
P.1.2. Convex ensembles and convex functions
If int U и J(u) C
2
(U ) , then for convexity of J (u ) on U
necessary and sufficiency that
J''(u) , 0, E n , u U .
132
1) J (u ) J v J' v , u-ν κ|u ν|2 , κ 0, u,ν U; (5)
(6)
J''(u) , 2 , 2 0, u, U .
133
Then function J (u ) ( u ) ( )d is convex on E 1 . To prove.
u
1
u 2 u1
7. Let J (u1 , u 2 ) ( u1 ) ( )d , u2 0, ( ) 0 ,
linear function.
10. Function J (u ) is convex on U , if and only if function
g ( ) J (u ( u )) by one variable , 0 1 is convex
under any u, U . If J (u ) is strictly convex on U , that
g ( ), 0 1 is strictly convex. To prove.
n
11. Function J (u ) determined on convex ensemble U from E is
identified quasiconvex, if J (u (1 ) ) maxJ (u ), J ( ),
u, U and under all , 0 1. Is any convex function
quasiconvex, conversely? To prove that J (u ) is quasiconvex on U if and
only if ensemble M ( ) u U J(u) J (ν) is convex under all
U .
12. To check in each of the following exercises whether the function
J (u ) convex (concave) on given ensemble U, or indicate such points from
U in the neighborhood of which J (u ) is not neither convex, nor concave:
a) J (u ) u16 u 22 u32 u 42 10u1 5u 2 3u 4 20, U E 4 ;
2u1u2
б ) J (u ) e , u E2 ;
134
в ) J (u ) u15 0,5u 32 7u1 u3 6,
U u (u1 , u 2 , u 3 ) E 3 u i 0, i 1,2,3 ;
г) J (u ) 6u12 u 23 6u32
12u1 8u 2 7, U u E 3 / u 0 .
2
13. Let J (u ) Au b Au b, Au b , where A is the
matrix of the order m n, b Е m - the vector. To prove that J (u ) is
convex on E n . If A*A is nondegenerate, that J (u ) is strictly convex on
E n . To find J u , J u .
J(u) 0 ,5 Au,u b, u , where A A is the matrix of
*
14. Let
u12 u 22 (u12 u 22 ) 2 .
19. To find on plane of the parameters , areas, where function
J (u1 ,u 2 ) u1α u 2β , u1 0, u 2 0 is convex (strictly convex) and
concave (strictly concave).
135
u1u2
20. To find on plane E2 areas, where function J (u1 , u 2 ) e is
convex, and areas in which it is concave.
21. Let J i (u ), i 1,m be convex, nonnegative, monotonous
m
increasing on E 1 functions. To prove that function J (u ) J i (u ).
i 1
possesses by these characteristics.
22. Let J (u ) be a function determined and convex on convex
ensemble U. To prove that: a) ensemble Г E / sup
n
,u
uU
136
P.1.3. Convex programming. KUNA-TUCKER’S theorem
J (u ) inf (1)
at condition
u U u E n u U 0 ; g i (u ) 0, i 1, m;
g i (u ) ai , u bi 0, i m 1, p;
g i (u ) ai , u bi 0, i p 1, s , (2)
is executed.
*
Lemma. In order to pair u* , U 0 0 be saddle point of the
Lagrange’s function (3), necessary and sufficiency performing the following
conditions:
137
L (u* ,* ) L (u,λ* ), u U 0 ,
(5)
*i g i (u* ) 0, i 1, s, u* U * U , * 0 ,
at condition
u1 u 2 1, 8u12 u 22 2, u1 0, u 2 0. (7)
at condition
u1 u 2 1, 8u12 u 22 2, u1 0, u 2 0. (9)
u U u u1 , u 2 E 2 /u U 0 ,g1 (u ) 8u12 u 22 2 0,
g 2 (u ) u1 u 2 1 0, U 0 u u1 , u 2
E 2 / u1 0, u 2 0 . (11)
138
We notice that ensemble U 0 E 2 is convex, function J (u ) is
convex on ensemble U 0 , since symmetrical matrix
16 - 12
J ' ' (u ) 0, J ' ' (u ) , 0, E 2 , u U 0 .
- 12 20
It is easy to check that functions g1 (u ), g 2 (u ) are convex on ensemble
U 0 . Finally, problem (10), (11) is a problem of the convex programming.
Entering ensembles
U1 u E 2 / g1 (u ) 0 , U 2 u E 2 / g 2 (u ) 0 ,
problem (10), (11) possible to write as
J (u ) inf, u U U 0 U1 U 2 .
U * u* U J (u* ) min J (u ) .
uU
In fact, ensembles U 0 , U 1 , U 2 are convex and closed, moreover
ensemble U is bounded, consequently, ensemble U is bounded and closed,
i.e. U - a compact ensemble in En. As follows from correlation (10) function
J (u ) is continuous (semicontinuous from below) on compact ensemble U.
Then according to the theorem 1 (lecture 2) ensemble U * . Now
problem (10), (11) possible to write as J (u ) min, u U E .
2
139
L (u , ) 8u12 10u 22 12u1u 2 50u1 80u 2
1 (8u12 u 22 2) 2 (u1 u 2 1), (12)
u (u1 , u 2 ) U 0 , λ ( λ1 ,λ2 ) 0 E 2 / 1 0
has saddle point.
30. Lagrange’s function has the form (12), area of its determination
U 0 0, , moreover 1 0, but 2 can be as positive, as negative.
u , U to Lagrange’s function
40. We define saddle point *
*
0 0
L u u* , * , u u* 0, u U 0, (13)
where
* (1* , *2 )
Conditions *i g i (u* ) 0, i 1,2 (the conditions of the
complementing elasticity) are written:
2 2
1* 8u1* u 2* 2 0, *2 (u1* u 2* 1) 0, 1* 0. (14)
140
2 2
1* (8u1* u 2* 2) 0, u1* u 2* 1 0, *2 0, 1* 0. (15)
We notice that from the second equation (14) and u* U follows that
u u 1 0, 0 . We solve the system of the algebraic equations
*
1
*
2
*
2
(15). It is possible several events:
1) 1* 0, *2 0. In this case the point u* is defined from solution of
2 2
the system 8u1* u 2* 2, u1* u 2* 1 . Thence we have u1* 0,47,
u 2* 0,53. Then value 1* is solution of the equation 126,2 6,461* 0 .
126,2
Thence we have 1 0.
*
6,46
It is impossible, since 1 0.
*
So the case is excluded. Thereby, conditions (13), (14) are not executed in
the internal points of the ensemble U 0 . It remains to consider the border
points of the ensemble U0 .
b) Since the pointu* U , that remains to check the conditions (13),
(14) in the border points (1,0), (0,1) of the ensemble U 0 . We notice that
border point (1,0) U , since in the point restriction 8u1 u 2 2 0.
2 2
is not executed. Then the conditions (13), (14) follows to check in the
singular point u* (0,1) U 0 . In the point g1 (u* ) 1 0,
consequently, value 1* 0. and equalities (14) are executed. Since
derivative L u (u* , ) (38 *2 ,60 *2 ), that inequality (13) is
*
141
*
50. Problem (10), (11) has the solutions u1 0, u 2* 1, J (u* ) 70 .
To solve the following problems:
1. J (u ) u12 u 22 6 sup; u1 u 2 5, u1 0.
2. J (u ) 2u1 2u1 4u 2 3u3 8 sup;
2
142
P.1.4. NONLINEAR PROGRAMMING
J (u ) inf, (1)
u U , (2)
U u E n / u U 0 , g i (u ) 0, i 1, m; g i (u ) 0, i m 1, s ,
where function J (u) C 1 (U 01 ), g i (u) C 1 (U 01 ), i 1, s, U 01 -
open ensemble contained the convex ensembles U 0 from E n ,
in particular, U 01 E n .
Theorem 1 (the necessary optimality condition). If the functions
J (u ) C 1 ( E n ), g i (u ) C 1 ( E n ), i 1, s, int U 0 , U 0 is a
convex ensemble, but ensemble U * , then for each point
u* U necessary exist the Lagrange’s coefficients
*
(*0 , 1* ,..., *s ) 0 E s 1 / 0 0, 1 0,..., m 0 ,
such the following conditions are executed:
*
0, *0 0, 1* 0, ..., *m 0, (3)
*
L u (u* , )u u *
S
*0 J (u* ) *i g i (u* ), u u* 0, u U 0 , (4)
i 1
E S 1 / 0 0, 1 0,..., m 0 .
143
At solution of the problem (1), (2) necessary to consider separately two
events: 1) *0 0 (degenerate problem); 2) *0 0 (nondegenerate
problem). In this case possible to take 1.
*
0
* *
g i (u* ) g (u )
y 0, i I , i * y 0, i m 1, s.
u u
We consider the solutions of the following example on base of the
theorems 1, 2 and solution algorithm of the nonlinear programming (the
lectures 13).
Example 1. To find the problem solution
at conditions
3u1 4u 2 8, u1 4u 22 2, u1 0, u 2 0.
(7)
Solution. Problem (6), (7) is tantamount to the problem
144
J (u ) 2u 22 2u1 u 2 6u1 2u 2 3 inf, (8)
u (u1 , u 2 ) U u E 2 / u U 0 , g1 (u ) 0, g 2 (u ) 0 , (9)
where U 0
u E 2 / u1 0, u2 0 ; g1 (u ) 3u1 4u2 8; g 2 (u)
) u1 4u 22 2 .
Function J (u ) is not a convex function on ensemble U 0 ,
consequently, we have the nonlinear programming problem.
10. We show that ensemble U * . Let U1 u E 2 / g1 (u ) 0,
U 2 u E 2 / g 2 (u ) 0 .
ThenU U 0 U1 U 2 - closed bounded ensemble, consequently,
it is convex. Function J (u ) is continuous on ensemble U . Thence in
effect Weierstrass theorems we have U * .
20. Generalized Lagrange’s function for problem (8), (9) has the form
L (u1 , u 2 , 0 , 12 ) 0 (2u 22 2u1u 2 6u1 2u 2 3)
1 (3u1 4u 2 8) 2 (u1 4u 22 2), u (u1 , u 2 ) U 0 ,
(0 , 1 , 2 ) 0 E / 0 0, 1 0, 2 0.
3
L u u* , * , u u* 0, u U 0 , (11)
where derivative
145
*0 (2u 2* 6) 31* *2
L u (u* , * ) * .
* *
0 (4u 2 2u1 2) 41 82 u 2
* * *
J (u1* , u2* ) 9,07 . However 1* 1,216, *2 0,5 0 , so the given
point u* (1,43; 0,926) can not be solution of the problem (8), (9).
1* 0, *2 0. In this case from (14) we have 2u 2* 6 *2 0,
2)
4u 2* 2u1* 2 8*2 u 2* 0, u1* 4u 2*2 2 0.
146
Thence we have u1 117, u 2 5,454,
* *
*2 4,908 0 . However
the point u* (117; 5,454) U , since the inequality 3u1 4u 2 8.
* *
aren’t executed.
3) 1* 0, *2 0. Equations (14) are written so: 2u 2* 6 31* 0,
4u 2* 2u1* 2 41* 0, 3u1* 4u 2* 8 0.
Thence we have u1* 28 / 9, u 2* 1 / 3 0, 1* 16 / 9 0 . The
point u* (28 / 9, 1 / 3) U .
4) 1 0, 2 0. In this case from equation (14) we have
* *
2u 2* 6 0, 4u 2* 2u1* 2 0.
Thence we have u 2* 3, u1* 5. However the point u* (5; 3) U
, since inequality u1 4u 2 2 0. are not executed. Thereby, the point
* *
u* intU 0 .
b) We expect that point u* ГрU 0 . Here possible the following
events:
1) u (0, u 2 ) ГрU 0 , u 2 0; 2) u u1 ,0 ГрU 0 , u1 0.
For the first type point of the restriction g1 4u 2 8 0,
g 2 4u22 2 0 , consequently, 0 u2 1/ 2 . Then.
u* u 0, u2 1 / 2 , J u* 3,4 . For the second type of the
*
1
147
3 2 1 2
S u1 , u 2 u1 u 2 .
36 16
Now optimization problem can be formulated so: to minimize the
function
3 2 1 2
J u1 , u 2 u1 u 2 inf (15)
36 16
at conditions
u1 u 2 l , u1 0, u 2 0. (16)
3 2 1 2
J u1 , u 2 u1 u 2 inf (17)
36 16
at conditions
u U u E 2 g1 (u ) u1 u 2 l 0, g 2 (u1 , u 2 ) u1 0,
g 3 (u1 , u 2 ) u 2 0 , (18)
148
3 2 1 2
L (u , ) L (u1 , u 2 , 0 , 1 , 2 , 3 ) 0 (u1 u 2 )
36 16
1 (u1 u 2 l ) 2 u1 3 u 2 , u u1 , u 2 E 2 ,
0 , 1 , 2 , 3 0 E 4 / 0 0, 2 0, 3 0.
3 * *
*
u1 0 1* *2
L (u* , ) 18 0, (19)
1 * * *
u 2 0 1 3
*
8
*
0, *0 0, *2 0, *3 0 . (21)
3 * 1
u1 1* *2 0, u 2* 1* *3 0, u1* u 2* l 0,
18 8
149
9l 4l 3
1. *2 0, *3 0 . Then u1* , u2* ,
94 3 94 3
l 3
1* .
18 8 3
2. *2 0, *3 0 . In this case we have u1* 0, u2* l ,
1* l / 8, *2 l / 8 .
3. 2 0, 3 0 . Then from systems of the equations (22) we get
* *
u u 2* l. is not executed.
*
1
3 2 1 2
40. Quadric form y L uu (u * , ) y y1 y 2 .
*
18 8
Hyperplane equations for events 1 - 3 are accordingly written:
g1 u* g u
1) y1 1 * y2 y1 y2 0 . Thence we have y1
u1 u 2
y 2 . Then y L uu (u* , * ) y 4 3 9 y12 0 , i.e. the point
( u 9l (9 4 3 ) , u 4l 3 (9 4 3 ) ) is not the point of the
* *
1 2
local minimum.
g1 (u* ) g (u )
y1 1 * y 2 y1 y 2 0,
u1 u 2
g 2 (u* ) g (u )
y1 2 * y 2 y1 0 .
u1 u 2
Thence we have y1 0, y 2 0; y L uu (u* , * ) y 0, y 0 . The
sufficient conditions of the local minimum are degenerated.
g 1 (u* ) g (u ) g 3 (u* )
3) y1 1 * y 2 y1 y 2 0 , y1
u1 u 2 u1
g (u )
3 * y2 y2 0 .
u 2
150
The data of the equation have the solutions y1 0, y2 0 . Again the
sufficient optimality conditions do not give the onedigit answer. Solution of
the problem is found by comparison of function values J (u1 , u 2 ) in the
last two events. In the second event value J (u1* , u 2* ) l 2 / 16 , but in the
third event J (u1* , u 2* ) 3l 2 / 36. Then the problem solution is the
point ( u1 0, u 2 l ), i.e. from the whole wire is made only square.
* *
151
b) J (u ) u1u 2 u 3 inf,
u1 u 2 u3 6, u1u 2 u1u3 u 2u3 12;
1 1 1 1
c) J (u ) inf, 2
2 1;
u1 u 2 u1 u 2
d) J u 2u1 3u 22 u 32 inf,
u1 u 2 u 3 8, u i 0, i 1,3;
e) J u u12 u 22 u 32 inf,
u1 u 2 u 3 12, u i 0, i 1,3 ;
f) J (u ) u1u 2 u1u 3 u 2 u 3 inf, u1 u 2 u 3 4;
g) J (u ) u12 u 2 u 22 u1 u1u 2 u 3 inf,
u1 u2 u3 15, ui 0, i 1,3;
l) J (u ) u12 2u1u 2 u 32 inf,
u1 2u 2 u 3 1, 2u1 u 2 u 3 5, u i 0, i 1,3.
10. To find the conditional extreme to function J (u ) (u1 2) 2
(u2 3) 2 at condition u12 u22 52 .
u U u E n / u 0, Au b , (2)
where
A a1 , a 2 ,..., a m , a m1 ,..., a n - matrix of the order m n;
a i E m , i 1, n ; vectors are identified by the condition vectors, but
b E m - by the restriction vector.
152
Simplex-method is a general method of the problem solution of the
linear programming in canonical form (1), (2). Since the general and the
basic problems of the linear programming are reduced to type (1), (2), it is
possible to consider that simplex-method is the general method of the
problem solution of the linear programming.
The base of the linear programming theory is stated in the lectures 15 -
17. We remind briefly a rule of the nondegenerate problem solution of the
linear programming in canonical form.
10. To build the initial extreme point u 0 U ensemble U . We
notice, if U * , that lower border to linear function (1) on U is
reached in the extreme point ensemble U , moreover in the nondegenerate
problem the extreme point has exactly m positive coordinates i.e.
u 0 u10 0,..., u m0 0,0,...,0 , and vectors a1 ,..., a m , corresponding
to the positive coordinates of the extreme point are linear independent. The
extreme point u U is defined in the general event by М-method
0
m
i.e. a
j
a u
i 1
i
ij AБ u j , where AБ a1 ,..., a m - nonsingular
153
30. The extreme point u1 U and corresponding to it simplex-table is
built on base of the simplex-table of the point u U . The index j0 is
0
ui0
defined from condition min i0 amongst uij0 0 . The extreme point
uij
u 1 (u10 0 u1 j0 ,..., u i00 1 0 u i0 1 j0 , 0, ui00 1 ui0 1 j0 ,..., um0 0umj0 ,
0,...,0, 0 ,0,...,0), 0 i0 .
The base consists of the vectors a1 ,..., a i0 1 , a j0 , a i0 1 ,..., a m .
Decomposition coefficients of the vectors a j , j 1, n on the base are
defined by the formula
u ij0 u i0 j
(u ij ) нов u ij , i i0 , j j 0 ;
u i0 j0
(3)
u i0 j
(u i0 j ) нов , j 1, n.
u i0 j0
solution of the problem (1), (2), but otherwise transition to the new extreme
point u U realized and etc.
2
154
u1 2u 2 u 4 3u 5 17, 4u 2 u 3 u 5 12,
u 2 8u 4 u 5 u 6 6,
u1 0, u 2 0, u 3 0, u 4 0, u 5 0, u 6 0 .
Matrix A , vectors b and c are equal to
1 2 0 1 3 0
A 0 4 1 0 1 0 a1 , a 2 , a 3 , a 4 , a 5 , a 6 ,
0 1 0 8 1 1
1 2 0 1 3 0
2 3 4 5 6
a 0 , a 4 , a 1 , a 0 , a 1 , a 0 ,
1
Base
с b
6 -3 5 -2 -4 2
0
a1 a2 a3 a4 a5 a6
I u0=(17, 0, 12, 0, 0, 6)
17
A1 6 17 1 2 0 1 3 0
2
12
A3 5 12 0 (4) 1 0 1 0 i0 3
4
6
A6 2 6 0 1 0 8 -1 1
1
155
0
zj cj 0 37 0 24 25 0 J (u ) 174
z j0 c j0 37,
j0 2
II u1=(11, 3, 0, 0, 0, 3)
a1 6 11 1 0 -1/2 1 5/2 0 11
a2 -3 3 0 1 1/4 0 1/4 0 -
a6 2 3 0 0 -1/4 (8) -5/4 1 3/8 i0 6
zj cj 0 0 - 37
4
24
63
4
0 1
J (u ) 63
z j0 c j0 24, j0 4
85
a1 6 8 1 0 - 15 0 ( 85 ) -1 4
32 32 8
1 1
a2 -3 3 0 1 4 0 4 0 -
a4 -2
3
0 0 - 1 1 - 5 1
- i0 1
8 32 32 8
2
J (u ) 54
zj cj 0 0 - 17
2 0
39
2 -3
39
z j0 c j0 , j0 5
2
IV u3=(0, 2, 0, 1, 4, 0)
zj cj 0
32
a5 -4 4 0 - 15 0 1 - 4
85 85 85 j 1,6
a2 -3 2 - 8 1
25
0 0
1
85 85 85
156
10
a4 -2 1
1
0 - 5 1 0
17 85 85
3
zj cj - 624
85
0 - 430
85
0 0 - 177
85
J ( u ) 24
Example 2.
J (u ) 2u1 3u2 5u3 6u4 4u5 inf;
2u1 u 2 u3 u 4 5; u1 3u 2 u3 u 4 2u5 8;
u1 4u 2 u 4 1, u j 0, j 1,5.
Corresponding М-problem has the form
J (u ) 2u1 3u 2 5u 3 6u 4 4u 5 Mu6 Mu7 inf;
2u1 u 2 u3 u 4 u6 5; u1 3u 2 u3 u 4 2u5 8;
u1 4u 2 u 4 u 7 1, u j 0, j 1,7.
2 1 1 1 0 1 0
,
A 1 3 1 1 2 0 0 (a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 )
1 4 0 1 0 0 1
5
b 8 , c (2,3,5,6,4, M , M ) .
1
Solution. The initial extreme point u0=(0, 0, 0, 0, 4, 5, 1). We present
value z j c j in the manner of z j c j j M j , j 1, n , i.e.
instead of one line for z j c j are entered two: in the first values j are
written, in the second - j . Since M – sufficiently great positive number,
then z j c j z k ck , if j k . But if j k , that j k .
-2 -3 5 -6 4 М М
Base с B
a1 a2 a3 a4 a5 а6 а7
157
I J (u ) 16 6M u 0 (0,0,0,0,4,5,1)
A6 M 5 2 1 -1 1 0 1 0
A5 4 4 1/2 3/2 1/2 -1/2 1 0 0
A7 M 1 -1 (4) 0 1 0 0 1
j 4 9 -3 4 0 0 0
zj cj
j 1 5 -1 2 0 0 0
j 9/4 0 -1 3/4 0 0
158
u1 u 2 u3 5; u j 0, j 1,4 .
2. J (u ) u1 u 2 u3 3u 4 u5 u6 3u7 inf;
3u3 u5 u6 6;
u 2 2u3 u 4 10;
u1 u6 0;
u3 u6 u7 6; u j 0, j 1,7 .
3. J (u ) u1 2u 2 u3 2u4 u5 inf;
u1 2u2 u3 2;
2u1 u2 u4 0;
u1 3u2 u5 6;
u1 0; u 2 0, u3 0, u 4 0.
4. J (u ) 2u1 3u 2 2u3 u 4 sup;
2u1 2u 2 3u3 u 4 6;
u 2 u3 u4 2;
u1 u2 2u3 5;
u1 0; u 2 0, u3 0.
5. J (u ) u1 2u 2 u3 sup;
3u1 u2 u3 4;
2u1 3u2 6;
u1 2 0 u 2 1, u j 0, j 1,3.
6. J (u ) u1 2u2 u3 inf;
u1 u2 1;
u1 2u2 u3 8;
u1 3u 2 3, u j 0, j 1,3.
7. J (u ) 2u1 u 2 u3 sup;
159
u1 2u2 u3 4;
u1 u2 2;
u1 u3 1, u j 0, j 1,3.
8. The steel twig by length 111 cm have gone into blanking shop. It is
necessary to cut them on stocking up on 19, 23 and 30 cm in amount 311,
215 and 190 pieces accordingly. To build the model on base of which it is
possible to formulate the extreme problem of the variant choice performing
the work under which the number cut twig are minimum.
9. There are two products which must pass processing on four
machines (I, II, III, IV) in process production. The time of the processing of
each product on each of these machines is specified in the following form:
Machine А В
I 2 ј
II 4 2
III 3 1
IV 1 4
The machines I, II, III and IV can be used accordingly during 45, 100,
300 and 50 hours. Price of the product А - 6 tenge per unit, product В - 4
tenge. In what correlation do follows to produce the product А and В to get
maximum profit? To solve the problem in expectation that product А is
required in amount not less 22 pieces.
10. Refinery disposes by two oil grade - А and В, moreover petrol and
fuel oil are got under processing. Three possible production processes are
characterized by the following scheme:
a) 1 unit of the sort A + 2 unit of the sort B 2 unit of the fuel oil + 3
unit of petrol;
b) 2 unit of the sort A + 1 unit of the sort B 5 unit of the fuel oil + 1
unit of petrol;
c) 2 unit of the sort A + 2 unit of the sort B 2 unit of the fuel oil + 1
unit of petrol;
We assume the price of the fuel oil 1 tenge per unit, but the price of the
petrol 10 tenge per unit. To find the most profitable production plan if there
are 10 units to oil of the sort А and 15 units oil of the sort B.
160
Appendix 2
Three term tasks for students of the 2-nd, 3-rd courses (5,6 terms
accordingly) educating on the specialties “Applied mathematics”,
“Mathematics”, “Mechanics”, “Informatics”, “Economical cybernetics” are
worked out by the following parts:
convex ensembles and convex function (1-th task),
convex and nonlinear functions (2-nd task),
linear programming (3-rd task),
Further we prefer variants of the tasks on course “Mathematical
programming”.
Task 1
2
5. J u u1 u 2 2u3 u1u 2 u1u3 u 2u3
2 2 2
5u2 25; U E 3 .
6. J u u2
1 2
2
5
u3 7u1 u3 6; U u E 3 : u 0 .
7. J u 3u1 u2 2u3 u1u2 3u1u3 u2u3 3
2 2 2
3u2 6; U E 3.
161
1 2
8. J u 5u1 u2 4u32 u1u2 2u1u3 2u2u3
2
2
u3 1; U E 3.
1 2 1
9. J u 2u1 u 2 5u 3 u1u 2 2u1u 3 u 2 u 3
2 2
2 2
3u1 2u 2 6; U E 3 .
10. J u u13 2u32 10u1 u 2 5u3 6;
U {u E 3 : u 0}.
11. J u 5u14 u26 u32 13u1 7u3 8; U E 3.
12. J u 3u12 2u22 u32 3u1u2 u1u3
2u2u3 17; U E 3.
1 4
13. J u 4u13 u24 u3 3u1 8u2 11;
2
U {u E 3 : u 0}.
14. J u 8u13 12u32 3u1u3 6u2 17;
U {u E 3 : u 0}.
15. J u 2u12 2u 22 4u32 2u1u 2 2u1u3
2u 2u3 16; U E 3 .
1 2
16. J u 2u1 u2 u3 2u1u2 8u3 12; U E 3 .
2 2
2
1 7 1 4
17. J u u2 u3 2u2u3 11u1 6;
2 2
U {u E 3 : u 0}.
5 3 1
18. J u u12 u 22 4u 32 u1u 2 2u1u 3 u 2 u 3
2 2 2
+ 8u3 13; U E 3 .
162
1 3
19. J u 3u1 u2 2u1u2 5u1u3 7u1 16;
2
2
U {u E 3 : u 0}.
3 2 1
20. J u 2u12 u22 u3 u1u2 u1u3
2 2
2u2u3 10; U E 3 .
3 5
21. J u 2u12 u32 u1u3 12u2 18; U E 3 .
2 2
22. J u 6u12 u23 6u32 12u1 8u2 7;
U {u E 3 : u 0}.
3 2
23. J u u1 u2 2u3 u1u2 u1u3 2
2 2
2
2u2u3 8u2 ; U E 3.
5 2
24. J u 4u1 u2 u3 4u1u2 11u3 14; U E .
2 2 3
2
7 2 4 3 1 3
25. J u u1 u 2 u 3 13u1 7u 3 9;
2 2 2
U {u E : u 0}.
3
5 3 1 5 3 3
26. J u u1 u2 u3 22u2 17;
6 4 2
U {u E : u 0}.
3
5 3 3 2
27. J u u1 2u 2 u 3 2u1u 2 3u1u 3 u 2 u 3 ;
2
6 2
U {u E : u 0}.
3
3 2 5 2 9 2
29. J u u1 u2 u3 3u1u3 7u2u3 ; U E 3 .
2 2 2
163
7 3 5 2 5 4 1 3
30. J u u1 u 2 u 3 u 4 3u 2 ;
6 2 12 2
U {u E : u 0}.
3
1 3
31. J u 3u1 u 2 2u1u 2 5u1u 3 7u1 16;
2
2
U {u E 3 : u 0}.
5 2 3
32. J u u1 u 22 4u 32 u1u 2 2u1u 3
2 2
1
u 2 u 3 8u 3 13; U E 3.
2
1 7 1 4
33. J u u2 u3 2u2u3 11u1 6;
2 2
U {u E 3 : u 0}.
1 2
34. J u 2u1 u2 u3 2u1u2 8u3 12; U E 3 .
2 2
2
35. J u 2u1 2u2 4u3 2u1u2 2u1u3
2 2 2
2u2u3 16; U E 3.
7 3 5 2 5 4 1 3
36. J u u1 u 2 u 3 u 4 3u 2 ;
6 2 12 2
U {u E : u 0}.
4
3 2 5 2 9 2
37. J u u1 u2 u3 3u1u3 7u2u3 ; U E .
3
2 2 2
38. J u 2u1 2u2 u3 u1u2 9u1u3 u2u3 9; U E .
2 2 2 3
5 3 3 2
39. J u u1 2u 2 u 3 2u1u 2 3u1u 3 u 2 u 3 ;
2
6 2
U {u E : u 0}.
3
164
5 3 1 5 3 3
40. J u u1 u2 u3 22u2 17;
6 4 2
U {u E : u 0}.
3
7 2 4 3 1 3
41. J u u1 u 2 u 3 13u1 7u 3 9;
2 3 2
U {u E : u 0}.
3
5 2
42. J u 4u1 u2 u3 4u1u2 11u3 14; U E .
2 2 3
2
3
43. J u u1 u2 2u3 u1u2 u1u3
2 2 2
2
2u2u3 8u2 ; U E 3.
44. J u 6u12 u32 6u32 12u1 8u2 7;
U {u E 3 : u 0}.
3 2 5
45. J u 2u1 u3 u1u3 12u2 18; U E .
2 3
2 2
3 1
46. J u 2u1 u2 u3 u1u2 u1u3
2 2 2
2 2
2u2u3 10; U E 3.
47. J u u16 u22 u32 u42 10u1 3u4 20; U E 4 .
48. J u e
2 u1 u 2
; U E 2.
49. J u u13 u23 u33 10u1 u2 15u3 10; U E 3.
1 2
50. J u u1 u2 u3 u1u2 u3 10; U E .
2 2 3
2
51. J u u1 u2 2u3 u1u2 u1u3 u2u3
2 2 2
5u2 25; U E 3.
1 2
52. J u u2 u3 7u1 u3 6; U {u E 3 : u 0}.
5
165
53. J u 3u12 u22 2u32 u1u2 3u1u3
u2u3 3u2 6; U E 3.
1 2
54. J u 5u1 u 2 4u 32 u1u 2 2u1u 3 2u 2 u 3 u 3 1;
2
2
U E 3.
55. J u u1 2u3 10u1 u2 5u3 6;
3 3
U {u E 3 : u 0}.
56. J u 5u14 u26 u32 13u1 7u3 8; U E 3.
57. J u 3u12 2u22 u32 3u1u2 u1u3
2u3u2 17; U E 3.
1 4
58. J u 4u1 u2 u3 3u1 8u2 11;
3 4
2
U {u E 3 : u 0}.
59. J u 3u13 12u32 3u1u3 6u2 17;
U {u E 3 : u 0}.
60. J u 2u12 2u22 4u32 2u1u2 2u1u3
2u2u3 16; U E 3.
61. J u u16 u22 u32 u42 10u1 3u4 20; U E 4 .
62. J u e
2 u1 u 2
; U E 2.
63. J u u13 u23 u33 10u1 u2 15u3 10; U E 3.
1 2
64. J u u1 u2 u3 u1u2 u3 10; U E .
2 2 3
2
65. J u u1 u2 2u3 u1u2 u1u3 u2u3
2 2 2
5u2 25; U E 3.
1 2
66. J u u2 u3 7u1 u3 6; U {u E 3 : u 0}.
5
166
67. J u 3u12 u22 2u32 u1u2 3u1u3 u2u3
3u2 6; U E 3.
1 2
68. J u 5u1 u2 u32 u1u2 2u1u3 2u2u3
2
2
u3 1; U E 3.
1 2 1
69. J u 2u1 u 2 5u 3 u1u 2 2u1u 3 u 2 u 3
2 2
2 2
3u1 2u 2 6; U E 3 .
70. J u u13 u33 10u1 u2 5u3 6;
U {u E 3 : u 0}.
71. J u 5u14 u26 u32 13u1 7u3 8; U E 3.
72. J u 3u12 2u22 2u32 3u1u2 u1u3
2u2u3 17; U E 3.
1 4
73. J u 4u3 u24 u3 3u1 8u2 11;
2
U {u E 3 : u 0}.
74. J u 8u13 12u32 3u1u3 6u2 17;
U {u E 3 : u 0}.
75. J u 2u12 2u22 4u32 u1u2 2u2u3 16; U E 3.
1 2
76. J u 2u1 u2 u3 2u1u2 8u3 12; U E .
2 2 3
2
1 7 1 4
77. J u u2 u3 2u2u3 11u1 6;
2 2
U {u E : u 0}.
3
167
5 2 3
78. J u u1 u 22 4u 32 u1u 2 2u1u 3
2 2
1
u 2 u 3 8u 3 13; U E 3.
2
1 3
79. J u u1 u 2 2u1u 2 5u1u 3 7u1 16;
2
2
U {u E 3 : u 0}.
3 2 1
80. J u 2u1 u2 u3 u1u2 u1u3
2 2
2 2
2u2u3 10; U E 3.
3 2 5
81. J u u1 u3 u1u3 12u2 18; U E .
2 3
2 2
82. J u 6u1 u2 6u3 12u1 8u2 17;
2 3 2
; U {u E 3 : u 0}.
3 2
83. J u u1 u2 2u2 u1u2 u1u3 2u2u3 8u2 ; U E .
2 2 3
2
7 2 4 3 1 3
84. J u u1 u 2 u 3 13u1 7u 3 9;
2 3 2
U {u E : u 0}.
3
5 3 1 5 3 3
85. J u u1 u2 u3 22u2 10;
6 4 2
U {u E : u 0}.
3
5 3 3 2
86. J u u1 2u 2 u 3 2u1u 2 3u1u 3 u 2 u 3 ;
2
6 2
U {u E : u 0}.
3
3 2 5 2 9 2
88. J u u1 u2 u3 3u1u3 7u2u3 ; U E .
3
2 2 2
168
7 3 5 2 5 4 1 3
89. J u u1 u2 u3 u4 3u2 ;
6 2 12 2
U {u E : u 0}.
3
Task 2
2
u1 3, u2 6, u1 0, u2 0.
3 2
3. J u 4u1 8u 2 u1 u 2 2u1u 2 max,
2
2
u1 u2 3, u1 0, u1 u2 1, u2 0.
3 2
4. J u 4u1 8u 2 u1 u 2 2u1u 2 max,
2
2
u1 u2 1, u1 0, u1 4, u2 0.
3 2
5. J u 4u1 8u 2 u1 u 2 2u1u 2 max,
2
2
3u1 5u2 15, u1 u2 1, u1 0, u2 0.
1 2
6. J u 3u1 2u 2 u1 u 2 u1u 2 max,
2
2
u1 2u2 2, u1 0, 2u1 u2 2, u2 0.
7. J u u1 6u 2 u1 3u 2 3u1u 2 max,
2 2
169
8. J u u1 6u 2 u1 3u 2 3u1u 2 max,
2 2
u1 u2 3, u1 0, 2u1 u2 2, u2 0.
9. J u u1 6u 2 u1 3u 2 3u1u 2 max,
2 2
u1 u2 0, u1 0, u2 5, u2 0.
3 2
10. J u 6u 2 u1 u 2 2u1u 2 max,
2
2
3u1 4u2 12, u1 0, u1 u2 2, u2 0.
3 2
11. J u 6u 2 u1 u 2 2u1u 2 max,
2
2
u1 2u2 2, u1 0, u1 2, u2 0.
3 2
12. J u 6u 2 u1 u 2 2u1u 2 max,
2
2
3u1 4u2 12, u1 0, u1 2u2 2, u2 0.
3 2
13. J u 8u1 12u 2 u1 u 2 max,
2
2
2u1 u2 4, u1 0, 2u1 5u2 10, u2 0.
3 2
14. J u 8u1 12u 2 u1 u 2 max,
2
2
u1 2u2 2, u1 0, u1 6, u2 0.
3 2
15. J u 8u1 12u 2 u1 u 2 max,
2
2
3u1 2u2 0, u1 0, 4u1 3u2 12, u2 0.
1 2
16. J u 3u1 2u 2 u1 u 2 u1u 2 max,
2
2
2u1 u2 2, u1 0, 2u1 3u2 6, u2 0.
1 2
17. J u 6u1 4u 2 u1 u 2 u1u 2 max,
2
2
u1 2u2 2, u1 0, 2u1 u2 0, u2 0.
170
1 2
18. J u 6u1 4u 2 u1 u 2 u1u 2 max,
2
2
2u1 u2 2, u1 0, u2 1, u2 0.
1 2
19. J u 6u1 4u 2 u1 u 2 u1u 2 max,
2
2
3u1 2u2 6, u1 0, 3u1 u2 3, u2 0.
20. J u 8u1 6u2 2u1 u2 max,
2 2
u1 u2 1, u1 0, 3u1 2u2 6, u2 0.
21. J u 8u1 6u2 2u1 u2 max,
2 2
u1 u2 1, u1 0, u1 3, u2 0.
22. J u 8u1 6u2 2u1 u2 max,
2 2
2u1 u2 4, u1 0, u1 u2 2, u2 0.
25. J u 2u1 2u 2 u1 2u 2 2u1u 2 max,
2 2
2u1 u2 2, u1 0, u 2 4, u 2 0.
26. J u 4u1 4u 2 3u1 u 2 2u1u 2 max,
2 2
1 1
u1 u2 6, u1 0, u1 u2 1, u2 0.
2 2
171
11 1 2 1
30. J u u1 u 2 u12 u 22 u1u 2 max,
2 6 3 2
2u1 u2 2, u1 0, u1 2u2 2, u2 0.
31. J u 18u1 12u 2 2u1 u 2 2u1u 2 max,
2 2
1
u1 u2 4, u1 0, u1 u2 1, u2 0.
2
1 2 5 2
32. J u 6u1 16u 2 u1 u 2 2u1u 2 max,
2 2
5u1 2u2 10, u1 0, 3u1 2u 2 6, u 2 0.
33. J u 11u1 8u 2 2u1 u 2 u1u 2 max,
2 2
u1 u2 2, u1 0, u1 5, u2 0.
13 1
39. J u u1 5u 2 2u12 u 22 u1u 2 max,
2 2
u1 u2 3, u1 0, u1 u2 1, u2 0.
9 27 1 1
40. J u u1 u 2 u12 2u 22 u1u 2 max,
2 2 2 2
u1 u2 2, u1 0, u1 u2 4, u2 0.
172
41. J u 2u1 8u 2 u1 5u 2 4u1u 2 max,
2 2
u1 u2 3, u1 0, 2u1 3u2 6, u2 0.
42. J u 8u1 18u 2 u1 2u 22 2u1u 2 max,
2
u1 u2 3, u1 0, u1 2, u2 0.
23 3
43. J u u1 u12 2u 22 u1u 2 max,
2 2
5u1 4u2 20, u1 0, u1 u2 2, u2 0.
44. J u 48u1 28u 2 4u1 2u 22 4u1u 2 max,
2
2u1 u2 6, u1 0, 2u1 u2 4, u2 0.
45. J u u1 10u 2 u1 2u 22 u1u 2 max,
2
3u1 2u2 6, u1 0, u1 u2 4, u2 0.
49. J u u1 u2 6 max, u1 u2 5, u2 0.
2 2
173
54. J u 10u1 5u 2 u1 2u 2 10 min,
2 2
2u12 u 2 4, u 2 0, u1 u2 8.
55. J u u3 u1u 2 6 min,
2
2u2 u3 3, u j 0, u1 u2 u3 2, j 1,2,3 .
56. J u 2u1 3u2 u1 6 max,
2 2
u12 u2 3, u2 0, 2u1 u2 5.
57. J u u1 3u1 5 min,
2
2
u12 4u1 u2 5, u12 6u1 u2 7.
2
2u1 9u 2 18, u 2 0, u1 u2 1.
2 2
3u12 6 u2 0, u12 u 12 9, u1 0, u 2 0.
5 2
64. J u 3u1 u 2 3u1u 2 7 min,
2
2
3u1 u 2 1, u22 2.
174
65. J u u1 2u1u2 u1 6 max,
6
u1 u2 1, u1 0 , 2u1 u2 5, u2 0.
70. J u 2u1 3u 2 u1u 2 6 max,
2 2
u1 u2 3. u12 u2 5, u1 0
71. J u u1 u 2 u1 5u 2 5 max, , u1 u2 5.
2 2
u12 u 2 u3 5, u1 5u 2 8, u1 0, u 2 0.
74. J u u12 u22 u32 min,
u1 u2 u3 3, 2u1 u2 u3 5.
75. J u 3u12 u22 u32 u1u2 6 min, 2
2u1 u2 u3 5, u2 3u3 8, u2 0.
76. J u u1u 2 u1u 3 2u 2 u 3 u1 5 max,
u1 u2 u3 3, 2u1 u 2 5, u1 0.
1 2 3 2
77. J u u1 u 2 u 3 12u1 13u 2 5u 3 min,
2
2 2
u1 5u2 4u3 16 , 2u1 7u 2 3u 3 2, u1 0.
175
78. J u u1 2u 2 30u1 16u 3 10 min,
2 2
2
u1 u2 2u3 3, 2u1 u 2 3u 3 11, u j 0, j 1,2,3.
176
Task 3
J u cu min; U u E n / Au b, u 0 .
3 1 1 1 1
1. c (5,1,1,0,0), b 5,4,11, A 2 1 3 0 0 .
0 5 6 1 0
1 2 1 6 1
2. c (6,1,1,2,0), b 4,1,9, A 3 1 1 1 0 .
1 3 5 0 0
3 1 1 6 1
3. c (0,6,1,1,0), b 6,6,6 , A 1 0 5 1 7 .
1 2 3 1 1
5 1 1 3 1
4. c (7,1,1,1,0), b 5,3,2 , A 0 2 4 1 1 .
1 3 5 0 0
1 1 1 2 1
5. c (8,1,3,0,0), b 4,3,6 , A 2 0 1 3 5 .
3 0 1 6 1
2 1 2 0 0
6. c (0,1,3,1,1), b 2,8,5, A 1 1 4 1 3 .
3 1 1 0 6
2 0 1 1 1
7. c (1,2,1,1,0), b 2,7,2 , A 4 1 3 1 2 .
1 0 1 2 1
177
6 1 1 2 1
8. c (0,1,6,1,3), b 9,14,3, A 1 0 1 7 8 .
1 0 2 1 1
2 0 3 1 1
9. c (8,1,1,1,0), b 5,9,3, A 3 1 1 6 2 .
1 0 2 1 2
2 0 3 1 0
10. c (1,3,1,1,0), b 4,4,15, A 1 0 1 2 3 .
3 3 6 3 6
4 1 1 0 1
11. c (0,2,0,1,3), b 6,1,24 , A 1 3 1 0 3 .
8 4 12 4 12
8 16 8 8 24
12. c (10,5,25,5,0), b 32,1,15, A 0 2 1 1 1 .
0 3 2 1 1
4 1 1 2 1
13. c (6,0,1,1,2), b 8,2,2 , A 2 1 0 1 0 .
1 1 0 0 1
1 2 3 4 1
14. c (5,1,3,1,0), b 7,7,12 , A 0 3 1 4 0 .
0 4 0 8 1
3 4 1 0 0
15. c (5,3,2,1,1), b 12,16,3, A 3 2 1 1 1 .
1 3 0 0 1
178
1 1 1 0 0
16. c (7,0,1,1,1), b 1,12,4, A 2 2 1 1 2 .
2 1 0 0 1
1 1 1 0 0
17. c (6,1,2,1,1), b 2,11,6 , A 5 2 1 1 1 .
3 2 0 0 1
2 1 1 1 3
18. c (0,0,3,2,1), b 5,7,2 , A 3 0 2 1 6 .
1 0 1 2 1
6 3 1 1 1
19. c (1,7,2,1,1), b 20,12,6 , A 4 3 0 1 0 .
3 2 0 0 1
1 2 1 0 0
20. c (2,0,1,1,1), b 2,14,1, A 3 5 1 1 2 .
1 1 0 0 1
1 2 1 0 0
21. c (6,1,0,1,2), b 2,18,2, A 2 6 2 1 1 .
1 2 0 0 1
1 2 1 0 0
22. c (0,3,1,1,1), b 2,2,6 , A 1 1 0 1 0 .
2 1 1 1 2
2 2 1 1 1
23. c (3,0,1,2,1), b 6,2,2 , A 2 1 0 1 0 .
1 1 0 0 1
179
1 1 1 0 0
24. c (0,5,1,1,1), b 2,2,10 , A 1 2 0 1 0 .
2 1 1 1 2
3 4 1 0 0
25. c (1,5,2,1,1), b 12,1,3, A 1 1 0 1 0 .
3 2 1 1 1
1 1 1 0 0
26. c (5,0,1,1,1), b 1,3,12 , A 3 1 0 1 0 .
2 2 1 1 2
1 1 1 0 0
27. c (7,0,2,1,1), b 2,3,11, A 3 1 0 1 0 .
5 2 1 1 1
5 5 1 2 1
28. c (1,4,1,1,1), b 28,2,12 , A 1 2 0 1 0 .
3 4 0 0 1
1 2 1 0 0
29. c (0,8,2,1,1), b 2,20,6 , A 6 3 1 1 1 .
3 2 0 0 1
3 5 1 1 2
30. c (0,2,1,1,1), b 14,10,1, A 2 5 0 1 0 .
1 1 0 0 1
1 2 1 0 0
31. c (7,2,0,1,2), b 2,12,18, A 3 4 0 1 0 .
2 6 2 1 1
180
1 2 1 0 0
32. c (1,3,1,1,1), b 2,6,1, A 2 1 1 1 2 .
1 1 0 0 1
1 2 1 0 0
33. c (5,1,1,1,2), b 2,8,2 , A 4 1 1 2 1 .
1 1 0 0 1
1 1 2 2 1
34. c (1,2,1,1,1), b 11,2,3, A 1 2 0 1 0 .
1 1 0 0 1
2 3 1 2 1
35. c (10,5,2,1,1), b 17,1,3, A 1 1 0 1 0 .
1 3 0 0 1
1 0 2 1 3
36. c (2,1,3,1,1), b 6,16,7 , A 2 2 4 8 4 .
1 0 1 7 1
1 1 1 0 0
37. c (4,1,1,2,1), b 2,13,16, A 4 3 2 1 0 .
3 2 0 0 1
4 3 1 0 0
38. c ( 2,2,1,2,1), b 12,2,26 , A 1 2 0 1 0 .
6 3 1 1 1
9 1 1 1 2
39. c (5,2,1,1,1), b 26,12,6 , A 4 3 0 1 0 .
3 2 0 0 1
181
2 6 1 1 1
40. c (1,11,1,2,1), b 13,10,1, A 2 5 0 1 0 .
1 1 0 0 1
3 1 3 1 0
41. c (5,1,1,2,0), b 1,6,2 , A 2 3 1 2 1 .
3 1 2 1 0
2 1 1 1 2
42. c (0,3,1,1,1), b 6,2,1, A 1 1 0 1 0 .
1 1 0 0 1
1 1 1 2 1
43. c (8,1,3,0,0), b 4,3,6, A 2 0 1 3 5 .
3 0 1 6 1
1 1 1 0 0
44. c (2,1,1,1,1), b 2,11,3, A 1 1 2 2 1 .
1 1 0 0 1
4 3 2 1 1
45. c (5,0,1,2,1), b 13,3,6 , A 3 1 0 1 0 .
3 2 0 0 1
1 1 1 0 0
46. c (1,3,1,1,1), b 1,17,4 , A 5 2 2 1 3 .
2 1 0 0 1
3 4 1 0 0
47. c (9,5,2,1,1), b 12,17,3, A 2 3 1 2 1 .
1 3 0 0 1
182
4 3 1 0 0
48. c (1,1,1,2,1), b 12,26,12 , A 6 3 1 1 1 .
3 4 0 0 1
1 2 1 0 0
49. c (0,7,1,1,1), b 2,26,6, A 9 1 1 1 2 .
3 2 0 0 1
1 2 1 0 0
50. c (4,8,1,2,1), b 2,13,1, A 2 6 1 1 1 .
1 1 0 0 1
1 1 0 1 2
51. c (3,1,1,1,0), b 3,6,5, A 2 1 1 2 3 .
3 2 0 3 8
1 2 1 0 0
52. c (1,3,1,2,1), b 2,2,5, A 1 1 0 1 0 .
1 2 1 1 1
1 2 1 0 0
53. c (0,1,1,2,1), b 2,2,6 , A 2 1 0 1 0 .
2 2 1 1 1
1 1 1 0 0
54. c (0,5,1,1,1), b 2,2,11, A 1 2 0 1 0 .
1 1 2 2 1
3 4 1 0 0
55. c (9,2,1,0,1), b 12,1,17 , A 1 1 0 1 0 .
2 3 1 2 1
183
1 1 1 0 0
56. c (1,0,1,1,1), b 1,3,17 , A 3 1 0 1 0 .
5 2 2 1 3
1 1 1 0 0
57. c (3,2,1,2,1), b 2,3,13, A 3 1 0 1 0 .
4 3 2 1 1
1 2 1 0 0
58. c (9,0,1,1,1), b 2,12,26 , A 4 3 0 1 0 .
9 1 1 1 2
6 3 1 1 1
59. c (5,5,1,2,1), b 26,2,12, A 1 2 0 1 0 .
3 4 0 0 1
1 2 1 0 0
60. c (0,10,1,2,1), b 2,10,1, A 2 5 0 1 0 .
2 6 1 1 1
3 1 3 1 2
61. c (3,2,1,1,0), b 5,5,5, A 3 2 1 1 1 .
7 2 2 0 1
5 10 5 15 10
62. c (1,2,1,1,0), b 25,3,5, A 0 1 1 6 2 .
0 6 1 1 1
2 1 0 1 1
63. c (1,1,2,1,0), b 4,7,9 , A 3 2 0 1 1 .
1 1 1 2 6
184
2 1 1 1 1
64. c (1,3,1,0,0), b 4,3,6 , A 1 0 2 1 3 .
3 0 3 1 2
3 1 1 2 3
65. c (2,1,1,5,0), b 7,1,9 , A 2 0 3 2 1.
3 0 1 1 6
4 1 1 2 1
66. c (6,0,1,1,2), b 8,2,2 , A 2 1 0 1 0 .
1 1 0 0 1
1 2 3 4 1
67. c (5,1,3,1,0), b 7,7,12, A 0 3 1 4 0 .
0 4 0 8 1
3 4 1 0 0
68. c (5,3,2,1,1), b 12,16,3, A 3 2 1 1 1 .
1 3 0 0 1
1 1 1 0 0
69. c (7,0,1,1,1), b 1,12,4 , A 2 2 1 1 2 .
2 1 0 0 1
1 1 1 0 0
70. c (6,1,2,1,1), b 2,11,6 , A 5 2 1 1 1 .
3 2 0 0 1
2 1 1 1 3
71. c (0,0,3,2,1), b 5,7,2 , A 3 0 2 1 6 .
1 0 1 2 1
185
6 3 1 1 1
72. c (1,7,2,1,1), b 20,12,6, A 4 3 0 1 0 .
3 2 0 0 1
1 2 1 0 0
73. c (2,0,1,1,1), b 2,14,1, A 3 5 1 1 2 .
1 1 0 0 1
1 2 1 0 0
74. c (6,1,0,1,2), b 2,18,2, A 2 6 2 1 1 .
1 2 0 0 1
1 2 1 0 0
75. c (0,3,1,1,1), b 2,2,6, A 1 1 0 1 0 .
2 1 1 1 2
2 2 1 1 1
76. c (3,0,1,2,1), b 6,2,2, A 2 1 0 1 0 .
1 1 0 0 1
1 1 1 0 0
77. c (0,5,1,1,1), b 2,2,10, A 1 2 0 1 0 .
2 2 1 1 2
3 4 1 0 0
78. c (1,5,2,1,1), b 12,1,3, A 1 1 0 1 0 .
3 2 1 1 1
1 1 1 0 0
79. c (5,0,1,1,1), b 1,3,12, A 3 1 0 1 0 .
2 2 1 1 2
186
1 1 1 0 0
80. c (7,0,2,1,1), b 2,3,11, A 3 1 0 1 0 .
5 2 1 1 1
5 5 1 2 1
81. c (1,4,1,1,1), b 28,2,12, A 1 2 0 1 0 .
3 4 0 0 1
1 2 1 0 0
82. c (0,8,2,1,1), b 2,20,6 , A 6 3 1 1 1 .
3 2 0 0 1
1 2 1 0 0
83. c (7,2,0,1,2), b 2,12,18, A 3 4 0 1 0 .
2 6 2 1 1
1 2 1 0 0
84. c (1,3,1,1,1), b 2,6,1, A 2 1 1 1 2 .
1 1 0 0 1
1 1 2 2 1
85. c (1,2,1,1,1), b 11,2,3, A 1 2 0 1 0 .
1 1 0 0 1
2 3 1 2 1
86. c (10,5,2,1,1), b 17,1,3, A 1 1 0 1 0 .
1 3 0 0 1
1 0 2 1 3
87. c (2,1,3,1,1), b 6,16,7 , A 2 2 4 8 4 .
1 0 1 7 1
187
2 6 1 1 1
88. c (1,11,1,2,1), b 13,10,1, A 2 5 0 1 0 .
7 4 0 0 1
1 1 1 0 0
89. c ( 2,1,1,1,1), b 2,11,3, A 1 1 2 2 1 .
1 1 0 0 1
188
Appendix 3
KEYS
189
To the 3-rd task
20
1. J (u* ) 5 ; 17. J (u* ) 15 ; 33. J (u * ) ;
3
4 22
2. J (u* ) 6 ; 18. J (u* ) ; 34. J (u* ) ;
3 3
J (u* ) 14 ; 437
3. J (u* ) 6 ; 19. 35. J (u* ) ;
13
14 43 31
4. J (u* ) ; 20. J (u* ) ; 36. J (u* ) ;
3 7 3
134 99
5. J (u* ) ; 21. J (u* ) ; 37. U ;
7 5
137 17 89
6. J (u* ) ; 22. J (u* ) ; 38. J (u* ) ;
33 3 5
7. J (u* ) 0 ; 23. J (u* ) 7 ; 39. J (u* ) 19 ;
11
8. J (u* ) ; 24. J (u* ) 6 ; 40. J (u* ) 21 ;
3
9. J (u* ) 0 ; 25. U ; 41. J (u* ) 3 ;
27 17
10. J (u* ) 7 ; 26. J (u* ) ; 42. J (u* ) ;
5 3
54 110
11. J (u* ) ; 27. J (u* ) 16 ; 43. J (u* ) ;
13 7
19
12. J (u* ) 10 ; 28. J (u* ) 26 ; 44. J (u* ) ;
2
13. J (u* ) 6 ; 29. J (u* ) 12 ; 45. J (u* ) 5 ;
7 3 77
14. J (u* ) ; 30. J (u* ) ; 46. J (u* ) ;
6 7 5
118 389
15. J (u* ) 26 ; 31. J (u* ) ; 47. J (u* ) ;
5 13
29 19 72
16. J (u* ) ; 32. J (u* ) ; 48. J (u* ) ;
5 3 5
438
49. J (u* ) 16 ; 55. J (u* ) ; 61. J (u* ) 5 ;
13
190
133
50. J (u* ) 23 ; 56. J (u* ) 8 ; 62. J (u* ) ;
37
23 24 21
51. J (u* ) ; 57. J (u* ) ; 63. J (u* ) ;
3 5 10
52. J (u* ) 6 ; 58. J (u* ) 22 ; 64. J (u* ) 1 ;
5
53. J (u* ) 4 ; 59. J (u* ) 28 ; 65. J (u* ) .
17
26
54. J (u * ) ; 60. U ;
5
191
Appendix 4
TESTS
* * * uU
A) Weierstrass’ theorems
B) Farkas’s theorems
C) Kuhn-Tucker’s theorems
D) Lagrange theorems
E) There isn’t any correct answer amongst A)-D).
192
D) u U,v U, number 0,1 for which
u u 1 v U
E) u U , v U , E the point 1
u u 1 v U
D) I u , u u 0, u U
* *
I u I v I u 1 v ,
E)
u U , v U , 0,1
8) Function I u determined on convex ensemble U E n is
identified strongly convex on U , if
A)
0, I u 1 v I u 1 I v 1 u v 2,
u U ,v U , 0,1
I u 1 v I u 1 I v ,
B)
u U , v U , 0,1
I u 1 v I u 1 I v ,
C)
u U , v U , 0,1
D)
0, I u 1 v I u 1 I v 1 u v 2 ,
u U , v U , 0,1
E) There isn’t any correct answer amongst A)-D).
194
C) I u , 0, u U , E n
D) I u u , v 0, u U , v U
E) I u u , v 0, u U , v U
195
C) Nondegenerate problem of linear programming
D) Problem of nonlinear programming
E) The simplest variational problem
I u c, u inf
13) Define type of the following problem
n
u U {u E / u 0, i 1, n, Au b}
j
A) General problem of linear programming
B) Canonical problem of linear programming
C) Nondegenerate problem of linear programming
D) Problem of nonlinear programming
E) The simplest variational problem
20) Choose correct statement for problem
I u inf, u U E n , U {u U / I u min I u }
* * * uU
197
A) If function I(u) is semicontinuous from below on compact ensemble
U, that ensemble
U
*
B) If ensemble U is convex, but function I u is continuous on U,
that ensemble
U
*
C) If ensemble U is limited, but function I u is convex on ensemble
U, that ensemble
U
*
D) If function I u is semicontinuous from below on ensemble U, but
ensemble U is limited, that
I u
U *
198
ensemble
U {u U / I u min I u }
* * * uU
consists of single point
U {u U / I u min I u }
E) define necessary and sufficient conditions that ensemble
* * * uU
23) For the convex programming problem of the type
I u inf, u U {u E n / u U , g u 0, i 1, m}
0 i
Lagrange’s function has the form
A) L u , I u m
g u , u U ,
i 1 i i 0
m
{ E / 0,..., m 0}
0 1
B) Lu, m g u , 0,
i 1 i i i
m
i 1, m, 1, u U
i 1 i 0
C)
Lu , I (u ), u U E1
0
D) m
Lu , I (u g u ), u U , E m
i 1 i i 0
m
E) Lu , I u g u , u U ,
i1 i i 0
{ E m / 0,..., m 0}
0 1
24) For the convex programming problem of the type I(u) →inf
u U {u E n / u U , g u 0, i 1, m,
0 i
g u a i , u b 0, i m 1, s}
i i
Lagrange’s function has the form
199
A) Lu, I u s g u , u U ,
i 1 i i 0
{u E s / u 0, u 0, ..., u m 0}
0 1 2
B) m m
Lu, g u , 0, i 1, m, 1, u U
i 1 i i i i 1 i 0
m
L u , I u g u , u U , E m
C)
i 1 i i 0
s
D) Lu , I u g u , u U ,
i 1 i i 0
{u E s / u 0, u 0, ..., u s 0}
0 1 2
E) Lu, I u , u U , E1
25) For the convex programming problem necessary and sufficient
conditions that for any point
* * * uU
u U {u U / I u min I u
exist Lagrange’s multipliers
* such that pair
(u , ) forms
0 * *
saddle point to Lagrange’s function are defined
A) Kuhn-Tucker’s theorems
B) Lagrange’s theorems
C) Weierstrass’ theorem
D) Farkas’ theorems
E) Bellman’s theorems
condition
A)
I u , u u 0, u U
* *
B)
I u , u u 0, u U
* *
200
C)
I u I u I u , u u , u U
* *
D)
*
I u , 0, E n
I u 0
E)
*
27) Pair
u* , * U 0 0 is identified saddle point to Lagrange’s
function s , if
Lu , I u g u
i1 i i
A)
*
L u , L u , * L u , * , u U ,
* 0 0
B)
*
L u , L u , , u U ,
* 0 0
C)
* *
L u , 0
D)
*
L u , L u , * L u , , u U ,
0 0
E)
* * * * uU
u U {u U / I u min I u }, * 0, i 1, s
j
28) If pair
u*, * U 0 0 is a saddle point to Lagrange’s
201
29) For the nonlinear programming problem
I u inf, u U {u E n / u U , g u 0, i 1, m, g u 0, i 1, m
0 i i
generalized Lagrange’s function has the form
A) L u ,
s
I u g u , u U ,
0 i 1 i i 0 0
{ E s 1 / 0 , 0 ,..., m 0}
0 1
B) m
Lu , I u g u , u U , E m
i1 i 0
m
Lu , I u g u , u U , E m
C)
0 i1 i i 0
s
Lu , I u g u , u U , E s 1
D)
0 i 1 i i 0
s
Lu , I u g u , u U , 0
E)
i 1 i 0
B) necessary and sufficient condition that the point
u U {u U / I u min I u }
* * * * uU
u U {u U / I u min I u }
C) sufficient condition that the point
* * * * uU
D) necessary and sufficient condition that the function I(u) is convex in
the point
u U
*
E) necessary and sufficient condition that ensemble
U {u U / I u min I u }
* * * uU
202
31) Indicate faithful statement
A) convex programming problem is a partial case of the nonlinear
programming problem
B) nondegenerate problem of nonlinear programming can be reduced to
the convex programming problem
C) convex programming problem is a partial case of the linear
programming problem
D) any nonlinear programming problem can be reduced to the convex
programming problem
E) nonlinear programming problem is a partial case of the convex
programming problem
33) What from enumerated methods can be used for solving of the
convex programming problem
A) Lagrange multipliers method
B) Simplex-method
C) Method of the least squares
D) Pontryagin maximum principle
E) Bellman’s dynamic programming method
34) What from enumerated methods can be used for solving of the
linear programming problem
A) Simplex-method
B) Method of the least squares
C) Pontryagin maximum principle
D) Bellman’s dynamic programming method
E) any method from A)- D)
203
D) only in the internal points ensemble U
E) in internal, border, isolated points ensemble U
U {u U / I u min I u }
* * * uU
contains single point u *
, that the
point is
A) an extreme point ensemble U
B) an isolated points ensemble U
C) an internal point ensemble U
D) an internal or extreme point ensemble U
E) extreme or isolated point ensemble U
204
D) I u c,u inf,
u U {u E n / u 0, j I , g u ai ,
j i
u b 0, i 1, m}
i
E) I u c,u inf,
u U {u E n / 0 g u b , i 1, m}
i i
39) The linear programming problem in canonical form
I u c,u inf, is
identified nondegenerate,
n
uU {u E / u 0, j 1, n, Au b}
j
if:
A) any point u U has not less than rang A positive coordinates
B) rangA = m, where A is a constant matrix of dimensionality
m n, m n
C) ensemble
U
* * *
{u U / I u min I u }
uU
D) rangA = m, where A is a constant matrix of dimensionality
m n, m n
E) any point u U has no more than rangA = m, positive
coordinates
205
CONTENTS
FOREWORD ..................................................................................................... 3
Introduction ........................................................................................................ 4
Lecture 1. THE MAIN DETERMINATIONS. STATEMENT
OF THE PROBLEM .......................................................................................... 4
Lecture 2. WEIERSTRASS’ THEOREM ......................................................... 12
206
У ч еб но е изда ние
MATHEMATICAL
PROGRAMMING
Учебное пособие
ИБ № 5074
Подписано в печать 16.03.11. Формат 60х84 1/16. Бумага офсетная.
Печать RISO. Объем 13,00 п.л. Тираж 500 экз. Заказ № 230.
Издательство «Қазақ университетi» Казахского национального
университета им. аль-Фараби. 050040, г. Алматы, пр. аль-Фараби, 71. КазНУ.
Отпечатано в типографии издательства «Қазақ университетi».
207