3.7 KiB
Отчет по теме 2
Мокеев Никита, А-03-24
1 Настройка текущего каталога
Установил в качестве текущей папки ТЕМА2:
2
Загрузка данных из текстового файла в матрицу:
>> XX=load('dan_vuz.txt');
Размерность:
>> size(XX)
ans =
290 15
Данные о скольких вузах России представлены в этой матрице? - 290
Данные о показателях результативности:
>> X=XX(:,3:13);
Матрица корреляций между показателями результативности:
>> R=corr(X);
Метод главных компонент, получение собственных значений и собственных векторов от квадратичной формы:
>> [vect,lambda]=eig(X'*X);
Извлечение диагональных элементов матрицы в вектор:
>> Sobst=diag(lambda)
Sobst =
2.2947e+01
1.9317e+03
2.5940e+03
3.4573e+03
5.6252e+03
8.6721e+03
1.8915e+04
4.7523e+04
5.7484e+04
2.2565e+05
7.4946e+06
Вывод собственных значений на экран с заголовком:
>> fprintf('Eigenvalues:n %f n',Sobst)
Eigenvalues:
22.946585
Eigenvalues:
1931.665464
Eigenvalues:
2593.979592
Eigenvalues:
3457.339562
Eigenvalues:
5625.151474
Eigenvalues:
8672.065947
Eigenvalues:
18914.627989
Eigenvalues:
47522.678185
Eigenvalues:
57483.681267
Eigenvalues:
225653.068540
Eigenvalues:
7494628.795394
Выделения собственного наибольшего значения и вектора:
>> SobMax=Sobst(end)
SobMax = 7.4946e+06
>> GlComp=vect(:,end)
GlComp =
3.5306e-02
4.6772e-02
4.8953e-02
6.1556e-01
2.4277e-01
7.3685e-01
9.5893e-02
1.6945e-04
1.7911e-02
5.9523e-02
1.7425e-02
Доля информации в главной компоненте:
>> Delt=100*SobMax/sum(Sobst);
>> fprintf('Delta= %d n ',round(Delt))
Delta= 95
Вектор оценки обобщенной результативности для каждого ВУЗа:
>> Res=X*GlComp;
Вывод заголовка результатов:
>> fprintf(' Results n ')
Results
Вывод таблицы результатов:
>> fprintf('%d %f n ',[XX(:,1),Res] ')
Сохранение вектора оценок в бинарном файле:
>> save res.mat Res -mat
Построение гистограммы:
>> hist(Res,20)
Обозначение осей и повторный вывод:
xlabel('Results ');
ylabel('Number of Unis ');
hist(Res,20)
Сохранение гистограммы:
saveas(gcf, 'Hist.jpg ', 'jpg ')
Расчет оценки коререляции:
>> CorFin=corr(Res,XX(:,2))
CorFin = 0.8437
И ее вывод:
>> fprintf('Correlation of Results and Money = %f n',CorFin)
Correlation of Results and Money = 0.843710
3
figure2 Prog2