Смекни!
smekni.com

Решение задач нелинейного программирования (стр. 2 из 3)

Второй этап – составление математической модели задачи

1. На основании полученных в первом этапе исходных данных и описания заданного производственного процесса составляется следующая таблица:

Продукты ресурсы 1 2 3
I 8 4 6 1000
II 160 240 200 280000

Через Х1 обозначим ресурсы I вида.

Через Х2 обозначим ресурсы II вида.

2. Обращаясь к условиям задачи, определяем все возможные ограничения, объединяя их в систему ограничений.

1 + 4Х2 + 6Х3 ≤ 1000

240Х1+ 200Х2 + 160Х3 ≤ 280000

Таким образом, получили задачу нелинейного программирования. Такие задачи называются задачами нелинейного программирования.

Решение задач нелинейного программирования осуществляется приведением их к задачам линейного программирования.

Для решения задачи линейного программирования применяется симплекс – метод.

Третий этап – выбор метода решения полученной математической задачи

Решение

1. Для решения задач линейного программирования симплекс – методом задача приводиться к каноническому виду:


1 + 4Х2 + 6Х3 + Х4= 1000

240Х1+ 200Х2 + 160Х3 + Х5= 280000

2. Составляем таблицу и определяем все неотрицательные базисные решения системы.

Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х4 8 4 6 1 0 1000
Х5 240 200 160 0 1 280000

А) Нашли некоторое неотрицательное базисное решение: Х4 =1300, Х5 = 190000. По заданию продолжаем искать базисные решения. Разрешающим элементом выбираем в 1 строке – Х2. Соответственно вся строка делится на 8, а все остальные элементы находятся по правилу прямоугольника.

Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х4 8 4 6 1 0 1000
Х5 240 200 160 0 1 280000
Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х2 ¾ 1 ½ 1/8 0 325/2
Х5 90 0 60 -25 1 157500

Б) Нашли некоторое неотрицательное базисное решение: Х2 =325/2, Х5 =157500. По заданию продолжаем искать базисные решения. Разрешающим элементом выбираем в 1 строке – Х1. Соответственно вся строка делится на 3/4, а все остальные элементы находятся по правилу прямоугольника.

Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х2 ¾ 1 ½ 1/8 0 325/2
Х5 90 0 60 -25 1 157500
Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х1 1 4/3 2/3 1/6 0 650/3
Х5 0 -120 0 -40 1 138000

В) Нашли некоторое неотрицательное базисное решение: Х1 =650/3, Х5 =138000. По заданию продолжаем искать базисные решения. Разрешающим элементом выбираем в 1 строке – Х3. Соответственно вся строка делится на 2/3, а все остальные элементы находятся по правилу прямоугольника.

Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х1 1 4/3 2/3 1/6 0 650/3
Х5 0 -120 0 -40 1 138000
Базисные переменные Х1 Х2 Х3 Х4 Х5 Свободный член
Х3 3/2 2 1 1/4 0 325
Х5 0 -120 0 -40 1 138000

Г) Нашли некоторое неотрицательное базисное решение: Х5 =138000, Х3 =325. Найдены все неотрицательные базисные решения.

2. Находим получаемый продукт.

Х1= 6*0+8*0+4*0=0

Х2=240*0+200*0+160*0=0

У1=3*00,4*00,6=0

У2=5*00,5*00,5=0

У3=8*00,6*00,4=0

F1=0+0+0=0

Х1= 6*0+8*325/2+4*0=1300

Х2=240*0+200*325/2+160*0=32500

У1=3*13000,4*325000,6=26904,728

У2=5*13000,5*325000,5=32500

У3=8*13000,6*325000,4=37688,542

F2=26904,728 +32500 +37688,542 = 97093,27

Х1= 6*650/3+8*0+4*0=1300

Х2=240*650/3+200*0+160*0=52000

У1=3*13000,4*520000,6=35699,794

У2=5*13000,5*520000,5=41109,610

У3=8*13000,6*520000,4=45483,862

F3= 35699,794+ 41109,610+ 45483,862= 122263,266

Х1= 6*0+8*0+4*325=1300

Х2=240*0+200*0+160*325=52000

У1=3*13000,4*520000,6=35699,794

У2=5*13000,5*520000,5=41109,610

У3=8*13000,6*520000,4=45483,862

F3= 35699,794+ 41109,610+ 45483,862= 122263,266

F1 < F2

F2 < F3

F3 = F4

Ответ: Fmax= 122263,266

Четвертый этап – подготовка словесного алгоритма решения задачи

1. Вводим данные в таблицу

2. Выбираем разрешающий элемент:

2.1. Берем каждый неотрицательный элемент первой строки и делим на свободный член первой строки.

2.2. Находим среди всех деленных элементов минимальный.

2.3. Берем каждый неотрицательный элемент второй строки и делим на свободный член второй строки.

2.4. Находим среди всех деленных элементов минимальный.

