Модель рыбной отрасли. Построение модели рынка. Математическая модель рынка с четырьмя потребителями и производителями, страница 5

- 1.    0.    0.    0.    0.    0.    0.    0.    0. 

-->a0=[1 -27/d1 -50/d1 0 0 0 0 0 0;

-->    1 0 0 -26/d2 -15/d2 0 0 0 0;

-->    1 0 0 0 0 -35/d3 -15/d3 0 0;

-->    1 0 0 0 0 0 0 -35/d4 -28/d4;

-->    1 -70/d5 0 0 0 -55/d5 0 0 0;

-->    1 0 -95/d6 0 0 0 -35/d6 0 0;

-->    1 0 0 -46/d7 0 0 0 -85/d7 0;

-->    1 0 0 0 -35/d8 0 0 0 -48/d8;

-->    0 70 0 0 0 55 0 0 0;

-->    0 -70 0 0 0 -55 0 0 0;

-->    0 0 95 0 0 0 35 0 0;

-->    0 0 -95 0 0 0 -35 0 0;

-->    0 0 0 46 0 0 0 85 0;

-->    0 0 0 -46 0 0 0 -85 0;

-->    0 0 0 0 35 0 0 0 48;

-->    0 0 0 0 -35 0 0 0 -48;

-->    0 43 45 0 0 0 0 0 0;

-->    0 0 0 20 20 0 0 0 0;

-->    0 0 0 0 0 20 20 0 0;

-->    0 0 0 0 0 0 0 50 20]

a0  =

column 1 to 7

1.  - 0.0023684  - 0.0043859    0.           0.           0.           0.        

1.    0.           0.         - 0.0028621  - 0.0016512    0.           0.        

1.    0.           0.           0.           0.         - 0.0026101  - 0.0011186 

1.    0.           0.           0.           0.           0.           0.        

1.    0.0116669    0.           0.           0.           0.0091669    0.        

1.    0.           0.0135718    0.           0.           0.           0.0050001 

1.    0.           0.           0.0070771    0.           0.           0.        

1.    0.           0.           0.           0.0050002    0.           0.        

0.    70.          0.           0.           0.           55.          0.        

0.  - 70.          0.           0.           0.         - 55.          0.        

0.    0.           95.          0.           0.           0.           35.       

0.    0.         - 95.          0.           0.           0.         - 35.       

0.    0.           0.           46.          0.           0.           0.        

0.    0.           0.         - 46.          0.           0.           0.        

0.    0.           0.           0.           35.          0.           0.         

0.    0.           0.           0.         - 35.          0.           0.        

0.    43.          45.          0.           0.           0.           0.        

0.    0.           0.           20.          20.          0.           0.        

0.    0.           0.           0.           0.           20.          20.       

0.    0.           0.           0.           0.           0.           0.         

column 8 to 9

0.           0.        

0.           0.        

0.           0.        

- 0.0025167  - 0.0020133 

0.           0.        

0.           0.        

0.0130773    0.        

0.           0.0068574 

0.           0.        

0.           0.        

0.           0.        

0.           0.        

85.          0.        

- 85.          0.        

0.           48.       

0.         - 48.       

0.           0.        

0.           0.        

0.           0.        

50.          20.       

-->b0=[-f1min/d1 -f2min/d2 -f3min/d3 -f4min/d4 f5max/d5 f6max/d6 f7max/d7 f8max/d8 14000 -8000 19000 -11000 16500 -10000 19000 -12000 14000 8000 19000 11000]

b0  =

column 1 to 7

- 0.0000011  - 0.3540009    0.0000178  - 0.0629176    2.3333714    2.4286209    2.38466 

column  8 to 15

2.714362    14000.  - 8000.    19000.  - 11000.    16500.  - 10000.    19000. 

column 16 to 20

- 12000.    14000.    8000.    19000.    11000. 

-->lb0=[0 0 0 0 0 0 0 0 0]

lb0  =

0.    0.    0.    0.    0.    0.    0.    0.    0. 

-->ub0=[1 300 300 300 300 300 300 300 300]

ub0  =

1.    300.    300.    300.    300.    300.    300.    300.    300. 

-->[x0,l0]=karmarkar([],[],L',[],[],[],[],[],a0,b0',lb0',ub0')

l0  =

- 0.5021956 

x0  =

0.5021956 

60.945263 

81.591913 

265.98529 

57.481546 

122.19219 

163.80883 

0.0033938 

280.68092 

-->Fx0=[f(1,:)*x0(2:9) f(2,:)*x0(2:9) f(3,:)*x0(2:9) f(4,:)*x0(2:9) f(5,:)*x0(2:9) f(6,:)*x0(2:9) f(7,:)*x0(2:9) f(8,:)*x0(2:9)]

Fx0  =

column 1 to 6

- 5725.1177  - 7777.8407  - 6733.8593  - 7859.1844    10986.739    13484.541 

column 7 to 8

12235.612    15484.538 

-->lx0=[(-Fx0(1)-f1min)/d1 (-Fx0(2)-f2min)/d2 (-Fx0(3)-f3min)/d3 (-Fx0(4)-f4min)/d4 (Fx0(5)+f5max)/d5 (Fx0(6)+f6max)/d6 (Fx0(7)+f7max)/d7 (Fx0(8)+f8max)/d8]

lx0  =

column 1 to 6

0.5021978    0.5022010    0.5021987    0.5021988    0.5022064    0.5022024 

column 7 to 8

0.5022079    0.5022064