2.9 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
Вывод собственных значений на экран с заголовком:
>> fprintf('Eigenvalues:\n %f \n',Sobst)
Eigenvalues:
22.946585
Eigenvalues:
1931.665464
Eigenvalues:
2593.979592
>> fprintf('\n')
Выделения собственного наибольшего значения и вектора:
>> SobMax=Sobst(end)
SobMax = 7.4946e+06
>> GlComp=vect(:,end)
GlComp =
3.5306e-02
4.6772e-02
4.8953e-02
Доля информации в главной компоненте:
> Delt=100*SobMax/sum(Sobst)
Delt = 95.273
>> fprintf('Delta= %d \n ',round(Delt))
Delta= 95
Вектор оценки обобщенной результативности для каждого ВУЗа:
>> Res=X*GlComp
Res =
9.2542e+01
7.3433e+01
5.8855e+00
Вывод таблицы результатов:
>> fprintf(' Results \n ')
Results
>> fprintf('%d %f \n ',[XX(:,1),Res] ')
197 92.541636
198 73.432513
199 5.885468
200 35.300393
Сохранение вектора оценок в бинарном файле:
>> save res.mat Res -mat
Построение гистограммы:
>> hist(Res,20)
Обозначение осей и повторный вывод:
>> hist(Res,20)
>> xlabel('Results ')
>> ylabel('Number of Unis ')
Сохранение гистограммы:
>> 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