2.5. Берем каждый неотрицательный элемент n-ой строки и делим на свободный член n-ой строки.

2.6. Находим среди всех деленных элементов минимальный.

2.7. Берем минимальные элементы первой, второй и n-ой строки и среди них находим минимальный (это и будет разрешающий элемент). При условии если минимальные элементы строк совпадают, берется элемент первой строки.

3. Вычисляем всю таблицу методом прямоугольника относительно разрешающего элемента:

3.1. Умножаем разрешающий элемент на элемент решаемой строки.

3.2. Отнимаем произведение соответствующего элемента решаемой строки на элемент разрешающего столбца решаемой строки

3.3. И делим ответ на разрешающий элемент.

3.4. Делим разрешающую строку на разрешающий элемент.

3.5. Берем каждый элемент разрешающей строки и делим на разрешающий элемент.

3.6. Всем элементам, кроме разрешающего элемента, разрешающего столбца присвоим (0)

3.7. Разрешающему элементу присвоим (1).

В индексе разрешающей строки присвоить индекс

4. Повторяем процедуру вычисления с 2 пункта.

5. В конечном результате находим все неотрицательные базисные решения. Подставляем значения и находим получаемый продукт.

6. Находим все F.

7. Выбираем наибольшую из них, которая будет являться оптимальным планом выпуска продукции.

Пятый этап – разработка программы для решения задачи

Private Sub Form_Load()

Left = (Screen. Width – Width) &bsol; 2

Top = (Screen. Height – Height) &bsol; 2

End Sub

Private Sub Timer1_Timer()

Unload Form1

Load Form2

Form2. Show

End Sub

‘Объявление переменных

Public a As Integer

Public b As Integer

Public c As Integer

Public d As Integer

Public e As Integer

Public f As Integer

Public aa As Integer

Public ab As Integer

Public ac As Integer

Public ad As Integer

Public ae As Integer

Public af As Integer

Public ba As Integer

Public bb As Integer

Public bc As Integer

Public bd As Integer

Public be As Integer

be = Text17. Text

bf = Text18. Text

ca = Text19. Text

cb = Text20. Text

cc = Text21. Text

cd = Text22. Text

ce = Text23. Text

cf = Text24. Text

X1 = Text25. Text

X2 = Text26. Text

X3 = Text27. Text

‘Проверка выполнения равенств

If a*x1+aa*x2+ba*x3=ca Then «Равенство выполняется» Else «Равенство не выполняется»

If b*x1+ab*x2+bb*x3=cb Then «Равенство выполняется» Else «Равенство не выполняется»

If c*x1+ac*x2+bc*x3=cc Then «Равенство выполняется» Else «Равенство не выполняется»

If d*x1+ad*x2+bd*x3=cd Then «Равенство выполняется» Else «Равенство не выполняется»

If e*x1+ae*x2+be*x3=ce Then «Равенство выполняется» Else «Равенство не выполняется»

F= f*x1+af*x2+bf*x3

If F<fmin Then «Решение не выполняется» Else «Решение выполняется, план является оптимальным»

Text28. Visible = True

Text29. Visible = True

End Sub

Private Sub Command2_Click()

‘очистка текстовых окон для следующего ввода данных

Text1. Text = «»

Text2. Text = «»

Text3. Text = «»

Text4. Text = «»

Text5. Text = «»

Text6. Text = «»

Text7. Text = «»

Text8. Text = «»

Text9. Text = «»

Text10. Text = «»

Text11. Text = «»

Text12. Text = «»

Text13. Text = «»

Text14. Text = «»

Text15. Text = «»

Text16. Text = «»

Text17. Text = «»

Text18. Text = «»

Text19. Text = «»

Text20. Text = «»

Text21. Text = «»

Text22. Text = «»

Text23. Text = «»

Text24. Text = «»

Text25. Text = «»

Text26. Text = «»

Text27. Text = «»

Text28. Visible = False

Text29. Visible = False

End Sub

Private Sub Command3_Click()

‘показать справку

Unload Form2

Load Form3

Form3. Show

End Sub

Private Sub Command4_Click()

Unload Form2

End Sub

Private Sub Form_Load()

Left = (Screen. Width – Width) &bsol; 2

Top = (Screen. Height – Height) &bsol; 2

‘подготовка текстовых окон к вводу данных при запуске рабочего окна

Text1. Text = «»

Text2. Text = «»

Text3. Text = «»

Text4. Text = «»

Text5. Text = «»

Text6. Text = «»

Text7. Text = «»

Text8. Text = «»

Text9. Text = «»

Text10. Text = «»

Text11. Text = «»

Text12. Text = «»

Text13. Text = «»

Text14. Text = «»

Text15. Text = «»

Text16. Text = «»

Text17. Text = «»

Text18. Text = «»

Text19. Text = «»

Text20. Text = «»

Text21. Text = «»

Text22. Text = «»

Text23. Text = «»

Text24. Text = «»

Text25. Text = «»

Text26. Text = «»