Добавил:
Studfiles2
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Программы к курсу методы оптимизации / LAB3 / SHELP
.TXT ‹ Ў®а в®а п а Ў®в N 3
‘ € Њ Џ ‹ … Љ ‘ Ќ › ‰
Њ … ’ Ћ „
ђ бб¬ ваЁў Ґвбп б«Ґ¤гой п § ¤ з «ЁҐ©®Ј® Їа®Ја ¬¬Ёа®ў Ёп .
Ќ ©вЁ ¬ЁЁ¬г¬ «ЁҐ©®© дгЄжЁЁ f(x1,x2,...,xn):
f = c[1]*x[1] + c[2]*x[2] +...+ c[n]*x[n] ,
Ј¤Ґ c[i] - Ї®бв®плҐ Є®нддЁжЁҐвл ,
¬®¦Ґб⢥ , § ¤ ®¬ Ў®а®¬ «ЁҐ©ле ®Ја ЁзҐЁ© :
a[1,1]*x[1] + ... + a[1,n]*x[n] >= b[1]
...
a[m,1]*x[1] + ... + a[m,n]*x[n] >= b[m]
x[1]>=0,...,x[n]>=0 ,
Ј¤Ґ a[i,j],b[i] - Ї®бв®плҐ Є®нддЁжЁҐвл .
‚ ¬ ваЁз®© д®а¬Ґ ®Ја ЁзҐЁп § ЇЁблў овбп б«Ґ¤гойЁ¬ ®Ўа §®¬ :
AX>=B , X>=0 .
–Ґ«Ґў п дгЄжЁп ¬®¦Ґв Ўлвм ЇаҐ¤бв ў«Ґ ў ўЁ¤Ґ бЄ «па®Ј®
Їа®Ё§ўҐ¤ҐЁп :
f = ( C,X ) .
ЉђЂ’Љ€… ЋЃ™€… ‘‚…„…Ќ€џ
------------------------
‘Ё¬Ї«ҐЄбл© ¬Ґв®¤ аҐиҐЁп § ¤ зЁ «ЁҐ©®Ј® Їа®Ја ¬¬Ёа®ў Ёп
б®бв®Ёв Ё§ ¤ўге нв Ї®ў :
1) Ї®ЁбЄ Єа ©Ґ© в®зЄЁ ¤®ЇгбвЁ¬®Ј® ¬®¦Ґбвў ,
2) Ї®ЁбЄ ®ЇвЁ¬ «м®© в®зЄЁ Їг⥬ Їа ў«Ґ®Ј® ЇҐаҐЎ®а
Єа ©Ёе в®зҐЄ .
Ља ©пп в®зЄ Ґ бгйҐбвўгҐв , Ґб«Ё ў в Ў«ЁжҐ бгйҐбвўгҐв бва®Є ,
-----------------
ўбҐ н«Ґ¬Ґвл Є®в®а®© ҐЇ®«®¦ЁвҐ«мл , Ї®б«Ґ¤Ё© н«Ґ¬Ґв -
®ваЁж ⥫мл© .
Ља ©пп в®зЄ © ¤ Ґ , Ґ«Ё ўбҐ н«Ґ¬Ґвл ўҐЄв®а -бв®«Ўж B
---------------
Ў®«миҐ г«п .
—в®Ўл ©вЁ Єа ©оо в®зЄг , ¤® :
1) ўлЎа вм бва®Єг i , ў Є®в®а®© b[i] < 0;
2) ўлЎа вм бв®«ЎҐж s , ў Є®в®а®¬ a[i,s]>=0;
3) ў бв®«ЎжҐ s § ¤ вм ®¬Ґа бва®ЄЁ r а §аҐи о饣® н«Ґ¬Ґв в Є ,
зв®Ўл ®ваЁж ⥫쮥 ®в®иҐЁҐ b[r]/a[r,s] Ўл«® ¬ ЄбЁ¬ «мл¬ .
4) Ї®¬Ґпвм ¬Ґбв ¬Ё Ё¬Ґ Є®®а¤Ё в ў в Ў«ЁжҐ Ё§ бва®ЄЁ r Ё бв®«Ўж s;
5) а бб¬ ваЁў п н«Ґ¬Ґв a[r,s] Є Є а §аҐи ойЁ© , Ґ®Ўе®¤Ё¬®
ЇаҐ®Ўа §®ў вм в Ў«Ёжг Ї® д®а¬г« ¬ :
ARS:= a[r,s];
z1[r,s]:= 1/ARS;
z1[r,j]:= -z[r,j]/ARS , j<>s;
z1[i,s]:= z[i,s]/ARS , i<>r;
z1[i,j]:= (z[i,j]*ARS - z[i,s]*z[r,j])/ARS , i<>r,j<>s;
z:=z1,
Ј¤Ґ Ї®¤ z Ё z1 Ї®Ё¬ Ґвбп ᮮ⢥вб⢥® ЇҐаў® з «м®Ґ Ё ЇаҐ®Ўа §®ў ®Ґ
§ 票Ґ в Ў«Ёжл ( Єа®¬Ґ «Ґў®Ј® бв®«Ўж Ё ўҐа奩 бва®ЄЁ ).
ЋЇвЁ¬ «м п в®зЄ © ¤ Ґ , Ґб«Ё ўбҐ н«Ґ¬Ґвл ўҐЄв®а-бва®ЄЁ ‘ >= 0
---------------
( ЇаЁ н⮬ ўбҐ н«Ґ¬Ґвл ўҐЄв®а-бв®«Ўж B >= 0 ).
ЋЇвЁ¬ «м п в®зЄ Ґ бгйҐбвўгҐв , Ґб«Ё ў в Ў«ЁжҐ Ґбвм бв®«ЎҐж j ,
-----------------
ў Є®в®а®¬ c[j] < 0 , ўбҐ a[i,j]>0 ЇаЁ «оЎ®¬ i .
—в®Ўл ©вЁ ®ЇвЁ¬ «мго в®зЄг , ¤® :
1) ўлЎа вм бв®«ЎҐж s , ў Є®в®а®¬ c[s] < 0;
2) ў бв®«ЎжҐ s § ¤ вм ®¬Ґа бва®ЄЁ r а §аҐи о饣® н«Ґ¬Ґв в Є ,
зв®Ўл ®ваЁж ⥫쮥 ®в®иҐЁҐ b[r]/a[r,s] Ўл«® ¬ ЄбЁ¬ «мл¬ ;
3) Ї®¬Ґпвм ¬Ґбв ¬Ё Ё¬Ґ Є®®а¤Ё в ў в Ў«ЁжҐ Ё§ бва®ЄЁ r Ё бв®«Ўж s;
4) а бб¬ ваЁў п н«Ґ¬Ґв a[r,s] Є Є а §аҐи ойЁ© , Ґ®Ўе®¤Ё¬®
ЇаҐ®Ўа §®ў вм в Ў«Ёжг Ї® д®а¬г« ¬ ( б¬.ўлиҐ ).
Љ®®а¤Ё вл ®ЇвЁ¬ «м®© в®зЄЁ ®ЇаҐ¤Ґ«повбп б«Ґ¤гойЁ¬ ®Ўа §®¬ :
1) Ґб«Ё x[j] 室Ёвбп i-¬ ¬Ґб⥠«Ґў®Ј® бв®«Ўж , в® ҐЈ® § 票Ґ
а ў® b[i];
2) Ґб«Ё x[i] 室Ёвбп j-¬ ¬Ґб⥠ўҐа奩 бва®ЄЁ , в® ҐЈ® § 票Ґ
а ў® 0 .
‘ € Њ Џ ‹ … Љ ‘ Ќ › ‰
Њ … ’ Ћ „
ђ бб¬ ваЁў Ґвбп б«Ґ¤гой п § ¤ з «ЁҐ©®Ј® Їа®Ја ¬¬Ёа®ў Ёп .
Ќ ©вЁ ¬ЁЁ¬г¬ «ЁҐ©®© дгЄжЁЁ f(x1,x2,...,xn):
f = c[1]*x[1] + c[2]*x[2] +...+ c[n]*x[n] ,
Ј¤Ґ c[i] - Ї®бв®плҐ Є®нддЁжЁҐвл ,
¬®¦Ґб⢥ , § ¤ ®¬ Ў®а®¬ «ЁҐ©ле ®Ја ЁзҐЁ© :
a[1,1]*x[1] + ... + a[1,n]*x[n] >= b[1]
...
a[m,1]*x[1] + ... + a[m,n]*x[n] >= b[m]
x[1]>=0,...,x[n]>=0 ,
Ј¤Ґ a[i,j],b[i] - Ї®бв®плҐ Є®нддЁжЁҐвл .
‚ ¬ ваЁз®© д®а¬Ґ ®Ја ЁзҐЁп § ЇЁблў овбп б«Ґ¤гойЁ¬ ®Ўа §®¬ :
AX>=B , X>=0 .
–Ґ«Ґў п дгЄжЁп ¬®¦Ґв Ўлвм ЇаҐ¤бв ў«Ґ ў ўЁ¤Ґ бЄ «па®Ј®
Їа®Ё§ўҐ¤ҐЁп :
f = ( C,X ) .
ЉђЂ’Љ€… ЋЃ™€… ‘‚…„…Ќ€џ
------------------------
‘Ё¬Ї«ҐЄбл© ¬Ґв®¤ аҐиҐЁп § ¤ зЁ «ЁҐ©®Ј® Їа®Ја ¬¬Ёа®ў Ёп
б®бв®Ёв Ё§ ¤ўге нв Ї®ў :
1) Ї®ЁбЄ Єа ©Ґ© в®зЄЁ ¤®ЇгбвЁ¬®Ј® ¬®¦Ґбвў ,
2) Ї®ЁбЄ ®ЇвЁ¬ «м®© в®зЄЁ Їг⥬ Їа ў«Ґ®Ј® ЇҐаҐЎ®а
Єа ©Ёе в®зҐЄ .
Ља ©пп в®зЄ Ґ бгйҐбвўгҐв , Ґб«Ё ў в Ў«ЁжҐ бгйҐбвўгҐв бва®Є ,
-----------------
ўбҐ н«Ґ¬Ґвл Є®в®а®© ҐЇ®«®¦ЁвҐ«мл , Ї®б«Ґ¤Ё© н«Ґ¬Ґв -
®ваЁж ⥫мл© .
Ља ©пп в®зЄ © ¤ Ґ , Ґ«Ё ўбҐ н«Ґ¬Ґвл ўҐЄв®а -бв®«Ўж B
---------------
Ў®«миҐ г«п .
—в®Ўл ©вЁ Єа ©оо в®зЄг , ¤® :
1) ўлЎа вм бва®Єг i , ў Є®в®а®© b[i] < 0;
2) ўлЎа вм бв®«ЎҐж s , ў Є®в®а®¬ a[i,s]>=0;
3) ў бв®«ЎжҐ s § ¤ вм ®¬Ґа бва®ЄЁ r а §аҐи о饣® н«Ґ¬Ґв в Є ,
зв®Ўл ®ваЁж ⥫쮥 ®в®иҐЁҐ b[r]/a[r,s] Ўл«® ¬ ЄбЁ¬ «мл¬ .
4) Ї®¬Ґпвм ¬Ґбв ¬Ё Ё¬Ґ Є®®а¤Ё в ў в Ў«ЁжҐ Ё§ бва®ЄЁ r Ё бв®«Ўж s;
5) а бб¬ ваЁў п н«Ґ¬Ґв a[r,s] Є Є а §аҐи ойЁ© , Ґ®Ўе®¤Ё¬®
ЇаҐ®Ўа §®ў вм в Ў«Ёжг Ї® д®а¬г« ¬ :
ARS:= a[r,s];
z1[r,s]:= 1/ARS;
z1[r,j]:= -z[r,j]/ARS , j<>s;
z1[i,s]:= z[i,s]/ARS , i<>r;
z1[i,j]:= (z[i,j]*ARS - z[i,s]*z[r,j])/ARS , i<>r,j<>s;
z:=z1,
Ј¤Ґ Ї®¤ z Ё z1 Ї®Ё¬ Ґвбп ᮮ⢥вб⢥® ЇҐаў® з «м®Ґ Ё ЇаҐ®Ўа §®ў ®Ґ
§ 票Ґ в Ў«Ёжл ( Єа®¬Ґ «Ґў®Ј® бв®«Ўж Ё ўҐа奩 бва®ЄЁ ).
ЋЇвЁ¬ «м п в®зЄ © ¤ Ґ , Ґб«Ё ўбҐ н«Ґ¬Ґвл ўҐЄв®а-бва®ЄЁ ‘ >= 0
---------------
( ЇаЁ н⮬ ўбҐ н«Ґ¬Ґвл ўҐЄв®а-бв®«Ўж B >= 0 ).
ЋЇвЁ¬ «м п в®зЄ Ґ бгйҐбвўгҐв , Ґб«Ё ў в Ў«ЁжҐ Ґбвм бв®«ЎҐж j ,
-----------------
ў Є®в®а®¬ c[j] < 0 , ўбҐ a[i,j]>0 ЇаЁ «оЎ®¬ i .
—в®Ўл ©вЁ ®ЇвЁ¬ «мго в®зЄг , ¤® :
1) ўлЎа вм бв®«ЎҐж s , ў Є®в®а®¬ c[s] < 0;
2) ў бв®«ЎжҐ s § ¤ вм ®¬Ґа бва®ЄЁ r а §аҐи о饣® н«Ґ¬Ґв в Є ,
зв®Ўл ®ваЁж ⥫쮥 ®в®иҐЁҐ b[r]/a[r,s] Ўл«® ¬ ЄбЁ¬ «мл¬ ;
3) Ї®¬Ґпвм ¬Ґбв ¬Ё Ё¬Ґ Є®®а¤Ё в ў в Ў«ЁжҐ Ё§ бва®ЄЁ r Ё бв®«Ўж s;
4) а бб¬ ваЁў п н«Ґ¬Ґв a[r,s] Є Є а §аҐи ойЁ© , Ґ®Ўе®¤Ё¬®
ЇаҐ®Ўа §®ў вм в Ў«Ёжг Ї® д®а¬г« ¬ ( б¬.ўлиҐ ).
Љ®®а¤Ё вл ®ЇвЁ¬ «м®© в®зЄЁ ®ЇаҐ¤Ґ«повбп б«Ґ¤гойЁ¬ ®Ўа §®¬ :
1) Ґб«Ё x[j] 室Ёвбп i-¬ ¬Ґб⥠«Ґў®Ј® бв®«Ўж , в® ҐЈ® § 票Ґ
а ў® b[i];
2) Ґб«Ё x[i] 室Ёвбп j-¬ ¬Ґб⥠ўҐа奩 бва®ЄЁ , в® ҐЈ® § 票Ґ
а ў® 0 .