Ткачев Ю.А. Открытые и закрытые числовые системы

67 end; Statisk(nf,m,X1,mx1,sigmax1,asimx1); Statisk(nf,m,X2,mx2,sigmax2,asimx2); Statisk(nf,m,X3,mx3,sigmax3,asimx3); {сохраним для сумм} sigmax1S:=sigmax1;sigmax2S:=sigmax2; sigmax3S:=sigmax3; {расчёт сумм компонентов по пробам} for i:=1 to nf do SumCompx[i]:=X1[i]+X2[i]+X3[i]; Correl(nf,mx1,mx2,sigmax1,sigmax2,X1,X2,rx12); Correl(nf,mx1,mx3,sigmax1,sigmax3,X1,X3,rx13); Correl(nf,mx2,mx3,sigmax2,sigmax3,X2,X3,rx23); skvox1:=sqrt(sigmax1);skvox2:=sqrt(sigmax2); skvox3:=sqrt(sigmax3); {печать статистик исходных данных} Writeln(f,’статистики исходных данных’); Writeln(f); Writeln(f,’Комп Средние Дисперсия Скво Asimx ’); Writeln(f,’ X1 ’,mx1:10:3, sigmax1:10:3, skvox1:10:3,asimx1:10:3); Writeln(f,’ X2 ’,mx2:10:3, sigmax2:10:3, skvox2:10:3,asimx2:10:3); Writeln(f,’ X3 ’,mx3:10:3, sigmax3:10:3, skvox3:10:3,asimx3:10:3); Writeln(f,’ rX rx12 rx13 rx23 ’); Writeln(f,’ ’, rx12:10:3,rx13:10:3, rx23:10:3); {расчёт компонентов закрытой системы} for i:=1 to nf do Y1[i]:=X1[i]/SumCompx[i]; { Y - доля комп.X1[i] во всей совокупности компонентов в пробе i} for i:=1 to nf do Y2[i]:=X2[i]/SumCompx[i]; for i:=1 to nf do Y3[i]:=X3[i]/SumCompx[i]; if nf>50 then i:=i else begin {печать данных по долям} Writeln(f); Writeln(f,’ массивы долей’); Writeln(f); Writeln(f,’ Y1 Y2 Y3 ’); for i:=1 to nf do Writeln(f,Y1[i]:10:3,Y2[i]:10:3, Y3[i]:10:3); Writeln(f);