j = j - 1; if j <= 1e-8

// Решение получено - выйти из цикла break

end end // if

x(j) = x(j) + 1; end // while 1

write(6,[' *********'

' Ответ ' ' *********'])

write(6, ' f_opt = ' + string(f_opt)) // печать в одну строку x_opt

10.2. Задача о загрузке

10.2.1. Динамическое программирование – вычислить C и R

//function RiukDiCR(c, a, b )

//Решение задачи о загрузке самолета методом ДП.

//==== 30 марта 2011 года =====Визгунов Н.П.======

clc, clear, mode(0), lines(0)

n= length(a);

55

u9 = 1 + floor(b / min(a)); // u9 - Кол. управлений

// 1) Очевидный способ вычисления С и R

C = -%inf * ones(u9, n); R = C;

for u = 1 : u9 // (u-1) - это количество предметов for j = 1 : n

temp = a(j) * (u -1); // j - номер предмета

if temp <= b // если предметы помещаются в самолет

C(u, j) = temp;

R(u, j) = c(j) * (u -1);

end

end

end C; R;

// 2) Вычисления С и R в духе SciLab

C = (0 : u9 -1)' * a

I = find(C > y1)

C(I) = -%inf

R = (0 : u9 -1)' * c;

R(I) = -%inf

end

end

end

[fs, index] = max(Fsu, 'c'); F(:,j) = fs;

X(:,j) = index - 1;

//Печать шапки и таблицы

//======================

u7 = 1 + floor(b / min(a(j)));

line50 = part('-', ones(1, 22)) + '|' + ...

part('-', ones(1, 3 + 6 * u7));

56

disp(Temp(:, 1 : 3 + u7))

write(6, line50)

end

//Вычислить и напечать ответ

//==========================

x_opt = zeros(1, n); s = s9;

for j = 1 : n

x_opt(j) = X(s, j);

s = s - a(j) * x_opt(j);

end

c, a, x_opt

write(6,' f_opt = ' + string(F(s9, 1)) + ...

'y1 = ' + string(y1))

10.2.2.Динамическое программирование – без вычисления C и R

//function RiukDi(c, a, b )

//Программа RiukDi.sce. Нет вычислений С и R !!!

//Решение задачи о загрузке самолета методом ДП.

//=====30 марта 2011 года=====Визгунов Н.П.=====

clc, clear, lines(0, 100), mode(0)

//доходы =

c = [ 65 80 30 ]; //вес предметов =

a = [ 2 3 1 ]; //грузоподъемность самолета = b = 5; y1 = b;

n= length(a);

end

end

[fs, index] = max(Fsu, 'c'); F(:,j) = fs;

X(:,j) = index - 1;

//Печать шапки и таблицы

//======================

u7 = 1 + floor(b / min(a(j)));

57

line50 = part('-', ones(1, 22)) + '|' + ...

part('-', ones(1, 3 + 6 * u7));

// line50 = code2str(str2code('-') * ones(1,58));

write(6, line50)

end

//Вычислить и напечать ответ

//==========================

x_opt = zeros(1, n); s = s9;

for j = 1 : n

x_opt(j) = X(s, j);

s = s - a(j) * x_opt(j);

end

c, a, x_opt

write(6,' f_opt = ' + string(F(s9,1)) + ...

' y1= ' + string(y1)) //endfunction

10.2.3. Полный перебор

//function RiukPe(c, a, b)

//программа RjukPe.sce Лабораторная номер 3

//=====30 марта 2011 года=====Визгунов Н.П.

clc, clear, lines(0), mode(0)

n= length(a);

if a * x' <= b // если предметы войдут в самолет f = c * x';

if f_opt < f

disp('f_opt < f <<<<<<<<<<<<<')

58