1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
| //Сортировка "Больше"
with MainForm.ActiveMDIChild as TChildForm do
begin
//строки в таблицах дочернего окна нумеруются с 1
//нулевая строка резервная
for i:=1 to SignsChild.RowCount-1 do
begin
if (SignsChild.Cells[0,i]='>') or (SignsChild.Cells[0,i]='>=') then
begin
inc(MoreCount);
SetLength(SimplexTable,LimChild.ColCount+2,MoreCount);
//коэффициенты
for j:=0 to LimChild.ColCount-1 do
SimplexTable[j+2,MoreCount-1]:=StrToFloat(LimChild.cells[j,i]);
//Пока нули (потом базис...)
SimplexTable[0,MoreCount-1]:=0;
//Значение (B i-ый)
SimplexTable[1,MoreCount-1]:=StrToFloat(BChild.cells[0,i]);
end;
end;
//Сортировка "Равно"
for i:=1 to SignsChild.RowCount-1 do
begin
if SignsChild.Cells[0,i]='=' then
begin
inc(EquallyCount);
SetLength(SimplexTable,LimChild.ColCount+2,MoreCount+EquallyCount);
//коэффициенты
for j:=0 to LimChild.ColCount-1 do
SimplexTable[j+2,MoreCount+EquallyCount-1]:=StrToFloat(LimChild.cells[j,i]);
//Пока нули (потом базис...)
SimplexTable[0,MoreCount+EquallyCount-1]:=0;
//Значение (B i-ый)
SimplexTable[1,MoreCount+EquallyCount-1]:=StrToFloat(BChild.cells[0,i]);
end;
end;
//Сортировка "Меньше"
for i:=1 to SignsChild.RowCount-1 do
begin
if (SignsChild.Cells[0,i]='<') or (SignsChild.Cells[0,i]='<=') then
begin
inc(LessCount);
SetLength(SimplexTable,LimChild.ColCount+2,MoreCount+EquallyCount+LessCount);
//коэффициенты
for j:=0 to LimChild.ColCount-1 do
SimplexTable[j+2,MoreCount+EquallyCount+LessCount-1]:=StrToFloat(LimChild.cells[j,i]);
//Пока нули (потом базис...)
SimplexTable[0,MoreCount+EquallyCount+LessCount-1]:=0;
//Значение (B i-ый)
SimplexTable[1,MoreCount+EquallyCount+LessCount-1]:=StrToFloat(BChild.cells[0,i]);
end;
end;
end;
{2}//---------------------------------------------------------------------------
//Порождение начального базиса
//2.1 Добавить коэф. -1 (>=)
for j:=0 to MoreCount-1 do
begin
Setlength(SimplexTable,length(SimplexTable)+1,MoreCount+EquallyCount+LessCount);
for i:=length(SimplexTable)-MoreCount+1 to length(SimplexTable)-1 do
SimplexTable[i,j]:=0;
SimplexTable[length(SimplexTable)-1,j]:=-1;
end;
end;
//Целевая функция GoalFun
GoalFun:=nil;
with MainForm.ActiveMDIChild as TChildForm do
begin
SetLength(GoalFun,GoalChild.ColCount+1);
for i:=1 to GoalChild.ColCount do
begin
if parametersForm.Min.Checked then GoalFun[i]:=StrToFloat(goalChild.Cells[i-1,1])
else GoalFun[i]:=-1*StrToFloat(goalChild.Cells[i-1,1]);
end;
end;
//Искусственная функция ArtFun
ArtFun:=nil;
SetLength(ArtFun,length(SimplexTable)-1-MoreCount);
//i=1 - Значение иск. функции
for i:=1 to length(SimplexTable)-3 do
for j:=0 to MoreCount-1 do ArtFun[i-1]:=ArtFun[i-1]-SimplexTable[i,j];
//------------------------------------------------------------------------------
//Минимизация искусственной функции
//Базис
if MoreCount>0 then
begin
for j:=0 to MoreCount-1 do
SimplexTable[0,j]:=length(simplexTable)-MoreCount+j-1;
for i:=MoreCount to length(simplexTable[0])-1 do
SimplexTable[0,i]:=length(simplexTable)-(LessCount+EquallyCount+MoreCount)+(i-MoreCount)-1;
end
else
for i:=0 to LessCount+EquallyCount-1 do
SimplexTable[0,i]:=length(simplexTable)-(LessCount+EquallyCount+MoreCount)+i-1;
//2 нижние строки для оценок
SetLength(SimplexTable,length(SimplexTable),length(SimplexTable[0])+2);
for i:=0 to length(GoalFun)-1 do SimplexTable[i+1,length(SimplexTable[0])-2]:=goalFun[i];
for i:=0 to length(ArtFun)-1 do SimplexTable[i+1,length(SimplexTable[0])-1]:=ArtFun[i];
SimplexTableNew:=nil;
SetLength(SimplexTableNew,length(SimplexTable),length(SimplexTable[0]));
//итерации...
up:
repeat
if not art then inc(IterCount);
if IterCount=Parametersform.CountIteration.Value then
begin
with MainForm.ActiveMDIChild as TChildForm do
begin
task.Items.Add('');
task.Items.Add('Достигнуто предельное число итераций. Решение не найдено');
exit;
end;
end;
{ |
(
