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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

76
REM Неявные ограничения записываются в строках 7000 - 7100
40 REM функция z=F(x1,x2,...,xn) вычисляется в подпрограмме
REM (строки 5000 - 5100)
60 REM Вычисление значений G1,G2,...,Gm и проверка
65 REM ограничений производится в подпрограмме
REM (строки 6000 - 8000)
80 PRINT "Введите количество неявных ограничений": INPUT M
100 PRINT "Введите количество переменных": INPUT N
120 DIM x(N), Y(N), l(N), u(N), xc(N), x0(N), xr(N), xh(N)
160 K = 2 * N: PP = 0
180 DIM c(K, N), F(K), g(M), Ic(M), Ec(2 * N)
200 PRINT "Введите начальные значения"
220 FOR J = 1 TO N
PRINT "x("; J; ")"; : INPUT x(J)
c(1, J) = x(J): xc(J) = x(J)
NEXT J
240 REM Прочитать значения нижних и верхних границ
REM l(J) - нижние границы переменных x(J)
REM u(J) - верхние границы переменных x(J)
260 FOR J = 1 TO N: READ l(J), u(J): NEXT J
nr = 0
280 REM Включить генератор случайных чисел
RANDOMIZE [1948]
290 x = 1
REM Построение начального комплекса
500 i = 1
520 GOSUB 5000: F(1) = z
600 i = i + 1
620 FOR J = 1 TO N: c(i, J) = l(J) + RND(x) * (u(J) - l(J))
x(J) = c(i, J): NEXT J
640 IM = 1: GOSUB 6000
660 IF Ic = 1 GOTO 720
670 REM Обновить значение "центра тяжести"
680 FOR J = 1 TO N: xc(J) = ((i - 1) * xc(J) + c(i, J)) / i: NEXT J
700 GOTO 760
Nitro PDF Trial
www.nitropdf.com

Страницы