幻灯片 1

Size: px

Start display at page:

Download "幻灯片 1"

禁墀窦
7 years ago
Views:

1 第 3 章 MATLAB 矩阵分析与处理 3.1 特殊矩阵 3.2 矩阵结构变换 3.3 矩阵求逆与线性方程组求解 3.4 矩阵求值 3.5 矩阵的特征值与特征向量 3.6 矩阵的超越函数

2 特殊矩阵通用的特殊矩阵常用的产生通用特殊矩阵的函数有 : zeros: 产生全 0 矩阵 ( 零矩阵 ) ones: 产生全 1 矩阵 ( 幺矩阵 ) eye: 产生单位矩阵 rand: 产生 0~1 间均匀分布的随机矩阵 randn: 产生均值为 0, 方差为 1 的标准正态分布随机矩阵

3 例 3.1 分别建立和与矩阵 A 同样大小的零矩阵 (1) 建立一个 3 3 零矩阵 zeros(3) (2) 建立一个 3 2 零矩阵 zeros(3,2) (3) 设 A 为 2 3 矩阵, 则可以用 zeros(size(a)) 建立一个与矩阵 A 同样大小零矩阵 A=[1 2 3;4 5 6]; % 产生一个 2 3 阶矩阵 A zeros(size(a)) % 产生一个与矩阵 A 同样大小的零矩阵

4 例 3.2 建立随机矩阵 : (1) 在区间 [20,50] 内均匀分布的 5 阶随机矩阵 (2) 均值为 0.6 方差为 0.1 的 5 阶正态分布随机矩阵命令如下 : x=20+(50-20)*rand(5) y=0.6+sqrt(0.1)*randn(5) 此外, 常用的函数还有 reshape(a,m,n), 它在矩阵总元素保持不变的前提下, 将矩阵 A 重新排成 m n 的二维矩阵

5 用于专门学科的特殊矩阵 (1) 魔方矩阵魔方矩阵有一个有趣的性质, 其每行每列及两条对角线上的元素和都相等对于 n 阶魔方阵, 其元素由 1,2,3,,n2 共 n2 个整数组成 MATLAB 提供了求魔方矩阵的函数 magic(n), 其功能是生成一个 n 阶魔方阵

6 例 3.3 将 101~125 等 25 个数填入一个 5 行 5 列的表格中, 使其每行每列及对角线的和均为 565 M=100+magic(5) %5 阶魔方矩阵每行每列及对角线的和均为 65

7 范得蒙矩阵范得蒙 (Vandermonde) 矩阵最后一列全为 1, 倒数第二列为一个指定的向量, 其他各列是其后列与倒数第二列的点乘积可以用一个指定向量生成一个范得蒙矩阵在 MATLAB 中, 函数 vander(v) 生成以向量 V 为基础向量的范得蒙矩阵例如,A=vander([1;2;3;5]) 即可得到上述范得蒙矩阵 A=[ ; ; ; ]

8 希尔伯特矩阵在 MATLAB 中, 生成希尔伯特矩阵的函数是 hilb(n),h ij =1/(i+j-1) 使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果 MATLAB 中, 有一个专门求希尔伯特矩阵的逆的函数 invhilb(n), 其功能是求 n 阶的希尔伯特矩阵的逆矩阵

9 例 3.4 求 4 阶希尔伯特矩阵及其逆矩阵命令如下 : format rat % 以有理形式输出 H=hilb(4) H=invhilb(4)

10 托普利兹矩阵托普利兹 (Toeplitz) 矩阵除第一行第一列外, 其他每个元素都与左上角的元素相同生成托普利兹矩阵的函数是 toeplitz(x,y), 它生成一个以 x 为第一列,y 为第一行的托普利兹矩阵这里 x, y 均为向量, 两者不必等长 toeplitz(x) 用向量 x 生成一个对称的托普利兹矩阵例如 T=toeplitz(1:6)

11 伴随矩阵 MATLAB 生成伴随矩阵的函数是 compan(p), 其中 p 是一个多项式的系数向量, 高次幂系数排在前, 低次幂排在后例如, 为了求多项式的 x3-7x+6 的伴随矩阵, 可使用命令 : p=[1,0,-7,6]; compan(p) ans=

12 帕斯卡矩阵二次项 (x+y) n 展开后的系数随 n 的增大组成一个三角形表, 称为杨辉三角形由杨辉三角形表组成的矩阵称为帕斯卡 (Pascal) 矩阵函数 pascal(n) 生成一个 n 阶帕斯卡矩阵

