2020-11-03 MATLAB学习小结（四）

特殊矩阵

通用性的特殊矩阵

常用的函数：

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

zeros函数的调用格式：
zeros(m)：产生mm零矩阵。
zeros(m,n)：产生mn零矩阵。
zeros(size(A))：产生与矩阵A同样大小的零矩阵。
例如：

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

rand函数：产生（0，1）开区间均匀分布的随机数x。
fix(a+(b-a+1)*x)：产生[a,b]区间上均匀分布的随机整数。
randn函数：产生均值为0、方差为1的标准正态分布随机x。
μ+σx：得到均值为μ、方差为σ²的随机数。

用于专门学科的特殊矩阵

1.魔方矩阵——Magic Square
例如：

解释：
（1）n阶魔方阵由1,2,3,…，n²共n²个整数组成，且每行、每列以及主、副对角线上各n个元素之和都相等。
（2）n阶魔方阵每行每列元素的和为（1+2+3+…+n²）/n = (n+n³)/2 。
（3）n>2时有很多不同的n阶魔方阵，MATLAB函数magic(n)产生一个特定的魔方阵。

例2 产生8阶魔方阵，求其每行每列元素的和。

2.范德蒙矩阵
对于向量v=[v1,v2,…,vn]，范德蒙矩阵的一般形式为：

范德蒙（Vandermonde）矩阵是法国数学家范德蒙提出的一种特殊矩阵。范得蒙矩阵的最后一列全为1，即向量v各元素的零次方，倒数第二列为指定的向量v，即向量v各元素的一次方， 其他各列是其后列与倒数第二列的点乘积。

在MATLAB中，函数vander(V)生成以向量V为基础的范得蒙矩阵。

范德蒙矩阵常用在各种通信系统的纠错编码中，例如，常用的Reed-Solomon编码即以范德蒙矩阵为基础。

3.希尔伯特矩阵
n阶希尔伯特（Hilbert）矩阵的一般形式为：

希尔伯特矩阵的元素为H(i,j)=1/(i+j-1)。

在MATLAB中，生成n阶希尔伯特矩阵的函数是hilb(n)。

format rat 是先将数字有理化。

希尔伯特矩阵是著名的病态矩阵，即任何一个元素发生较小的变动，整个矩阵的值和逆矩阵都会发生很大变化。病态程度和矩阵的阶数相关，随着阶数的增加病态越严重。

4.伴随矩阵
设多项式p(x)为anxⁿ+an-1xⁿ﹣¹+…+a1x+a0，则多项式的伴随矩阵是：

p(x)称为A的特征多项式，方程p(x)=0的根称为A的特征值。

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

可以求出伴随矩阵的特征值，该特征值等于多项式方程的根。

5.帕斯卡矩阵
根据二项式定理，(x+y)ⁿ展开后的系数随着n的增大组成一个三角形表，这个三角形称为杨辉三角形。
把二项式系数依次填写在矩阵的左侧对角线上，然后提取左侧的n行n列元素即为n阶帕斯卡（Pascal）矩阵。

帕斯卡矩阵的第一行元素和第一列元素都为1，其余位置的元素是该元素的左边元素与上面元素相加，P(i,j)=P(i,j-1)+P(i-1,j)，且P(i,1)=1,P(1,j)=1。
函数pascal(n)生成一个n阶帕斯卡矩阵。

例3 生成5阶帕斯卡矩阵，验证它的逆矩阵的所有元素也为整数。