MATLAB之特殊矩阵

通用的特殊矩阵

zeros函数 产生全0函数，即零矩阵
ones函数 产生全1函数，即幺矩阵
eye函数 产生对角线为1的函数。当矩阵是方阵时，得到一个单位矩阵
rand函数 产生（0，1）区间均匀分布的随机矩阵
randn函数 产生均值为0，方差为1的标准正态分布随机矩阵X μ+σX可以得到均值为μ方差为σ方的随机数

zeros(m)
zeros(m,n)
zeros(size(A))  %产生与A矩阵同型的零矩阵

例子5阶两位随机整数矩阵A，在产生均值为0.6，方差为0.1的5阶正态分布随机矩阵B，最后验证（A+B）I=IA+BI I为单位矩阵

A=rand(5);
A=fix(10+(99-10+1).*A)
B=randn(5);
B=0.6+sqrt(0.1).*B
I=eye(5);
(A+B)*I
I*A+B*I

用于专门学科特殊矩阵

魔方矩阵

>> M=magic(3)

M =

     8     1     6
     3     5     7
     4     9     2

n阶魔方矩阵由1，2，3，。。。，n方个整数组成，且每行每列以及主副对角线上各点n个元素之和都相等
和为（1+2+3+……+n方）/n=（n+n立方）/2

范德蒙矩阵 常用在各种通讯系统的纠错编码中，例如常用的Reed-Solomon

vander（V）为以向量v为基础的范德蒙矩阵

>> vander(1:5)

ans =

     1     1     1     1     1
    16     8     4     2     1
    81    27     9     3     1
   256    64    16     4     1
   625   125    25     5     1

希尔伯特矩阵 hilb（n）

   >> hilb(4)

ans =

    1.0000    0.5000    0.3333    0.2500
    0.5000    0.3333    0.2500    0.2000
    0.3333    0.2500    0.2000    0.1667
    0.2500    0.2000    0.1667    0.1429

伴随矩阵

compan§其中p是一个多项式的系数向量，高次幂系数排在前，低次幂系数排在后。
例如，生成一个x3-2x2-5x+6的伴随矩阵

>> p=[1,-2,-5,6];A=compan(p)

A =

     2     5    -6
     1     0     0
     0     1     0

帕斯卡矩阵

二项式定理，(x+y)^n展开后的系数随着n的增大组成一个三角形表，杨辉三角，
把二项式系数，依次填在左对角线上，然后提取左侧的n行n列元素即为n阶帕斯卡矩阵
P(i,j)=P(i,j-1)+P(i-1,j),且P(i,1)=1,P(1,j)=1
生成5阶帕斯卡矩阵，验证他的逆矩阵的所有元素也为整数

>> P=pascal(5)

P =

     1     1     1     1     1
     1     2     3     4     5
     1     3     6    10    15
     1     4    10    20    35
     1     5    15    35    70

>> inv(P)

ans =

    5.0000  -10.0000   10.0000   -5.0000    1.0000
  -10.0000   30.0000  -35.0000   19.0000   -4.0000
   10.0000  -35.0000   46.0000  -27.0000    6.0000
   -5.0000   19.0000  -27.0000   17.0000   -4.0000
    1.0000   -4.0000    6.0000   -4.0000    1.0000