13 例 3.5 求 (x+y) 5 的展开式在 MATLAB 命令窗口, 输入命令 : pascal(6) 矩阵次对角线上的元素 1,5,10,10,5,1 即为展开式的系数

14 矩阵结构调整变换对角阵 1. 对角阵只有对角线上有非 0 元素的矩阵称为对角矩阵, 对角线上的元素相等的对角矩阵称为数量矩阵, 对角线上的元素都为 1 的对角矩阵称为单位矩阵

15 提取矩阵的对角线元素设 A 为 m n 矩阵,diag(A) 函数用于提取矩阵 A 主对角线元素, 产生一个具有 min(m,n) 个元素的列向量 diag(a) 函数还有一种形式 diag(a,k), 其功能是提取第 k 条对角线的元素

16 >>A = [1 2 3;4 5 6; 7 8 9] ; >>D=diag(A) D = >> diag(a,1) ans = 2 6

17 构造对角矩阵设 V 为具有 m 个元素的向量,diag(V) 将产生一个 m m 对角矩阵, 其主对角线元素即为向量 V 的元素 diag(v) 函数也有另一种形式 diag(v,k), 其功能是产生一个 n n(n=m+ k ) 对角阵, 其第 k 条对角线的元素即为向量 V 的元素 >>diag(d) ans =

18 例 3.6 先建立 5 5 矩阵 A, 然后将 A 的第一行元素乘以 1, 第二行乘以 2,, 第五行乘以 5 A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12, 19,21,3;... 11,18,25,2,19]; D=diag(1:5); D*A % 用 D 左乘 A, 对 A 的每行乘以一个指定常数

19 三角阵三角阵又进一步分为上三角阵和下三角阵, 所谓上三角阵, 即矩阵的对角线以下的元素全为 0 的一种矩阵, 而下三角阵则是对角线以上的元素全为 0 的一种矩阵

20 上三角矩阵求矩阵 A 的上三角阵的 MATLAB 函数是 triu(a) triu(a) 函数也有另一种形式 triu(a,k), 其功能是求矩阵 A 的第 k 条对角线以上的元素例如, 提取矩阵 A 的第 2 条对角线以上的元素, 形成新的矩阵 B triu(ones(4,4),-1) ans =

21 下三角矩阵在 MATLAB 中, 提取矩阵 A 的下三角矩阵的函数是 tril(a) 和 tril(a,k), 其用法与提取上三角矩阵的函数 triu(a) 和 triu(a,k) 完全相同 tril(ones(4,4),-1) ans =

22 矩阵的转置转置运算符是单撇号 ( ) >> A = [1 2 3;4 5 6; 7 8 9] A = >> A' ans =

23 矩阵的旋转利用函数 rot90(a,k) 将矩阵 A 按逆时针旋转 90º 的 k 倍, 当 k 为 1 时可省略 X = Y = rot90(x) Y =

24 矩阵的左右翻转对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换, 第二列和倒数第二列调换,, 依次类推 MATLAB 对矩阵 A 实施左右翻转的函数是 fliplr(a) If A is the 3-by-2 matrix, A = then fliplr(a) produces If A is a row vector, A = then fliplr(a) produces

25 矩阵的上下翻转 MATLAB 对矩阵 A 实施上下翻转的函数是 flipud(a) If A is the 3-by-2 matrix, A = then flipud(a) produces If A is a column vector, A = then flipud(a) produces A = 7 5 3

26 矩阵的逆与伪逆对于一个方阵 A, 如果存在一个与其同阶的方阵 B, 使得 :A B=B A=I (I 为单位矩阵 ) 则称 B 为 A 的逆矩阵, 当然,A 也是 B 的逆矩阵求一个矩阵的逆是一件非常烦琐的工作, 容易出错, 但在 MATLAB 中, 求一个矩阵的逆非常容易求方阵 A 的逆矩阵可调用函数 inv(a) 例 3.7 求方阵 A 的逆矩阵, 且验证 A 与 A^-1 是否是互逆. A=[1,-1,1;5,-4,3;2,1,1]; B=inv(A); A*B B*A

27 矩阵的逆与伪逆如果矩阵 A 不是一个方阵, 或者 A 是一个非满秩的方阵时, 矩阵 A 没有逆矩阵, 但可以找到一个与 A 的转置矩阵 A 同型的矩阵 B, 使得 : A B A=A B A B=B 此时称矩阵 B 为矩阵 A 的伪逆, 也称为广义逆矩阵在 MATLAB 中, 求一个矩阵伪逆的函数是 pinv(a) >> A=[ ; ; ]; >> B=pinv(A) B =

