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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

82
ELSE
i = x3 - x1: x0 = x1: x1 = x2: x2 = x0 + t2 * i: F1 = F2
CALL Func(x2, F2)
END IF
LOOP WHILE i > E * (B - A)
x = x1: F = F1
END SUB
SUB Func (x, F)
SHARED k, SA, So, CA
x1 = 1.34 - .34 * x - SQR((1.34 - .34 * x) ^ 2 - x)
tau = 2.68 * x1 / (k * (1 - x1))
F = SA / x + tau * So / (CA * x)
END SUB
2. Программа нахождения минимума функции(Turbo Pascal)
Program Gold;
var A,B,E,x,F:real;
{****************************************************************}
Procedure GoldMin (A, B, E:real; var x,F:real);
{ "Метод золотого сечения"
Процедура GoldMin производит поиск интервала [X1,X2],
в котором лежит точка минимума функции F(x) }
{F(x) определяется в процедуре Fх
А,В - интервал возможного варьирования переменной х
А - нижний предел; В - верхний предел
Е - относительная погрешность вычислений
(относительно интервала А - В) (Х2-Х1)<=Е*(В-А) }
var i,t1,t2,x0,x1,x2,x3,F0,F1,F2,F3 :real;
{****************************************************************}
Procedure Fx (var F:real; var x:real);
begin
F:=(x-2)*(x-2)+1;
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);
i:=x3 - x1;
While (i>E*(B-A)) do begin
Nitro PDF Trial
www.nitropdf.com

Страницы