{"id":446015,"date":"2025-01-24T21:00:29","date_gmt":"2025-01-24T21:00:29","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=446015"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=446015","title":{"rendered":"<span>Armadillo: \u043c\u0430\u0442\u0440\u0438\u0446\u044b, \u043a\u0443\u0431\u044b \u0438 \u0440\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u0435 \u0434\u0430\u043d\u043d\u044b\u0435 \u043d\u0430 C++<\/span>"},"content":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p><em>\u041f\u0440\u0438\u0432\u0435\u0442, \u0425\u0430\u0431\u0440!<\/em><\/p>\n<p>\u0420\u0430\u0431\u043e\u0442\u0430 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438 \u0432\u00a0\u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u043e\u043c C++\u00a0\u2014 \u044d\u0442\u043e \u0431\u043e\u043b\u044c, \u0441\u0442\u0440\u0430\u0434\u0430\u043d\u0438\u0435 \u0438 \u043c\u0430\u0441\u0441\u0430 \u043d\u0435\u043d\u0443\u0436\u043d\u043e\u0433\u043e \u043a\u043e\u0434\u0430. \u041f\u043e\u0447\u0435\u043c\u0443? \u041f\u043e\u0442\u043e\u043c\u0443 \u0447\u0442\u043e\u00a0\u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u0430\u044f \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 STL \u0432\u043e\u043e\u0431\u0449\u0435 \u043d\u0435\u00a0\u043f\u0440\u0435\u0434\u043d\u0430\u0437\u043d\u0430\u0447\u0435\u043d\u0430 \u0434\u043b\u044f\u00a0\u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0430\u043b\u0433\u0435\u0431\u0440\u044b. \u041d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u0435\u0441\u043b\u0438 \u0432\u044b \u0445\u043e\u0442\u0438\u0442\u0435 \u0443\u043c\u043d\u043e\u0436\u0438\u0442\u044c \u0434\u0432\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b, \u0432\u0430\u043c \u043f\u0440\u0438\u0434\u0451\u0442\u0441\u044f \u043f\u0438\u0441\u0430\u0442\u044c \u043a\u0430\u0441\u0442\u043e\u043c\u043d\u044b\u0439 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c, \u043e\u0442\u043b\u0430\u0436\u0438\u0432\u0430\u0442\u044c \u0435\u0433\u043e, \u0430\u00a0\u0437\u0430\u0442\u0435\u043c \u0435\u0449\u0451 \u0440\u0430\u0437 \u043f\u0438\u0441\u0430\u0442\u044c \u0435\u0433\u043e, \u043a\u043e\u0433\u0434\u0430 \u043f\u043e\u0439\u043c\u0451\u0442\u0435, \u0447\u0442\u043e\u00a0\u043e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043d\u0435\u00a0\u0442\u0430.<\/p>\n<p>\u041d\u043e\u00a0\u0435\u0441\u0442\u044c \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 Armadillo \u2014 \u0438 \u0441\u00a0\u043d\u0435\u0439 \u0432\u0441\u0451 \u0438\u043d\u0430\u0447\u0435. \u0412\u043e\u0442 \u0435\u0451 \u0433\u043b\u0430\u0432\u043d\u044b\u0435 \u0444\u0438\u0447\u0438:<\/p>\n<ol>\n<li>\n<p><strong>\u041f\u0440\u043e\u0441\u0442\u043e\u0442\u0430:<\/strong> \u0432\u044b \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442\u0435 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438, \u043a\u0430\u043a\u00a0\u0441\u00a0\u043e\u0431\u044b\u0447\u043d\u044b\u043c\u0438 \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u043c\u0438. \u041d\u0438\u043a\u0430\u043a\u0438\u0445 \u00ab\u0434\u0436\u0443\u043d\u0433\u043b\u0435\u0439\u00bb \u0438\u0437\u00a0\u0446\u0438\u043a\u043b\u043e\u0432.<\/p>\n<\/li>\n<li>\n<p><strong>\u0421\u043a\u043e\u0440\u043e\u0441\u0442\u044c:<\/strong> \u043f\u043e\u0434\u0434\u0435\u0440\u0436\u043a\u0430 OpenBLAS, MKL \u0438 LAPACK \u0434\u0435\u043b\u0430\u0435\u0442 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u043e\u043b\u043d\u0438\u0435\u043d\u043e\u0441\u043d\u044b\u043c\u0438.<\/p>\n<\/li>\n<li>\n<p><strong>\u0412\u043e\u0437\u043c\u043e\u0436\u043d\u043e\u0441\u0442\u0438:<\/strong> \u0432\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u044b\u0435 \u043c\u0435\u0442\u043e\u0434\u044b \u0434\u043b\u044f\u00a0\u043a\u043b\u0430\u0441\u0442\u0435\u0440\u0438\u0437\u0430\u0446\u0438\u0438, PCA, \u0440\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u0445 \u043c\u0430\u0442\u0440\u0438\u0446 \u0438 \u0434\u0430\u0436\u0435 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u0433\u043e \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u044f.<\/p>\n<\/li>\n<li>\n<p><strong>\u0423\u043d\u0438\u0432\u0435\u0440\u0441\u0430\u043b\u044c\u043d\u043e\u0441\u0442\u044c:<\/strong> \u043c\u043e\u0436\u043d\u043e \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u043a\u0430\u043a\u00a0\u0441\u00a0\u043f\u043b\u043e\u0442\u043d\u044b\u043c\u0438, \u0442\u0430\u043a \u0438 \u0441\u00a0\u0440\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u043c\u0438 \u0434\u0430\u043d\u043d\u044b\u043c\u0438.<\/p>\n<\/li>\n<\/ol>\n<p>\u0422\u0435\u043f\u0435\u0440\u044c \u043f\u0435\u0440\u0435\u0439\u0434\u0451\u043c \u0441\u0438\u043d\u0442\u0430\u043a\u0441\u0438\u0441\u0443.<\/p>\n<h3>\u0411\u0430\u0437\u043e\u0432\u044b\u0439 \u0441\u0438\u043d\u0442\u0430\u043a\u0441\u0438\u0441 Armadillo<\/h3>\n<p>\u0414\u043b\u044f\u00a0\u0440\u0430\u0431\u043e\u0442\u044b \u0441\u00a0Armadillo \u0432\u0430\u043c \u043d\u0443\u0436\u043d\u043e \u043f\u043e\u0434\u043a\u043b\u044e\u0447\u0438\u0442\u044c \u0435\u0451 \u0437\u0430\u0433\u043e\u043b\u043e\u0432\u043e\u0447\u043d\u044b\u0439 \u0444\u0430\u0439\u043b:<\/p>\n<pre><code class=\"cpp\">#include &lt;armadillo&gt;<\/code><\/pre>\n<p>\u0427\u0442\u043e\u0431\u044b \u0438\u0437\u0431\u0435\u0436\u0430\u0442\u044c \u043f\u043e\u0441\u0442\u043e\u044f\u043d\u043d\u043e\u0433\u043e \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u0438\u044f <code>arma::<\/code>, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u043e \u0438\u043c\u0451\u043d:<\/p>\n<pre><code class=\"cpp\">using namespace arma;<\/code><\/pre>\n<p>\u0412\u043e\u0442 \u0438 \u0432\u0441\u0451. \u0422\u0435\u043f\u0435\u0440\u044c \u0432\u0441\u0435 \u0433\u043e\u0442\u043e\u0432\u043e \u043a\u00a0\u0440\u0430\u0431\u043e\u0442\u0435 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438, \u0432\u0435\u043a\u0442\u043e\u0440\u0430\u043c\u0438 \u0438 \u0434\u0440\u0443\u0433\u0438\u043c\u0438 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c\u0438 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438.<\/p>\n<h4>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446<\/h4>\n<p>\u0421\u0430\u043c\u044b\u0439 \u043f\u043e\u043f\u0443\u043b\u044f\u0440\u043d\u044b\u0439 \u043a\u043b\u0430\u0441\u0441 \u0432\u00a0Armadillo\u00a0\u2014 \u044d\u0442\u043e <code>mat<\/code>. \u041e\u043d \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0434\u043b\u044f\u00a0\u0440\u0430\u0431\u043e\u0442\u044b \u0441\u00a0\u043f\u043b\u043e\u0442\u043d\u044b\u043c\u0438 \u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043f\u0443\u0441\u0442\u043e\u0439 \u043c\u0430\u0442\u0440\u0438\u0446\u044b:<\/p>\n<pre><code class=\"cpp\">mat A;  \/\/ \u041f\u043e\u043a\u0430 \u0431\u0435\u0437 \u0440\u0430\u0437\u043c\u0435\u0440\u043e\u0432<\/code><\/pre>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b \u0444\u0438\u043a\u0441\u0438\u0440\u043e\u0432\u0430\u043d\u043d\u043e\u0433\u043e \u0440\u0430\u0437\u043c\u0435\u0440\u0430:<\/p>\n<pre><code class=\"cpp\">mat A(3, 4);  \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 3x4, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u043c\u0443\u0441\u043e\u0440\u043e\u043c (\u0434\u043e Armadillo 10.4) \u0438\u043b\u0438 \u043d\u0443\u043b\u044f\u043c\u0438 (\u0441 \u0432\u0435\u0440\u0441\u0438\u0438 10.5).<\/code><\/pre>\n<p><strong>\u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0438\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438:<\/strong><\/p>\n<ul>\n<li>\n<p>\u0421\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0435 \u0447\u0438\u0441\u043b\u0430 (\u0440\u0430\u0432\u043d\u043e\u043c\u0435\u0440\u043d\u043e\u0435 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435):<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randu);  \/\/ \u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0430 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u043c\u0438 \u0447\u0438\u0441\u043b\u0430\u043c\u0438 \u043e\u0442 0 \u0434\u043e 1<\/code><\/pre>\n<ul>\n<li>\n<p>\u0421\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0435 \u0447\u0438\u0441\u043b\u0430 (\u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435):<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randn);  \/\/ \u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0430 \u0447\u0438\u0441\u043b\u0430\u043c\u0438 \u0441 \u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u044b\u043c \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435\u043c<\/code><\/pre>\n<ul>\n<li>\n<p>\u0415\u0434\u0438\u043d\u0438\u0447\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat I = eye&lt;mat&gt;(4, 4);  \/\/ \u0415\u0434\u0438\u043d\u0438\u0447\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430 4x4<\/code><\/pre>\n<ul>\n<li>\n<p>\u041c\u0430\u0442\u0440\u0438\u0446\u0430 \u0438\u0437\u00a0\u043d\u0443\u043b\u0435\u0439 \u0438\u043b\u0438\u00a0\u0435\u0434\u0438\u043d\u0438\u0446:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat Z = zeros&lt;mat&gt;(5, 5);  \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 5x5, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u043d\u0443\u043b\u044f\u043c\u0438 mat O = ones&lt;mat&gt;(5, 5);   \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 5x5, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u0435\u0434\u0438\u043d\u0438\u0446\u0430\u043c\u0438<\/code><\/pre>\n<h4>\u0414\u043e\u0441\u0442\u0443\u043f \u043a \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c<\/h4>\n<p>\u041c\u043e\u0436\u043d\u043e \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u044c \u0434\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c \u043c\u0430\u0442\u0440\u0438\u0446\u044b \u0442\u0430\u043a\u00a0\u0436\u0435, \u043a\u0430\u043a\u00a0\u0432\u00a0Python \u0438\u043b\u0438\u00a0MATLAB:<\/p>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randu); double value = A(0, 1);  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442 \u0438\u0437 \u043f\u0435\u0440\u0432\u043e\u0439 \u0441\u0442\u0440\u043e\u043a\u0438 \u0438 \u0432\u0442\u043e\u0440\u043e\u0433\u043e \u0441\u0442\u043e\u043b\u0431\u0446\u0430 A(1, 2) = 3.14;          \/\/ \u041f\u0440\u0438\u0441\u0432\u0430\u0438\u0432\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435<\/code><\/pre>\n<p>\u0414\u043b\u044f\u00a0\u0434\u043e\u0441\u0442\u0443\u043f\u0430 \u0431\u0435\u0437\u00a0\u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0433\u0440\u0430\u043d\u0438\u0446 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 <code>.at()<\/code>:<\/p>\n<pre><code class=\"cpp\">double value = A.at(0, 1);  \/\/ \u0427\u0443\u0442\u044c \u0431\u044b\u0441\u0442\u0440\u0435\u0435, \u043d\u043e \u0431\u0435\u0437 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u043d\u0430 \u0432\u044b\u0445\u043e\u0434 \u0437\u0430 \u043f\u0440\u0435\u0434\u0435\u043b\u044b<\/code><\/pre>\n<h4>\u041e\u0441\u043d\u043e\u0432\u043d\u044b\u0435 \u043e\u043f\u0435\u0440\u0430\u0446\u0438\u0438 \u0441 \u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438<\/h4>\n<p>\u0422\u0440\u0430\u043d\u0441\u043f\u043e\u043d\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435:<\/p>\n<pre><code class=\"cpp\">mat B = A.t();<\/code><\/pre>\n<p>\u0421\u043b\u043e\u0436\u0435\u043d\u0438\u0435\/\u0432\u044b\u0447\u0438\u0442\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446:<\/p>\n<pre><code class=\"cpp\">mat C = A + B;  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442\u044b \u0441\u043a\u043b\u0430\u0434\u044b\u0432\u0430\u044e\u0442\u0441\u044f \u043f\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043d\u043e mat D = A - B;  \/\/ \u0410\u043d\u0430\u043b\u043e\u0433\u0438\u0447\u043d\u043e<\/code><\/pre>\n<p>\u0423\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446:<\/p>\n<pre><code class=\"cpp\">mat E = A * B;  \/\/ \u041a\u043b\u0430\u0441\u0441\u0438\u0447\u0435\u0441\u043a\u043e\u0435 \u0443\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446<\/code><\/pre>\n<p>\u042d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u0440\u043d\u043e\u0435 (\u043f\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043d\u043e\u0435) \u0443\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435:<\/p>\n<pre><code class=\"cpp\">mat F = A % B;  \/\/ \u041a\u0430\u0436\u0434\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 A(i, j) \u0443\u043c\u043d\u043e\u0436\u0430\u0435\u0442\u0441\u044f \u043d\u0430 B(i, j)<\/code><\/pre>\n<p>\u0414\u0435\u043b\u0435\u043d\u0438\u0435 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432:<\/p>\n<pre><code class=\"cpp\">mat G = A \/ 2.0;  \/\/ \u0414\u0435\u043b\u0435\u043d\u0438\u0435 \u0432\u0441\u0435\u0445 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043d\u0430 2<\/code><\/pre>\n<h4>\u0412\u0435\u043a\u0442\u043e\u0440\u044b<\/h4>\n<p>\u0412\u00a0\u0432\u0435\u043a\u0442\u043e\u0440\u044b \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u044b \u043a\u043b\u0430\u0441\u0441\u0430\u043c\u0438 <code>vec<\/code> \u0438 <code>rowvec<\/code>.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u0432\u0435\u043a\u0442\u043e\u0440\u0430:<\/p>\n<pre><code class=\"cpp\">vec x = {1.0, 2.0, 3.0};  \/\/ \u0421\u0442\u043e\u043b\u0431\u0435\u0446 rowvec y = {4.0, 5.0, 6.0};  \/\/ \u0421\u0442\u0440\u043e\u043a\u0430<\/code><\/pre>\n<p>\u0414\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c:<\/p>\n<pre><code class=\"cpp\">x(1) = 42.0;  \/\/ \u041f\u0440\u0438\u0441\u0432\u0430\u0438\u0432\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0432\u0442\u043e\u0440\u043e\u043c\u0443 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0443<\/code><\/pre>\n<p><strong>\u041e\u0441\u043d\u043e\u0432\u043d\u044b\u0435 \u043e\u043f\u0435\u0440\u0430\u0446\u0438\u0438:<\/strong><\/p>\n<ul>\n<li>\n<p>\u041d\u043e\u0440\u043c\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u044f:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">x = normalise(x);  \/\/ \u041f\u0440\u0438\u0432\u043e\u0434\u0438\u0442 \u0434\u043b\u0438\u043d\u0443 \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u043a \u0435\u0434\u0438\u043d\u0438\u0446\u0435<\/code><\/pre>\n<ul>\n<li>\n<p>\u0421\u043a\u0430\u043b\u044f\u0440\u043d\u043e\u0435 \u043f\u0440\u043e\u0438\u0437\u0432\u0435\u0434\u0435\u043d\u0438\u0435:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">double result = dot(x, y);<\/code><\/pre>\n<h4>\u0420\u0430\u0431\u043e\u0442\u0430 \u0441 \u043a\u0443\u0431\u0430\u043c\u0438<\/h4>\n<p><code>cube<\/code>\u2014 \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0435 \u043c\u0430\u0441\u0441\u0438\u0432\u044b, \u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0442\u0441\u044f \u0432\u00a0\u0437\u0430\u0434\u0430\u0447\u0430\u0445 \u043e\u0431\u0440\u0430\u0431\u043e\u0442\u043a\u0438 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0439 \u0438\u043b\u0438\u00a0\u0432\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u0445 \u0440\u044f\u0434\u043e\u0432.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043a\u0443\u0431\u0430:<\/p>\n<pre><code class=\"cpp\">cube C(4, 4, 3, fill::randu);  \/\/ 4x4x3 \u043c\u0430\u0441\u0441\u0438\u0432 \u0441\u043e \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u043c\u0438 \u0447\u0438\u0441\u043b\u0430\u043c\u0438<\/code><\/pre>\n<p>\u0414\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c:<\/p>\n<pre><code class=\"cpp\">C(0, 1, 2) = 3.14;  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442 \u0432 0-\u0439 \u0441\u0442\u0440\u043e\u043a\u0435, 1-\u043c \u0441\u0442\u043e\u043b\u0431\u0446\u0435, 2-\u043c \"\u0441\u0440\u0435\u0437\u0435\"<\/code><\/pre>\n<p>\u0421\u0440\u0435\u0437\u044b \u043a\u0443\u0431\u0430:<\/p>\n<pre><code class=\"cpp\">mat slice1 = C.slice(0);  \/\/ \u0418\u0437\u0432\u043b\u0435\u043a\u0430\u0435\u043c \u043f\u0435\u0440\u0432\u044b\u0439 \"\u0441\u0440\u0435\u0437\" \u043a\u0443\u0431\u0430 (\u043c\u0430\u0442\u0440\u0438\u0446\u0430 4x4)<\/code><\/pre>\n<h4>\u0420\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b<\/h4>\n<p>\u0414\u043b\u044f\u00a0\u044d\u043a\u043e\u043d\u043e\u043c\u0438\u0438 \u043f\u0430\u043c\u044f\u0442\u0438, \u0435\u0441\u043b\u0438 \u043c\u0430\u0442\u0440\u0438\u0446\u0430 \u0441\u043e\u0434\u0435\u0440\u0436\u0438\u0442 \u043c\u043d\u043e\u0433\u043e \u043d\u0443\u043b\u0435\u0439, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 <code>sp_mat<\/code>:<\/p>\n<pre><code class=\"cpp\">sp_mat S(1000, 1000);  \/\/ \u0420\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430 1000x1000 S(0, 1) = 3.14;        \/\/ \u0414\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435<\/code><\/pre>\n<h3>Armadillo \u0432 \u0437\u0430\u0434\u0430\u0447\u0430\u0445<\/h3>\n<p><strong>\u041b\u0438\u043d\u0435\u0439\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f<\/strong><\/p>\n<p>\u041b\u0438\u043d\u0435\u0439\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f \u0440\u0435\u0448\u0430\u0435\u0442\u0441\u044f \u0447\u0435\u0440\u0435\u0437 \u043c\u0435\u0442\u043e\u0434 \u043d\u0430\u0438\u043c\u0435\u043d\u044c\u0448\u0438\u0445 \u043a\u0432\u0430\u0434\u0440\u0430\u0442\u043e\u0432:<\/p>\n<pre><code class=\"cpp\">mat X = {{1, 2}, {3, 4}, {5, 6}};  \/\/ \u041f\u0440\u0438\u0437\u043d\u0430\u043a\u0438 vec y = {7, 8, 9};                 \/\/ \u0426\u0435\u043b\u0435\u0432\u0430\u044f \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u0430\u044f  vec beta = solve(X.t() * X, X.t() * y);  \/\/ \u0420\u0435\u0448\u0435\u043d\u0438\u0435 \u0443\u0440\u0430\u0432\u043d\u0435\u043d\u0438\u044f  beta.print(\"Coefficients:\");<\/code><\/pre>\n<p><strong>K\u2011Means<\/strong><\/p>\n<p>\u041a\u043b\u0430\u0441\u0442\u0435\u0440\u0438\u0437\u0430\u0446\u0438\u044f \u0441\u00a0\u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043d\u0438\u0435\u043c K\u2011Means:<\/p>\n<pre><code class=\"cpp\">mat data(100, 2, fill::randu);  \/\/ 100 \u0442\u043e\u0447\u0435\u043a \u0432 2D \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435 mat centroids;  kmeans(centroids, data.t(), 3);  \/\/ \u0420\u0430\u0437\u0434\u0435\u043b\u044f\u0435\u043c \u043d\u0430 3 \u043a\u043b\u0430\u0441\u0442\u0435\u0440\u0430 centroids.print(\"Cluster centroids:\"); <\/code><\/pre>\n<p><strong>PCA<\/strong><\/p>\n<p>\u0414\u043b\u044f\u00a0\u0443\u043c\u0435\u043d\u044c\u0448\u0435\u043d\u0438\u044f \u0440\u0430\u0437\u043c\u0435\u0440\u043d\u043e\u0441\u0442\u0438 \u0434\u0430\u043d\u043d\u044b\u0445:<\/p>\n<pre><code class=\"cpp\">mat data(100, 5, fill::randu);  \/\/ 100 \u0442\u043e\u0447\u0435\u043a \u0432 5-\u043c\u0435\u0440\u043d\u043e\u043c \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435  vec eigval; mat eigvec;  princomp(eigvec, eigval, data);  \/\/ PCA  eigval.print(\"Eigenvalues:\"); eigvec.print(\"Eigenvectors:\");<\/code><\/pre>\n<p>\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0430\u044f \u0434\u043e\u043a\u0443\u043c\u0435\u043d\u0442\u0430\u0446\u0438\u044f <a href=\"https:\/\/arma.sourceforge.net\/docs.html\">\u0434\u043e\u0441\u0442\u0443\u043f\u043d\u0430 \u043f\u043e\u00a0\u0441\u0441\u044b\u043b\u043a\u0435<\/a>.<\/p>\n<hr\/>\n<p>\u0412 \u0437\u0430\u043a\u043b\u044e\u0447\u0435\u043d\u0438\u0435 \u0440\u0435\u043a\u043e\u043c\u0435\u043d\u0434\u0443\u044e \u043a \u043f\u043e\u0441\u0435\u0449\u0435\u043d\u0438\u044e \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b\u0435 \u043e\u0442\u043a\u0440\u044b\u0442\u044b\u0435 \u0443\u0440\u043e\u043a\u0438 \u0432 Otus:<\/p>\n<ul>\n<li>\n<p>12 \u0444\u0435\u0432\u0440\u0430\u043b\u044f: Hadoop \u0432 Docker. <a href=\"https:\/\/otus.pw\/xO92\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<li>\n<p>27 \u044f\u043d\u0432\u0430\u0440\u044f: SVD \u0438 ALS \u043d\u0430 \u0441\u043b\u0443\u0436\u0431\u0435 \u0440\u0435\u043a\u043e\u043c\u0435\u043d\u0434\u0430\u0442\u0435\u043b\u044c\u043d\u044b\u0445 \u0441\u0438\u0441\u0442\u0435\u043c. <a href=\"https:\/\/otus.pw\/4HrS\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<li>\n<p>30 \u044f\u043d\u0432\u0430\u0440\u044f: \u041e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438 \u043d\u0430 C++. <a href=\"https:\/\/otus.pw\/yxpS\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<\/ul>\n<p>\u041f\u043e\u043b\u043d\u044b\u0439 \u0441\u043f\u0438\u0441\u043e\u043a \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b\u0445 \u0443\u0440\u043e\u043a\u043e\u0432 \u043f\u043e Data Science \u0438 Machine Learning, \u0430 \u0442\u0430\u043a\u0436\u0435 \u043f\u043e \u0434\u0440\u0443\u0433\u0438\u043c IT-\u043d\u0430\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u0438\u044f\u043c \u043c\u043e\u0436\u043d\u043e \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0435\u0442\u044c <a href=\"https:\/\/otus.pw\/yhKh\/\">\u0432 \u043a\u0430\u043b\u0435\u043d\u0434\u0430\u0440\u0435<\/a>.<\/p>\n<\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!----><!----><\/div>\n<p><!----><!----><br \/> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"https:\/\/habr.com\/ru\/articles\/876190\/\"> https:\/\/habr.com\/ru\/articles\/876190\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p><em>\u041f\u0440\u0438\u0432\u0435\u0442, \u0425\u0430\u0431\u0440!<\/em><\/p>\n<p>\u0420\u0430\u0431\u043e\u0442\u0430 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438 \u0432\u00a0\u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u043e\u043c C++\u00a0\u2014 \u044d\u0442\u043e \u0431\u043e\u043b\u044c, \u0441\u0442\u0440\u0430\u0434\u0430\u043d\u0438\u0435 \u0438 \u043c\u0430\u0441\u0441\u0430 \u043d\u0435\u043d\u0443\u0436\u043d\u043e\u0433\u043e \u043a\u043e\u0434\u0430. \u041f\u043e\u0447\u0435\u043c\u0443? \u041f\u043e\u0442\u043e\u043c\u0443 \u0447\u0442\u043e\u00a0\u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u0430\u044f \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 STL \u0432\u043e\u043e\u0431\u0449\u0435 \u043d\u0435\u00a0\u043f\u0440\u0435\u0434\u043d\u0430\u0437\u043d\u0430\u0447\u0435\u043d\u0430 \u0434\u043b\u044f\u00a0\u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0430\u043b\u0433\u0435\u0431\u0440\u044b. \u041d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u0435\u0441\u043b\u0438 \u0432\u044b \u0445\u043e\u0442\u0438\u0442\u0435 \u0443\u043c\u043d\u043e\u0436\u0438\u0442\u044c \u0434\u0432\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b, \u0432\u0430\u043c \u043f\u0440\u0438\u0434\u0451\u0442\u0441\u044f \u043f\u0438\u0441\u0430\u0442\u044c \u043a\u0430\u0441\u0442\u043e\u043c\u043d\u044b\u0439 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c, \u043e\u0442\u043b\u0430\u0436\u0438\u0432\u0430\u0442\u044c \u0435\u0433\u043e, \u0430\u00a0\u0437\u0430\u0442\u0435\u043c \u0435\u0449\u0451 \u0440\u0430\u0437 \u043f\u0438\u0441\u0430\u0442\u044c \u0435\u0433\u043e, \u043a\u043e\u0433\u0434\u0430 \u043f\u043e\u0439\u043c\u0451\u0442\u0435, \u0447\u0442\u043e\u00a0\u043e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043d\u0435\u00a0\u0442\u0430.<\/p>\n<p>\u041d\u043e\u00a0\u0435\u0441\u0442\u044c \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 Armadillo \u2014 \u0438 \u0441\u00a0\u043d\u0435\u0439 \u0432\u0441\u0451 \u0438\u043d\u0430\u0447\u0435. \u0412\u043e\u0442 \u0435\u0451 \u0433\u043b\u0430\u0432\u043d\u044b\u0435 \u0444\u0438\u0447\u0438:<\/p>\n<ol>\n<li>\n<p><strong>\u041f\u0440\u043e\u0441\u0442\u043e\u0442\u0430:<\/strong> \u0432\u044b \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442\u0435 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438, \u043a\u0430\u043a\u00a0\u0441\u00a0\u043e\u0431\u044b\u0447\u043d\u044b\u043c\u0438 \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u043c\u0438. \u041d\u0438\u043a\u0430\u043a\u0438\u0445 \u00ab\u0434\u0436\u0443\u043d\u0433\u043b\u0435\u0439\u00bb \u0438\u0437\u00a0\u0446\u0438\u043a\u043b\u043e\u0432.<\/p>\n<\/li>\n<li>\n<p><strong>\u0421\u043a\u043e\u0440\u043e\u0441\u0442\u044c:<\/strong> \u043f\u043e\u0434\u0434\u0435\u0440\u0436\u043a\u0430 OpenBLAS, MKL \u0438 LAPACK \u0434\u0435\u043b\u0430\u0435\u0442 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u043e\u043b\u043d\u0438\u0435\u043d\u043e\u0441\u043d\u044b\u043c\u0438.<\/p>\n<\/li>\n<li>\n<p><strong>\u0412\u043e\u0437\u043c\u043e\u0436\u043d\u043e\u0441\u0442\u0438:<\/strong> \u0432\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u044b\u0435 \u043c\u0435\u0442\u043e\u0434\u044b \u0434\u043b\u044f\u00a0\u043a\u043b\u0430\u0441\u0442\u0435\u0440\u0438\u0437\u0430\u0446\u0438\u0438, PCA, \u0440\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u0445 \u043c\u0430\u0442\u0440\u0438\u0446 \u0438 \u0434\u0430\u0436\u0435 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u0433\u043e \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u044f.<\/p>\n<\/li>\n<li>\n<p><strong>\u0423\u043d\u0438\u0432\u0435\u0440\u0441\u0430\u043b\u044c\u043d\u043e\u0441\u0442\u044c:<\/strong> \u043c\u043e\u0436\u043d\u043e \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u043a\u0430\u043a\u00a0\u0441\u00a0\u043f\u043b\u043e\u0442\u043d\u044b\u043c\u0438, \u0442\u0430\u043a \u0438 \u0441\u00a0\u0440\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u043c\u0438 \u0434\u0430\u043d\u043d\u044b\u043c\u0438.<\/p>\n<\/li>\n<\/ol>\n<p>\u0422\u0435\u043f\u0435\u0440\u044c \u043f\u0435\u0440\u0435\u0439\u0434\u0451\u043c \u0441\u0438\u043d\u0442\u0430\u043a\u0441\u0438\u0441\u0443.<\/p>\n<h3>\u0411\u0430\u0437\u043e\u0432\u044b\u0439 \u0441\u0438\u043d\u0442\u0430\u043a\u0441\u0438\u0441 Armadillo<\/h3>\n<p>\u0414\u043b\u044f\u00a0\u0440\u0430\u0431\u043e\u0442\u044b \u0441\u00a0Armadillo \u0432\u0430\u043c \u043d\u0443\u0436\u043d\u043e \u043f\u043e\u0434\u043a\u043b\u044e\u0447\u0438\u0442\u044c \u0435\u0451 \u0437\u0430\u0433\u043e\u043b\u043e\u0432\u043e\u0447\u043d\u044b\u0439 \u0444\u0430\u0439\u043b:<\/p>\n<pre><code class=\"cpp\">#include &lt;armadillo&gt;<\/code><\/pre>\n<p>\u0427\u0442\u043e\u0431\u044b \u0438\u0437\u0431\u0435\u0436\u0430\u0442\u044c \u043f\u043e\u0441\u0442\u043e\u044f\u043d\u043d\u043e\u0433\u043e \u043d\u0430\u043f\u0438\u0441\u0430\u043d\u0438\u044f <code>arma::<\/code>, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u043e \u0438\u043c\u0451\u043d:<\/p>\n<pre><code class=\"cpp\">using namespace arma;<\/code><\/pre>\n<p>\u0412\u043e\u0442 \u0438 \u0432\u0441\u0451. \u0422\u0435\u043f\u0435\u0440\u044c \u0432\u0441\u0435 \u0433\u043e\u0442\u043e\u0432\u043e \u043a\u00a0\u0440\u0430\u0431\u043e\u0442\u0435 \u0441\u00a0\u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438, \u0432\u0435\u043a\u0442\u043e\u0440\u0430\u043c\u0438 \u0438 \u0434\u0440\u0443\u0433\u0438\u043c\u0438 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c\u0438 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438.<\/p>\n<h4>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446<\/h4>\n<p>\u0421\u0430\u043c\u044b\u0439 \u043f\u043e\u043f\u0443\u043b\u044f\u0440\u043d\u044b\u0439 \u043a\u043b\u0430\u0441\u0441 \u0432\u00a0Armadillo\u00a0\u2014 \u044d\u0442\u043e <code>mat<\/code>. \u041e\u043d \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442\u0441\u044f \u0434\u043b\u044f\u00a0\u0440\u0430\u0431\u043e\u0442\u044b \u0441\u00a0\u043f\u043b\u043e\u0442\u043d\u044b\u043c\u0438 \u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043f\u0443\u0441\u0442\u043e\u0439 \u043c\u0430\u0442\u0440\u0438\u0446\u044b:<\/p>\n<pre><code class=\"cpp\">mat A;  \/\/ \u041f\u043e\u043a\u0430 \u0431\u0435\u0437 \u0440\u0430\u0437\u043c\u0435\u0440\u043e\u0432<\/code><\/pre>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b \u0444\u0438\u043a\u0441\u0438\u0440\u043e\u0432\u0430\u043d\u043d\u043e\u0433\u043e \u0440\u0430\u0437\u043c\u0435\u0440\u0430:<\/p>\n<pre><code class=\"cpp\">mat A(3, 4);  \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 3x4, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u043c\u0443\u0441\u043e\u0440\u043e\u043c (\u0434\u043e Armadillo 10.4) \u0438\u043b\u0438 \u043d\u0443\u043b\u044f\u043c\u0438 (\u0441 \u0432\u0435\u0440\u0441\u0438\u0438 10.5).<\/code><\/pre>\n<p><strong>\u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0438\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438:<\/strong><\/p>\n<ul>\n<li>\n<p>\u0421\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0435 \u0447\u0438\u0441\u043b\u0430 (\u0440\u0430\u0432\u043d\u043e\u043c\u0435\u0440\u043d\u043e\u0435 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435):<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randu);  \/\/ \u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0430 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u043c\u0438 \u0447\u0438\u0441\u043b\u0430\u043c\u0438 \u043e\u0442 0 \u0434\u043e 1<\/code><\/pre>\n<ul>\n<li>\n<p>\u0421\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0435 \u0447\u0438\u0441\u043b\u0430 (\u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435):<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randn);  \/\/ \u0417\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u0430 \u0447\u0438\u0441\u043b\u0430\u043c\u0438 \u0441 \u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u044b\u043c \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u0435\u043c<\/code><\/pre>\n<ul>\n<li>\n<p>\u0415\u0434\u0438\u043d\u0438\u0447\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat I = eye&lt;mat&gt;(4, 4);  \/\/ \u0415\u0434\u0438\u043d\u0438\u0447\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430 4x4<\/code><\/pre>\n<ul>\n<li>\n<p>\u041c\u0430\u0442\u0440\u0438\u0446\u0430 \u0438\u0437\u00a0\u043d\u0443\u043b\u0435\u0439 \u0438\u043b\u0438\u00a0\u0435\u0434\u0438\u043d\u0438\u0446:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">mat Z = zeros&lt;mat&gt;(5, 5);  \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 5x5, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u043d\u0443\u043b\u044f\u043c\u0438 mat O = ones&lt;mat&gt;(5, 5);   \/\/ \u041c\u0430\u0442\u0440\u0438\u0446\u0430 5x5, \u0437\u0430\u043f\u043e\u043b\u043d\u0435\u043d\u043d\u0430\u044f \u0435\u0434\u0438\u043d\u0438\u0446\u0430\u043c\u0438<\/code><\/pre>\n<h4>\u0414\u043e\u0441\u0442\u0443\u043f \u043a \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c<\/h4>\n<p>\u041c\u043e\u0436\u043d\u043e \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u044c \u0434\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c \u043c\u0430\u0442\u0440\u0438\u0446\u044b \u0442\u0430\u043a\u00a0\u0436\u0435, \u043a\u0430\u043a\u00a0\u0432\u00a0Python \u0438\u043b\u0438\u00a0MATLAB:<\/p>\n<pre><code class=\"cpp\">mat A(3, 3, fill::randu); double value = A(0, 1);  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442 \u0438\u0437 \u043f\u0435\u0440\u0432\u043e\u0439 \u0441\u0442\u0440\u043e\u043a\u0438 \u0438 \u0432\u0442\u043e\u0440\u043e\u0433\u043e \u0441\u0442\u043e\u043b\u0431\u0446\u0430 A(1, 2) = 3.14;          \/\/ \u041f\u0440\u0438\u0441\u0432\u0430\u0438\u0432\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435<\/code><\/pre>\n<p>\u0414\u043b\u044f\u00a0\u0434\u043e\u0441\u0442\u0443\u043f\u0430 \u0431\u0435\u0437\u00a0\u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0433\u0440\u0430\u043d\u0438\u0446 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 <code>.at()<\/code>:<\/p>\n<pre><code class=\"cpp\">double value = A.at(0, 1);  \/\/ \u0427\u0443\u0442\u044c \u0431\u044b\u0441\u0442\u0440\u0435\u0435, \u043d\u043e \u0431\u0435\u0437 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u043d\u0430 \u0432\u044b\u0445\u043e\u0434 \u0437\u0430 \u043f\u0440\u0435\u0434\u0435\u043b\u044b<\/code><\/pre>\n<h4>\u041e\u0441\u043d\u043e\u0432\u043d\u044b\u0435 \u043e\u043f\u0435\u0440\u0430\u0446\u0438\u0438 \u0441 \u043c\u0430\u0442\u0440\u0438\u0446\u0430\u043c\u0438<\/h4>\n<p>\u0422\u0440\u0430\u043d\u0441\u043f\u043e\u043d\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435:<\/p>\n<pre><code class=\"cpp\">mat B = A.t();<\/code><\/pre>\n<p>\u0421\u043b\u043e\u0436\u0435\u043d\u0438\u0435\/\u0432\u044b\u0447\u0438\u0442\u0430\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446:<\/p>\n<pre><code class=\"cpp\">mat C = A + B;  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442\u044b \u0441\u043a\u043b\u0430\u0434\u044b\u0432\u0430\u044e\u0442\u0441\u044f \u043f\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043d\u043e mat D = A - B;  \/\/ \u0410\u043d\u0430\u043b\u043e\u0433\u0438\u0447\u043d\u043e<\/code><\/pre>\n<p>\u0423\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446:<\/p>\n<pre><code class=\"cpp\">mat E = A * B;  \/\/ \u041a\u043b\u0430\u0441\u0441\u0438\u0447\u0435\u0441\u043a\u043e\u0435 \u0443\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435 \u043c\u0430\u0442\u0440\u0438\u0446<\/code><\/pre>\n<p>\u042d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u0440\u043d\u043e\u0435 (\u043f\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043d\u043e\u0435) \u0443\u043c\u043d\u043e\u0436\u0435\u043d\u0438\u0435:<\/p>\n<pre><code class=\"cpp\">mat F = A % B;  \/\/ \u041a\u0430\u0436\u0434\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 A(i, j) \u0443\u043c\u043d\u043e\u0436\u0430\u0435\u0442\u0441\u044f \u043d\u0430 B(i, j)<\/code><\/pre>\n<p>\u0414\u0435\u043b\u0435\u043d\u0438\u0435 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432:<\/p>\n<pre><code class=\"cpp\">mat G = A \/ 2.0;  \/\/ \u0414\u0435\u043b\u0435\u043d\u0438\u0435 \u0432\u0441\u0435\u0445 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043d\u0430 2<\/code><\/pre>\n<h4>\u0412\u0435\u043a\u0442\u043e\u0440\u044b<\/h4>\n<p>\u0412\u00a0\u0432\u0435\u043a\u0442\u043e\u0440\u044b \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u044b \u043a\u043b\u0430\u0441\u0441\u0430\u043c\u0438 <code>vec<\/code> \u0438 <code>rowvec<\/code>.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u0432\u0435\u043a\u0442\u043e\u0440\u0430:<\/p>\n<pre><code class=\"cpp\">vec x = {1.0, 2.0, 3.0};  \/\/ \u0421\u0442\u043e\u043b\u0431\u0435\u0446 rowvec y = {4.0, 5.0, 6.0};  \/\/ \u0421\u0442\u0440\u043e\u043a\u0430<\/code><\/pre>\n<p>\u0414\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c:<\/p>\n<pre><code class=\"cpp\">x(1) = 42.0;  \/\/ \u041f\u0440\u0438\u0441\u0432\u0430\u0438\u0432\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0432\u0442\u043e\u0440\u043e\u043c\u0443 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0443<\/code><\/pre>\n<p><strong>\u041e\u0441\u043d\u043e\u0432\u043d\u044b\u0435 \u043e\u043f\u0435\u0440\u0430\u0446\u0438\u0438:<\/strong><\/p>\n<ul>\n<li>\n<p>\u041d\u043e\u0440\u043c\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u044f:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">x = normalise(x);  \/\/ \u041f\u0440\u0438\u0432\u043e\u0434\u0438\u0442 \u0434\u043b\u0438\u043d\u0443 \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u043a \u0435\u0434\u0438\u043d\u0438\u0446\u0435<\/code><\/pre>\n<ul>\n<li>\n<p>\u0421\u043a\u0430\u043b\u044f\u0440\u043d\u043e\u0435 \u043f\u0440\u043e\u0438\u0437\u0432\u0435\u0434\u0435\u043d\u0438\u0435:<\/p>\n<\/li>\n<\/ul>\n<pre><code class=\"cpp\">double result = dot(x, y);<\/code><\/pre>\n<h4>\u0420\u0430\u0431\u043e\u0442\u0430 \u0441 \u043a\u0443\u0431\u0430\u043c\u0438<\/h4>\n<p><code>cube<\/code>\u2014 \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0435 \u043c\u0430\u0441\u0441\u0438\u0432\u044b, \u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0442\u0441\u044f \u0432\u00a0\u0437\u0430\u0434\u0430\u0447\u0430\u0445 \u043e\u0431\u0440\u0430\u0431\u043e\u0442\u043a\u0438 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0439 \u0438\u043b\u0438\u00a0\u0432\u0440\u0435\u043c\u0435\u043d\u043d\u044b\u0445 \u0440\u044f\u0434\u043e\u0432.<\/p>\n<p>\u0421\u043e\u0437\u0434\u0430\u043d\u0438\u0435 \u043a\u0443\u0431\u0430:<\/p>\n<pre><code class=\"cpp\">cube C(4, 4, 3, fill::randu);  \/\/ 4x4x3 \u043c\u0430\u0441\u0441\u0438\u0432 \u0441\u043e \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u043c\u0438 \u0447\u0438\u0441\u043b\u0430\u043c\u0438<\/code><\/pre>\n<p>\u0414\u043e\u0441\u0442\u0443\u043f \u043a\u00a0\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u0430\u043c:<\/p>\n<pre><code class=\"cpp\">C(0, 1, 2) = 3.14;  \/\/ \u042d\u043b\u0435\u043c\u0435\u043d\u0442 \u0432 0-\u0439 \u0441\u0442\u0440\u043e\u043a\u0435, 1-\u043c \u0441\u0442\u043e\u043b\u0431\u0446\u0435, 2-\u043c \"\u0441\u0440\u0435\u0437\u0435\"<\/code><\/pre>\n<p>\u0421\u0440\u0435\u0437\u044b \u043a\u0443\u0431\u0430:<\/p>\n<pre><code class=\"cpp\">mat slice1 = C.slice(0);  \/\/ \u0418\u0437\u0432\u043b\u0435\u043a\u0430\u0435\u043c \u043f\u0435\u0440\u0432\u044b\u0439 \"\u0441\u0440\u0435\u0437\" \u043a\u0443\u0431\u0430 (\u043c\u0430\u0442\u0440\u0438\u0446\u0430 4x4)<\/code><\/pre>\n<h4>\u0420\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u044b\u0435 \u043c\u0430\u0442\u0440\u0438\u0446\u044b<\/h4>\n<p>\u0414\u043b\u044f\u00a0\u044d\u043a\u043e\u043d\u043e\u043c\u0438\u0438 \u043f\u0430\u043c\u044f\u0442\u0438, \u0435\u0441\u043b\u0438 \u043c\u0430\u0442\u0440\u0438\u0446\u0430 \u0441\u043e\u0434\u0435\u0440\u0436\u0438\u0442 \u043c\u043d\u043e\u0433\u043e \u043d\u0443\u043b\u0435\u0439, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 <code>sp_mat<\/code>:<\/p>\n<pre><code class=\"cpp\">sp_mat S(1000, 1000);  \/\/ \u0420\u0430\u0437\u0440\u0435\u0436\u0435\u043d\u043d\u0430\u044f \u043c\u0430\u0442\u0440\u0438\u0446\u0430 1000x1000 S(0, 1) = 3.14;        \/\/ \u0414\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435<\/code><\/pre>\n<h3>Armadillo \u0432 \u0437\u0430\u0434\u0430\u0447\u0430\u0445<\/h3>\n<p><strong>\u041b\u0438\u043d\u0435\u0439\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f<\/strong><\/p>\n<p>\u041b\u0438\u043d\u0435\u0439\u043d\u0430\u044f \u0440\u0435\u0433\u0440\u0435\u0441\u0441\u0438\u044f \u0440\u0435\u0448\u0430\u0435\u0442\u0441\u044f \u0447\u0435\u0440\u0435\u0437 \u043c\u0435\u0442\u043e\u0434 \u043d\u0430\u0438\u043c\u0435\u043d\u044c\u0448\u0438\u0445 \u043a\u0432\u0430\u0434\u0440\u0430\u0442\u043e\u0432:<\/p>\n<pre><code class=\"cpp\">mat X = {{1, 2}, {3, 4}, {5, 6}};  \/\/ \u041f\u0440\u0438\u0437\u043d\u0430\u043a\u0438 vec y = {7, 8, 9};                 \/\/ \u0426\u0435\u043b\u0435\u0432\u0430\u044f \u043f\u0435\u0440\u0435\u043c\u0435\u043d\u043d\u0430\u044f  vec beta = solve(X.t() * X, X.t() * y);  \/\/ \u0420\u0435\u0448\u0435\u043d\u0438\u0435 \u0443\u0440\u0430\u0432\u043d\u0435\u043d\u0438\u044f  beta.print(\"Coefficients:\");<\/code><\/pre>\n<p><strong>K\u2011Means<\/strong><\/p>\n<p>\u041a\u043b\u0430\u0441\u0442\u0435\u0440\u0438\u0437\u0430\u0446\u0438\u044f \u0441\u00a0\u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043d\u0438\u0435\u043c K\u2011Means:<\/p>\n<pre><code class=\"cpp\">mat data(100, 2, fill::randu);  \/\/ 100 \u0442\u043e\u0447\u0435\u043a \u0432 2D \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435 mat centroids;  kmeans(centroids, data.t(), 3);  \/\/ \u0420\u0430\u0437\u0434\u0435\u043b\u044f\u0435\u043c \u043d\u0430 3 \u043a\u043b\u0430\u0441\u0442\u0435\u0440\u0430 centroids.print(\"Cluster centroids:\"); <\/code><\/pre>\n<p><strong>PCA<\/strong><\/p>\n<p>\u0414\u043b\u044f\u00a0\u0443\u043c\u0435\u043d\u044c\u0448\u0435\u043d\u0438\u044f \u0440\u0430\u0437\u043c\u0435\u0440\u043d\u043e\u0441\u0442\u0438 \u0434\u0430\u043d\u043d\u044b\u0445:<\/p>\n<pre><code class=\"cpp\">mat data(100, 5, fill::randu);  \/\/ 100 \u0442\u043e\u0447\u0435\u043a \u0432 5-\u043c\u0435\u0440\u043d\u043e\u043c \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435  vec eigval; mat eigvec;  princomp(eigvec, eigval, data);  \/\/ PCA  eigval.print(\"Eigenvalues:\"); eigvec.print(\"Eigenvectors:\");<\/code><\/pre>\n<p>\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0430\u044f \u0434\u043e\u043a\u0443\u043c\u0435\u043d\u0442\u0430\u0446\u0438\u044f <a href=\"https:\/\/arma.sourceforge.net\/docs.html\">\u0434\u043e\u0441\u0442\u0443\u043f\u043d\u0430 \u043f\u043e\u00a0\u0441\u0441\u044b\u043b\u043a\u0435<\/a>.<\/p>\n<hr\/>\n<p>\u0412 \u0437\u0430\u043a\u043b\u044e\u0447\u0435\u043d\u0438\u0435 \u0440\u0435\u043a\u043e\u043c\u0435\u043d\u0434\u0443\u044e \u043a \u043f\u043e\u0441\u0435\u0449\u0435\u043d\u0438\u044e \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b\u0435 \u043e\u0442\u043a\u0440\u044b\u0442\u044b\u0435 \u0443\u0440\u043e\u043a\u0438 \u0432 Otus:<\/p>\n<ul>\n<li>\n<p>12 \u0444\u0435\u0432\u0440\u0430\u043b\u044f: Hadoop \u0432 Docker. <a href=\"https:\/\/otus.pw\/xO92\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<li>\n<p>27 \u044f\u043d\u0432\u0430\u0440\u044f: SVD \u0438 ALS \u043d\u0430 \u0441\u043b\u0443\u0436\u0431\u0435 \u0440\u0435\u043a\u043e\u043c\u0435\u043d\u0434\u0430\u0442\u0435\u043b\u044c\u043d\u044b\u0445 \u0441\u0438\u0441\u0442\u0435\u043c. <a href=\"https:\/\/otus.pw\/4HrS\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<li>\n<p>30 \u044f\u043d\u0432\u0430\u0440\u044f: \u041e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438 \u043d\u0430 C++. <a href=\"https:\/\/otus.pw\/yxpS\/\">\u041f\u043e\u0434\u0440\u043e\u0431\u043d\u0435\u0435<\/a><\/p>\n<\/li>\n<\/ul>\n<p>\u041f\u043e\u043b\u043d\u044b\u0439 \u0441\u043f\u0438\u0441\u043e\u043a \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b\u0445 \u0443\u0440\u043e\u043a\u043e\u0432 \u043f\u043e Data Science \u0438 Machine Learning, \u0430 \u0442\u0430\u043a\u0436\u0435 \u043f\u043e \u0434\u0440\u0443\u0433\u0438\u043c IT-\u043d\u0430\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u0438\u044f\u043c \u043c\u043e\u0436\u043d\u043e \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0435\u0442\u044c <a href=\"https:\/\/otus.pw\/yhKh\/\">\u0432 \u043a\u0430\u043b\u0435\u043d\u0434\u0430\u0440\u0435<\/a>.<\/p>\n<\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!----><!----><\/div>\n<p><!----><!----><br \/> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"https:\/\/habr.com\/ru\/articles\/876190\/\"> https:\/\/habr.com\/ru\/articles\/876190\/<\/a><br \/><\/br><\/br><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-446015","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/446015","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=446015"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/446015\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=446015"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=446015"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=446015"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}