diff --git a/ТЕМА2/Hist.jpg b/ТЕМА2/Hist.jpg new file mode 100644 index 0000000..80a336d Binary files /dev/null and b/ТЕМА2/Hist.jpg differ diff --git a/ТЕМА2/Hist2.jpg b/ТЕМА2/Hist2.jpg new file mode 100644 index 0000000..80a336d Binary files /dev/null and b/ТЕМА2/Hist2.jpg differ diff --git a/ТЕМА2/Prog1.m b/ТЕМА2/Prog1.m new file mode 100644 index 0000000..0671f38 --- /dev/null +++ b/ТЕМА2/Prog1.m @@ -0,0 +1,21 @@ +XX=load('dan_vuz.txt'); +size(XX); +X=XX(:,3:13);; +R=corr(X); +[vect,lambda]=eig(X'*X); +Sobst=diag(lambda); +fprintf('Eigenvalues:\n %f \n',Sobst); +fprintf('\n'); +SobMax=Sobst(end); +GlComp=vect(:,end); +Delt=100*SobMax/sum(Sobst); +fprintf('Delta= %d \n ',round(Delt)); +Res=X*GlComp; +fprintf(' Results \n '); +fprintf('%d %f \n ',[XX(:,1),Res] '); +save res.mat Res -mat; +hist(Res,20) +xlabel('Results '); +ylabel('Number of Unis '); +hist(Res,20) + diff --git a/ТЕМА2/Prog2.m b/ТЕМА2/Prog2.m new file mode 100644 index 0000000..1ec0aac --- /dev/null +++ b/ТЕМА2/Prog2.m @@ -0,0 +1,23 @@ +XX=load('dan_vuz.txt'); +size(XX); +X=XX(:,3:13);; +R=corr(X); +[vect,lambda]=eig(X'*X); +Sobst=diag(lambda); +fprintf('Eigenvalues:\n %f \n',Sobst); +fprintf('\n'); +SobMax=Sobst(end); +GlComp=vect(:,end); +Delt=100*SobMax/sum(Sobst); +fprintf('Delta= %d \n ',round(Delt)); +Res=X*GlComp; +fprintf(' Results \n '); +fprintf('%d %f \n ',[XX(:,1),Res] '); +save res.mat Res -mat; +hist(Res,20) +xlabel('Results ') +ylabel('Number of Unis ') +saveas(gcf, 'Hist2.jpg ', 'jpg '); +CorFin=corr(Res,XX(:,2)); +fprintf('Correlation of Results and Money = %f \n',CorFin); + diff --git a/ТЕМА2/assets/figure0.png b/ТЕМА2/assets/figure0.png new file mode 100644 index 0000000..d7f5e2e Binary files /dev/null and b/ТЕМА2/assets/figure0.png differ diff --git a/ТЕМА2/assets/figure1.png b/ТЕМА2/assets/figure1.png new file mode 100644 index 0000000..b933886 Binary files /dev/null and b/ТЕМА2/assets/figure1.png differ diff --git a/ТЕМА2/assets/figure2.png b/ТЕМА2/assets/figure2.png new file mode 100644 index 0000000..4471b03 Binary files /dev/null and b/ТЕМА2/assets/figure2.png differ diff --git a/ТЕМА2/assets/figure3.png b/ТЕМА2/assets/figure3.png new file mode 100644 index 0000000..35f2d99 Binary files /dev/null and b/ТЕМА2/assets/figure3.png differ diff --git a/ТЕМА2/report.md b/ТЕМА2/report.md new file mode 100644 index 0000000..e78c070 --- /dev/null +++ b/ТЕМА2/report.md @@ -0,0 +1,322 @@ +# Отчет по теме 2 + + +Мокеев Никита, А-03-24 + + + +## 1 Настройка текущего каталога + + + +Установил в качестве текущей папки ТЕМА2: + + + +![Скриншот выбора текущей папки](assets/figure0.png) + + + +## 2 + + + +Загрузка данных из текстового файла в матрицу: + +```matlab +>> XX=load('dan_vuz.txt'); + +``` + +Размерность: + +```matlab + +>> size(XX) + +ans = + + + +  290 15 + +``` + + + +*Данные о скольких вузах России представлены в этой матрице?* - 290 + + + +Данные о показателях результативности: + +```matlab + +>> X=XX(:,3:13); + +``` + + + +Матрица корреляций между показателями результативности: + + + +```matlab + +>> R=corr(X); + +``` + + + +Метод главных компонент, получение собственных значений и собственных векторов от квадратичной формы: + +```matlab + +>> [vect,lambda]=eig(X'*X); + +``` + + + +Извлечение диагональных элементов матрицы в вектор: + +```matlab + +>> 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 + +``` + + + +Вывод собственных значений на экран с заголовком: + +```matlab + +>> 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 + +``` + + + +Выделения собственного наибольшего значения и вектора: + +```matlab + +>> 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 + +``` + + + +Доля информации в главной компоненте: + +```matlab + +>> Delt=100*SobMax/sum(Sobst); + +>> fprintf('Delta= %d n ',round(Delt)) + +Delta= 95 + +``` + + + +Вектор оценки обобщенной результативности для каждого ВУЗа: + +```matlab + +>> Res=X*GlComp; + +``` + + + +Вывод заголовка результатов: + +```matlab + +>> fprintf(' Results n ') + + Results + +``` + + + +Вывод таблицы результатов: + +```matlab + +>> fprintf('%d %f n ',[XX(:,1),Res] ') + +``` + + + +Сохранение вектора оценок в бинарном файле: + +```matlab + +>> save res.mat Res -mat + +``` + + + +Построение гистограммы: + +```matlab + +>> hist(Res,20) + +``` + +![Скриншот гистограммы](assets/figure1.png) + +Обозначение осей и повторный вывод: +```matlab +xlabel('Results '); +ylabel('Number of Unis '); +hist(Res,20) +``` + +![Скриншот гистограммы](assets/figure2.png) + +Сохранение гистограммы: +```matlab +saveas(gcf, 'Hist.jpg ', 'jpg ') +``` + +Расчет оценки коререляции: +```matlab +>> CorFin=corr(Res,XX(:,2)) +CorFin = 0.8437 +``` + +И ее вывод: +```matlab +>> fprintf('Correlation of Results and Money = %f n',CorFin) +Correlation of Results and Money = 0.843710 +``` + +## 3 + + +figure2 +Prog2 + + + + + + + + + + + + + diff --git a/ТЕМА2/res.mat b/ТЕМА2/res.mat new file mode 100644 index 0000000..d4cf7fe Binary files /dev/null and b/ТЕМА2/res.mat differ