Курсовой проект - Динамическое программирование. Динамическая задача выбора объема партий. Вариант 8.2

Проверим правильность графика, подставив

i

y26

в целевую функцию.







N

i

N

n

ii

yxNxf

11

i

26)416()y ,(

0,

8,15,10

1122132





ii

yx

yxyxyxy

1.

10

32

xy

3

x

10

560*26)10*416()y ,(

i



i

xf

2

y

0

3

x

0

12,1610*26)0*40()y ,(

i



i

xf

2

y

10

2.

15

221

 yxy

3

x

10

1320*26))1510(42*16()y ,(

i



i

xf

2

y

0

2

x

15

1

y

0

3

x

10

75,7515*26)10*416()y ,(

i



i

xf

2

y

0

2

x

0

1

y

15

3

x

0

11

12,13210*26)25*416()y ,(

i



i

xf

2

y

10

2

x

25

1

y

0

3

x

0

61,4110*2625*26)0*40()y ,(

i



i

xf

2

y

10

2

x

0

1

y

25

3.

8

11

 yx

3

x

10

1800*26))81510(43*16()y ,(

i



i

xf

2

y

0

2

x

15

1

y

0

1

x

8

3

x

10

75,18315*26))1023(42*16()y ,(

i



i

xf

2

y

0

2

x

0

1

y

15

1

x

23

3

x

0

12,18010*26))825(42*16()y ,(

i



i

xf

2

y

10

2

x

25

1

y

0

1

x

8

3

x

0

61,189)10*2625*2()33*416()y ,(

i



i

xf

2

y

10

2

x

0

1

y

25

1

x

33

Рассчитав 3 шага, можно проверить программу, убрав ее ограничения на

количество отрезков N и поля «Затраты на хранения». Подставим значения в

программу:

12

Получаем, что по оптимальной программе хранить продукцию на складе не

выгодно и целевая функция не зависит от затрат на хранение S, следовательно,

график построен верно.

5. Алгоритм решения задачи

13

1. Проверка правильности ввода данных.

2. Начиная с последнего дня, подсчет на основе системы ограничений

всех возможных оптимальных для каждого k-го шага

комбинаций удовлетворения спроса.

3. Подсчет значений целевой функции для каждой

полученной комбинации удовлетворения спроса.

4. Поиск минимального значения целевой функции

и соответствующей ей

комбинации удовлетворения спроса.

5. Вывод в окно программы минимального значения целевой функции и

соответствующей ей комбинации удовлетворения спроса в виде удобном для

просмотра.

14

6. Листинг программы, реализующий алгоритм задачи

unit Unit1;

interface

uses

Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,

Dialogs, StdCtrls, ComCtrls, Grids, Menus, Math;

type

Tmas = array of array of byte;

TRec=record

days: integer;

day1: string[5];

day2: string[5];

day3: string[5];

day4: string[5];

day5: string[5];

day6: string[5];

day7: string[5];

day8: string[5];

day9: string[5];

end;

TDays=array[1..9] of integer;

var

Form1: TForm1;

F: File of Trec;

implementation

{$R *.dfm}

function GeneralFunction(XY: Tmas; Variation: integer):real;

var

SummXi, SummYi: real;

i,j,N: integer;

Ykoef: Tdays;

begin

Ykoef[1]:=StrToInt(Form1.Edit11.Text);

Ykoef[2]:=StrToInt(Form1.Edit12.Text);

Ykoef[3]:=StrToInt(Form1.Edit13.Text);

Ykoef[4]:=StrToInt(Form1.Edit14.Text);

Ykoef[5]:=StrToInt(Form1.Edit15.Text);

Ykoef[6]:=StrToInt(Form1.Edit16.Text);

Ykoef[7]:=StrToInt(Form1.Edit17.Text);

Ykoef[8]:=StrToInt(Form1.Edit18.Text);

Ykoef[9]:=StrToInt(Form1.Edit19.Text);

SummXi:=0;

SummYi:=0;

j:=Form1.TrackBar1.Position;

N:=0;

for i:=0 to 2*Form1.TrackBar1.Position-1 do

begin

if not odd(i) then

begin

SummXi:=SummXi+sqrt(XY[Variation][i]);

if XY[Variation][i]<>0 then

15

inc(N);

end

else

begin

SummYi:=SummYi+sqrt(XY[Variation][i] * Ykoef[j]);

dec(j);

end;

GeneralFunction:=RoundTo((StrToInt(Form1.Edit20.Text)*N +

StrToInt(Form1.Edit10.Text) * SummXi) + SummYi, -2);

end;

function mudoxren(range:integer; s:integer; w:integer; d:integer; var mas:Tmas;

err:boolean):tmas;

var i:integer;

begin

if err then

for i:=s to s+range-1 do

mas[i][w]:=d-mas[i][w+1]+mas[i][w-1]

else

for i:=s to s+range-1 do

mas[i][w+1]:=d+mas[i][w-1]+mas[i][w];

end;

procedure TForm1.Button1Click(Sender: TObject);

var

d: Tdays;

matr: Tmas;

i,j,current_day,range,s,w,u,max_i,i1,i2: integer;

err:boolean;

res: array of real;

temp: real;

begin

d[1]:=strtoint(Edit1.Text);

d[2]:=strtoint(Edit2.Text);

d[3]:=strtoint(Edit3.Text);

d[4]:=strtoint(Edit4.Text);

d[5]:=strtoint(Edit5.Text);

d[6]:=strtoint(Edit6.Text);

d[7]:=strtoint(Edit7.Text);

d[8]:=strtoint(Edit8.Text);

d[9]:=strtoint(Edit9.Text);

SetLength(matr, Round(Power(2,TrackBar1.Position-1)), TrackBar1.Position*2-1);

for i:=0 to Round(Power(2,TrackBar1.Position-1))-1 do

for j:=0 to TrackBar1.Position*2-1 do

matr[i,j]:=0;

range:=Round(Power(2, TrackBar1.Position-1)/2);

current_day:=TrackBar1.Position;

w:=0;

err:=true;

u:=1;

repeat

s:=0;

j:=1;

for i:=0 to u do

begin

mudoxren(range,s,w,d[current_day],matr,err);

err:=not(err);

s:=j*range;

16

inc(j);

end;

range:=Round(range/2);

current_day:=current_day-1;

w:=w+2;

u:=(u+1)*2-1;

until u>Round(Power(2, TrackBar1.Position)-1)-1;

for i:=0 to Round(Power(2, TrackBar1.Position-1))-1 do

matr[i][TrackBar1.Position*2-2]:=matr[i][TrackBar1.Position*2-3]+d[1];

SetLength(res,Round(Power(2, TrackBar1.Position-1)));

for i:=0 to Round(Power(2, TrackBar1.Position-1))-1 do

res[i]:=GeneralFunction(matr,i);

max_i:=0;

temp:=res[0];

for i:=1 to Round(Power(2, TrackBar1.Position-1))-1 do

if temp > res[i] then

begin

temp:=res[i];

max_i:=i;

end;

i1:=TrackBar1.Position;

i2:=TrackBar1.Position;

for j:=0 to 2*TrackBar1.Position-2 do

begin

if not odd(j) then

begin

if matr[max_i][j]<>0 then

StringGrid1.Cells[2,i1]:=' '+FloatToStr(matr[max_i][j])

else

StringGrid1.Cells[2,i1]:=' - ';

StringGrid1.Cells[0,i1]:=' '+IntToStr(i1);

StringGrid1.Cells[1,i1]:=' '+IntToStr(d[i1]);

dec(i1);

end

else

begin

if matr[max_i][j]<>0 then

begin

StringGrid1.Cells[4,i2-1]:=' '+FloatToStr(matr[max_i][j]);

StringGrid1.Cells[3,i2]:=' '+FloatToStr(matr[max_i][j]);

end

else

begin

StringGrid1.Cells[4,i2-1]:=' - ';

StringGrid1.Cells[3,i2]:=' - ';

end;

dec(i2);

end;

StringGrid1.Cells[4,TrackBar1.Position]:=' - ';

StringGrid1.Cells[3,1]:=' - ';

Label18.Caption:='Общие затраты: '+FloatToStr(temp)+' у.е.';

end;

end.

procedure TForm1.Button2Click(Sender: TObject);

var

d: Tdays;

matr: Tmas;

17

i,j,current_day,range,s,w,u,Ykoef: integer;

err:boolean;

res: array of real;

temp,temp2,min: real;

p: boolean;

begin

form2.Show;

Form2.Chart1.Series[0].Clear;

Form2.Chart1.Series[1].Clear;

Form2.Chart1.Series[2].Clear;

d[1]:=strtoint(Edit1.Text);

d[2]:=strtoint(Edit2.Text);

d[3]:=strtoint(Edit3.Text);

d[4]:=strtoint(Edit4.Text);

d[5]:=strtoint(Edit5.Text);

d[6]:=strtoint(Edit6.Text);

d[7]:=strtoint(Edit7.Text);

d[8]:=strtoint(Edit8.Text);

d[9]:=strtoint(Edit9.Text);

Ykoef:=1;

temp:=0;

p:=true;

SetLength(matr, Round(Power(2,TrackBar1.Position-1)), TrackBar1.Position*2-1);

repeat

temp2:=temp;

for i:=0 to Round(Power(2,TrackBar1.Position-1))-1 do

for j:=0 to TrackBar1.Position*2-1 do

matr[i,j]:=0;

range:=Round(Power(2, TrackBar1.Position-1)/2);

current_day:=TrackBar1.Position;

w:=0;

err:=true;

u:=1;

repeat

s:=0;

j:=1;

for i:=0 to u do

begin

mudoxren(range,s,w,d[current_day],matr,err);

err:=not(err);

s:=j*range;

inc(j);

end;

range:=Round(range/2);

current_day:=current_day-1;

w:=w+2;

u:=(u+1)*2-1;

until u>Round(Power(2, TrackBar1.Position-1))-1;

for i:=0 to Round(Power(2, TrackBar1.Position-1))-1 do

matr[i][TrackBar1.Position*2-2]:=matr[i][TrackBar1.Position*2-3]+d[1];

SetLength(res,Round(Power(2, TrackBar1.Position-1)));

for i:=0 to Round(Power(2, TrackBar1.Position-1))-1 do

res[i]:=GeneralFunction(matr,Ykoef,i);

temp:=res[0];

for i:=1 to Round(Power(2, TrackBar1.Position-1))-1 do

if temp > res[i] then

temp:=res[i];

if p then

min:=temp;

p:=false;

Form2.Chart1.Series[0].AddXY(Ykoef,temp);

inc(Ykoef);

until temp=temp2;

Form2.Chart1.Series[0].AddXY(Ykoef+8,temp);

Form2.Chart1.Series[2].AddXY(Ykoef,temp+3);

Form2.Chart1.Series[1].AddXY(0,temp);

18

Form2.Chart1.Series[1].AddXY(Ykoef-2,temp);

Form2.Chart1.Series[1].AddXY(Ykoef-2,min);

end;

7. Руководство пользователя

7.1 Системные требования

Процессор: Pentium I или аналогичный AMD 400 MHz и выше

ОЗУ: 64 Мб и более

ОС: Windows 98, 2000, ХР

7.2 Описание возможностей

Данная программа предназначена для расчета календарной программы

выпуска некоторого вида изделий на плановый период, состоящий из N отрезков в

течение 7,8 или 9 последовательных дней. Выходными данными для программы

являются: минимальная стоимость затрат на выпуск и изготовление деталей.

Данные, для которых разрабатывалась данная программа, введены по

умолчанию. Поэтому для расчета оптимальной программы выпуска продукции на

основе этих данных достаточно нажать на кнопку "Рассчитать".

19

7.3 Основное окно программы

Главное меню

программы

Область начальных данных

Область вывода результата Кнопка расчета результата

Кнопка вывода графика зависимости

целевой функции от затрат на хранение

Кнопка выхода из программы

20

Курсовой проект - Динамическое программирование. Динамическая задача выбора объема партий. Вариант 8.2

Подождите немного. Документ загружается.