Оптимизация химико-технологических процессов органического синтеза. Бочкарев В.В - 48 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

48
{$I C:\pascal\Procedur\noLineUr.pas}
{***********************************************************}
Procedure GoldMin (A, B, E:real; var xGold,FGold:real);
{"Метод золотого сечения"
ПроцедураGoldMin производит поиск интервала[A,B],
в котором лежит решение уравненияF(x)=0}
var interval,t1,t2,x0,x1,x2,x3,F0,F1,F2,F3 : real;
{***********************************************************}
Procedure Fx (var FGold:real; var xGold:real);
begin
alfa:=xGold;
noLineUraw(n,e,intX,f,x);
Ratio:=xh2/xco;
writeln(alfa:5:3,xco:7:3,xh2:7:3,xch4:7:3,xco2:7:3,xh2o:7:3,Ratio:7:3);
FGold:=sqr(Ratio-RatioOpt);
End;
{***********************************************************}
begin
t1:=0.3819660113; t2:=1 - t1;
x0:=A; x1:=A + t1*(B - A); x2:=A + t2*(B - A); x3:=B;
Fx(F1, x1); Fx(F2, x2);
interval:=x3 - x1;
While interval>E do begin
If F2 < F1 then
begin
interval:=x3 - x1; x0:=x1; x1:=x2; x2:=x0 +
t2*interval;
F1:=F2; Fx(F2,x2);
end {end then}
else
begin
interval:=x2-x0; x3:=x2; x2:=x1; x1:=x0+t1*interval;
F2:=F1; Fx(F1,x1);
end; {end if}
end; {end while}
xGold:=x1; FGold:=F1;
END;
************************************************************}
begin
ClrScr;
n:=3;
e:=0.00001;
Nitro PDF Trial
www.nitropdf.com

Страницы