Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

лекции по спец. курсу(лубочкин)

.pdf
Скачиваний:
19
Добавлен:
25.03.2015
Размер:
464.24 Кб
Скачать

n #%&&' '"+ %"4+& ")'"'#'( ( ')3% % ! +& ( )$"4 *%5%, n 5,3'4 ( )3+/

%5 7, +%&"+(&'"#- 4,, '"' ")'"'#, )3+&+5, # (3%5! 3,#'# 5% -&+* &, 4,, '+ +7 n ('75', & & ')%3,*++

.#'7&, +5 %3%7 zj ( j = 1, n $4,7,&& % +%&"+(&'"#+ 6 :

 

 

! ')+",&+! 4,, '+ ')%3,*++ +")' -7$ #"! &'35,#+(&

% ,&& %8

 

αj

6j = 1, n: ( )$"4 *%5%, j/+ 5,34+ )3+ % +&+ &'+ +%&"+(&'/

"#+ 64'

+ %"#(' # " #'&& *%5%, j/+ 5,34+ ( ( )$"4,%5'"' 7, # " ,"'(

3,

#'#

5%

-&+* :

! 6)' (%"$: 4 +&4%3,(

 

 

 

 

 

 

 

γ1j

6j = 1, n: '

 

 

γ2j

6j =

1, n

: '

