第五讲线性代数中的数值计算问题.pptVIP

第五讲 线性代数中的 数值计算问题【引 例 】求下列三阶线性代数方程组的近似解 2 x5 x4 x ?51 2 3x5 x2 x ?61 2 3? x2 x4 x ?5 1 2 3MATLAB程序为: A[2 -5 4;1 5 -2;-1 2 4]; b[5;6;5]; xA\b在MATLAB命令窗口,先输入下列命令构造 系数矩阵A和右端向量b: A[2 -5 4;1 5 -2;-1 2 4] A 2 -5 41 5 -2 -1 2 4 b[5;6;5] b 565 然后只需输入命令xA\b即可求得解x: xA\b x 2.7674 1.1860 1.3488一、 特殊矩阵的实现一、 特殊矩阵的实现 常见的特殊矩阵有零矩阵、幺矩阵、单位矩阵、三角 形矩阵等,这类特殊矩阵在线性代数中具有通用性;还 有一类特殊矩阵在专门学科中有用,如有名的希尔伯特 Hilbert矩阵、范德蒙Vandermonde 矩阵等。 1.零矩阵:所有元素值为零的矩阵称为零矩阵。零矩阵可 以用zeros函数实现。zeros是MATLAB内部函数,使用格式 如下: zerosm:产生m m阶零矩阵; zerosm,n:产生m n阶零矩阵,当mn时等同于 zerosm; zerossizeA:产生与矩阵A同样大小的零矩阵。一、 特殊矩阵的实现 2.幺矩阵:所有元素值为1的矩阵称为幺矩阵。 幺矩阵可以用ones函数实现。它的调用格式与 zeros函数一样。 【例1】 试用ones分别建立3?2阶幺矩阵、和与 前例矩阵A同样大小的幺矩阵。 用ones3,2 建立一个3 2阶幺阵: ones3,2 % 一个3?2阶幺阵 ans 1 11 11 1一、 特殊矩阵的实现 3.单位矩阵:主对角线的元素值为1、其余元素值为0的 矩阵称为单位矩阵。它可以用MATLAB内部函数eye建立, 使用格式与zeros相同。 4.数量矩阵:主对角线的元素值为一常数d、其余元素 值为0的矩阵称为数量矩阵。显然,当d1时,即为单位 矩阵,故数量矩阵可以用eyem*d或eyem,n*d建立。 5.对角阵:对角线的元素值为常数、其余元素值为0的 矩阵称为对角阵。我们可以通过MATLAB内部函数diag, 利用一个向量构成对角阵;或从矩阵中提取某对角线构 成一个向量。使用 格式为diagV和diagV,k两种。一、 特殊矩阵的实现 6.用一个向量V构成一个对角阵 设V为具有m个元素的向量,diagV将产生一个m ?m 阶对角阵,其主对角线的元素值即为向量的元素值; diagV,k将产生一个n?nnm+|k|,k为一整数阶对角阵, 其第k条对角线的元素值即为向量的元素值。注意:当k 0,则该对角线位于主对角线的上方第k条;当k0, 该对角线位于主对角线的下方第|k|条;当k0,则等同于 diagV。用diag建立的对角阵必是方阵。 【例2】已知向量v,试建立以向量v作为主对角线的对 角阵A;建立分别以向量v作为主对角线两侧的对角线的 对角阵B和C。 MATLAB程序如下:一、 特殊矩阵的实现 % 按各种对角线情况构成相应的对角阵A、B和C v [1;2;3]; % 建立一个已知的向量A Adiagv A1 0 00 2 00 0 3 Bdiagv,1 B 0 1 0 00 0 2 00 0 0 30 0 0 0 Cdiagv,-1 C 0 0 0 01 0 0 00 2 0 00 0 3 一、 特殊矩阵的实现 7.从矩阵中提取某对角线 我们也可以用diag从矩阵中提取某对角线构成 一个向量。设A为m n阶矩阵,diagA将从矩 阵A中提取其主对角线产生一个具有minm,n个 元素的向量。diagA,k的功能是: 当k0,则将从矩阵A中提取位于主对角线的 上方第k条对角线构成一个具有n-k个元素的向 量;当k0,则将从矩阵A中提取位于主对角线 的下方第|k|条对角线构成一个具有m+k个元素的 向量;当k0,则等同于diagA。一、 特殊矩阵的实现 【例3】 已知矩阵A,试从矩阵A分别提取主 对角线及它两侧的对角线构成向量B、C和D。 MATLAB程序如下: A[1 2 3;4 5 6]; % 建立一个已知的2 ?3阶矩阵A % 按各种对角线情况构成向量B、C和D BdiagA B 15 CdiagA,1 C 26 DdiagA,-1 D 4一、 特殊矩阵的实现 8.上三角阵:使用格式为triuA、triuA,k 设A为m ?n阶矩阵,triuA将从矩阵A中提取主 对角线之上的上三角部分构成一个m n阶上三 角阵;triuA,k将从矩阵A中提取主对角线第|k|条 对角线之上的上三角部分构成一个