线性代数中的数值计算.pptVIP

第4章线性代数中的数值计算 4.1 特殊矩阵的生成 4.1.1 通用的特殊矩阵 ●zeros函数：产生全0矩阵，即零矩阵。 ●ones函数：产生全1矩阵，即幺矩阵。 ●eye函数：产生单位矩阵，即对角线上的元素为1、其余元素为0的矩阵。 ●rand函数：产生0～1均匀分布的随机矩阵。 ●randn函数：产生均值为0、方差为1的标准正态分布随机矩阵。 这几个函数的调用格式相似： 如果参数是一个，那么创建一个方阵,行数和列数均为这个数； 如果参数有两个，第一个参数代表行数，第二个参数代表列数。 zeros函数的调用格式如下。 ●zeros(m)：产生m?×?m零矩阵。 ●zeros(m,n)：产生m?×?n零矩阵。当m?=?n时，等同于zeros(m)。 ●zeros(size(A))：产生与矩阵A同样大小的零矩阵。 【例4.1】分别建立3?×?3、3?×?2和与矩阵A同样大小的零矩阵。 （1）建立一个3?×?3的零矩阵。 zeros(3) ans= 0 0 0 0 0 0 0 0 0 （2）建立一个2?×?3的零矩阵。 zeros(2,3) （3）设A为2?×?3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小的零矩阵。 A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵A zeros(size(A)) %产生一个与矩阵A同样大小的零矩阵 【例4.2】建立随机矩阵： （1）在区间[10,?30]内均匀分布的4X4随机矩阵。 （2）均值为0.6、方差为0.1的4X4正态分布随机矩阵。 若想得到在任意[a,?b]区间上均匀分布的随机数，可用yi?=?a?+?(b???a)xi 计算. xi为(0,1)区间均匀分布的随机数矩阵，可用rand函数产生。 要得到均值为μ、方差为σ的随机数，可用yi?=?μ?+?σxi计算出来。 xi为均值为0、方差为1的标准正态分布随机数矩阵，可用randn 函数。 a=10; b=30; x=a+(b-a)*rand(4) y=0.6+sqrt(0.1)*randn(4) 4.2 矩阵分析 4.2.1 矩阵结构变换 1．对角阵 （1）提取矩阵的对角线元素 设A为m?×?n矩阵，函数diag(A)用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。 A=[1,2,3;4,5,6] A=? 1 2 3 4 5 6 D=diag(A) D= 1 5 diag(A,k) 提取第k条对角线的元素。 主对角线为第0条对角线；与主对角线平行， 往上为第1条，第2条，…，第n条对角线， 往下为第???1条，第???2条，…，第???n条对角线。 （2）构造对角矩阵 diag(V,k)的功能是产生一个n?×?n（n?=?m?+?|k|）对角阵，其第k条对角线的元素即为向量V的元素。 例如： diag(1:3,-1) ans =?0 0 0 0 1 0 0 0 0 2 0 0 0 0 3 0 省略k时，相当于k为0，其主对角线元素即为向量V的元素。 2．三角阵 上三角阵是对角线以下的元素全为0的一种矩阵. 下三角阵是对角线以上的元素全为0的一种矩阵。 求矩阵A的上三角阵的MATLAB函数是triu(A)。 triu(A,k)，其功能是求矩阵A的第k条对角线以上的元素。 提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k) 3．矩阵的转置 转置运算符是单撇号()。 4．矩阵的旋转 利用函数rot90(A,k)将矩阵A旋转90°的k倍，当k为负整数时，对矩阵A按顺时针方向进行旋转；当k为1时可省略。 A=[9,37,38;-2,31,8;0,84,5]; B=rot90(A) B= 38 8 5 37 31 84 9 -2 0 ? rot90(A,4) ans =9 37 38 -2 31 8 0 84 5 5．矩阵的翻转 左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换……依此类推。 对矩阵A实施左右翻转的函数是fliplr(A)。 对矩阵A实施上下翻转的函数是flipud(A)。 4.2.2 矩阵求值 1．方阵的行列式值 求方阵A所对应的行列式的值