Стохастикалық программалаудың өзара екі жақтылық есебі

Математиканың көптеген салаларында екі жақтылық теоремалары деп аталатын теоремалар кездеседі. Олардың әрқайсысы берілген теорияның кез-келген тұжырымы белгілі-бір ережелер стандарты бойынша екінші бір тұжырым алуға мүмкіндік береді және біріншісінің дұрыстығынан автоматты түрде екіншісінің де дұрыстығы алынады. Екі жақтылық есебінің мысалдары стохастикалық бағдарлау есебінде де кездеседі.

Өзара екі жақтылық есебін сипаттайық. Бастапқы немесе І-ші есепдеген атпен  стохастикалық бағдарлаудың қандай да бір болмасын есебі берілсін де шектеулер теңсіздіктер түрінде болсын. Сонымен

a11х112х2+…

а211222+… 

п белгісізді т теңсіздікте жүйесі мен сызықты форма

f= с1х12х2+… + с„х„(2)

берілсін.

(1)-ші жүйенің барлық теріс емес шешімдерінің ішінен (2)-ші форманы минимумдейтінін табу керек. Қойылған есеппен екінші бір есепті, яғни бастапқы есептен екі жақтылық байланыста болатын есепті қарастырайық. Оның қойылуы мынадай: т — белгісізі бар п теңдеулер жүйесі берілген:

а х  + а х  + … + а  х  ≤ b

…   …   …   …   …   …   … …

a x  + a x  + … + a x  ≤ b

a x  + a x  + … + a x = b

…   …   …   …   …   …   …  …

a x + a x  + … + a x  = b

және стохастикалық форма

 

cр = Ь1у12у2+… + Ьтут

түрінде жазылады.

(Г) жүйесінің барлық теріс емес шешімдерінің ішінен    (2′) формасын

максимумдейтінін таңдап алу керек. Бұл есепті (П) — есеп деп алайық.

Екі есепті салыстыру арқылы мынаны байқауға болады.

1. Бастапқы  есептегі  белгісіздер жанындағы  коэффициенттерден құрылған  матрицасы мен екі жақтылық есебіндегі матрицалары бір — бірін «транспонирлеуден», яғни жатық жолды бағаналармен ауыстырудан алынады. Олардың реттері сақталады.

  1. Әр есептің шектеулер жүйесінің оң жағында екінші есептің сызықты формасының коэффициенттері орналасқан.
  2. Бастапқы есептегі теңсіздіктер > — типті, і- формасының минимумы
    табылу керек те, ал екі жақтылық есебіндегі шектеулерде < — типті теңсіздіктер,
    яғни қарама-қарсы мағыналы, ал  ср     формасы максимум мэн, яғни қарсы
    мағыналы мэн қабылдайды.
  3. Біреуінің белгісіздерінің саны, екіншісінің шектеулер санына тең.
  4. Осы бес шартты қанағаттандыратын стохастикалық бағдарлау есебі симметриялы деп аталады. Оның біреуі негізгі, ал екіншісі — екі жақтылы есептер деп аталады.
  5. Стохастикалық  бағдарлау есептерінде симметриялы екі жақты парлармен қатар, симметриялы емес те екі жақты парлар кездеседі. Онда:
  6. а) Негізгі есеп: Теңдік сақталады:

а х + а х + … + а  х  ≤ b

…   …   …   …   …   …   … …

a x + a x + … + a x  ≤ b

a x  + a x  + … + a x = b

…   …   …   …   …   …   …  …

a x + a x  + … + a x = b

б) Екі жақтылық есебі:

а у  + а у  + … + а  у  ≤с

…   …   …   …   …   …   … …

a у  + a у  + … + a у  ≤ с

a у  + a у  + … + a у = с

…   …   …   …   …   …   …  …

a у + a у  + … + a у  = с

Бұл   есептер   симметриялы   пар   жағдайынан   мынадай   екі   жолмен ерекшеленеді:

  1. (3) — (4)-ші есептердегі шектеулер теңдікпен берілген.
  2. (3′)-(4′)-есептерінде  у{   айнымалыларының  і = 1,т) теріс еместік

белгілері жоқ.

2. Сонымен қатар, екі жақтылық есебіндегі стохастикалық форма үшін

