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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

78
1610 REM Если явные ограничения нарушены,то
1615 REM поместить точку внутри границ
1620 FOR J = 1 TO N
1640 IF Ec(J) = 1 THEN
xr(J) = l(J) + .00001
x(J) = xr(J)
END IF
1160 IF Ec(J + N) = 1 THEN
xr(J) = u(J) - .00001
x(J) = xr(J)
END IF
1680 NEXT J
1800 IF Ic = 0 GOTO 2000
1810 REM Если неявные ограничения нарушены, то
1815 REM переместить на половину расстояния к "центру тяжести"
1820 FOR l = 1 TO N: xr(l) = (xr(l) + x0(l)) / 2: x(l) = xr(l): NEXT l
1840 GOTO 1490
2000 GOSUB 5000: Fr = z
2010 REM Если новое значение - наихудшее,
2015 REM то переместиться на половину расстояния к точке x0
2018 REM и вычислить новое значение Функции
2020 IF Fr < F(K) GOTO 2400
2040 nr = nr + 1
IF nr >= 10 THEN
nr = 0: GOTO 2400
END IF
FOR l = 1 TO N: xr(l) = (xr(l) + x0(l)) / 2: x(l) = xr(l): NEXT l
2060 GOTO 1490
2400 REM Обновить xc и заменить наихудшую точку новой точкой
2410 F(K) = Fr
2420 FOR l = 1 TO N
2440 xc(l) = K * xc(l) - c(K, l) + xr(l)
2460 xc(l) = xc(l) / K: c(K, l) = xr(l)
2480 NEXT l
2490 REM Упорядочить значения Функции и точки,
2495 REM в которых она вычислена
2500 FOR J = 1 TO K - 1
2520 FOR i = J + 1 TO K
2540 IF F(J) <= F(i) GOTO 2600
Nitro PDF Trial
www.nitropdf.com

Страницы