1. 标准数组:全1数组,全0数组,单位矩阵,随机矩阵,对角矩阵以及元素为指定常数的数组。
2.全1数组用ones函数,全0数组用zeros函数。
对于ones和zeros函数,当只有一个输入参数时,即ones(n)或zeros(n),Matlab就分别生成一个n×n的全1或者全0数组。当有两个输入参数时,即ones (r,c)或者zeros(r,c),Matlab就分别生成r行c列的全1或者全0数组。要想生成一个与其他数组相同维数的全1或者全0数组,用户只要在ones或者zeros的参数中调用size函数就可以了。测试数组:ones(4),m = ones(4,8)zeros(4),zeros(3,5),size(m),zeros(size(m))。
3.单位矩阵用eye函数。该函数用与ones和zeros函数相同的语法格式来生成单位矩阵。单位矩阵或数组是具有如下取值的矩阵或数组:除A(i,i)之外,所有其他元素都为0,其中i=min(r,c),min(r,c)是矩阵A中的行数和列数的最小数。
4.随机矩阵用rand函数。函数rand生成均匀分布的随机数组,其元素取值介于0-1之间。直接调用rand产生一个随机数,随机数组用rand(n)。另外randn 函数将生成均值为0,方差为1的正态分布矩阵。rand和randn用法和ones相同。
5.对角矩阵用diag函数。在该数组中,一个向量可以被放在与数组的主对角线平行的任何位置。验证:a = 1:5 diag(a) diag(a,1) diag(a,-2)
6.几种生成所有元素都相同的数组的方法,先令d=pi
(1)d*one(3,4) slowest method
(2)d+zeros(3,4) slower method
(3)d(ones(3,4)) fast method
(4)repmat(d,3,4) fastest method
数组的数据量较小时,4种方法都可以。随着数组维数的增大,含有标量乘法的方法d*one(3,4)就会使矩阵生成过程变慢。因为加法通常都比乘法运算速度快,较好的办法就是将用到的标量加到加到一个全0数组zeros(3,4)上。
后两种方法虽不直观,它们却是生成大数组的最快方法,因为都用到了数组索引。(3)d(ones(r,c))先生成一个r×c的全1数组,然后用这个数组来索引和复制标量d。尽管这种方法没有用到浮点运算,但生成一个全1临时数组将会占用内存,并且消耗时间,因此使得这种方法的速度变慢。
方法repmat(d,r,c)调用函数repmat,即replicate matrix(复制矩阵)的缩写。对于标量,该函数执行如下操作步骤:
D(r*c) = d;-->>创建数组D,其第r*c个数为d,前面的数都为0
D( : )=d;
D=reshape(D,r,c);
上面的代码首先采用标量扩展方法生成一个具有r×c个元素的向量(特别注意:生成的为行向量而非矩阵)。然后将数组所有值赋值为d,再用函数reshape将这个向量变换成一个r×c的数组。