mах ср = —mіп(- ср)

болатынын көрсетуге болады. Олай болса, екіншісі үшін екі жақтылық есебі, 1-ші есеппен бірдей болады, яғни симметриялы пардағы екі есеп, өзара екі жақты деп аталады.

Екі жақтылық теоремасы.Егер бастапқы есептің тиімді шешімі бар болса, онда екінші екі жақтылық есебінің де тиімді шешімі бар болады. Әрі / -формасының минимумы, ср -формасының максимумына тең болады.

тіп / = тах ф.

Егер бастапқы есепте стохастикалық форма төменнен шектелмеген болса, онда екі жақтылық есебінде шектеулер жүйесінің бірде-бір теріс емес шешімі болмайды.

 Екі жақтылық әдісі

Екі жақтылық теоремасымен бағдарлау есебін шығарудың әдістерінің бірі «екі жақтылық әдісі»деп аталатын әдіспен байланысты. Оның мағынасы мынадай:

Стохастикалық бағдарлау есебі теңсіздік түріндегі шектеулермен берілсін. Оны біз Г-есебі ретінде қабылдайық та, ол үшін екі жақтылық есебі болып табылатын І-ші есепті құрып, оның шешімін симплекс-әдісті пайдаланып табайық. Бұл жағдайда екі жақтылық теоремасына мынадай қосымша енгізуге болады. Өзара екі жақты (I) жэне (II) есептері (1)-(2) жэне (1’)-(2′) түрлерінде берілсін. Егер І-ші есепте біз бастапқы базистен тиімді базиске көшетін болсақ, ал одан кейін (II) есебіндегі сәйкес базистерді ауыстырсақ, онда (II) есебіндегі жаңа базис те тиімді болады. Әрі біреуінің базисті шешімі екі жақтылық есебіндегі стохастикалық формасының сәйкес базисті емес айнымалыларының жанындағы коэффициенттерінің теріс таңбамен алынған жаңа базистік айнымалыларына теңестіру арқылы алынады.

Сәйкестік былай орнатылады:

Х        Х„       Х  …    Хп+т

.    .    .    .    .   .    .   .   .

Ут+1′»  Ут+п   У  . . .      Ут

Осылай тұжырымдау нәтижесінде (II) есебінің шешімі табылады. Екі    жақты    симплекс    эдіске    бастапқы    есептегі    шектеулер    саны, айнымалылар санынан көп болған кезде қолдану ыңғайлы.

Мысал 1. ф = -ух — 2у2 — Зу3формасын максимумдеу қажет болсын жэне ол келесі шектеулермен берілген

— у, + 2у7 -Зу. 1

2у -у2ъ<-1

Бұл қойылған есепті ІІ деп алдық. Оған екі жақтылық есебі болатын І-ші есепті құрайық. Ол былай қойылады:

f = x1 — x2формасын минимумдеу керек жэне ол төменгі шектеулермен жүйесімен берілген.

х + 2x2 >-1

2x1 — x2 > -2

— 3×1 — x2 > -3

Енді бұл форма мен шектеулердің қалай алынғанын түсіндірейік. Алдындағы келтірілген теория бойынша (1)-(2) жүйесі мен І»-формасы, (1′)-(2′) жүйесі және ср — формасымен мынадай байланыста болады.

ІІ есепте айнымалылар саны үшеу және де: сх = 1-ге, с2 = —1 —ге, т = 3 —тең. Теңдеудің саны екеу, ал жалпы түрде

ср = Ьхх2 2+… + Ът т.

Салыстыру арқылы bх =, b2 =-2, bъ=-3- терге тең және (Г)-тің бастапқы екі теңдеуінен

ап =-1,    а =2, аъх =-3, ап =2,      а22 =-1, аЪ2 =-.ә

-екендіктеріне көз жеткіземіз.

Осы айтылғандарды пайдаланып І-есепті құрамыз.

/ = схх22+… + сп -х„(2)

және (1)-ді табамыз. т = 3 — болғандықтан шектеулер жүйесі 3 теңдеуден тұруы керек.

сх = 1, с2 = -2 — болып табылды, сондықтан

/ = х+(-)-х2х2.

