Методы оптимизации и расчеты на ЭВМ технико-экономических задач. Ромашова О.Ю. - 144 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

144
Для расчета целевой функции воспользуемся готовой программой
расчета термического КПД. Для этого из приложения 1 выберем модуль
KPDMOD3
, содержащий функцию KPD3(P0, tpw, Pk, X0). Первые
два параметра функции – оптимизируемые, остальные
(Pk и X0)
должны быть заданы в основной программе в качестве исходных данных.
4. Моделирование оптимизационной задачи
Алгоритм расчета оптимальных значений P
0
и t
п.в.
методом координатного спуска
Алгоритм расчета методом координатного спуска приведен на
рис. 4.30. Для одномерной минимизации используется метод «золотого»
сечения.
Программа расчета максимального значения КПД
методом координатного спуска
Программа составлена в соответствии с алгоритмом на рис. 4.30.
Ввод исходных данных для простоты приведен в тексте основной про-
граммы.
Program KOORDIN;
Uses TABL1, KPDMOD3; {Подключение модуля для расчета КПД}
Var x0, Pk, … : real;
Begin { Начало основной программы }
PK:=0.005; {МПа} X0:=0.995;
E1:=
P0max:=22; {МПа}
tks:=TSP(PK); tpwmin:=tks;
P0opt:=0; tpwopt:=0;
tpwf:=100; { Зафиксировали второй параметр }
Writeln(‘ tpwf= ‘,tpwf:6:1);
REPEAT { Координатный спуск }
{ Метод «золотого» сечения для оптимизации P0 }
P0min:= PST(tpwf);
a:=P0min;
b:=P0max;
Writeln(‘ Границы изменения tpw ‘);
Writeln(‘ a= ‘, a:6:3, ‘ b= ’, b:6:3);
REPEAT { «Золотое сечение» }
P01:=a+0.3820*(b-a);
P02:=a+0.6180*(b-a);
f1:=KPD3(P01,tpwf,Pk,x0);
f2:= KPD3(P02,tpwf,Pk,x0);
If f1<=f2 then b:=x2;

Страницы