第3部分MATLAB矩阵分析与处理.pptVIP

第3章 MATLAB矩阵分析与处理 第3章 MATLAB矩阵分析与处理 特殊矩阵 矩阵结构变换 矩阵求逆与线性方程组求解 矩阵求值 3.1 特殊矩阵 3.1.1 通用的特殊矩阵 常用的产生通用特殊矩阵的函数有： zeros：产生全0矩阵(零矩阵)。 ones：产生全1矩阵(幺矩阵)。 eye：产生单位矩阵。 rand：产生0～1间均匀分布的随机矩阵。 randn：产生均值为0，方差为1的标准正态分布随机 矩阵。 这几个函数的调用格式相似，下面以产生零矩 阵的zeros函数为例进行说明。其调用格式为： zeros(m): 产生m×m零矩阵 zeros(m,n) 产生m×n零矩阵 zeros(size(A)) 产生与矩阵A同样大小的零矩阵。 例3.1 分别建立3×3、3×2和与矩阵A同样大小的零矩阵 (1) 建立一个3×3零矩阵。zeros(3)ans= 0 0 0 0 0 0 0 0 0 (2) 建立一个3×2零矩阵。 zeros(3,2) ans= 0 0 0 0 0 0 (3) 设A为2×3矩阵，则可以用zeros(size(A))建立一个与 矩阵A同样大小零矩阵。A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵Azeros(size(A)) %产生一个与矩阵A同样大小的零矩阵 ans= 0 0 0 0 0 0 例3.2 建立随机矩阵： (1) 在区间[20,50]内均匀分布的5阶随机矩阵。 (2) 均值为0.6、方差为0.1的5阶正态分布随机矩阵。 rand：产生0～1间均匀分布的随机矩阵。 要得到[a,b]区间上均匀分布的随机数，需用yi=a+(b-a)xi randn：产生均值为0，方差为1的标准正态分布随机矩阵。 命令如下：x=20+(50-20)*rand(5)x = 48.5039 42.8629 38.4630 32.1712 21.7367 26.9342 33.6940 43.7581 48.0641 30.5860 38.2053 20.5551 47.6544 47.5071 44.3950 34.5795 44.6422 42.1462 32.3081 20.2958 46.7390 33.3411 25.2880 46.8095 24.1667 y=0.6+sqrt(0.1)*randn(5) y = 0.8713 0.4735 0.8114 0.0927 0.7672 0.9966 0.8182 0.9766 0.6814 0.6694 0.0960 0.8579 0.2197 0.2659 0.3085 0.1443 0.8251 0.5937 1.0475 -0.0864 0.7806 1.0080 0.5504 0.3454 0.5813 3.1.2 用于专门学科的特殊矩阵 (1) 魔方矩阵 魔方矩阵有一个有趣的性质，其每行、每列及 两条对角线上的元素和都相等。对于n阶魔方阵，其 元素由1,2,3,…,n2共n2个整数组成。MATLAB提供 了求魔方矩阵的函数magic(n)，其功能是生成一个n 阶魔方阵。 magic(3) ans = 8 1 6 3 5 7 4 9 2 例3.3 将101-125等25个数填入一个5行5列的表格 中，使其每行每列及对角线的和均为565。 一个5姐魔方矩阵的每行、每列及对角线的和 均为65，对其每个元素都加100后，这些和变成565. (2) 范得蒙德矩阵 范得蒙德(Vandermonde)矩阵最后一列全为1， 倒数第二列为一个指定的向量，其他各列是其后列 与倒数第二列的点乘积。可以用一个指定向量生成 一个范得蒙矩阵。在MATLAB中，函数vander(V)生 成以向量V为基础向量的范得蒙矩阵。例如， A=vander([1;2;3;5])即可得到上述范得蒙矩阵。 例3.4 求4阶希尔伯特矩阵及其逆矩阵。 命令如下： (4) 托普利兹矩阵 托普利兹(Toeplitz)矩阵