SHARE
TWEET

Untitled

a guest May 6th, 2019 100 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. clc
  2. clear
  3. printf("------------task_1----------\n")
  4. c=[1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0]';
  5. A=[0 2.2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  6.    0 0 -0.2 0 0 2.2 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  7.    0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0 0 0 0 0 0 0 0 0;
  8.    0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0 0 0 0 0;
  9.    0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0;
  10.    0.588 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  11.    0 0 0 -0.08 0.588 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  12.    0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0 0 0 0 0 0 0 0 0;
  13.    0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0 0 0 0 0;
  14.    0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0;
  15.    1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  16.    0 -2.2 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0;
  17.    0 0 0 0 0 -2.2 0 0 1 1 1 1 0 0 0 0 0 0 0 0;
  18.    0 0 0 0 0 0 0 0 0 -2.2 0 0 1 1 1 1 0 0 0 0;
  19.    0 0 0 0 0 0 0 0 0 0 0 0 0 -2.2 0 0 1 1 1 1;
  20.    1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
  21.    0 -0.5 0 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0;
  22.    0 0 0 0 0 -0.5 0 0 1 0 0 1 0 0 0 0 0 0 0 0;
  23.    0 0 0 0 0 0 0 0 0 -0.5 0 0 1 0 0 1 0 0 0 0;
  24.    0 0 0 0 0 0 0 0 0 0 0 0 0 -0.5 0 0 1 0 0 1];
  25. b=[5600 5600 5600 5600 5600 800 800 800 800 800 5500 0 0 0 0 2750 0 0 0 0]';
  26. [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1)
  27. printf("------------task_2_1----------\n")
  28. j=1;
  29. k=1000
  30. b(11)=0;
  31. b(16)=0;
  32. while(b(11)!=[10000])
  33.   b(11)=b(11)+k;
  34.   b(16)=b(16)+k/2;
  35.   [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  36.   y(j)=f_max;
  37.   j=j+1;
  38.  endwhile
  39. disp(y);
  40. x=1000:1000:10000;
  41. plot(x,y);
  42. title('Кол-во станков')
  43. legend('boxoff');
  44. legend('Кол-во стали, потраченное на производство станков',1)
  45. xlabel('Изначальное кол-во станков')
  46. ylabel('Потраченное кол-во стали')
  47. grid on;
  48. printf("------------task_2_2----------\n")
  49. j=1;
  50. k=1000
  51. b(11)=0;
  52. b(16)=0;
  53. while(b(11)!=[10000])
  54.   b(11)=b(11)+k;
  55.   b(16)=b(16)+k/2;
  56.   [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  57.   y(j)=f_max;
  58.   y1(j)=x_max(1)+x_max(5)+x_max(9)+x_max(13)+x_max(17);
  59.   y2(j)=x_max(2)+x_max(6)+x_max(10)+x_max(14)+x_max(18);
  60.   y3(j)=x_max(3)+x_max(7)+x_max(11)+x_max(15)+x_max(19);
  61.   y4(j)=x_max(4)+x_max(8)+x_max(12)+x_max(16)+x_max(20);
  62.   j=j+1;
  63.  endwhile
  64.  x=1000:1000:10000;
  65.  figure()
  66.  plot(x,y1,x,y2,x,y3,x,y4);
  67.  title('Кол-во станков')
  68. legend('boxoff');
  69. legend('Кол-во стали, потраченное на производство станков','Кол-во стали, потраченное на производство стали','Кол-во стали, потраченное на улучшение сталелинейного завода','Кол-во стали, потраченное на улучшение станкостроительного завода',4)
  70. xlabel('Изначальное кол-во станков')
  71. ylabel('Потраченное кол-во стали')
  72. grid on;
  73. printf("------------task_3_1----------\n")
  74. j=1;
  75. k=100
  76. b(6)=0;
  77. b(7)=0;
  78. b(8)=0;
  79. b(9)=0;
  80. b(10)=0;
  81. while(b(6)!=[5000]);
  82.   b(6)=b(6)+k;
  83.   b(7)=b(7)+k;
  84.   b(8)=b(8)+k;
  85.   b(9)=b(9)+k;
  86.   b(10)=b(10)+k;
  87.   [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  88.   y(j)=f_max;
  89.   j=j+1;
  90.  endwhile
  91.  disp(y);
  92.  x=100:100:5000;
  93.  figure()
  94.  plot(x,y);
  95.  title('Кол-во станков')
  96. legend('boxoff');
  97. legend('Кол-во стали, потраченное на производство станков',1)
  98. xlabel('Мощность производства')
  99. ylabel('Кол-во стали')
  100. grid on;
  101. printf("------------task_3_1----------\n")
  102. j=1;
  103. k=100
  104. b(6)=0;
  105. b(7)=0;
  106. b(8)=0;
  107. b(9)=0;
  108. b(10)=0;
  109. while(b(6)!=[5000]);
  110.   b(6)=b(6)+k;
  111.   b(7)=b(7)+k;
  112.   b(8)=b(8)+k;
  113.   b(9)=b(9)+k;
  114.   b(10)=b(10)+k;
  115.   [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  116.   y(j)=f_max;
  117.   y1(j)=x_max(1)+x_max(5)+x_max(9)+x_max(13)+x_max(17);
  118.   y2(j)=x_max(2)+x_max(6)+x_max(10)+x_max(14)+x_max(18);
  119.   y3(j)=x_max(3)+x_max(7)+x_max(11)+x_max(15)+x_max(19);
  120.   y4(j)=x_max(4)+x_max(8)+x_max(12)+x_max(16)+x_max(20);
  121.   j=j+1;
  122.  endwhile
  123.  disp(y);
  124.  x=100:100:5000;
  125.  figure()
  126.  plot(x,y1,x,y2,x,y3,x,y4);
  127.  title('Кол-во станков')
  128. legend('boxoff');
  129. legend('Кол-во стали, потраченное на производство станков','Кол-во стали, потраченное на производство стали','Кол-во стали, потраченное на улучшение сталелинейного завода','Кол-во стали, потраченное на улучшение станкостроительного завода',4)
  130. xlabel('Мощность производства')
  131. ylabel('Кол-во стали')
  132. grid on;
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top