28 用矩阵求逆方法求解线性方程组在线性方程组 Ax=b 两边各左乘 A -1, 有由于 A -1 A=I, 故得例 3.8 用求逆矩阵的方法解线性方程组命令如下 : A=[1,2,3;1,4,9;1,8,27]; b=[5,-2,6]'; x=inv(a)*b 也可以运用左除运算符 \ 求解线性代数方程组 A=[1,2,3;1,4,9;1,8,27]; b=[5,-2,6]'; x=a\b x 2y 3z 5 x 4y 9z 2 x 8y 27z 6 A -1 Ax=A -1 b x=a -1 b Matlab 与数学实验

29 矩阵求值方阵的行列式把一个方阵看作一个行列式, 并对其按行列式的规则求值, 这个值就称为所对应的行列式的值在 MATLAB 中, 求方阵 A 所对应的行列式的值的函数是 det(a) A=rand(5) B=det(A)

30 矩阵的秩与迹矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩在 MATLAB 中, 求矩阵秩的函数是 rank(a) 矩阵的迹矩阵的迹等于矩阵的对角线元素之和, 也等于矩阵的特征值之和在 MATLAB 中, 求矩阵的迹的函数是 trace(a)

31 向量和矩阵的范数矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度范数有多种方法定义, 其定义不同, 范数值也就不同向量的 3 种常用范数及其计算函数在 MATLAB 中, 求向量范数的函数为 : (1) norm(v) 或 norm(v,2): 计算向量 V 的 2 范数 (2) norm(v,1): 计算向量 V 的 1 范数 (3) norm(v,inf): 计算向量 V 的范数

32 向量和矩阵的范数矩阵的范数及其计算函数 MATLAB 提供了求 3 种矩阵范数的函数, 其函数调用格式与求向量的范数的函数完全相同矩阵的条件数在 MATLAB 中, 计算矩阵 A 的 3 种条件数的函数是 : (1) cond(a,1) 计算 A 的 1 范数下的条件数 (2) cond(a) 或 cond(a,2) 计算 A 的 2 范数数下的条件数 (3) cond(a,inf) 计算 A 的范数下的条件数

33 矩阵的特征值与特征向量在 MATLAB 中, 计算矩阵 A 的特征值和特征向量的函数是 eig(a), 常用的调用格式有 3 种 : (1) E=eig(A): 求矩阵 A 的全部特征值, 构成向量 E (2) [V,D]=eig(A): 求矩阵 A 的全部特征值, 构成对角阵 D, 并求 A 的特征向量构成 V 的列向量 (3) [V,D]=eig(A, nobalance ): 与第 2 种格式类似, 但第 2 种格式中先对 A 作相似变换后求矩阵 A 的特征值和特征向量, 而格式 3 直接求矩阵 A 的特征值和特征向量

34 矩阵的特征值与特征向量例 3.9 用求特征值的方法解方程 3x 5-7x 4 +5x 2 +2x-18=0 p=[3,-7,0,5,2,-18]; A=compan(p); %A 的伴随矩阵 x1=eig(a) % 求 A 的特征值 x2=roots(p) % 直接求多项式 p 的零点

35 矩阵的超越函数矩阵平方根 sqrtm sqrtm(a) 计算矩阵 A 的平方根矩阵对数 logm logm(a) 计算矩阵 A 的自然对数此函数输入参数的条件与输出结果间的关系和函数 sqrtm(a) 完全一样

36 矩阵的超越函数矩阵指数 expm expm1 expm2 expm3 expm(a) expm1(a) expm2(a) expm3(a) 的功能都求矩阵指数 e A 普通矩阵函数 funm funm(a, fun ) 用来计算直接作用于矩阵 A 的由 fun 指定的超越函数值当 fun 取 sqrt 时, funm(a, sqrt ) 可以计算矩阵 A 的平方根, 与 sqrtm(a) 的计算结果一样

Matlab和应用数学2-2.ppt [兼容模式]

Matlab和应用数学2-2.ppt [兼容模式] MATLAB 矩阵运算上课日期 :2007.7.12 2007.7.18 教学方法 : 授课上机练习和自学相结合授课教师 : 梁克维 (matlkw@zju.edu.cn) 上上课 : 上午 8:30 11:30, 西 1-204 机 : 下午 1:00-4:00, 计算机中心 MATLAB 矩阵矩阵的建立 1. 直