ответвлено от main/it-labs
33 строки
814 B
Matlab
33 строки
814 B
Matlab
fp=fopen('prtcl.txt', 'w');
|
|
XX=load('dan_vuz.txt');
|
|
size(XX);
|
|
X=XX(:,3:13);
|
|
R = corr(X);
|
|
[vect,lambda]=eig(X'*X);
|
|
Sobst=diag(lambda);
|
|
fprintf(fp, 'Eigenvalues:\n%f\n', Sobst);
|
|
fprintf(fp, '\n');
|
|
SobMax = Sobst(end);
|
|
GlComp = vect(:, end);
|
|
Delt = 100*SobMax/sum(Sobst);
|
|
fprintf('Delta=%d\n', round(Delt));
|
|
Res = X * GlComp;
|
|
fprintf(fp, 'Results\n');
|
|
fprintf(fp, '%d %f \n', [XX(:, 1), Res]');
|
|
MinRes=min(Res);
|
|
MaxRes=max(Res);
|
|
MeanRes=mean(Res);
|
|
StdRes=std(Res);
|
|
fprintf(fp, 'Min=%f\n', MinRes);
|
|
fprintf(fp, 'Max=%f\n', MaxRes);
|
|
fprintf(fp, 'Mean=%f\n', MeanRes);
|
|
fprintf(fp, 'Std=%f\n', StdRes);
|
|
save res.mat Res -mat;
|
|
hist(Res, 20);
|
|
xlabel('Results');
|
|
ylabel('Number of Unis');
|
|
saveas(gcf,'Hist.jpg','jpg');
|
|
CorFin = corr(Res, XX(:,2));
|
|
fprintf(fp, 'Correlation of Results and Money = %f\n', CorFin);
|
|
fclose(fp);
|