Енді (1)-ші шектеулерді жазайық

х2>-1 (=bХ) (2)-х1+(-)-х2>-2 (=b2) (-3)-х1+(-1)-х2>-3 (=b3

—    x1 + 2x2 >-1
2x1 — x2 > -2

—    3x1 — x2 > -3

F = x12.

—    Шектеулер жүйесі және мақсатты функция алынды. Екі жақтылы есебі шықты. Жаңадан    x3,x4,x5  айнымалыларын қосу арқылы І-ші есепті мына

—    түрде жазуға болады.

—    1-ші есеп.

—    — x1 + 2x2 — x3 = -1

—    5=-3

—    Минимумдеу

—    / = x1 — x2 + 0 • x3 + 0 • x4 + 0 • x

— және барлық

—    айнымалылары теріс болмауы тиіс.

—    Есепті  симплекс-әдіспен  шығаруға  болады.   Тек,   есепті  келесі  түрге келтірген ыңғайлы

x2 =2 + 2x4 — x5,

x2 = 2 + 2x1 — x4,

f= x 4    — x

    (II) — екінші  қарастырылған  есептің  экономикалық  мағынасына тоқталайық.

Ресурстары bг і = ,т) өндіріс орны, бұл ресурстарды өнім шығару үшін

пайдаланады және оны өткізеді. Мұндағы қойылатын мақсат ресурстарды артық жұмсамай, өнім бағасын максимум пайда түсіретіндей етіп анықтау керек. Бастапқы (1)-(2) есебі теріс еместік шартымен қарама-қарсыға өзгерген түрінде алайық.

Ах<b,         х>0,     /тах=тах<с,х>,    ху>0      (у=м)                               (3) берілсін.

мұндағы   x} — } -ші түрлі  өнім  шамасы,   С} — } -ші түрлі  өнімнің  бір

бірлігінің шамасы, ац — } -ші өнім бірлігіне жұмсалатын / -ші ресурс шығыны.

Енді қандай бір болмасын өндіріс орны ресурстарын өнім дайындауға шығындамай-ақ, тікелей сатуға жіберуге үйғарды дейік.

Ресурсты қандай бағамен сату керек деген сұрақ туады?

Баға сатушыны да, сатып алушыны да қанағаттандыруы керек.

Сатып алушы ресурстар үшін неғұрлым аз төлегенді қалайды. Ал сатушы ресурстар үшін алғаны, бұрынғы дайын тауар үшін тапқан пайдасынан кем болмауын қалайды.

Есептің математикалыц моделі

1. Мақсатты функция

     —             (р = bг у1+b2у2 +. . . + bтуп—>mіп  (4)

сатып алатын жақтың мақсатын сипаттайды. ут/-ші ресурс бірлігінің бағасы.

2.   Шектеулер жүйесі сатып алушылардың қызығушылығын сипаттайды.

Әр өнім бірлігін дайындауға кеткен ресурстарды бағалау қажет. Сонымен қатар ол ресурстар бағасын өткізілетін өнім бірліктері бағасымен шектеу керек

3. Баға теріс бола алмайды.

Уі>0   (і = ^п)                                                                                                      

Сонда екі жақтылық белгісі негізгі есеппен жаңадан қойылған есептердің модельдерін салыстыру арқылы алынады.

  Негізгі есептің үлгісі Жаңадан қойылған есептің үлгісі
1. f =>mах cр=>mіп
2. Шектеулер жүйесі < — типті Шектеулер жүйесі > — типті
3. С} -лер    мақсатты    функцияның коэффициенттері С}бос мүшелер
4. bгбос мүшелер bг-лер      мақсатты      функцияның коэффициенттері
5. Шектеулер                         жүйесінің коэффициенттері => бағаналар Шектеулер                            жүйесінің коэффициенттері => жатық жол

Осы екі жақты есеп болудың бес белгісі бізде бұрын айтылған.

І — тура есеп, ІІ — оған екі жақтылық есебі.

Егер осы есептердің біреуінің стохастикалық функциясы шектеусіз болса, онда есеп қарама-қайшылы.

Екі жақтылықтың екінші теоремасы:

