RcppArmadillo幾乎和C++中的Armadillo一樣，因此本文主要參考Armadillo主頁上的手冊(cè)。

http://arma.sourceforge.net/docs.html

1. 常用數(shù)據(jù)類型

Mat<type>為模板類，其中type可以是：float, double, std::complex, std::complex, short, int, long, and unsigned versions of short, int, long等。為方便起見，Armadillo C++已經(jīng)預(yù)定義了以下類型。

在Armadillo中，矩陣是按照一列一列(column by column)存在內(nèi)存中的(column-major ordering)。

2. 數(shù)學(xué)運(yùn)算

3. 向量、矩陣和域的創(chuàng)建

基本創(chuàng)建

//Matrix: X
mat X(n_rows, n_cols); 
mat X(n_rows, n_cols, fill_type); 
mat X(size(Y)); 
mat X(size(Y), fill_type);
mat X("1 0;2 0"); \\create a matrix from a string. (first row is "1 0", second row is "2 0")
cx_mat X(mat,mat);   \\for constructing a complex matrix out of two real matrices

//Vector: X
vec X(n_elem);
vec X(n_elem, fill_type);
vec X(size(Y)); 
vec X(size(Y), fill_type); 
cx_vec X(vec,vec);   \\for constructing a complex vector out of two real vectors

//Cube or 3-D matrix: X
cube X(n_rows, n_cols, n_slices); 
cube X(n_rows, n_cols, n_slices, fill_type); 
cube X(size(Y));
cube X(size(Y), fill_type); 
cx_cube X(cube, cube);   \\for constructing a complex cube out of two real cubes

//Field: F
field<object_type> F(n_elem) 
field<object_type> F(n_rows, n_cols) 
field<object_type> F(n_rows, n_cols, n_slices) 
field<object_type> F(size(X)) 
//似乎object_type只能是mat或者cube,且不能存儲(chǔ)復(fù)數(shù)類型，比如cx_mat.

//An example of field
mat A(3,4,fill::randu);
vec B(5,fill::randn);
field <mat> F(2,1);
F(0)=A;
F(1)=B;

其中fill_type是可選的，可以是以下選擇

用函數(shù)創(chuàng)建

注釋：

rand_type可以是randi()、randu()、randn()和randg()，分別代表 [a,b] 區(qū)間中的整數(shù)隨機(jī)值， U[0,1] 分布中的隨機(jī)浮點(diǎn)值，從標(biāo)準(zhǔn)正態(tài)分布中抽取的隨機(jī)值，從參數(shù)為a,b的Gamma分布中抽取的隨機(jī)值。distr_param(a,b) 只適用于randi()和randg()。
e.g.

vec v=randu<vec>(5);

regspace()中delta默認(rèn)是1或-1。

4. 初始化，元素訪問，屬性和成員函數(shù)

4.1. 元素初始化 Element initialization

// C++11
vec v = { 1, 2, 3 };
mat A = { {1, 3, 5},
          {2, 4, 6} };

4.2. 元素訪問 Element access

無論是向量vec，矩陣mat，立方體cube，還是域field，每個(gè)維數(shù)均是從0開始的。

4.3. 子矩陣訪問 Submatrix view

矩陣X的連續(xù)子集訪問

注釋:
(1) span(start,end)可以被span::all代替，意味著這一維上所有的元素。
(2) X.diag(k)可以改變第k個(gè)對(duì)角線的值。

mat X=randn<mat>(4,4);
vec v={1,2,3,4};
X.diag()=v;

向量V的連續(xù)子集訪問

向量或矩陣X的間斷子集訪問

立方體（三維矩陣）Q 的切片 slice

域F的子集訪問

4.4. 屬性 Attribute

注釋:

• 返回值是無符號(hào)整數(shù)(uword)
• 返回值是read-only的；如果要改變大?。ňS數(shù)），用成員函數(shù).set_size(), .copy_size(), .zeros(), .ones(), 或者.reset()。

.set_size()

.set_size( n_elem )
.set_size( n_rows, n_cols )
.set_size( n_rows, n_cols, n_slices )
.set_size( size(X) )

.copy_size(A)
把維數(shù)設(shè)置成和A一樣。

.zeros()

.zeros( n_elem )
.zeros( n_rows, n_cols )
.zeros( n_rows, n_cols, n_slices )  
.zeros( size(X) )   

.ones()
參見.zeros()。

.reset()
把維數(shù)設(shè)置成0，意味著無元素。

4.5. 其他成員函數(shù) Other member function

注釋：

.transform(lambda_function)

// C++11 only example
mat A = ones<mat>(4,5);
// add 123 to every element
A.transform( [](double val) { return (val + 123.0); } );

.reshape()和.resize()的區(qū)別在于前者不會(huì)保存原對(duì)象的布局,而后者會(huì)保留原對(duì)象的布局，且后者更快。例如，如果新對(duì)象的維數(shù)大于原對(duì)象的維數(shù)，則新對(duì)象中原維數(shù)外的元素會(huì)被設(shè)置成0。
e.g.

mat A = randu<mat>(2,3);
A.reshape(4,4);
           [,1]      [,2] [,3] [,4]
[1,] 0.02567623 0.8880936    0    0
[2,] 0.12546129 0.6520889    0    0
[3,] 0.52724939 0.0000000    0    0
[4,] 0.30407942 0.0000000    0    0

mat A = randu<mat>(2,3);
A.resize(4,4);
          [,1]      [,2]      [,3] [,4]
[1,] 0.5451790 0.2632051 0.6375933    0
[2,] 0.3753245 0.8050394 0.1627499    0
[3,] 0.0000000 0.0000000 0.0000000    0
[4,] 0.0000000 0.0000000 0.0000000    0

.memptr()
可被用于和一些庫交互，比如FFTW。

      mat A = randu<mat>(5,5);
const mat B = randu<mat>(5,5);
      double* A_mem = A.memptr();
const double* B_mem = B.memptr();

5. 常用函數(shù)

5.1. 向量、矩陣和立方體的一般函數(shù)

5.2. 其他一些數(shù)學(xué)函數(shù)

miscellaneous element-wise functions:

三角函數(shù) Trigonometric element-wise functions (cos, sin, tan, …)

cos, acos, cosh, acosh

sin, asin, sinh, asinh

tan, atan, tanh, atanh

atan2, hypot

5.3. 矩陣的分解、因子化、逆矩陣和線性方程的解

Dense Matrix

Sparse Matrix

5.4. 信號(hào)和圖像處理

5.5. 統(tǒng)計(jì)和聚類

以上就是R語言科學(xué)計(jì)算RcppArmadillo簡(jiǎn)明手冊(cè)的詳細(xì)內(nèi)容，更多關(guān)于RcppArmadillo簡(jiǎn)明手冊(cè)的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章！

最新評(píng)論