//Зан. млн ч.0,2 0,5 0,76 1,4
//ПРИМЕР 2
disp ("***Primer 2***")
Aish=[6.1 12.5 11.2;
3.1 7.3 14.9;
8.6 7.1 5.3]
Y=[33.2 21.4 50.8]
X=[63 46.7 71.8]
A=[Aish(:,1)/X(1) Aish(:,2)/X(2) Aish(:,3)/X(3)]
Ax=A*X'
Xras=A*X'+Y'
disp("***Primer 3***")
SumA=[sum(A(:,1)) sum(A(:,2)) sum(A(:,3))]
disp("***Primer 4***")
I=[1 0 0; 0 1 0; 0 0 1]
B=(I-A)^-1
XrasB=B*Y'
YX=Y'./X'
Xprom=X-Y
Vipusk_KS=[B(1.1)*Y(1)/X(1)*100 B(1.2)*Y(2)/X(1)*100 B(1.3)*Y(3)/X(1)*100;
B(2.1)*Y(1)/X(2)*100 B(2.2)*Y(2)/X(2)*100 B(2.3)*Y(3)/X(2)*100;
B(3.1)*Y(1)/X(3)*100 B(3.2)*Y(2)/X(3)*100 B(3.3)*Y(3)/X(3)*100]
SumB=[sum(B(:,1)) sum(B(:,2)) sum(B(:,3))]
A1=A(1.1)+A(1.2)+A(1.3)
A1o=[A(1.1)/A1 A(1.2)/A1 A(1.3)/A1]
disp("***Primer 5***")
tz=[0.2 0.5 0.76]
t= tz./X
T=t*B
Tk=T-t
Tkc=T./t
disp("***Primer 6***")
Y=[33.2 21.4 50.8]
SumY=sum(Y')
al=Y/SumY
bet=B*al'
Sumbet=sum(bet')
disp("***Primer 7***")
gam=T*al'
disp("***Primer 8***")
Ys=[11.9 1.2 20.1;
24.1 4.9 -7.6;
34.7 4.2 11.9]
Xs=(I-A)^-1*Ys
SumXs=[sum(Xs(:,1)) sum(Xs(:,2)) sum(Xs(:,3))]
SumXss=[sum(Xs(1,:)) sum(Xs(2,:)) sum(Xs(3,:))]
Qs=T*Ys
Xs_pr=[Xs(1,1)/SumXss(1)*100 Xs(1,2)/SumXss(1)*100 Xs(1,3)/SumXss(1)*100; Xs(2,1)/SumXss(2)*100 Xs(2,2)/SumXss(2)*100 Xs(2,3)/SumXss(2)*100; Xs(3,1)/SumXss(3)*100 Xs(3,2)/SumXss(3)*100 Xs(3,3)/SumXss(3)*100]
disp("***Primer 9***")
M=[70 55 80]
TrudR=t*M'
TrudRfakt=t*X'
Y1=[34 23 51]
Y2=[35 22 52]
Y3=[36 24 54]
Y4=[34 23 52]
Y5=[33.2 23 53]
Y6=[33.2 22 52]
Y7=[33.2 24 53]
Y8=[33/2 23 52]
X1=B*Y1'
X2=B*Y2'
X3=B*Y3'
X4=B*Y4'
X5=B*Y5'
X6=B*Y6'
X7=B*Y7'
X8=B*Y8'
disp("***Primer 10***")
A1=[0.05 0.25 0.1;
0.2 0.15 0.2;
0.25 0.2 0.05]
B1=[1.19 0.4 0.21;
0.37 1.36 0.326;
0.39 0.39 1.176]
Xlg=[20 50 30]
Ylg=(I-A1)*Xlg'
Xlg=[22 60 40]
Ylg=(I-A1)*Xlg'
Yr=(I-A)*M'
TrudRfakt=t*X'
dTrud=t*M'-t*X'
X1=[55.22 48.52 72.48]
X2=[66.27 47.65 73.57]
X3=[68.91 50.92 76.65]
X4=[65.52 48.81 73.65]
X5=[64.87 49.02 74.66]
X6=[64.14 47.44 73.22]
X7=[65.30 50.29 74.93]
X8=[64.57 48.72 73.49]
Y1=(I-A)*X1'
Y2=(I-A)*X2'
Y3=(I-A)*X3'
Y4=(I-A)*X4'
Y5=(I-A)*X5'
Y6=(I-A)*X6'
Y7=(I-A)*X7'
Y8=(I-A)*X8'
disp("***Primer 11***")
a=[-1 -0.253 -0.150;
0 0.832 -0.188;
0 -0.16 0.962]
b=[-0.927*49100 14000+0.048*49100 35500+0.121*49100]
x=inv(a)*b'
Primer 2
Aish =
6.1 12.5 11.2
3.1 7.3 14.9
8.6 7.1 5.3
Y=
33.2 21.4 50.8
X=
63 46.7 71.8
A =
0.0968254 0.2676660 0.1559889
0.0492063 0.1563169 0.2075209
0.1365079 0.1520343 0.0738162
AX =
29.8
25.3
21
Xras =
63
46.7
71.8
***Primer 3***
SumA =
0.2825397 0.5760171 0.4373259
***Primer 4***
SumA =
1.2958827
***Primer 4***
I =
1. 0. 0.
0. 1. 0.
0. 0. 1.
B =
1.1747013 0.4255174 0.2931859
0.1157733 1.2770871 0.3056430
0.1921406 0.2723514 1.1730829
XrasB =
63
46.7
71.8
YX =
0.5269841
0.4582441
0.7075209
Xprom =
29.8 25.3 21
Vipisk_KS =
61.904893 39.902551 94.721945
8.2305629 5.3052424 12.593753
8.8844974 5.7267543 13.594351
SumB =
1.4826152 1.9749559 1.7719118
A1 =
0.2904762
A1 =
0.3333333 0.3333333 0.3333333
***Primer 5***
tz =
0.2 0.5 0.76
t =
0.0031746 0.0107066 0.0105850
T =
0.0070026 0.0179070 0.0166202
Tk =
0.0038280 0.0072003 0.0060352
Tkc =
2.2058044 1.6725126 1.5701707
***Primer 6***
Y =
33.2 21.4 50.8
SumY =
105.4
a1 =
0.3149905 0.2030361 0.4819734
bet =
0.5977230
0.443074
0.6812144
Sumbet =
1.7220114
***Primer 7***
gam =
0.0138520
***Primer 8***
Ys =
11.9 1.2 20.1
24.1 4.9 -7.6
34.7 4.2 11.9
Xs =
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.