! "+)",(

 

 

γ3j

6j

=

1, n

: ' ! "3,& ,4, ( *%5%% j/+ 5,34+

6("% (% + +&

#%3$#"! +7 3%*%)#$3 *%5%, j/+ 5,34+:

 

a1j 6j = 1, n: 7,#3,# 4 +&4%3, )3+ % +&+ &'+ +%&"+(&'"#+ j/"'

")'"'#,(

 

 

 

 

 

 

= 1, n: 7,#3,# "+)",(

a2j

6j

a3j

6j =

1, n

: 7,#3,# "3,& ,4, 6(% + +& 3,"" +# (, #"! )' αj +

γij )' !'35$ ,5 aij = αj γij :

a4j

6j =

1, n

: "$55,3& % &'35,#+(& % 7,#3,# 6( "#'+5'"#&'5 ( 3,/

,%&++: %4#3' &%3"++ + 5% +&

#%

T 3%"$3" (3%5%&+ 5% -&+*

 

63'5% #'"' ( 7, ,& 8 '#)$"4& % *%& *%5%'( dj ( j = 1, n 7,#3,# &,

+")' -7$%5'% " 3-% pi ( i =

1, 4

64 +&4%3 ( "+)"( "3,& ,4 (

%4#3' &%3"+!:

'"3,&+ %&+!

"&+7$ )' ( )$"4$ q j ( j =

1, n

)3% % -& +

'# %5 ( )$"4,

4 +&4%3, r1

 

 

 

 

 

 

 

0(% %5 "

% $ +% )%3%5%&& %8 qj '# %5 ( )$"4, *%5%, j/+ 5,34+

(

j =

1, n

ri

 

7,#3,# "''#(%#"#(%&&' "+)",( "3,& ,4, +

%4#3' &%3"++

(

i= 2, 4

#,(+#"! 7, , , 5,4"+5+7,*++ '&' , *%5%&'"' 7,(' , /' '&' '5 ( ,&&'5 " $ ,% )'&+5,%#"! 3,7&'"#- 5%, $ '# %5'5 ( )$"4, ( "#'+5'"#/ &'5 ( 3,,%&++ + +7 %3, 4,5+ &, " 3-% $ %#'5 ((% %&& & '#'7&, %&++ 5,#%5,#+ %"4,! 5' % - #'+ 7, , + +5%%# (+

n

4

 

 

 

 

 

 

 

X

X

− p1r1 → max, qj = αj zj

 

 

 

dj qj

piri

≥ q j , j =

1, n

;

j=1

i=2

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

X

 

 

≤ r1 ;

X

 

ri =

aij zj , i =

1, 4; r1

zj ≤ T.

 

j=1

 

j=1

6. ':

"4 +5 +7 7, , + 6. ': )%3%5%&& % qj ( j = 1, n( + ri ( i = 2, 4 ! #'/

"' ('")' -7$%5"! ( 3,,%&+!5+ ! #+& )%3%5%&& & %3%7 )%3%5%&& % zj ( j = 1, n +7 #'+ 7, , + ( , #,4,% ((% %5 '#'7&, %&+% cj = dj αj P4i=2 piaij (

7

j =

1, n

0% + +&, cj '7&, ,%# '&' ( )' $ ,%5+ '# j/"' #%&&' '"+ %/

"4'

"' ")'"'#, )3+ %"' % +&+ &'+

+%&"+(&'"#+ 6#%7 $ %#, 7,#3,# &, +7"'/

#'( %&+% 4 +&4%3,: /'"4' -4$ (% + +&, p1r1

! +4"+3'(,&,( #' %% 5', &'

')$"#+#- + $ %"#- ( +#'"'('5 3%7$

-#,#% )'"

% 3% %&+! 7, ,

+ ,4+5 '#/

3,7'5 ( 7, , , 6.':

4(+(, %&, "

% $ %+

4"#3%5,

-&'+ 7, , %8

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

→ max,

 

 

 

X

 

 

 

≤ r1 ;

 

 

 

 

 

 

 

 

 

 

 

 

cj zj

 

 

 

 

 

 

a1j zj

 

 

 

6. 01:

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

zj ≤ T ; zj

, j =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, n.

 

 

 

 

 

 

 

 

 

 

j=1

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 5

 

0 )

. . #

+ )'"#3'%&

3,7

+

&

% 5' %

+

 

4'&'5+ %"4+

& 7, , ( )3% "#,(/

! +% "'

#'+ 3,7 + & % 7, ,

 

+

1 / 8

 

 

6

. )

:(

6

. 9:(

6. -

:(

6. .1:

.# 4,, '+

+7 &+& 5', &' )%3%+#+ 4

 

4(+(,

%&'

+

 

5' %

+ 3

$"'

"

' (+ , 0 ,"#&'"#+ (

'#

6. ):

( 6. -:( 6.

.1: 5'

, &'

)%3%

+#+ 4

 

 

4(+(,

%&'

+ 7, , %

1 / ( 4,&'/

&+ %"4'

+ ! '35%(

# %

4 7, ,

% (+ ,

6

. 9:

,""5'#3+5 3,7 + & % " $ ,+

%5%'( 7, , 1 / (

'#

+ &

&

'#

 

%5%'( 7, , + 6. 9:8

 

.: " + ( +"&' &'+ 7, , % *%

%($

 

 

!$&4*+ &$, &' 5+&+5+7+3'(,#-

 

 

 

 

 

 

 

 

 

 

 

c 0 x → min,

 

 

 

 

 

 

 

 

 

 

#' )%3%+#+ 4 7, , % &, 5,4"+5$5 5', &' "

% $ +5 '#3,7'5 8

 

 

 

 

 

 

 

 

 

 

−c 0 x → max.

 

 

 

 

 

 

 

 

 

 

0(' ! &'( % ),3,5%#3 c˜ = −c( )'

$ +5

%5% 7, , + 6. 9:8

 

 

 

 

 

 

 

 

 

 

 

0 x → max.

 

 

 

 

 

 

 

 

 

 

: "

+ ( 3%7$

-#,#% )'"#3'%&+! 5' %

+ )'

$ %& (, '"3,&+ %&+! (+ ,

 

 

 

 

 

 

 

 

 

A1x = b1, A2x = b2,

 

 

 

 

 

 

#' ( ((' ! &'( % 5,#3+*$ + (%4#'3

, b = b1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

b2

 

 

 

 

 

 

 

#+ '"3,&+

%&+! 5',

&' 7,)+",#- ( (+ % ' &'"' '"3,&+ %&+!

 

 

 

 

 

 

 

A1

x =

 

b1

,

 

 

 

# %

Ax = b

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

 

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6+ ( #'5 "

$ ,% )' $ %& %5% 7, , + 6. 9::

 

 

 

 

 

 

 

77

0: " + )' $ %&' '"3,&+ %&+% ( (+ % &%3,(%&"#(,

 

 

n

 

X

 

 

Ax ≤ b, # %

aij xj ≤ bi, i =

1, m,

 

 

j=1

#' ((' !# &'( % ')' &+#% -

& % )%3%5%&& % xn+i ≥ 0( i =

1, m

( + '#,/

( ! # +& 6)' ' &'+: 4 '"3,&+

%&+!5/&%3,(%&"#(,5 8

 

 

n

X aij xj + xn+i = bi; xn+i ≥ 0, i = 1, m.

j=1

'" , &'( % )%3%5%&& % )3+&+5, # 7&, %&+!

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

xn+i = bi

aij xj ≥ 0, i =

1, m,

 

 

 

 

 

 

 

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

+ '7&, , #( &,)3+5%3 ( (%

+ +&$ &%+")'

-7'(,&&'"' 3%"$3", Ri ( i =

1, m

 

0 (%4#'3&'/5,#3+ &'+

! '35% ( %

# + )3' % ,&

"

% $ +%

%+/

"#(+! 8 #

((% %& ')' &+#%

-& + m/(%4#'3 x = (xn+i,

 

i = 1, m)

≥ 0(

4'#'3 + #

'#,( %& 4

%('+ ,"#+ '"3,&+ %&+!/&%3,(%&"#(,8

 

 

 

Ax + x = b, x ≥ 0.

0(' ! % +&+ &$ m × m/5,#3+*$ E = diag {1, i = 1, m}( m × (n + m)/

5,#3+*$

˜

 

(

 

 

/(%4#'3

 

 

 

 

(

'"3,&+ %&+! 7,)+ %5 (

 

A = [A, E]

 

(n + m)

 

 

 

x˜ = (x, x )

 

 

 

 

(+ %

 

 

 

 

 

˜

 

 

≥ 0

 

 

 

 

 

 

 

 

 

 

Ax˜ = b,

 

 

 

 

 

6# % + (

#'5 "

$ ,% )' $

%&

 

%5%

7, , + 6

. 9::

 

 

&, '"+ &' (

#' -4' ( +#,!

')'

&+#%

-& % )%3%5%&& %( )'"#$),%5 (

" $ ,% '"3,&+ %&+! (+ ,

Ax ≥ b

8 )' ,",%5

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

A = [A, −E]

 

 

): %3(,

-& % '"3,&+ %&+! b ≤ Ax ≤ b 7,)+ %5 ( (+ % ($& '"3,/

&+ %&++ 8 Ax ≤ b ( Ax ≥ b 0(' ! ')'

&+#% -& % )%3%5%&& % x(1) ≥ 0(

x(2) ≥ 0( "'" ,"&' " $ , 0: )'

$ +5

 

 

 

 

 

 

 

 

 

 

 

Ax + x(1) = b ,

Ax − x(2) = b .

 

 

0(' ! 5,#3+*$ + (%4#'3

 

 

 

 

 

(1)

 

 

 

b

 

 

 

A

E

0

 

 

 

 

 

 

 

 

A˜ = A

0

−E

, x˜ =

x

,

 

˜b = b

,

 

x(2)

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"'" ,"&' " $ , : )' $ +5 '"3,&+ %&+!

˜

˜

(1)

≥ 0 x

(2)

≥ 0.

Ax˜ = b, x

 

 

7

: ($&"#'3'&&+% '"3,&+ %&+! d ≤ x ≤ d

)%3%)+ %5 ( (+ % x ≤ d (

x d ( # % Ex d ( Ex d , %% )'"#$),%5 )' "&%5% " $ ,! ):

 

9:

 

 

 

 

 

 

)%3%5%&&(1)

&

&% &,

'

,

%&' '

"

3,&+

 

%&+% &%'#3+*,#%

-/

 

" + &,(1) ,"#(2)-

 

 

 

 

&'"#+8 x = (x

, x

)( x

≥ 0 6 #' 5',%# " $ +#"! ( &,)3+5%3 ( )3+ ( )'

/

&%&++ )3%'

#

 

 

 

+ (

 

&

& )

 

(

 

,

 

 

 

3,7'(,&+(2) 7,5%&!%5')+",&&3,7&'"#(- " $($,!& &%'#3+*,#%: + :: #'-4,& & )%3%5%&&$ )%3%5%&/& 8

&$ (%4#'3, x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xj(2) = xˇj(2) − xˆj(2)

( j(2) ≥ 0( j(2)

≥ 0( # %

x(2)

= xˇ(2) − xˆ(2) ( (2)

≥ 0( (2) ≥ 0

/'" % &%% 3,(%&"#(' 5', &' )%3%)+",#- ( (+ % Ex(2) − Exˇ(2) + Exˆ(2) = 0

, %% )'"#$),%5 ,&, '"+ &' )3% $ %5$

7

 

 

0

,&&'+ " ,(% '#'"&'( (,%#"! ' +& +7 "'(3%5%&& & !!%4#+(&& 5%#'/

'( 1 / , ,)#+(& + 5%#'

 

2

 

4

 

2

 

4

 

 

 

4 4

 

,""5'#3+5 4,&'&+ %"4$ 7, , $ 1 / " ($&"#'3'&&+5+ )3!5 5+ '"3,&+/

%&+!5+

0

x → max,

6 .:

 

f (x) = c

 

Ax = b,

6 :

 

d ≤ x ≤ d .

6 0:

% %"-8 c = c(J ) = (cj , j J ) = (c1, c2, . . . , cn)( x = x(J ) ( d = d (J )( d = d (J )n/(%4#'3 b = b(I) = (bi, i I) = (b1, b2, . . . , bm) m/(%4#'3 A = A(I, J ) = (aij , i I, j J ) m × n/5,#3+*,( m < n I = {1, 2, . . . , m}( J = {1, 2, . . . , n}

5&',%"#(, +& %4"'( ),3,5%#3'( b1 ( b2 ( ( bm + )%3%5%&& & 7, , + x1 ( x2 ( ( xn "''#(%#"#(%&&' $ %5 " +#,#-( #' rankA = m

+"

, m + n &,3,4#%3+7$ # 3,75%3 7, , + &',%"#(' {c, A, b, d , d }

'#3,7$%# "'('4$)&'"#- %% ),3,5%#3'( ,"#' +5 )3+"(,+(,

# ")%*+, -& %

&,7(,&+! 8 c (%4#'3 "#'+5'"#+ ( b (%4#'3 3%"$3"'( ( A 5,#3+*, 7,#3,#

6%% "#'

#* aj = A(I, j) ( j J (%4#'3 $" '(++:

 

&'

,%"#(' ')$"#+5 & (%4#'3'( 7, , + 6 .:2 6 0:

 

 

X = {x Rn: Ax = b, d ≤ x ≤ d }

 

7, ,%#"! "+"#%5'+ +7 m +&%+&& $3,(&%&++ 6 : 6'"&'(&'% '"3,&+ %&+%:

+ 2n

+&%+&& &%3,(%&"#( 6 0: 6)3!5 % '"3,&+ %&+!:

%'5%#3+ %"4+

5&',%"#(' X )3% "#,( !%# "'#'+ )%3%"% %&+% 4'&% &'/

"'

+" , "+)%3)

'"4'"#%

+

+ )'

$)3'"#3,&"#( ,4+%

5&',%"#(, &,7 (,/

# 5&'"'"3,&&

5+ ( , (

"

$ ,%

+& '"3,&+ %&&'"#+

5&'"'"3,&&+4,5+

6"5 3+" .( :

 

 

 

 

 

 

 

 

 

7

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

 

 

 

 

 

 

 

 

7

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7 7

x

0

, "%'5%#3+

%"4'5 !7

4% 7, , , 6 .:2 6 0: "'"#'+# ( )'+"4% #,4'+

#' 4+

5&'"'"3,&&'"' 5&',%"#(, X ( ( 4'#'3'+

 

+&%+&,! ! '35, 6 .: '"#+",%#

5,4"+5,

-&'"' 7&, %&+! 6"5 3+" 0:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

f (x) = c 0x0 = max c 0x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

!!

!!

!!T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

6

 

f (x) = c 0x = α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BBB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

%5% 5&', %"#(, X &,7 (, # ) ,&,5+ ( # % ) ,& #'

#'+ n/

(%4#'3 x( $ '(

%#('3! ++ '"&'(&'5$ 6 : + )3!5 5 6 0: '"3,&+ %&+!5

7, , +

x

0 ( &,

4'#'3'5

*%

%(,! !$&4*+!

 

6 .: '"#+",%# 5,4"+5,

-&'"'

 

/ ,&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7&, %&+! (

&,7 (, # 3% %&+%5 7, ,

+ + + %% ')#+5,

-& 5 )

,&'5

#')#+/

 

 

, )3,4#+4% 43'5% ')#+5, -&&

) ,&'(ε

(,, &$

3' - +"

3,

 

# "$

5,

-& % 6ε/')#+5, -& %: )

,& /

,& x

)3+ 7, ,&&'5 +"

% ε

≥ 0 &,7 /

7

(, # ε/')#+5,

-& 5 ( %"

+ '& $ '( %#('3!%# &%3,(%&"#($

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x0) − f (xε) ≤ ε.

 

 

 

6 ):

&',%"#(' "$#')#+5,

 

-& & ) ,&'( +7'#3,,%&' &, 3+" )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

 

!

 

x0

 

 

 

 

f (x) = c 0x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

 

!

!

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!

!

 

 

 

 

 

 

X ε

 

T

 

T

 

 

 

f (x) = c 0x0

ε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TTTT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BBB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

2

2

4 4

 

 

 

 

 

'

"

,"&' 7,)+"+ 7, ,

 

+

 

6

.:2

6

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0: '"&'(&'0 #' (')3'"4+ 73%&+!)3+ )3,4#+4+%% 3% %&++(,,"'"#'/& 5

+# ( )'"#3'%&++ ')#+5,

 

 

-&'"' )

 

,&, x

 

 

 

 

 

!( !%#"! #,4,% )'+"4 "$#')#+5, -&'"' )

,&, xε

" ",3,+3'(,&&'+ '*%&/

4'+ 4, %"#(, ε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/3+ +""

% '(,&++ 7, , 1 / ('7&+4, # + 3$"+% (')3'" ( (,, & % !

)3+

'

,%&+

+ 8

 

 

 

 

 

 

 

 

+ '"3,&+ %&+! 7, , + ( # % +5%%#"!

+ &'#! # ' +&

1) &% )3'#+('3% +(

 

 

 

 

 

)

,&

6,4 %"' )'"#3'+#-

6)3'# %5 )%3('+ ! ,7 :

 

 

 

2) %"#-

+ "3% + '"3,&+ %&++ +7##' & %( +"4 %&+% 4'#'3& &% ( +!%#

 

 

&, 3%7

$

-#,#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%(,! !$&4*+!

 

 

 

3) '"3,&+ %&,

+ &, 5&',%"#(% ) ,&'( *%

 

 

 

4) !(

!%#"! + ')#+5,

 

-& + )

 

 

 

,& % +&"#(%&& 5

 

 

 

5) 4,4 +75%&+#"! ')#+5,

-&+ )

 

,& + ')#+5, -&'% 7&, %&+% *%

%('+ !$&4/

 

 

*++ )3+ +75%&%&++ ),3,5%#3'( 7, , +

6)3'#

%5, $("#(+#% -&'"#+:

6) 4,4 4'&'5&' "4'33%4#+3'(,#- ')#+5,

-&+ )

,& )3+ &%#'

- +& +75%/

 

 

&%&+!& ),3,5%#3'(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

$

#

')#+5, -& % ) ,&

 

 

5',

&' "#3'+#- 3,7 + & 5+ 5%#' ,5+

6+7 3,7 + /

&&

4

,""'(: .)+" (,%5 + ,

 

%% 5%#' )3+&,

%, +# 4

,""$ 8

 

1) )3!5& (

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

2)#' & & (

3)3% ,4",*+'&& &

1)

#' '7&, ,%#( #'

 

 

 

 

 

 

 

 

 

 

 

( )3'*%""% 3,#'# 5%#' , )3%'#3,7$%#"! +&!'35,*+! '# +"&' &'+ 6, &%

2)

4,4'+/ +#' (")'5'",#% -&'+: 7, , %(

 

 

 

 

 

 

 

 

 

k (

( )3'*%""% %"' 3,#'# "%&%3+3$%#"! )'" % '(,#% -&'"#- (%4#'3'( x

 

 

k = 1, 2, . . . ( #' &' 6, &% )3+#

+,%&&': $ '(

%#('3! +& ("%5 '"3,&+/

3)

%&+!5 7, , + (

&!%#"! &%3,(%&"#(' c

0

x

k+1

≥ c

0

x

k

k ( " %

'# +#%3,*++ 4 +#%3,*++ ( )'

 

 

 

 

 

k &'5%3 +#%3,*++ 5%#' ,

 

 

 

 

 

 

 

 

 

 

 

 

* 4,7,&& % "('+"#(, (,, &

" #' 4+ 73%&+! )3+4 , & & 7, , ,#%5,/

#+ %"4,! 5' %

- 6 .:2 6 0:(

4,4 )3,(+ ' ( &%)'"3% "#(%&&' "(!7,&, " )3,4/

#+ %"4'+ 7, , %+ /' #'5$

 

1)

#,! +&!'35,*+! ' &%+ 5', %# # #- +")' -7'(,&, ")%*+, +"#,5+ ,&/

2)

&'+

)3,4#+ %"4'+ '# ,"#+ (

"'#

%&+% ( )3'*%""% 3,#'# 5%#' , ("%& '"3,&+ %&++ )'7(' !%# ")%/

 

*+,

+"#,5 (

#'+ 5'5% +5%#- #' &$ 4,3#+&$ +"" % $%5'+ )3,4#+/

%"4'+ 7, , + 6)3+ &,3$ %&++ '"3,&+ %&++ '"#$)&,! +&!'35,*+! 5'/

,%# &% +5%#- &+4,4'"' )3,4#+ %"4'"' "5 " ,:(

3)"%&%3+3'(,&+% $ $ , %+"! #' &'+ +&!'35,*++ ' )3,4#+ %"4'+ 7, , %

)'7('

!%# &, %, &' 4'' +3'(,#- )3'*%"" 3%

%&+!

 

 

/' #'5$ '"#,&'( )3'*%"", 3% %&+!

6)'

# 5 )3+

+&,5: )'7(' !%# (

#,4'5 4

 

,""% 5%#' '( )'

$

+#- 4'5)%&",*+

7, +73,"&' '(,&& % 3%"$3"

0 3$"

+

& 5%#' ,& #,4'+

4'5)%&",*++ 5',%# + &%

# #- (

# 4

1) '"#$)&,! 4 5'5%$ '"#,&'(, +&!'35,*+! 5', %# &% +5%#- &+4,4'"'

'#&' %&+! 4 7, , % 6 .:2 6 0:(

 

 

 

 

 

 

 

 

2) #, +&!'35,*+! 5', %# "

+ 4'5 "3$#' ')+" (,#- "+#$,*+ (

3) 4, %"#(' +&! '35,*++ 5',%# '4,7,#-"! &$, % +"&' &'+

6'&% &' ( '"#+,%&+%

# & )' %7&

&

"('+"#( '# &' "')3!,%&' "

')3% %

(

%&& 5+ 3,"&' ,5+

%#' ( '#

, , +% ')3% %

%&& 5+ '"#'+&/

"#(,5+

 

5'"$# ( 3$"+&

'#&' %&+!&

)3'+"3 (,#- 5%#' ,5 ( #+5+ '"#'/

+&"#(,5+ &% '# , , +5 ,4 ( )3!5

%

( #' & %(

3% ,4",*+'&& % 5%#'

5'"$# 3% ,#- 7, , $ '

- %( %5 &%)3!5 %

( )3+#

+,%&&

% + &%3% ,4",*+/

'&& % 0 #'3 4 ,"", 5%#' , 3% %&+! 7,(+"# '# 4'&% &'+

*% + 3,(&+(,#-

)' !!%4#+(&'"#+ 5', &'

+ - 5%#'

 

' &'"' 4

,"", 6/3!5'+ "+5) %4"/

5%#'

, ,&*+", '#&'"+#"! 4 #'5$

 

,% 4

,""$:

 

 

 

0 "+ $ ( #3,&&'"' 4

,"", 5%#' , 3% %&+! #

$ %5 " +#,#-( #' ( "'"#,(

+"&' &'+ +&! '35,*++ (&' +# &%4'#'3 + 6&, ,

-& +: )

 

,& x

7

2

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/3+ 7, ,&&'5 ) ,&% x 5%#' ' ,%& +5%#-

"3% "#(,

3,")'7&,(,&+!

%"'

')#+5, -&'"#+ + + ε/')#+5,

-&'"#+ " + )

,& '4,,%#"! &% ε/')#+5,

-/

& 5 ( #' &$, &' $5%#- )'"#3'+#- 3$"'+ )

,& x( 4'#'3 + &% &$,% +"&' &'"'

 

#'# 3%

+#- )%3% +"

%&& % 7, , + (

&$, &' $5%#- !!%4#+(&' "3,(&+/

(,#- ) ,& x " 3$"+5+ )

 

,&,5+ ( # % "#3'+#- (%4#'3 ( $ '(

%#('3! +%

("%5 '"3,&+ %&+!5 7, ,

+

 

"&' ( #'

"

,(&,! " ', &'"#- "'"#'+# ( "'#

/

%&++ '"&'(&'"' '"

3,&+ %&+!

6 :

( )'"4' -4$ &,

+#+ (%4#'3 ( $ '(

%#('3!/

 

+

+

#'

-4' )3!5

5 '"3,&+

%&+!5

6 0:

( &% "'"#,( !%# &+4,4'"

' #3$ ,

,4+5 '#

3,7'5 (

)%3('+

)3+ 3,73,#'#4% )3!5

&

5%#' '( 3% %&+! 7, ,

+

6 .:2

6 0: ('7&+4,%# )3'

#

 

%5, $)3,(

%&+! '"&'(&

5 %% '"3,&+

%&+%5

 

!

!!%4#+(&'"

' 3,73%

%&+!

#'+ )3'#

%5 ,

%% )3% ,",%#"! ")%*+,

-&'%

"3% "#(' (

&,7

(,%5'% ')'3'+

 

 

6

.:2 6 0: &,7'(%5 $)3,( !%5

 

5 (

%"/

 

."&'(&'% '

"3,&+

%&+%

 

6 : 7, , +

+

!

 

# & m/(%4#'3, b + n/(%4#'3, x &,+ %#"! #,4'+ n/(%4#'3

 

x( #'

&, (%4#'3% x = x +

x #$ %# ( )'

&!#-"! '"&'(&'% '"3,&+ %&+%8 Ax = b

 

7 )'"

% &%"' 3,(%&"#(, +5%%5 Ax = A(x+

x) = Ax+A

x = b( '#4$ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A x = b − Ax.

 

 

 

 

 

 

 

 

6 :

,4+5 '#3,7'5 (

7, , , $)3,(

!%5'"#+ '"&'(& 5 '

"3,&+ %&+%5 "(%

,"- 4

)3'#

 

%5% 3,73% +5'"#+ '#&'"+#%

-&'

x "+"#%5

6 :( "'"#'! %+

+7 m

+&%+& & , "%#3,+ %"4+&

$3,(&%&++ " n &%+7(%"#& 5+

xj

( j J

7

4$3",

 

+&%+&'+ , "%

#3

+7(%"#&' ( #' "+"#%5,

6 : 3,73% +5,

6)3+

 

# &

)3,(&

 

,"#!& : #'" , + #'

-4' #'" ,( 4'" , 5,#3+*, A )'

&'"' 3,&",( # %

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rankA = m.

 

 

 

 

 

 

 

 

6 9:

 

,4+5 '#3,7'5 ( '4,7,&,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

643+#%3++

$)3,(

!%5'"#+ '"&'(&'"'

'"3,&+ %&+!:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7 " % ,&&&

)3% )'

',%&++ )3+ )'"#,&'(4% 7, , + 6 .:2 6 0: "

% $%#(

#' %% '"&'(&'% '"3,&+ %&+% $)3,(

!%5'

 

 

 

 

 

 

 

 

 

3,&+/

 

6,4 " % $%# +7 #%'3%5

.( ( "

$ ,% $)3,( !%5'"#+ '"&'(&'"' '"

%&+! 6 #' + )3% )'

,",%#"!: ( 5,#3+*% A 5', &' ( % +#- "+"#%5$ +7 m

+&%

+

&' &%7,(+"+5

& "#'

 

#*'( /$"#- #,4,! "+"#%5, (

#3,&, &'

,%"#('

+& %4"'( #+& "#'

#*'( '#'7&, +5 %3%7 J ( , "''#(%#"#($ $ ( 3%74

$

+7 5,#3+* A %3%7 P = A = A(I, J ) = (aj = A(I, j), j J ) '" ,

+5%%5

detP

6= 0.

6 $:

 

 

7

 

 

&',%"#(' +& %4"'( J J ( |J | = |I| = m( &,7'(%5 ')'3'+ ( %" + !

')'3&'+ 5,#3+* P ( )' &!%#"! $"

'(+% 6 $:

 

 

 

 

 

 

 

 

 

 

 

* +#

(,! )'"

% &%% ')3% %

%&+% ( )' $ +5

 

 

 

 

 

 

 

 

 

 

 

 

2

2

6(#'3,! !'35, 43+#%3+! $)3,(

 

!%5'"#+ '"&'(&'"' '"3,/

&+ %&+!:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%#

7 " % ,&& &

)3% )' ', %&++ "

% $%#( #' ( 7, , % 6 .:2

6 0: "$ %"#($/

&'#!

# ' &, ')'3,

 

 

&'"#+ ')'3

 

"'"#'+# ( #'5

( #' " %% )'5'/

 

3

$#'% )3% "#,(

 

%&+% ' "$

-

5'

, &' $)3,(

 

!#- '"&'(&

5 '"3,&+ %&+%5

' %% #'&4++ "5 " ')'3

#$ %# , %% )'"#%)%&&' ( !( !#-"!

 

5%#' '( '"$ %"#(

 

!%#"! )%3%&'

 

/'"4' -4$ &, 4,

, '+ +#%3,*++ )3!5 &

 

'# )

,&, x 7, , + 4 &%4'#'3'5$ 3

$"'5$ )

,&$ x = x +

 

 

 

x( #' )3+ 3,7/

3,

#'#4% )3!5 &

5%#' '( ('7&+4,%#

 

,"#& +

" $

,+

$)3,(

 

%&+! '"&'(& 5

'"3,&+ %&+%5 ( , +5%&&' ( )'"4'

-4$ x )

,& 7, , + ( #' Ax = b( + +7 6 :

)'

$ ,%5

 

 

 

 

 

 

 

 

A

 

x = 0.

 

 

 

 

 

 

 

 

 

 

 

6 -:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

")'

-7$! ')'3$(

3,(%&"#('

6 -: 7,)+ %5 ( (+ % P

x + A

x = 0(

" % A

 

 

 

)(

 

 

 

 

 

 

 

x

 

x(J )( J

 

 

 

 

 

 

 

= A(I, J

 

x

=

x(J ) (

=

= J \ J ( '#4$ ,

" % $%#( #'

 

 

 

 

 

x = −QA x , Q = P

−1

.

 

 

 

 

 

 

6 ':

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,4+5 '#3,7'5 ( 4,4+5 #

&+ (7!#- (n − m)/(%4#'3

x

 

( )' " +#,( )'

! '35$ % 6 ': m/(%4#'3

x ( )' $ +5 n/(%4#'3

 

x = (Δx ,

x

)( '#,(/

%&+% 4'#'3'"' 4

 

#'5$ ) ,&$ x ,%# n/(%4#'3 x = x +

 

x( !(

! ++"!

)"%( ') ,&'5 ( # % (%4#'3'5 (

&, 4'#'3'5 ( )' &!%#"! '"&'(&'% '"3,&+ %/

&+%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%,5%#+5 #,4,%( #' ( #'3'5 4'5)'&% x = x(J ) (%4#'3, x )' ! '3/

5$ %

 

 

 

 

 

 

 

 

x = Q(b − A x )

 

 

 

 

 

 

 

 

 

6 .1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)3+

#'5 ( #'3% (%4#'3, b + '"#,

-& & 4'5)'&% x

= x(J ) ( 5', &'

'#%")% +#- ( )' &%&+% &, (%4#'3% x '"&'(&'"' '"3,&+ %&+! ( # % (%4#'3 x

#$ %# )"%( ') ,&'5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,&,5+ /'/

 

0 )3'*%""% 3,

#'# 5%#' , ')'3, #$ %# 5%&!#-"! &,3! $ " )

#'5$ '"&'(& 5 '#

%4#'5 )3%'#3,7'(,&+! "#,&'(+#"! ),3, +7 )

,&, + ')'/

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,&'5 +#,%5 %"'

 

/,3$ {x, J } +7 ) ,&, + ')'3 &,7'(%5 ')'3& 5 )

&%( 3', %&& 5 ( %"

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d j < xj < dj , j J .

 

 

 

 

 

 

 

 

 

6 ..: