Смекни!
smekni.com

Численные методы решения задач управления технологическими процессами (стр. 4 из 4)

a=[-4 -6 -15;-2 -2 0;-5 -3 -4;-7 -3 -12]; »b=[-250-60-100-220];

»f=[44 35 100]

f =

44 35 100

» [x,y,z]=linprog(f,a,b) Optimization terminated. x= 13.2143

16.7857

6.4286

y =

1.8118e+003 z= 1

Задача №3 Вариант№13

N1 N2 N3 N4 N5 В1
А1 2 5 7 9 12 4000
А2 4 9 16 21 28 6000
АЗ 7 12 18 27 34 8000
А4 10 18 26 34 42 10000
А5 14 23 32 44 51 12000
А6 18 29 39 50 63 14000
С 10 11 12 13 14

Расчет производится в приложении МАТLАВ.

a=[2,5,7,9,12;4,9,16,21,28;7,12/18,27,34;10,18,26,34,42;14l23,32,44,51;18,29,39,50,63]

b=[4000;60000;8000;10000;12000;14000]

lb=zeros(5,l)

f=[-10-ll -12-13-14]

»[x,y,z]=linprog(f,a,b,[],[],lb)

Optimization terminated.

x =

777.7778

0.0000

0.0000

0.0000

0.0000

y =

-7.7778e+003

z =

1

Задача №4 Вариант №8

12.

Уваривание рецептурной смеси в змеевиковом уварочном аппарате.

Общий вид зависимости: τвар=f(P пара)

Результаты моделирования: y=178-865,5*x+1629,4*x^2-1025*x^3

Минимизировать время уваривания смеси до влажности 16% при ограничении на давлении в магистрали пара, обусловленного экономической эффективность.

Ограничения: 0,2 <Р пара<0,7 Мпа; f`(Р пара)≤-44

Расчет производится в приложении МАТLАВ.

function f=myfun(x)

function[c,ceq]=ogr(x); c=-865,5+3258,8*x-3075*x^2+44;

ceq=[];

»[x,y,z]=fmincon(@myfun,[1],[],[],[],[],[0,7],[0,2],@ogr)

x =

0

Y =

-687

z =

1

13.

Уваривание рецептурной смеси в змеевиковом уварочном аппарате.

Общий вид зависимости: y отк=f(P апп)

Результаты моделирования: y=47,5-2,01*x+0,03*x^2-0,000015*x^3

Повысить качество смеси, выраженное в заданной взякости 20 Па*с в зависимости от ограничения для температуры в аппарате и частотой вращения мешалки.

Ограничения: 60 <Т апп<76 °С; 15<w лоп<35 об/мин;

Расчет производится в приложении МАТLАВ.

function f=myfun(x)

f=47,5-2,01*x+0,03*x^2-0,00015*x^3;

function[c,ceq]=ogr(x); c=-2,01+0,06*x-0,00045*x^2; ceq=[];

»[x,y;z]=fmincon(@myfun,[l],[],[],[],[],[65],[93],@ogr)

x =

65.0001

y =

47

z =

1