Теорема 2 . Егер тура есептегі  стохастикалықфункция шектеусіз болса, онда оның екі жацтылыц есебінің шешімі болмайды.

Мысал 1.

Тура есеп.

f = 6x1 + 4x2-» mах

2x1 + 4x2 < 8

 2x1 + x2 < 6

Екі жақтылық есебі.

ср = 8x + 6x2 —> mіп

2y + 2y2 > 6

4у, +2у2 >4

х2>0.»  

 Екіжақтылықесебінқұрудыңмысалдары

Мысал 1. Стохастикалықпрограммалау есебі қойылған: cызықты функция

Ғ = 14x1 +10х2 +14х3 +11х

  келесі шектеулермен берілген.

4x1 + 2x2 + 2x4 + 3x4 < 35 x12 +2х3 +3х4 <30 3x1 + x2 + 2x3 + x4 > 40 x1 >0,х2 >0,х3 >0,х4 >0.

Қойылған есепке екі жақтылық есебін құрыңдар.

Шешуі:Есептің шешу жолын келесі кесте түрінде береміз.

Кесте 1

Алгоритм Қойылған алгоритмге сэйкестік
        Орнату

1

2     3

1

Шектеулер жүйесіндегі барлық Үшінші теңсіздікті  аталған түрге
  теңсіздіктерді < түріне келтіру

келтіреміз: —

3x1 — x2 — 2x3 — x4 < -40

2

Берілген жүйенің коффициенттері-      
  нен құрылған Аі матрицасын   4

2     2

3     35  
  құрамыз. Оның құрамына: Л 1

1     2

3    30  
  а) айнымалылар жанындағы коэф- А — -3

-1-2

-1-40  
  фициенттер енеді;   ^14 10 14 11     Ғ  
  б) шектеулер жүйесінің оң жағын- ч    
  дағы бос мүшелер бағанасы енеді;      
  в) стохастикалық функцияның қүрамын-      
  дағы айнымалылары жанындағы      
  коэффициенттер енеді.      

3

Аі матрицаның транспонирленген   4

і -:

3        14 л  
  А| матрицасын табу.   2

1   —

1        10  
    4 = 2

2 —

2     -14  
      3

3    —

1        11  
      V35

30 —

40       Ғу  

4

Алынған А| матрицасының негізінде 2 = 35у,

+ 30_у2

+ ? 40^3 -^ шіп
  екі жақтылық есебін құру. 4Уг+У2

~з.у3 ^

14,
    2Уі+У°2

«3^3   *

ю,
    х+2у2

-2Уъ

>14,
     

-Уз^

11
    У^0,у2

^0,у3

>0.

Мысал 2. Стохастикалық программалау есебі қойылған:

Ғ = 5x1 — x2 + 8x3 — x4 -^ тах

2x1 + 5x2 — x3 + 7x4 < 2,

■5x3 — x4 > 3,

■3x3 + 7x4 < 5

X

x2

x1 >0,х2 >0,х3 >0,х4 >0.

Қойылған есепке екі жақтылық есебін құрыңдар Екі жақтылық есебіне тағы бір мысал қарастырайық. Мысал -3. Берілген стохастикалық бағдарлау есебіне

+10х2 —» шіп

+3х2 > 6  + x2 > 9

I x — 8x < 8

x1 > 0, x2 > 0

екі жақтылық есебін құрыңдар.

  Шешуі: Келесі кестені құрайық.

Кесте 2

Алгоритм Қойылған алгоритмге сәйкестік орнату

1

2 3

1

Шектеулер жүйесіндегі барлық теңсіздіктерді > түріне келтіру Үшінші теңсіздікті  г+8x2 > -8 түріне келтіріп жазамыз

2

Кеңейтілген Аі матрицасын құрамыз. А = 13       6 3      1      9 -1   8   -8 ч80 10     2у  

3

Аі матрицаның транспонирленген 1′ матрицасын табу. 4 =      3    -1    80 3      1      8     10 6     9    -8     Ғ

V                                   )

 

4

Табылған       А|        матрицасының негізінде    екі    жақтылық    есебін құрыңдар. Ғ = вух + 9у2 — 8у3 —» тах ^+3^2-^3>80 Ъух2 +8у3 >10 Уі >0,у2 >0,у3 >0.