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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

141
Программа расчета оптимальных давлений промперегрева
Программа составлена на языке ПАСКАЛЬ и реализована на
ПЭВМ IBM PC AT. Текст программы приведен ниже.
Program GRADIENT;
Uses KPDmod;
Var A,A1,A2,Amax1,Amax2,GR1,GR2,E1,E2,dP1,dP2,KPD0,KPD1,
KPDmax,P0,PK,P10,P20,P11,PP11,P21,PP21,P1opt,P2opt : real;
IMF: string[20];
F : text;
Begin IMF:='GRAD.ISH';
Assign(F,IMF); reset(F);
Readln(F,P0); close (F);
Writeln('P0=',P0:4:1,'МПа',' PK=',PK:6:4);
Writeln('Введите нач. приближения Рпп1 и Рпп2 в МПа');
Readln(P10,P20);
E1:=0.0001; E2:=0.0001; {Погрешности определения Р1 и Р2}
dP1:=0.005*P10;
dP2:=0.005*P20;
Repeat
GR1:=GRAD1(P10,P20,dP1);
GR2:=GRAD2(P10,P20,dP2);
{ Выбираем пределы изменения шага A для P1 и P2 }
P11:=P10; P21:=P20;
Repeat
if GR1>0 then Amax1:=0.95*((P10-P21)/GR1);
if GR1<0 then Amax1:=0.95*(P10-P0)/GR1;
if GR2>0 then Amax2:=0.95*((P20-PK)/GR2);
if GR2<0 then Amax2:=0.95*(P20-P11)/GR2;
{В качестве максимального шага из Amax1 и Amax2 выбираем
наименьший }
If Amax1<Amax2 then A:=Amax1
else A:=Amax2;
A1:=A*GR1;
PP11:=P11;
P11:=P10-A1;
A2:=A*GR2;
PP21:=P21;
P21:=P20-A2;
Until (ABS((P11-PP11)/P11)<=0.001) and (ABS((P21-P21)/P21)<=0.01);
Repeat
A1:=A*GR1;
A2:=A*GR2;
P11:=P10-A1;
P21:=P20-A2;

Страницы